Cell-type specialization is encoded by specific chromatin topologies.

The three-dimensional (3D) structure of chromatin is intrinsically associated with gene regulation and cell function1-3. Methods based on chromatin conformation capture have mapped chromatin structures in neuronal systems such as in vitro differentiated neurons, neurons isolated through fluorescence-activated cell sorting from cortical tissues pooled from different animals and from dissociated whole hippocampi4-6. However, changes in chromatin organization captured by imaging, such as the relocation of Bdnf away from the nuclear periphery after activation7, are invisible with such approaches8. Here we developed immunoGAM, an extension of genome architecture mapping (GAM)2,9, to map 3D chromatin topology genome-wide in specific brain cell types, without tissue disruption, from single animals. GAM is a ligation-free technology that maps genome topology by sequencing the DNA content from thin (about 220 nm) nuclear cryosections. Chromatin interactions are identified from the increased probability of co-segregation of contacting loci across a collection of nuclear slices. ImmunoGAM expands the scope of GAM to enable the selection of specific cell types using low cell numbers (approximately 1,000 cells) within a complex tissue and avoids tissue dissociation2,10. We report cell-type specialized 3D chromatin structures at multiple genomic scales that relate to patterns of gene expression. We discover extensive 'melting' of long genes when they are highly expressed and/or have high chromatin accessibility. The contacts most specific of neuron subtypes contain genes associated with specialized processes, such as addiction and synaptic plasticity, which harbour putative binding sites for neuronal transcription factors within accessible chromatin regions. Moreover, sensory receptor genes are preferentially found in heterochromatic compartments in brain cells, which establish strong contacts across tens of megabases. Our results demonstrate that highly specific chromatin conformations in brain cells are tightly related to gene regulation mechanisms and specialized functions.

Comparison of the Hi-C, GAM and SPRITE methods using polymer models of chromatin.

Hi-C, split-pool recognition of interactions by tag extension (SPRITE) and genome architecture mapping (GAM) are powerful technologies utilized to probe chromatin interactions genome wide, but how faithfully they capture three-dimensional (3D) contacts and how they perform relative to each other is unclear, as no benchmark exists. Here, we compare these methods in silico in a simplified, yet controlled, framework against known 3D structures of polymer models of murine and human loci, which can recapitulate Hi-C, GAM and SPRITE experiments and multiplexed fluorescence in situ hybridization (FISH) single-molecule conformations. We find that in silico Hi-C, GAM and SPRITE bulk data are faithful to the reference 3D structures whereas single-cell data reflect strong variability among single molecules. The minimal number of cells required in replicate experiments to return statistically similar contacts is different across the technologies, being lowest in SPRITE and highest in GAM under the same conditions. Noise-to-signal levels follow an inverse power law with detection efficiency and grow with genomic distance differently among the three methods, being lowest in GAM for genomic separations >1 Mb.

Polymer models are a versatile tool to study chromatin 3D organization.

The development of new experimental technologies is opening the way to a deeper investigation of the three-dimensional organization of chromosomes inside the cell nucleus. Genome architecture is linked to vital functional purposes, yet a full comprehension of the mechanisms behind DNA folding is still far from being accomplished. Theoretical approaches based on polymer physics have been employed to understand the complexity of chromatin architecture data and to unveil the basic mechanisms shaping its structure. Here, we review some recent advances in the field to discuss how Polymer Physics, combined with numerical Molecular Dynamics simulation and Machine Learning based inference, can capture important aspects of genome organization, including the description of tissue-specific structural rearrangements, the detection of novel, regulatory-linked architectural elements and the structural variability of chromatin at the single-cell level.

Mega-Enhancer Bodies Organize Neuronal Long Genes in the Cerebellum.

Dynamic regulation of gene expression plays a key role in establishing the diverse neuronal cell types in the brain. Recent findings in genome biology suggest that three-dimensional (3D) genome organization has important, but mechanistically poorly understood functions in gene transcription. Beyond local genomic interactions between promoters and enhancers, we find that cerebellar granule neurons undergoing differentiation in vivo exhibit striking increases in long-distance genomic interactions between transcriptionally active genomic loci, which are separated by tens of megabases within a chromosome or located on different chromosomes. Among these interactions, we identify a nuclear subcompartment enriched for near-megabase long enhancers and their associated neuronal long genes encoding synaptic or signaling proteins. Neuronal long genes are differentially recruited to this enhancer-dense subcompartment to help shape the transcriptional identities of granule neuron subtypes in the cerebellum. SPRITE analyses of higher-order genomic interactions, together with IGM-based 3D genome modeling and imaging approaches, reveal that the enhancer-dense subcompartment forms prominent nuclear structures, which we term mega-enhancer bodies. These novel nuclear bodies reside in the nuclear periphery, away from other transcriptionally active structures, including nuclear speckles located in the nuclear interior. Together, our findings define additional layers of higher-order 3D genome organization closely linked to neuronal maturation and identity in the brain.

A Polymer Physics Model to Dissect Genome Organization in Healthy and Pathological Phenotypes.

Novel technologies revealed a nontrivial spatial organization of the chromosomes within the cell nucleus, which includes different levels of compartmentalization and architectural patterns. Notably, such complex three-dimensional structure plays a crucial role in vital biological functions and its alterations can produce extensive rewiring of genomic regulatory regions, thus leading to gene misexpression and disease. Here, we show that theoretical and computational approaches, based on polymer physics, can be employed to dissect chromatin contacts in three-dimensional space and to predict the effects of pathogenic structural variants on the genome architecture. In particular, we discuss the folding of the human EPHA4 and the murine Pitx1 loci as case studies.

Physical mechanisms of chromatin spatial organization.

In higher eukaryotes, chromosomes have a complex three-dimensional (3D) conformation in the cell nucleus serving vital functional purposes, yet their folding principles remain poorly understood at the single-molecule level. Here, we summarize recent approaches from polymer physics to comprehend the physical mechanisms underlying chromatin architecture. In particular, we focus on two models that have been supported by recent, growing experimental evidence, the Loop Extrusion model and the Strings&Binders phase separation model. We discuss their key ingredients, how they compare to experimental data and some insight they provide on chromatin architecture and gene regulation. Progress in that research field are opening the possibility to predict how genomic mutations alter the network of contacts between genes and their regulators and how that is linked to genetic diseases, such as congenital disorders and cancer.

Dynamic and equilibrium properties of finite-size polymer models of chromosome folding.

Novel technologies are revealing that chromosomes have a complex three-dimensional organization within the cell nucleus that serves functional purposes. Models from polymer physics have been developed to quantitively understand the molecular principles controlling their structure and folding mechanisms. Here, by using massive molecular-dynamics simulations we show that classical scaling laws combined with finite-size effects of a simple polymer model can effectively explain the scaling behavior that chromatin exhibits at the topologically associating domains level, as revealed by experimental observations. Model results are then validated against recently published high-resolution in situ Hi-C data.

Computational approaches from polymer physics to investigate chromatin folding.

Microscopy and sequencing-based technologies are providing increasing insights into chromatin architecture. Nevertheless, a full comprehension of chromosome folding and its link with vital cell functions is far from accomplished at the molecular level. Recent theoretical and computational approaches are providing important support to experiments to dissect the three-dimensional structure of chromosomes and its organizational mechanisms. Here, we review, in particular, the String&Binders polymer model of chromatin that describes the textbook scenario where contacts between distal DNA sites are established by cognate binders. It has been shown to recapitulate key features of chromosome folding and to be able at predicting how phenotypes causing structural variants rewire the interactions between genes and regulators.

Higher-order Chromosome Structures Investigated by Polymer Physics in Cellular Morphogenesis and Differentiation.

Experimental advances in Molecular Biology demonstrated that chromatin architecture and gene regulation are deeply related. Hi-C data, for instance, returned a scenario where chromosomes form a complex pattern of interactions, including TADs, metaTADs, and compartments, correlated with genomic and epigenomic features. Here, we discuss the emerging hierarchical organization of chromatin and show how it remains partially conserved during mouse neuronal differentiation with changes highly related to modifications in gene expression. In this scenario, models of polymer physics, such as the Strings & Binders (SBS) model, can be a crucial instrument to understand the molecular mechanisms underlying the formation of such a higher order 3D structure. In particular, we focus on the case study of the murine Pitx1 genomic region. At this locus, two alternative spatial conformations take place in the hindlimb and forelimb tissues, corresponding to two different transcriptional states of Pitx1. We finally show how the structural variants can affect the locus 3D organization leading to ectopic gene expression and limb malformations.

A Dynamic Folded Hairpin Conformation Is Associated with α-Globin Activation in Erythroid Cells.

We investigate the three-dimensional (3D) conformations of the α-globin locus at the single-allele level in murine embryonic stem cells (ESCs) and erythroid cells, combining polymer physics models and high-resolution Capture-C data. Model predictions are validated against independent fluorescence in situ hybridization (FISH) data measuring pairwise distances, and Tri-C data identifying three-way contacts. The architecture is rearranged during the transition from ESCs to erythroid cells, associated with the activation of the globin genes. We find that in ESCs, the spatial organization conforms to a highly intermingled 3D structure involving non-specific contacts, whereas in erythroid cells the α-globin genes and their enhancers form a self-contained domain, arranged in a folded hairpin conformation, separated from intermingling flanking regions by a thermodynamic mechanism of micro-phase separation. The flanking regions are rich in convergent CTCF sites, which only marginally participate in the erythroid-specific gene-enhancer contacts, suggesting that beyond the interaction of CTCF sites, multiple molecular mechanisms cooperate to form an interacting domain.