Hi-C analyses with GENOVA: a case study with cohesin variants

Noncoding function of super enhancer derived Cpox pre-mRNA in modulating neighbouring gene expression and chromatin interactions

Super enhancers are important regulators of gene expression that often overlap with protein-coding genes. However, it is unclear whether the overlapping protein-coding genes and the RNA derived from them contribute to enhancer activity. Using an erythroid-specific super enhancer that overlaps the Cpox gene as a model, Cpox pre-mRNA is found to have a non-coding function in regulating neighbouring protein-coding genes, eRNA expression and TAD interactions. Depletion of Cpox pre-mRNA leads to accumulation of H3K27me3 and release of p300 from the Cpox locus, activating an intra-TAD enhancer and gene expression. Additionally, a head-to-tail interaction between the TAD boundary genes Cpox and Dcbld2 is identified, facilitated by a novel type of repressive loop anchored by p300 and PRC2/H3K27me3. These results uncover a regulatory role for pre-mRNA transcribed within a super enhancer context and provide insight into head-to-tail inter-gene interaction in the regulation of gene expression and oncogene activation.

Endometrial tumorigenesis involves epigenetic plasticity demarcating non-coding somatic mutations and 3D-genome alterations

BACKGROUND

The incidence and mortality of endometrial cancer (EC) is on the rise. Eighty-five percent of ECs depend on estrogen receptor alpha (ERα) for proliferation, but little is known about its transcriptional regulation in these tumors.

RESULTS

We generate epigenomics, transcriptomics, and Hi-C datastreams in healthy and tumor endometrial tissues, identifying robust ERα reprogramming and profound alterations in 3D genome organization that lead to a gain of tumor-specific enhancer activity during EC development. Integration with endometrial cancer risk single-nucleotide polymorphisms and whole-genome sequencing data from primary tumors and metastatic samples reveals a striking enrichment of risk variants and non-coding somatic mutations at tumor-enriched ERα sites. Through machine learning-based predictions and interaction proteomics analyses, we identify an enhancer mutation which alters 3D genome conformation, impairing recruitment of the transcriptional repressor EHMT2/G9a/KMT1C, thereby alleviating transcriptional repression of ESR1 in EC.

CONCLUSIONS

In summary, we identify a complex genomic-epigenomic interplay in EC development and progression, altering 3D genome organization to enhance expression of the critical driver ERα.

H3K27me3 and the PRC1-H2AK119ub pathway cooperatively maintain heterochromatin and transcriptional silencing after the loss of H3K9 methylation

BACKGROUND

Heterochromatin is a fundamental component of eukaryotic chromosome architecture, crucial for genome stability and cell type-specific gene regulation. In mammalian nuclei, heterochromatin forms condensed B compartments, distinct from the transcriptionally active euchromatic A compartments. Histone H3 lysine 9 and lysine 27 trimethylation (H3K9me3 and H3K27me3) are two major epigenetic modifications that enrich constitutive and facultative heterochromatin, respectively. Previously, we found that the redistribution of H3K27me3 following the loss of H3K9 methylation contributes to heterochromatin maintenance, while the simultaneous loss of both H3K27me3 and H3K9 methylation induces heterochromatin decondensation in mouse embryonic fibroblasts. However, the spatial positioning of B compartments largely persists, suggesting additional mechanisms are involved.

RESULTS

In this study, we investigated the role of H2AK119 monoubiquitylation (uH2A), a repressive chromatin mark deposited by Polycomb Repressive Complex 1 (PRC1), in maintaining heterochromatin structure following the loss of H3K9 and H3K27 methylation. We observed that uH2A and H3K27me3 are independently enriched in B compartments after H3K9 methylation loss. Despite the absence of H3K9me3 and H3K27me3, uH2A remained localized and contributed to heterochromatin retention. These results suggest that PRC1-mediated uH2A functions independently and cooperatively with H3K27me3 to maintain heterochromatin organization originally created by H3K9 methylation.

CONCLUSION

Our findings highlight a compensatory role for uH2A in preserving heterochromatin structure after the loss of other repressive chromatin modifications. The PRC1-uH2A pathway plays a critical role in maintaining the integrity of B compartments and suggests that heterochromatin architecture is supported by a network of redundant epigenetic mechanisms in mammalian cells.

Hi-C profiling in tissues reveals 3D chromatin-regulated breast tumor heterogeneity informing a looping-mediated therapeutic avenue

The limitations of Hi-C (high-throughput chromosome conformation capture) profiling in in vitro cell culture include failing to recapitulate disease-specific physiological properties and lacking a clinically relevant disease microenvironment. In this study, we conduct Hi-C profiling in a pilot cohort of 12 breast tissues comprising two normal tissues, five ER+ breast primary tumors, and five tamoxifen-treated recurrent tumors. We demonstrate 3D chromatin-regulated breast tumor heterogeneity and identify a looping-mediated target gene, CA2, which might play a role in driving tamoxifen resistance. The inhibition of CA2 impedes tumor growth both in vitro and in vivo and reverses chromatin looping. The disruption of CA2 looping reduces tamoxifen-resistant cancer cell proliferation, decreases CA2 mRNA and protein expression, and weakens the looping interaction. Our study thus provides mechanistic and functional insights into the role of 3D chromatin architecture in regulating breast tumor heterogeneity and informs a new looping-mediated therapeutic avenue for treating breast cancer.

Evolutionary divergence in CTCF-mediated chromatin topology drives transcriptional innovation in humans

Chromatin topology can impact gene regulation, but how evolutionary divergence in chromatin topology has shaped gene regulatory landscapes for distinctive human traits remains poorly understood. CTCF sites determine chromatin topology by forming domains and loops. Here, we show evolutionary divergence in CTCF-mediated chromatin topology at the domain and loop scales during primate evolution, elucidating distinct mechanisms for shaping regulatory landscapes. Human-specific divergent domains lead to a broad rewiring of transcriptional landscapes. Divergent CTCF loops concord with species-specific enhancer activity, influencing enhancer connectivity to target genes in a concordant yet constrained manner. Under this concordant mechanism, we establish the role of human-specific CTCF loops in shaping transcriptional isoform diversity, with functional implications for disease susceptibility. Furthermore, we validate the function of these human-specific CTCF loops using human forebrain organoids. This study advances our understanding of genetic evolution from the perspective of genome architecture.

Binding domain mutations provide insight into CTCF's relationship with chromatin and its contribution to gene regulation

Here we used a series of CTCF mutations to explore CTCF's relationship with chromatin and its contribution to gene regulation. CTCF's impact depends on the genomic context of bound sites and the unique binding properties of WT and mutant CTCF proteins. Specifically, CTCF's signal strength is linked to changes in accessibility, and the ability to block cohesin is linked to its binding stability. Multivariate modelling reveals that both CTCF and accessibility contribute independently to cohesin binding and insulation, however CTCF signal strength has a stronger effect. CTCF and chromatin have a bidirectional relationship such that at CTCF sites, accessibility is reduced in a cohesin-dependent, mutant specific fashion. In addition, each mutant alters TF binding and accessibility in an indirect manner, changes which impart the most influence on rewiring transcriptional networks and the cell's ability to differentiate. Collectively, the mutant perturbations provide a rich resource for determining CTCF's site-specific effects.

Post-mitotic transcriptional activation and 3D regulatory interactions show locus- and differentiation-specific sensitivity to cohesin depletion

Prior studies showed that structural loops collapse upon acute cohesin depletion, while regulatory enhancer-promoter (E-P) loops largely persist, consistent with minimal transcriptional changes. However, these studies, conducted in asynchronous cells, could not resolve whether cohesin is required for the establishment of regulatory interactions and transcriptional activation during cell division or cell state transitions. To address this gap, we degraded RAD21, a core cohesin subunit, in naïve mouse embryonic stem cells (ESCs) transitioning from mitosis to G1 either in self-renewal condition or during differentiation toward formative pluripotency. Although most structural loops failed to be re-established without cohesin, about 35% of regulatory loops reformed at normal or higher frequencies. Cohesin-independent loops showed characteristics of strong active enhancers and promoters and a significant association with H3K27ac mitotic bookmarks. However, inhibition of CBP/p300 during mitotic exit did not impact these cohesin-independent interactions, suggesting the presence of complex compensatory mechanisms. At the transcriptional level, cohesin depletion induced only minor changes, supporting that post-mitotic transcriptional reactivation is largely independent of cohesin. The few genes with impaired reactivation were directly bound by RAD21 at their promoters, engaged in many structural loops, and located within strongly insulated TADs with low gene density. Importantly, degrading cohesin during the M-to-G1 transition in the presence of EpiLC differentiation signals revealed a larger group of susceptible genes, including key signature genes and transcription factors. Impaired activation of these genes was partly due to the failure to establish de novo EpiLC-specific interactions in the absence of cohesin. These experiments revealed locus-specific and context-specific dependencies between cohesin, E-P interactions, and transcription.

B cell stimulation changes the structure and higher-order organization of the inactive X chromosome

X Chromosome Inactivation (XCI) equalizes X-linked gene expression between sexes. B cells exhibit dynamic XCI, with Xist RNA/heterochromatic marks absent on the inactive X (Xi) in naive B cells but returning following mitogenic stimulation. The impact of dynamic XCI on Xi structure and maintenance was previously unknown. Here, we find dosage compensation of the Xi with state-specific XCI escape genes in naive and in vitro activated B cells. Allele-specific OligoPaints indicate similar Xi and Xa territories in B cells that are less compact than in fibroblasts. Allele-specific Hi-C reveals a lack of TAD-like structures on the Xi of naive B cells, and stimulation-induced alterations in TAD-like boundary strength independent of gene expression. Notably, Xist deletion in B cells changes TAD boundaries and large-scale Xi compaction. Altogether, our results uncover B cell-specific Xi plasticity which could underlie sex-biased biological mechanisms.

Chromatin conformation, gene transcription, and nucleosome remodeling as an emergent system

In single cells, variably sized nanoscale chromatin structures are observed, but it is unknown whether these form a cohesive framework that regulates RNA transcription. Here, we demonstrate that the human genome is an emergent, self-assembling, reinforcement learning system. Conformationally defined heterogeneous, nanoscopic packing domains form by the interplay of transcription, nucleosome remodeling, and loop extrusion. We show that packing domains are not topologically associated domains. Instead, packing domains exist across a structure-function life cycle that couples heterochromatin and transcription in situ, explaining how heterochromatin enzyme inhibition can produce a paradoxical decrease in transcription by destabilizing domain cores. Applied to development and aging, we show the pairing of heterochromatin and transcription at myogenic genes that could be disrupted by nuclear swelling. In sum, packing domains represent a foundation to explore the interactions of chromatin and transcription at the single-cell level in human health.

Epigenomic and 3D genomic mapping reveals developmental dynamics and subgenomic asymmetry of transcriptional regulatory architecture in allotetraploid cotton

Although epigenetic modification has long been recognized as a vital force influencing gene regulation in plants, the dynamics of chromatin structure implicated in the intertwined transcriptional regulation of duplicated genes in polyploids have yet to be understood. Here, we document the dynamic organization of chromatin structure in two subgenomes of allotetraploid cotton (Gossypium hirsutum) by generating 3D genomic, epigenomic and transcriptomic datasets from 12 major tissues/developmental stages covering the life cycle. We systematically identify a subset of genes that are closely associated with specific tissue functions. Interestingly, these genes exhibit not only higher tissue specificity but also a more pronounced homoeologous bias. We comprehensively elucidate the intricate process of subgenomic collaboration and divergence across various tissues. A comparison among subgenomes in the 12 tissues reveals widespread differences in the reorganization of 3D genome structures, with the Dt subgenome exhibiting a higher extent of dynamic chromatin status than the At subgenome. Moreover, we construct a comprehensive atlas of putative functional genome elements and discover that 37 cis-regulatory elements (CREs) have selection signals acquired during domestication and improvement. These data and analyses are publicly available to the research community through a web portal. In summary, this study provides abundant resources and depicts the regulatory architecture of the genome, which thereby facilitates the understanding of biological processes and guides cotton breeding.

Inter-chromosomal transcription hubs shape the 3D genome architecture of African trypanosomes

The eukaryotic nucleus exhibits a highly organized 3D genome architecture, with RNA transcription and processing confined to specific nuclear structures. While intra-chromosomal interactions, such as promoter-enhancer dynamics, are well-studied, the role of inter-chromosomal interactions remains poorly understood. Investigating these interactions in mammalian cells is challenging due to large genome sizes and the need for deep sequencing. Additionally, transcription-dependent 3D topologies in mixed cell populations further complicate analyses. To address these challenges, we used high-resolution DNA-DNA contact mapping (Micro-C) in Trypanosoma brucei, a parasite with continuous RNA polymerase II (RNAPII) transcription and polycistronic transcription units (PTUs). With approximately 300 transcription start sites (TSSs), this genome organization simplifies data interpretation. To minimize scaffolding artifacts, we also generated a highly contiguous phased genome assembly using ultra-long sequencing reads. Our Micro-C analysis revealed an intricate 3D genome organization. While the T. brucei genome displays features resembling chromosome territories, its chromosomes are arranged around polymerase-specific transcription hubs. RNAPI-transcribed genes cluster, as expected from their localization to the nucleolus. However, we also found that RNAPII TSSs form distinct inter-chromosomal transcription hubs with other RNAPII TSSs. These findings highlight the evolutionary significance of inter-chromosomal transcription hubs and provide new insights into genome organization in T. brucei.

E-box independent chromatin recruitment turns MYOD into a transcriptional repressor

MYOD is an E-box sequence-specific basic Helix-Loop-Helix (bHLH) transcriptional activator that, when expressed in non-muscle cells, induces nuclear reprogramming toward skeletal myogenesis by promoting chromatin accessibility at previously silent loci. Here, we report on the identification of a previously unrecognized property of MYOD as repressor of gene expression, via E-box-independent chromatin binding within accessible genomic elements, which invariably leads to reduced chromatin accessibility. MYOD-mediated repression requires the integrity of functional domains previously implicated in MYOD-mediated activation of gene expression. Repression of mitogen- and growth factor-responsive genes occurs through promoter binding and requires a highly conserved domain within the first helix. Repression of cell-of-origin/alternative lineage genes occurs via binding and decommissioning of distal regulatory elements, such as super-enhancers (SE), which requires the N-terminal activation domain as well as two chromatin-remodeling domains and leads to reduced strength of CTCF-mediated chromatin interactions. Surprisingly, MYOD-mediated chromatin compaction and repression of transcription do not associate with reduction of H3K27ac, the conventional histone mark of enhancer or promoter activation, but with reduced levels of the recently discovered histone H4 acetyl-methyl lysine modification (Kacme). These results extend MYOD biological properties beyond the current dogma that restricts MYOD function to a monotone transcriptional activator and reveal a previously unrecognized functional versatility arising from an alternative chromatin recruitment through E-box or non-E-box sequences. The E-box independent repression of gene expression by MYOD might provide a promiscuous mechanism to reduce chromatin accessibility and repress cell-of-origin/alternative lineage and growth factor/mitogen-responsive genes to safeguard the integrity of cell identity during muscle progenitor commitment toward the myogenic lineage.

Gene mobility elements mediate cell type specific genome organization and radial gene movement in vivo

Understanding the level of genome organization that governs gene regulation remains a challenge despite advancements in chromatin profiling techniques. Cell type specific chromatin architectures may be obscured by averaging heterogeneous cell populations. Here we took a reductionist perspective, starting with the relocation of the hunchback gene to the nuclear lamina in Drosophila neuroblasts. We previously found that this event terminates competence to produce early-born neurons and is mediated by an intronic 250 base-pair element, which we term gene mobility element (GME). Here we found over 800 putative GMEs globally that are chromatin accessible and are Polycomb (PcG) target sites. GMEs appear to be distinct from PcG response elements, however, which are largely chromatin inaccessible in neuroblasts. Performing in situ Hi-C of purified neuroblasts, we found that GMEs form megabase-scale chromatin interactions, spanning multiple topologically associated domain borders, preferentially contacting other GMEs. These interactions are cell type and stage-specific. Notably, GMEs undergo developmentally- timed mobilization to/from the neuroblast nuclear lamina, and domain swapping a GFP reporter transgene intron with a GME relocates the transgene to the nuclear lamina in embryos. We propose that GMEs constitute a genome organizational framework and mediate gene-to-lamina mobilization during progenitor competence state transitions in vivo .

STAG2 loss in Ewing sarcoma alters enhancer-promoter contacts dependent and independent of EWS::FLI1

Cohesin complexes carrying STAG1 or STAG2 organize the genome into chromatin loops. STAG2 loss-of-function mutations promote metastasis in Ewing sarcoma, a pediatric cancer driven by the fusion transcription factor EWS::FLI1. We integrated transcriptomic data from patients and cellular models to identify a STAG2-dependent gene signature associated with worse prognosis. Subsequent genomic profiling and high-resolution chromatin interaction data from Capture Hi-C indicated that cohesin-STAG2 facilitates communication between EWS::FLI1-bound long GGAA repeats, presumably acting as neoenhancers, and their target promoters. Changes in CTCF-dependent chromatin contacts involving signature genes, unrelated to EWS::FLI1 binding, were also identified. STAG1 is unable to compensate for STAG2 loss and chromatin-bound cohesin is severely decreased, while levels of the processivity factor NIPBL remain unchanged, likely affecting DNA looping dynamics. These results illuminate how STAG2 loss modifies the chromatin interactome of Ewing sarcoma cells and provide a list of potential biomarkers and therapeutic targets.

PWOs repress gene transcription by regulating chromatin structures in Arabidopsis

PWWP-DOMAIN INTERACTOR OF POLYCOMBS (PWO) family proteins play a vital role in regulating plant development. However, the molecular mechanisms of how PWOs regulate chromatin structure is elusive. Our data show that the PWO1 binding sites are enriched with positive modifications but exclusive with H3K27me3. Moreover, PWO1 binds to the H3K27me3-enriched compartment domain (H3K27me3-CD) boundary regions, and functions to maintain the boundary strength. Meanwhile, we found that PWOs and Polycomb repressive complex 2 (PRC2) function parallelly in maintaining H3K27me3-CDs' structure. Loss of either PWOs or PRC2 leads to H3K27me3-CD strength reduction, B to A compartment switching as well as the H3K27me3-CD relocating away from the nuclear periphery. Additionally, PWOs and lamin-like proteins collaborate to regulate multiple chromatin structures to repress gene transcription within H3K27me3-CDs. We conclude that PWOs maintain H3K27me3-CDs' repressive state and regulate their spatial position in the nucleus.

A genome wide code to define cell-type specific CTCF binding and chromatin organization

CTCF-mediated chromatin folding plays a key role in gene regulation, however the mechanisms underlying cell type-specific control are not fully elucidated. Comprehensive analyses reveal that CTCF binding stability and cohesin overlap in mice and humans is regulated by species specific differences in accessibility, the presence of CTCF extended Up and Downstream binding site sequences and motifs corresponding to expressed TFs enriched at most bound sites. CTCF repositions nucleosomes asymmetrically, with linker spacing altered by accessibility, while cohesin-mediated nucleosome phasing is affected by surrounding motifs that similarly impact chromatin insulation. Importantly, cell type-specific transcriptional programs determine motif enrichment at CTCF bound sites, reflecting the stabilizing/destabilizing effect of individual TFs. These studies identify mechanisms underlying cell type-specific CTCF profiles, linked to local and long-range chromatin organization.

Multiscale 3D genome rewiring during PTF1A-mediated somatic cell reprogramming into neural stem cells

The genome is intricately folded into chromatin compartments, topologically associating domains (TADs) and loops unique to each cell type. How this higher-order genome organization regulates cell fate transition remains elusive. Here we show how a single non-neural progenitor transcription factor, PTF1A, reorchestrates the 3D genome during fibroblast transdifferentiation into neural stem cells (NSCs). Multiomics analyses integrating Hi-C data, PTF1A and CTCF DNA-binding profiles, H3K27ac modification, and gene expression, demonstrate that PTF1A binds to subTAD boundaries subsequently associated with elevated CTCF binding and enhanced boundary insulation, and reorganizes chromatin loops, leading to gene expression changes that drive transdifferentiation into NSCs. Moreover, PTF1A activates enhancers and super-enhancers near low-insulation boundaries and modulates H3K27ac deposition, promoting cell fate transitions. Together, our data implicate an involvement of 3D genome in transcriptional and cell fate alterations, and highlight an essential role for PTF1A in gene expression control and multiscale 3D genome remodeling during cell reprogramming.

Heterochromatin protein 1 alpha (HP1α) undergoes a monomer to dimer transition that opens and compacts live cell genome architecture

Our understanding of heterochromatin nanostructure and its capacity to mediate gene silencing in a living cell has been prevented by the diffraction limit of optical microscopy. Thus, here to overcome this technical hurdle, and directly measure the nucleosome arrangement that underpins this dense chromatin state, we coupled fluorescence lifetime imaging microscopy (FLIM) of Förster resonance energy transfer (FRET) between histones core to the nucleosome, with molecular editing of heterochromatin protein 1 alpha (HP1α). Intriguingly, this super-resolved readout of nanoscale chromatin structure, alongside fluorescence fluctuation spectroscopy (FFS) and FLIM-FRET analysis of HP1α protein-protein interaction, revealed nucleosome arrangement to be differentially regulated by HP1α oligomeric state. Specifically, we found HP1α monomers to impart a previously undescribed global nucleosome spacing throughout genome architecture that is mediated by trimethylation on lysine 9 of histone H3 (H3K9me3) and locally reduced upon HP1α dimerisation. Collectively, these results demonstrate HP1α to impart a dual action on chromatin that increases the dynamic range of nucleosome proximity. We anticipate that this finding will have important implications for our understanding of how live cell heterochromatin structure regulates genome function.

Identifying genetic variants associated with chromatin looping and genome function

Here we present a comprehensive HiChIP dataset on naïve CD4 T cells (nCD4) from 30 donors and identify QTLs that associate with genotype-dependent and/or allele-specific variation of HiChIP contacts defining loops between active regulatory regions (iQTLs). We observe a substantial overlap between iQTLs and previously defined eQTLs and histone QTLs, and an enrichment for fine-mapped QTLs and GWAS variants. Furthermore, we describe a distinct subset of nCD4 iQTLs, for which the significant variation of chromatin contacts in nCD4 are translated into significant eQTL trends in CD4 T cell memory subsets. Finally, we define connectivity-QTLs as iQTLs that are significantly associated with concordant genotype-dependent changes in chromatin contacts over a broad genomic region (e.g., GWAS SNP in the RNASET2 locus). Our results demonstrate the importance of chromatin contacts as a complementary modality for QTL mapping and their power in identifying previously uncharacterized QTLs linked to cell-specific gene expression and connectivity.

Clusters of lineage-specific genes are anchored by ZNF274 in repressive perinucleolar compartments

Long known as the site of ribosome biogenesis, the nucleolus is increasingly recognized for its role in shaping three-dimensional (3D) genome organization. Still, the mechanisms governing the targeting of selected regions of the genome to nucleolus-associated domains (NADs) remain enigmatic. Here, we reveal the essential role of ZNF274, a SCAN-bearing member of the Krüppel-associated box (KRAB)-containing zinc finger protein (KZFP) family, in sequestering lineage-specific gene clusters within NADs. Ablation of ZNF274 triggers transcriptional activation across entire genomic neighborhoods-encompassing, among others, protocadherin and KZFP-encoding genes-with loss of repressive chromatin marks, altered the 3D genome architecture and de novo CTCF binding. Mechanistically, ZNF274 anchors target DNA sequences at the nucleolus and facilitates their compartmentalization via a previously uncharted function of the SCAN domain. Our findings illuminate the mechanisms underlying NAD organization and suggest that perinucleolar entrapment into repressive hubs constrains the activation of tandemly arrayed genes to enable selective expression and modulate cell differentiation programs during development.

Reversible acetylation of HDAC8 regulates cell cycle

HDAC8, a member of class I HDACs, plays a pivotal role in cell cycle regulation by deacetylating the cohesin subunit SMC3. While cyclins and CDKs are well-established cell cycle regulators, our knowledge of other regulators remains limited. Here we reveal the acetylation of K202 in HDAC8 as a key cell cycle regulator responsive to stress. K202 acetylation in HDAC8, primarily catalyzed by Tip60, restricts HDAC8 activity, leading to increased SMC3 acetylation and cell cycle arrest. Furthermore, cells expressing the mutant form of HDAC8 mimicking K202 acetylation display significant alterations in gene expression, potentially linked to changes in 3D genome structure, including enhanced chromatid loop interactions. K202 acetylation impairs cell cycle progression by disrupting the expression of cell cycle-related genes and sister chromatid cohesion, resulting in G2/M phase arrest. These findings indicate the reversible acetylation of HDAC8 as a cell cycle regulator, expanding our understanding of stress-responsive cell cycle dynamics.

Condensin I folds the Caenorhabditis elegans genome

The structural maintenance of chromosome (SMC) complexes-cohesin and condensins-are crucial for chromosome separation and compaction during cell division. During the interphase, mammalian cohesins additionally fold the genome into loops and domains. Here we show that, in Caenorhabditis elegans, a species with holocentric chromosomes, condensin I is the primary, long-range loop extruder. The loss of condensin I and its X-specific variant, condensin IDC, leads to genome-wide decompaction, chromosome mixing and disappearance of X-specific topologically associating domains, while reinforcing fine-scale epigenomic compartments. In addition, condensin I/IDC inactivation led to the upregulation of X-linked genes and unveiled nuclear bodies grouping together binding sites for the X-targeting loading complex of condensin IDC. C. elegans condensin I/IDC thus uniquely organizes holocentric interphase chromosomes, akin to cohesin in mammals, as well as regulates X-chromosome gene expression.

Brain and cancer associated binding domain mutations provide insight into CTCF's relationship with chromatin and its ability to act as a chromatin organizer

Although only a fraction of CTCF motifs are bound in any cell type, and approximately half of the occupied sites overlap cohesin, the mechanisms underlying cell-type specific attachment and ability to function as a chromatin organizer remain unknown. To investigate the relationship between CTCF and chromatin we applied a combination of imaging, structural and molecular approaches, using a series of brain and cancer associated CTCF mutations that act as CTCF perturbations. We demonstrate that binding and the functional impact of WT and mutant CTCF depend not only on the unique properties of each protein, but also on the genomic context of bound sites. Our studies also highlight the reciprocal relationship between CTCF and chromatin, demonstrating that the unique binding properties of WT and mutant proteins have a distinct impact on accessibility, TF binding, cohesin overlap, chromatin interactivity and gene expression programs, providing insight into their cancer and brain related effects.

SMCHD1 maintains heterochromatin and genome compartments in human myoblasts

Mammalian genomes are subdivided into euchromatic A compartments that contain mostly active chromatin, and inactive, heterochromatic B compartments. However, it is unknown how A and B genome compartments are established and maintained. Here we studied SMCHD1, an SMC-like protein in human male myoblasts. SMCHD1 colocalizes with Lamin B1 and the heterochromatin mark H3K9me3. Loss of SMCHD1 leads to extensive heterochromatin depletion at the nuclear lamina and acquisition of active chromatin states along all chromosomes. In absence of SMCHD1, long range intra-chromosomal and inter-chromosomal contacts between B compartments are lost while many new TADs and loops are formed. Inactivation of SMCHD1 promotes numerous B to A compartment transitions accompanied by activation of silenced genes. SMCHD1 functions as an anchor for heterochromatin domains ensuring that these domains are inaccessible to epigenome modification enzymes that typically operate in active chromatin. Therefore, A compartments are formed by default when not prevented by SMCHD1.

An atlas of the tomato epigenome reveals that KRYPTONITE shapes TAD-like boundaries through the control of H3K9ac distribution

In recent years, the exploration of genome three-dimensional (3D) conformation has yielded profound insights into the regulation of gene expression and cellular functions in both animals and plants. While animals exhibit a characteristic genome topology defined by topologically associating domains (TADs), plants display similar features with a more diverse conformation across species. Employing advanced high-throughput sequencing and microscopy techniques, we investigated the landscape of 26 histone modifications and RNA polymerase II distribution in tomato (Solanum lycopersicum). Our study unveiled a rich and nuanced epigenetic landscape, shedding light on distinct chromatin states associated with heterochromatin formation and gene silencing. Moreover, we elucidated the intricate interplay between these chromatin states and the overall topology of the genome. Employing a genetic approach, we delved into the role of the histone modification H3K9ac in genome topology. Notably, our investigation revealed that the ectopic deposition of this chromatin mark triggered a reorganization of the 3D chromatin structure, defining different TAD-like borders. Our work emphasizes the critical role of H3K9ac in shaping the topology of the tomato genome, providing valuable insights into the epigenetic landscape of this agriculturally significant crop species.

PHF2 regulates genome topology and DNA replication in neural stem cells via cohesin

Cohesin plays a crucial role in the organization of topologically-associated domains (TADs), which influence gene expression and DNA replication timing. Whether epigenetic regulators may affect TADs via cohesin to mediate DNA replication remains elusive. Here, we discover that the histone demethylase PHF2 associates with RAD21, a core subunit of cohesin, to regulate DNA replication in mouse neural stem cells (NSC). PHF2 loss impairs DNA replication due to the activation of dormant replication origins in NSC. Notably, the PHF2/RAD21 co-bound genomic regions are characterized by CTCF enrichment and epigenomic features that resemble efficient, active replication origins, and can act as boundaries to separate adjacent domains. Accordingly, PHF2 loss weakens TADs and chromatin loops at the co-bound loci due to reduced RAD21 occupancy. The observed topological and DNA replication defects in PHF2 KO NSC support a cohesin-dependent mechanism. Furthermore, we demonstrate that the PHF2/RAD21 complex exerts little effect on gene regulation, and that PHF2's histone-demethylase activity is dispensable for normal DNA replication and proliferation of NSC. We propose that PHF2 may serve as a topological accessory to cohesin for cohesin localization to TADs and chromatin loops, where cohesin represses dormant replication origins directly or indirectly, to sustain DNA replication in NSC.

SETDB1 regulates short interspersed nuclear elements and chromatin loop organization in mouse neural precursor cells

BACKGROUND

Transposable elements play a critical role in maintaining genome architecture during neurodevelopment. Short Interspersed Nuclear Elements (SINEs), a major subtype of transposable elements, are known to harbor binding sites for the CCCTC-binding factor (CTCF) and pivotal in orchestrating chromatin organization. However, the regulatory mechanisms controlling the activity of SINEs in the developing brain remains elusive.

RESULTS

In our study, we conduct a comprehensive genome-wide epigenetic analysis in mouse neural precursor cells using ATAC-seq, ChIP-seq, whole genome bisulfite sequencing, in situ Hi-C, and RNA-seq. Our findings reveal that the SET domain bifurcated histone lysine methyltransferase 1 (SETDB1)-mediated H3K9me3, in conjunction with DNA methylation, restricts chromatin accessibility on a selective subset of SINEs in neural precursor cells. Mechanistically, loss of Setdb1 increases CTCF access to these SINE elements and contributes to chromatin loop reorganization. Moreover, de novo loop formation contributes to differential gene expression, including the dysregulation of genes enriched in mitotic pathways. This leads to the disruptions of cell proliferation in the embryonic brain after genetic ablation of Setdb1 both in vitro and in vivo.

CONCLUSIONS

In summary, our study sheds light on the epigenetic regulation of SINEs in mouse neural precursor cells, suggesting their role in maintaining chromatin organization and cell proliferation during neurodevelopment.

Estrogen receptor 1 chromatin profiling in human breast tumors reveals high inter-patient heterogeneity with enrichment of risk SNPs and enhancer activity at most-conserved regions

Estrogen Receptor 1 (ESR1; also known as ERα, encoded by ESR1 gene) is the main driver and prime drug target in luminal breast cancer. ESR1 chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ESR1 chromatin action, along with its biological implications. Here, we use a large set of ESR1 ChIP-seq data from 70 ESR1+ breast cancers to explore inter-patient heterogeneity in ESR1 DNA binding to reveal a striking inter-tumor heterogeneity of ESR1 action. Of note, commonly shared ESR1 sites show the highest estrogen-driven enhancer activity and are most engaged in long-range chromatin interactions. In addition, the most commonly shared ESR1-occupied enhancers are enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ESR1 and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we can confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ESR1-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ESR1 landscape, with the most common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.

Open2C, Abdennur N, Abraham S, Fudenberg G, Flyamer IM, Galitsyna AA, Goloborodko A, Imakaev M, Oksuz BA, Venev SV, Xiao Y. Cooltools: Enabling high-resolution Hi-C analysis in Python

Chromosome conformation capture (3C) technologies reveal the incredible complexity of genome organization. Maps of increasing size, depth, and resolution are now used to probe genome architecture across cell states, types, and organisms. Larger datasets add challenges at each step of computational analysis, from storage and memory constraints to researchers' time; however, analysis tools that meet these increased resource demands have not kept pace. Furthermore, existing tools offer limited support for customizing analysis for specific use cases or new biology. Here we introduce cooltools (https://github.com/open2c/cooltools), a suite of computational tools that enables flexible, scalable, and reproducible analysis of high-resolution contact frequency data. Cooltools leverages the widely-adopted cooler format which handles storage and access for high-resolution datasets. Cooltools provides a paired command line interface (CLI) and Python application programming interface (API), which respectively facilitate workflows on high-performance computing clusters and in interactive analysis environments. In short, cooltools enables the effective use of the latest and largest genome folding datasets.

Tet-mediated DNA methylation dynamics affect chromosome organization

DNA Methylation is a significant epigenetic modification that can modulate chromosome states, but its role in orchestrating chromosome organization has not been well elucidated. Here we systematically assessed the effects of DNA Methylation on chromosome organization with a multi-omics strategy to capture DNA Methylation and high-order chromosome interaction simultaneously on mouse embryonic stem cells with DNA methylation dioxygenase Tet triple knock-out (Tet-TKO). Globally, upon Tet-TKO, we observed weakened compartmentalization, corresponding to decreased methylation differences between CpG island (CGI) rich and poor domains. Tet-TKO could also induce hypermethylation for the CTCF binding peaks in TAD boundaries and chromatin loop anchors. Accordingly, CTCF peak generally weakened upon Tet-TKO, which results in weakened TAD structure and depletion of long-range chromatin loops. Genes that lost enhancer-promoter looping upon Tet-TKO showed DNA hypermethylation in their gene bodies, which may compensate for the disruption of gene expression. We also observed distinct effects of Tet1 and Tet2 on chromatin organization and increased DNA methylation correlation on spatially interacted fragments upon Tet inactivation. Our work showed the broad effects of Tet inactivation and DNA methylation dynamics on chromosome organization.

The PTM profiling of CTCF reveals the regulation of 3D chromatin structure by O-GlcNAcylation

CCCTC-binding factor (CTCF), a ubiquitously expressed and highly conserved protein, is known to play a critical role in chromatin structure. Post-translational modifications (PTMs) diversify the functions of protein to regulate numerous cellular processes. However, the effects of PTMs on the genome-wide binding of CTCF and the organization of three-dimensional (3D) chromatin structure have not been fully understood. In this study, we uncovered the PTM profiling of CTCF and demonstrated that CTCF can be O-GlcNAcylated and arginine methylated. Functionally, we demonstrated that O-GlcNAcylation inhibits CTCF binding to chromatin. Meanwhile, deficiency of CTCF O-GlcNAcylation results in the disruption of loop domains and the alteration of chromatin loops associated with cellular development. Furthermore, the deficiency of CTCF O-GlcNAcylation increases the expression of developmental genes and negatively regulates maintenance and establishment of stem cell pluripotency. In conclusion, these results provide key insights into the role of PTMs for the 3D chromatin structure.

Chromatin organization and behavior in HRAS-transformed mouse fibroblasts

In higher eukaryotic cells, a string of nucleosomes, where long genomic DNA is wrapped around core histones, are rather irregularly folded into a number of condensed chromatin domains, which have been revealed by super-resolution imaging and Hi-C technologies. Inside these domains, nucleosomes fluctuate and locally behave like a liquid. The behavior of chromatin may be highly related to DNA transaction activities such as transcription and repair, which are often upregulated in cancer cells. To investigate chromatin behavior in cancer cells and compare those of cancer and non-cancer cells, we focused on oncogenic-HRAS (Gly12Val)-transformed mouse fibroblasts CIRAS-3 cells and their parental 10T1/2 cells. CIRAS-3 cells are tumorigenic and highly metastatic. First, we found that HRAS-induced transformation altered not only chromosome structure, but also nuclear morphology in the cell. Using single-nucleosome imaging/tracking in live cells, we demonstrated that nucleosomes are locally more constrained in CIRAS-3 cells than in 10T1/2 cells. Consistently, heterochromatin marked with H3K27me3 was upregulated in CIRAS-3 cells. Finally, Hi-C analysis showed enriched interactions of the B-B compartment in CIRAS-3 cells, which likely represents transcriptionally inactive chromatin. Increased heterochromatin may play an important role in cell migration, as they have been reported to increase during metastasis. Our study also suggests that single-nucleosome imaging provides new insights into how local chromatin is structured in living cells.

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.

Hi-C profiling in tissues reveals 3D chromatin-regulated breast tumor heterogeneity and tumor-specific looping-mediated biological pathways

Current knowledge in three-dimensional (3D) chromatin regulation in normal and disease states was mostly accumulated through Hi-C profiling in in vitro cell culture system. The limitations include failing to recapitulate disease-specific physiological properties and often lacking clinically relevant disease microenvironment. In this study, we conduct tissue-specific Hi-C profiling in a pilot cohort of 12 breast tissues comprising of two normal tissues (NTs) and ten ER+ breast tumor tissues (TTs) including five primary tumors (PTs), and five tamoxifen-treated recurrent tumors (RTs). We find largely preserved compartments, highly heterogeneous topological associated domains (TADs) and intensively variable chromatin loops among breast tumors, demonstrating 3D chromatin-regulated breast tumor heterogeneity. Further cross-examination identifies RT-specific looping-mediated biological pathways and suggests CA2, an enhancer-promoter looping (EPL)-mediated target gene within the bicarbonate transport metabolism pathway, might play a role in driving the tamoxifen resistance. Remarkably, the inhibition of CA2 not only impedes tumor growth both in vitro and in vivo , but also reverses chromatin looping. Our study thus yields significant mechanistic insights into the role and clinical relevance of 3D chromatin architecture in breast cancer endocrine resistance.

The potential of epigenetic therapy to target the 3D epigenome in endocrine-resistant breast cancer

Three-dimensional (3D) epigenome remodeling is an important mechanism of gene deregulation in cancer. However, its potential as a target to counteract therapy resistance remains largely unaddressed. Here, we show that epigenetic therapy with decitabine (5-Aza-mC) suppresses tumor growth in xenograft models of pre-clinical metastatic estrogen receptor positive (ER+) breast tumor. Decitabine-induced genome-wide DNA hypomethylation results in large-scale 3D epigenome deregulation, including de-compaction of higher-order chromatin structure and loss of boundary insulation of topologically associated domains. Significant DNA hypomethylation associates with ectopic activation of ER-enhancers, gain in ER binding, creation of new 3D enhancer-promoter interactions and concordant up-regulation of ER-mediated transcription pathways. Importantly, long-term withdrawal of epigenetic therapy partially restores methylation at ER-enhancer elements, resulting in a loss of ectopic 3D enhancer-promoter interactions and associated gene repression. Our study illustrates the potential of epigenetic therapy to target ER+ endocrine-resistant breast cancer by DNA methylation-dependent rewiring of 3D chromatin interactions, which are associated with the suppression of tumor growth.

SATB2 organizes the 3D genome architecture of cognition in cortical neurons

The DNA-binding protein SATB2 is genetically linked to human intelligence. We studied its influence on the three-dimensional (3D) epigenome by mapping chromatin interactions and accessibility in control versus SATB2-deficient cortical neurons. We find that SATB2 affects the chromatin looping between enhancers and promoters of neuronal-activity-regulated genes, thus influencing their expression. It also alters A/B compartments, topologically associating domains, and frequently interacting regions. Genes linked to SATB2-dependent 3D genome changes are implicated in highly specialized neuronal functions and contribute to cognitive ability and risk for neuropsychiatric and neurodevelopmental disorders. Non-coding DNA regions with a SATB2-dependent structure are enriched for common variants associated with educational attainment, intelligence, and schizophrenia. Our data establish SATB2 as a cell-type-specific 3D genome modulator, which operates both independently and in cooperation with CCCTC-binding factor (CTCF) to set up the chromatin landscape of pyramidal neurons for cognitive processes.

Orchestrating chromosome conformation capture analysis with Bioconductor

Genome-wide chromatin conformation capture assays provide formidable insights into the spatial organization of genomes. However, due to the complexity of the data structure, their integration in multi-omics workflows remains challenging. We present data structures, computational methods and visualization tools available in Bioconductor to investigate Hi-C, micro-C and other 3C-related data, in R. An online book ( https://bioconductor.org/books/OHCA/ ) further provides prospective end users with a number of workflows to process, import, analyze and visualize any type of chromosome conformation capture data.

TDP-43 chronic deficiency leads to dysregulation of transposable elements and gene expression by affecting R-loop and 5hmC crosstalk

TDP-43 is an RNA/DNA-binding protein that forms aggregates in various brain disorders. TDP-43 engages in many aspects of RNA metabolism, but its molecular roles in regulating genes and transposable elements (TEs) have not been extensively explored. Chronic TDP-43 knockdown impairs cell proliferation and cellular responses to DNA damage. At the molecular level, TDP-43 chronic deficiency affects gene expression either locally or distally by concomitantly altering the crosstalk between R-loops and 5-hydroxymethylcytosine (5hmC) in gene bodies and long-range enhancer/promoter interactions. Furthermore, TDP-43 knockdown induces substantial disease-relevant TE activation by influencing their R-loop and 5hmC homeostasis in a locus-specific manner. Together, our findings highlight the genomic roles of TDP-43 in modulating R-loop-5hmC coordination in coding genes, distal regulatory elements, and TEs, presenting a general and broad molecular mechanism underlying the contributions of proteinopathies to the etiology of neurodegenerative disorders.

An Atlas of Promoter Chromatin Modifications and HiChIP Regulatory Interactions in Human Subcutaneous Adipose-Derived Stem Cells

The genome of human adipose-derived stem cells (ADSCs) from abdominal and gluteofemoral adipose tissue depots are maintained in depot-specific stable epigenetic conformations that influence cell-autonomous gene expression patterns and drive unique depot-specific functions. The traditional approach to explore tissue-specific transcriptional regulation has been to correlate differential gene expression to the nearest-neighbor linear-distance regulatory region defined by associated chromatin features including open chromatin status, histone modifications, and DNA methylation. This has provided important information; nonetheless, the approach is limited because of the known organization of eukaryotic chromatin into a topologically constrained three-dimensional network. This network positions distal regulatory elements in spatial proximity with gene promoters which are not predictable based on linear genomic distance. In this work, we capture long-range chromatin interactions using HiChIP to identify remote genomic regions that influence the differential regulation of depot-specific genes in ADSCs isolated from different adipose depots. By integrating these data with RNA-seq results and histone modifications identified by ChIP-seq, we uncovered distal regulatory elements that influence depot-specific gene expression in ADSCs. Interestingly, a subset of the HiChIP-defined chromatin loops also provide previously unknown connections between waist-to-hip ratio GWAS variants with genes that are known to significantly influence ADSC differentiation and adipocyte function.

Stage-specific dynamic reorganization of genome topology shapes transcriptional neighborhoods in developing human retinal organoids

We have generated a high-resolution Hi-C map of developing human retinal organoids to elucidate spatiotemporal dynamics of genomic architecture and its relationship with gene expression patterns. We demonstrate progressive stage-specific alterations in DNA topology and correlate these changes with transcription of cell-type-restricted gene markers during retinal differentiation. Temporal Hi-C reveals a shift toward A compartment for protein-coding genes and B compartment for non-coding RNAs, displaying high and low expression, respectively. Notably, retina-enriched genes are clustered near lost boundaries of topologically associated domains (TADs), and higher-order assemblages (i.e., TAD cliques) localize in active chromatin regions with binding sites for eye-field transcription factors. These genes gain chromatin contacts at their transcription start site as organoid differentiation proceeds. Our study provides a global view of chromatin architecture dynamics associated with diversification of cell types during retinal development and serves as a foundational resource for in-depth functional investigations of retinal developmental traits.

3D genome remodeling and homologous pairing during meiotic prophase of mouse oogenesis and spermatogenesis

During meiosis, the chromatin and transcriptome undergo prominent switches. Although recent studies have explored the genome reorganization during spermatogenesis, the chromatin remodeling in oogenesis and characteristics of homologous pairing remain largely elusive. We comprehensively compared chromatin structures and transcriptomes at successive substages of meiotic prophase in both female and male mice using low-input high-through chromosome conformation capture (Hi-C) and RNA sequencing (RNA-seq). Compartments and topologically associating domains (TADs) gradually disappeared and slowly recovered in both sexes. We found that homologs adopted different sex-conserved pairing strategies prior to and after the leptotene-to-zygotene transition, changing from long interspersed nuclear element (LINE)-enriched compartments B to short interspersed nuclear element (SINE)-enriched compartments A. We complemented marker genes and predicted the sex-specific meiotic sterile genes for each substage. This study provides valuable insights into the similarities and distinctions between sexes in chromosome architecture, homologous pairing, and transcriptome during meiotic prophase of both oogenesis and spermatogenesis.

Loss of cohesin regulator PDS5A reveals repressive role of Polycomb loops

Polycomb Repressive Complexes 1 and 2 (PRC1, PRC2) are conserved epigenetic regulators that promote transcriptional gene silencing. PRC1 and PRC2 converge on shared targets, catalyzing repressive histone modifications. Additionally, a subset of PRC1/PRC2 targets engage in long-range interactions whose functions in gene silencing are poorly understood. Using a CRISPR screen in mouse embryonic stem cells, we found that the cohesin regulator PDS5A links transcriptional silencing by Polycomb and 3D genome organization. PDS5A deletion impairs cohesin unloading and results in derepression of a subset of endogenous PRC1/PRC2 target genes. Importantly, derepression is not linked to loss of Polycomb chromatin domains. Instead, PDS5A removal causes aberrant cohesin activity leading to ectopic insulation sites, which disrupt the formation of ultra-long Polycomb loops. We show that these loops are important for robust silencing at a subset of PRC1/PRC2 target genes and that maintenance of cohesin-dependent genome architecture is critical for Polycomb regulation.

Noncoding function of super enhancer derived mRNA in modulating neighboring gene expression and TAD interaction

Super enhancers are important regulators of gene expression that often overlap with protein-coding genes. However, it is unclear whether the overlapping protein-coding genes and the mRNA derived from them contribute to enhancer activity. Using an erythroid-specific super enhancer that overlaps the Cpox gene as a model, we found that Cpox mRNA has a non-coding function in regulating neighboring protein-coding genes, eRNA expression and TAD interactions. Depletion of Cpox mRNA leads to accumulation of H3K27me3 and release of p300 from the Cpox locus, activating an intra-TAD enhancer and gene expression. Additionally, we identified a head-to-tail interaction between the TAD boundary genes Cpox and Dcbld2 that is facilitated by a novel type of repressive loop anchored by p300 and PRC2/H3K27me3. Our results uncover a regulatory role for mRNA transcribed within a super enhancer context and provide insight into head-to-tail inter-gene interaction in the regulation of gene expression and oncogene activation.

Lineage-specific 3D genome organization is assembled at multiple scales by IKAROS

A generic level of chromatin organization generated by the interplay between cohesin and CTCF suffices to limit promiscuous interactions between regulatory elements, but a lineage-specific chromatin assembly that supersedes these constraints is required to configure the genome to guide gene expression changes that drive faithful lineage progression. Loss-of-function approaches in B cell precursors show that IKAROS assembles interactions across megabase distances in preparation for lymphoid development. Interactions emanating from IKAROS-bound enhancers override CTCF-imposed boundaries to assemble lineage-specific regulatory units built on a backbone of smaller invariant topological domains. Gain of function in epithelial cells confirms IKAROS' ability to reconfigure chromatin architecture at multiple scales. Although the compaction of the Igκ locus required for genome editing represents a function of IKAROS unique to lymphocytes, the more general function to preconfigure the genome to support lineage-specific gene expression and suppress activation of extra-lineage genes provides a paradigm for lineage restriction.

Breast cancer risk SNPs converge on estrogen receptor binding sites commonly shared between breast tumors to locally alter estrogen signalling output

Estrogen Receptor alpha (ERα) is the main driver and prime drug target in luminal breast. ERα chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ERα chromatin action, along with its biological implications. Here, we use a large set of ERα ChIP-seq data from 70 ERα+ breast cancers to explore inter-patient heterogeneity in ERα DNA binding, to reveal a striking inter-tumor heterogeneity of ERα action. Interestingly, commonly-shared ERα sites showed the highest estrogen-driven enhancer activity and were most-engaged in long-range chromatin interactions. In addition, the most-commonly shared ERα-occupied enhancers were enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ERα and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we could confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ERα-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ERα landscape, with the most-common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.

Super-enhancer hijacking drives ectopic expression of hedgehog pathway ligands in meningiomas

Hedgehog signaling mediates embryologic development of the central nervous system and other tissues and is frequently hijacked by neoplasia to facilitate uncontrolled cellular proliferation. Meningiomas, the most common primary brain tumor, exhibit Hedgehog signaling activation in 6.5% of cases, triggered by recurrent mutations in pathway mediators such as SMO. In this study, we find 35.6% of meningiomas that lack previously known drivers acquired various types of somatic structural variations affecting chromosomes 2q35 and 7q36.3. These cases exhibit ectopic expression of Hedgehog ligands, IHH and SHH, respectively, resulting in Hedgehog signaling activation. Recurrent tandem duplications involving IHH permit de novo chromatin interactions between super-enhancers within DIRC3 and a locus containing IHH. Our work expands the landscape of meningioma molecular drivers and demonstrates enhancer hijacking of Hedgehog ligands as a route to activate this pathway  in neoplasia.

Three-way contact analysis characterizes the higher order organization of the Tcra locus

The generation of highly diverse antigen receptors in T and B lymphocytes relies on V(D)J recombination. The enhancer Eα has been implicated in regulating the accessibility of Vα and Jα genes through long-range interactions during rearrangements of the T-cell antigen receptor gene Tcra. However, direct evidence for Eα physically mediating the interaction of Vα and Jα genes is still lacking. In this study, we utilized the 3C-HTGTS assay, a chromatin interaction technique based on 3C, to analyze the higher order chromatin structure of the Tcra locus. Our analysis revealed the presence of sufficient information in the 3C-HTGTS data to detect multiway contacts. Three-way contact analysis of the Tcra locus demonstrated the co-occurrence of the proximal Jα genes, Vα genes and Eα in CD4+CD8+ double-positive thymocytes. Notably, the INT2-TEAp loop emerged as a prominent structure likely to be responsible for bringing the proximal Jα genes and the Vα genes into proximity. Moreover, the enhancer Eα utilizes this loop to establish physical proximity with the proximal Vα gene region. This study provides insights into the higher order chromatin structure of the Tcra locus, shedding light on the spatial organization of chromatin and its impact on V(D)J recombination.

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.

Simultaneous profiling of chromatin architecture and transcription in single cells

The three-dimensional structure of chromatin plays a crucial role in development and disease, both of which are associated with transcriptional changes. However, given the heterogeneity in single-cell chromatin architecture and transcription, the regulatory relationship between the three-dimensional chromatin structure and gene expression is difficult to explain based on bulk cell populations. Here we develop a single-cell, multimodal, omics method allowing the simultaneous detection of chromatin architecture and messenger RNA expression by sequencing (single-cell transcriptome sequencing (scCARE-seq)). Applying scCARE-seq to examine chromatin architecture and transcription from 2i to serum single mouse embryonic stem cells, we observe improved separation of cell clusters compared with single-cell chromatin conformation capture. In addition, after defining the cell-cycle phase of each cell through chromatin architecture extracted by scCARE-seq, we find that periodic changes in chromatin architecture occur in parallel with transcription during the cell cycle. These findings highlight the potential of scCARE-seq to facilitate comprehensive analyses that may boost our understanding of chromatin architecture and transcription in the same single cell.

Replication dynamics identifies the folding principles of the inactive X chromosome

Chromosome-wide late replication is an enigmatic hallmark of the inactive X chromosome (Xi). How it is established and what it represents remains obscure. By single-cell DNA replication sequencing, here we show that the entire Xi is reorganized to replicate rapidly and uniformly in late S-phase during X-chromosome inactivation (XCI), reflecting its relatively uniform structure revealed by 4C-seq. Despite this uniformity, only a subset of the Xi became earlier replicating in SmcHD1-mutant cells. In the mutant, these domains protruded out of the Xi core, contacted each other and became transcriptionally reactivated. 4C-seq suggested that they constituted the outermost layer of the Xi even before XCI and were rich in escape genes. We propose that this default positioning forms the basis for their inherent heterochromatin instability in cells lacking the Xi-binding protein SmcHD1 or exhibiting XCI escape. These observations underscore the importance of 3D genome organization for heterochromatin stability and gene regulation.