Formation of Chromosomal Domains by Loop Extrusion

Machine and Deep Learning Methods for Predicting 3D Genome Organization

Three-dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, topologically associating domains (TADs), and A/B compartments, play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers and transcription factor binding site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, and TAD boundaries) and analyze their pros and cons. We also point out obstacles to the computational prediction of 3D interactions and suggest future research directions.

Long non-coding RNAs: roles in cellular stress responses and epigenetic mechanisms regulating chromatin

Most of the genome is transcribed into RNA but only 2% of the sequence codes for proteins. Non-coding RNA transcripts include a very large number of long noncoding RNAs (lncRNAs). A growing number of identified lncRNAs operate in cellular stress responses, for example in response to hypoxia, genotoxic stress, and oxidative stress. Additionally, lncRNA plays important roles in epigenetic mechanisms operating at chromatin and in maintaining chromatin architecture. Here, we address three lncRNA topics that have had significant recent advances. The first is an emerging role for many lncRNAs in cellular stress responses. The second is the development of high throughput screening assays to develop causal relationships between lncRNAs across the genome with cellular functions. Finally, we turn to recent advances in understanding the role of lncRNAs in regulating chromatin architecture and epigenetics, advances that build on some of the earliest work linking RNA to chromatin architecture.

DNA packaging by molecular motors: from bacteriophage to human chromosomes

Dense packaging of genomic DNA is crucial for organismal survival, as DNA length always far exceeds the dimensions of the cells that contain it. Organisms, therefore, use sophisticated machineries to package their genomes. These systems range across kingdoms from a single ultra-powerful rotary motor that spools the DNA into a bacteriophage head, to hundreds of thousands of relatively weak molecular motors that coordinate the compaction of mitotic chromosomes in eukaryotic cells. Recent technological advances, such as DNA proximity-based sequencing approaches, polymer modelling and in vitro reconstitution of DNA loop extrusion, have shed light on the biological mechanisms driving DNA organization in different systems. Here, we discuss DNA packaging in bacteriophage, bacteria and eukaryotic cells, which, despite their extreme variation in size, structure and genomic content, all rely on the action of molecular motors to package their genomes.

Temperature regulates negative supercoils to modulate meiotic crossovers and chromosome organization

Crossover recombination is a hallmark of meiosis that holds the paternal and maternal chromosomes (homologs) together for their faithful segregation, while promoting genetic diversity of the progeny. The pattern of crossover is mainly controlled by the architecture of the meiotic chromosomes. Environmental factors, especially temperature, also play an important role in modulating crossovers. However, it is unclear how temperature affects crossovers. Here, we examined the distribution of budding yeast axis components (Red1, Hop1, and Rec8) and the crossover-associated Zip3 foci in detail at different temperatures, and found that both increased and decreased temperatures result in shorter meiotic chromosome axes and more crossovers. Further investigations showed that temperature changes coordinately enhanced the hyperabundant accumulation of Hop1 and Red1 on chromosomes and the number of Zip3 foci. Most importantly, temperature-induced changes in the distribution of axis proteins and Zip3 foci depend on changes in DNA negative supercoils. These results suggest that yeast meiosis senses temperature changes by increasing the level of negative supercoils to increase crossovers and modulate chromosome organization. These findings provide a new perspective on understanding the effect and mechanism of temperature on meiotic recombination and chromosome organization, with important implications for evolution and breeding.

The chromatin tapestry as a framework for neurodevelopment

The neuronal nucleus houses a meticulously organized genome. Within this structure, genetic material is not simply compacted but arranged into a precise and functional 3D chromatin landscape essential for cellular regulation. This mini-review highlights the importance of this chromatin landscape in healthy neurodevelopment, as well as the diseases that occur with aberrant chromatin architecture. We discuss insights into the fundamental mechanistic relationship between histone modifications, DNA methylation, and genome organization. We then discuss findings that reveal how these epigenetic features change throughout normal neurodevelopment. Finally, we highlight single-gene neurodevelopmental disorders that illustrate the interdependence of epigenetic features, showing how disruptions in DNA methylation or genome architecture can ripple across the entire epigenome. As such, we emphasize the importance of measuring multiple chromatin architectural aspects, as the disruption of one mechanism can likely impact others in the intricate epigenetic network. This mini-review underscores the vast gaps in our understanding of chromatin structure in neurodevelopmental diseases and the substantial research needed to understand the interplay between chromatin features and neurodevelopment.

An integrated view of the structure and function of the human 4D nucleome

The dynamic three-dimensional (3D) organization of the human genome (the "4D Nucleome") is closely linked to genome function. Here, we integrate a wide variety of genomic data generated by the 4D Nucleome Project to provide a detailed view of human 3D genome organization in widely used embryonic stem cells (H1-hESCs) and immortalized fibroblasts (HFFc6). We provide extensive benchmarking of 3D genome mapping assays and integrate these diverse datasets to annotate spatial genomic features across scales. The data reveal a rich complexity of chromatin domains and their sub-nuclear positions, and over one hundred thousand structural loops and promoter-enhancer interactions. We developed 3D models of population-based and individual cell-to-cell variation in genome structure, establishing connections between chromosome folding, nuclear organization, chromatin looping, gene transcription, and DNA replication. We demonstrate the use of computational methods to predict genome folding from DNA sequence, uncovering potential effects of genetic variants on genome structure and function. Together, this comprehensive analysis contributes insights into human genome organization and enhances our understanding of connections between the regulation of genome function and 3D genome organization in general.

3D genome topology distinguishes molecular subgroups of medulloblastoma

Four main medulloblastoma (MB) molecular subtypes have been identified based on transcriptional, DNA methylation, and genetic profiles. However, it is currently not known whether 3D genome architecture differs between MB subtypes. To address this question, we performed in situ Hi-C to reconstruct the 3D genome architecture of MB subtypes. In total, we generated Hi-C and matching transcriptome data for 28 surgical specimens and Hi-C data for one patient-derived xenograft. The average resolution of the Hi-C maps was 6,833 bp. Using these data, we found that insulation scores of topologically associating domains (TADs) were effective at distinguishing MB molecular subgroups. TAD insulation score differences between subtypes were globally not associated with differential gene expression, although we identified few exceptions near genes expressed in the lineages of origin of specific MB subtypes. Our study therefore supports the notion that TAD insulation scores can distinguish MB subtypes independently of their transcriptional differences.

Targeting the 3D genome by anthracyclines for chemotherapeutic effects

The chromatins are folded into three-dimensional (3D) structures inside cells, which coordinates the regulation of gene transcription by the non-coding regulatory elements. Aberrant chromatin 3D folding has been shown in many diseases, such as acute myeloid leukemia (AML), and may contribute to tumorigenesis. The anthracycline topoisomerase II inhibitors can induce histone eviction and DNA damage. We performed genome-wide high-resolution mapping of the chemotherapeutic effects of various clinically used anthracycline drugs. ATAC-seq was used to profile the histone eviction effects of different anthracyclines. TOP2A ChIP-seq was used to profile the potential DNA damage regions. Integrated analyses show that different anthracyclines have distinct target selectivity on epigenomic regions, based on their respective ATAC-seq and ChIP-seq profiles. We identified the underlying molecular mechanism that unique anthracycline variants selectively target chromatin looping anchors via disrupting CTCF binding, suggesting an additional potential therapeutic effect on the 3D genome. We further performed Hi-C experiments, and data from K562 cells treated with the selective anthracycline drugs indicate that the 3D chromatin organization is disrupted. Furthermore, AML patients receiving anthracycline drugs showed altered chromatin structures around potential looping anchors, which linked to distinct clinical outcomes. Our data indicate that anthracyclines are potent and selective epigenomic targeting drugs and can target the 3D genome for anticancer therapy, which could be used for personalized medicine to treat tumors with aberrant 3D chromatin structures.

DFF-ChIP: a method to detect and quantify complex interactions between RNA polymerase II, transcription factors, and chromatin

Recently, we introduced a chromatin immunoprecipitation (ChIP) technique utilizing the human DNA Fragmentation Factor (DFF) to digest the DNA prior to immunoprecipitation (DFF-ChIP) that provides the precise location of transcription complexes and their interactions with neighboring nucleosomes. Here we expand the technique to new targets and provide useful information concerning purification of DFF, digestion conditions, and the impact of crosslinking. DFF-ChIP analysis was performed individually for subunits of Mediator, DSIF, and NELF that that do not interact with DNA directly, but rather interact with RNA polymerase II (Pol II). We found that Mediator was associated almost exclusively with preinitiation complexes (PICs). DSIF and NELF were associated with engaged Pol II and, in addition, potential intermediates between PICs and early initiation complexes. DFF-ChIP was then used to analyze the occupancy of a tight binding transcription factor, CTCF, and a much weaker binding factor, glucocorticoid receptor (GR), with and without crosslinking. These results were compared to those from standard ChIP-Seq that employs sonication and to CUT&RUN which utilizes MNase to fragment the genomic DNA. Our findings indicate that DFF-ChIP reveals details of occupancy that are not available using other methods including information revealing pertinent protein:protein interactions.

Dynamic barriers modulate cohesin positioning and genome folding at fixed occupancy

In mammalian interphase cells, genomes are folded by cohesin loop extrusion limited by directional CTCF barriers. This interplay leads to the enrichment of cohesin at barriers, isolation between neighboring topologically associating domains, and elevated contact frequency between convergent CTCF barriers across the genome. However, recent in vivo measurements present a puzzle: reported residence times for CTCF on chromatin are in the range of a few minutes, while lifetimes for cohesin are much longer. Can the observed features of genome folding result from the action of relatively transient barriers? To address this question, we developed a dynamic barrier model, where CTCF sites switch between bound and unbound states with rates that can be directly compared with biophysical measurements. Using this model, we investigated how barrier dynamics would impact observables for a range of experimental genomic and imaging datasets, including ChIP-seq, Hi-C, and microscopy. We found the interplay of CTCF and cohesin binding timescales influence the strength of each of these features, leaving a signature of barrier dynamics even in the population-averaged snapshots offered by genomic datasets. First, in addition to barrier occupancy, barrier bound times are crucial for instructing features of genome folding. Second, the ratio of boundary to extruder lifetime greatly alters simulated ChIP-seq and simulated Hi-C. Third, large-scale changes in chromosome morphology observed experimentally after increasing extruder lifetime require dynamic barriers. By integrating multiple sources of experimental data, our biophysical model argues that CTCF barrier bound times effectively approach those of cohesin extruder lifetimes. Together, we demonstrate how models that are informed by biophysically measured protein dynamics broaden our understanding of genome folding.

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.

The MVPs (masterful versatile players): Chromatin factors as pivotal mediators between 3D genome organization and the response to environment

In recent years, the study of genome dynamics has become a prominent research field due to its influence on understanding the control of gene expression. The study of 3D genome organization has unveiled multiple mechanisms in orchestrating chromosome folding. Growing evidence reveals that these mechanisms are not only important for genome organization, but play a pivotal role in enabling plants to adapt to environmental stimuli. In this review, we provide an overview of the current knowledge concerning epigenetic factors and regulatory elements driving 3D genome dynamics and their responses to external stimuli. We discuss the most recent findings, previous evidence, and explore their implications for future research.

Factors that determine cell type-specific CTCF binding in health and disease

A number of factors contribute to cell type-specific CTCF chromatin binding, but how they act in concert to determine binding stability and functionality has not been fully elucidated. In this review, we tie together different layers of regulation to provide a holistic view of what is known. What emerges from these studies is a multifaceted system in which DNA sequence, DNA and chromatin accessibility, and cell type-specific transcription factors together contribute to CTCF binding profile and function. We discuss these findings in the light of disease settings in which changes in the chromatin landscape and transcriptional programming can disrupt CTCF's binding profile and involvement in looping.

Structure and dynamics of nuclear A/B compartments and subcompartments

Mammalian chromosomes form a hierarchical structure within the cell nucleus, from chromatin loops, megabase (Mb)-sized topologically associating domains (TADs) to larger-scale A/B compartments. The molecular basis of the structures of loops and TADs has been actively studied. However, the A and B compartments, which correspond to early-replicating euchromatin and late-replicating heterochromatin, respectively, are still relatively unexplored. In this review, we focus on the A/B compartments, discuss their close relationship to DNA replication timing (RT), and introduce recent findings on the features of subcompartments revealed by detailed classification of the A/B compartments. In doing so, we speculate on the structure, potential function, and developmental dynamics of A/B compartments and subcompartments in mammalian cells.

Spandrels of the cell nucleus

S.J. Gould and R. Lewontin in their famous "Spandrels paper" (1979) argued that many anatomical elements arise in evolution not due to their "current utility" but rather due to other "reasons for origin", such as other developmental processes, physical constraints and mechanical forces. Here, in the same spirit, we argue that a variety of molecular processes, physical constraints, and mechanical forces, alone or together, generate structures that are detectable in the cell nucleus, yet these structures themselves may not carry any specific function, being a mere reflection of processes that produced them.

Permeable TAD boundaries and their impact on genome-associated functions

TAD boundaries are genomic elements that separate biological processes in neighboring domains by blocking DNA loops that are formed through Cohesin-mediated loop extrusion. Most TAD boundaries consist of arrays of binding sites for the CTCF protein, whose interaction with the Cohesin complex blocks loop extrusion. TAD boundaries are not fully impermeable though and allow a limited amount of inter-TAD loop formation. Based on the reanalysis of Nano-C data, a multicontact Chromosome Conformation Capture assay, we propose a model whereby clustered CTCF binding sites promote the successive stalling of Cohesin and subsequent dissociation from the chromatin. A fraction of Cohesin nonetheless achieves boundary read-through. Due to a constant rate of Cohesin dissociation elsewhere in the genome, the maximum length of inter-TAD loops is restricted though. We speculate that the DNA-encoded organization of stalling sites regulates TAD boundary permeability and discuss implications for enhancer-promoter loop formation and other genomic processes.

Pushing the TAD boundary: Decoding insulator codes of clustered CTCF sites in 3D genomes

Topologically associating domain (TAD) boundaries are the flanking edges of TADs, also known as insulated neighborhoods, within the 3D structure of genomes. A prominent feature of TAD boundaries in mammalian genomes is the enrichment of clustered CTCF sites often with mixed orientations, which can either block or facilitate enhancer-promoter (E-P) interactions within or across distinct TADs, respectively. We will discuss recent progress in the understanding of fundamental organizing principles of the clustered CTCF insulator codes at TAD boundaries. Specifically, both inward- and outward-oriented CTCF sites function as topological chromatin insulators by asymmetrically blocking improper TAD-boundary-crossing cohesin loop extrusion. In addition, boundary stacking and enhancer clustering facilitate long-distance E-P interactions across multiple TADs. Finally, we provide a unified mechanism for RNA-mediated TAD boundary function via R-loop formation for both insulation and facilitation. This mechanism of TAD boundary formation and insulation has interesting implications not only on how the 3D genome folds in the Euclidean nuclear space but also on how the specificity of E-P interactions is developmentally regulated.

Crosslinking intensity modulates the reliability and sensitivity of chromatin conformation detection at different structural levels

Formaldehyde (FA) is a chemical that facilitates crosslinking between DNA and proteins. It is widely used in various biochemical assays, such as chromosome conformation capture (3C) and Chromatin Immunoprecipitation (ChIP). While the concentration and temperature of FA treatment are recognized as crucial factors in crosslinking, their quantitative effects have largely remained unexplored. In this study, we employed 3C as a model system to systematically assess the impacts of these two factors on crosslinking. Our findings indicate that the strength of crosslinking significantly influences chromatin conformation detection at nearly all known structural levels. Specifically, a delicate balance between sensitivity and reliability is required when detecting higher-level structures, such as chromosome compartments. Conversely, intense crosslinking is preferred when targeting lower-level structures, such as topologically associated domains (TADs) or chromatin loops. Based on our data, we propose a conceptual molecular thermal motion model to elucidate the roles of these two factors in restricting FA crosslinking. Our results not only shed light on the previously overlooked confounding factor in FA crosslinking but also highlight the need for caution in new technology developments that rely on FA crosslinking.

Local volume concentration, packing domains, and scaling properties of chromatin

We propose the Self Returning Excluded Volume (SR-EV) model for the structure of chromatin based on stochastic rules and physical interactions. The SR-EV rules of return generate conformationally defined domains observed by single-cell imaging techniques. From nucleosome to chromosome scales, the model captures the overall chromatin organization as a corrugated system, with dense and dilute regions alternating in a manner that resembles the mixing of two disordered bi-continuous phases. This particular organizational topology is a consequence of the multiplicity of interactions and processes occurring in the nuclei, and mimicked by the proposed return rules. Single configuration properties and ensemble averages show a robust agreement between theoretical and experimental results including chromatin volume concentration, contact probability, packing domain identification and size characterization, and packing scaling behavior. Model and experimental results suggest that there is an inherent chromatin organization regardless of the cell character and resistant to an external forcing such as RAD21 degradation.

Brain and cancer associated 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 be reprogrammed. Collectively, the mutant perturbations provide a rich resource for determining CTCF's site-specific effects.

Prediction and functional interpretation of inter-chromosomal genome architecture from DNA sequence with TwinC

Three-dimensional nuclear DNA architecture comprises well-studied intra-chromosomal (cis) folding and less characterized inter-chromosomal (trans) interfaces. Current predictive models of 3D genome folding can effectively infer pairwise cis-chromatin interactions from the primary DNA sequence but generally ignore trans contacts. There is an unmet need for robust models of trans-genome organization that provide insights into their underlying principles and functional relevance. We present TwinC, an interpretable convolutional neural network model that reliably predicts trans contacts measurable through genome-wide chromatin conformation capture (Hi-C). TwinC uses a paired sequence design from replicate Hi-C experiments to learn single base pair relevance in trans interactions across two stretches of DNA. The method achieves high predictive accuracy (AUROC=0.80) on a cross-chromosomal test set from Hi-C experiments in heart tissue. Mechanistically, the neural network learns the importance of compartments, chromatin accessibility, clustered transcription factor binding and G-quadruplexes in forming trans contacts. In summary, TwinC models and interprets trans genome architecture, shedding light on this poorly understood aspect of gene regulation.

Polymer Physics Models Reveal Structural Folding Features of Single-Molecule Gene Chromatin Conformations

Here, we employ polymer physics models of chromatin to investigate the 3D folding of a 2 Mb wide genomic region encompassing the human LTN1 gene, a crucial DNA locus involved in key cellular functions. Through extensive Molecular Dynamics simulations, we reconstruct in silico the ensemble of single-molecule LTN1 3D structures, which we benchmark against recent in situ Hi-C 2.0 data. The model-derived single molecules are then used to predict structural folding features at the single-cell level, providing testable predictions for super-resolution microscopy experiments.

Genetic-epigenetic interplay in the determination of plant 3D genome organization

The 3D chromatin organization plays a major role in the control of gene expression. However, our comprehension of the governing principles behind nuclear organization remains incomplete. Particularly, the spatial segregation of loci with similar repressive transcriptional states in plants poses a significant yet poorly understood puzzle. In this study, employing a combination of genetics and advanced 3D genomics approaches, we demonstrated that a redistribution of facultative heterochromatin marks in regions usually occupied by constitutive heterochromatin marks disrupts the 3D genome compartmentalisation. This disturbance, in turn, triggers novel chromatin interactions between genic and transposable element (TE) regions. Interestingly, our results imply that epigenetic features, constrained by genetic factors, intricately mold the landscape of 3D genome organisation. This study sheds light on the profound genetic-epigenetic interplay that underlies the regulation of gene expression within the intricate framework of the 3D genome. Our findings highlight the complexity of the relationships between genetic determinants and epigenetic features in shaping the dynamic configuration of the 3D genome.

Cohesin-mediated chromatin remodeling controls the differentiation and function of conventional dendritic cells

The cohesin protein complex extrudes chromatin loops, stopping at CTCF-bound sites, to organize chromosomes into topologically associated domains, yet the biological implications of this process are poorly understood. We show that cohesin is required for the post-mitotic differentiation and function of antigen-presenting dendritic cells (DCs), particularly for antigen cross-presentation and IL-12 secretion by type 1 conventional DCs (cDC1s) in vivo. The chromatin organization of DCs was shaped by cohesin and the DC-specifying transcription factor IRF8, which controlled chromatin looping and chromosome compartmentalization, respectively. Notably, optimal expression of IRF8 itself required CTCF/cohesin-binding sites demarcating the Irf8 gene. During DC activation, cohesin was required for the induction of a subset of genes with distal enhancers. Accordingly, the deletion of CTCF sites flanking the Il12b gene reduced IL-12 production by cDC1s. Our data reveal an essential role of cohesin-mediated chromatin regulation in cell differentiation and function in vivo, and its bi-directional crosstalk with lineage-specifying transcription factors.

Three-dimensional chromatin mapping of sensory neurons reveals that promoter-enhancer looping is required for axonal regeneration

The in vivo three-dimensional genomic architecture of adult mature neurons at homeostasis and after medically relevant perturbations such as axonal injury remains elusive. Here, we address this knowledge gap by mapping the three-dimensional chromatin architecture and gene expression program at homeostasis and after sciatic nerve injury in wild-type and cohesin-deficient mouse sensory dorsal root ganglia neurons via combinatorial Hi-C, promoter-capture Hi-C, CUT&Tag for H3K27ac and RNA-seq. We find that genes involved in axonal regeneration form long-range, complex chromatin loops, and that cohesin is required for the full induction of the regenerative transcriptional program. Importantly, loss of cohesin results in disruption of chromatin architecture and severely impaired nerve regeneration. Complex enhancer-promoter loops are also enriched in the human fetal cortical plate, where the axonal growth potential is highest, and are lost in mature adult neurons. Together, these data provide an original three-dimensional chromatin map of adult sensory neurons in vivo and demonstrate a role for cohesin-dependent long-range promoter interactions in nerve regeneration.

Dynamics of microcompartment formation at the mitosis-to-G1 transition

As cells exit mitosis and enter G1, mitotic chromosomes decompact and transcription is reestablished. Previously, Hi-C studies showed that essentially all interphase 3D genome features including A/B-compartments, TADs, and CTCF loops, are lost during mitosis. However, Hi-C remains insensitive to features such as microcompartments, nested focal interactions between cis-regulatory elements (CREs). We therefore applied Region Capture Micro-C to cells from mitosis to G1. Unexpectedly, we observe microcompartments in prometaphase, which further strengthen in ana/telophase before gradually weakening in G1. Loss of loop extrusion through condensin depletion differentially impacts microcompartments and large A/B-compartments, suggesting that they are partially distinct. Using polymer modeling, we show that microcompartment formation is favored by chromatin compaction and disfavored by loop extrusion activity, explaining why ana/telophase likely provides a particularly favorable environment. Our results suggest that CREs exhibit intrinsic homotypic affinity leading to microcompartment formation, which may explain transient transcriptional spiking observed upon mitotic exit.

Interphase chromosome conformation is specified by distinct folding programs inherited via mitotic chromosomes or through the cytoplasm

Identity-specific interphase chromosome conformation must be re-established each time a cell divides. To understand how interphase folding is inherited, we developed an experimental approach that physically segregates mediators of G1 folding that are intrinsic to mitotic chromosomes from cytoplasmic factors. Proteins essential for nuclear transport, RanGAP1 and Nup93, were degraded in pro-metaphase arrested DLD-1 cells to prevent the establishment of nucleo-cytoplasmic transport during mitotic exit and isolate the decondensing mitotic chromatin of G1 daughter cells from the cytoplasm. Using this approach, we discover a transient folding intermediate entirely driven by chromosome-intrinsic factors. In addition to conventional compartmental segregation, this chromosome-intrinsic folding program leads to prominent genome-scale microcompartmentalization of mitotically bookmarked and cell type-specific cis-regulatory elements. This microcompartment conformation is formed during telophase and subsequently modulated by a second folding program driven by factors inherited through the cytoplasm in G1. This nuclear import-dependent folding program includes cohesin and factors involved in transcription and RNA processing. The combined and inter-dependent action of chromosome-intrinsic and cytoplasmic inherited folding programs determines the interphase chromatin conformation as cells exit mitosis.

HYENA detects oncogenes activated by distal enhancers in cancer

Somatic structural variations (SVs) in cancer can shuffle DNA content in the genome, relocate regulatory elements, and alter genome organization. Enhancer hijacking occurs when SVs relocate distal enhancers to activate proto-oncogenes. However, most enhancer hijacking studies have only focused on protein-coding genes. Here, we develop a computational algorithm 'HYENA' to identify candidate oncogenes (both protein-coding and non-coding) activated by enhancer hijacking based on tumor whole-genome and transcriptome sequencing data. HYENA detects genes whose elevated expression is associated with somatic SVs by using a rank-based regression model. We systematically analyze 1146 tumors across 25 types of adult tumors and identify a total of 108 candidate oncogenes including many non-coding genes. A long non-coding RNA TOB1-AS1 is activated by various types of SVs in 10% of pancreatic cancers through altered 3-dimensional genome structure. We find that high expression of TOB1-AS1 can promote cell invasion and metastasis. Our study highlights the contribution of genetic alterations in non-coding regions to tumorigenesis and tumor progression.

Spectral-based detection of chromatin loops in multiplexed super-resolution FISH data

Involved in mitotic condensation, interaction of transcriptional regulatory elements and isolation of structural domains, loop formation has become a paradigm in the deciphering of chromatin architecture and its functional role. Despite the emergence of increasingly powerful genome visualization techniques, the high variability in cell populations and the randomness of conformations still make loop detection a challenge. We introduce an approach for determining the presence and frequency of loops in a collection of experimental conformations obtained by multiplexed super-resolution imaging. Based on a spectral approach, in conjunction with neural networks, this method offers a powerful tool to detect loops in large experimental data sets, both at the population and single-cell levels. The method's performance is confirmed on experimental FISH data where Hi-C and other loop detection results are available. The method is then applied to recently published experimental data, where it provides a detailed and statistically quantified description of the global architecture of the chromosomal region under study.

YY1-controlled regulatory connectivity and transcription are influenced by the cell cycle

Few transcription factors have been examined for their direct roles in physically connecting enhancers and promoters. Here acute degradation of Yin Yang 1 (YY1) in erythroid cells revealed its requirement for the maintenance of numerous enhancer-promoter loops, but not compartments or domains. Despite its reported ability to interact with cohesin, the formation of YY1-dependent enhancer-promoter loops does not involve stalling of cohesin-mediated loop extrusion. Integrating mitosis-to-G1-phase dynamics, we observed partial retention of YY1 on mitotic chromatin, predominantly at gene promoters, followed by rapid rebinding during mitotic exit, coinciding with enhancer-promoter loop establishment. YY1 degradation during the mitosis-to-G1-phase interval revealed a set of enhancer-promoter loops that require YY1 for establishment during G1-phase entry but not for maintenance in interphase, suggesting that cell cycle stage influences YY1's architectural function. Thus, as revealed here for YY1, chromatin architectural functions of transcription factors can vary in their interplay with CTCF and cohesin as well as by cell cycle stage.

A Comparison of Two Versions of the CRISPR-Sirius System for the Live-Cell Visualization of the Borders of Topologically Associating Domains

In recent years, various technologies have emerged for the imaging of chromatin loci in living cells via catalytically inactive Cas9 (dCas9). These technologies facilitate a deeper understanding of the mechanisms behind the chromatin dynamics and provide valuable kinetic data that could not have previously been obtained via FISH applied to fixed cells. However, such technologies are relatively complicated, as they involve the expression of several chimeric proteins as well as sgRNAs targeting the visualized loci, a process that entails many technical subtleties. Therefore, the effectiveness in visualizing a specific target locus may be quite low. In this study, we directly compared two versions of a previously published CRISPR-Sirius method based on the use of sgRNAs containing eight MS2 or PP7 stem loops and the expression of MCP or PCP fused to fluorescent proteins. We assessed the visualization efficiency for several unique genomic loci by comparing the two approaches in delivering sgRNA genes (transient transfection and lentiviral transduction), as well as two CRISPR-Sirius versions (with PCP and with MCP). The efficiency of visualization varied among the loci, and not all loci could be visualized. However, the MCP-sfGFP version provided more efficient visualization in terms of the number of cells with signals than PCP-sfGFP for all tested loci. We also showed that lentiviral transduction was more efficient in locus imaging than transient transfection for both CRISPR-Sirius systems. Most of the target loci in our study were located at the borders of topologically associating domains, and we defined a set of TAD borders that could be effectively visualized using the MCP-sfGFP version of the CRISPR-Sirius system. Altogether, our study validates the use of the CRISPR-Sirius technology for live-cell visualization and highlights various technical details that should be considered when using this method.

Unraveling the three-dimensional (3D) genome architecture in Neurodevelopmental Disorders (NDDs)

The human genome, comprising millions of pairs of bases, serves as the blueprint of life, encoding instructions for cellular processes. However, genomes are not merely linear sequences; rather, the complex of DNA and histones, known as chromatin, exhibits complex organization across various levels, which profoundly influence gene expression and cellular function. Central to understanding genome organization is the emerging field of three-dimensional (3D) genome studies. Utilizing advanced techniques such as Hi-C, researchers have unveiled non-random dispositions of genomic elements, highlighting their importance in transcriptional regulation and disease mechanisms. Topologically Associating Domains (TADs), that demarcate regions of chromatin with preferential internal interactions, play crucial roles in gene regulation and are increasingly implicated in various diseases such as cancer and schizophrenia. However, their role in Neurodevelopmental Disorders (NDDs) remains poorly understood. Here, we focus on TADs and 3D conservation across the evolution and between cell types in NDDs. The investigation into genome organization and its impact on disease has led to significant breakthroughs in understanding NDDs etiology such ASD (Autism Spectrum Disorder). By elucidating the wide spectrum of ASD manifestations, researchers aim to uncover the underlying genetic and epigenetic factors contributing to its heterogeneity. Moreover, studies linking TAD disruption to NDDs underscore the importance of spatial genome organization in maintaining proper brain development and function. In summary, this review highlights the intricate interplay between genome organization, transcriptional control, and disease pathology, shedding light on fundamental biological processes and offering insights into the mechanisms underlying NDDs like ASD.

The physical chemistry of interphase loop extrusion

Loop extrusion constitutes a universal mechanism of genome organization, whereby structural maintenance of chromosomes (SMC) protein complexes load onto the chromatin fiber and generate DNA loops of increasingly-larger sizes until their eventual release. In mammalian interphase cells, loop extrusion is mediated by the cohesin complex, which is dynamically regulated by the interchange of multiple accessory proteins. Although these regulators bind the core cohesin complex only transiently, their disruption can dramatically alter cohesin dynamics, gene expression, chromosome morphology and contact patterns. Still, a theory of how cohesin regulators and their molecular interplay with the core complex modulate genome folding remains at large. Here we derive a model of cohesin loop extrusion from first principles, based on in vivo measurements of the abundance and dynamics of cohesin regulators. We systematically evaluate potential chemical reaction networks that describe the association of cohesin with its regulators and with the chromatin fiber. Remarkably, experimental observations are consistent with only a single biochemical reaction cycle, which results in a unique minimal model that may be fully parameterized by quantitative protein measurements. We demonstrate how distinct roles for cohesin regulators emerge simply from the structure of the reaction network, and how their dynamic exchange can regulate loop extrusion kinetics over time-scales that far exceed their own chromatin residence times. By embedding our cohesin biochemical reaction network within biophysical chromatin simulations, we evidence how variations in regulatory protein abundance can alter chromatin architecture across multiple length- and time-scales. Predictions from our model are corroborated by biophysical and biochemical assays, optical microscopy observations, and Hi-C conformation capture techniques. More broadly, our theoretical and numerical framework bridges the gap between in vitro observations of extrusion motor dynamics at the molecular scale and their structural consequences at the genome-wide level.

LDB1 establishes multi-enhancer networks to regulate gene expression

How specific enhancer-promoter pairing is established is still mostly unclear. Besides the CTCF/cohesin machinery, only a few nuclear factors have been studied for a direct role in physically connecting regulatory elements. Here, we show via acute degradation experiments that LDB1 directly and broadly promotes enhancer-promoter loops. Most LDB1-mediated contacts, even those spanning hundreds of kb, can form in the absence of CTCF, cohesin, or YY1 as determined via the use of multiple degron systems. Moreover, an engineered LDB1-driven chromatin loop is cohesin independent. Cohesin-driven loop extrusion does not stall at LDB1 occupied sites but may aid the formation of a subset of LDB1 anchored loops. Leveraging the dynamic reorganization of nuclear architecture during the transition from mitosis to G1-phase, we establish a relationship between LDB1-dependent interactions in the context of TAD organization and gene activation. Lastly, Tri-C and Region Capture Micro-C reveal that LDB1 organizes multi-enhancer networks to activate transcription. This establishes LDB1 as a direct driver of regulatory network inter-connectivity.

IVEA: an integrative variational Bayesian inference method for predicting enhancer-gene regulatory interactions

Motivation

Enhancers play critical roles in cell-type-specific transcriptional control. Despite the identification of thousands of candidate enhancers, unravelling their regulatory relationships with their target genes remains challenging. Therefore, computational approaches are needed to accurately infer enhancer-gene regulatory relationships.

Results

In this study, we propose a new method, IVEA, that predicts enhancer-gene regulatory interactions by estimating promoter and enhancer activities. Its statistical model is based on the gene regulatory mechanism of transcriptional bursting, which is characterized by burst size and frequency controlled by promoters and enhancers, respectively. Using transcriptional readouts, chromatin accessibility, and chromatin contact data as inputs, promoter and enhancer activities were estimated using variational Bayesian inference, and the contribution of each enhancer-promoter pair to target gene transcription was calculated. Our analysis demonstrates that the proposed method can achieve high prediction accuracy and provide biologically relevant enhancer-gene regulatory interactions.

Availability and implementation

The IVEA code is available on GitHub at https://github.com/yasumasak/ivea. The publicly available datasets used in this study are described in Supplementary Table S4.

OpenNucleome for high-resolution nuclear structural and dynamical modeling

The intricate structural organization of the human nucleus is fundamental to cellular function and gene regulation. Recent advancements in experimental techniques, including high-throughput sequencing and microscopy, have provided valuable insights into nuclear organization. Computational modeling has played significant roles in interpreting experimental observations by reconstructing high-resolution structural ensembles and uncovering organization principles. However, the absence of standardized modeling tools poses challenges for furthering nuclear investigations. We present OpenNucleome-an open-source software designed for conducting GPU-accelerated molecular dynamics simulations of the human nucleus. OpenNucleome offers particle-based representations of chromosomes at a resolution of 100 KB, encompassing nuclear lamina, nucleoli, and speckles. This software furnishes highly accurate structural models of nuclear architecture, affording the means for dynamic simulations of condensate formation, fusion, and exploration of non-equilibrium effects. We applied OpenNucleome to uncover the mechanisms driving the emergence of 'fixed points' within the nucleus-signifying genomic loci robustly anchored in proximity to specific nuclear bodies for functional purposes. This anchoring remains resilient even amidst significant fluctuations in chromosome radial positions and nuclear shapes within individual cells. Our findings lend support to a nuclear zoning model that elucidates genome functionality. We anticipate OpenNucleome to serve as a valuable tool for nuclear investigations, streamlining mechanistic explorations and enhancing the interpretation of experimental observations.

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.

3D chromatin architecture, BRD4, and Mediator have distinct roles in regulating genome-wide transcriptional bursting and gene network

Discontinuous transcription is evolutionarily conserved and a fundamental feature of gene regulation; yet, the exact mechanisms underlying transcriptional bursting are unresolved. Analyses of bursting transcriptome-wide have focused on the role of cis-regulatory elements, but other factors that regulate this process remain elusive. We applied mathematical modeling to single-cell RNA sequencing data to infer bursting dynamics transcriptome-wide under multiple conditions to identify possible molecular mechanisms. We found that Mediator complex subunit 26 (MED26) primarily regulates frequency, MYC regulates burst size, while cohesin and Bromodomain-containing protein 4 (BRD4) can modulate both. Despite comparable effects on RNA levels among these perturbations, acute depletion of MED26 had the most profound impact on the entire gene regulatory network, acting downstream of chromatin spatial architecture and without affecting TATA box-binding protein (TBP) recruitment. These results indicate that later steps in the initiation of transcriptional bursts are primary nodes for integrating gene networks in single cells.

Chromosome structure in Drosophila is determined by boundary pairing not loop extrusion

Two different models have been proposed to explain how the endpoints of chromatin looped domains ('TADs') in eukaryotic chromosomes are determined. In the first, a cohesin complex extrudes a loop until it encounters a boundary element roadblock, generating a stem-loop. In this model, boundaries are functionally autonomous: they have an intrinsic ability to halt the movement of incoming cohesin complexes that is independent of the properties of neighboring boundaries. In the second, loops are generated by boundary:boundary pairing. In this model, boundaries are functionally non-autonomous, and their ability to form a loop depends upon how well they match with their neighbors. Moreover, unlike the loop-extrusion model, pairing interactions can generate both stem-loops and circle-loops. We have used a combination of MicroC to analyze how TADs are organized, and experimental manipulations of the even skipped TAD boundary, homie, to test the predictions of the 'loop-extrusion' and the 'boundary-pairing' models. Our findings are incompatible with the loop-extrusion model, and instead suggest that the endpoints of TADs in flies are determined by a mechanism in which boundary elements physically pair with their partners, either head-to-head or head-to-tail, with varying degrees of specificity. Although our experiments do not address how partners find each other, the mechanism is unlikely to require loop extrusion.

Structural perturbation of chromatin domains with multiple developmental regulators can severely impact gene regulation and development

Chromatin domain boundaries delimited by CTCF motifs can restrict the range of enhancer action. However, disruption of domain structure often results in mild gene dysregulation and thus predicting the impact of boundary rearrangements on animal development remains challenging. Here, we tested whether structural perturbation of a chromatin domain with multiple developmental regulators can result in more acute gene dysregulation and severe developmental phenotypes. We targeted clusters of CTCF motifs in a domain of the mouse genome containing three FGF ligand genes-Fgf3, Fgf4, and Fgf15-that regulate several developmental processes. Deletion of the 23.9kb cluster that defines the centromeric boundary of this domain resulted in ectopic interactions of the FGF genes with enhancers located across the deleted boundary that are active in the developing brain. This caused strong induction of FGF expression and perinatal lethality with encephalocele and orofacial cleft phenotypes. Heterozygous boundary deletion was sufficient to cause these fully penetrant phenotypes, and strikingly, loss of a single CTCF motif within the cluster also recapitulated ectopic FGF expression and caused encephalocele. However, such phenotypic sensitivity to perturbation of domain structure did not extend to all CTCF clusters of this domain, nor to all developmental processes controlled by these three FGF genes-for example, the ability to undergo lineage specification in the blastocyst and pre-implantation development were not affected. By tracing the impact of different chromosomal rearrangements throughout mouse development, we start to uncover the determinants of phenotypic robustness and sensitivity to perturbation of chromatin boundaries. Our data show how small sequence variants at certain domain boundaries can have a surprisingly outsized effect and must be considered as potential sources of gene dysregulation during development and disease.

The relationship between nucleosome positioning and higher-order genome folding

Eukaryotic genomes are organized into 3D structures, which range from small-scale nucleosome arrays to large-scale chromatin domains. These structures have an important role in the regulation of transcription and other nuclear processes. Despite advances in our understanding of the properties, functions, and underlying mechanisms of genome structures, there are many open questions about the interplay between these structures across scales. In particular, it is not well understood if and how 1D features of nucleosome arrays influence large-scale 3D genome folding patterns. In this review, we discuss recent studies that address these questions and summarize our current understanding of the relationship between nucleosome positioning and higher-order genome folding.

Cohesin prevents cross-domain gene coactivation

The contrast between the disruption of genome topology after cohesin loss and the lack of downstream gene expression changes instigates intense debates regarding the structure-function relationship between genome and gene regulation. Here, by analyzing transcriptome and chromatin accessibility at the single-cell level, we discover that, instead of dictating population-wide gene expression levels, cohesin supplies a general function to neutralize stochastic coexpression tendencies of cis-linked genes in single cells. Notably, cohesin loss induces widespread gene coactivation and chromatin co-opening tens of million bases apart in cis. Spatial genome and protein imaging reveals that cohesin prevents gene co-bursting along the chromosome and blocks spatial mixing of transcriptional hubs. Single-molecule imaging shows that cohesin confines the exploration of diverse enhancer and core promoter binding transcriptional regulators. Together, these results support that cohesin arranges nuclear topology to control gene coexpression in single cells.

Cohesin composition and dosage independently affect early development in zebrafish

Cohesin, a chromatin-associated protein complex with four core subunits (Smc1a, Smc3, Rad21 and either Stag1 or 2), has a central role in cell proliferation and gene expression in metazoans. Human developmental disorders termed 'cohesinopathies' are characterized by germline variants of cohesin or its regulators that do not entirely eliminate cohesin function. However, it is not clear whether mutations in individual cohesin subunits have independent developmental consequences. Here, we show that zebrafish rad21 or stag2b mutants independently influence embryonic tailbud development. Both mutants have altered mesoderm induction, but only homozygous or heterozygous rad21 mutation affects cell cycle gene expression. stag2b mutants have narrower notochords and reduced Wnt signaling in neuromesodermal progenitors as revealed by single-cell RNA sequencing. Stimulation of Wnt signaling rescues transcription and morphology in stag2b, but not rad21, mutants. Our results suggest that mutations altering the quantity versus composition of cohesin have independent developmental consequences, with implications for the understanding and management of cohesinopathies.

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.

Unique territorial and compartmental organization of chromosomes in the holocentric silkmoth

The hallmarks of chromosome organization in multicellular eukaryotes are chromosome territories (CT), chromatin compartments, and insulated domains, including topologically associated domains (TADs). Yet, most of these elements of chromosome organization are derived from analyses of a limited set of model organisms, while large eukaryotic groups, including insects, remain mostly unexplored. Here we combine Hi-C, biophysical modeling, and microscopy to characterize the 3D genome architecture of the silkmoth, Bombyx mori. In contrast to other eukaryotes, B. mori chromosomes form highly separated territories. Similar to other eukaryotes, B. mori chromosomes segregate into active A and inactive B compartments, yet unlike in vertebrate systems, contacts between euchromatic A regions appear to be a strong driver of compartmentalization. Remarkably, we also identify a third compartment, called secluded "S," with a unique contact pattern. Each S region shows prominent short-range self-contacts and is remarkably devoid of contacts with the rest of the chromosome, including other S regions. Compartment S hosts a unique combination of genetic and epigenetic features, localizes towards the periphery of CTs, and shows developmental plasticity. Biophysical modeling reveals that the formation of such secluded domains requires highly localized loop extrusion within them, along with a low level of extrusion in A and B. Our Hi-C data supports predicted genome-wide and localized extrusion. Such a broad, non-uniform distribution of extruders has not been seen in other organisms. Overall, our analyses support loop extrusion in insects and highlight the evolutionary plasticity of 3D genome organization, driven by a new combination of known processes.

Exploring the roles of RNAs in chromatin architecture using deep learning

Recent studies have highlighted the impact of both transcription and transcripts on 3D genome organization, particularly its dynamics. Here, we propose a deep learning framework, called AkitaR, that leverages both genome sequences and genome-wide RNA-DNA interactions to investigate the roles of chromatin-associated RNAs (caRNAs) on genome folding in HFFc6 cells. In order to disentangle the cis- and trans-regulatory roles of caRNAs, we have compared models with nascent transcripts, trans-located caRNAs, open chromatin data, or DNA sequence alone. Both nascent transcripts and trans-located caRNAs improve the models' predictions, especially at cell-type-specific genomic regions. Analyses of feature importance scores reveal the contribution of caRNAs at TAD boundaries, chromatin loops and nuclear sub-structures such as nuclear speckles and nucleoli to the models' predictions. Furthermore, we identify non-coding RNAs (ncRNAs) known to regulate chromatin structures, such as MALAT1 and NEAT1, as well as several new RNAs, RNY5, RPPH1, POLG-DT and THBS1-IT1, that might modulate chromatin architecture through trans-interactions in HFFc6. Our modeling also suggests that transcripts from Alus and other repetitive elements may facilitate chromatin interactions through trans R-loop formation. Our findings provide insights and generate testable hypotheses about the roles of caRNAs in shaping chromatin organization.

CatLearning: highly accurate gene expression prediction from histone mark

Histone modifications, known as histone marks, are pivotal in regulating gene expression within cells. The vast array of potential combinations of histone marks presents a considerable challenge in decoding the regulatory mechanisms solely through biological experimental approaches. To overcome this challenge, we have developed a method called CatLearning. It utilizes a modified convolutional neural network architecture with a specialized adaptation Residual Network to quantitatively interpret histone marks and predict gene expression. This architecture integrates long-range histone information up to 500Kb and learns chromatin interaction features without 3D information. By using only one histone mark, CatLearning achieves a high level of accuracy. Furthermore, CatLearning predicts gene expression by simulating changes in histone modifications at enhancers and throughout the genome. These findings help comprehend the architecture of histone marks and develop diagnostic and therapeutic targets for diseases with epigenetic changes.

SnapHiC-G: identifying long-range enhancer-promoter interactions from single-cell Hi-C data via a global background model

Harnessing the power of single-cell genomics technologies, single-cell Hi-C (scHi-C) and its derived technologies provide powerful tools to measure spatial proximity between regulatory elements and their target genes in individual cells. Using a global background model, we propose SnapHiC-G, a computational method, to identify long-range enhancer-promoter interactions from scHi-C data. We applied SnapHiC-G to scHi-C datasets generated from mouse embryonic stem cells and human brain cortical cells. SnapHiC-G achieved high sensitivity in identifying long-range enhancer-promoter interactions. Moreover, SnapHiC-G can identify putative target genes for noncoding genome-wide association study (GWAS) variants, and the genetic heritability of neuropsychiatric diseases is enriched for single-nucleotide polymorphisms (SNPs) within SnapHiC-G-identified interactions in a cell-type-specific manner. In sum, SnapHiC-G is a powerful tool for characterizing cell-type-specific enhancer-promoter interactions from complex tissues and can facilitate the discovery of chromatin interactions important for gene regulation in biologically relevant cell types.

An enhancer RNA recruits KMT2A to regulate transcription of Myb

The Myb proto-oncogene encodes the transcription factor c-MYB, which is critical for hematopoiesis. Distant enhancers of Myb form a hub of interactions with the Myb promoter. We identified a long non-coding RNA (Myrlin) originating from the -81-kb murine Myb enhancer. Myrlin and Myb are coordinately regulated during erythroid differentiation. Myrlin TSS deletion using CRISPR-Cas9 reduced Myrlin and Myb expression and LDB1 complex occupancy at the Myb enhancers, compromising enhancer contacts and reducing RNA Pol II occupancy in the locus. In contrast, CRISPRi silencing of Myrlin left LDB1 and the Myb enhancer hub unperturbed, although Myrlin and Myb expressions were downregulated, decoupling transcription and chromatin looping. Myrlin interacts with the KMT2A/MLL1 complex. Myrlin CRISPRi compromised KMT2A occupancy in the Myb locus, decreasing CDK9 and RNA Pol II binding and resulting in Pol II pausing in the Myb first exon/intron. Thus, Myrlin directly participates in activating Myb transcription by recruiting KMT2A.

Cohesin-mediated 3D contacts tune enhancer-promoter regulation

Enhancers are key drivers of gene regulation thought to act via 3D physical interactions with the promoters of their target genes. However, genome-wide depletions of architectural proteins such as cohesin result in only limited changes in gene expression, despite a loss of contact domains and loops. Consequently, the role of cohesin and 3D contacts in enhancer function remains debated. Here, we developed CRISPRi of regulatory elements upon degron operation (CRUDO), a novel approach to measure how changes in contact frequency impact enhancer effects on target genes by perturbing enhancers with CRISPRi and measuring gene expression in the presence or absence of cohesin. We systematically perturbed all 1,039 candidate enhancers near five cohesin-dependent genes and identified 34 enhancer-gene regulatory interactions. Of 26 regulatory interactions with sufficient statistical power to evaluate cohesin dependence, 18 show cohesin-dependent effects. A decrease in enhancer-promoter contact frequency upon removal of cohesin is frequently accompanied by a decrease in the regulatory effect of the enhancer on gene expression, consistent with a contact-based model for enhancer function. However, changes in contact frequency and regulatory effects on gene expression vary as a function of distance, with distal enhancers (e.g., >50Kb) experiencing much larger changes than proximal ones (e.g., <50Kb). Because most enhancers are located close to their target genes, these observations can explain how only a small subset of genes - those with strong distal enhancers - are sensitive to cohesin. Together, our results illuminate how 3D contacts, influenced by both cohesin and genomic distance, tune enhancer effects on gene expression.