CHROMATIX: computing the functional landscape of many-body chromatin interactions in transcriptionally active loci from deconvolved single cells

SEgene identifies links between super enhancers and gene expression across cell types

Enhancers are non-coding DNA regions that facilitate gene transcription, with a specialized subset, super-enhancers, known to exert exceptionally strong transcriptional activation effects. Super-enhancers have been implicated in oncogenesis, and their identification is achievable through histone mark chromatin immunoprecipitation followed by sequencing data using existing analytical tools. However, conventional super-enhancer detection methodologies often do not accurately reflect actual gene expression levels, and the large volume of identified super-enhancers complicates comprehensive analysis. To address these limitations, we developed the super-enhancer to gene links (SE-to-gene Links) analysis, a platform named "SEgene" which incorporates the peak-to-gene links approach-a statistical method designed to reveal correlations between genes and peak regions ( https://github.com/hamamoto-lab/SEgene ). This platform enables a targeted evaluation of super-enhancer regions in relation to gene expression, facilitating the identification of super-enhancers that are functionally linked to transcriptional activity. Here, we demonstrate the application of SE-to-gene Links analysis to public datasets, confirming its efficacy in accurately detecting super-enhancers and identifying functionally associated genes. Additionally, SE-to-gene Links analysis identified ERBB2 as a significant gene of interest in the lung adenocarcinoma dataset from the National Cancer Center Japan cohort, suggesting a potential impact across multiple patient samples. Thus, the SE-to-gene Links analysis provides an analytical tool for evaluating super-enhancers as potential therapeutic targets, supporting the identification of clinically significant super-enhancer regions and their functionally associated genes.

Organization and Dynamics of Chromosomes

How long thread-like eukaryotic chromosomes fit tidily in the small volume of the nucleus without significant entanglement is just beginning to be understood, thanks to major advances in experimental techniques. Several polymer models, which reproduce contact maps that measure the probabilities that two loci are in spatial contact, have predicted the 3D structures of interphase chromosomes. Data-driven approaches, using contact maps as input, predict that mitotic helical chromosomes are characterized by a switch in handedness, referred to as perversion. By using experimentally derived effective interactions between chromatin loci in simulations, structures of conventional and inverted nuclei have been accurately predicted. Polymer theory and simulations show that the dynamics of individual loci in chromatin exhibit subdiffusive behavior but the diffusion exponents are broadly distributed, which accords well with experiments. Although coarse-grained models are successful, many challenging problems remain, which require the creation of new experimental and computational tools to understand genome biology.

Effects of Lamina-Chromatin Attachment on Super Long-Range Chromatin Interactions

The interactions between chromatin and lamin proteins localized on the nuclear envelope play a crucial role in the three-dimensional (3D) organization of the genome. This study investigates the influence of lamin associated domains (LADs) on genome organization at the chromosome level using 3D polymer models of mouse embryonic fibroblasts (MEFs) and embryonic stem cells (mESCs). By integrating genome-wide LAD maps from DamID assays, we simulated chromatin conformations with and without LAD attachment to the nuclear envelope. Our results show that incorporating LAD-lamin interactions yields a radial chromatin distribution consistent with experimental observations. Moreover, LAD-lamin interactions induce significant super long-range chromatin contacts across distant genomic regions. These findings suggest two distinct mechanisms driving induction of chromatin interactions by LAD-lamin attachment.

Structural basis of differential gene expression at eQTLs loci from high-resolution ensemble models of 3D single-cell chromatin conformations

MOTIVATION

Techniques such as high-throughput chromosome conformation capture (Hi-C) have provided a wealth of information on nucleus organization and genome important for understanding gene expression regulation. Genome-Wide Association Studies have identified numerous loci associated with complex traits. Expression quantitative trait loci (eQTL) studies have further linked the genetic variants to alteration in expression levels of associated target genes across individuals. However, the functional roles of many eQTLs in noncoding regions remain unclear. Current joint analyses of Hi-C and eQTLs data lack advanced computational tools, limiting what can be learned from these data.

RESULTS

We developed a computational method for simultaneous analysis of Hi-C and eQTL data, capable of identifying a small set of nonrandom interactions from all Hi-C interactions. Using these nonrandom interactions, we reconstructed large ensembles (×105) of high-resolution single-cell 3D chromatin conformations with thorough sampling, accurately replicating Hi-C measurements. Our results revealed many-body interactions in chromatin conformation at the single-cell level within eQTL loci, providing a detailed view of how 3D chromatin structures form the physical foundation for gene regulation, including how genetic variants of eQTLs affect the expression of associated eGenes. Furthermore, our method can deconvolve chromatin heterogeneity and investigate the spatial associations of eQTLs and eGenes at subpopulation level, revealing their regulatory impacts on gene expression. Together, ensemble modeling of thoroughly sampled single-cell chromatin conformations combined with eQTL data, helps decipher how 3D chromatin structures provide the physical basis for gene regulation, expression control, and aid in understanding the overall structure-function relationships of genome organization.

AVAILABILITY AND IMPLEMENTATION

It is available at https://github.com/uic-liang-lab/3DChromFolding-eQTL-Loci.

Advancements and future directions in single-cell Hi-C based 3D chromatin modeling

Single-cell Hi-C data provides valuable insights into the three-dimensional organization of chromatin within individual cells, yet modeling this data poses significant challenges due to its inherent sparsity and variability. This review comprehensively explores the predominant approaches to reconstructing 3D chromatin structures from single-cell Hi-C data, positioning these methods within the broader contexts of single-cell Hi-C research and bulk Hi-C data modeling. We categorize the modeling strategies based on their objective functions, which are framed in terms of force fields, potentials, cost functions, or likelihood probabilities. Despite their diverse methodologies, these approaches exhibit deep underlying similarities. We further dissect the basic components of these models, such as attractive restraint forces and repulsive forces, and discuss additional terms like fluid viscosity and variation penalties. The review also critically evaluates the current state of model validation, highlighting the inconsistencies across various studies and emphasizing the need for a comprehensive validation framework. We detail common validation techniques, including the comparison of distance matrices and the assessment of contact violations. We argue that the future of single-cell Hi-C modeling lies in integrating multiple data modalities and incorporating cell cycle trajectory information. Such integration could significantly advance our understanding of chromatin conformation dynamics during cell cycle progression and cell differentiation. We also foresee the continued growth of optimization-based and molecular dynamics approaches, supported by general molecular dynamics toolkits.

Computational methods for analysing multiscale 3D genome organization

Recent progress in whole-genome mapping and imaging technologies has enabled the characterization of the spatial organization and folding of the genome in the nucleus. In parallel, advanced computational methods have been developed to leverage these mapping data to reveal multiscale three-dimensional (3D) genome features and to provide a more complete view of genome structure and its connections to genome functions such as transcription. Here, we discuss how recently developed computational tools, including machine-learning-based methods and integrative structure-modelling frameworks, have led to a systematic, multiscale delineation of the connections among different scales of 3D genome organization, genomic and epigenomic features, functional nuclear components and genome function. However, approaches that more comprehensively integrate a wide variety of genomic and imaging datasets are still needed to uncover the functional role of 3D genome structure in defining cellular phenotypes in health and disease.

Efficient Hi-C inversion facilitates chromatin folding mechanism discovery and structure prediction

Genome-wide chromosome conformation capture (Hi-C) experiments have revealed many structural features of chromatin across multiple length scales. Further understanding genome organization requires relating these discoveries to the mechanisms that establish chromatin structures and reconstructing these structures in three dimensions, but both objectives are difficult to achieve with existing algorithms that are often computationally expensive. To alleviate this challenge, we present an algorithm that efficiently converts Hi-C data into contact energies, which measure the interaction strength between genomic loci brought into proximity. Contact energies are local quantities unaffected by the topological constraints that correlate Hi-C contact probabilities. Thus, extracting contact energies from Hi-C contact probabilities distills the biologically unique information contained in the data. We show that contact energies reveal the location of chromatin loop anchors, support a phase separation mechanism for genome compartmentalization, and parameterize polymer simulations that predict three-dimensional chromatin structures. Therefore, we anticipate that contact energy extraction will unleash the full potential of Hi-C data and that our inversion algorithm will facilitate the widespread adoption of contact energy analysis.

Techniques for and challenges in reconstructing 3D genome structures from 2D chromosome conformation capture data

Chromosome conformation capture technologies that provide frequency information for contacts between genomic regions have been crucial for increasing our understanding of genome folding and regulation. However, such data do not provide direct evidence of the spatial 3D organization of chromatin. In this opinion article, we discuss the development and application of computational methods to reconstruct chromatin 3D structures from experimental 2D contact data, highlighting how such modeling provides biological insights and can suggest mechanisms anchored to experimental data. By applying different reconstruction methods to the same contact data, we illustrate some state-of-the-art of these techniques and discuss our gene resolution approach based on Brownian dynamics and Monte Carlo sampling.

Efficient Hi-C inversion facilitates chromatin folding mechanism discovery and structure prediction

Genome-wide chromosome conformation capture (Hi-C) experiments have revealed many structural features of chromatin across multiple length scales. Further understanding genome organization requires relating these discoveries to the mechanisms that establish chromatin structures and reconstructing these structures in three dimensions, but both objectives are difficult to achieve with existing algorithms that are often computationally expensive. To alleviate this challenge, we present an algorithm that efficiently converts Hi-C data into contact energies, which measure the interaction strength between genomic loci brought into proximity. Contact energies are local quantities unaffected by the topological constraints that correlate Hi-C contact probabilities. Thus, extracting contact energies from Hi-C contact probabilities distills the biologically unique information contained in the data. We show that contact energies reveal the location of chromatin loop anchors, support a phase separation mechanism for genome compartmentalization, and parameterize polymer simulations that predict three-dimensional chromatin structures. Therefore, we anticipate that contact energy extraction will unleash the full potential of Hi-C data and that our inversion algorithm will facilitate the widespread adoption of contact energy analysis.

Predicting Pathology of Missense Mutations through Protein-Specific Evolutionary Pattern

Missense mutations, which are single base pair genetic alternation resulting in a different amino acid, are among the most common occurring variants in exon regions of the human genome and may lead to diseases. Thus to assess the effects of missense mutations, it is essential to investigate the evolutionary history of the protein under selection pressures. In this study, we employ a continuous-time Markov model to investigate the evolutionary patterns in protein sequences and a Bayesian Markov chain Monte Carlo method to estimate the substitution rates for protein of interest, from which we obtain scoring matrices. Specifically, we examined the evolutionary patterns of protein sequences containing missense mutations using a species tree to define the phylogeny of the protein of interest. We thoroughly studied the evolutionary pattern of human muscle glycogen phosphorylase containing 127 known missense mutations, and identified characteristic evolutionary patterns in 63 proteins with 2,238 missense mutations, including both deleterious and neutral effects. Our results show that the estimated protein-specific evolutionary pattern-based scoring matrices (PSM) lead to higher sensitivity in detecting the pathological effects of missense mutations, compared to the general evolutionary pattern-based scoring matrix of Blosum62 (BL62) matrix. By incorporating PSM, the performance of a recently released structure-based model SPRI for evaluating missense mutations is further improved.

Analysis of super-enhancer using machine learning and its application to medical biology

The analysis of super-enhancers (SEs) has recently attracted attention in elucidating the molecular mechanisms of cancer and other diseases. SEs are genomic structures that strongly induce gene expression and have been reported to contribute to the overexpression of oncogenes. Because the analysis of SEs and integrated analysis with other data are performed using large amounts of genome-wide data, artificial intelligence technology, with machine learning at its core, has recently begun to be utilized. In promoting precision medicine, it is important to consider information from SEs in addition to genomic data; therefore, machine learning technology is expected to be introduced appropriately in terms of building a robust analysis platform with a high generalization performance. In this review, we explain the history and principles of SE, and the results of SE analysis using state-of-the-art machine learning and integrated analysis with other data are presented to provide a comprehensive understanding of the current status of SE analysis in the field of medical biology. Additionally, we compared the accuracy between existing machine learning methods on the benchmark dataset and attempted to explore the kind of data preprocessing and integration work needed to make the existing algorithms work on the benchmark dataset. Furthermore, we discuss the issues and future directions of current SE analysis.

The magic of unraveling genome architecture and function

Over the last decades, technological breakthroughs in super-resolution microscopy have allowed us to reach molecular resolution and design experiments of unprecedented complexity. Investigating how chromatin is folded in 3D, from the nucleosome level up to the entire genome, is becoming possible by "magic" (imaging genomic), i.e., the combination of imaging and genomic approaches. This offers endless opportunities to delve into the relationship between genome structure and function. Here, we review recently achieved objectives and the conceptual and technical challenges the field of genome architecture is currently undertaking. We discuss what we have learned so far and where we are heading. We elucidate how the different super-resolution microscopy approaches and, more specifically, live-cell imaging have contributed to the understanding of genome folding. Moreover, we discuss how future technical developments could address remaining open questions.

An Igh distal enhancer modulates antigen receptor diversity by determining locus conformation

The mouse Igh locus is organized into a developmentally regulated topologically associated domain (TAD) that is divided into subTADs. Here we identify a series of distal VH enhancers (EVHs) that collaborate to configure the locus. EVHs engage in a network of long-range interactions that interconnect the subTADs and the recombination center at the DHJH gene cluster. Deletion of EVH1 reduces V gene rearrangement in its vicinity and alters discrete chromatin loops and higher order locus conformation. Reduction in the rearrangement of the VH11 gene used in anti-PtC responses is a likely cause of the observed reduced splenic B1 B cell compartment. EVH1 appears to block long-range loop extrusion that in turn contributes to locus contraction and determines the proximity of distant VH genes to the recombination center. EVH1 is a critical architectural and regulatory element that coordinates chromatin conformational states that favor V(D)J rearrangement.

Modeling and simulation of cell nuclear architecture reorganization process

We develop a special phase field/diffusive interface method to model the nuclear architecture reorganization process. In particular, we use a Lagrange multiplier approach in the phase field model to preserve the specific physical and geometrical constraints for the biological events. We develop several efficient and robust linear and weakly nonlinear schemes for this new model. To validate the model and numerical methods, we present ample numerical simulations which in particular reproduce several processes of nuclear architecture reorganization from the experiment literature.

Extracting multi-way chromatin contacts from Hi-C data

There is a growing realization that multi-way chromatin contacts formed in chromosome structures are fundamental units of gene regulation. However, due to the paucity and complexity of such contacts, it is challenging to detect and identify them using experiments. Based on an assumption that chromosome structures can be mapped onto a network of Gaussian polymer, here we derive analytic expressions for n-body contact probabilities (n > 2) among chromatin loci based on pairwise genomic contact frequencies available in Hi-C, and show that multi-way contact probability maps can in principle be extracted from Hi-C. The three-body (triplet) contact probabilities, calculated from our theory, are in good correlation with those from measurements including Tri-C, MC-4C and SPRITE. Maps of multi-way chromatin contacts calculated from our analytic expressions can not only complement experimental measurements, but also can offer better understanding of the related issues, such as cell-line dependent assemblies of multiple genes and enhancers to chromatin hubs, competition between long-range and short-range multi-way contacts, and condensates of multiple CTCF anchors.

Perspectives for the reconstruction of 3D chromatin conformation using single cell Hi-C data

Construction of chromosomes 3D models based on single cell Hi-C data constitute an important challenge. We present a reconstruction approach, DPDchrom, that incorporates basic knowledge whether the reconstructed conformation should be coil-like or globular and spring relaxation at contact sites. In contrast to previously published protocols, DPDchrom can naturally form globular conformation due to the presence of explicit solvent. Benchmarking of this and several other methods on artificial polymer models reveals similar reconstruction accuracy at high contact density and DPDchrom advantage at low contact density. To compare 3D structures insensitively to spatial orientation and scale, we propose the Modified Jaccard Index. We analyzed two sources of the contact dropout: contact radius change and random contact sampling. We found that the reconstruction accuracy exponentially depends on the number of contacts per genomic bin allowing to estimate the reconstruction accuracy in advance. We applied DPDchrom to model chromosome configurations based on single-cell Hi-C data of mouse oocytes and found that these configurations differ significantly from a random one, that is consistent with other studies.

Single-cell Hi-C data analysis: safety in numbers

Over the past decade, genome-wide assays for chromatin interactions in single cells have enabled the study of individual nuclei at unprecedented resolution and throughput. Current chromosome conformation capture techniques survey contacts for up to tens of thousands of individual cells, improving our understanding of genome function in 3D. However, these methods recover a small fraction of all contacts in single cells, requiring specialised processing of sparse interactome data. In this review, we highlight recent advances in methods for the interpretation of single-cell genomic contacts. After discussing the strengths and limitations of these methods, we outline frontiers for future development in this rapidly moving field.

Minimalistic 3D chromatin models: Sparse interactions in single cells drive the chromatin fold and form many-body units

Computational three-dimensional chromatin modeling has helped uncover principles of genome organization. Here, we discuss methods for modeling three-dimensional chromatin structures, with focus on a minimalistic polymer model which inverts population Hi-C into single-cell conformations. Utilizing only basic physical properties, this model reveals that a few specific Hi-C interactions can fold chromatin into conformations consistent with single-cell imaging, Dip-C, and FISH measurements. Aggregated single-cell chromatin conformations also reproduce Hi-C frequencies. This approach allows quantification of structural heterogeneity and discovery of many-body interaction units and has revealed additional insights, including (1) topologically associating domains as a byproduct of folding driven by specific interactions, (2) cell subpopulations with different structural scaffolds are developmental stage dependent, and (3) the functional landscape of many-body units within enhancer-rich regions. We also discuss these findings in relation to the genome structure-function relationship.

Multiscale modeling of genome organization with maximum entropy optimization

Three-dimensional (3D) organization of the human genome plays an essential role in all DNA-templated processes, including gene transcription, gene regulation, and DNA replication. Computational modeling can be an effective way of building high-resolution genome structures and improving our understanding of these molecular processes. However, it faces significant challenges as the human genome consists of over 6 × 109 base pairs, a system size that exceeds the capacity of traditional modeling approaches. In this perspective, we review the progress that has been made in modeling the human genome. Coarse-grained models parameterized to reproduce experimental data via the maximum entropy optimization algorithm serve as effective means to study genome organization at various length scales. They have provided insight into the principles of whole-genome organization and enabled de novo predictions of chromosome structures from epigenetic modifications. Applications of these models at a near-atomistic resolution further revealed physicochemical interactions that drive the phase separation of disordered proteins and dictate chromatin stability in situ. We conclude with an outlook on the opportunities and challenges in studying chromosome dynamics.

High-resolution single-cell 3D-models of chromatin ensembles during Drosophila embryogenesis

Single-cell chromatin studies provide insights into how chromatin structure relates to functions of individual cells. However, balancing high-resolution and genome wide-coverage remains challenging. We describe a computational method for the reconstruction of large 3D-ensembles of single-cell (sc) chromatin conformations from population Hi-C that we apply to study embryogenesis in Drosophila. With minimal assumptions of physical properties and without adjustable parameters, our method generates large ensembles of chromatin conformations via deep-sampling. Our method identifies specific interactions, which constitute 5-6% of Hi-C frequencies, but surprisingly are sufficient to drive chromatin folding, giving rise to the observed Hi-C patterns. Modeled sc-chromatins quantify chromatin heterogeneity, revealing significant changes during embryogenesis. Furthermore, >50% of modeled sc-chromatin maintain topologically associating domains (TADs) in early embryos, when no population TADs are perceptible. Domain boundaries become fixated during development, with strong preference at binding-sites of insulator-complexes upon the midblastula transition. Overall, high-resolution 3D-ensembles of sc-chromatin conformations enable further in-depth interpretation of population Hi-C, improving understanding of the structure-function relationship of genome organization.

Exploring chromosomal structural heterogeneity across multiple cell lines

Using computer simulations, we generate cell-specific 3D chromosomal structures and compare them to recently published chromatin structures obtained through microscopy. We demonstrate using machine learning and polymer physics simulations that epigenetic information can be used to predict the structural ensembles of multiple human cell lines. Theory predicts that chromosome structures are fluid and can only be described by an ensemble, which is consistent with the observation that chromosomes exhibit no unique fold. Nevertheless, our analysis of both structures from simulation and microscopy reveals that short segments of chromatin make two-state transitions between closed conformations and open dumbbell conformations. Finally, we study the conformational changes associated with the switching of genomic compartments observed in human cell lines. The formation of genomic compartments resembles hydrophobic collapse in protein folding, with the aggregation of denser and predominantly inactive chromatin driving the positioning of active chromatin toward the surface of individual chromosomal territories.