Chromatin Structural Domains Research Articles

Three-dimensional genome landscape comprehensively reveals patterns of spatial gene regulation in papillary and anaplastic thyroid cancers: a study using representative cell lines for each cancer type

BackgroundSpatial chromatin structure is intricately linked with somatic aberrations, and somatic mutations of various cancer-related genes, termed co-mutations (CoMuts), occur in certain patterns during cancer initiation and progression. The functional mechanisms underlying these genetic events remain largely unclear in thyroid cancer (TC). With discrepant differentiation, papillary thyroid cancer (PTC) and anaplastic thyroid cancer (ATC) differ greatly in characteristics and prognosis. We aimed to reveal the spatial gene alterations and regulations between the two TC subtypes.MethodsWe systematically investigated and compared the spatial co-mutations between ATC (8305C), PTC (BCPAP and TPC-1), and normal thyroid cells (Nthy-ori-3–1). We constructed a framework integrating whole-genome sequencing (WGS), high-throughput chromosome conformation capture (Hi-C), and transcriptome sequencing, to systematically detect the associations between the somatic co-mutations of cancer-related genes, structural variations (SVs), copy number variations (CNVs), and high-order chromatin conformation.ResultsSpatial co-mutation hotspots were enriched around topologically associating domains (TADs) in TC. A common set of 227 boundaries were identified in both ATC and PTC, with significant overlaps between them. The spatial proximities of the co-mutated gene pairs in the two TC types were significantly greater than in the gene-level and overall backgrounds, and ATC cells had higher TAD contact frequency with CoMuts > 10 compared with PTC cells. Compared with normal thyroid cells, in ATC the number of the created novel three-dimensional chromatin structural domains increased by 10%, and the number of shifted TADs decreased by 7%. We found five TAD blocks with CoMut genes/events specific to ATC with certain mutations in genes including MAST-NSUN4, AM129B/TRUB2, COL5A1/PPP1R26, PPP1R26/GPSM1/CCDC183, and PRAC2/DLX4. For the majority of ATC and PTC cells, the HOXA10 and HIF2α signals close to the transcription start sites of CoMut genes within TADs were significantly stronger than those at the background. CNV breakpoints significantly overlapped with TAD boundaries in both TC subtypes. ATCs had more CNV losses overlapping with TAD boundaries, and noncoding SVs involved in intrachromosomal SVs, amplified inversions, and tandem duplication differed between ATC and PTC. TADs with short range were more abundant in ATC than PTC. More switches of A/B compartment types existed in ATC cells compared with PTC. Gene expression was significantly synchronized, and orchestrated by complex epigenetics and regulatory elements.ConclusionChromatin interactions and gene alterations and regulations are largely heterogeneous in TC. CNVs and complex SVs may function in the TC genome by interplaying with TADs, and are largely different between ATC and PTC. Complexity of TC genomes, which are highly organized by 3D genome-wide interactions mediating mutational and structural variations and gene activation, may have been largely underappreciated. Our comprehensive analysis may provide key evidence and targets for more customized diagnosis and treatment of TC.

MCI-frcnn: A deep learning method for topological micro-domain boundary detection.

Chromatin structural domains, or topologically associated domains (TADs), are a general organizing principle in chromatin biology. RNA polymerase II (RNAPII) mediates multiple chromatin interactive loops, tethering together as RNAPII-associated chromatin interaction domains (RAIDs) to offer a framework for gene regulation. RAID and TAD alterations have been found to be associated with diseases. They can be further dissected as micro-domains (micro-TADs and micro-RAIDs) by clustering single-molecule chromatin-interactive complexes from next-generation three-dimensional (3D) genome techniques, such as ChIA-Drop. Currently, there are few tools available for micro-domain boundary identification. In this work, we developed the MCI-frcnn deep learning method to train a Faster Region-based Convolutional Neural Network (Faster R-CNN) for micro-domain boundary detection. At the training phase in MCI-frcnn, 50 images of RAIDs from Drosophila RNAPII ChIA-Drop data, containing 261 micro-RAIDs with ground truth boundaries, were trained for 7days. Using this well-trained MCI-frcnn, we detected micro-RAID boundaries for the input new images, with a fast speed (5.26 fps), high recognition accuracy (AUROC = 0.85, mAP = 0.69), and high boundary region quantification (genomic IoU = 76%). We further applied MCI-frcnn to detect human micro-TADs boundaries using human GM12878 SPRITE data and obtained a high region quantification score (mean gIoU = 85%). In all, the MCI-frcnn deep learning method which we developed in this work is a general tool for micro-domain boundary detection.

Protein Receptors Evolved from Homologous Cohesion Modules That Self-Associated and Are Encoded by Interactive Networked Genes.

Previously, it was proposed that protein receptors evolved from self-binding peptides that were encoded by self-interacting gene segments (inverted repeats) widely dispersed in the genome. In addition, self-association of the peptides was thought to be mediated by regions of amino acid sequence similarity. To extend these ideas, special features of receptors have been explored, such as their degree of homology to other proteins, and the arrangement of their genes for clues about their evolutionary origins and dynamics in the genome. As predicted, BLASTP searches for homologous proteins detected a greater number of unique hits for queries with receptor sequences than for sequences of randomly-selected, non-receptor proteins. This suggested that the building blocks (cohesion modules) for receptors were duplicated, dispersed, and maintained in the genome, due to structure/function relationships discussed here. Furthermore, the genes coding for a representative panel of receptors participated in a larger number of gene–gene interactions than for randomly-selected genes. This could conceivably reflect a greater evolutionary conservation of the receptor genes, with their more extensive integration into networks along with inherent properties of the genes themselves. In support of the latter possibility, some receptor genes were located in active areas of adaptive gene relocation/amalgamation to form functional blocks of related genes. It is suggested that adaptive relocation might allow for their joint regulation by common promoters and enhancers, and affect local chromatin structural domains to facilitate or repress gene expression. Speculation is included about the nature of the coordinated communication between receptors and the genes that encode them.

Immunochemical Detection of Modified Cytosine Species in Lampbrush Chromatin.

The lampbrush chromosomes found in the giant nucleus or germinal vesicle (GV) of amphibian oocytes provide unique opportunities for discrete closed and open chromatin structural domains to be directly observable by simple light microscopy. Moreover, the method described here for preparing spreads of lampbrush chromatin for immunostaining enables a straightforward approach to establishing the distributions of modified nucleotides within and between structurally and functionally distinctive chromatin domains.

Integrative detection and analysis of structural variation in cancer genomes.

Structural variants (SVs) can contribute to oncogenesis through a variety of mechanisms. Despite their importance, the identification of SVs in cancer genomes remains challenging. Here, we present a framework that integrates optical mapping, high-throughput chromosome conformation capture (Hi-C), and whole-genome sequencing to systematically detect SVs in a variety of normal or cancer samples and cell lines. We identify the unique strengths of each method and demonstrate that only integrative approaches can comprehensively identify SVs in the genome. By combining Hi-C and optical mapping, we resolve complex SVs and phase multiple SV events to a single haplotype. Furthermore, we observe widespread structural variation events affecting the functions of noncoding sequences, including the deletion of distal regulatory sequences, alteration of DNA replication timing, and the creation of novel three-dimensional chromatin structural domains. Our results indicate that noncoding SVs may be underappreciated mutational drivers in cancer genomes.

Transcription instability in high-risk neuroblastoma is associated with a global perturbation of chromatin domains.

Chromosome instability has a pivotal role among the hallmarks of cancer, but its transcriptional counterpart is rarely considered a relevant factor in cell destabilization. To examine transcription instability (TIN), we first devised a metric we named TIN index and used it to evaluate TIN on a dataset containing more than 500 neuroblastoma samples. We found that metastatic tumors from high‐risk (HR) patients are characterized by significantly different TIN index values compared to low/intermediate‐risk patients. Our results indicate that the TIN index is a good predictor of neuroblastoma patient's outcome, and a related TIN index gene signature (TIN‐signature) is also able to predict the neuroblastoma patient's outcome with high confidence. Interestingly, we find that TIN‐signature genes have a strong positional association with superenhancers in neuroblastoma tumors. Finally, we show that TIN is linked to chromatin structural domains and interferes with their integrity in HR neuroblastoma patients. This novel approach to gene expression analysis broadens the perspective of genome instability investigations to include functional aspects.

Association of modified cytosines and the methylated DNA-binding protein MeCP2 with distinctive structural domains of lampbrush chromatin

The lampbrush chromosomes (LBCs) of the amphibian oocyte are the largest known chromosomes and offer the unique opportunity to visualize the distinctive closed and open chromatin structural domains with an unprecedented spatial resolution. In this study we investigated the association of DNA methylation and proteins interpreting methylation state, such as MeCP2, with LBCs. We expressed HA-tagged MeCP2 in Xenopus laevis oocytes, and we observed that while it predominantly targeted the transcriptionally inactive chromomeres, a minor fraction of HA-MeCP2 also associated with many of the transcriptionally active lateral loops of LBCs. We further demonstrated that the association of MeCP2 with LBCs is directly determined by its 5-methylcytosine-binding domain. We then defined the distribution of 5-methylcytosines (5mC) by immunostaining Xenopus and axolotl LBCs and confirmed the pattern suggested by the expression of HA-MeCP2 targeting of intense staining of the chromomeres and of many loop bases. In addition, we found that short interstitial regions of the active transcriptional units could be clearly stained for 5mC. Interestingly, these 5mC-positive regions corresponded precisely to segments of active transcription units from which RNA polymerase II and nascent transcripts were simultaneously absent. Together, these data support a model presenting MeCP2 as both a transcriptional repressor and activator.

Association of modified cytosines and the methylated DNA-binding protein MeCP2 with distinctive structural domains of lampbrush chromatin

We have investigated the association of DNA methylation and proteins interpreting methylation state with the distinctive closed and open chromatin structural domains that are directly observable in the lampbrush chromosomes (LBCs) of amphibian oocytes. To establish the distribution in LBCs of MeCP2, one of the key proteins binding 5-methylcytosine-modified DNA (5mC), we expressed HA-tagged MeCP2 constructs in Xenopus laevis oocytes. Full-length MeCP2 was predominantly targeted to the closed, transcriptionally inactive chromomere domains in a pattern proportional to chromomeric DNA density and consistent with a global role in determining chromatin state. A minor fraction of HA-MeCP2 was also found to associate with a distinctive structural domain, namely a short region at the bases of some of the extended lateral loops. Expression in oocytes of deleted constructs and of point mutants derived from Rett syndrome patients demonstrated that the association of MeCP2 with LBCs was determined by its 5mC-binding domain. We also examined more directly the distribution of 5mC by immunostaining Xenopus and axolotl LBCs and confirmed the pattern suggested by MeCP2 targeting of intense staining of the chromomeres and of some loop bases. In addition, we found in the longer loops of axolotl LBCs that short interstitial regions could also be clearly stained for 5mC. These 5mC regions corresponded precisely to unusual segments of active transcription units from which RNA polymerase II (pol II) and nascent transcripts were simultaneously absent. We also examined by immunostaining the distribution in lampbrush chromatin of the oxidized 5mC derivative, 5-hydroxymethylcytosine (5hmC). Although in general, the pattern resembled that obtained for 5mC, one antibody against 5hmC produced intense staining of restricted chromosomal foci. These foci corresponded to a third type of lampbrush chromatin domain, the transcriptionally active but less extended structures formed by clusters of genes transcribed by pol III. This raises the possibility that 5hmC may play a role in establishing the distinctive patterns of gene repression and activation that characterize specific pol III-transcribed gene families in amphibian genomes.

Visualization of chromatin domains created by the gypsy insulator of Drosophila.

Insulators might regulate gene expression by establishing and maintaining the organization of the chromatin fiber within the nucleus. Biochemical fractionation and in situ high salt extraction of lysed cells show that two known protein components of the gypsy insulator are present in the nuclear matrix. Using FISH with DNA probes located between two endogenous Su(Hw) binding sites, we show that the intervening DNA is arranged in a loop, with the two insulators located at the base. Mutations in insulator proteins, subjecting the cells to a brief heat shock, or destruction of the nuclear matrix lead to disruption of the loop. Insertion of an additional gypsy insulator in the center of the loop results in the formation of paired loops through the attachment of the inserted sequences to the nuclear matrix. These results suggest that the gypsy insulator might establish higher-order domains of chromatin structure and regulate nuclear organization by tethering the DNA to the nuclear matrix and creating chromatin loops.

Chromatin insulators and boundaries: effects on transcription and nuclear organization.

Chromatin boundaries and insulators are transcriptional regulatory elements that modulate interactions between enhancers and promoters and protect genes from silencing effects by the adjacent chromatin. Originally discovered in Drosophila, insulators have now been found in a variety of organisms, ranging from yeast to humans. They have been found interspersed with regulatory sequences in complex genes and at the boundaries between active and inactive chromatin. Insulators might modulate transcription by organizing the chromatin fiber within the nucleus through the establishment of higher-order domains of chromatin structure.

Genetic and molecular analysis of the gypsy chromatin insulator of Drosophila.

Boundary or insulator elements set up independent territories of gene activity by establishing higher order domains of chromatin structure. The gypsy retrotransposon of Drosophila contains an insulator element that represses enhancer-promoter interactions and is responsible for the mutant phenotypes caused by insertion of this element. The gypsy insulator inhibits the interaction of promoter-distal enhancers with the transcription complex without affecting the functionality of promoter-proximal enhancers; in addition, these sequences can buffer a transgene from chromosomal position effects. Two proteins have been identified that bind gypsy insulator sequences and are responsible for their effects on transcription. The suppressor of Hairy-wing [su(Hw)] protein affects enhancer function both upstream and downstream of its binding site by causing a silencing effect similar to that of heterochromatin. The modifier of mdg4 [mod(mdg4)] protein interacts with su(Hw) to transform this bi-directional repression into the polar effect characteristic of insulators. These effects seem to be modulated by changes in chromatin structure.

Large-scale chromatin structural domains within mitotic and interphase chromosomes in vivo and in vitro

Higher-order chromatin structural domains approximately 130 nm in width are observed as prominent components of both Drosophila melanogaster and human mitotic chromosomes using buffer conditions which preserve chromosome morphology as determined by light microscopic comparison with chromosomes within living cells. Spatially discrete chromatin structural domains of similar size also exist as prominent components within interphase nuclei prepared under equivalent conditions. Examination of chromosomes during the anaphase-telophase transition suggests that chromosomes decondense largely through the progressive straightening or uncoiling of these large-scale chromatin domains. A quantitative analysis of the size distribution of these higher-order domains in telophase nuclei indicated a mean width of 126 +/- 36 nm. Three-dimensional views using stereopairs of chromosomes and interphase nuclei from 0.5 micron thick sections suggest that these large-scale chromatin domains consist of 30 nm fibers packed by tight folding into larger, linear, fiber-like elements. Reduction in vitro of either polyamine or divalent cation concentrations within two different buffer systems results in a loss of these large-scale domains, with no higher-order chromatin organization evident above the 20-30 nm fiber. Under these conditions the DNA distribution within mitotic chromosomes and interphase nuclei appears significantly diffuse relative to the appearance by light microscopy within living cells, or, by electron microscopy, within cells fixed directly without permeabilization in buffer. These results suggest that these large-scale chromatin structural domains are fundamental elements of chromosome architecture in vivo.