Search

The coil-globule transition of hetero-polymer chains is studied here. By means of extensive Molecular Dynamics simulations, we show that the transition is directly linked to the complexity of the chain, which depends on the number of chemical species defined in the environment and the location of the binding sites along the polymer. In addition, when the number of species increases, we find that the distribution of binding sites plays an important role in triggering the transition, beyond the standard control parameters of the polymer model, i.e. binders concentration and binding affinity. Overall, our results show that by increasing the system complexity new organizational layers emerge, thus allowing a more structured control on the polymer thermodynamic state. This can be potentially applied to the study of chromatin architecture, as such polymer models have been broadly used to understand the molecular mechanisms of genome folding.

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.

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. Progresses 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.

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.

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.

SummaryHow interphase chromatin compacts into mitotic chromosomes has eluded researchers for over a century. Here we show that PHF13, a tightly regulated H3K4me epigenetic reader, and cohesin, a mediator of chromatin architecture, cooperate to drive polymer-polymer-phase-separation (PPPS) and higher order compaction of chromatin into chromosomes. PHF13 interacts with cohesin, shows similar dynamics during mitosis, and their co-depletion dramatically impairs mitotic condensation. Mechanistically, we demonstrate that PHF13 stabilizes cohesin chromatin interactions and that itself oligomerizes, resulting in a polymer with increased chromatin avidity and the ability to bridge neighboring and distant nucleosomes. Consistently, molecular dynamic simulations modelling the ability of PHF13 and cohesin to drive chromatin phase separations recapitulated our in vivo observations and are in line with 2-step condensation models.

8-Oxo-7,8-dihydro-2′-deoxyguanosine (8-oxodG), a major product of the DNA oxidization process, has been proposed to have an epigenetic function in gene regulation and has been associated with genome instability. NGS-based methodologies are contributing to the characterization of the 8-oxodG function in the genome. However, the 8-oxodG epigenetic role at a genomic level and the mechanisms controlling the genomic 8-oxodG accumulation/maintenance have not yet been fully characterized. In this study, we report the identification and characterization of a set of enhancer regions accumulating 8-oxodG in human epithelial cells. We found that these oxidized enhancers are mainly super-enhancers and are associated with bidirectional-transcribed enhancer RNAs and DNA Damage Response activation. Moreover, using ChIA-PET and HiC data, we identified specific CTCF-mediated chromatin loops in which the oxidized enhancer and promoter regions physically associate. Oxidized enhancers and their associated chromatin loops accumulate endogenous double-strand breaks which are in turn repaired by NHEJ pathway through a transcription-dependent mechanism. Our work suggests that 8-oxodG accumulation in enhancers–promoters pairs occurs in a transcription-dependent manner and provides novel mechanistic insights on the intrinsic fragility of chromatin loops containing oxidized enhancers-promoters interactions.

The development of novel experimental technologies able to map genome-wide chromatin contacts, as Hi-C, GAM or SPRITE, allowed to derive detailed information about the spatial structure of chromosomes in the cell nucleus. They revealed that the genome has a complex spatial organisation, which is highly connected with its activity. In the last years, such an abundance of experimental data prompted the development of quantitative models based on Polymer Physics to describe the chromatin architecture, clarifying many aspects about the molecular mechanisms underlying genome folding. Efficient algorithms are thus fundamental to perform massive numerical simulations for testing the accuracy of these models and provide a good description for small genomic regions or for whole chromosomes. Here, we consider the performances of Molecular Dynamics (MD) implementation of commonly used polymer physics models. Such models can be combined with Machine Learning approaches informed with experimental data to produce more accurate descriptions of real genomic regions. However, the execution times increase as a power-law with the size of the input data, which ultimately reflects the complexity of the investigated system. The best strategy is therefore a convenient trade-off between the accuracy in the description and the availability of computational resources. The combination of innovative experimental data and polymer physics theories allow to reconstruct the 3D genome structure. This is achieved by the use of machine learning approaches and massive parallel computing. Efficient algorithms and computational resources are then fundamental to produce models of increasingly high accuracy.

Cooks syndrome (CS) is an ultrarare limb malformation due to in tandem microduplications involving KCNJ2 and extending to the 5' regulatory element of SOX9. To date, six CS families were resolved at the molecular level. Subsequent studies explored the evolutionary and pathological complexities of the SOX9-KCNJ2/Sox9-Kcnj2 locus, and suggested a key role for the formation of novel topologically associating domain (TAD) by inter-TAD duplications in causing CS. Here, we report a unique case of CS associated with a de novo 1;17 translocation affecting the KCNJ2 locus. On chromosome 17, the breakpoint mapped between KCNJ16 and KCNJ2, and combined with a ~ 5 kb deletion in the 5' of KCNJ2. Based on available capture Hi-C data, the breakpoint on chromosome 17 separated KCNJ2 from a putative enhancer. Gene expression analysis demonstrated downregulation of KCNJ2 in both patient's blood cells and cultured skin fibroblasts. Our findings suggest that a complex rearrangement falling in the 5' of KCNJ2 may mimic the developmental consequences of in tandem duplications affecting the SOX9-KCNJ2/Sox9-Kcnj2 locus. This finding adds weight to the notion of an intricate role of gene regulatory regions and, presumably, the related three-dimensional chromatin structure in normal and abnormal human morphology.

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.

Regulatory landscapes drive complex developmental gene expression, but it remains unclear how their integrity is maintained when incorporating novel genes and functions during evolution. Here, we investigated how a placental mammal-specific gene, Zfp42, emerged in an ancient vertebrate topologically associated domain (TAD) without adopting or disrupting the conserved expression of its gene, Fat1. In ESCs, physical TAD partitioning separates Zfp42 and Fat1 with distinct local enhancers that drive their independent expression. This separation is driven by chromatin activity and not CTCF/cohesin. In contrast, in embryonic limbs, inactive Zfp42 shares Fat1’s intact TAD without responding to active Fat1 enhancers. However, neither Fat1 enhancer-incompatibility nor nuclear envelope-attachment account for Zfp42’s unresponsiveness. Rather, Zfp42’s promoter is rendered inert to enhancers by context-dependent DNA methylation. Thus, diverse mechanisms enabled the integration of independent Zfp42 regulation in the Fat1 locus. Critically, such regulatory complexity appears common in evolution as, genome wide, most TADs contain multiple independently expressed genes., We thank the Montpellier Ressources Imagerie facility (BioCampus Montpellier, Centre National de la Recherche Scientifique [CNRS], INSERM, University of Montpellier) and for computer resources from CINECA (ISCRA grant thanks to computer resources from INFN and CINECA [ISCRA Grant HP10C8JWU7]). G.C., Q.S., and F.B. were supported by a grant from the European Research Council (Advanced Grant 3DEpi, 788972) and by the CNRS. This work was funded by EMBO and the Wellcome Trust (ALTF1554-2016 and 206475/Z/17/Z; to M.I.R.) as well as the Deutsche Forschungsgemeinschaft (KR3985/7-3 and MU 880/16-1 to S.M.).

SUMMARYCohesin loop extrusion facilitates precise gene expression by continuously driving promoters to sample all enhancers located within the same topologically-associated domain (TAD). However, many TADs contain multiple genes with divergent expression patterns, thereby indicating additional forces further refine how enhancer activities are utilised. Here, we unravel the mechanisms enabling a new gene, Rex1, to emerge with divergent expression within the ancient Fat1 TAD in placental mammals. We show that such divergent expression is not determined by a strict enhancer-promoter compatibility code, intra-TAD position or nuclear envelope-attachment. Instead, TAD-restructuring in embryonic stem cells (ESCs) separates Rex1 and Fat1 with distinct proximal enhancers that independently drive their expression. By contrast, in later embryonic tissues, DNA methylation renders the inactive Rex1 promoter profoundly unresponsive to Fat1 enhancers within the intact TAD. Combined, these features adapted an ancient regulatory landscape during evolution to support two entirely independent Rex1 and Fat1 expression programs. Thus, rather than operating only as rigid blocks of co-regulated genes, TAD-regulatory landscapes can orchestrate complex divergent expression patterns in evolution.HIGHLIGHTSNew genes can emerge in evolution without taking on the expression pattern of their surrounding pre-existing TAD.Compartmentalisation can restructure seemingly evolutionarily stable TADs to control a promoter’s access to enhancers.Lamina-associated domains neither prevent transcriptional activation nor enhancer-promoter communication.Repression rather than promoter-specificity refines when genes respond to promiscuous enhancer activities in specific tissues.

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.

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.

The mammalian genome has a complex 3D organization, serving vital functional purposes, yet it remains largely unknown how the multitude of specific DNA contacts, e.g., between transcribed and regulatory regions, is orchestrated by chromatin organizers, such as Transcription Factors. Here, we implement a method combining machine learning and polymer physics to infer from only Hi-C data the genomic 1D arrangement of the minimal set of binding sites sufficient to recapitulate, through only physics, 3D contact patterns genome-wide in human and mouse cells. The inferred binding sites are validated by their predictions on how chromatin refolds in a set of duplications at the Sox9 locus against available independent cHi-C data, showing that their different phenotypes originate from distinct enhancer hijackings in their 3D structure. Albeit derived from only Hi-C, our binding sites fall in epigenetic classes that well match chromatin states from epigenetic segmentation studies, such as active, poised and repressed states. However, the inferred binding domains have an overlapping, combinatorial organization along chromosomes, missing in epigenetic segmentations, which is required to explain Hi-C contact specificity with high accuracy. In a reverse approach, the epigenetic profile of binding domains provides a code to derive from only epigenetic marks the DNA binding sites and, hence, the 3D architecture, as validated by successful predictions of Hi-C matrices in an independent set of chromosomes. Overall, our results shed light on how complex 3D architectural information is encrypted in 1D epigenetics via the related, combinatorial arrangement of specific binding sites along the genome.

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. This Analysis reports a computational approach to implement Hi-C, SPRITE and GAM, which allows researchers to assess the performances of the three technologies to capture DNA contacts in chromatin three-dimensional models.

Understanding the mechanisms driving the folding of chromosomes in nuclei is a major goal of modern Molecular Biology. Recent technological advances in microscopy (FISH, STORM) and sequencing approaches (Hi-C, GAM, SPRITE) enabled to collect quantitative data about chromatin 3D architecture, revealing a non-random and highly specific organization. To transform such tremendous amount of data into valuable insights on genome folding, heavy computational analyses are required. Here, we study the performances of PRISMR, a computational tool based on Machine Learning strategies and Polymer Physics principles, to explore genome 3D structure from Genome Architecture Mapping (GAM) data. Using such data, we show that PRISMR can successfully reconstruct the 3D structure of real genomic regions at various length scales, from mega-base sized loci to whole chromosomes. Importantly, the inferred structures are validated against independent Hi-C data. Finally, we show how PRISMR can be effectively employed to explore differences between experimental methods.

Insulators play a critical role in spatiotemporal gene regulation in animals. The evolutionarily conserved CCCTC-binding factor (CTCF) is required for insulator function in mammals, but not all of its binding sites act as insulators. Here we explore the sequence requirements of CTCF-mediated transcriptional insulation using a sensitive insulator reporter in mouse embryonic stem cells. We find that insulation potency depends on the number of CTCF-binding sites in tandem. Furthermore, CTCF-mediated insulation is dependent on upstream flanking sequences at its binding sites. CTCF-binding sites at topologically associating domain boundaries are more likely to function as insulators than those outside topologically associating domain boundaries, independently of binding strength. We demonstrate that insulators form local chromatin domain boundaries and weaken enhancer-promoter contacts. Taken together, our results provide genetic, molecular and structural evidence connecting chromatin topology to the action of insulators in the mammalian genome.

SummaryTechnologies for measuring 3D genome topology are increasingly important for studying mechanisms of gene regulation, for genome assembly and for mapping of genome rearrangements. Hi-C and other ligation-based methods have become routine but have specific biases. Here, we develop multiplex-GAM, a faster and more affordable version of Genome Architecture Mapping (GAM), a ligation-free technique to map chromatin contacts genome-wide. We perform a detailed comparison of contacts obtained by multiplex-GAM and Hi-C using mouse embryonic stem (mES) cells. We find that both methods detect similar topologically associating domains (TADs). However, when examining the strongest contacts detected by either method, we find that only one third of these are shared. The strongest contacts specifically found in GAM often involve “active” regions, including many transcribed genes and super-enhancers, whereas in Hi-C they more often contain “inactive” regions. Our work shows that active genomic regions are involved in extensive complex contacts that currently go under-estimated in genome-wide ligation-based approaches, and highlights the need for orthogonal advances in genome-wide contact mapping technologies.

Insulators play a critical role in spatiotemporal gene expression in metazoans by separating active and repressive chromatin domains and preventing inappropriate enhancer-promoter contacts. The evolutionarily conserved CCCTC-binding factor (CTCF) is required for insulator function in mammals, but not all of its binding sites act as insulators. Here, we explore the sequence requirements of CTCF-mediated transcriptional insulation with the use of a sensitive insulator reporter assay in mouse embryonic stem cells. We find that insulation potency depends on the number of CTCF binding sites in tandem. Furthermore, CTCF-mediated insulation is dependent on DNA sequences flanking its core binding motifs, and CTCF binding sites at topologically associating domain(TAD) boundaries are more likely to function as insulators than those outside TAD boundaries, independent of binding strength. Using chromosomal conformation capture assays and high-resolution chromatin imaging techniques, we demonstrate that insulators form local chromatin domain boundaries and reduce enhancer-promoter contacts. Taken together, our results provide strong genetic, molecular, and structural evidence connecting chromatin topology to the action of insulators in the mammalian genome.

Chromosome spatial organization controls functional interactions between genes and regulators, yet the molecular and physical mechanisms underlying folding at the single DNA molecule level remain to be understood. Here we employ models of polymer physics to investigate the conformations of two 2Mb-wide DNA loci in human HCT116 and IMR90 wild-type and cohesin depleted cells. Model predictions on the 3D structure of single-molecules are consistently validated against super-resolution single-cell imaging data, providing evidence that the architecture of the studied loci is controlled by a thermodynamics mechanism of polymer phase separation whereby chromatin self-assembles in segregated globules. The process is driven by interactions between distinct types of cognate binding sites, correlating each with a different combination of chromatin factors, including CTCF, cohesin and histone marks. The intrinsic thermodynamics degeneracy of conformations results in a broad structural and time variability of single-molecules, reflected in their varying TAD-like contact patterns. Globules breathe in time, inducing stochastic unspecific interactions, yet they produce stable, compact environments where specific contacts become highly favored between regions enriched for cognate binding sites, albeit characterized by weak biochemical affinities. Cohesin depletion tends to reverse globule phase separation into a coil, randomly folded state, resulting in much more variable contacts across single-molecules, hence erasing population-averaged patterns. Overall, globule phase separation appears to be a robust, reversible mechanism of chromatin organization, where stochasticity and specificity coexist.

Neurons and oligodendrocytes are terminally differentiated cells that perform highly specialized functions, which depend on cascades of gene activation and repression to retain homeostatic control over a lifespan. Gene expression is regulated by three-dimensional (3D) genome organisation, from local levels of chromatin compaction to the organisation of topological domains and chromosome compartments. Whereas our understanding of 3D genome architecture has vastly increased in the past decade, it remains difficult to study specialized cells in their native environment without disturbing their activity. To develop the application of Genome Architecture Mapping (GAM) in small numbers of specialized cells in complex tissues, we combined GAM with immunoselection. We applied immunoGAM to map the genome architecture of specific cell populations in the juvenile/adult mouse brain: dopaminergic neurons (DNs) from the midbrain, pyramidal glutamatergic neurons (PGNs) from the hippocampus, and oligodendrocyte lineage cells (OLGs) from the cortex. We integrate 3D genome organisation with single-cell transcriptomics data, and find specific chromatin structures that relate with cell-type specific patterns of gene expression. We discover abundant changes in compartment organisation, especially a strengthening of heterochromatin compartments which establish strong contacts spanning tens of megabases, especially in brain cells. These compartments contain olfactory and taste receptor genes, which are de-repressed in a subpopulation of PGNs with molecular signatures of long-term potentiation (LTP). We also show extensive reorganisation of topological domains where activation of neuronal or oligodendrocyte genes coincides with formation of new TAD borders. Finally, we discover loss of TAD organisation, or ‘TAD melting’, at long (>1Mb) neuronal genes when they are most highly expressed. Our work shows that the 3D organisation of the genome is highly cell-type specific in terminally differentiated cells of the brain, and essential to better understand brain-specific mechanisms of gene regulation.

Summary: 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. : Chiariello et al. use polymer physics models to infer the 3D conformations of the murine α-globin locus. In the transition from ESCs to erythroid cells, the locus rearranges from an inactive highly intermingled conformation to an active, hairpin-shaped structure, marked by cell-specific three-way contacts accurately predicted by the models. Keywords: globin loci, higher-order chromatin organization, polymer physics, gene regulation

Summary Enhancer RNAs (eRNAs) are a subset of long noncoding RNA generated from genomic enhancers: they are thought to act as potent promoters of the expression of nearby genes through interaction with the transcriptional and epigenomic machineries. In the present work, we describe two eRNAs transcribed from the enhancer of Nkx2-5—a gene specifying a master cardiomyogenic lineage transcription factor (TF)—which we call Intergenic Regulatory Element Nkx2-5 Enhancers (IRENEs). The IRENEs are encoded, respectively, on the same strand (SS) and in the divergent direction (div) respect to the nearby gene. Of note, these two eRNAs have opposing roles in the regulation of Nkx2-5: IRENE-SS acts as a canonical promoter of transcription, whereas IRENE-div represses the activity of the enhancer through recruitment of the histone deacetylase sirtuin 1. Thus, we have identified an autoregulatory loop controlling expression of the master cardiac TF NKX2-5, in which one eRNA represses transcription., Graphical Abstract, Highlights • Two eRNAs (IRENE-SS, IRENE-div) with opposing functions are found upstream of Nkx2-5 • IRENE-SS works as a classical eRNA, acting as a transcriptional activator • IRENE-div acts unconventionally, functioning as a transcriptional repressor • IRENEs epigenetically control enhancer status and, subsequently, locus architecture, Biological Sciences; Molecular Biology; Molecular Mechanism of Gene Regulation

Innovative experimental protocols from Molecular Biology provided in recent years quantitative data about the structure of the cell nucleus. These technologies, such as Hi-C, GAM or SPRITE, revealed that the genome has a non-random three-dimensional (3D) spatial organization, which serves functional purposes. In order to dissect the complexity of chromosome folding, models from Polymer Physics have been employed, highlighting many key aspects of large-scale chromatin organization. A deep understanding of the molecular mechanisms underlying the genome architecture is currently a crucial problem in Biology, since chromatin misfolding or structural variants can reconfigure chromatin domains, thereby resulting in pathogenic phenotypes and disease. Here, we discuss a numerical Polymer-Physics-based approach (PRISMR), able to model 3D chromatin folding by using Machine Learning strategies informed with experimental data. Using as a case study the Pitx1 locus, a genomic region critically involved in hindlimb development, we show that the PRISMR algorithm reproduces in silico with high accuracy the experimental contact data, thus providing a powerful computational tool for analyzing and predicting the 3D chromatin structure.

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.

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.

The constant development of sophisticated technologies is allowing to dissect three-dimensional chromatin structure at high resolution level. The tremendous amount of quantitative experimental data available today requires a conceptual framework able to make sense of them. In this perspective, polymer physics offers a key tool to interpret chromatin architecture data and to unveil the basic mechanisms shaping its structure. In the very last years, several polymer models have been proposed and have allowed to capture complex features emerging from the data. The major peculiarity distinguishing the different models is represented by the more or less complicated physical mechanism used to explain chromatin folding. Here, we review very popular models which have been recently developed and which represent brilliant examples from this interdisciplinary research field. In order to highlight the wide range of practical applications they have, we discuss the cases of the murine Pitx1 and the human EPHA4 loci, showing that polymer physics allows to effectively study chromatin structure in different cell lines and to predict the impact of pathogenic structural variants on the genome three-dimensional architecture. This article is categorized under: Structure and Mechanism > Computational Biochemistry and Biophysics Molecular and Statistical Mechanics > Molecular Dynamics and Monte-Carlo Methods Theoretical and Physical Chemistry > Statistical Mechanics.

In recent years the development of novel technologies, as Hi-C or GAM, allowed to investigate the spatial structure of chromatin in the cell nucleus with a constantly increasing level of accuracy. Polymer physics models have been developed and improved to better interpret the wealth of complex information coming from the experimental data, providing highly accurate understandings on chromatin architecture and on the mechanisms regulating genome folding. To investigate the capability of the models to explain the experiments and to test their agreement with the data, massive parallel simulations are needed and efficient algorithms are fundamental. In this work, we consider general computational Molecular Dynamics (MD) techniques commonly used to implement such models, with a special focus on the Strings & Binders Switch polymer model. By combining this model with machine learning computational approaches, it is possible to give an accurate description of real genomic loci. In addition, it is also possible to make predictions about the impact of structural variants of the genomic sequence, which are known to be linked to severe congenital diseases.

It is not clear how spontaneous DNA double-strand breaks (DSBs) form and are processed in normal cells, and whether they predispose to cancer-associated translocations. We show that DSBs in normal mammary cells form upon release of paused RNA polymerase II (Pol II) at promoters, 5' splice sites and active enhancers, and are processed by end-joining in the absence of a canonical DNA-damage response. Logistic and causal-association models showed that Pol II pausing at long genes is the main predictor and determinant of DSBs. Damaged introns with paused Pol II-pS5, TOP2B and XRCC4 are enriched in translocation breakpoints, and map at topologically associating domain boundary-flanking regions showing high interaction frequencies with distal loci. Thus, in unperturbed growth conditions, release of paused Pol II at specific loci and chromatin territories favors DSB formation, leading to chromosomal translocations.

The depth and complexity of data now available on chromosome 3D architecture, derived by new technologies such as Hi-C, have triggered the development of models based on polymer physics to explain the observed patterns and the underlying molecular folding mechanisms. Here, we give an overview of some of the ideas and models from physics introduced to date, along with their progresses and limitations in the description of experimental data. In particular, we focus on the Strings&Binders and the Loop Extrusion model of chromatin architecture. This article is categorized under: Analytical and Computational Methods > Computational Methods.

Long-range gene regulation involves physical proximity between enhancers and promoters to generate precise patterns of gene expression in space and time. However, in some cases, proximity coincides with gene activation, whereas, in others, preformed topologies already exist before activation. In this study, we investigate the preformed configuration underlying the regulation of the Shh gene by its unique limb enhancer, the ZRS , in vivo during mouse development. Abrogating the constitutive transcription covering the ZRS region led to a shift within the Shh–ZRS contacts and a moderate reduction in Shh transcription. Deletion of the CTCF binding sites around the ZRS resulted in the loss of the Shh–ZRS preformed interaction and a 50% decrease in Shh expression but no phenotype, suggesting an additional, CTCF-independent mechanism of promoter–enhancer communication. This residual activity, however, was diminished by combining the loss of CTCF binding with a hypomorphic ZRS allele, resulting in severe Shh loss of function and digit agenesis. Our results indicate that the preformed chromatin structure of the Shh locus is sustained by multiple components and acts to reinforce enhancer–promoter communication for robust transcription.

Complex architectural rearrangements are associated to the control of the HoxD genes in different cell types; yet, how they are implemented in single cells remains unknown. By use of polymer models, we dissect the locus 3D structure at the single DNA molecule level in mouse embryonic stem and cortical neuronal cells, as the HoxD cluster changes from a poised to a silent state. Our model describes published Hi-C, 3-way 4C, and FISH data with high accuracy and is validated against independent 4C data on the Nsi-SB 0.5-Mb duplication and on triple contacts. It reveals the mode of action of compartmentalization on the regulation of the HoxD genes that have gene- and cell-type-specific multi-way interactions with their regulatory elements and high cell-to-cell variability. It shows that TADs and higher-order 3D structures, such as metaTADs, associate with distinct combinations of epigenetic factors, including but not limited to CCCTC-binding factor (CTCF) and histone marks. Bianco et al. reconstruct the 3D structure of the murine HoxD locus by using polymer models at the single DNA molecule level. The locus architecture rearranges upon differentiation from embryonic stem cells to cortical neurons in connection to epigenetic changes, as cell-type- and gene-specific multi-way contacts are established with regulatory elements.

The promoters of mammalian genes are commonly regulated by multiple distal enhancers, which physically interact within discrete chromatin domains. How such domains form and how the regulatory elements within them interact within single cells is not understood. To address this we developed Tri-C, a new Chromosome Conformation Capture (3C) approach to identify concurrent chromatin interactions at individual alleles within single cells. The heterogeneity of interactions observed between such cells shows that CTCF-mediated formation of chromatin domains and interactions within them are dynamic processes. Importantly, our analyses reveal higher-order structures involving simultaneous interactions between multiple enhancers and promoters within individual cells. This provides a structural basis for understanding how multiple cis-elements act together to establish robust regulation of gene expression.

The spatial organisation of interphase chromosomes is known to affect genomic function, yet the principles behind such organisation remain elusive. Here, we first compare and then combine two well-known biophysical models, the transcription factor (TF) and loop extrusion (LE) models, and dissect their respective roles in organising the genome. Our results suggest that extrusion and transcription factors play complementary roles in folding the genome: the former are necessary to compact gene deserts or “inert chromatin” regions, the latter are sufficient to explain most of the structure found in transcriptionally active or repressed domains. Finally, we find that to reproduce interaction patterns found in HiC experiments we do not need to postulate an explicit motor activity of cohesin (or other extruding factors): a model where co-hesin molecules behave as molecular slip-links sliding diffusively along chromatin works equally well.

The regulatory specificity of enhancers and their interaction with gene promoters is thought to be controlled by their sequence and the binding of transcription factors. By studying Pitx1, a regulator of hindlimb development, we show that dynamic changes in chromatin conformation can restrict the activity of enhancers. Inconsistent with its hindlimb-restricted expression, Pitx1 is controlled by an enhancer (Pen) that shows activity in forelimbs and hindlimbs. By Capture Hi-C and three-dimensional modeling of the locus, we demonstrate that forelimbs and hindlimbs have fundamentally different chromatin configurations, whereby Pen and Pitx1 interact in hindlimbs and are physically separated in forelimbs. Structural variants can convert the inactive into the active conformation, thereby inducing Pitx1 misexpression in forelimbs, causing partial arm-to-leg transformation in mice and humans. Thus, tissue-specific three-dimensional chromatin conformation can contribute to enhancer activity and specificity in vivo and its disturbance can result in gene misexpression and disease.