Search

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., A new technique called immunoGAM, which combines genome architecture mapping (GAM) with immunoselection, enabled the discovery of specialized chromatin conformations linked to gene expression in specific cell populations from mouse brain tissues.

Loop-extrusion and phase-separation have been proposed as mechanisms that shape chromosome large-scale spatial organization. It is unclear, however, how they perform relative to each other in explaining chromatin architecture data and whether they compete or co-exist at the single-molecule level. Here, we compare models of polymer physics based on loop-extrusion and phase-separation, as well as models where both mechanisms act simultaneously in a single molecule, against multiplexed FISH data available in human loci in IMR90 and HCT116 cells. We find that the different models recapitulate bulk Hi-C and average microscopy data. Single-molecule chromatin conformations are also well captured, especially by phase-separation based models that better reflect the experimentally reported segregation in globules of the considered genomic loci and their cell-to-cell structural variability. Such a variability is consistent with two main concurrent causes: single-cell epigenetic heterogeneity and an intrinsic thermodynamic conformational degeneracy of folding. Overall, the model combining loop-extrusion and polymer phase-separation provides a very good description of the data, particularly higher-order contacts, showing that the two mechanisms can co-exist in shaping chromatin architecture in single cells.

The mammalian genome has a complex, functional 3D organization. However, it remains largely unknown how DNA contacts are orchestrated by chromatin organizers. Here, we infer from only Hi-C the cell-type-specific arrangement of DNA binding sites sufficient to recapitulate, through polymer physics, contact patterns genome wide. Our model is validated by its predictions in a set of duplications at Sox9 against available independent data. The binding site types fall in classes that well match chromatin states from segmentation studies, yet they have an overlapping, combinatorial organization along chromosomes necessary to accurately explain contact specificity. The chromatin signatures of the binding site types return a code linking chromatin states to 3D architecture. The code is validated by extensive de novo predictions of Hi-C maps in an independent set of chromosomes. Overall, our results shed light on how 3D information is encrypted in 1D chromatin via the specific combinatorial arrangement of binding sites.

The combination of modelling and experimental advances can provide deep insights for understanding chromatin 3D organization and ultimately its underlying mechanisms. In particular, models of polymer physics can help comprehend the complexity of genomic contact maps, as those emerging from technologies such as Hi-C, GAM or SPRITE. Here we discuss a method to reconstruct 3D structures from Genome Architecture Mapping (GAM) data, based on PRISMR, a computational approach introduced to find the minimal polymer model best describing Hi-C input data from only polymer physics. After recapitulating the PRISMR procedure, we describe how we extended it for treating GAM data. We successfully test the method on a 6 Mb region around the Sox9 gene and, at a lower resolution, on the whole chromosome 7 in mouse embryonic stem cells. The PRISMR derived 3D structures from GAM co-segregation data are finally validated against independent Hi-C contact maps. The method results to be versatile and robust, hinting that it can be similarly applied to different experimental data, such as SPRITE or microscopy distance data.

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.

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.

Structural variants (SVs) can result in changes in gene expression due to abnormal chromatin folding and cause disease. However, the prediction of such effects remains a challenge. Here we present a polymer-physics-based approach (PRISMR) to model 3D chromatin folding and to predict enhancer-promoter contacts. PRISMR predicts higher-order chromatin structure from genome-wide chromosome conformation capture (Hi-C) data. Using the EPHA4 locus as a model, the effects of pathogenic SVs are predicted in silico and compared to Hi-C data generated from mouse limb buds and patient-derived fibroblasts. PRISMR deconvolves the folding complexity of the EPHA4 locus and identifies SV-induced ectopic contacts and alterations of 3D genome organization in homozygous or heterozygous states. We show that SVs can reconfigure topologically associating domains, thereby producing extensive rewiring of regulatory interactions and causing disease by gene misexpression. PRISMR can be used to predict interactions in silico, thereby providing a tool for analyzing the disease-causing potential of SVs.

New Hi-C technologies have revealed that chromosomes have a complex network of spatial contacts in the cell nucleus of higher organisms, whose organisation is only partially understood. Here, we investigate the structure of such a network in human GM12878 cells, to derive a large scale picture of nuclear architecture. We find that the intensity of intra-chromosomal interactions is power-law distributed. Inter-chromosomal interactions are two orders of magnitude weaker and exponentially distributed, yet they are not randomly arranged along the genomic sequence. Intra-chromosomal contacts broadly occur between epigenomically homologous regions, whereas inter-chromosomal contacts are especially associated with regions rich in highly expressed genes. Overall, genomic contacts in the nucleus appear to be structured as a network of networks where a set of strongly individual chromosomal units, as envisaged in the 'chromosomal territory' scenario derived from microscopy, interact with each other via on average weaker, yet far from random and functionally important interactions.