Chromatin Folding Research Articles

Type II topoisomerases shape multi-scale 3D chromatin folding in regions of positive supercoils.

Type II topoisomerases (TOP2s) resolve torsional stress accumulated during various cellular processes and are enriched at chromatin loop anchors and topologically associated domain (TAD) boundaries, where, when trapped, can lead to genomic instability promoting the formation of oncogenic fusions. Whether TOP2s relieve topological constraints at these positions and/or participate in 3D chromosome folding remains unclear. Here, we combine 3D genomics, imaging, and GapRUN, a method for the genome-wide profiling of positive supercoiling, to assess the role of TOP2s in shaping chromosome organization in human cells. Acute TOP2 depletion led to the emergence of new, large-scale contacts at the boundaries between active, positively supercoiled, and lamina-associated domains. TOP2-dependent changes at the higher-order chromatin folding were accompanied by remodeling of chromatin-nuclear lamina interactions and of gene expression, while at the chromatin loop level, TOP2 depletion predominantly remodeled transcriptionally anchored, positively supercoiled loops. We propose that TOP2s act as a fine regulator of chromosome folding at multiple scales.

Regulation of cholesterol biosynthesis by CTCF and H3K27 methylation is critical for cell migration

CTCF is a key factor in three-dimensional chromatin folding and transcriptional control that was found to affect cancer cell migration by a mechanism that is still poorly understood. To identify this mechanism, we used mouse melanoma cells with a partial loss of function (pLoF) of CTCF. We found that CTCF pLoF inhibits cell migration rate while leading to an increase in the expression of multiple enzymes in the cholesterol biosynthesis pathway along with an elevation in the cellular cholesterol level. In agreement with the cholesterol change we detected altered membrane dynamics in CTCF pLoF cells as measured by reduced formation of migrasomes, extracellular vesicles formed at the rear side of migrating cells. Inhibition of cholesterol synthesis in CTCF pLoF cells restored the cellular migration rate and migrasome formation, suggesting that CTCF supports cell migration by suppressing cholesterol synthesis. Detailed analysis of the promoter of Hmgcs1, an early enzyme in the cholesterol synthesis pathway, revealed that CTCF prevents formation of a loop between that promoter and another promoter 200 kb away. CTCF also supports PRC2 recruitment to the promoter and deposition of H3K27me3. H3K27me3 at the promoter of Hmgcs1 prevents SREBP2 binding and activation of transcription. By this mechanism, CTCF fine-tunes cholesterol levels to support cell migration. Notably, genome wide association studies suggest a link between CTCF and cholesterol-associated diseases, thus CTCF emerges as a new regulator of cholesterol biosynthesis.

TAD border deletion at the Kit locus causes tissue-specific ectopic activation of a neighboring gene

Topologically associated domains (TADs) restrict promoter-enhancer interactions, thereby maintaining the spatiotemporal pattern of gene activity. However, rearrangements of the TADs boundaries do not always lead to significant changes in the activity pattern. Here, we investigated the consequences of the TAD boundaries deletion on the expression of developmentally important genes encoding tyrosine kinase receptors: Kit, Kdr, Pdgfra. We used genome editing in mice to delete the TADs boundaries at the Kit locus and characterized chromatin folding and gene expression in pure cultures of fibroblasts, mast cells, and melanocytes. We found that although Kit is highly active in both mast cells and melanocytes, deletion of the TAD boundary between the Kit and Kdr genes results in ectopic activation only in melanocytes. Thus, the epigenetic landscape, namely the mutual arrangement of enhancers and actively transcribing genes, is important for predicting the consequences of the TAD boundaries removal. We also found that mice without a TAD border between the Kit and Kdr genes have a phenotypic manifestation of the mutation — a lighter coloration. Thus, the data obtained shed light on the principles of interaction between the 3D chromatin organization and epigenetic marks in the regulation of gene activity.

Dynamical phase transition in models that couple chromatin folding with histone modifications.

Genomic regions can acquire heritable epigenetic states through unique histone modifications, which lead to stable gene expression patterns without altering the underlying DNA sequence. However, the relationship between chromatin conformational dynamics and epigenetic stability is poorly understood. In this paper, we propose kinetic models to investigate the dynamic fluctuations of histone modifications and the spatial interactions between nucleosomes. Our model explicitly incorporates the influence of chemical modifications on the structural stability of chromatin and the contribution of chromatin contacts to the cooperative nature of chemical reactions. Through stochastic simulations and analytical theory, we have discovered distinct steady-state outcomes in different kinetic regimes, resembling a dynamical phase transition. Importantly, we have validated that the emergence of this transition, which occurs on biologically relevant timescales, is robust against variations in model design and parameters. Our findings suggest that the viscoelastic properties of chromatin and the timescale at which it transitions from a gel-like to a liquidlike state significantly impact dynamic processes that occur along the one-dimensional DNA sequence.

Imprinted lncRNA KCNQ1OT1 regulates CDKN1C expression through promoter binding and chromatin folding in pigs

Long noncoding RNAs (lncRNAs) are implicated in a number of regulatory functions in eukaryotic genomes. In humans, KCNQ1OT1 is a 91 kb imprinted lncRNA that inhibits multiple surrounding genes in cis. Among them, CDKN1C is closely related to KCNQ1OT1 and is involved in multiple epigenetic disorders. Here, we found that pigs also had a relatively conserved paternal allele expressing KCNQ1OT1 and had a shorter 5′ end (∼27 kb) compared to human KCNQ1OT1. Knockdown of KCNQ1OT1 using antisense oligonucleotides (ASO) showed that upregulation of CDKN1C expression in pigs. However, porcine KCNQ1OT1 did not affect the DNA methylation status of the CpG islands in the promoters of KCNQ1OT1 and CDKN1C. Inhibition of DNA methyltransferase using Decitabine treatment resulted in a significant increase in both KCNQ1OT1 and CDKN1C expression, suggesting that the regulation between KCNQ1OT1 and CDKN1C may not be dependent on RNA interference. Further use of chromosome conformation capture and reverse transcription-associated trap detection in the region where CDKN1C was located revealed that KCNQ1OT1 bound to the CDKN1C promoter and affected chromosome folding. Phenotypically, inhibition of KCNQ1OT1 at the cumulus-oocyte complex promoted cumulus cell transformation, and to upregulated the expression of ALPL at the early stage of osteogenic differentiation of porcine bone marrow mesenchymal stem cells. Our results confirm that the expression of KCNQ1OT1 imprinting in pigs as well as porcine KCNQ1OT1 regulates the expression of CDKN1C through direct promoter binding and chromatin folding alteration. And this regulatory mechanism played an important role in cell differentiation.

Pervasive structural heterogeneity rewires glioblastoma chromosomes to sustain patient-specific transcriptional programs

Glioblastoma multiforme (GBM) encompasses brain malignancies marked by phenotypic and transcriptional heterogeneity thought to render these tumors aggressive, resistant to therapy, and inevitably recurrent. However, little is known about how the spatial organization of GBM genomes underlies this heterogeneity and its effects. Here, we compile a cohort of 28 patient-derived glioblastoma stem cell-like lines (GSCs) known to reflect the properties of their tumor-of-origin; six of these were primary-relapse tumor pairs from the same patient. We generate and analyze 5 kbp-resolution chromosome conformation capture (Hi-C) data from all GSCs to systematically map thousands of standalone and complex structural variants (SVs) and the multitude of neoloops arising as a result. By combining Hi-C, histone modification, and gene expression data with chromatin folding simulations, we explain how the pervasive, uneven, and idiosyncratic occurrence of neoloops sustains tumor-specific transcriptional programs via the formation of new enhancer-promoter contacts. We also show how even moderately recurrent neoloops can relate to patient-specific vulnerabilities. Together, our data provide a resource for dissecting GBM biology and heterogeneity, as well as for informing therapeutic approaches.

High resolution maps of chromatin reorganization through mouse meiosis reveal novel features of the 3D meiotic structure.

When germ cells transition from the mitotic cycle into meiotic prophase I (MPI), chromosomes condense into an array of chromatin loops that are required to promote homolog pairing and genetic recombination. To identify the changes in chromosomal conformation, we isolated nuclei on a trajectory from spermatogonia to the end of MPI. At each stage along this trajectory, we built genomic interaction maps with the highest temporal and spatial resolution to date. The changes in chromatin folding coincided with a concurrent decline in mitotic cohesion and a rise in meiotic cohesin complexes. We found that the stereotypical large-scale A and B compartmentalization was lost during meiotic prophase I alongside the loss of topological associating domains (TADs). Still, local subcompartments were detected and maintained throughout meiosis. The enhanced Micro-C resolution revealed that, despite the loss of TADs, higher frequency contact sites between two loci were detectable during meiotic prophase I coinciding with CTCF bound sites. The pattern of interactions around these CTCF sites with their neighboring loci showed that CTCF sites were often anchoring the meiotic loops. Additionally, the localization of CTCF to the meiotic axes indicated that these anchors were at the base of loops. Strikingly, even in the face of the dramatic reconfiguration of interphase chromatin into a condensed loop-array, the interactions between regulatory elements remained well preserved. This establishes a potential mechanism for how the meiotic chromatin maintains active transcription within a highly structured genome. In summary, the high temporal and spatial resolution of these data revealed previously unappreciated aspects of mammalian meiotic chromatin organization.

Depletion of lamins B1 and B2 promotes chromatin mobility and induces differential gene expression by a mesoscale-motion-dependent mechanism

BackgroundB-type lamins are critical nuclear envelope proteins that interact with the three-dimensional genomic architecture. However, identifying the direct roles of B-lamins on dynamic genome organization has been challenging as their joint depletion severely impacts cell viability. To overcome this, we engineered mammalian cells to rapidly and completely degrade endogenous B-type lamins using Auxin-inducible degron technology.ResultsUsing live-cell Dual Partial Wave Spectroscopic (Dual-PWS) microscopy, Stochastic Optical Reconstruction Microscopy (STORM), in situ Hi-C, CRISPR-Sirius, and fluorescence in situ hybridization (FISH), we demonstrate that lamin B1 and lamin B2 are critical structural components of the nuclear periphery that create a repressive compartment for peripheral-associated genes. Lamin B1 and lamin B2 depletion minimally alters higher-order chromatin folding but disrupts cell morphology, significantly increases chromatin mobility, redistributes both constitutive and facultative heterochromatin, and induces differential gene expression both within and near lamin-associated domain (LAD) boundaries. Critically, we demonstrate that chromatin territories expand as upregulated genes within LADs radially shift inwards. Our results indicate that the mechanism of action of B-type lamins comes from their role in constraining chromatin motion and spatial positioning of gene-specific loci, heterochromatin, and chromatin domains.ConclusionsOur findings suggest that, while B-type lamin degradation does not significantly change genome topology, it has major implications for three-dimensional chromatin conformation at the single-cell level both at the lamina-associated periphery and the non-LAD-associated nuclear interior with concomitant genome-wide transcriptional changes. This raises intriguing questions about the individual and overlapping roles of lamin B1 and lamin B2 in cellular function and disease.

Epigenetic Regulation in Schizophrenia: Focus on Methylation and Histone Modifications in Human Studies.

Despite extensive research over the last few decades, the etiology of schizophrenia (SZ) remains unclear. SZ is a pathological disorder that is highly debilitating and deeply affects the lifestyle and minds of those affected. Several factors (one or in combination) have been reported as contributors to SZ pathogenesis, including neurodevelopmental, environmental, genetic and epigenetic factors. Deoxyribonucleic acid (DNA) methylation and post-translational modification (PTM) of histone proteins are potentially contributing epigenetic processes involved in transcriptional activity, chromatin folding, cell division and apoptotic processes, and DNA damage and repair. After establishing a summary of epigenetic processes in the context of schizophrenia, this review aims to highlight the current understanding of the role of DNA methylation and histone PTMs in this disorder and their potential roles in schizophrenia pathophysiology and pathogenesis.

From Nucleosomes to Compartments: Physicochemical Interactions Underlying Chromatin Organization.

Chromatin organization plays a critical role in cellular function by regulating access to genetic information. However, understanding chromatin folding is challenging due to its complex, multiscale nature. Significant progress has been made in studying in vitro systems, uncovering the structure of individual nucleosomes and their arrays, and elucidating the role of physicochemical forces in stabilizing these structures. Additionally, remarkable advancements have been achieved in characterizing chromatin organization in vivo, particularly at the whole-chromosome level, revealing important features such as chromatin loops, topologically associating domains, and nuclear compartments. However, bridging the gap between in vitro and in vivo studies remains challenging. The resemblance between in vitro and in vivo chromatin conformations and the relevance of internucleosomal interactions for chromatin folding in vivo are subjects of debate. This article reviews experimental and computational studies conducted at various length scales, highlighting the significance of intrinsic interactions between nucleosomes and their roles in chromatin folding in vivo.

A De Novo Noncoding RARB Variant Associated with Complex Microphthalmia Alters a Putative Regulatory Element.

Retinoic acid receptor beta (RARB) is a transcriptional regulator crucial for coordinating retinoic acid- (RA-) mediated morphogenic movements, cell growth, and differentiation during eye development. Loss- or gain-of-function RARB coding variants have been associated with microphthalmia, coloboma, and anterior segment defects. We identified a de novo variant c.157+1895G>A located within a conserved region (CR1) in the first intron of RARB in an individual with complex microphthalmia and significant global developmental delay. Based on the phenotypic overlap, we further investigated the possible effects of the variant on mRNA splicing and/or transcriptional regulation through in silico and functional studies. In silico analysis identified the possibility of alternative splicing, suggested by one out of three (HSF, SpliceAI, and MaxEntScan) splicing prediction programs, and a strong indication of regulatory function based on publicly available DNase hypersensitivity, histone modification, chromatin folding, and ChIP-seq data sets. Consistent with the predictions of SpliceAI and MaxEntScan, in vitro minigene assays showed no effect on RARB mRNA splicing. Evaluation of CR1 for a regulatory role using luciferase reporter assays in human lens epithelial cells demonstrated a significant increase in the activity of the RARB promoter in the presence of wild-type CR1. This activity was further significantly increased in the presence of CR1 carrying the c.157+1895G>A variant, suggesting that the variant may promote RARB overexpression in human cells. Induction of RARB overexpression in human lens epithelial cells resulted in increased cell proliferation and elevated expression of FOXC1, a known downstream target of RA signaling and a transcription factor whose down- and upregulation is associated with ocular phenotypes overlapping the RARB spectrum. These results support a regulatory role for the CR1 element and suggest that the de novo c.157+1895G>A variant affecting this region may alter the proper regulation of RARB and, as a result, its downstream genes, possibly leading to abnormal development.

Bottom-up data integration in polymer models of chromatin organisation

Cellular functions crucially depend on the precise execution of complex biochemical reactions taking place on the chromatin fiber in the tightly packed environment of the cell nucleus. Despite the availability of large data sets probing this process from multiple angles, bottom-up frameworks which allow the incorporation of the sequence-specific nature of biochemistry in a unified model of 3D chromatin structure remain scarce. Here we propose SEMPER (Sequence Enhanced Magnetic PolymER), a novel stochastic polymer model which naturally incorporates observational data about sequence-driven biochemical processes, such as binding of transcription factor proteins, in a 3D model of chromatin structure. We introduce a novel approximate Bayesian algorithm to quantify a posteriori the relative importance of various factors, including the polymeric nature of DNA, in determining chromatin epigenetic state, thus providing a transparent way to generate biological hypotheses. While accurate prediction of contact frequencies (a problem already extensively studied in the literature) is not our main aim, as a byproduct of the inference procedure and without additional input from the genome 3D structure, our model can predict with reasonable accuracy some notable and nontrivial conformational features of chromatin folding within the nucleus. Our work highlights the importance of introducing physically realistic statistical models for predicting chromatin states from epigenetic data and opens the way to a new class of more systematic approaches to interpreting epigenomic data.

The 3D Genome in Brain Development: An Exploration of Molecular Mechanisms and Experimental Methods.

The human brain contains multiple cell types that are spatially organized into functionally distinct regions. The proper development of the brain requires complex gene regulation mechanisms in both neurons and the non-neuronal cell types that support neuronal function. Studies across the last decade have discovered that the 3D nuclear organization of the genome is instrumental in the regulation of gene expression in the diverse cell types of the brain. In this review, we describe the fundamental biochemical mechanisms that regulate the 3D genome, and comprehensively describe in vitro and ex vivo studies on mouse and human brain development that have characterized the roles of the 3D genome in gene regulation. We highlight the significance of the 3D genome in linking distal enhancers to their target promoters, which provides insights on the etiology of psychiatric and neurological disorders, as the genetic variants associated with these disorders are primarily located in noncoding regulatory regions. We also describe the molecular mechanisms that regulate chromatin folding and gene expression in neurons. Furthermore, we describe studies with an evolutionary perspective, which have investigated features that are conserved from mice to human, as well as human gained 3D chromatin features. Although most of the insights on disease and molecular mechanisms have been obtained from bulk 3C based experiments, we also highlight other approaches that have been developed recently, such as single cell 3C approaches, as well as non-3C based approaches. In our future perspectives, we highlight the gaps in our current knowledge and emphasize the need for 3D genome engineering and live cell imaging approaches to elucidate mechanisms and temporal dynamics of chromatin interactions, respectively.

Enhancer-driven 3D chromatin domain folding modulates transcription in human mammary tumor cells.

The genome is organized in functional compartments and structural domains at the sub-megabase scale. How within these domains interactions between numerous cis-acting enhancers and promoters regulate transcription remains an open question. Here, we determined chromatin folding and composition over several hundred kb around estrogen-responsive genes in human breast cancer cell lines after hormone stimulation. Modeling of 5C data at 1.8 kb resolution was combined with quantitative 3D analysis of multicolor FISH measurements at 100 nm resolution and integrated with ChIP-seq data on transcription factor binding and histone modifications. We found that rapid estradiol induction of the progesterone gene expression occurs in the context of preexisting, cell type-specific chromosomal architectures encompassing the 90 kb progesterone gene coding region and an enhancer-spiked 5' 300 kb upstream genomic region. In response to estradiol, interactions between estrogen receptor α (ERα) bound regulatory elements are reinforced. Whereas initial enhancer-gene contacts coincide with RNA Pol 2 binding and transcription initiation, sustained hormone stimulation promotes ERα accumulation creating a regulatory hub stimulating transcript synthesis. In addition to implications for estrogen receptor signaling, we uncover that preestablished chromatin architectures efficiently regulate gene expression upon stimulation without the need for de novo extensive rewiring of long-range chromatin interactions.

Multiscale modeling reveals the ion-mediated phase separation of nucleosome core particles

Due to the vast length scale inside the cell nucleus, multiscale models are required to understand chromatin folding, structure, and dynamics and how they regulate genomic activities such as DNA transcription, replication, and repair. We study the interactions and structure of condensed phases formed by the universal building block of chromatin, the nucleosome core particle (NCP), using bottom-up multiscale coarse-grained (CG) simulations with a model extracted from all-atom MD simulations. In the presence of the multivalent cations Mg(H2O)62+ or CoHex3+, we analyze the internal structures of the NCP aggregates and the contributions of histone tails and ions to the aggregation patterns. We then derive a "super" coarse-grained (SCG) NCP model to study the macroscopic scale phase separation of NCPs. The SCG simulations show the formation of NCP aggregates with Mg(H2O)62+ concentration-dependent densities and sizes. Variation of the CoHex3+ concentrations results in highly ordered lamellocolumnar and hexagonal columnar phases in agreement with experimental data. The results give detailed insights into nucleosome interactions and for understanding chromatin folding in the cell nucleus.

Context-dependent perturbations in chromatin folding and the transcriptome by cohesin and related factors

Cohesin regulates gene expression through context-specific chromatin folding mechanisms such as enhancer–promoter looping and topologically associating domain (TAD) formation by cooperating with factors such as cohesin loaders and the insulation factor CTCF. We developed a computational workflow to explore how three-dimensional (3D) structure and gene expression are regulated collectively or individually by cohesin and related factors. The main component is CustardPy, by which multi-omics datasets are compared systematically. To validate our methodology, we generated 3D genome, transcriptome, and epigenome data before and after depletion of cohesin and related factors and compared the effects of depletion. We observed diverse effects on the 3D genome and transcriptome, and gene expression changes were correlated with the splitting of TADs caused by cohesin loss. We also observed variations in long-range interactions across TADs, which correlated with their epigenomic states. These computational tools and datasets will be valuable for 3D genome and epigenome studies.

Cryoelectron tomography reveals the multiplex anatomy of condensed native chromatin and its unfolding by histone citrullination

Nucleosome chains fold and self-associate to form higher-order structures whose internal organization is unknown. Here, cryoelectron tomography (cryo-ET) of native human chromatin reveals intrinsic folding motifs such as (1) non-uniform nucleosome stacking, (2) intermittent parallel and perpendicular orientations of adjacent nucleosome planes, and (3) a regressive nucleosome chain path, which deviates from the direct zigzag topology seen in reconstituted nucleosomal arrays. By examining the self-associated structures, we observed prominent nucleosome stacking in cis and anti-parallel nucleosome interactions, which are consistent with partial nucleosome interdigitation in trans. Histone citrullination strongly inhibits nucleosome stacking and self-association with a modest effect on chromatin folding, whereas the reconstituted arrays undergo a dramatic unfolding into open zigzag chains induced by histone citrullination. This study sheds light on the internal structure of compact chromatin nanoparticles and suggests a mechanism for how epigenetic changes in chromatin folding are retained across both open and condensed forms.

Toward single-cell epigenome analysis: A microfluidic device for isolating, stretching, and imaging individual chromosomes

We developed a microfluidic device with microstructures that can moor the ends of chromatin fibers obtained by unfolding single intact chromosomes. This mooring involves trapping both ends of the chromosome with antibody-conjugated microspheres and placing each microsphere on the microstructure. In principle, the device is structured in such a way that mooring of chromatin fibers, i.e., partially unfolded chromosomes, is achieved under uniform tension, which is necessary to obtain an accurate distribution of folded/unfolded regions and fluorescence probes on the chromatin fibers. In addition, changes in the distribution of chromatin folding structures were observed by controlling the solution conditions, and specific histone chemical modification regions were visualized by immunofluorescence staining of chromatin fibers in the elongated form with uniform tension. This method is advantageous for the epigenetic analysis of non-fragmented natural chromatin fibers.

Advances in biological functions and mechanisms of histone variants in plants.

Nucleosome is the basic subunit of chromatin, consisting of approximately 147bp DNA wrapped around a histone octamer, containing two copies of H2A, H2B, H3 and H4. A linker histone H1 can bind nucleosomes through its conserved GH1 domain, which may promote chromatin folding into higher-order structures. Therefore, the complexity of histones act importantly for specifying chromatin and gene activities. Histone variants, encoded by separate genes and characterized by only a few amino acids differences, can affect nucleosome packaging and stability, and then modify the chromatin properties. Serving as carriers of pivotal genetic and epigenetic information, histone variants have profound significance in regulating plant growth and development, response to both biotic and abiotic stresses. At present, the biological functions of histone variants in plant have become a research hotspot. Here, we summarize recent researches on the biological functions, molecular chaperons and regulatory mechanisms of histone variants in plant, and propose some novel research directions for further study of plant histone variants research field. Our study will provide some enlightens for studying and understanding the epigenetic regulation and chromatin specialization mediated by histone variant in plant.

Efficient Hi-C inversion facilitates chromatin folding mechanism discovery and structure prediction

Genome-wide chromosome conformation capture (Hi-C) experiments have revealed many structural features of chromatin across multiple length scales. Further understanding genome organization requires relating these discoveries to the mechanisms that establish chromatin structures and reconstructing these structures in three dimensions, but both objectives are difficult to achieve with existing algorithms that are often computationally expensive. To alleviate this challenge, we present an algorithm that efficiently converts Hi-C data into contact energies, which measure the interaction strength between genomic loci brought into proximity. Contact energies are local quantities unaffected by the topological constraints that correlate Hi-C contact probabilities. Thus, extracting contact energies from Hi-C contact probabilities distills the biologically unique information contained in the data. We show that contact energies reveal the location of chromatin loop anchors, support a phase separation mechanism for genome compartmentalization, and parameterize polymer simulations that predict three-dimensional chromatin structures. Therefore, we anticipate that contact energy extraction will unleash the full potential of Hi-C data and that our inversion algorithm will facilitate the widespread adoption of contact energy analysis.