Complex small-world regulatory networks emerge from the 3D organisation of the human genome

Biomarker Selection for Adaptive Systems

Biomarkers enable objective monitoring of a given cell or state in a biological system and are widely used in research, biomanufacturing, and clinical practice. However, identifying appropriate biomarkers that are both robustly measurable and capture a state accurately remains challenging. We present a framework for biomarker identification based upon observability guided sensor selection. Our methods, Dynamic Sensor Selection (DSS) and Structure-Guided Sensor Selection (SGSS), utilize temporal models and experimental data, offering a template for applying observability theory to data from biological systems. Unlike conventional methods that assume well-known, fixed dynamics, DSS adaptively select biomarkers or sensors that maximize observability while accounting for the time-varying nature of biological systems. Additionally, SGSS incorporates structural information and diverse data to identify sensors which are resilient against inaccuracies in our model of the underlying system. We validate our approaches by performing estimation on high dimensional systems derived from temporal gene expression data from partial observations. Our algorithms reliably identify known biomarkers and uncover new ones within our datasets. Additionally, integrating chromosome conformation and gene expression data addresses noise and uncertainty, enhancing the reliability of our biomarker selection approach for the genome.

Topological Spectra and Entropy of Chromatin Loop Networks

The 3D folding of a mammalian gene can be studied by a polymer model, where the chromatin fiber is represented by a semiflexible polymer which interacts with multivalent proteins, representing complexes of DNA-binding transcription factors and RNA polymerases. This physical model leads to the natural emergence of clusters of proteins and binding sites, accompanied by the folding of chromatin into a set of topologies, each associated with a different network of loops. Here, we combine numerics and analytics to first classify these networks and then find their relative importance or statistical weight, when the properties of the underlying polymer are those relevant to chromatin. Unlike polymer networks previously studied, our chromatin networks have finite average distances between successive binding sites, and this leads to giant differences between the weights of topologies with the same number of edges and nodes but different wiring. These weights strongly favor rosettelike structures with a local cloud of loops with respect to more complicated nonlocal topologies. Our results suggest that genes should overwhelmingly fold into a small fraction of all possible 3D topologies, which can be robustly characterized by the framework we propose here.

The scales, mechanisms, and dynamics of the genome architecture

Even when split into several chromosomes, DNA molecules that make up our genome are too long to fit into the cell nuclei unless massively folded. Such folding must accommodate the need for timely access to selected parts of the genome by transcription factors, RNA polymerases, and DNA replication machinery. Here, we review our current understanding of the genome folding inside the interphase nuclei. We consider the resulting genome architecture at three scales with a particular focus on the intermediate (meso) scale and summarize the insights gained from recent experimental observations and diverse computational models.

Chromatin phase separated nanoregions explored by polymer cross-linker models and reconstructed from single particle trajectories

Phase separated domains (PSDs) are ubiquitous in cell biology, representing nanoregions of high molecular concentration. PSDs appear at diverse cellular domains, such as neuronal synapses but also in eukaryotic cell nucleus, limiting the access of transcription factors and thus preventing gene expression. We develop a generalized cross-linker polymer model, to study PSDs: we show that increasing the number of cross-linkers induces a polymer condensation, preventing access of diffusing molecules. To investigate how the PSDs restrict the motion of diffusing molecules, we compute the mean residence and first escaping times. Finally, we develop a method based on mean-square-displacement of single particle trajectories to reconstruct the properties of PSDs from the continuum range of anomalous exponents. We also show here that PSD generated by polymers do not induces a long-range attracting field (potential well), in contrast with nanodomains at neuronal synapses. To conclude, PSDs can result from condensed chromatin organization, where the number of cross-linkers controls molecular access.

Transcription modulates chromatin dynamics and locus configuration sampling

In living cells, the 3D structure of gene loci is dynamic, but this is not revealed by 3C and FISH experiments in fixed samples, leaving a notable gap in our understanding. To overcome these limitations, we applied the highly predictive heteromorphic polymer (HiP-HoP) model to determine chromatin fiber mobility at the Pax6 locus in three mouse cell lines with different transcription states. While transcriptional activity minimally affects movement of 40-kbp regions, we observed that motion of smaller 1-kbp regions depends strongly on local disruption to chromatin fiber structure marked by H3K27 acetylation. This also substantially influenced locus configuration dynamics by modulating protein-mediated promoter-enhancer loops. Importantly, these simulations indicate that chromatin dynamics are sufficiently fast to sample all possible locus conformations within minutes, generating wide dynamic variability within single cells. This combination of simulation and experimental validation provides insight into how transcriptional activity influences chromatin structure and gene dynamics.

Ultra-long-range interactions between active regulatory elements

Contacts between enhancers and promoters are thought to relate to their ability to activate transcription. Investigating factors that contribute to such chromatin interactions is therefore important for understanding gene regulation. Here, we have determined contact frequencies between millions of pairs of cis-regulatory elements from chromosome conformation capture data sets and analyzed a collection of hundreds of DNA-binding factors for binding at regions of enriched contacts. This analysis revealed enriched contacts at sites bound by many factors associated with active transcription. We show that active regulatory elements, independent of cohesin and polycomb, interact with each other across distances of tens of megabases in vertebrate and invertebrate genomes and that interactions correlate and change with activity. However, these ultra-long-range interactions are not dependent on RNA polymerase II transcription or individual transcription cofactors. Using simulations, we show that a model of chromatin and multivalent binding factors can give rise to long-range interactions via bridging-induced clustering. We propose that long-range interactions between cis-regulatory elements are driven by at least three distinct processes: cohesin-mediated loop extrusion, polycomb contacts, and clustering of active regions.

Esearch3D: propagating gene expression in chromatin networks to illuminate active enhancers

Most cell type-specific genes are regulated by the interaction of enhancers with their promoters. The identification of enhancers is not trivial as enhancers are diverse in their characteristics and dynamic in their interaction partners. We present Esearch3D, a new method that exploits network theory approaches to identify active enhancers. Our work is based on the fact that enhancers act as a source of regulatory information to increase the rate of transcription of their target genes and that the flow of this information is mediated by the folding of chromatin in the three-dimensional (3D) nuclear space between the enhancer and the target gene promoter. Esearch3D reverse engineers this flow of information to calculate the likelihood of enhancer activity in intergenic regions by propagating the transcription levels of genes across 3D genome networks. Regions predicted to have high enhancer activity are shown to be enriched in annotations indicative of enhancer activity. These include: enhancer-associated histone marks, bidirectional CAGE-seq, STARR-seq, P300, RNA polymerase II and expression quantitative trait loci (eQTLs). Esearch3D leverages the relationship between chromatin architecture and transcription, allowing the prediction of active enhancers and an understanding of the complex underpinnings of regulatory networks. The method is available at: https://github.com/InfOmics/Esearch3D and https://doi.org/10.5281/zenodo.7737123.

Network models of chromatin structure

Increasing numbers of datasets and experimental assays that capture the organization of chromatin inside the nucleus warrant an effort to develop tools to visualize and analyze these structures. Alongside polymer physics or constraint-based modeling, network theory approaches to describe 3D epigenome organization have gained in popularity. Representing genomic regions as nodes in a network enables visualization of 1D epigenomics datasets in the context of chromatin structure maps, while network theory metrics can be used to describe 3D epigenome organization and dynamics. In this review, we summarize the most salient applications of network theory to the study of chromatin contact maps, demonstrating its potential in revealing epigenomic patterns and relating them to cellular phenotypes.

Synthetic analysis of chromatin tracing and live-cell imaging indicates pervasive spatial coupling between genes

The role of the spatial organization of chromosomes in directing transcription remains an outstanding question in gene regulation. Here, we analyze two recent single-cell imaging methodologies applied across hundreds of genes to systematically analyze the contribution of chromosome conformation to transcriptional regulation. Those methodologies are (1) single-cell chromatin tracing with super-resolution imaging in fixed cells; and (2) high-throughput labeling and imaging of nascent RNA in living cells. Specifically, we determine the contribution of physical distance to the coordination of transcriptional bursts. We find that individual genes adopt a constrained conformation and reposition toward the centroid of the surrounding chromatin upon activation. Leveraging the variability in distance inherent in single-cell imaging, we show that physical distance - but not genomic distance - between genes on individual chromosomes is the major factor driving co-bursting. By combining this analysis with live-cell imaging, we arrive at a corrected transcriptional correlation of [Formula: see text] for genes separated by < 400 nm. We propose that this surprisingly large correlation represents a physical property of human chromosomes and establishes a benchmark for future experimental studies.

Deciphering multi-way interactions in the human genome

Chromatin architecture, a key regulator of gene expression, can be inferred using chromatin contact data from chromosome conformation capture, or Hi-C. However, classical Hi-C does not preserve multi-way contacts. Here we use long sequencing reads to map genome-wide multi-way contacts and investigate higher order chromatin organization in the human genome. We use hypergraph theory for data representation and analysis, and quantify higher order structures in neonatal fibroblasts, biopsied adult fibroblasts, and B lymphocytes. By integrating multi-way contacts with chromatin accessibility, gene expression, and transcription factor binding, we introduce a data-driven method to identify cell type-specific transcription clusters. We provide transcription factor-mediated functional building blocks for cell identity that serve as a global signature for cell types.

Mapping higher order chromatin architecture is important. Here the authors use long sequencing reads to map genome-wide multi-way contacts and investigate higher order chromatin organisation; they use hypergraph theory for data representation and analysis, and apply this to different cell types.

High-resolution simulations of chromatin folding at genomic rearrangements in malignant B cells provide mechanistic insights into proto-oncogene deregulation

Genomic rearrangements are known to result in proto-oncogene deregulation in many cancers, but the link to 3D genome structure remains poorly understood. Here, we used the highly predictive heteromorphic polymer (HiP-HoP) model to predict chromatin conformations at the proto-oncogene CCND1 in healthy and malignant B cells. After confirming that the model gives good predictions of Hi-C data for the nonmalignant human B cell-derived cell line GM12878, we generated predictions for two cancer cell lines, U266 and Z-138. These possess genome rearrangements involving CCND1 and the immunoglobulin heavy locus (IGH), which we mapped using targeted genome sequencing. Our simulations showed that a rearrangement in U266 cells where a single IGH super-enhancer is inserted next to CCND1 leaves the local topologically associated domain (TAD) structure intact. We also observed extensive changes in enhancer-promoter interactions within the TAD, suggesting that it is the downstream chromatin remodeling which gives rise to the oncogene activation, rather than the presence of the inserted super-enhancer DNA sequence per se. Simulations of the IGH-CCND1 reciprocal translocation in Z-138 cells revealed that an oncogenic fusion TAD is created, encompassing CCND1 and the IGH super-enhancers. We predicted how the structure and expression of CCND1 changes in these different cell lines, validating this using qPCR and fluorescence in situ hybridization microscopy. Our work demonstrates the power of polymer simulations to predict differences in chromatin interactions and gene expression for different translocation breakpoints.

The 3D architecture of the pepper genome and its relationship to function and evolution

The organization of chromatin into self-interacting domains is universal among eukaryotic genomes, though how and why they form varies considerably. Here we report a chromosome-scale reference genome assembly of pepper (Capsicum annuum) and explore its 3D organization through integrating high-resolution Hi-C maps with epigenomic, transcriptomic, and genetic variation data. Chromatin folding domains in pepper are as prominent as TADs in mammals but exhibit unique characteristics. They tend to coincide with heterochromatic regions enriched with retrotransposons and are frequently embedded in loops, which may correlate with transcription factories. Their boundaries are hotspots for chromosome rearrangements but are otherwise depleted for genetic variation. While chromatin conformation broadly affects transcription variance, it does not predict differential gene expression between tissues. Our results suggest that pepper genome organization is explained by a model of heterochromatin-driven folding promoted by transcription factories and that such spatial architecture is under structural and functional constraints.

Predicting genome organisation and function with mechanistic modelling

Fitting-free mechanistic models based on polymer simulations predict chromatin folding in 3D by focussing on the underlying biophysical mechanisms. This class of models has been increasingly used in conjunction with experiments to study the spatial organisation of eukaryotic chromosomes. Feedback from experiments to models leads to successive model refinement and has previously led to the discovery of new principles for genome organisation. Here, we review the basis of mechanistic polymer simulations, explain some of the more recent approaches and the contexts in which they have been useful to explain chromosome biology, and speculate on how they might be used in the future.