Spatially coordinated heterochromatinization of long synaptic genes in fragile X syndrome.

Short tandem repeat (STR) instability causes transcriptional silencing in several repeat expansion disorders. In fragile X syndrome (FXS), mutation-length expansion of a CGG STR represses FMR1 via local DNA methylation. Here, we find megabase-scale H3K9me3 domains on autosomes and encompassing FMR1 on the X chromosome in FXS patient-derived iPSCs, iPSC-derived neural progenitors, EBV-transformed lymphoblasts, and brain tissue with mutation-length CGG expansion. H3K9me3 domains connect via inter-chromosomal interactions and demarcate severe misfolding of TADs and loops. They harbor long synaptic genes replicating at the end of S phase, replication-stress-induced double-strand breaks, and STRs prone to stepwise somatic instability. CRISPR engineering of the mutation-length CGG to premutation length reverses H3K9me3 on the X chromosome and multiple autosomes, refolds TADs, and restores gene expression. H3K9me3 domains can also arise in normal-length iPSCs created with perturbations linked to genome instability, suggesting their relevance beyond FXS. Our results reveal Mb-scale heterochromatinization and trans interactions among loci susceptible to instability.

Multi-feature clustering of CTCF binding creates robustness for loop extrusion blocking and Topologically Associating Domain boundaries.

Topologically Associating Domains (TADs) separate vertebrate genomes into insulated regulatory neighborhoods that focus genome-associated processes. TADs are formed by Cohesin-mediated loop extrusion, with many TAD boundaries consisting of clustered binding sites of the CTCF insulator protein. Here we determine how this clustering of CTCF binding contributes to the blocking of loop extrusion and the insulation between TADs. We identify enrichment of three features of CTCF binding at strong TAD boundaries, consisting of strongly bound and closely spaced CTCF binding peaks, with a further enrichment of DNA-binding motifs within these peaks. Using multi-contact Nano-C analysis in cells with normal and perturbed CTCF binding, we establish that individual CTCF binding sites contribute to the blocking of loop extrusion, but in an incomplete manner. When clustered, individual CTCF binding sites thus create a stepwise insulation between neighboring TADs. Based on these results, we propose a model whereby multiple instances of temporal loop extrusion blocking create strong insulation between TADs.

High-throughput Oligopaint screen identifies druggable 3D genome regulators.

The human genome functions as a three-dimensional chromatin polymer, driven by a complex collection of chromosome interactions1-3. Although the molecular rules governing these interactions are being quickly elucidated, relatively few proteins regulating this process have been identified. Here, to address this gap, we developed high-throughput DNA or RNA labelling with optimized Oligopaints (HiDRO)-an automated imaging pipeline that enables the quantitative measurement of chromatin interactions in single cells across thousands of samples. By screening the human druggable genome, we identified more than 300 factors that influence genome folding during interphase. Among these, 43 genes were validated as either increasing or decreasing interactions between topologically associating domains. Our findings show that genetic or chemical inhibition of the ubiquitous kinase GSK3A leads to increased long-range chromatin looping interactions in a genome-wide and cohesin-dependent manner. These results demonstrate the importance of GSK3A signalling in nuclear architecture and the use of HiDRO for identifying mechanisms of spatial genome organization.

Co-depletion of NIPBL and WAPL balance cohesin activity to correct gene misexpression.

The relationship between cohesin-mediated chromatin looping and gene expression remains unclear. NIPBL and WAPL are two opposing regulators of cohesin activity; depletion of either is associated with changes in both chromatin folding and transcription across a wide range of cell types. However, a direct comparison of their individual and combined effects on gene expression in the same cell type is lacking. We find that NIPBL or WAPL depletion in human HCT116 cells each alter the expression of ~2,000 genes, with only ~30% of the genes shared between the conditions. We find that clusters of differentially expressed genes within the same topologically associated domain (TAD) show coordinated misexpression, suggesting some genomic domains are especially sensitive to both more or less cohesin. Finally, co-depletion of NIPBL and WAPL restores the majority of gene misexpression as compared to either knockdown alone. A similar set of NIPBL-sensitive genes are rescued following CTCF co-depletion. Together, this indicates that altered transcription due to reduced cohesin activity can be functionally offset by removal of either its negative regulator (WAPL) or the physical barriers (CTCF) that restrict loop-extrusion events.

CTCF blocks antisense transcription initiation at divergent promoters.

Transcription at most promoters is divergent, initiating at closely spaced oppositely oriented core promoters to produce sense transcripts along with often unstable upstream antisense transcripts (uasTrx). How antisense transcription is regulated and to what extent it is coordinated with sense transcription is not well understood. Here, by combining acute degradation of the multi-functional transcription factor CTCF and nascent transcription measurements, we find that CTCF specifically suppresses antisense but not sense transcription at hundreds of divergent promoters. Primary transcript RNA-FISH shows that CTCF lowers burst fraction but not burst intensity of uasTrx and that co-bursting of sense and antisense transcripts is disfavored. Genome editing, chromatin conformation studies and high-resolution transcript mapping revealed that precisely positioned CTCF directly suppresses the initiation of uasTrx, in a manner independent of its architectural function. In sum, CTCF shapes the transcriptional landscape in part by suppressing upstream antisense transcription.

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.

Cohesin promotes stochastic domain intermingling to ensure proper regulation of boundary-proximal genes.

The human genome can be segmented into topologically associating domains (TADs), which have been proposed to spatially sequester genes and regulatory elements through chromatin looping. Interactions between TADs have also been suggested, presumably because of variable boundary positions across individual cells. However, the nature, extent and consequence of these dynamic boundaries remain unclear. Here, we combine high-resolution imaging with Oligopaint technology to quantify the interaction frequencies across both weak and strong boundaries. We find that chromatin intermingling across population-defined boundaries is widespread but that the extent of permissibility is locus-specific. Cohesin depletion, which abolishes domain formation at the population level, does not induce ectopic interactions but instead reduces interactions across all boundaries tested. In contrast, WAPL or CTCF depletion increases inter-domain contacts in a cohesin-dependent manner. Reduced chromatin intermingling due to cohesin loss affects the topology and transcriptional bursting frequencies of genes near boundaries. We propose that cohesin occasionally bypasses boundaries to promote incorporation of boundary-proximal genes into neighboring domains.

Single cell analysis pushes the boundaries of TAD formation and function.

Eukaryotic genomes encode genetic information in their linear sequence, but appropriate expression of their genes requires chromosomes to fold into complex three-dimensional structures. Fueled by a growing collection of sequencing and imaging-based technologies, studies have uncovered a hierarchy of DNA interactions, from small chromatin loops that connect genes and enhancers to larger topologically associated domains (TADs) and compartments. However, despite the remarkable conservation of these organizational features, we have a very limited understanding of how this organization influences gene expression. This issue is further complicated in the context of single-cell heterogeneity, as has recently been revealed at both the level of gene activation and chromatin topology. Here, we provide a perspective on recent studies that address cell-to-cell variability and the relationship between structural heterogeneity and gene expression. We propose that transcription is regulated by variable 3D structures driven by at least two independent and partially redundant mechanisms. Collectively, this may provide flexibility to transcriptional regulation at the level of individual cells as well as reproducibility across whole tissues.

ARID1A spatially partitions interphase chromosomes.

ARID1A, a subunit of the SWItch/Sucrose Non-Fermentable (SWI/SNF) chromatin-remodeling complex, localizes to both promoters and enhancers to influence transcription. However, the role of ARID1A in higher-order spatial chromosome partitioning and genome organization is unknown. Here, we show that ARID1A spatially partitions interphase chromosomes and regulates higher-order genome organization. The SWI/SNF complex interacts with condensin II, and they display significant colocalizations at enhancers. ARID1A knockout drives the redistribution of condensin II preferentially at enhancers, which positively correlates with changes in transcription. ARID1A and condensin II contribute to transcriptionally inactive B-compartment formation, while ARID1A weakens the border strength of topologically associated domains. Condensin II redistribution induced by ARID1A knockout positively correlates with chromosome sizes, which negatively correlates with interchromosomal interactions. ARID1A loss increases the trans interactions of small chromosomes, which was validated by three-dimensional interphase chromosome painting. These results demonstrate that ARID1A is important for large-scale genome folding and spatially partitions interphase chromosomes.

Hotspots of aberrant enhancer activity punctuate the colorectal cancer epigenome.

In addition to mutations in genes, aberrant enhancer element activity at non-coding regions of the genome is a key driver of tumorigenesis. Here, we perform epigenomic enhancer profiling of a cohort of more than forty genetically diverse human colorectal cancer (CRC) specimens. Using normal colonic crypt epithelium as a comparator, we identify enhancers with recurrently gained or lost activity across CRC specimens. Of the enhancers highly recurrently activated in CRC, most are constituents of super enhancers, are occupied by AP-1 and cohesin complex members, and originate from primed chromatin. Many activate known oncogenes, and CRC growth can be mitigated through pharmacologic inhibition or genome editing of these loci. Nearly half of all GWAS CRC risk loci co-localize to recurrently activated enhancers. These findings indicate that the CRC epigenome is defined by highly recurrent epigenetic alterations at enhancers which activate a common, aberrant transcriptional programme critical for CRC growth and survival.

Modeling disease risk through analysis of physical interactions between genetic variants within chromatin regulatory circuitry.

SNPs associated with disease susceptibility often reside in enhancer clusters, or super-enhancers. Constituents of these enhancer clusters cooperate to regulate target genes and often extend beyond the linkage disequilibrium (LD) blocks containing risk SNPs identified in genome-wide association studies (GWAS). We identified 'outside variants', defined as SNPs in weak LD with GWAS risk SNPs that physically interact with risk SNPs as part of a target gene's regulatory circuitry. These outside variants further explain variation in target gene expression beyond that explained by GWAS-associated SNPs. Additionally, the clinical risk associated with GWAS SNPs is considerably modified by the genotype of outside variants. Collectively, these findings suggest a potential model in which outside variants and GWAS SNPs that physically interact in 3D chromatin collude to influence target transcript levels as well as clinical risk. This model offers an additional hypothesis for the source of missing heritability for complex traits.