Resolving the 3D Landscape of Transcription-Linked Mammalian Chromatin Folding

BRD2 compartmentalizes the accessible genome

Mammalian chromosomes are organized into megabase-sized compartments that are further subdivided into topologically associated domains (TADs). While the formation of TADs is dependent on Cohesin, the mechanism behind compartmentalization remains enigmatic. Here, we show that the bromodomain and extraterminal (BET) family scaffold protein BRD2 promotes spatial mixing and compartmentalization of active chromatin after Cohesin loss. This activity is independent of transcription but requires BRD2 to recognize acetylated targets through its double bromodomain and interact with binding partners with its low complexity domain. Notably, genome compartmentalization mediated by BRD2 is antagonized on one hand by Cohesin and on the other by the BET homolog protein BRD4, both of which inhibit BRD2 binding to chromatin. Polymer simulation of our data supports a BRD2-Cohesin interplay model of nuclear topology, where genome compartmentalization results from a competition between loop extrusion and chromatin state-specific affinity interactions.

ProbC: joint modeling of epigenome and transcriptome effects in 3D genome

BACKGROUND

Hi-C and its high nucleosome resolution variant Micro-C provide a window into the spatial packing of a genome in 3D within the cell. Even though both techniques do not directly depend on the binding of specific antibodies, previous work has revealed enriched interactions and domain structures around multiple chromatin marks; epigenetic modifications and transcription factor binding sites. However, the joint impact of chromatin marks in Hi-C and Micro-C interactions have not been globally characterized, which limits our understanding of 3D genome characteristics. An emerging question is whether it is possible to deduce 3D genome characteristics and interactions by integrative analysis of multiple chromatin marks and associate interactions to functionality of the interacting loci.

RESULT

We come up with a probabilistic method PROBC to decompose Hi-C and Micro-C interactions by known chromatin marks. PROBC is based on convex likelihood optimization, which can directly take into account both interaction existence and nonexistence. Through PROBC, we discover histone modifications (H3K27ac, H3K9me3, H3K4me3, H3K4me1) and CTCF as particularly predictive of Hi-C and Micro-C contacts across cell types and species. Moreover, histone modifications are more effective than transcription factor binding sites in explaining the genome's 3D shape through these interactions. PROBC can successfully predict Hi-C and Micro-C interactions in given species, while it is trained on different cell types or species. For instance, it can predict missing nucleosome resolution Micro-C interactions in human ES cells trained on mouse ES cells only from these 5 chromatin marks with above 0.75 AUC. Additionally, PROBC outperforms the existing methods in predicting interactions across almost all chromosomes.

CONCLUSION

Via our proposed method, we optimally decompose Hi-C interactions in terms of these chromatin marks at genome and chromosome levels. We find a subset of histone modifications and transcription factor binding sites to be predictive of both Hi-C and Micro-C interactions and TADs across human, mouse, and different cell types. Through learned models, we can predict interactions on species just from chromatin marks for which Hi-C data may be limited.

Pentad: a tool for distance-dependent analysis of Hi-C interactions within and between chromatin compartments

Background

Understanding the role of various factors in 3D genome organization is essential to determine their impact on shaping large-scale chromatin units such as euchromatin (A) and heterochromatin (B) compartments. At this level, chromatin compaction is extensively modulated when transcription and epigenetic profiles change upon cell differentiation and response to various external impacts. However, detailed analysis of chromatin contact patterns within and between compartments is complicated because of a lack of suitable computational methods.

Results

We developed a tool, Pentad, to perform calculation, visualisation and quantitative analysis of the average chromatin compartment from the Hi-C matrices in cis, trans, and specified genomic distances. As we demonstrated by applying Pentad to publicly available Hi-C datasets, it helps to reliably detect redistribution of contact frequency in the chromatin compartments and assess alterations in the compartment strength.

Conclusions

Pentad is a simple tool for the analysis of changes in chromatin compartmentalization in various biological conditions. Pentad is freely available at https://github.com/magnitov/pentad.

The glucocorticoid receptor associates with the cohesin loader NIPBL to promote long-range gene regulation

The cohesin complex is central to chromatin looping, but mechanisms by which these long-range chromatin interactions are formed and persist remain unclear. We demonstrate that interactions between a transcription factor (TF) and the cohesin loader NIPBL regulate enhancer-dependent gene activity. Using mass spectrometry, genome mapping, and single-molecule tracking methods, we demonstrate that the glucocorticoid (GC) receptor (GR) interacts with NIPBL and the cohesin complex at the chromatin level, promoting loop extrusion and long-range gene regulation. Real-time single-molecule experiments show that loss of cohesin markedly diminishes the concentration of TF molecules at specific nuclear confinement sites, increasing TF local concentration and promoting gene regulation. Last, patient-derived acute myeloid leukemia cells harboring cohesin mutations exhibit a reduced response to GCs, suggesting that the GR-NIPBL-cohesin interaction is defective in these patients, resulting in poor response to GC treatment.

Mammalian Micro-C-XL

Chromosome Conformation Capture (3C) methods are a family of sequencing-based assays to measure the three-dimensional structure of genomes, with Hi-C as the most prominent method in widespread use. The Micro-C-XL protocol is technical variant that improves the resolution and signal-to-noise ratio of the Hi-C protocol and therefore offers enhanced detection of chromatin features such as chromosome loops and fine-grained resolution of topologically associated domains. Here we describe a detailed step-by-step protocol for Micro-C-XL in mammalian cells.

Stripenn detects architectural stripes from chromatin conformation data using computer vision

Architectural stripes tend to form at genomic regions harboring genes with salient roles in cell identity and function. Therefore, the accurate identification and quantification of these features are essential for understanding lineage-specific gene regulation. Here, we present Stripenn, an algorithm rooted in computer vision to systematically detect and quantitate architectural stripes from chromatin conformation measurements using various technologies. We demonstrate that Stripenn outperforms existing methods and highlight its biological applications in the context of B and T lymphocytes. By comparing stripes across distinct cell types and different species, we find that these chromatin features are highly conserved and form at genes with prominent roles in cell-type-specific processes. In summary, Stripenn is a computational method that borrows concepts from widely used image processing techniques to demarcate and quantify architectural stripes.

Stripenn detects architectural stripes from chromatin conformation data using computer vision

Architectural stripes tend to form at genomic regions harboring genes with salient roles in cell identity and function. Therefore, the accurate identification and quantification of these features are essential for understanding lineage-specific gene regulation. Here, we present Stripenn, an algorithm rooted in computer vision to systematically detect and quantitate architectural stripes from chromatin conformation measurements using various technologies. We demonstrate that Stripenn outperforms existing methods and highlight its biological applications in the context of B and T lymphocytes. By comparing stripes across distinct cell types and different species, we find that these chromatin features are highly conserved and form at genes with prominent roles in cell-type-specific processes. In summary, Stripenn is a computational method that borrows concepts from widely used image processing techniques to demarcate and quantify architectural stripes.

CTCF supports preferentially short lamina-associated domains

More than one third of the mammalian genome is in a close association with the nuclear lamina, thus these genomic regions were termed lamina-associated domains (LADs). This association is fundamental for many aspects of chromatin biology including transcription, replication, and DNA damage repair. LADs association with the nuclear envelope is thought to be dependent on two major mechanisms: The first mechanism is the interaction between nuclear membrane proteins such as LBR with heterochromatin modifications that are enriched in LADs chromatin. The second mechanism is based on proteins that bind the borders of the LADs and support the association of the LADs with the nuclear envelope. Two factors were suggested to support the second mechanism: CCCTC-binding factor (CTCF) and YY1 based on their enriched binding to LADs borders. However, this mechanism has not been proven yet at a whole genome level. Here, to test if CTCF supports the LADs landscape, we generated melanoma cells with a partial loss of function (pLoF) of CTCF by the CRISPR-Cas9 system and determined the LADs landscape by lamin B ChIP-seq analysis. We found that under regular growth conditions, CTCF pLoF led to modest changes in the LADs landscape that included an increase in the signal of 2% of the LADs and a decrease in the signal of 8% of the LADs. However, CTCF importance for the LADs landscape was much higher upon induction of a chromatin stress. We induced chromatin stress by inhibiting RNA polymerase II, an intervention that is known to alter chromatin compaction and supercoiling. Notably, only in CTCF pLoF cells, the chromatin stress led to the dissociation of 7% of the LADs from the lamina. The CTCF-dependent LADs had almost three times shorter median length than the non-affected LADs, were enriched in CTCF binding at their borders, and were higher in their facultative-status (cell-type specific). Thus, it appears that CTCF is a key factor in facilitating the association of short facultative LADs with the nuclear lamina upon chromatin stress.

A distinct core regulatory module enforces oncogene expression in KMT2A-rearranged leukemia

In this study, Harada et al. identified the transcription factors MEF2D and IRF8 as selective transcriptional dependencies of KMT2A-rearranged AML, where MEF2D displays partially redundant functions with its paralog, MEF2C. This study illustrates a mechanism of context-specific transcriptional addiction whereby a specific AML subclass depends on a highly specialized core regulatory module to directly enforce expression of common leukemia oncogenes.

Acute myeloid leukemia with KMT2A (MLL) rearrangements is characterized by specific patterns of gene expression and enhancer architecture, implying unique core transcriptional regulatory circuitry. Here, we identified the transcription factors MEF2D and IRF8 as selective transcriptional dependencies of KMT2A-rearranged AML, where MEF2D displays partially redundant functions with its paralog, MEF2C. Rapid transcription factor degradation followed by measurements of genome-wide transcription rates and superresolution microscopy revealed that MEF2D and IRF8 form a distinct core regulatory module with a narrow direct transcriptional program that includes activation of the key oncogenes MYC, HOXA9, and BCL2. Our study illustrates a mechanism of context-specific transcriptional addiction whereby a specific AML subclass depends on a highly specialized core regulatory module to directly enforce expression of common leukemia oncogenes.

Regulating specificity in enhancer-promoter communication

Enhancers are cis-regulatory elements that can activate transcription remotely to regulate a specific pattern of a gene's expression. Genes typically have many enhancers that are often intermingled in the loci of other genes. To regulate expression, enhancers must therefore activate their correct promoter while ignoring others that may be in closer linear proximity. In this review, we discuss mechanisms by which enhancers engage with promoters, including recent findings on the role of cohesin and the Mediator complex, and how this specificity in enhancer-promoter communication is encoded. Genetic dissection of model loci, in addition to more recent findings using genome-wide approaches, highlight the core promoter sequence, its accessibility, cofactor-promoter preference, in addition to the surrounding genomic context, as key components.

Whole-genome methods to define DNA and histone accessibility and long-range interactions in chromatin

Defining the genome-wide chromatin landscape has been a goal of experimentalists for decades. Here we review highlights of these efforts, from seminal experiments showing discontinuities in chromatin structure related to gene activation to extensions of these methods elucidating general features of chromatin related to gene states by exploiting deep sequencing methods. We also review chromatin conformational capture methods to identify patterns in long-range interactions between genomic loci.

Spatial organization of transcribed eukaryotic genes

Despite the well-established role of nuclear organization in the regulation of gene expression, little is known about the reverse: how transcription shapes the spatial organization of the genome. Owing to the small sizes of most previously studied genes and the limited resolution of microscopy, the structure and spatial arrangement of a single transcribed gene are still poorly understood. Here we study several long highly expressed genes and demonstrate that they form open-ended transcription loops with polymerases moving along the loops and carrying nascent RNAs. Transcription loops can span across micrometres, resembling lampbrush loops and polytene puffs. The extension and shape of transcription loops suggest their intrinsic stiffness, which we attribute to decoration with multiple voluminous nascent ribonucleoproteins. Our data contradict the model of transcription factories and suggest that although microscopically resolvable transcription loops are specific for long highly expressed genes, the mechanisms underlying their formation could represent a general aspect of eukaryotic transcription.

Dynamic Runx1 chromatin boundaries affect gene expression in hematopoietic development

The transcription factor RUNX1 is a critical regulator of developmental hematopoiesis and is frequently disrupted in leukemia. Runx1 is a large, complex gene that is expressed from two alternative promoters under the spatiotemporal control of multiple hematopoietic enhancers. To dissect the dynamic regulation of Runx1 in hematopoietic development, we analyzed its three-dimensional chromatin conformation in mouse embryonic stem cell (ESC) differentiation cultures. Runx1 resides in a 1.1 Mb topologically associating domain (TAD) demarcated by convergent CTCF motifs. As ESCs differentiate to mesoderm, chromatin accessibility, Runx1 enhancer-promoter (E-P) interactions, and CTCF-CTCF interactions increase in the TAD, along with initiation of Runx1 expression from the P2 promoter. Differentiation to hematopoietic progenitor cells is associated with the formation of tissue-specific sub-TADs over Runx1, a shift in E-P interactions, P1 promoter demethylation, and robust expression from both Runx1 promoters. Deletion of promoter-proximal CTCF sites at the sub-TAD boundaries has no obvious effects on E-P interactions but leads to partial loss of domain structure, mildly affects gene expression, and delays hematopoietic development. Together, our analysis of gene regulation at a large multi-promoter developmental gene reveals that dynamic sub-TAD chromatin boundaries play a role in establishing TAD structure and coordinated gene expression.

Genome organization controls transcriptional dynamics during development

Past studies offer contradictory claims for the role of genome organization in the regulation of gene activity. Here, we show through high-resolution chromosome conformation analysis that the Drosophila genome is organized by two independent classes of regulatory sequences, tethering elements and insulators. Quantitative live imaging and targeted genome editing demonstrate that this two-tiered organization is critical for the precise temporal dynamics of Hox gene transcription during development. Tethering elements mediate long-range enhancer-promoter interactions and foster fast activation kinetics. Conversely, the boundaries of topologically associating domains (TADs) prevent spurious interactions with enhancers and silencers located in neighboring TADs. These two levels of genome organization operate independently of one another to ensure precision of transcriptional dynamics and the reliability of complex patterning processes.

Tethering gene regulation to chromatin organization

[Figure: see text].

NicE-C efficiently reveals open chromatin-associated chromosome interactions at high resolution

Enhancer–promoter communication is known to regulate spatiotemporal dynamics of gene expression. Several methods are available to capture enhancer–promoter interactions, but they either require large amounts of starting materials and are costly, or provide a relative low resolution in chromatin contact maps. Here, we present nicking enzyme-assisted open chromatin interaction capture (NicE-C), a method that leverages nicking enzyme–mediated open chromatin profiling and chromosome conformation capture to enable robust and cost-effective detection of open chromatin interactions at high resolution, especially enhancer–promoter interactions. Using TNF stimulation and mouse kidney aging as models, we applied NicE-C to reveal characteristics of dynamic enhancer–promoter interactions.

The stochastic nature of genome organization and function

Genomes have complex three-dimensional structures. High-resolution population-based biochemical studies over the last decade have painted a mostly static picture of the genome characterized by universal organizational features, such as chromatin domains and compartments. Yet, when analyzed at the single cell level, these architectural elements are highly variable. The heterogeneity in genome organization is in line with the inherent stochasticity of transcription that shows high variation between individual cells. We highlight recent findings on single-cell variability in genome organization and describe a framework for how the stochastic nature of chromatin organization may relate to transcription dynamics.

Pre-configuring chromatin architecture with histone modifications guides hematopoietic stem cell formation in mouse embryos

The gene activity underlying cell differentiation is regulated by a diverse set of transcription factors (TFs), histone modifications, chromatin structures and more. Although definitive hematopoietic stem cells (HSCs) are known to emerge via endothelial-to-hematopoietic transition (EHT), how the multi-layered epigenome is sequentially unfolded in a small portion of endothelial cells (ECs) transitioning into the hematopoietic fate remains elusive. With optimized low-input itChIP-seq and Hi-C assays, we performed multi-omics dissection of the HSC ontogeny trajectory across early arterial ECs (eAECs), hemogenic endothelial cells (HECs), pre-HSCs and long-term HSCs (LT-HSCs) in mouse embryos. Interestingly, HSC regulatory regions are already pre-configurated with active histone modifications as early as eAECs, preceding chromatin looping dynamics within topologically associating domains. Chromatin looping structures between enhancers and promoters only become gradually strengthened over time. Notably, RUNX1, a master TF for hematopoiesis, enriched at half of these loops is observed early from eAECs through pre-HSCs but its enrichment further increases in HSCs. RUNX1 and co-TFs together constitute a central, progressively intensified enhancer-promoter interactions. Thus, our study provides a framework to decipher how temporal epigenomic configurations fulfill cell lineage specification during development.

Enhancer-promoter interactions and transcription are largely maintained upon acute loss of CTCF, cohesin, WAPL or YY1

It remains unclear why acute depletion of CTCF (CCCTC-binding factor) and cohesin only marginally affects expression of most genes despite substantially perturbing three-dimensional (3D) genome folding at the level of domains and structural loops. To address this conundrum, we used high-resolution Micro-C and nascent transcript profiling in mouse embryonic stem cells. We find that enhancer-promoter (E-P) interactions are largely insensitive to acute (3-h) depletion of CTCF, cohesin or WAPL. YY1 has been proposed as a structural regulator of E-P loops, but acute YY1 depletion also had minimal effects on E-P loops, transcription and 3D genome folding. Strikingly, live-cell, single-molecule imaging revealed that cohesin depletion reduced transcription factor (TF) binding to chromatin. Thus, although CTCF, cohesin, WAPL or YY1 is not required for the short-term maintenance of most E-P interactions and gene expression, our results suggest that cohesin may facilitate TFs to search for and bind their targets more efficiently.

Advances in profiling chromatin architecture shed light on the regulatory dynamics underlying brain disorders

Understanding the exquisitely complex nature of the three-dimensional organization of the genome and how it affects gene regulation remains a central question in biology. Recent advances in sequencing- and imaging-based approaches in decoding the three-dimensional chromatin landscape have enabled a systematic characterization of gene regulatory architecture. In this review, we outline how chromatin architecture provides a reference atlas to predict the functional consequences of non-coding variants associated with human traits and disease. High-throughput perturbation assays such as massively parallel reporter assays (MPRA) and CRISPR-based genome engineering in combination with a reference atlas opened an avenue for going beyond observational studies to experimentally validating the regulatory principles of the genome. We conclude by providing a suggested path forward by calling attention to barriers that can be addressed for a more complete understanding of the regulatory landscape of the human brain.

The transcription factor activity gradient (TAG) model: contemplating a contact-independent mechanism for enhancer-promoter communication

How distal cis-regulatory elements (e.g., enhancers) communicate with promoters remains an unresolved question of fundamental importance. Although transcription factors and cofactors are known to mediate this communication, the mechanism by which diffusible molecules relay regulatory information from one position to another along the chromosome is a biophysical puzzle-one that needs to be revisited in light of recent data that cannot easily fit into previous solutions. Here we propose a new model that diverges from the textbook enhancer-promoter looping paradigm and offer a synthesis of the literature to make a case for its plausibility, focusing on the coactivator p300.

Mapping Mammalian 3D Genomes by Micro-C

3D genome mapping aims at connecting the physics of chromatin folding to the underlying biological events, and applications of various chromosomal conformation capture (3C) assays continue to discover critical roles of genome folding in regulating nuclear functions. To interrogate the full spectrum of chromatin folding ranging from the level of nucleosomes to full chromosomes in mammals, we developed an enhanced 3C-based method called Micro-C. The protocol employs Micrococcal nuclease (MNase) to fragment the genome, which overcomes the resolution limit of restriction enzyme-based methods, enabling the estimation of contact frequencies between proximal nucleosomes. Such improvements successfully resolve the fine-scale level of chromatin folding, including enhancer-promoter or promoter-promoter interactions, genic and nucleosomal folding, and boost the signal-to-noise ratio in detecting loops and substructures underlying TADs. In this chapter, we will thoroughly discuss the details of the Micro-C protocol and critical parameters to consider for generating high-quality Micro-C maps.

cLoops2: a full-stack comprehensive analytical tool for chromatin interactions

Investigating chromatin interactions between regulatory regions such as enhancer and promoter elements is vital for understanding the regulation of gene expression. Compared to Hi-C and its variants, the emerging 3D mapping technologies focusing on enriched signals, such as TrAC-looping, reduce the sequencing cost and provide higher interaction resolution for cis-regulatory elements. A robust pipeline is needed for the comprehensive interpretation of these data, especially for loop-centric analysis. Therefore, we have developed a new versatile tool named cLoops2 for the full-stack analysis of these 3D chromatin interaction data. cLoops2 consists of core modules for peak-calling, loop-calling, differentially enriched loops calling and loops annotation. It also contains multiple modules for interaction resolution estimation, data similarity estimation, features quantification, feature aggregation analysis, and visualization. cLoops2 with documentation and example data are open source and freely available at GitHub: https://github.com/KejiZhaoLab/cLoops2.

Extracting multi-way chromatin contacts from Hi-C data

There is a growing realization that multi-way chromatin contacts formed in chromosome structures are fundamental units of gene regulation. However, due to the paucity and complexity of such contacts, it is challenging to detect and identify them using experiments. Based on an assumption that chromosome structures can be mapped onto a network of Gaussian polymer, here we derive analytic expressions for n-body contact probabilities (n > 2) among chromatin loci based on pairwise genomic contact frequencies available in Hi-C, and show that multi-way contact probability maps can in principle be extracted from Hi-C. The three-body (triplet) contact probabilities, calculated from our theory, are in good correlation with those from measurements including Tri-C, MC-4C and SPRITE. Maps of multi-way chromatin contacts calculated from our analytic expressions can not only complement experimental measurements, but also can offer better understanding of the related issues, such as cell-line dependent assemblies of multiple genes and enhancers to chromatin hubs, competition between long-range and short-range multi-way contacts, and condensates of multiple CTCF anchors.

Predicting genome organisation and function with mechanistic modelling

Fitting-free mechanistic models based on polymer simulations predict chromatin folding in 3D by focussing on the underlying biophysical mechanisms. This class of models has been increasingly used in conjunction with experiments to study the spatial organisation of eukaryotic chromosomes. Feedback from experiments to models leads to successive model refinement and has previously led to the discovery of new principles for genome organisation. Here, we review the basis of mechanistic polymer simulations, explain some of the more recent approaches and the contexts in which they have been useful to explain chromosome biology, and speculate on how they might be used in the future.

Enhancer architecture-dependent multilayered transcriptional regulation orchestrates RA signaling-induced early lineage differentiation of ESCs

Signaling pathway-driven target gene transcription is critical for fate determination of embryonic stem cells (ESCs), but enhancer-dependent transcriptional regulation in these processes remains poorly understood. Here, we report enhancer architecture-dependent multilayered transcriptional regulation at the Halr1–Hoxa1 locus that orchestrates retinoic acid (RA) signaling-induced early lineage differentiation of ESCs. We show that both homeobox A1 (Hoxa1) and Hoxa adjacent long non-coding RNA 1 (Halr1) are identified as direct downstream targets of RA signaling and regulated by RARA/RXRA via RA response elements (RAREs). Chromosome conformation capture-based screens indicate that RA signaling promotes enhancer interactions essential for Hoxa1 and Halr1 expression and mesendoderm differentiation of ESCs. Furthermore, the results also show that HOXA1 promotes expression of Halr1 through binding to enhancer; conversely, loss of Halr1 enhances interaction between Hoxa1 chromatin and four distal enhancers but weakens interaction with chromatin inside the HoxA cluster, leading to RA signaling-induced Hoxa1 overactivation and enhanced endoderm differentiation. These findings reveal complex transcriptional regulation involving synergistic regulation by enhancers, transcription factors and lncRNA. This work provides new insight into intrinsic molecular mechanisms underlying ESC fate determination during RA signaling-induced early differentiation.

Multiple stages of evolutionary change in anthrax toxin receptor expression in humans

The advent of animal husbandry and hunting increased human exposure to zoonotic pathogens. To understand how a zoonotic disease may have influenced human evolution, we study changes in human expression of anthrax toxin receptor 2 (ANTXR2), which encodes a cell surface protein necessary for Bacillus anthracis virulence toxins to cause anthrax disease. In immune cells, ANTXR2 is 8-fold down-regulated in all available human samples compared to non-human primates, indicating regulatory changes early in the evolution of modern humans. We also observe multiple genetic signatures consistent with recent positive selection driving a European-specific decrease in ANTXR2 expression in multiple tissues affected by anthrax toxins. Our observations fit a model in which humans adapted to anthrax disease following early ecological changes associated with hunting and scavenging, as well as a second period of adaptation after the rise of modern agriculture.

Local chromatin fiber folding represses transcription and loop extrusion in quiescent cells

A longstanding hypothesis is that chromatin fiber folding mediated by interactions between nearby nucleosomes represses transcription. However, it has been difficult to determine the relationship between local chromatin fiber compaction and transcription in cells. Further, global changes in fiber diameters have not been observed, even between interphase and mitotic chromosomes. We show that an increase in the range of local inter-nucleosomal contacts in quiescent yeast drives the compaction of chromatin fibers genome-wide. Unlike actively dividing cells, inter-nucleosomal interactions in quiescent cells require a basic patch in the histone H4 tail. This quiescence-specific fiber folding globally represses transcription and inhibits chromatin loop extrusion by condensin. These results reveal that global changes in chromatin fiber compaction can occur during cell state transitions, and establish physiological roles for local chromatin fiber folding in regulating transcription and chromatin domain formation.

The macro and micro of chromosome conformation capture

The 3D organization of the genome facilitates gene regulation, replication, and repair, making it a key feature of genomic function and one that remains to be properly understood. Over the past two decades, a variety of chromosome conformation capture (3C) methods have delineated genome folding from megabase-scale compartments and topologically associating domains (TADs) down to kilobase-scale enhancer-promoter interactions. Understanding the functional role of each layer of genome organization is a gateway to understanding cell state, development, and disease. Here, we discuss the evolution of 3C-based technologies for mapping 3D genome organization. We focus on genomics methods and provide a historical account of the development from 3C to Hi-C. We also discuss ChIP-based techniques that focus on 3D genome organization mediated by specific proteins, capture-based methods that focus on particular regions or regulatory elements, 3C-orthogonal methods that do not rely on restriction digestion and proximity ligation, and methods for mapping the DNA-RNA and RNA-RNA interactomes. We consider the biological discoveries that have come from these methods, examine the mechanistic contributions of CTCF, cohesin, and loop extrusion to genomic folding, and detail the 3D genome field's current understanding of nuclear architecture. Finally, we give special consideration to Micro-C as an emerging frontier in chromosome conformation capture and discuss recent Micro-C findings uncovering fine-scale chromatin organization in unprecedented detail. This article is categorized under: Gene Expression and Transcriptional Hierarchies > Regulatory Mechanisms Gene Expression and Transcriptional Hierarchies > Gene Networks and Genomics.

RNA polymerase II is required for spatial chromatin reorganization following exit from mitosis

Mammalian chromosomes are three-dimensional entities shaped by converging and opposing forces. Mitotic cell division induces marked chromosome condensation, but following reentry into the G₁ phase of the cell cycle, chromosomes reestablish their interphase organization. Here, we tested the role of RNA polymerase II (RNAPII) in this transition using a cell line that allows its auxin-mediated degradation. In situ Hi-C showed that RNAPII is required for both compartment and loop establishment following mitosis. RNAPs often counteract loop extrusion, and in their absence, longer and more prominent loops arose. Evidence from chromatin binding, super-resolution imaging, and in silico modeling allude to these effects being a result of RNAPII-mediated cohesin loading upon G₁ reentry. Our findings reconcile the role of RNAPII in gene expression with that in chromatin architecture.

Spatial Organization of Chromatin: Emergence of Chromatin Structure During Development

Nuclei are central hubs for information processing in eukaryotic cells. The need to fit large genomes into small nuclei imposes severe restrictions on genome organization and the mechanisms that drive genome-wide regulatory processes. How a disordered polymer such as chromatin, which has vast heterogeneity in its DNA and histone modification profiles, folds into discernibly consistent patterns is a fundamental question in biology. Outstanding questions include how genomes are spatially and temporally organized to regulate cellular processes with high precision and whether genome organization is causally linked to transcription regulation. The advent of next-generation sequencing, super-resolution imaging, multiplexed fluorescent in situ hybridization, and single-molecule imaging in individual living cells has caused a resurgence in efforts to understand the spatiotemporal organization of the genome. In this review, we discuss structural and mechanistic properties of genome organization at different length scales and examine changes in higher-order chromatin organization during important developmental transitions.

The Stochastic Genome and Its Role in Gene Expression

Mammalian genomes have distinct levels of spatial organization and structure that have been hypothesized to play important roles in transcription regulation. Although much has been learned about these architectural features with ensemble techniques, single-cell studies are showing a new universal trend: Genomes are stochastic and dynamic at every level of organization. Stochastic gene expression, on the other hand, has been studied for years. In this review, we probe whether there is a causative link between the two phenomena. We specifically discuss the functionality of chromatin state, topologically associating domains (TADs), and enhancer biology in light of their stochastic nature and their specific roles in stochastic gene expression. We highlight persistent fundamental questions in this area of research.

BRD4 orchestrates genome folding to promote neural crest differentiation

Higher-order chromatin structure regulates gene expression, and mutations in proteins mediating genome folding underlie developmental disorders known as cohesinopathies. However, the relationship between three-dimensional genome organization and embryonic development remains unclear. Here we define a role for bromodomain-containing protein 4 (BRD4) in genome folding, and leverage it to understand the importance of genome folding in neural crest progenitor differentiation. Brd4 deletion in neural crest results in cohesinopathy-like phenotypes. BRD4 interacts with NIPBL, a cohesin agonist, and BRD4 depletion or loss of the BRD4-NIPBL interaction reduces NIPBL occupancy, suggesting that BRD4 stabilizes NIPBL on chromatin. Chromatin interaction mapping and imaging experiments demonstrate that BRD4 depletion results in compromised genome folding and loop extrusion. Finally, mutation of individual BRD4 amino acids that mediate an interaction with NIPBL impedes neural crest differentiation into smooth muscle. Remarkably, loss of WAPL, a cohesin antagonist, rescues attenuated smooth muscle differentiation resulting from BRD4 loss. Collectively, our data reveal that BRD4 choreographs genome folding and illustrates the relevance of balancing cohesin activity for progenitor differentiation.

Complex small-world regulatory networks emerge from the 3D organisation of the human genome

The discovery that overexpressing one or a few critical transcription factors can switch cell state suggests that gene regulatory networks are relatively simple. In contrast, genome-wide association studies (GWAS) point to complex phenotypes being determined by hundreds of loci that rarely encode transcription factors and which individually have small effects. Here, we use computer simulations and a simple fitting-free polymer model of chromosomes to show that spatial correlations arising from 3D genome organisation naturally lead to stochastic and bursty transcription as well as complex small-world regulatory networks (where the transcriptional activity of each genomic region subtly affects almost all others). These effects require factors to be present at sub-saturating levels; increasing levels dramatically simplifies networks as more transcription units are pressed into use. Consequently, results from GWAS can be reconciled with those involving overexpression. We apply this pan-genomic model to predict patterns of transcriptional activity in whole human chromosomes, and, as an example, the effects of the deletion causing the diGeorge syndrome.

Transcriptional enhancers and their communication with gene promoters

Transcriptional enhancers play a key role in the initiation and maintenance of gene expression programmes, particularly in metazoa. How these elements control their target genes in the right place and time is one of the most pertinent questions in functional genomics, with wide implications for most areas of biology. Here, we synthesise classic and recent evidence on the regulatory logic of enhancers, including the principles of enhancer organisation, factors that facilitate and delimit enhancer-promoter communication, and the joint effects of multiple enhancers. We show how modern approaches building on classic insights have begun to unravel the complexity of enhancer-promoter relationships, paving the way towards a quantitative understanding of gene control.

tagHi-C Reveals 3D Chromatin Architecture Dynamics during Mouse Hematopoiesis

Spatiotemporal chromatin reorganization during hematopoietic differentiation has not been comprehensively characterized, mainly because of the large numbers of starting cells required for current chromatin conformation capture approaches. Here, we introduce a low-input tagmentation-based Hi-C (tagHi-C) method to capture the chromatin structures of hundreds of cells. Using tagHi-C, we are able to map the spatiotemporal dynamics of chromatin structure in ten primary hematopoietic stem, progenitor, and differentiated cell populations from mouse bone marrow. Our results reveal that changes in compartment dynamics and the Rabl configuration occur during hematopoietic cell differentiation. We identify gene-body-associating domains (GADs) as general structures for highly expressed genes. Moreover, we extend the body of knowledge regarding genes influenced by genome-wide association study (GWAS) loci through spatial chromatin looping. Our study provides the tagHi-C method for studying the three-dimensional (3D) genome of a small number of cells and maps the comprehensive 3D chromatin landscape of bone marrow hematopoietic cells.

Systematic evaluation of chromosome conformation capture assays

Chromosome conformation capture (3C) assays are used to map chromatin interactions genome-wide. Chromatin interaction maps provide insights into the spatial organization of chromosomes and the mechanisms by which they fold. Hi-C and Micro-C are widely used 3C protocols that differ in key experimental parameters including cross-linking chemistry and chromatin fragmentation strategy. To understand how the choice of experimental protocol determines the ability to detect and quantify aspects of chromosome folding we have performed a systematic evaluation of 3C experimental parameters. We identified optimal protocol variants for either loop or compartment detection, optimizing fragment size and cross-linking chemistry. We used this knowledge to develop a greatly improved Hi-C protocol (Hi-C 3.0) that can detect both loops and compartments relatively effectively. In addition to providing benchmarked protocols, this work produced ultra-deep chromatin interaction maps using Micro-C, conventional Hi-C and Hi-C 3.0 for key cell lines used by the 4D Nucleome project.

This analysis systematically evaluates cross-linking chemistry and chromatin fragmentation strategies commonly used in 3C assays and introduces an improved Hi-C protocol for detecting loops and compartments.

Cohesin: behind dynamic genome topology and gene expression reprogramming

Beyond its originally discovered role tethering replicated sister chromatids, cohesin has emerged as a master regulator of gene expression. Recent advances in chromatin topology resolution and single-cell studies have revealed that cohesin has a pivotal role regulating highly dynamic chromatin interactions linked to transcription control. The dynamic association of cohesin with chromatin and its capacity to perform loop extrusion contribute to the heterogeneity of chromatin contacts. Additionally, different cohesin subcomplexes, with specific properties and regulation, control gene expression across the cell cycle and during developmental cell commitment. Here, we discuss the most recent literature in the field to highlight the role of cohesin in gene expression regulation during transcriptional shifts and its relationship with human diseases.

Seeing the forest through the trees: prioritising potentially functional interactions from Hi-C

Eukaryotic genomes are highly organised within the nucleus of a cell, allowing widely dispersed regulatory elements such as enhancers to interact with gene promoters through physical contacts in three-dimensional space. Recent chromosome conformation capture methodologies such as Hi-C have enabled the analysis of interacting regions of the genome providing a valuable insight into the three-dimensional organisation of the chromatin in the nucleus, including chromosome compartmentalisation and gene expression. Complicating the analysis of Hi-C data, however, is the massive amount of identified interactions, many of which do not directly drive gene function, thus hindering the identification of potentially biologically functional 3D interactions. In this review, we collate and examine the downstream analysis of Hi-C data with particular focus on methods that prioritise potentially functional interactions. We classify three groups of approaches: structural-based discovery methods, e.g. A/B compartments and topologically associated domains, detection of statistically significant chromatin interactions, and the use of epigenomic data integration to narrow down useful interaction information. Careful use of these three approaches is crucial to successfully identifying potentially functional interactions within the genome.

CTCF and transcription influence chromatin structure re-configuration after mitosis

During mitosis, transcription is globally attenuated and chromatin architecture is dramatically reconfigured. We exploited the M- to G1-phase progression to interrogate the contributions of the architectural factor CTCF and the process of transcription to genome re-sculpting in newborn nuclei. Depletion of CTCF during the M- to G1-phase transition alters short-range compartmentalization after mitosis. Chromatin domain boundary re-formation is impaired upon CTCF loss, but a subset of boundaries, characterized by transitions in chromatin states, is established normally. Without CTCF, structural loops fail to form, leading to illegitimate contacts between cis-regulatory elements (CREs). Transient CRE contacts that are normally resolved after telophase persist deeply into G1-phase in CTCF-depleted cells. CTCF loss-associated gains in transcription are often linked to increased, normally illegitimate enhancer-promoter contacts. In contrast, at genes whose expression declines upon CTCF loss, CTCF seems to function as a conventional transcription activator, independent of its architectural role. CTCF-anchored structural loops facilitate formation of CRE loops nested within them, especially those involving weak CREs. Transcription inhibition does not significantly affect global architecture or transcription start site-associated boundaries. However, ongoing transcription contributes considerably to the formation of gene domains, regions of enriched contacts along gene bodies. Notably, gene domains emerge in ana/telophase prior to completion of the first round of transcription, suggesting that epigenetic features in gene bodies contribute to genome reconfiguration prior to transcription. The focus on the de novo formation of nuclear architecture during G1 entry yields insights into the contributions of CTCF and transcription to chromatin architecture dynamics during the mitosis to G1-phase progression.

Acute depletion of CTCF rewires genome-wide chromatin accessibility

BACKGROUND

The transcription factor CTCF appears indispensable in defining topologically associated domain boundaries and maintaining chromatin loop structures within these domains, supported by numerous functional studies. However, acute depletion of CTCF globally reduces chromatin interactions but does not significantly alter transcription.

RESULTS

Here, we systematically integrate multi-omics data including ATAC-seq, RNA-seq, WGBS, Hi-C, Cut&Run, and CRISPR-Cas9 survival dropout screens, and time-solved deep proteomic and phosphoproteomic analyses in cells carrying auxin-induced degron at endogenous CTCF locus. Acute CTCF protein degradation markedly rewires genome-wide chromatin accessibility. Increased accessible chromatin regions are frequently located adjacent to CTCF-binding sites at promoter regions and insulator sites associated with enhanced transcription of nearby genes. In addition, we use CTCF-associated multi-omics data to establish a combinatorial data analysis pipeline to discover CTCF co-regulatory partners. We successfully identify 40 candidates, including multiple established partners. Interestingly, many CTCF co-regulators that have alterations of their respective downstream gene expression do not show changes of their own expression levels across the multi-omics measurements upon acute CTCF loss, highlighting the strength of our system to discover hidden co-regulatory partners associated with CTCF-mediated transcription.

CONCLUSIONS

This study highlights that CTCF loss rewires genome-wide chromatin accessibility, which plays a critical role in transcriptional regulation.

Linker histone H1.8 inhibits chromatin binding of condensins and DNA topoisomerase II to tune chromosome length and individualization

DNA loop extrusion by condensins and decatenation by DNA topoisomerase II (topo II) are thought to drive mitotic chromosome compaction and individualization. Here, we reveal that the linker histone H1.8 antagonizes condensins and topo II to shape mitotic chromosome organization. In vitro chromatin reconstitution experiments demonstrate that H1.8 inhibits binding of condensins and topo II to nucleosome arrays. Accordingly, H1.8 depletion in Xenopus egg extracts increased condensins and topo II levels on mitotic chromatin. Chromosome morphology and Hi-C analyses suggest that H1.8 depletion makes chromosomes thinner and longer through shortening the average loop size and reducing the DNA amount in each layer of mitotic loops. Furthermore, excess loading of condensins and topo II to chromosomes by H1.8 depletion causes hyper-chromosome individualization and dispersion. We propose that condensins and topo II are essential for chromosome individualization, but their functions are tuned by the linker histone to keep chromosomes together until anaphase.

Mitochondrial nucleoid in cardiac homeostasis: bidirectional signaling of mitochondria and nucleus in cardiac diseases

Metabolic function and energy production in eukaryotic cells are regulated by mitochondria, which have been recognized as the intracellular 'powerhouses' of eukaryotic cells for their regulation of cellular homeostasis. Mitochondrial function is important not only in normal developmental and physiological processes, but also in a variety of human pathologies, including cardiac diseases. An emerging topic in the field of cardiovascular medicine is the implication of mitochondrial nucleoid for metabolic reprogramming. This review describes the linear/3D architecture of the mitochondrial nucleoid (e.g., highly organized protein-DNA structure of nucleoid) and how it is regulated by a variety of factors, such as noncoding RNA and its associated R-loop, for metabolic reprogramming in cardiac diseases. In addition, we highlight many of the presently unsolved questions regarding cardiac metabolism in terms of bidirectional signaling of mitochondrial nucleoid and 3D chromatin structure in the nucleus. In particular, we explore novel techniques to dissect the 3D structure of mitochondrial nucleoid and propose new insights into the mitochondrial retrograde signaling, and how it regulates the nuclear (3D) chromatin structures in mitochondrial diseases.

Understanding 3D genome organization by multidisciplinary methods

Understanding how chromatin is folded in the nucleus is fundamental to understanding its function. Although 3D genome organization has been historically difficult to study owing to a lack of relevant methodologies, major technological breakthroughs in genome-wide mapping of chromatin contacts and advances in imaging technologies in the twenty-first century considerably improved our understanding of chromosome conformation and nuclear architecture. In this Review, we discuss methods of 3D genome organization analysis, including sequencing-based techniques, such as Hi-C and its derivatives, Micro-C, DamID and others; microscopy-based techniques, such as super-resolution imaging coupled with fluorescence in situ hybridization (FISH), multiplex FISH, in situ genome sequencing and live microscopy methods; and computational and modelling approaches. We describe the most commonly used techniques and their contribution to our current knowledge of nuclear architecture and, finally, we provide a perspective on up-and-coming methods that open possibilities for future major discoveries.

The essential but enigmatic regulatory role of HERVH in pluripotency

Human specific endogenous retrovirus H (HERVH) is highly expressed in both naive and primed stem cells and is essential for pluripotency. Despite the proven relationship between HERVH expression and pluripotency, there is no single definitive model for the function of HERVH. Instead, several hypotheses of a regulatory function have been put forward including HERVH acting as enhancers, long noncoding RNAs (lncRNAs), and most recently as markers of topologically associating domain (TAD) boundaries. Recently several enhancer-associated lncRNAs have been characterized, which bind to Mediator and are necessary for promoter-enhancer folding interactions. We propose a synergistic model of HERVH function combining relevant findings and discuss the current limitations for its role in regulation, including the lack of evidence for a pluripotency-associated target gene.

Inducible Protein Degradation to Understand Genome Architecture

We review exciting recent advances in protein degradation, with a focus on chromatin structure. In our analysis of the literature, we highlight studies of kinetic control of protein stability for cohesin, condensin, ATP-dependent chromatin remodeling, and pioneer transcription factors. With new connections emerging between chromatin remodeling and genome structure, we anticipate exciting developments at the intersection of these topics to be revealed in the coming years. Moreover, we pay special attention to the 20-year anniversary of PROTACs, with an overview of E3 ligase/target pairings and central questions that might lead to the next generation of PROTACs with an expanded scope and generality. While steady-state experimental measurements with constitutive genome editing are impactful, we highlight complementary approaches for rapid kinetic protein degradation to uncover early targeting functions and to understand the central determinants of genome structure-function relationships.

The 3D Genome Structure of Single Cells

The spatial organization of the genome in the cell nucleus is pivotal to cell function. However, how the 3D genome organization and its dynamics influence cellular phenotypes remains poorly understood. The very recent development of single-cell technologies for probing the 3D genome, especially single-cell Hi-C (scHi-C), has ushered in a new era of unveiling cell-to-cell variability of 3D genome features at an unprecedented resolution. Here, we review recent developments in computational approaches to the analysis of scHi-C, including data processing, dimensionality reduction, imputation for enhancing data quality, and the revealing of 3D genome features at single-cell resolution. While much progress has been made in computational method development to analyze single-cell 3D genomes, substantial future work is needed to improve data interpretation and multimodal data integration, which are critical to reveal fundamental connections between genome structure and function among heterogeneous cell populations in various biological contexts.

Enhancers navigate the three-dimensional genome to direct cell fate decisions

The activity and selectivity of transcriptional enhancers determine gene expression patterns that enable a zygote to become a complex organism. How enhancers convey regulatory information is a central conundrum in biology. Here, we discuss recent progress provided by rapidly evolving technologies in understanding enhancer-promoter interactions in the context of overall nuclear genome organization.

Cooperative DNA looping by PRC2 complexes

Polycomb repressive complex 2 (PRC2) is an essential protein complex that silences gene expression via post-translational modifications of chromatin. This paper combined homology modeling, atomistic and coarse-grained molecular dynamics simulations, and single-molecule force spectroscopy experiments to characterize both its full-length structure and PRC2-DNA interactions. Using free energy calculations with a newly parameterized protein-DNA force field, we studied a total of three potential PRC2 conformations and their impact on DNA binding and bending. Consistent with cryo-EM studies, we found that EZH2, a core subunit of PRC2, provides the primary interface for DNA binding, and its curved surface can induce DNA bending. Our simulations also predicted the C2 domain of the SUZ12 subunit to contact DNA. Multiple PRC2 complexes bind with DNA cooperatively via allosteric communication through the DNA, leading to a hairpin-like looped configuration. Single-molecule experiments support PRC2-mediated DNA looping and the role of AEBP2 in regulating such loop formation. The impact of AEBP2 can be partly understood from its association with the C2 domain, blocking C2 from DNA binding. Our study suggests that accessory proteins may regulate the genomic location of PRC2 by interfering with its DNA interactions.