Chromatin alternates between A and B compartments at kilobase scale for subgenic organization

Divergent 3D genome organization in livers of cave and surface morphs of Astyanax mexicanus as a potential driver of unique metabolic adaptations in cave environment

The cave morphs of Astyanax mexicanus have evolved a suite of distinct adaptations to life in perpetual darkness, including the loss of eyes and pigmentation loss, as well as profound metabolic changes such as hyperphagia and starvation resilience, traits that sharply contrast with those of their river-dwelling surface counterparts. While changed gene expression is a primary driver of these adaptations, the underlying role of 3D genome organization - a key regulator of gene expression - remains unexplored. Here, we investigate the 3D genome architecture of the livers of surface fish and two cavefish morphs (Pachón and Tinaja) using Hi-C, performing the first comparative 3D genomic analysis in this species. We analyzed and identified cave-specific 3D genomic features, such as genomic compartments and loops, which were conserved in both the cave populations but absent in surface fish. Integrating the 3D genome data with transcriptomic and epigenetic datasets, linked these changes to differential expression of metabolically relevant genes, such as Arhgef19 and Endog . Additionally, our study also uncovered genomic inversions unique to cavefish, potentially tied to cave adaptation. Our findings suggest that 3D genome organization contributes to transcriptomic shifts underlying cavefish phenotypes, providing a novel intra-species and morph specific perspective on 3D chromatin evolution. This study establishes a foundation for exploring how genome architecture potentially facilitates adaptation to new environments. Comparison of morphs within the same species also establishes a foundation for better understanding of how 3D genome reorganization may drive speciation and phenotypic diversity.

Toward decoding the mechanisms that shape sub-megabase-scale genome organization

Understanding genome organization at the kilobase to megabase scale is critical, as it encompasses genes and regulatory elements. Improvements in the resolution of experimental techniques have revealed novel structural motifs at this scale, including micro-compartments, nucleosome clutches, microdomains, and packing domains. Here we review recent progress on developing theories to explain these observations. Key advances include elucidating the role of nucleosome positioning and epigenetic modifications, the role and mechanisms of compartmentalization in local structure, and the interplay between loop extrusion and phase separation. This work has revealed probable mechanisms by which the observed structures emerge, but it remains unclear how these factors act together in the cell. To this end, recent studies have used chromatin conformation capture data in concert with diverse genomics datasets to create native-like models of chromatin at nucleosome resolution and below. While several roadblocks remain, this strategy promises to decode how molecular forces sum to shape chromatin structure and ultimately regulate transcription.

Chromatin loops are an ancestral hallmark of the animal regulatory genome

In bilaterian animals, gene regulation is shaped by a combination of linear and spatial regulatory information. Regulatory elements along the genome are integrated into gene regulatory landscapes through chromatin compartmentalization1,2, insulation of neighbouring genomic regions3,4 and chromatin looping that brings together distal cis-regulatory sequences5. However, the evolution of these regulatory features is unknown because the three-dimensional genome architecture of most animal lineages remains unexplored6,7. To trace the evolutionary origins of animal genome regulation, here we characterized the physical organization of the genome in non-bilaterian animals (sponges, ctenophores, placozoans and cnidarians)8,9 and their closest unicellular relatives (ichthyosporeans, filastereans and choanoflagellates)10 by combining high-resolution chromosome conformation capture11,12 with epigenomic marks and gene expression data. Our comparative analysis showed that chromatin looping is a conserved feature of genome architecture in ctenophores, placozoans and cnidarians. These sequence-determined distal contacts involve both promoter-enhancer and promoter-promoter interactions. By contrast, chromatin loops are absent in the unicellular relatives of animals. Our findings indicate that spatial genome regulation emerged early in animal evolution. This evolutionary innovation introduced regulatory complexity, ultimately facilitating the diversification of animal developmental programmes and cell type repertoires.

The homie insulator has sub-elements with different insulating and long-range pairing properties

Chromatin insulators are major determinants of chromosome architecture. Specific architectures induced by insulators profoundly influence nuclear processes, including how enhancers and promoters interact over long distances and between homologous chromosomes. Insulators can pair with copies of themselves in trans to facilitate homolog pairing. They can also pair with other insulators, sometimes with great specificity, inducing long-range chromosomal loops. Contrary to their canonical function of enhancer blocking, these loops can bring distant enhancers and promoters together to activate gene expression, while at the same time blocking other interactions in cis. The details of these effects depend on the choice of pairing partner, and on the orientation specificity of pairing, implicating the 3D architecture as a major functional determinant. Here, we dissect the homie insulator from the Drosophila even skipped (eve) locus, to understand its substructure. We test pairing function based on homie-carrying transgenes interacting with endogenous eve. The assay is sensitive to both pairing strength and orientation. Using this assay, we found that a Su(Hw) binding site in homie is required for efficient long-range interaction, although some activity remains without it. This binding site also contributes to the canonical insulator activities of enhancer blocking and barrier function. Based on this and other results from our functional dissection, each of the canonical insulator activities, chromosomal loop formation, enhancer blocking, and barrier activity, are partially separable. Our results show the complexity inherent in insulator functions, which can be provided by an array of different proteins with both shared and distinct properties.

pC-SAC: A method for high-resolution 3D genome reconstruction from low-resolution Hi-C data

The three-dimensional (3D) organization of the genome is crucial for gene regulation, with disruptions linked to various diseases. High-throughput Chromosome Conformation Capture (Hi-C) and related technologies have advanced our understanding of 3D genome organization by mapping interactions between distal genomic regions. However, capturing enhancer-promoter interactions at high resolution remains challenging due to the high sequencing depth required. We introduce pC-SAC (probabilistically Constrained Self-Avoiding Chromatin), a novel computational method for producing accurate high-resolution Hi-C matrices from low-resolution data. pC-SAC uses adaptive importance sampling with sequential Monte Carlo to generate ensembles of 3D chromatin chains that satisfy physical constraints derived from low-resolution Hi-C data. Our method achieves over 95% accuracy in reconstructing high-resolution chromatin maps and identifies novel interactions enriched with candidate cis-regulatory elements (cCREs) and expression quantitative trait loci (eQTLs). Benchmarking against state-of-the-art deep learning models demonstrates pC-SAC's performance in both short- and long-range interaction reconstruction. pC-SAC offers a cost-effective solution for enhancing the resolution of Hi-C data, thus enabling deeper insights into 3D genome organization and its role in gene regulation and disease. Our tool can be found at https://github.com/G2Lab/pCSAC.

3D genome folding in epigenetic regulation and cellular memory

The 3D folding of the genome is tightly linked to its epigenetic state which maintains gene expression programmes. Although the relationship between gene expression and genome organisation is highly context dependent, 3D genome organisation is emerging as a novel epigenetic layer to reinforce and stabilise transcriptional states. Whether regulatory information carried in genome folding could be transmitted through mitosis is an area of active investigation. In this review, we discuss the relationship between epigenetic state and nuclear organisation, as well as the interplay between transcriptional regulation and epigenetic genome folding. We also consider the architectural remodelling of nuclei as cells enter and exit mitosis, and evaluate the potential of the 3D genome to contribute to cellular memory.

Chromatin as a Coevolutionary Graph: Modeling the Interplay of Replication with Chromatin Dynamics

Modeling DNA replication poses significant challenges due to the intricate interplay of biophysical processes and the need for precise parameter optimization. In this study, we explore the interactions among three key biophysical factors that influence chromatin folding: replication, loop extrusion, and compartmentalization. Replication forks, known to act as barriers to the motion of loop extrusion factors, also correlate with the phase separation of chromatin into A and B compartments. Our approach integrates three components: (1) a numerical model that takes into advantage single-cell replication timing data to simulate replication fork propagation; (2) a stochastic Monte Carlo simulation that captures the interplay between the biophysical factors, with loop extrusion factors binding, unbinding, and extruding dynamically, while CTCF barriers and replication forks act as static and moving barriers, and a Potts Hamiltonian governs the spreading of epigenetic states driving chromatin compartmentalization; and (3) a 3D OpenMM simulation that reconstructs the chromatin's 3D structure based on the states generated by the stochastic model. To our knowledge, this is the first framework to dynamically integrate and simulate these three biophysical factors, enabling insights into chromatin behavior during replication. Furthermore, we investigate how replication stress alters these dynamics and affects chromatin structure.

The synergy between compartmentalization and motorization in chromatin architecture

High-resolution techniques capable of manipulating from single molecules to millions of cells are combined with three-dimensional modeling followed by simulation to comprehend the specific aspects of chromosomes. From the theoretical perspective, the energy landscape theory from protein folding inspired the development of the minimal chromatin model (MiChroM). In this work, two biologically relevant MiChroM energy terms were minimized under different conditions, revealing a competition between loci compartmentalization and motor-driven activity mechanisms in chromatin folding. Enhancing the motor activity energy baseline increased the lengthwise compaction and reduced the polymer entanglement. Concomitantly, decreasing compartmentalization-related interactions reduced the overall polymer collapse, although compartmentalization given by the microphase separation remained almost intact. For multiple chromosome simulations, increased motorization intensified the territory formation of the different chains and reduced compartmentalization strength lowered the probability of contact formation of different loci between multiple chains, approximating to the experimental inter-contacts of the human chromosomes. These findings have direct implications for experimental data-driven chromosome modeling, specially those involving multiple chromosomes. The interplay between phase-separation and territory formation mechanisms should be properly implemented in order to recover the genome architecture and dynamics, features that might play critical roles in regulating nuclear functions.

In silico nanoscope to study the interplay of genome organization and transcription regulation

In eukaryotic genomes, regulated access and communication between cis-regulatory elements (CREs) are necessary for enhancer-mediated transcription of genes. The molecular framework of the chromatin organization underlying such communication remains poorly understood. To better understand it, we develop a multiscale modeling pipeline to build near-atomistic models of the 200 kb Nanog gene locus in mouse embryonic stem cells comprising nucleosomes, transcription factors, co-activators, and RNA polymerase II-mediator complexes. By integrating diverse experimental data, including protein localization, genomic interaction frequencies, cryo-electron microscopy, and single-molecule fluorescence studies, our model offers novel insights into chromatin organization and its role in enhancer-promoter communication. The models equilibrated by high-performance molecular dynamics simulations span a scale of ∼350 nm, revealing an experimentally consistent local and global organization of chromatin and transcriptional machinery. Our models elucidate that the sequence-regulated chromatin accessibility facilitates the recruitment of transcription regulatory proteins exclusively at CREs, guided by the contrasting nucleosome organization compared to other regions. By constructing an experimentally consistent near-atomic model of chromatin in the cellular environment, our approach provides a robust framework for future studies on nuclear compartmentalization, chromatin organization, and transcription regulation.

Novel AI-powered computational method using tensor decomposition for identification of common optimal bin sizes when integrating multiple Hi-C datasets

Identifying the optimal bin sizes (or resolutions) for the integration of multiple Hi-C datasets is a challenge due to the fact that bin sizes must be common over multiple datasets. By contrast, the dependence of quality upon bin sizes can vary from dataset to dataset. Moreover, common structures should not be sought in bin sizes smaller than the optimal bin sizes, below which common structure cannot be the primary structure any more even after increasing the number of mapped short reads per bin. In this case, there are no common structures at finer resolutions, suggesting that individual Hi-C datasets may have to be analyzed separately in the bin sizes smaller than the optimal one. Thus, quality assessments of individual datasets have a limited ability to determine the best bin size for all datasets. In this study, we propose a novel application of tensor decomposition (TD) based unsupervised feature extraction (FE) to choose the optimal bin sizes for the integration of multiple Hi-C datasets. TD-based unsupervised FE exhibit phase transition-like phenomena through which the smallest possible bin size (or the highest resolution) can be automatically estimated empirically, without the need to manually set a threshold value for the integration of multiple Hi-C datasets, retrieved from GEO with GEO ID, GSE260760 and GSE255264. To our knowledge, ours is the first one that can optimize bin sizes over multiple Hi-C profiles without any tunable parameters.

Memory CD4+ T cells sequentially restructure their 3D genome during stepwise activation

Background

CD4+ T cells are a highly differentiated cell type that maintain enough transcriptomic plasticity to cycle between activated and memory statuses. How the 1D chromatin state and 3D chromatin architecture support this plasticity is under intensive investigation.

Methods

Here, we wished to test a commercially available in situ Hi-C kit (Arima Genomics Inc.) to establish whether published performance on limiting cell numbers from clonal cell lines copies across to a primary immune cell type. We achieved comparable contact matrices from 50,000, 250,000, and 1,000,000 memory CD4+ T-cell inputs. We generated multiple Hi-C and RNA-seq libraries from the same biological blood donors under three separate conditions: unstimulated fresh ex vivo, IL-2-only stimulated, and T cell receptor (TCR)+CD28+IL-2-stimulated, conferring increasingly stronger activation signals. We wished to capture the magnitude and progression of 3D chromatin shifts and correlate these to expression changes under the two stimulations.

Results

Although some genome organization changes occurred concomitantly with changes in gene expression, at least as many changes occurred without corresponding changes in expression. Counter to the hypothesis that topologically associated domains (TADs) are largely invariant structures providing a scaffold for dynamic looping contacts between enhancers and promotors, we found that there were at least as many dynamic TAD changes. Stimulation with IL-2 alone triggered many changes in genome organization, and many of these changes were strengthened by additional TCR and CD28 co-receptor stimulation.

Conclusions

This suggests a stepwise process whereby mCD4+ T cells undergo sequential buildup of 3D architecture induced by distinct or combined stimuli likely to "prime" or "deprime" them for expression responses to subsequent TCR-antigen ligation or additional cytokine stimulation.

Rethinking chromatin accessibility: from compaction to dynamic interactions

The genome is traditionally divided into condensed heterochromatin and open euchromatin. However, recent findings challenge this binary classification and the notion that chromatin condensation solely governs the accessibility of transcription factors (TFs) and, consequently, gene expression. Instead, chromatin accessibility is emerging as a factor-specific property that is influenced by multiple determinants. These include the mobility of the chromatin fiber, the capacity of TFs to engage repeatedly with it through multivalent interactions, and the four-dimensional organization of its surrounding diffusible space. Unraveling the molecular and biophysical principles that render a genomic target truly accessible remains a significant challenge, but innovative methods for locally perturbing chromatin, coupled with microscopy techniques that offer single-molecule sensitivity, provide an exciting experimental playground to test new hypotheses.

COCOA: A Framework for Fine-scale Mapping of Cell-type-specific Chromatin Compartments Using Epigenomic Information

Chromatin compartmentalization and epigenomic modifications play crucial roles in cell differentiation and disease development. However, precise mapping of chromatin compartment patterns requires Hi-C or Micro-C data at high sequencing depth. Exploring the systematic relationship between epigenomic modifications and compartment patterns remains challenging. To address these issues, we present COCOA, a deep neural network framework using convolution and attention mechanisms to infer fine-scale chromatin compartment patterns from six histone modification signals. COCOA extracts 1D track features through bidirectional feature reconstruction after resolution-specific binning of epigenomic signals. These track features are then cross-fused with contact features using an attention mechanism and transformed into chromatin compartment patterns through residual feature reduction. COCOA demonstrates accurate inference of chromatin compartmentalization at a fine-scale resolution and exhibits stable performance on test sets. Additionally, we explored the impact of histone modifications on chromatin compartmentalization prediction through in silico epigenomic perturbation experiments. Unlike obscure compartments observed in high-depth experimental data at 1-kb resolution, COCOA generates clear and detailed compartment patterns, highlighting its superior performance. Finally, we demonstrate that COCOA enables cell-type-specific prediction of unrevealed chromatin compartment patterns in various biological processes, making it an effective tool for gaining insights into chromatin compartmentalization from epigenomics in diverse biological scenarios. The COCOA Python code is publicly available at https://github.com/onlybugs/COCOA and https://ngdc.cncb.ac.cn/biocode/tools/BT007498.

Genome-wide absolute quantification of chromatin looping

3D genomics methods such as Hi-C and Micro-C have uncovered chromatin loops across the genome and linked these loops to gene regulation. However, these methods only measure 3D interaction probabilities on a relative scale. Here, we overcome this limitation by using live imaging data to calibrate Micro-C in mouse embryonic stem cells, thus obtaining absolute looping probabilities for 36,804 chromatin loops across the genome. We find that the looped state is generally rare, with a mean probability of 2.3% and a maximum of 26% across the quantified loops. On average, CTCF-CTCF loops are stronger than loops between cis-regulatory elements (3.2% vs. 1.1%). Our findings can be extended to human stem cells and differentiated cells under certain assumptions. Overall, we establish an approach for genome-wide absolute loop quantification and report that loops generally occur with low probabilities, generalizing recent live imaging results to the whole genome.

Hi-GDT: A Hi-C-based 3D gene domain analysis tool for analyzing local chromatin contacts in plants

BACKGROUND

Three-dimensional (3D) chromatin organization is emerging as a key factor in gene regulation in eukaryotes. Recent studies using high-resolution Hi-C analysis have explored fine-scale local chromatin contact domains in plants, as exemplified by the basic contact domains established at accessible gene border regions in Arabidopsis (Arabidopsis thaliana). However, we lack effective tools to identify these contact domains and examine their structural dynamics.

RESULTS

We developed the Hi-C-based 3D Gene Domain analysis Tool (Hi-GDT) to identify fine-scale local chromatin contact domains in plants, with a particular focus on gene borders. Hi-GDT successfully identifies local contact domains, including single-gene and multigene domains, with high reproducibility. Hi-GDT can also be used to discover local contact domains that are differentially organized in association with differences in gene expression between tissue types, genotypes, or in response to environmental stimuli.

CONCLUSIONS

Hi-GDT is a valuable tool for identifying genes regulated by dynamic 3D conformational changes, expanding our understanding of the structural and functional relevance of local 3D chromatin organization in plants. Hi-GDT is publicly available at https://github.com/CDL-HongwooLee/Hi-GDT.

Widespread 3D genome reorganization precedes programmed DNA rearrangement in Oxytricha trifallax

Genome organization recapitulates function, yet ciliates like Oxytricha trifallax possess highly-specialized germline genomes, which are largely transcriptionally silent. During post-zygotic development, Oxytricha's germline undergoes large-scale genome editing, rearranging precursor genome elements into a transcriptionally-active genome with thousands of gene-sized nanochromosomes. Transgenerationally-inherited RNAs, derived from the parental somatic genome, program the retention and reordering of germline fragments. Retained and eliminated DNA must be distinguished and processed separately, but the role of chromatin organization in this process is unknown. We developed tools for studying Oxytricha nuclei and apply them to map the 3D organization of precursor and developmental states using Hi-C. We find that the precursor conformation primes the germline for development, while a massive spatial reorganization during development differentiates retained from eliminated regions before DNA rearrangement. Further experiments suggest a role for RNA-DNA interactions and chromatin remodeling in this process, implying a critical role for 3D architecture in programmed genome rearrangement.

ChromatinHD connects single-cell DNA accessibility and conformation to gene expression through scale-adaptive machine learning

Gene regulation is inherently multiscale, but scale-adaptive machine learning methods that fully exploit this property in single-nucleus accessibility data are still lacking. Here, we develop ChromatinHD, a pair of scale-adaptive models that uses the raw accessibility data, without peak-calling or windows, to link regions to gene expression and determine differentially accessible chromatin. We show how ChromatinHD consistently outperforms existing peak and window-based approaches and find that this is due to a large number of uniquely captured, functional accessibility changes within and outside of putative cis-regulatory regions. Furthermore, ChromatinHD can delineate collaborating regulatory regions, including their preferential genomic conformations, that drive gene expression. Finally, our models also use changes in ATAC-seq fragment lengths to identify dense binding of transcription factors, a feature not captured by footprinting methods. Altogether, ChromatinHD, available at https://chromatinhd.org , is a suite of computational tools that enables a data-driven understanding of chromatin accessibility at various scales and how it relates to gene expression.

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.

Topoisomerase-modulated genome-wide DNA supercoiling domains colocalize with nuclear compartments and regulate human gene expression

DNA supercoiling is a biophysical feature of the double helix with a pivotal role in biological processes. However, understanding of DNA supercoiling in the chromatin remains limited. Here, we developed azide-trimethylpsoralen sequencing (ATMP-seq), a DNA supercoiling assay offering quantitative accuracy while minimizing genomic bias and background noise. Using ATMP-seq, we directly visualized transcription-dependent negative and positive twin-supercoiled domains around genes and mapped kilobase-resolution DNA supercoiling throughout the human genome. Remarkably, we discovered megabase-scale supercoiling domains (SDs) across all chromosomes that are modulated mainly by topoisomerases I and IIβ. Transcription activities, but not the consequent supercoiling accumulation in the local region, contribute to SD formation, indicating the long-range propagation of transcription-generated supercoiling. Genome-wide SDs colocalize with A/B compartments in both human and Drosophila cells but are distinct from topologically associating domains (TADs), with negative supercoiling accumulation at TAD boundaries. Furthermore, genome-wide DNA supercoiling varies between cell states and types and regulates human gene expression, underscoring the importance of supercoiling dynamics in chromatin regulation and function.

DDX18 coordinates nucleolus phase separation and nuclear organization to control the pluripotency of human embryonic stem cells

Pluripotent stem cells possess a unique nuclear architecture characterized by a larger nucleus and more open chromatin, which underpins their ability to self-renew and differentiate. Here, we show that the nucleolus-specific RNA helicase DDX18 is essential for maintaining the pluripotency of human embryonic stem cells. Using techniques such as Hi-C, DNA/RNA-FISH, and biomolecular condensate analysis, we demonstrate that DDX18 regulates nucleolus phase separation and nuclear organization by interacting with NPM1 in the granular nucleolar component, driven by specific nucleolar RNAs. Loss of DDX18 disrupts nucleolar substructures, impairing centromere clustering and perinucleolar heterochromatin (PNH) formation. To probe this further, we develop NoCasDrop, a tool enabling precise nucleolar targeting and controlled liquid condensation, which restores centromere clustering and PNH integrity while modulating developmental gene expression. This study reveals how nucleolar phase separation dynamics govern chromatin organization and cell fate, offering fresh insights into the molecular regulation of stem cell pluripotency.

An integrative TAD catalog in lymphoblastoid cell lines discloses the functional impact of deletions and insertions in human genomes

The human genome is packaged within a three-dimensional (3D) nucleus and organized into structural units known as compartments, topologically associating domains (TADs), and loops. TAD boundaries, separating adjacent TADs, have been found to be well conserved across mammalian species and more evolutionarily constrained than TADs themselves. Recent studies show that structural variants (SVs) can modify 3D genomes through the disruption of TADs, which play an essential role in insulating genes from outside regulatory elements' aberrant regulation. However, how SV affects the 3D genome structure and their association among different aspects of gene regulation and candidate cis-regulatory elements (cCREs) have rarely been studied systematically. Here, we assess the impact of SVs intersecting with TAD boundaries by developing an integrative Hi-C analysis pipeline, which enables the generation of an in-depth catalog of TADs and TAD boundaries in human lymphoblastoid cell lines (LCLs) to fill the gap of limited resources. Our catalog contains 18,865 TADs, including 4596 sub-TADs, with 185 SVs (TAD-SVs) that alter chromatin architecture. By leveraging the ENCODE registry of cCREs in humans, we determine that 34 of 185 TAD-SVs intersect with cCREs and observe significant enrichment of TAD-SVs within cCREs. This study provides a database of TADs and TAD-SVs in the human genome that will facilitate future investigations of the impact of SVs on chromatin structure and gene regulation in health and disease.

HiCrayon reveals distinct layers of multi-state 3D chromatin organization

Chromatin contact maps are often shown as 2D heatmaps and visually compared to 1D genomic data by simple juxtaposition. While common, this strategy is imprecise, placing the onus on the reader to align features with each other. To remedy this, we developed HiCrayon, an interactive tool that facilitates the integration of 3D chromatin organization maps and 1D datasets. This visualization method integrates data from genomic assays directly into the chromatin contact map by coloring interactions according to 1D signal. HiCrayon is implemented using R shiny and python to create a graphical user interface application, available in both web and containerized format to promote accessibility. We demonstrate the utility of HiCrayon in visualizing the effectiveness of compartment calling and the relationship between ChIP-seq and various features of chromatin organization. We also demonstrate the improved visualization of other 3D genomic phenomena, such as differences between loops associated with CTCF/cohesin versus those associated with H3K27ac. We then demonstrate HiCrayon's visualization of organizational changes that occur during differentiation and use HiCrayon to detect compartment patterns that cannot be assigned to either A or B compartments, revealing a distinct third chromatin compartment.

Multiscale molecular modeling of chromatin with MultiMM: From nucleosomes to the whole genome

Motivation: We present a user-friendly 3D chromatin simulation model for the human genome based on OpenMM, addressing the challenges posed by existing models with use-specific implementations. Our approach employs a multi-scale energy minimization strategy, capturing chromatin's hierarchical structure. Initiating with a Hilbert curve-based structure, users can input files specifying nucleosome positioning, loops, compartments, or subcompartments. Results: The model utilizes an energy minimization approach with a large choice of numerical integrators, providing the entire genome's structure within minutes. Output files include the generated structures for each chromosome, offering a versatile and accessible tool for chromatin simulation in bioinformatics studies. Furthermore, MultiMM is capable of producing nucleosome-resolution structures by making simplistic geometric assumptions about the structure and the density of nucleosomes on the DNA. Code availability: Open-source software and the manual are freely available on https://github.com/SFGLab/MultiMM or via pip https://pypi.org/project/MultiMM/.

Transcriptional machinery as an architect of genome structure

Chromatin organization, facilitated by compartmentalization and loop extrusion, is crucial for proper gene expression and cell viability. Transcription has long been considered important for shaping genome architecture due to its pervasive activity across the genome and impact on the local chromatin environment. Although earlier studies suggested a minimal contribution of transcription to shaping global genome structure, recent insights from high-resolution chromatin contact mapping, polymer simulations, and acute perturbations have revealed its critical role in dynamic chromatin organization at the level of active genes and enhancer-promoter interactions. In this review, we discuss these latest advances, highlighting the direct interplay between transcriptional machinery and loop extrusion. Finally, we explore how transcription of genes and non-coding regulatory elements may contribute to the specificity of gene regulation, focusing on enhancers as sites of targeted cohesin loading.

Uncovering the functions and mechanisms of regulatory elements-associated non-coding RNAs

Over the past decade, regulatory non-coding RNAs (ncRNAs) produced by RNA Pol II have been revealed as meaningful players in various essential cellular functions. In particular, thousands of ncRNAs are produced at transcriptional regulatory elements such as enhancers and promoters, where they may exert multiple functions to regulate proper development, cellular programming, transcription or genomic stability. Here, we review the mechanisms involving these regulatory element-associated ncRNAs, and particularly enhancer RNAs (eRNAs) and PROMoter uPstream Transcripts (PROMPTs). We contextualize the mechanisms described to the processing and degradation of these short lived RNAs. We summarize recent findings explaining how ncRNAs operate locally at promoters and enhancers, or further away, either shortly after their production by RNA Pol II, or through post-transcriptional stabilization. Such discoveries lead to a converging model accounting for how ncRNAs influence cellular fate, by acting on transcription and chromatin structure, which may further involve factors participating to 3D nuclear organization.

Functional associations of evolutionarily recent human genes exhibit sensitivity to the 3D genome landscape and disease

Genome organization is intricately tied to regulating genes and associated cell fate decisions. Here, we examine the positioning and functional significance of human genes, grouped by their lineage restriction level, within the 3D organization of the genome. We reveal that genes of different lineage restriction levels have distinct positioning relationships with both domains and loop anchors, and remarkably consistent relationships with boundaries across cell types. While the functional associations of each group of genes are primarily cell type-specific, associations of conserved genes maintain greater stability across 3D genomic features and disease than recently evolved genes. Furthermore, the expression of these genes across various tissues follows an evolutionary progression, such that RNA levels increase from young lineage restricted genes to ancient genes present in most species. Thus, the distinct relationships of gene evolutionary age, function, and positioning within 3D genomic features contribute to tissue-specific gene regulation in development and disease.

Type II topoisomerases shape multi-scale 3D chromatin folding in regions of positive supercoils

Type II topoisomerases (TOP2s) resolve torsional stress accumulated during various cellular processes and are enriched at chromatin loop anchors and topologically associated domain (TAD) boundaries, where, when trapped, can lead to genomic instability promoting the formation of oncogenic fusions. Whether TOP2s relieve topological constraints at these positions and/or participate in 3D chromosome folding remains unclear. Here, we combine 3D genomics, imaging, and GapRUN, a method for the genome-wide profiling of positive supercoiling, to assess the role of TOP2s in shaping chromosome organization in human cells. Acute TOP2 depletion led to the emergence of new, large-scale contacts at the boundaries between active, positively supercoiled, and lamina-associated domains. TOP2-dependent changes at the higher-order chromatin folding were accompanied by remodeling of chromatin-nuclear lamina interactions and of gene expression, while at the chromatin loop level, TOP2 depletion predominantly remodeled transcriptionally anchored, positively supercoiled loops. We propose that TOP2s act as a fine regulator of chromosome folding at multiple scales.

The chromosome folding problem and how cells solve it

Every cell must solve the problem of how to fold its genome. We describe how the folded state of chromosomes is the result of the combined activity of multiple conserved mechanisms. Homotypic affinity-driven interactions lead to spatial partitioning of active and inactive loci. Molecular motors fold chromosomes through loop extrusion. Topological features such as supercoiling and entanglements contribute to chromosome folding and its dynamics, and tethering loci to sub-nuclear structures adds additional constraints. Dramatically diverse chromosome conformations observed throughout the cell cycle and across the tree of life can be explained through differential regulation and implementation of these basic mechanisms. We propose that the first functions of chromosome folding are to mediate genome replication, compaction, and segregation and that mechanisms of folding have subsequently been co-opted for other roles, including long-range gene regulation, in different conditions, cell types, and species.

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.

Systematic identification of interchromosomal interaction networks supports the existence of specialized RNA factories

Most studies of genome organization have focused on intrachromosomal (cis) contacts because they harbor key features such as DNA loops and topologically associating domains. Interchromosomal (trans) contacts have received much less attention, and tools for interrogating potential biologically relevant trans structures are lacking. Here, we develop a computational framework that uses Hi-C data to identify sets of loci that jointly interact in trans This method, trans-C, initiates probabilistic random walks with restarts from a set of seed loci to traverse an input Hi-C contact network, thereby identifying sets of trans-contacting loci. We validate trans-C in three increasingly complex models of established trans contacts: the Plasmodium falciparum var genes, the mouse olfactory receptor "Greek islands," and the human RBM20 cardiac splicing factory. We then apply trans-C to systematically test the hypothesis that genes coregulated by the same trans-acting element (i.e., a transcription or splicing factor) colocalize in three dimensions to form "RNA factories" that maximize the efficiency and accuracy of RNA biogenesis. We find that many loci with multiple binding sites of the same DNA-binding proteins interact with one another in trans, especially those bound by factors with intrinsically disordered domains. Similarly, clustered binding of a subset of RNA-binding proteins correlates with trans interaction of the encoding loci. We observe that these trans-interacting loci are close to nuclear speckles. These findings support the existence of trans- interacting chromatin domains (TIDs) driven by RNA biogenesis. Trans-C provides an efficient computational framework for studying these and other types of trans interactions, empowering studies of a poorly understood aspect of genome architecture.

The 3D Genome in Brain Development: An Exploration of Molecular Mechanisms and Experimental Methods

The human brain contains multiple cell types that are spatially organized into functionally distinct regions. The proper development of the brain requires complex gene regulation mechanisms in both neurons and the non-neuronal cell types that support neuronal function. Studies across the last decade have discovered that the 3D nuclear organization of the genome is instrumental in the regulation of gene expression in the diverse cell types of the brain. In this review, we describe the fundamental biochemical mechanisms that regulate the 3D genome, and comprehensively describe in vitro and ex vivo studies on mouse and human brain development that have characterized the roles of the 3D genome in gene regulation. We highlight the significance of the 3D genome in linking distal enhancers to their target promoters, which provides insights on the etiology of psychiatric and neurological disorders, as the genetic variants associated with these disorders are primarily located in noncoding regulatory regions. We also describe the molecular mechanisms that regulate chromatin folding and gene expression in neurons. Furthermore, we describe studies with an evolutionary perspective, which have investigated features that are conserved from mice to human, as well as human gained 3D chromatin features. Although most of the insights on disease and molecular mechanisms have been obtained from bulk 3C based experiments, we also highlight other approaches that have been developed recently, such as single cell 3C approaches, as well as non-3C based approaches. In our future perspectives, we highlight the gaps in our current knowledge and emphasize the need for 3D genome engineering and live cell imaging approaches to elucidate mechanisms and temporal dynamics of chromatin interactions, respectively.

LAP2alpha facilitates myogenic gene expression by preventing nucleoplasmic lamin A/C from spreading to active chromatin regions

A-type lamins form a filamentous meshwork beneath the nuclear membrane that anchors large heterochromatic genomic regions at the nuclear periphery. A-type lamins also exist as a dynamic, non-filamentous pool in the nuclear interior, where they interact with lamin-associated polypeptide 2 alpha (LAP2α). Both proteins associate with largely overlapping euchromatic genomic regions in the nucleoplasm, but the functional significance of this interaction is poorly understood. Here, we report that LAP2α relocates towards regions containing myogenic genes in the early stages of muscle differentiation, possibly facilitating efficient gene regulation, while lamins A and C mostly associate with genomic regions away from these genes. Strikingly, upon depletion of LAP2α, A-type lamins spread across active chromatin and accumulate at regions of active H3K27ac and H3K4me3 histone marks in the vicinity of myogenic genes whose expression is impaired in the absence of LAP2α. Reorganization of A-type lamins on chromatin is accompanied by depletion of the active chromatin mark H3K27ac and a significantly impaired myogenic differentiation. Thus, the interplay of LAP2α and A-type lamins is crucial for proper positioning of intranuclear lamin A/C on chromatin to allow efficient myogenic differentiation.

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.

Genome organization and stability in mammalian pre-implantation development

A largely stable genome is required for normal development, even as genetic change is an integral aspect of reproduction, genetic adaptation and evolution. Recent studies highlight a critical window of mammalian development with intrinsic DNA replication stress and genome instability in the first cell divisions after fertilization. Patterns of DNA replication and genome stability are established very early in mammals, alongside patterns of nuclear organization, and before the emergence of gene expression patterns, and prior to cell specification and germline formation. The study of DNA replication and genome stability in the mammalian embryo provides a unique cellular system due to the resetting of the epigenome to a totipotent state, and the de novo establishment of the patterns of nuclear organization, gene expression, DNA methylation, histone modifications and DNA replication. Studies on DNA replication and genome stability in the early mammalian embryo is relevant for understanding both normal and disease-causing genetic variation, and to uncover basic principles of genome regulation.

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.

Genome structural dynamics: insights from Gaussian network analysis of Hi-C data

Characterization of the spatiotemporal properties of the chromatin is essential to gaining insights into the physical bases of gene co-expression, transcriptional regulation and epigenetic modifications. The Gaussian network model (GNM) has proven in recent work to serve as a useful tool for modeling chromatin structural dynamics, using as input high-throughput chromosome conformation capture data. We focus here on the exploration of the collective dynamics of chromosomal structures at hierarchical levels of resolution, from single gene loci to topologically associating domains or entire chromosomes. The GNM permits us to identify long-range interactions between gene loci, shedding light on the role of cross-correlations between distal regions of the chromosomes in regulating gene expression. Notably, GNM analysis performed across diverse cell lines highlights the conservation of the global/cooperative movements of the chromatin across different types of cells. Variations driven by localized couplings between genomic loci, on the other hand, underlie cell differentiation, underscoring the significance of the four-dimensional properties of the genome in defining cellular identity. Finally, we demonstrate the close relation between the cell type-dependent mobility profiles of gene loci and their gene expression patterns, providing a clear demonstration of the role of chromosomal 4D features in defining cell-specific differential expression of genes.

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.

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.

Loss of multi-level 3D genome organization during breast cancer progression

Breast cancer entails intricate alterations in genome organization and expression. However, how three-dimensional (3D) chromatin structure changes in the progression from a normal to a breast cancer malignant state remains unknown. To address this, we conducted an analysis combining Hi-C data with lamina-associated domains (LADs), epigenomic marks, and gene expression in an in vitro model of breast cancer progression. Our results reveal that while the fundamental properties of topologically associating domains (TADs) are overall maintained, significant changes occur in the organization of compartments and subcompartments. These changes are closely correlated with alterations in the expression of oncogenic genes. We also observe a restructuring of TAD-TAD interactions, coinciding with a loss of spatial compartmentalization and radial positioning of the 3D genome. Notably, we identify a previously unrecognized interchromosomal insertion event, wherein a locus on chromosome 8 housing the MYC oncogene is inserted into a highly active subcompartment on chromosome 10. This insertion is accompanied by the formation of de novo enhancer contacts and activation of MYC, illustrating how structural genomic variants can alter the 3D genome to drive oncogenic states. In summary, our findings provide evidence for the loss of genome organization at multiple scales during breast cancer progression revealing novel relationships between genome 3D structure and oncogenic processes.

Stem-loop and circle-loop TADs generated by directional pairing of boundary elements have distinct physical and regulatory properties

The chromosomes in multicellular eukaryotes are organized into a series of topologically independent loops called TADs. In flies, TADs are formed by physical interactions between neighboring boundaries. Fly boundaries exhibit distinct partner preferences, and pairing interactions between boundaries are typically orientation-dependent. Pairing can be head-to-tail or head-to-head. The former generates a stem-loop TAD, while the latter gives a circle-loop TAD. The TAD that encompasses the Drosophila even skipped (eve) gene is formed by the head-to-tail pairing of the nhomie and homie boundaries. To explore the relationship between loop topology and the physical and regulatory landscape, we flanked the nhomie boundary region with two attP sites. The attP sites were then used to generate four boundary replacements: λ DNA, nhomie forward (WT orientation), nhomie reverse (opposite of WT orientation), and homie forward (same orientation as WT homie). The nhomie forward replacement restores the WT physical and regulatory landscape: in MicroC experiments, the eve TAD is a 'volcano' triangle topped by a plume, and the eve gene and its regulatory elements are sequestered from interactions with neighbors. The λ DNA replacement lacks boundary function: the endpoint of the 'new' eve TAD on the nhomie side is ill-defined, and eve stripe enhancers activate a nearby gene, eIF3j. While nhomie reverse and homie forward restore the eve TAD, the topology is a circle-loop, and this changes the local physical and regulatory landscape. In MicroC experiments, the eve TAD interacts with its neighbors, and the plume at the top of the eve triangle peak is converted to a pair of 'clouds' of contacts with the next-door TADs. Consistent with the loss of isolation afforded by the stem-loop topology, the eve enhancers weakly activate genes in the neighboring TADs. Conversely, eve function is partially disrupted.

Non-coding variants impact cis-regulatory coordination in a cell type-specific manner

BACKGROUND

Interactions among cis-regulatory elements (CREs) play a crucial role in gene regulation. Various approaches have been developed to map these interactions genome-wide, including those relying on interindividual epigenomic variation to identify groups of covariable regulatory elements, referred to as chromatin modules (CMs). While CM mapping allows to investigate the relationship between chromatin modularity and gene expression, the computational principles used for CM identification vary in their application and outcomes.

RESULTS

We comprehensively evaluate and streamline existing CM mapping tools and present guidelines for optimal utilization of epigenome data from a diverse population of individuals to assess regulatory coordination across the human genome. We showcase the effectiveness of our recommended practices by analyzing distinct cell types and demonstrate cell type specificity of CRE interactions in CMs and their relevance for gene expression. Integration of genotype information revealed that many non-coding disease-associated variants affect the activity of CMs in a cell type-specific manner by affecting the binding of cell type-specific transcription factors. We provide example cases that illustrate in detail how CMs can be used to deconstruct GWAS loci, assess variable expression of cell surface receptors in immune cells, and reveal how genetic variation can impact the expression of prognostic markers in chronic lymphocytic leukemia.

CONCLUSIONS

Our study presents an optimal strategy for CM mapping and reveals how CMs capture the coordination of CREs and its impact on gene expression. Non-coding genetic variants can disrupt this coordination, and we highlight how this may lead to disease predisposition in a cell type-specific manner.

Replication timing alterations are associated with mutation acquisition during breast and lung cancer evolution

During each cell cycle, the process of DNA replication timing is tightly regulated to ensure the accurate duplication of the genome. The extent and significance of alterations in this process during malignant transformation have not been extensively explored. Here, we assess the impact of altered replication timing (ART) on cancer evolution by analysing replication-timing sequencing of cancer and normal cell lines and 952 whole-genome sequenced lung and breast tumours. We find that 6%-18% of the cancer genome exhibits ART, with regions with a change from early to late replication displaying an increased mutation rate and distinct mutational signatures. Whereas regions changing from late to early replication contain genes with increased expression and present a preponderance of APOBEC3-mediated mutation clusters and associated driver mutations. We demonstrate that ART occurs relatively early during cancer evolution and that ART may have a stronger correlation with mutation acquisition than alterations in chromatin structure.

Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample

Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth's death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.

Emerging roles of nuclear bodies in genome spatial organization

Nuclear bodies (NBs) are biomolecular condensates that participate in various cellular processes and respond to cellular stimuli in the nucleus. The assembly and function of these protein- and RNA-rich bodies, such as nucleoli, nuclear speckles, and promyelocytic leukemia (PML) NBs, contribute to the spatial organization of the nucleus, regulating chromatin activities locally and globally. Recent technological advancements, including spatial multiomics approaches, have revealed novel roles of nucleoli in modulating ribosomal DNA (rDNA) and adjacent non-rDNA chromatin activity, nuclear speckles in scaffolding active genome architecture, and PML NBs in maintaining genome stability during stress conditions. In this review, we summarize emerging functions of these important NBs in the spatial organization of the genome, aided by recently developed spatial multiomics approaches toward this direction.

DNA choreography: correlating mobility and organization of DNA across different resolutions from loops to chromosomes

The dynamics of DNA in the cell nucleus plays a role in cellular processes and fates but the interplay of DNA mobility with the hierarchical levels of DNA organization is still underexplored. Here, we made use of DNA replication to directly label genomic DNA in an unbiased genome-wide manner. This was followed by live-cell time-lapse microscopy of the labeled DNA combining imaging at different resolutions levels simultaneously and allowing one to trace DNA motion across organization levels within the same cells. Quantification of the labeled DNA segments at different microscopic resolution levels revealed sizes comparable to the ones reported for DNA loops using 3D super-resolution microscopy, topologically associated domains (TAD) using 3D widefield microscopy, and also entire chromosomes. By employing advanced chromatin tracking and image registration, we discovered that DNA exhibited higher mobility at the individual loop level compared to the TAD level and even less at the chromosome level. Additionally, our findings indicate that chromatin movement, regardless of the resolution, slowed down during the S phase of the cell cycle compared to the G1/G2 phases. Furthermore, we found that a fraction of DNA loops and TADs exhibited directed movement with the majority depicting constrained movement. Our data also indicated spatial mobility differences with DNA loops and TADs at the nuclear periphery and the nuclear interior exhibiting lower velocity and radius of gyration than the intermediate locations. On the basis of these insights, we propose that there is a link between DNA mobility and its organizational structure including spatial distribution, which impacts cellular processes.

Single-molecule imaging reveals a direct role of CTCF's zinc fingers in SA interaction and cluster-dependent RNA recruitment

CTCF is a zinc finger protein associated with transcription regulation that also acts as a barrier factor for topologically associated domains (TADs) generated by cohesin via loop extrusion. These processes require different properties of CTCF-DNA interaction, and it is still unclear how CTCF's structural features may modulate its diverse roles. Here, we employ single-molecule imaging to study both full-length CTCF and truncation mutants. We show that CTCF enriches at CTCF binding sites (CBSs), displaying a longer lifetime than observed previously. We demonstrate that the zinc finger domains mediate CTCF clustering and that clustering enables RNA recruitment, possibly creating a scaffold for interaction with RNA-binding proteins like cohesin's subunit SA. We further reveal a direct recruitment and an increase of SA residence time by CTCF bound at CBSs, suggesting that CTCF-SA interactions are crucial for cohesin stability on chromatin at TAD borders. Furthermore, we establish a single-molecule T7 transcription assay and show that although a transcribing polymerase can remove CTCF from CBSs, transcription is impaired. Our study shows that context-dependent nucleic acid binding determines the multifaceted CTCF roles in genome organization and transcription regulation.

Multiscale Bayesian simulations reveal functional chromatin condensation of gene loci

Chromatin, the complex assembly of DNA and associated proteins, plays a pivotal role in orchestrating various genomic functions. To aid our understanding of the principles underlying chromatin organization, we introduce Hi-C metainference, a Bayesian approach that integrates Hi-C contact frequencies into multiscale prior models of chromatin. This approach combines both bottom-up (the physics-based prior) and top-down (the data-driven posterior) strategies to characterize the 3D organization of a target genomic locus. We first demonstrate the capability of this method to accurately reconstruct the structural ensemble and the dynamics of a system from contact information. We then apply the approach to investigate the Sox2, Pou5f1, and Nanog loci of mouse embryonic stem cells using a bottom-up chromatin model at 1 kb resolution. We observe that the studied loci are conformationally heterogeneous and organized as crumpled globules, favoring contacts between distant enhancers and promoters. Using nucleosome-resolution simulations, we then reveal how the Nanog gene is functionally organized across the multiple scales of chromatin. At the local level, we identify diverse tetranucleosome folding motifs with a characteristic distribution along the genome, predominantly open at cis-regulatory elements and compact in between. At the larger scale, we find that enhancer-promoter contacts are driven by the transient condensation of chromatin into compact domains stabilized by extensive internucleosome interactions. Overall, this work highlights the condensed, but dynamic nature of chromatin in vivo, contributing to a deeper understanding of gene structure-function relationships.

Spatial orchestration of the genome: topological reorganisation during X-chromosome inactivation

Genomes are organised through hierarchical structures, ranging from local kilobase-scale cis-regulatory contacts to large chromosome territories. Most notably, (sub)-compartments partition chromosomes according to transcriptional activity, while topologically associating domains (TADs) define cis-regulatory landscapes. The inactive X chromosome in mammals has provided unique insights into the regulation and function of the three-dimensional (3D) genome. Concurrent with silencing of the majority of genes and major alterations of its chromatin state, the X chromosome undergoes profound spatial rearrangements at multiple scales. These include the emergence of megadomains, alterations of the compartment structure and loss of the majority of TADs. Moreover, the Xist locus, which orchestrates X-chromosome inactivation, has provided key insights into regulation and function of regulatory domains. This review provides an overview of recent insights into the control of these structural rearrangements and contextualises them within a broader understanding of 3D genome organisation.

Dynamic microenvironments shape nuclear organization and gene expression

Live imaging has revealed that the regulation of gene expression is largely driven by transient interactions. For example, many regulatory proteins bind chromatin for just seconds, and loop-like genomic contacts are rare and last only minutes. These discoveries have been difficult to reconcile with our canonical models that are predicated on stable and hierarchical interactions. Proteomic microenvironments that concentrate nuclear factors may explain how brief interactions can still mediate gene regulation by creating conditions where reactions occur more frequently. Here, we summarize new imaging technologies and recent discoveries implicating microenvironments as a potential driver of nuclear function. Finally, we propose that key properties of proteomic microenvironments, such as their size, enrichment, and lifetimes, are directly linked to regulatory function.