3DeFDR: statistical methods for identifying cell type-specific looping interactions in 5C and Hi-C data

Cell-type-specific loops linked to RNA polymerase II elongation in human neural differentiation

DNA is folded into higher-order structures that shape and are shaped by genome function. The role of long-range loops in the establishment of new gene expression patterns during cell fate transitions remains poorly understood. Here, we investigate the link between cell-specific loops and RNA polymerase II (RNA Pol II) during neural lineage commitment. We find thousands of loops decommissioned or gained de novo upon differentiation of human induced pluripotent stem cells (hiPSCs) to neural progenitor cells (NPCs) and post-mitotic neurons. During hiPSC-to-NPC and NPC-to-neuron transitions, genes changing from RNA Pol II initiation to elongation are >4-fold more likely to anchor cell-specific loops than repressed genes. Elongated genes exhibit significant mRNA upregulation when connected in cell-specific promoter-enhancer loops but not invariant promoter-enhancer loops or promoter-promoter loops or when unlooped. Genes transitioning from repression to RNA Pol II initiation exhibit a slight mRNA increase independent of loop status. Our data link cell-specific loops and robust RNA Pol II-mediated elongation during neural cell fate transitions.

A multi-looping chromatin signature predicts dysregulated gene expression in neurons with familial Alzheimer's disease mutations

Mammalian genomes fold into tens of thousands of long-range loops, but their functional role and physiologic relevance remain poorly understood. Here, using human post-mitotic neurons with rare familial Alzheimer's disease (FAD) mutations, we identify hundreds of reproducibly dysregulated genes and thousands of miswired loops prior to amyloid accumulation and tau phosphorylation. Single loops do not predict expression changes; however, the severity and direction of change in mRNA levels and single-cell burst frequency strongly correlate with the number of FAD-gained or -lost promoter-enhancer loops. Classic architectural proteins CTCF and cohesin do not change occupancy in FAD-mutant neurons. Instead, we unexpectedly find TAATTA motifs amenable to binding by DLX homeodomain transcription factors and changing noncoding RNAPolII signal at FAD-dynamic promoter-enhancer loops. DLX1/5/6 mRNA levels are strongly upregulated in FAD-mutant neurons coincident with a shift in excitatory-to-inhibitory gene expression and miswiring of multi-loops connecting enhancers to neural subtype genes. DLX1 overexpression is sufficient for loop miswiring in wildtype neurons, including lost and gained loops at enhancers with tandem TAATTA arrays and singular TAATTA motifs, respectively. Our data uncover a genome structure-function relationship between multi-loop miswiring and dysregulated excitatory and inhibitory transcriptional programs during lineage commitment of human neurons homozygously-engineered with rare FAD mutations.

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.

Cell type-specific loops linked to RNA polymerase II elongation in human neural differentiation

DNA is folded into higher-order structures that shape and are shaped by genome function. The role for long-range loops in the establishment of new gene expression patterns during cell fate transitions remains poorly understood. Here, we investigate the link between cell-specific loops and RNA polymerase II (RNAPolII) during neural lineage commitment. We find thousands of loops decommissioned or gained de novo upon differentiation of human induced pluripotent stem cells (hiPSCs) to neural progenitors (NPCs) and post-mitotic neurons. During hiPSC-to-NPC and NPC-to-neuron transitions, genes changing from RNAPolII initiation to elongation are >4-fold more likely to anchor cell-specific loops than repressed genes. Elongated genes exhibit significant mRNA upregulation when connected in cell-specific promoter-enhancer loops but not invariant promoter-enhancer loops, promoter-promoter loops, or unlooped. Genes transitioning from repression to RNAPolII initiation exhibit slight mRNA increase independent of loop status. Our data link cell-specific loops and robust RNAPolII-mediated elongation during neural cell fate transitions.

A deep learning method for replicate-based analysis of chromosome conformation contacts using Siamese neural networks

The organisation of the genome in nuclear space is an important frontier of biology. Chromosome conformation capture methods such as Hi-C and Micro-C produce genome-wide chromatin contact maps that provide rich data containing quantitative and qualitative information about genome architecture. Most conventional approaches to genome-wide chromosome conformation capture data are limited to the analysis of pre-defined features, and may therefore miss important biological information. One constraint is that biologically important features can be masked by high levels of technical noise in the data. Here we introduce a replicate-based method for deep learning from chromatin conformation contact maps. Using a Siamese network configuration our approach learns to distinguish technical noise from biological variation and outperforms image similarity metrics across a range of biological systems. The features extracted from Hi-C maps after perturbation of cohesin and CTCF reflect the distinct biological functions of cohesin and CTCF in the formation of domains and boundaries, respectively. The learnt distance metrics are biologically meaningful, as they mirror the density of cohesin and CTCF binding. These properties make our method a powerful tool for the exploration of chromosome conformation capture data, such as Hi-C capture Hi-C, and Micro-C.

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.

Cohesin-mediated loop anchors confine the locations of human replication origins

DNA replication occurs through an intricately regulated series of molecular events and is fundamental for genome stability1,2. At present, it is unknown how the locations of replication origins are determined in the human genome. Here we dissect the role of topologically associating domains (TADs)3-6, subTADs7 and loops8 in the positioning of replication initiation zones (IZs). We stratify TADs and subTADs by the presence of corner-dots indicative of loops and the orientation of CTCF motifs. We find that high-efficiency, early replicating IZs localize to boundaries between adjacent corner-dot TADs anchored by high-density arrays of divergently and convergently oriented CTCF motifs. By contrast, low-efficiency IZs localize to weaker dotless boundaries. Following ablation of cohesin-mediated loop extrusion during G1, high-efficiency IZs become diffuse and delocalized at boundaries with complex CTCF motif orientations. Moreover, G1 knockdown of the cohesin unloading factor WAPL results in gained long-range loops and narrowed localization of IZs at the same boundaries. Finally, targeted deletion or insertion of specific boundaries causes local replication timing shifts consistent with IZ loss or gain, respectively. Our data support a model in which cohesin-mediated loop extrusion and stalling at a subset of genetically encoded TAD and subTAD boundaries is an essential determinant of the locations of replication origins in human S phase.

Cohesin-dependence of neuronal gene expression relates to chromatin loop length

Cohesin and CTCF are major drivers of 3D genome organization, but their role in neurons is still emerging. Here, we show a prominent role for cohesin in the expression of genes that facilitate neuronal maturation and homeostasis. Unexpectedly, we observed two major classes of activity-regulated genes with distinct reliance on cohesin in mouse primary cortical neurons. Immediate early genes (IEGs) remained fully inducible by KCl and BDNF, and short-range enhancer-promoter contacts at the IEGs Fos formed robustly in the absence of cohesin. In contrast, cohesin was required for full expression of a subset of secondary response genes characterized by long-range chromatin contacts. Cohesin-dependence of constitutive neuronal genes with key functions in synaptic transmission and neurotransmitter signaling also scaled with chromatin loop length. Our data demonstrate that key genes required for the maturation and activation of primary cortical neurons depend on cohesin for their full expression, and that the degree to which these genes rely on cohesin scales with the genomic distance traversed by their chromatin contacts.

Assessing Specific Networks of Chromatin Interactions with HiChIP

The introduction of chromosome conformation capture (3C)-based technologies coupled with next-generation sequencing have significantly advanced our understanding of how the genetic material is organized within the eukaryotic nucleus. Three-dimensional (3D) genomic organization occurs at hierarchical levels, ranging from chromosome territories and subnuclear compartments to smaller self-associated domains and fine-scale chromatin interactions. The latter can be further categorized into different subtypes, such as structural or regulatory, based either on their presumed functionality and/or the factors that mediate their formation. Various enrichment strategies coupled with 3C-based technologies have been developed to prospectively isolate and quantify chromatin interactions around regions occupied by specific proteins or marks of interest. These approaches not only enable high-resolution characterization of the selected chromatin contacts at a cost-effective manner, but also offer important biological insights into their organizational principles and regulatory function. In this chapter, we will focus on the recently developed HiChIP technology with an emphasis on the discovery of putative active enhancers and promoter interactions in cell types of interest. We will describe the specific steps for designing, performing and analyzing successful HiChIP experiments as well as important limitations and considerations.