Multiple parameters shape the 3D chromatin structure of single nuclei

The spatial organization of chromatin at the scale of topologically associating domains (TADs) and below displays large cell-to-cell variations. Up until now, how this heterogeneity in chromatin conformation is shaped by chromatin condensation, TAD insulation, and transcription has remained mostly elusive. Here, we used Hi-M, a multiplexed DNA-FISH imaging technique providing developmental timing and transcriptional status, to show that the emergence of TADs at the ensemble level partially segregates the conformational space explored by single nuclei during the early development of Drosophila embryos. Surprisingly, a substantial fraction of nuclei displayed strong insulation even before TADs emerged. Moreover, active transcription within a TAD led to minor changes to the local inter- and intra-TAD chromatin conformation in single nuclei and only weakly affected insulation to the neighboring TAD. Overall, our results indicate that multiple parameters contribute to shaping the chromatin architecture of single nuclei at the TAD scale.

Aberrant GATA2 Activation in Pediatric B-Cell Acute Lymphoblastic Leukemia

GATA2 is a transcription factor that is critical for the generation and survival of hematopoietic stem cells (HSCs). It also plays an important role in the regulation of myeloid differentiation. Accordingly, GATA2 expression is restricted to HSCs and hematopoietic progenitors as well as early erythroid cells and megakaryocytic cells. Here we identified aberrant GATA2 expression in B-cell acute lymphoblastic leukemia (B-ALL) by analyzing transcriptome sequencing data obtained from St. Jude Cloud. Differentially expressed genes upon GATA2 activation showed significantly myeloid-like transcription signature. Further analysis identified several tumor-associated genes as targets of GATA2 activation including BAG3 and EPOR. In addition, the correlation between KMT2A-USP2 fusion and GATA2 activation not only indicates a potential trans-activating mechanism of GATA2 but also suggests that GATA2 is a target of KMT2A-USP2. Furthermore, by integrating whole-genome and transcriptome sequencing data, we showed that GATA2 is also cis activated. A somatic focal deletion located in the GATA2 neighborhood that disrupts the boundaries of topologically associating domains was identified in one B-ALL patient with GATA2 activation. These evidences support the hypothesis that GATA2 could be involved in leukemogenesis of B-ALL and can be transcriptionally activated through multiple mechanisms. The findings of aberrant activation of GATA2 and its molecular function extend our understanding of transcriptional factor dysregulation in B-ALL.

Rad21 is the core subunit of the cohesin complex involved in directing genome organization

The ring-shaped cohesin complex is an important factor regulating genome structure. It is thought to mediate the formation of chromatin loops and topologically associating domains (TADs) by loop extrusion. However, the regulation of association between cohesin and chromatin is poorly understood. In this study, we directly visualized cohesin loading after up-regulation of cohesin subunit Rad21 by identifying the formation of vermicelli-like structures via live cell super-resolution imaging. We also reveal that cohesin loading can be promoted by Rad21-loader interactions and accumulated contacts were shown at TAD corners while inter-TAD interactions increased after vermicelli formation, indicating that Rad21 is an important determinant of chromatin structure. Moreover, we find that cohesin saddle on topologically associating domains by FISH assay, which is consistent with the CTCF/cohesin-anchored chromatin loop model. Importantly, expression of Rad21 is strictly controlled, and aberrant expression of Rad21 leads to the formation of Rad21 “beads” in the nucleus. In summary, our observations provided important new biological insights into the mechanism of cohesin loading and its functions.

Dynamics of CTCF and cohesin mediated chromatin looping revealed by live-cell imaging

Animal genomes are folded into loops and topologically associating domains (TADs) by CTCF and cohesin, but whether these loops are stable or dynamic is unknown. Here, we directly visualize chromatin looping at the Fbn2 TAD in mouse embryonic stem cells using super-resolution live-cell imaging and quantify looping dynamics by Bayesian inference. Our results are consistent with cohesin-mediated loop extrusion in cells, and with CTCF both stopping and stabilizing cohesin. Surprisingly, the Fbn2 loop is both rare and dynamic, with a looped fraction of ~3-6.5% and a median loop lifetime of ~10-30 minutes. Instead of a stable loop, our results establish a highly dynamic view of TADs and loops where the Fbn2 TAD exists predominantly in a partially extruded conformation. This dynamic and quantitative view of TADs may facilitate a mechanistic understanding of their functions.

Deciphering the species-level structure of topologically associating domains

Chromosome conformation capture technologies such as Hi-C have revealed a rich hierarchical structure of chromatin, with topologically associating domains (TADs) as a key organizational unit, but experimentally reported TAD architectures, currently determined separately for each cell type, are lacking for many cell/tissue types. A solution to address this issue is to integrate existing epigenetic data across cells and tissue types to develop a species-level consensus map relating genes to TADs. Here, we introduce the TAD Map, a bag-of-genes representation that we use to infer, or "impute", TAD architectures for those cells/tissues with limited Hi-C experimental data. The TAD Map enables a systematic analysis of gene coexpression induced by chromatin structure. By overlaying transcriptional data from hundreds of bulk and single-cell assays onto the TAD Map, we assess gene coexpression in TADs and find that expressed genes cluster into fewer TADs than would be expected by chance, and show that time-course and RNA velocity studies further reveal this clustering to be strongest in the early stages of cell differentiation; it is also strong in tumor cells. We provide a probabilistic model to summarize any scRNA-seq transcriptome in terms of its TAD activation profile, which we term a TAD signature, and demonstrate its value for cell type inference, cell fate prediction, and multimodal synthesis. More broadly, our work indicates that the TAD Map's comprehensive, quantitative integration of chromatin structure and scRNA-seq data should play a key role in epigenetic and transcriptomic analyses. Software availability: https://tadmap.csail.mit.edu

Complete lung agenesis caused by complex genomic rearrangements with neo-TAD formation at the SHH locus

During human organogenesis, lung development is a timely and tightly regulated developmental process under the control of a large number of signaling molecules. Understanding how genetic variants can disturb normal lung development causing different lung malformations is a major goal for dissecting molecular mechanisms during embryogenesis. Here, through exome sequencing (ES), array CGH, genome sequencing (GS) and Hi-C, we aimed at elucidating the molecular basis of bilateral isolated lung agenesis in three fetuses born to a non-consanguineous family. We detected a complex genomic rearrangement containing duplicated, triplicated and deleted fragments involving the SHH locus in fetuses presenting complete agenesis of both lungs and near-complete agenesis of the trachea, diagnosed by ultrasound screening and confirmed at autopsy following termination. The rearrangement did not include SHH itself, but several regulatory elements for lung development, such as MACS1, a major SHH lung enhancer, and the neighboring genes MNX1 and NOM1. The rearrangement incorporated parts of two topologically associating domains (TADs) including their boundaries. Hi-C of cells from one of the affected fetuses showed the formation of two novel TADs each containing SHH enhancers and the MNX1 and NOM1 genes. Hi-C together with GS indicate that the new 3D conformation is likely causative for this condition by an inappropriate activation of MNX1 included in the neo-TADs by MACS1 enhancer, further highlighting the importance of the 3D chromatin conformation in human disease.

Analysis of sub-kilobase chromatin topology reveals nano-scale regulatory interactions with variable dependence on cohesin and CTCF

Enhancers and promoters predominantly interact within large-scale topologically associating domains (TADs), which are formed by loop extrusion mediated by cohesin and CTCF. However, it is unclear whether complex chromatin structures exist at sub-kilobase-scale and to what extent fine-scale regulatory interactions depend on loop extrusion. To address these questions, we present an MNase-based chromosome conformation capture (3C) approach, which has enabled us to generate the most detailed local interaction data to date and precisely investigate the effects of cohesin and CTCF depletion on chromatin architecture. Our data reveal that cis-regulatory elements have distinct internal nano-scale structures, within which local insulation is dependent on CTCF, but which are independent of cohesin. In contrast, we find that depletion of cohesin causes a subtle reduction in longer-range enhancer-promoter interactions and that CTCF depletion can cause rewiring of regulatory contacts. Together, our data show that loop extrusion is not essential for enhancer-promoter interactions, but contributes to their robustness and specificity and to precise regulation of gene expression.

The role of epigenetic modifications, long-range contacts, enhancers and topologically associating domains in the regulation of glioma grade-specific genes

Genome-wide studies have uncovered specific genetic alterations, transcriptomic patterns and epigenetic profiles associated with different glioma types. We have recently created a unique atlas encompassing genome-wide profiles of open chromatin, histone H3K27ac and H3Kme3 modifications, DNA methylation and transcriptomes of 33 glioma samples of different grades. Here, we intersected genome-wide atlas data with topologically associating domains (TADs) and demonstrated that the chromatin organization and epigenetic landscape of enhancers have a strong impact on genes differentially expressed in WHO low grade versus high grade gliomas. We identified TADs enriched in glioma grade-specific genes and/or epigenetic marks. We found the set of transcription factors, including REST, E2F1 and NFKB1, that are most likely to regulate gene expression in multiple TADs, containing specific glioma-related genes. Moreover, many genes associated with the cell–matrix adhesion Gene Ontology group, in particular 14 PROTOCADHERINs, were found to be regulated by long-range contacts with enhancers. Presented results demonstrate the existence of epigenetic differences associated with chromatin organization driving differential gene expression in gliomas of different malignancy.