Computational 3D genome modeling using Chrom3D

Multiscale molecular modeling of chromatin with MultiMM: From nucleosomes to the whole genome

Motivation: We present a user-friendly 3D chromatin simulation model for the human genome based on OpenMM, addressing the challenges posed by existing models with use-specific implementations. Our approach employs a multi-scale energy minimization strategy, capturing chromatin's hierarchical structure. Initiating with a Hilbert curve-based structure, users can input files specifying nucleosome positioning, loops, compartments, or subcompartments. Results: The model utilizes an energy minimization approach with a large choice of numerical integrators, providing the entire genome's structure within minutes. Output files include the generated structures for each chromosome, offering a versatile and accessible tool for chromatin simulation in bioinformatics studies. Furthermore, MultiMM is capable of producing nucleosome-resolution structures by making simplistic geometric assumptions about the structure and the density of nucleosomes on the DNA. Code availability: Open-source software and the manual are freely available on https://github.com/SFGLab/MultiMM or via pip https://pypi.org/project/MultiMM/.

Loss of multi-level 3D genome organization during breast cancer progression

Breast cancer entails intricate alterations in genome organization and expression. However, how three-dimensional (3D) chromatin structure changes in the progression from a normal to a breast cancer malignant state remains unknown. To address this, we conducted an analysis combining Hi-C data with lamina-associated domains (LADs), epigenomic marks, and gene expression in an in vitro model of breast cancer progression. Our results reveal that while the fundamental properties of topologically associating domains (TADs) are overall maintained, significant changes occur in the organization of compartments and subcompartments. These changes are closely correlated with alterations in the expression of oncogenic genes. We also observe a restructuring of TAD-TAD interactions, coinciding with a loss of spatial compartmentalization and radial positioning of the 3D genome. Notably, we identify a previously unrecognized interchromosomal insertion event, wherein a locus on chromosome 8 housing the MYC oncogene is inserted into a highly active subcompartment on chromosome 10. This insertion is accompanied by the formation of de novo enhancer contacts and activation of MYC, illustrating how structural genomic variants can alter the 3D genome to drive oncogenic states. In summary, our findings provide evidence for the loss of genome organization at multiple scales during breast cancer progression revealing novel relationships between genome 3D structure and oncogenic processes.

SMCHD1 maintains heterochromatin and genome compartments in human myoblasts

Mammalian genomes are subdivided into euchromatic A compartments that contain mostly active chromatin, and inactive, heterochromatic B compartments. However, it is unknown how A and B genome compartments are established and maintained. Here we studied SMCHD1, an SMC-like protein in human male myoblasts. SMCHD1 colocalizes with Lamin B1 and the heterochromatin mark H3K9me3. Loss of SMCHD1 leads to extensive heterochromatin depletion at the nuclear lamina and acquisition of active chromatin states along all chromosomes. In absence of SMCHD1, long range intra-chromosomal and inter-chromosomal contacts between B compartments are lost while many new TADs and loops are formed. Inactivation of SMCHD1 promotes numerous B to A compartment transitions accompanied by activation of silenced genes. SMCHD1 functions as an anchor for heterochromatin domains ensuring that these domains are inaccessible to epigenome modification enzymes that typically operate in active chromatin. Therefore, A compartments are formed by default when not prevented by SMCHD1.

Active transcription and epigenetic reactions synergistically regulate meso-scale genomic organization

In interphase nuclei, chromatin forms dense domains of characteristic sizes, but the influence of transcription and histone modifications on domain size is not understood. We present a theoretical model exploring this relationship, considering chromatin-chromatin interactions, histone modifications, and chromatin extrusion. We predict that the size of heterochromatic domains is governed by a balance among the diffusive flux of methylated histones sustaining them and the acetylation reactions in the domains and the process of loop extrusion via supercoiling by RNAPII at their periphery, which contributes to size reduction. Super-resolution and nano-imaging of five distinct cell lines confirm the predictions indicating that the absence of transcription leads to larger heterochromatin domains. Furthermore, the model accurately reproduces the findings regarding how transcription-mediated supercoiling loss can mitigate the impacts of excessive cohesin loading. Our findings shed light on the role of transcription in genome organization, offering insights into chromatin dynamics and potential therapeutic targets.

Near telomere-to-telomere genome of the model plant Physcomitrium patens

The model plant Physcomitrium patens has played a pivotal role in enhancing our comprehension of plant evolution and development. However, the current genome harbours numerous regions that remain unfinished and erroneous. To address these issues, we generated an assembly using Oxford Nanopore reads and Hi-C mapping. The assembly incorporates telomeric and centromeric regions, thereby establishing it as a near telomere-to-telomere genome except a region in chromosome 1 that is not fully assembled due to its highly repetitive nature. This near telomere-to-telomere genome resolves the chromosome number at 26 and provides a gap-free genome assembly as well as updated gene models to aid future studies using this model organism.

DiffDomain enables identification of structurally reorganized topologically associating domains

Topologically associating domains (TADs) are critical structural units in three-dimensional genome organization of mammalian genome. Dynamic reorganizations of TADs between health and disease states are associated with essential genome functions. However, computational methods for identifying reorganized TADs are still in the early stages of development. Here, we present DiffDomain, an algorithm leveraging high-dimensional random matrix theory to identify structurally reorganized TADs using high-throughput chromosome conformation capture (Hi-C) contact maps. Method comparison using multiple real Hi-C datasets reveals that DiffDomain outperforms alternative methods for false positive rates, true positive rates, and identifying a new subtype of reorganized TADs. Applying DiffDomain to Hi-C data from different cell types and disease states demonstrates its biological relevance. Identified reorganized TADs are associated with structural variations and epigenomic changes such as changes in CTCF binding sites. By applying to a single-cell Hi-C data from mouse neuronal development, DiffDomain can identify reorganized TADs between cell types with reasonable reproducibility using pseudo-bulk Hi-C data from as few as 100 cells per condition. Moreover, DiffDomain reveals differential cell-to-population variability and heterogeneous cell-to-cell variability in TADs. Therefore, DiffDomain is a statistically sound method for better comparative analysis of TADs using both Hi-C and single-cell Hi-C data.

Large-scale data-driven and physics-based models offer insights into the relationships among the structures, dynamics, and functions of chromosomes

The organized three-dimensional chromosome architecture in the cell nucleus provides scaffolding for precise regulation of gene expression. When the cell changes its identity in the cell-fate decision-making process, extensive rearrangements of chromosome structures occur accompanied by large-scale adaptations of gene expression, underscoring the importance of chromosome dynamics in shaping genome function. Over the last two decades, rapid development of experimental methods has provided unprecedented data to characterize the hierarchical structures and dynamic properties of chromosomes. In parallel, these enormous data offer valuable opportunities for developing quantitative computational models. Here, we review a variety of large-scale polymer models developed to investigate the structures and dynamics of chromosomes. Different from the underlying modeling strategies, these approaches can be classified into data-driven ('top-down') and physics-based ('bottom-up') categories. We discuss their contributions to offering valuable insights into the relationships among the structures, dynamics, and functions of chromosomes and propose the perspective of developing data integration approaches from different experimental technologies and multidisciplinary theoretical/simulation methods combined with different modeling strategies.

Fundamental insights into the correlation between chromosome configuration and transcription

Eukaryotic chromosomes exhibit a hierarchical organization that spans a spectrum of length scales, ranging from sub-regions known as loops, which typically comprise hundreds of base pairs, to much larger chromosome territories that can encompass a few mega base pairs. Chromosome conformation capture experiments that involve high-throughput sequencing methods combined with microscopy techniques have enabled a new understanding of inter- and intra-chromosomal interactions with unprecedented details. This information also provides mechanistic insights on the relationship between genome architecture and gene expression. In this article, we review the recent findings on three-dimensional interactions among chromosomes at the compartment, topologically associating domain, and loop levels and the impact of these interactions on the transcription process. We also discuss current understanding of various biophysical processes involved in multi-layer structural organization of chromosomes. Then, we discuss the relationships between gene expression and genome structure from perturbative genome-wide association studies. Furthermore, for a better understanding of how chromosome architecture and function are linked, we emphasize the role of epigenetic modifications in the regulation of gene expression. Such an understanding of the relationship between genome architecture and gene expression can provide a new perspective on the range of potential future discoveries and therapeutic research.

True-to-scale DNA-density maps correlate with major accessibility differences between active and inactive chromatin

Chromatin compaction differences may have a strong impact on accessibility of individual macromolecules and macromolecular assemblies to their DNA target sites. Estimates based on fluorescence microscopy with conventional resolution, however, suggest only modest compaction differences (∼2-10×) between the active nuclear compartment (ANC) and inactive nuclear compartment (INC). Here, we present maps of nuclear landscapes with true-to-scale DNA densities, ranging from <5 to >300 Mbp/μm3. Maps are generated from individual human and mouse cell nuclei with single-molecule localization microscopy at ∼20 nm lateral and ∼100 nm axial optical resolution and are supplemented by electron spectroscopic imaging. Microinjection of fluorescent nanobeads with sizes corresponding to macromolecular assemblies for transcription into nuclei of living cells demonstrates their localization and movements within the ANC and exclusion from the INC.

Strong interactions between highly dynamic lamina-associated domains and the nuclear envelope stabilize the 3D architecture of Drosophila interphase chromatin

BACKGROUND

Interactions among topologically associating domains (TADs), and between the nuclear envelope (NE) and lamina-associated domains (LADs) are expected to shape various aspects of three-dimensional (3D) chromatin structure and dynamics; however, relevant genome-wide experiments that may provide statistically significant conclusions remain difficult.

RESULTS

We have developed a coarse-grained dynamical model of D. melanogaster nuclei at TAD resolution that explicitly accounts for four distinct epigenetic classes of TADs and LAD-NE interactions. The model is parameterized to reproduce the experimental Hi-C map of the wild type (WT) nuclei; it describes time evolution of the chromatin over the G1 phase of the interphase. The simulations include an ensemble of nuclei, corresponding to the experimentally observed set of several possible mutual arrangements of chromosomal arms. The model is validated against multiple structural features of chromatin from several different experiments not used in model development. Predicted positioning of all LADs at the NE is highly dynamic-the same LAD can attach, detach and move far away from the NE multiple times during interphase. The probabilities of LADs to be in contact with the NE vary by an order of magnitude, despite all having the same affinity to the NE in the model. These probabilities are mostly determined by a highly variable local linear density of LADs along the genome, which also has the same strong effect on the predicted positioning of individual TADs -- higher probability of a TAD to be near NE is largely determined by a higher linear density of LADs surrounding this TAD. The distribution of LADs along the chromosome chains plays a notable role in maintaining a non-random average global structure of chromatin. Relatively high affinity of LADs to the NE in the WT nuclei substantially reduces sensitivity of the global radial chromatin distribution to variations in the strength of TAD-TAD interactions compared to the lamin depleted nuclei, where a small (0.5 kT) increase of cross-type TAD-TAD interactions doubles the chromatin density in the central nucleus region.

CONCLUSIONS

A dynamical model of the entire fruit fly genome makes multiple genome-wide predictions of biological interest. The distribution of LADs along the chromatin chains affects their probabilities to be in contact with the NE and radial positioning of highly mobile TADs, playing a notable role in creating a non-random average global structure of the chromatin. We conjecture that an important role of attractive LAD-NE interactions is to stabilize global chromatin structure against inevitable cell-to-cell variations in TAD-TAD interactions.

The interplay of 3D genome organization with UV-induced DNA damage and repair

The 3D organization of the eukaryotic genome is crucial for various cellular processes such as gene expression and epigenetic regulation, as well as for maintaining genome integrity. However, the interplay between UV-induced DNA damage and repair with the 3D structure of the genome is not well understood. Here, we used state-of-the-art Hi-C, Damage-seq, and XR-seq datasets and in silico simulations to investigate the synergistic effects of UV damage and 3D genome organization. Our findings demonstrate that the peripheral 3D organization of the genome shields the central regions of genomic DNA from UV-induced damage. Additionally, we observed that potential damage sites of pyrimidine-pyrimidone (6-4) photoproducts are more prevalent in the nucleus center, possibly indicating an evolutionary pressure against those sites at the periphery. Interestingly, we found no correlation between repair efficiency and 3D structure after 12 min of irradiation, suggesting that UV radiation alters the genome's 3D organization in a short period of time. Interestingly, however, 2 h after UV induction, we observed more efficient repair levels in the center of the nucleus relative to the periphery. These results have implications for understanding the etiology of cancer and other diseases, as the interplay between UV radiation and the 3D genome may play a role in the development of genetic mutations and genomic instability.

Drosophila Models Reveal Properties of Mutant Lamins That Give Rise to Distinct Diseases

Mutations in the LMNA gene cause a collection of diseases known as laminopathies, including muscular dystrophies, lipodystrophies, and early-onset aging syndromes. The LMNA gene encodes A-type lamins, lamins A/C, intermediate filaments that form a meshwork underlying the inner nuclear membrane. Lamins have a conserved domain structure consisting of a head, coiled-coil rod, and C-terminal tail domain possessing an Ig-like fold. This study identified differences between two mutant lamins that cause distinct clinical diseases. One of the LMNA mutations encodes lamin A/C p.R527P and the other codes lamin A/C p.R482W, which are typically associated with muscular dystrophy and lipodystrophy, respectively. To determine how these mutations differentially affect muscle, we generated the equivalent mutations in the Drosophila Lamin C (LamC) gene, an orthologue of human LMNA. The muscle-specific expression of the R527P equivalent showed cytoplasmic aggregation of LamC, a reduced larval muscle size, decreased larval motility, and cardiac defects resulting in a reduced adult lifespan. By contrast, the muscle-specific expression of the R482W equivalent caused an abnormal nuclear shape without a change in larval muscle size, larval motility, and adult lifespan compared to controls. Collectively, these studies identified fundamental differences in the properties of mutant lamins that cause clinically distinct phenotypes, providing insights into disease mechanisms.

Robust 3D modeling reveals spatiosyntenic properties of animal genomes

Animal genomes are organized into chromosomes that are remarkably conserved in their gene content, forming distinct evolutionary units (synteny). Using versatile chromosomal modeling, we infer three-dimensional topology of genomes from representative clades spanning the earliest animal diversification. We apply a partitioning approach using interaction spheres to compensate for varying quality of topological data. Using comparative genomics approaches, we test whether syntenic signal at gene pair, local, and whole chromosomal scale is reflected in the reconstructed spatial organization. We identify evolutionarily conserved three-dimensional networks at all syntenic scales revealing novel evolutionarily maintained interactors associated with known conserved local gene linkages (such as hox). We thus present evidence for evolutionary constraints that are associated with three-, rather than just two-, dimensional animal genome organization, which we term spatiosynteny. As more accurate topological data become available, together with validation approaches, spatiosynteny may become relevant in understanding the functionality behind the observed conservation of animal chromosomes.

3D Genome Plasticity in Normal and Diseased Neurodevelopment

Non-random spatial organization of the chromosomal material inside the nuclei of brain cells emerges as an important regulatory layer of genome organization and function in health and disease. Here, we discuss how integrative approaches assessing chromatin in context of the 3D genome is providing new insights into normal and diseased neurodevelopment. Studies in primate (incl. human) and rodent brain have confirmed that chromosomal organization in neurons and glia undergoes highly dynamic changes during pre- and early postnatal development, with potential for plasticity across a much wider age window. For example, neuronal 3D genomes from juvenile and adult cerebral cortex and hippocampus undergo chromosomal conformation changes at hundreds of loci in the context of learning and environmental enrichment, viral infection, and neuroinflammation. Furthermore, locus-specific structural DNA variations, such as micro-deletions, duplications, repeat expansions, and retroelement insertions carry the potential to disrupt the broader epigenomic and transcriptional landscape far beyond the boundaries of the site-specific variation, highlighting the critical importance of long-range intra- and inter-chromosomal contacts for neuronal and glial function.

Genome organization in cardiomyocytes expressing mutated A-type lamins

Cardiomyopathy is a myocardial disorder, in which the heart muscle is structurally and functionally abnormal, often leading to heart failure. Dilated cardiomyopathy is characterized by a compromised left ventricular function and contributes significantly to the heart failure epidemic, which represents a staggering clinical and public health problem worldwide. Gene mutations have been identified in 35% of patients with dilated cardiomyopathy. Pathogenic variants in LMNA, encoding nuclear A-type lamins, are one of the major causative causes of dilated cardiomyopathy (i.e. CardioLaminopathy). A-type lamins are type V intermediate filament proteins, which are the main components of the nuclear lamina. The nuclear lamina is connected to the cytoskeleton on one side, and to the chromatin on the other side. Among the models proposed to explain how CardioLaminopathy arises, the “chromatin model” posits an effect of mutated A-type lamins on the 3D genome organization and thus on the transcription activity of tissue-specific genes. Chromatin contacts with the nuclear lamina via specific genomic regions called lamina-associated domains lamina-associated domains. These LADs play a role in the chromatin organization and gene expression regulation. This review focuses on the identification of LADs and chromatin remodeling in cardiac muscle cells expressing mutated A-type lamins and discusses the methods and relevance of these findings in disease.

Uncovering the Principles of Genome Folding by 3D Chromatin Modeling

Our understanding of how genomic DNA is tightly packed inside the nucleus, yet is still accessible for vital cellular processes, has grown dramatically over recent years with advances in microscopy and genomics technologies. Computational methods have played a pivotal role in the structural interpretation of experimental data, which helped unravel some organizational principles of genome folding. Here, we give an overview of current computational efforts in mechanistic and data-driven 3D chromatin structure modeling. We discuss strengths and limitations of different methods and evaluate the added value and benefits of computational approaches to infer the 3D structural and dynamic properties of the genome and its underlying mechanisms at different scales and resolution, ranging from the dynamic formation of chromatin loops and topological associated domains to nuclear compartmentalization of chromatin and nuclear bodies.

Biology and Model Predictions of the Dynamics and Heterogeneity of Chromatin-Nuclear Lamina Interactions

Associations of chromatin with the nuclear lamina, at the nuclear periphery, help shape the genome in 3 dimensions. The genomic landscape of lamina-associated domains (LADs) is well characterized, but much remains unknown on the physical and mechanistic properties of chromatin conformation at the nuclear lamina. Computational models of chromatin folding at, and interactions with, a surface representing the nuclear lamina are emerging in attempts to characterize these properties and predict chromatin behavior at the lamina in health and disease. Here, we highlight the heterogeneous nature of the nuclear lamina and LADs, outline the main 3-dimensional chromatin structural modeling methods, review applications of modeling chromatin-lamina interactions and discuss biological insights inferred from these models in normal and disease states. Lastly, we address perspectives on future developments in modeling chromatin interactions with the nuclear lamina.

Chromatin domain alterations linked to 3D genome organization in a large cohort of schizophrenia and bipolar disorder brains

Chromosomal organization, scaling from the 147-base pair (bp) nucleosome to megabase-ranging domains encompassing multiple transcriptional units, including heritability loci for psychiatric traits, remains largely unexplored in the human brain. In this study, we constructed promoter- and enhancer-enriched nucleosomal histone modification landscapes for adult prefrontal cortex from H3-lysine 27 acetylation and H3-lysine 4 trimethylation profiles, generated from 388 controls and 351 individuals diagnosed with schizophrenia (SCZ) or bipolar disorder (BD) (n = 739). We mapped thousands of cis-regulatory domains (CRDs), revealing fine-grained, 104-106-bp chromosomal organization, firmly integrated into Hi-C topologically associating domain stratification by open/repressive chromosomal environments and nuclear topography. Large clusters of hyper-acetylated CRDs were enriched for SCZ heritability, with prominent representation of regulatory sequences governing fetal development and glutamatergic neuron signaling. Therefore, SCZ and BD brains show coordinated dysregulation of risk-associated regulatory sequences assembled into kilobase- to megabase-scaling chromosomal domains.

Gene architecture directs splicing outcome in separate nuclear spatial regions

How the splicing machinery defines exons or introns as the spliced unit has remained a puzzle for 30 years. Here, we demonstrate that peripheral and central regions of the nucleus harbor genes with two distinct exon-intron GC content architectures that differ in the splicing outcome. Genes with low GC content exons, flanked by long introns with lower GC content, are localized in the periphery, and the exons are defined as the spliced unit. Alternative splicing of these genes results in exon skipping. In contrast, the nuclear center contains genes with a high GC content in the exons and short flanking introns. Most splicing of these genes occurs via intron definition, and aberrant splicing leads to intron retention. We demonstrate that the nuclear periphery and center generate different environments for the regulation of alternative splicing and that two sets of splicing factors form discrete regulatory subnetworks for the two gene architectures. Our study connects 3D genome organization and splicing, thus demonstrating that exon and intron definition modes of splicing occur in different nuclear regions.

Modeling the 3D Genome Using Hi-C and Nuclear Lamin-Genome Contacts

The three-dimensional (3D) organization of the genome is shaped by interactions with multiple structures within the nucleus, affecting gene expression outcomes. Technological breakthroughs in recent years have generated vast data reflecting various aspects of nuclear genome architecture in space and time. Integrating these datasets into comprehensive 3D genome models can reveal new insights into genome structure and regulation in normal and disease states. In this chapter, we provide a step-by-step guide on how to generate publication-ready integrated 3D genome models from (raw) Hi-C data and from lamin-genome (LAD) contacts.

Non-coding regulatory elements: Potential roles in disease and the case of epilepsy

Non-coding DNA (ncDNA) refers to the portion of the genome that does not code for proteins and accounts for the greatest physical proportion of the human genome. ncDNA includes sequences that are transcribed into RNA molecules, such as ribosomal RNAs (rRNAs), microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and un-transcribed sequences that have regulatory functions, including gene promoters and enhancers. Variation in non-coding regions of the genome have an established role in human disease, with growing evidence from many areas, including several cancers, Parkinson's disease and autism. Here, we review the features and functions of the regulatory elements that are present in the non-coding genome and the role that these regions have in human disease. We then review the existing research in epilepsy and emphasise the potential value of further exploring non-coding regulatory elements in epilepsy. In addition, we outline the most widely used techniques for recognising regulatory elements throughout the genome, current methodologies for investigating variation and the main challenges associated with research in the field of non-coding DNA.

Time-dependent effect of 1,6-hexanediol on biomolecular condensates and 3D chromatin organization

BACKGROUND

Biomolecular condensates have been implicated in multiple cellular processes. However, the global role played by condensates in 3D chromatin organization remains unclear. At present, 1,6-hexanediol (1,6-HD) is the only available tool to globally disrupt condensates, yet the conditions of 1,6-HD vary considerably between studies and may even trigger apoptosis.

RESULTS

In this study, we first analyzed the effects of different concentrations and treatment durations of 1,6-HD and found that short-term exposure to 1.5% 1,6-HD dissolved biomolecular condensates whereas long-term exposure caused aberrant aggregation without affecting cell viability. Based on this condition, we drew a time-resolved map of 3D chromatin organization and found that short-term treatment with 1.5% 1,6-HD resulted in reduced long-range interactions, strengthened compartmentalization, homogenized A-A interactions, B-to-A compartment switch and TAD reorganization, whereas longer exposure had the opposite effects. Furthermore, the long-range interactions between condensate-component-enriched regions were markedly weakened following 1,6-HD treatment.

CONCLUSIONS

In conclusion, our study finds a proper 1,6-HD condition and provides a resource for exploring the role of biomolecular condensates in 3D chromatin organization.

Genome-wide profiles of UV lesion susceptibility, repair, and mutagenic potential in melanoma

Exposure to the ultraviolet (UV) radiation in sunlight creates DNA lesions, which if left unrepaired can induce mutations and contribute to skin cancer. The two most common UV-induced DNA lesions are the cis-syn cyclobutane pyrimidine dimers (CPDs) and pyrimidine (6-4) pyrimidone photoproducts (6-4PPs), both of which can initiate mutations. Interestingly, mutation frequency across the genomes of many cancers is heterogenous with significant increases in heterochromatin. Corresponding increases in UV lesion susceptibility and decreases in repair are observed in heterochromatin versus euchromatin. However, the individual contributions of CPDs and 6-4PPs to mutagenesis have not been systematically examined in specific genomic and epigenomic contexts. In this study, we compared genome-wide maps of 6-4PP and CPD lesion abundances in primary cells and conducted comprehensive analyses to determine the genetic and epigenetic features associated with susceptibility. Overall, we found a high degree of similarity between 6-4PP and CPD formation, with an enrichment of both in heterochromatin regions. However, when examining the relative levels of the two UV lesions, we found that bivalent and Polycomb-repressed chromatin states were uniquely more susceptible to 6-4PPs. Interestingly, when comparing UV susceptibility and repair with melanoma mutation frequency in these regions, disparate patterns were observed in that susceptibility was not always inversely associated with repair and mutation frequency. Functional enrichment analysis hint at mechanisms of negative selection for these regions that are essential for cell viability, immune function and induce cell death when mutated. Ultimately, these results reveal both the similarities and differences between UV-induced lesions that contribute to melanoma.

Computational methods for predicting 3D genomic organization from high-resolution chromosome conformation capture data

The advent of high-resolution chromosome conformation capture assays (such as 5C, Hi-C and Pore-C) has allowed for unprecedented sequence-level investigations into the structure-function relationship of the genome. In order to comprehensively understand this relationship, computational tools are required that utilize data generated from these assays to predict 3D genome organization (the 3D genome reconstruction problem). Many computational tools have been developed that answer this need, but a comprehensive comparison of their underlying algorithmic approaches has not been conducted. This manuscript provides a comprehensive review of the existing computational tools (from November 2006 to September 2019, inclusive) that can be used to predict 3D genome organizations from high-resolution chromosome conformation capture data. Overall, existing tools were found to use a relatively small set of algorithms from one or more of the following categories: dimensionality reduction, graph/network theory, maximum likelihood estimation (MLE) and statistical modeling. Solutions in each category are far from maturity, and the breadth and depth of various algorithmic categories have not been fully explored. While the tools for predicting 3D structure for a genomic region or single chromosome are diverse, there is a general lack of algorithmic diversity among computational tools for predicting the complete 3D genome organization from high-resolution chromosome conformation capture data.

HMGB1 coordinates SASP-related chromatin folding and RNA homeostasis on the path to senescence

Spatial organization and gene expression of mammalian chromosomes are maintained and regulated in conjunction with cell cycle progression. This is perturbed once cells enter senescence and the highly abundant HMGB1 protein is depleted from nuclei to act as an extracellular proinflammatory stimulus. Despite its physiological importance, we know little about the positioning of HMGB1 on chromatin and its nuclear roles. To address this, we mapped HMGB1 binding genome-wide in two primary cell lines. We integrated ChIP-seq and Hi-C with graph theory to uncover clustering of HMGB1-marked topological domains that harbor genes involved in paracrine senescence. Using simplified Cross-Linking and Immuno-Precipitation and functional tests, we show that HMGB1 is also a bona fide RNA-binding protein (RBP) binding hundreds of mRNAs. It presents an interactome rich in RBPs implicated in senescence regulation. The mRNAs of many of these RBPs are directly bound by HMGB1 and regulate availability of SASP-relevant transcripts. Our findings reveal a broader than hitherto assumed role for HMGB1 in coordinating chromatin folding and RNA homeostasis as part of a regulatory loop controlling cell-autonomous and paracrine senescence.

Application of Hi-C and other omics data analysis in human cancer and cell differentiation research

With the development of 3C (chromosome conformation capture) and its derivative technology Hi-C (High-throughput chromosome conformation capture) research, the study of the spatial structure of the genomic sequence in the nucleus helps researchers understand the functions of biological processes such as gene transcription, replication, repair, and regulation. In this paper, we first introduce the research background and purpose of Hi-C data visualization analysis. After that, we discuss the Hi-C data analysis methods from genome 3D structure, A/B compartment, TADs (topologically associated domain), and loop detection. We also discuss how to apply genome visualization technologies to the identification of chromosome feature structures. We continue with a review of correlation analysis differences among multi-omics data, and how to apply Hi-C and other omics data analysis into cancer and cell differentiation research. Finally, we summarize the various problems in joint analyses based on Hi-C and other multi-omics data. We believe this review can help researchers better understand the progress and applications of 3D genome technology.

Polymer modelling unveils the roles of heterochromatin and nucleolar organizing regions in shaping 3D genome organization in Arabidopsis thaliana

The 3D genome is characterized by a complex organization made of genomic and epigenomic layers with profound implications on gene regulation and cell function. However, the understanding of the fundamental mechanisms driving the crosstalk between nuclear architecture and (epi)genomic information is still lacking. The plant Arabidopsis thaliana is a powerful model organism to address these questions owing to its compact genome for which we have a rich collection of microscopy, chromosome conformation capture (Hi-C) and ChIP-seq experiments. Using polymer modelling, we investigate the roles of nucleolus formation and epigenomics-driven interactions in shaping the 3D genome of A. thaliana. By validation of several predictions with published data, we demonstrate that self-attracting nucleolar organizing regions and repulsive constitutive heterochromatin are major mechanisms to regulate the organization of chromosomes. Simulations also suggest that interphase chromosomes maintain a partial structural memory of the V-shapes, typical of (sub)metacentric chromosomes in anaphase. Additionally, self-attraction between facultative heterochromatin regions facilitates the formation of Polycomb bodies hosting H3K27me3-enriched gene-clusters. Since nucleolus and heterochromatin are highly-conserved in eukaryotic cells, our findings pave the way for a comprehensive characterization of the generic principles that are likely to shape and regulate the 3D genome in many species.

Parsing the Functional Impact of Noncoding Genetic Variants in the Brain Epigenome

The heritability of common psychiatric disorders has motivated global efforts to identify risk-associated genetic variants and elucidate molecular pathways connecting DNA sequence to disease-associated brain dysfunction. The overrepresentation of risk variants among gene regulatory loci instead of protein-coding loci, however, poses a unique challenge in discerning which among the many thousands of variants identified contribute functionally to disease etiology. Defined broadly, psychiatric epigenomics seeks to understand the effects of disease-associated genetic variation on functional readouts of chromatin in an effort to prioritize variants in terms of their impact on gene expression in the brain. Here, we provide an overview of epigenomic mapping in the human brain and highlight findings of particular relevance to psychiatric genetics. Computational methods, including convolutional neuronal networks, and other machine learning approaches hold great promise for elucidating the functional impact of both common and rare genetic variants, thereby refining the epigenomic architecture of psychiatric disorders and enabling integrative analyses of regulatory noncoding variants in the context of large population-level genome and phenome databases.

Hi-C implementation of genome structure for in silico models of radiation-induced DNA damage

Developments in the genome organisation field has resulted in the recent methodology to infer spatial conformations of the genome directly from experimentally measured genome contacts (Hi-C data). This provides a detailed description of both intra- and inter-chromosomal arrangements. Chromosomal intermingling is an important driver for radiation-induced DNA mis-repair. Which is a key biological endpoint of relevance to the fields of cancer therapy (radiotherapy), public health (biodosimetry) and space travel. For the first time, we leverage these methods of inferring genome organisation and couple them to nano-dosimetric radiation track structure modelling to predict quantities and distribution of DNA damage within cell-type specific geometries. These nano-dosimetric simulations are highly dependent on geometry and are benefited from the inclusion of experimentally driven chromosome conformations. We show how the changes in Hi-C contract maps impact the inferred geometries resulting in significant differences in chromosomal intermingling. We demonstrate how these differences propagate through to significant changes in the distribution of DNA damage throughout the cell nucleus, suggesting implications for DNA repair fidelity and subsequent cell fate. We suggest that differences in the geometric clustering for the chromosomes between the cell-types are a plausible factor leading to changes in cellular radiosensitivity. Furthermore, we investigate changes in cell shape, such as flattening, and show that this greatly impacts the distribution of DNA damage. This should be considered when comparing in vitro results to in vivo systems. The effect may be especially important when attempting to translate radiosensitivity measurements at the experimental in vitro level to the patient or human level.

The Interchromatin Compartment Participates in the Structural and Functional Organization of the Cell Nucleus

This article focuses on the role of the interchromatin compartment (IC) in shaping nuclear landscapes. The IC is connected with nuclear pore complexes (NPCs) and harbors splicing speckles and nuclear bodies. It is postulated that the IC provides routes for imported transcription factors to target sites, for export routes of mRNA as ribonucleoproteins toward NPCs, as well as for the intranuclear passage of regulatory RNAs from sites of transcription to remote functional sites (IC hypothesis). IC channels are lined by less-compacted euchromatin, called the perichromatin region (PR). The PR and IC together form the active nuclear compartment (ANC). The ANC is co-aligned with the inactive nuclear compartment (INC), comprising more compacted heterochromatin. It is postulated that the INC is accessible for individual transcription factors, but inaccessible for larger macromolecular aggregates (limited accessibility hypothesis). This functional nuclear organization depends on still unexplored movements of genes and regulatory sequences between the two compartments.

Interplay of lamin A and lamin B LADs on the radial positioning of chromatin

Immunosuppressive drugs such as cyclosporin A (CsA) can elicit hepatotoxicity by affecting gene expression. Here, we address the link between CsA and large-scale chromatin organization in HepG2 hepatocarcinoma cells. We show the existence of lamina-associated domains (LADs) interacting with lamin A, lamin B, or both. These ‘A-B’, ‘A-only’ and ‘B-only’ LADs display distinct fates after CsA treatment: A-B LADs remain constitutive or lose A, A-only LADs mainly lose A or switch to B, and B-only LADs remain B-only or acquire A. LAD rearrangement is overall uncoupled from changes in gene expression. Three-dimensional (3D) genome modeling predicts changes in radial positioning of LADs as LADs switch identities, which are corroborated by fluorescence in situ hybridization. Our results reveal interplay between A- and B-type lamins on radial locus positioning, suggesting complementary contributions to large-scale genome architecture. The data also unveil a hitherto unsuspected impact of cytotoxic drugs on genome conformation.Abbreviations: ChIP-seq: chromatin immunoprecipitation sequencing; CsA: cyclosporin A; FISH; fluorescence in situ hybridization; ICMT: isoprenylcysteine methyltransferase; LAD: lamina-associated domain; TAD: topologically-associated domain

Computing 3D Chromatin Configurations from Contact Probability Maps by Inverse Brownian Dynamics

The three-dimensional (3D) organization of chromatin, on the length scale of a few genes, is crucial in determining the functional state-accessibility and amount of gene expression-of the chromatin. Recent advances in chromosome conformation capture experiments provide partial information on the chromatin organization in a cell population, namely the contact count between any segment pairs, but not on the interaction strength that leads to these contact counts. However, given the contact matrix, determining the complete 3D organization of the whole chromatin polymer is an inverse problem. In this work, a novel inverse Brownian dynamics method based on a coarse-grained bead-spring chain model has been proposed to compute the optimal interaction strengths between different segments of chromatin such that the experimentally measured contact count probability constraints are satisfied. Applying this method to the α-globin gene locus in two different cell types, we predict the 3D organizations corresponding to active and repressed states of chromatin at the locus. We show that the average distance between any two segments of the region has a broad distribution and cannot be computed as a simple inverse relation based on the contact probability alone. The results presented for multiple normalization methods suggest that all measurable quantities may crucially depend on the nature of normalization. We argue that by experimentally measuring predicted quantities, one may infer the appropriate form of normalization.

A chromosomal connectome for psychiatric and metabolic risk variants in adult dopaminergic neurons

BACKGROUND

Midbrain dopaminergic neurons (MDN) represent 0.0005% of the brain's neuronal population and mediate cognition, food intake, and metabolism. MDN are also posited to underlay the neurobiological dysfunction of schizophrenia (SCZ), a severe neuropsychiatric disorder that is characterized by psychosis as well as multifactorial medical co-morbidities, including metabolic disease, contributing to markedly increased morbidity and mortality. Paradoxically, however, the genetic risk sequences of psychosis and traits associated with metabolic disease, such as body mass, show very limited overlap.

METHODS

We investigated the genomic interaction of SCZ with medical conditions and traits, including body mass index (BMI), by exploring the MDN's "spatial genome," including chromosomal contact landscapes as a critical layer of cell type-specific epigenomic regulation. Low-input Hi-C protocols were applied to 5-10 × 103 dopaminergic and other cell-specific nuclei collected by fluorescence-activated nuclei sorting from the adult human midbrain.

RESULTS

The Hi-C-reconstructed MDN spatial genome revealed 11 "Euclidean hot spots" of clustered chromatin domains harboring risk sequences for SCZ and elevated BMI. Inter- and intra-chromosomal contacts interconnecting SCZ and BMI risk sequences showed massive enrichment for brain-specific expression quantitative trait loci (eQTL), with gene ontologies, regulatory motifs and proteomic interactions related to adipogenesis and lipid regulation, dopaminergic neurogenesis and neuronal connectivity, and reward- and addiction-related pathways.

CONCLUSIONS

We uncovered shared nuclear topographies of cognitive and metabolic risk variants. More broadly, our PsychENCODE sponsored Hi-C study offers a novel genomic approach for the study of psychiatric and medical co-morbidities constrained by limited overlap of their respective genetic risk architectures on the linear genome.

MOCHI enables discovery of heterogeneous interactome modules in 3D nucleome

The composition of the cell nucleus is highly heterogeneous, with different constituents forming complex interactomes. However, the global patterns of these interwoven heterogeneous interactomes remain poorly understood. Here we focus on two different interactomes, chromatin interaction network and gene regulatory network, as a proof of principle to identify heterogeneous interactome modules (HIMs), each of which represents a cluster of gene loci that is in spatial contact more frequently than expected and that is regulated by the same group of transcription factors. HIM integrates transcription factor binding and 3D genome structure to reflect “transcriptional niche” in the nucleus. We develop a new algorithm, MOCHI, to facilitate the discovery of HIMs based on network motif clustering in heterogeneous interactomes. By applying MOCHI to five different cell types, we found that HIMs have strong spatial preference within the nucleus and show distinct functional properties. Through integrative analysis, this work shows the utility of MOCHI to identify HIMs, which may provide new perspectives on the interplay between transcriptional regulation and 3D genome organization.

Finding Friends in the Crowd: Three-Dimensional Cliques of Topological Genomic Domains

The mammalian genome is intricately folded in a three-dimensional topology believed to be important for the orchestration of gene expression regulating development, differentiation and tissue homeostasis. Important features of spatial genome conformation in the nucleus are promoter-enhancer contacts regulating gene expression within topologically-associated domains (TADs), short- and long-range interactions between TADs and associations of chromatin with nucleoli and nuclear speckles. In addition, anchoring of chromosomes to the nuclear lamina via lamina-associated domains (LADs) at the nuclear periphery is a key regulator of the radial distribution of chromatin. To what extent TADs and LADs act in concert as genomic organizers to shape the three-dimensional topology of chromatin has long remained unknown. A new study addressing this key question provides evidence of (i) preferred long-range associations between TADs forming TAD “cliques” which organize large heterochromatin domains, and (ii) stabilization of TAD cliques by LADs at the nuclear periphery after induction of terminal differentiation. Here, we review these findings, address the issue of whether TAD cliques exist in single cells and discuss the extent of cell-to-cell heterogeneity in higher-order chromatin conformation. The recent observations provide a first appreciation of changes in 4-dimensional higher-order genome topologies during differentiation.

Long-range interactions between topologically associating domains shape the four-dimensional genome during differentiation

Genomic information is selectively used to direct spatial and temporal gene expression during differentiation. Interactions between topologically associating domains (TADs) and between chromatin and the nuclear lamina organize and position chromosomes in the nucleus. However, how these genomic organizers together shape genome architecture is unclear. Here, using a dual-lineage differentiation system, we report long-range TAD-TAD interactions that form constitutive and variable TAD cliques. A differentiation-coupled relationship between TAD cliques and lamina-associated domains suggests that TAD cliques stabilize heterochromatin at the nuclear periphery. We also provide evidence of dynamic TAD cliques during mouse embryonic stem-cell differentiation and somatic cell reprogramming and of inter-TAD associations in single-cell high-resolution chromosome conformation capture (Hi-C) data. TAD cliques represent a level of four-dimensional genome conformation that reinforces the silencing of repressed developmental genes.

Lamin A, Chromatin and FPLD2: Not Just a Peripheral Ménage-à-Trois

At the nuclear periphery, the genome is anchored to A- and B-type nuclear lamins in the form of heterochromatic lamina-associated domains. A-type lamins also associate with chromatin in the nuclear interior, away from the peripheral nuclear lamina. This nucleoplasmic lamin A environment tends to be euchromatic, suggesting distinct roles of lamin A in the regulation of gene expression in peripheral and more central regions of the nucleus. The hot-spot lamin A R482W mutation causing familial partial lipodystrophy of Dunnigan-type (FPLD2), affects lamin A association with chromatin at the nuclear periphery and in the nuclear interior, and is associated with 3-dimensional (3D) rearrangements of chromatin. Here, we highlight features of nuclear lamin association with the genome at the nuclear periphery and in the nuclear interior. We address recent data showing a rewiring of such interactions in cells from FPLD2 patients, and in adipose progenitor and induced pluripotent stem cell models of FPLD2. We discuss associated epigenetic and genome conformation changes elicited by the lamin A R482W mutation at the gene level. The findings argue that the mutation adversely impacts both global and local genome architecture throughout the nucleus space. The results, together with emerging new computational modeling tools, mark the start of a new era in our understanding of the 3D genomics of laminopathies.