The role of 3D chromatin domains in gene regulation: a multi-facetted view on genome organization

Unveiling chromatin dynamics with virtual epigenome

The three-dimensional organization of chromatin is essential for gene regulation and cellular function, with epigenome playing a key role. Hi-C methods have expanded our understanding of chromatin interactions, but their high cost and complexity limit their use. Existing models for predicting chromatin interactions rely on limited ChIP-seq inputs, reducing their accuracy and generalizability. In this work, we present a computational approach, EpiVerse, which leverages imputed epigenetic signals and advanced deep learning techniques. EpiVerse significantly improves the accuracy of cross-cell-type Hi-C prediction, while also enhancing model interpretability by incorporating chromatin state prediction within a multitask learning framework. Moreover, EpiVerse predicts Hi-C contact maps across an array of 39 human tissues, which provides a comprehensive view of the complex relationship between chromatin structure and gene regulation. Furthermore, EpiVerse facilitates unprecedented in silico perturbation experiments at the "epigenome-level" to unveil the chromatin architecture under specific conditions. EpiVerse is available on GitHub: https://github.com/jhhung/EpiVerse .

Tracing the Shared Foundations of Gene Expression and Chromatin Structure

The three-dimensional organization of chromatin into topologically associating domains (TADs) may impact gene regulation by bringing distant genes into contact. However, many questions about TADs' function and their influence on transcription remain unresolved due to technical limitations in defining TAD boundaries and measuring the direct effect that TADs have on gene expression. Here, we develop consensus TAD maps for human and mouse with a novel "bag-of-genes" approach for defining the gene composition within TADs. This approach enables new functional interpretations of TADs by providing a way to capture species-level differences in chromatin organization. We also leverage a generative AI foundation model computed from 33 million transcriptomes to define contextual similarity, an embedding-based metric that is more powerful than co-expression at representing functional gene relationships. Our analytical framework directly leads to testable hypotheses about chromatin organization across cellular states. We find that TADs play an active role in facilitating gene co-regulation, possibly through a mechanism involving transcriptional condensates. We also discover that the TAD-linked enhancement of transcriptional context is strongest in early developmental stages and systematically declines with aging. Investigation of cancer cells show distinct patterns of TAD usage that shift with chemotherapy treatment, suggesting specific roles for TAD-mediated regulation in cellular development and plasticity. Finally, we develop "TAD signatures" to improve statistical analysis of single-cell transcriptomic data sets in predicting cancer cell-line drug response. These findings reshape our understanding of cellular plasticity in development and disease, indicating that chromatin organization acts through probabilistic mechanisms rather than deterministic rules.

Software availability

https://singhlab.net/tadmap.

Loop Catalog: a comprehensive HiChIP database of human and mouse samples

HiChIP enables cost-effective and high-resolution profiling of chromatin loops. To leverage the increasing number of HiChIP datasets, we developed Loop Catalog (https://loopcatalog.lji.org), a web-based database featuring loop calls from 1000+ distinct human and mouse HiChIP samples from 152 studies plus 44 high-resolution Hi-C samples. We demonstrate its utility for interpreting GWAS and eQTL variants through SNP-to-gene linking, identifying enriched sequence motifs and motif pairs, and generating regulatory networks and 2D representations of chromatin structure. Our catalog spans over 4.19M unique loops, and with embedded analysis modules, constitutes an important resource for the field.

Interphase chromatin biophysics and mechanics: new perspectives and open questions

The physical organization and properties of chromatin within the interphase nucleus are intimately linked to a wide range of functional DNA-based processes. In this context, interphase chromatin mechanics - that is, how chromatin, physically, responds to forces - is gaining increasing attention. Recent methodological advances for probing the force-response of chromatin in cellulo open new avenues for research, as well as new questions. This review discusses emerging views from these approaches and others, including recent in vitro single-molecule studies of cohesin and condensin motor activities, providing insights into physical and material aspects of chromatin, its plasticity in the context of functional processes, its nonequilibrium or 'active matter' properties, and the importance of factors such as chromatin fiber tension and stiffness. This growing field offers exciting opportunities to better understand the interplay between interphase chromosome structure, dynamics, mechanics, and functions.

Sequence-Based Machine Learning Reveals 3D Genome Differences between Bonobos and Chimpanzees

The 3D structure of the genome is an important mediator of gene expression. As phenotypic divergence is largely driven by gene regulatory variation, comparing genome 3D contacts across species can further understanding of the molecular basis of species differences. However, while experimental data on genome 3D contacts in humans are increasingly abundant, only a handful of 3D genome contact maps exist for other species. Here, we demonstrate that human experimental data can be used to close this data gap. We apply a machine learning model that predicts 3D genome contacts from DNA sequence to the genomes from 56 bonobos and chimpanzees and identify species-specific patterns of genome folding. We estimated 3D divergence between individuals from the resulting contact maps in 4,420 1 Mb genomic windows, of which ∼17% were substantially divergent in predicted genome contacts. Bonobos and chimpanzees diverged at 89 windows, overlapping genes associated with multiple traits implicated in Pan phenotypic divergence. We discovered 51 bonobo-specific variants that individually produce the observed bonobo contact pattern in bonobo-chimpanzee divergent windows. Our results demonstrate that machine learning methods can leverage human data to fill in data gaps across species, offering the first look at population-level 3D genome variation in nonhuman primates. We also identify loci where changes in 3D folding may contribute to phenotypic differences in our closest living relatives.

ENT3C: an entropy-based similarity measure for Hi-C and micro-C derived contact matrices

Hi-C and micro-C sequencing have shed light on the profound importance of 3D genome organization in cellular function by probing 3D contact frequencies across the linear genome. The resulting contact matrices are extremely sparse and susceptible to technical- and sequence-based biases, making their comparison challenging. The development of reliable, robust and efficient methods for quantifying similarity between contact matrices is crucial for investigating variations in the 3D genome organization in different cell types or under different conditions, as well as evaluating experimental reproducibility. We present a novel method, ENT3C, which measures the change in pattern complexity in the vicinity of contact matrix diagonals to quantify their similarity. ENT3C provides a robust, user-friendly Hi-C or micro-C contact matrix similarity metric and a characteristic entropy signal that can be used to gain detailed biological insights into 3D genome organization.

Establishment and maintenance of random monoallelic expression

This Review elucidates the regulatory principles of random monoallelic expression by focusing on two well-studied examples: the X-chromosome inactivation regulator Xist and the olfactory receptor gene family. Although the choice of a single X chromosome or olfactory receptor occurs in different developmental contexts, common gene regulatory principles guide monoallelic expression in both systems. In both cases, an event breaks the symmetry between genetically and epigenetically identical copies of the gene, leading to the expression of one single random allele, stabilized through negative feedback control. Although many regulatory steps that govern the establishment and maintenance of monoallelic expression have been identified, key pieces of the puzzle are still missing. We provide an overview of the current knowledge and models for the monoallelic expression of Xist and olfactory receptors. We discuss their similarities and differences, and highlight open questions and approaches that could guide the study of other monoallelically expressed genes.

3D epigenomics and 3D epigenopathies

Mammalian genomes are intricately compacted to form sophisticated 3-dimensional structures within the tiny nucleus, so called 3D genome folding. Despite their shapes reminiscent of an entangled yarn, the rapid development of molecular and next-generation sequencing technologies (NGS) has revealed that mammalian genomes are highly organized in a hierarchical order that delicately affects transcription activities. An increasing amount of evidence suggests that 3D genome folding is implicated in diseases, giving us a clue on how to identify novel therapeutic approaches. In this review, we will study what 3D genome folding means in epigenetics, what types of 3D genome structures there are, how they are formed, and how the technologies have developed to explore them. We will also discuss the pathological implications of 3D genome folding. Finally, we will discuss how to leverage 3D genome folding and engineering for future studies. [BMB Reports 2024; 57(5): 216-231].

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

Boundary stacking interactions enable cross-TAD enhancer-promoter communication during limb development

Although promoters and their enhancers are frequently contained within a topologically associating domain (TAD), some developmentally important genes have their promoter and enhancers within different TADs. Hypotheses about molecular mechanisms enabling cross-TAD interactions remain to be assessed. To test these hypotheses, we used optical reconstruction of chromatin architecture to characterize the conformations of the Pitx1 locus on single chromosomes in developing mouse limbs. Our data support a model in which neighboring boundaries are stacked as a result of loop extrusion, bringing boundary-proximal cis-elements into contact. This stacking interaction also contributes to the appearance of architectural stripes in the population average maps. Through molecular dynamics simulations, we found that increasing boundary strengths facilitates the formation of the stacked boundary conformation, counter-intuitively facilitating border bypass. This work provides a revised view of the TAD borders' function, both facilitating and preventing cis-regulatory interactions, and introduces a framework to distinguish border-crossing from border-respecting enhancer-promoter pairs.

Accessible gene borders establish a core structural unit for chromatin architecture in Arabidopsis

Three-dimensional (3D) chromatin structure is linked to transcriptional regulation in multicellular eukaryotes including plants. Taking advantage of high-resolution Hi-C (high-throughput chromatin conformation capture), we detected a small structural unit with 3D chromatin architecture in the Arabidopsis genome, which lacks topologically associating domains, and also in the genomes of tomato, maize, and Marchantia polymorpha. The 3D folding domain unit was usually established around an individual gene and was dependent on chromatin accessibility at the transcription start site (TSS) and transcription end site (TES). We also observed larger contact domains containing two or more neighboring genes, which were dependent on accessible border regions. Binding of transcription factors to accessible TSS/TES regions formed these gene domains. We successfully simulated these Hi-C contact maps via computational modeling using chromatin accessibility as input. Our results demonstrate that gene domains establish basic 3D chromatin architecture units that likely contribute to higher-order 3D genome folding in plants.

Sequence-based machine learning reveals 3D genome differences between bonobos and chimpanzees

Phenotypic divergence between closely related species, including bonobos and chimpanzees (genus Pan), is largely driven by variation in gene regulation. The 3D structure of the genome mediates gene expression; however, genome folding differences in Pan are not well understood. Here, we apply machine learning to predict genome-wide 3D genome contact maps from DNA sequence for 56 bonobos and chimpanzees, encompassing all five extant lineages. We use a pairwise approach to estimate 3D divergence between individuals from the resulting contact maps in 4,420 1 Mb genomic windows. While most pairs were similar, ∼17% were predicted to be substantially divergent in genome folding. The most dissimilar maps were largely driven by single individuals with rare variants that produce unique 3D genome folding in a region. We also identified 89 genomic windows where bonobo and chimpanzee contact maps substantially diverged, including several windows harboring genes associated with traits implicated in Pan phenotypic divergence. We used in silico mutagenesis to identify 51 3D-modifying variants in these bonobo-chimpanzee divergent windows, finding that 34 or 66.67% induce genome folding changes via CTCF binding motif disruption. Our results reveal 3D genome variation at the population-level and identify genomic regions where changes in 3D folding may contribute to phenotypic differences in our closest living relatives.

CEBPA phase separation links transcriptional activity and 3D chromatin hubs

Cell identity is orchestrated through an interplay between transcription factor (TF) action and genome architecture. The mechanisms used by TFs to shape three-dimensional (3D) genome organization remain incompletely understood. Here we present evidence that the lineage-instructive TF CEBPA drives extensive chromatin compartment switching and promotes the formation of long-range chromatin hubs during induced B cell-to-macrophage transdifferentiation. Mechanistically, we find that the intrinsically disordered region (IDR) of CEBPA undergoes in vitro phase separation (PS) dependent on aromatic residues. Both overexpressing B cells and native CEBPA-expressing cell types such as primary granulocyte-macrophage progenitors, liver cells, and trophectoderm cells reveal nuclear CEBPA foci and long-range 3D chromatin hubs at CEBPA-bound regions. In short, we show that CEBPA can undergo PS through its IDR, which may underlie in vivo foci formation and suggest a potential role of PS in regulating CEBPA function.

Shifting landscapes: the role of 3D genomic organizations in gene regulatory strategies

3D genome folding enables the physical storage of chromosomes into the compact volume of a cell's nucleus, allows for the accurate segregation of chromatin to daughter cells, and has been shown to be tightly coupled to the way in which genetic information is converted into transcriptional programs [1-3]. Importantly, this link between chromatin architecture and gene regulation is a selectable feature in which modifications to chromatin organization accompany, or perhaps even drive the establishment of new regulatory strategies with enduring impacts on animal body plan complexity. Here, we discuss the nature of different 3D genome folding systems found across the tree of life, with particular emphasis on metazoans, and the relative influence of these systems on gene regulation. We suggest how the properties of these folding systems have influenced regulatory strategies employed by different lineages and may have catalyzed the partitioning and specialization of genetic programs that enabled multicellularity and organ-grade body plan complexity.

DeDoc2 Identifies and Characterizes the Hierarchy and Dynamics of Chromatin TAD-Like Domains in the Single Cells

Topologically associating domains (TADs) are functional chromatin units with hierarchical structure. However, the existence, prevalence, and dynamics of such hierarchy in single cells remain unexplored. Here, a new generation TAD-like domain (TLD) detection algorithm, named deDoc2, to decode the hierarchy of TLDs in single cells, is reported. With dynamic programming, deDoc2 seeks genome partitions with global minimal structure entropy for both whole and local contact matrix. Notably, deDoc2 outperforms state-of-the-art tools and is one of only two tools able to identify the hierarchy of TLDs in single cells. By applying deDoc2, it is showed that the hierarchy of TLDs in single cells is highly dynamic during cell cycle, as well as among human brain cortex cells, and that it is associated with cellular identity and functions. Thus, the results demonstrate the abundance of information potentially encoded by TLD hierarchy for functional regulation. The deDoc2 can be freely accessed at https://github.com/zengguangjie/deDoc2.

Congenital aniridia beyond black eyes: From phenotype and novel genetic mechanisms to innovative therapeutic approaches

Congenital PAX6-aniridia, initially characterized by the absence of the iris, has progressively been shown to be associated with other developmental ocular abnormalities and systemic features making congenital aniridia a complex syndromic disorder rather than a simple isolated disease of the iris. Moreover, foveal hypoplasia is now recognized as a more frequent feature than complete iris hypoplasia and a major visual prognosis determinant, reversing the classical clinical picture of this disease. Conversely, iris malformation is also a feature of various anterior segment dysgenesis disorders caused by PAX6-related developmental genes, adding a level of genetic complexity for accurate molecular diagnosis of aniridia. Therefore, the clinical recognition and differential genetic diagnosis of PAX6-related aniridia has been revealed to be much more challenging than initially thought, and still remains under-investigated. Here, we update specific clinical features of aniridia, with emphasis on their genotype correlations, as well as provide new knowledge regarding the PAX6 gene and its mutational spectrum, and highlight the beneficial utility of clinically implementing targeted Next-Generation Sequencing combined with Whole-Genome Sequencing to increase the genetic diagnostic yield of aniridia. We also present new molecular mechanisms underlying aniridia and aniridia-like phenotypes. Finally, we discuss the appropriate medical and surgical management of aniridic eyes, as well as innovative therapeutic options. Altogether, these combined clinical-genetic approaches will help to accelerate time to diagnosis, provide better determination of the disease prognosis and management, and confirm eligibility for future clinical trials or genetic-specific therapies.

Loops and the activity of loop extrusion factors constrain chromatin dynamics

The chromosomes-DNA polymers and their binding proteins-are compacted into a spatially organized, yet dynamic, three-dimensional structure. Recent genome-wide chromatin conformation capture experiments reveal a hierarchical organization of the DNA structure that is imposed, at least in part, by looping interactions arising from the activity of loop extrusion factors. The dynamics of chromatin reflects the response of the polymer to a combination of thermal fluctuations and active processes. However, how chromosome structure and enzymes acting on chromatin together define its dynamics remains poorly understood. To gain insight into the structure-dynamics relationship of chromatin, we combine high-precision microscopy in living Schizosaccharomyces pombe cells with systematic genetic perturbations and Rouse model polymer simulations. We first investigated how the activity of two loop extrusion factors, the cohesin and condensin complexes, influences chromatin dynamics. We observed that deactivating cohesin, or to a lesser extent condensin, increased chromatin mobility, suggesting that loop extrusion constrains rather than agitates chromatin motion. Our corresponding simulations reveal that the introduction of loops is sufficient to explain the constraining activity of loop extrusion factors, highlighting that the conformation adopted by the polymer plays a key role in defining its dynamics. Moreover, we find that the number of loops or residence times of loop extrusion factors influence the dynamic behavior of the chromatin polymer. Last, we observe that the activity of the INO80 chromatin remodeler, but not the SWI/SNF or RSC complexes, is critical for ATP-dependent chromatin mobility in fission yeast. Taking the data together, we suggest that thermal and INO80-dependent activities exert forces that drive chromatin fluctuations, which are constrained by the organization of the chromosome into loops.

When 3D genome changes cause disease: the impact of structural variations in congenital disease and cancer

Large structural variations (SV) are a class of mutations that have long been known to cause a wide range of genetic diseases, from rare congenital disease to cancer. Many of these SVs do not directly disrupt disease-related genes and determining causal genotype-phenotype relationships has been challenging to disentangle in the past. This has started to change with our increased understanding of the 3D genome folding. The pathophysiologies of the different types of genetic diseases influence the type of SVs observed and their genetic consequences, and how these are connected to 3D genome folding. We propose guiding principles for interpreting disease-associated SVs based on our current understanding of 3D chromatin architecture and the gene-regulatory and physiological mechanisms disrupted in disease.

Super-enhancer landscape rewiring in cancer: The epigenetic control at distal sites

Super-enhancers evolve as elements at the top of the hierarchical control of gene expression. They are important end-gatherers of signaling pathways that control stemness, differentiation or adaptive responses. Many epigenetic regulations focus on these regions, and not surprisingly, during the process of tumorigenesis, various alterations can account for their dysfunction. Super-enhancers are emerging as key drivers of the aberrant gene expression landscape that sustain the aggressiveness of cancer cells. In this review, we will describe and discuss about the structure of super-enhancers, their epigenetic regulation, and the major changes affecting their functionality in cancer.

Recent progress and challenges in single-cell imaging of enhancer-promoter interaction

In the past two years, approaches relying on high-resolution microscopy and live-cell imaging have increasingly contributed to our understanding of the 3D genome organization and its importance for transcriptional control. Here, we describe recent progress that has highlighted how flexible and heterogeneous 3D chromatin structure is, on the length scales relevant to transcriptional control. We describe work that has investigated how robust transcriptional outcomes may be derived from such flexible organization without the need for clearly distinct structures in active and silent cells. We survey the latest state of the art in directly observing the dynamics of chromatin interactions, and suggest how some recent, apparently contradictory conclusions may be reconciled.

Function and Evolution of the Loop Extrusion Machinery in Animals

Structural maintenance of chromosomes (SMC) complexes are essential proteins found in genomes of all cellular organisms. Essential functions of these proteins, such as mitotic chromosome formation and sister chromatid cohesion, were discovered a long time ago. Recent advances in chromatin biology showed that SMC proteins are involved in many other genomic processes, acting as active motors extruding DNA, which leads to the formation of chromatin loops. Some loops formed by SMC proteins are highly cell type and developmental stage specific, such as SMC-mediated DNA loops required for VDJ recombination in B-cell progenitors, or dosage compensation in Caenorhabditis elegans and X-chromosome inactivation in mice. In this review, we focus on the extrusion-based mechanisms that are common for multiple cell types and species. We will first describe an anatomy of SMC complexes and their accessory proteins. Next, we provide biochemical details of the extrusion process. We follow this by the sections describing the role of SMC complexes in gene regulation, DNA repair, and chromatin topology.

The spatial organization of transcriptional control

In animals, the sequences for controlling gene expression do not concentrate just at the transcription start site of genes, but are frequently thousands to millions of base pairs distal to it. The interaction of these sequences with one another and their transcription start sites is regulated by factors that shape the three-dimensional (3D) organization of the genome within the nucleus. Over the past decade, indirect tools exploiting high-throughput DNA sequencing have helped to map this 3D organization, have identified multiple key regulators of its structure and, in the process, have substantially reshaped our view of how 3D genome architecture regulates transcription. Now, new tools for high-throughput super-resolution imaging of chromatin have directly visualized the 3D chromatin organization, settling some debates left unresolved by earlier indirect methods, challenging some earlier models of regulatory specificity and creating hypotheses about the role of chromatin structure in transcriptional regulation.

Drosophila insulator proteins exhibit in vivo liquid-liquid phase separation properties

Drosophila insulator proteins and the cohesin subunit Rad21 coalesce in vivo to form liquid-droplet condensates, suggesting that liquid–liquid phase separation mediates their function in 3D genome organization.

Mounting evidence implicates liquid–liquid phase separation (LLPS), the condensation of biomolecules into liquid-like droplets in the formation and dissolution of membraneless intracellular organelles (MLOs). Cells use MLOs or condensates for various biological processes, including emergency signaling and spatiotemporal control over steady-state biochemical reactions and heterochromatin formation. Insulator proteins are architectural elements involved in establishing independent domains of transcriptional activity within eukaryotic genomes. In Drosophila, insulator proteins form nuclear foci known as insulator bodies in response to osmotic stress. However, the mechanism through which insulator proteins assemble into bodies is yet to be investigated. Here, we identify signatures of LLPS by insulator bodies, including high disorder tendency in insulator proteins, scaffold–client–dependent assembly, extensive fusion behavior, sphericity, and sensitivity to 1,6-hexanediol. We also show that the cohesin subunit Rad21 is a component of insulator bodies, adding to the known insulator protein constituents and γH2Av. Our data suggest a concerted role of cohesin and insulator proteins in insulator body formation and under physiological conditions. We propose a mechanism whereby these architectural proteins modulate 3D genome organization through LLPS.

Mode and Tempo of 3D Genome Evolution in Drosophila

Topologically associating domains (TADs) are thought to play an important role in preventing gene misexpression by spatially constraining enhancer-promoter contacts. The deleterious nature of gene misexpression implies that TADs should, therefore, be conserved among related species. Several early studies comparing chromosome conformation between species reported high levels of TAD conservation; however, more recent studies have questioned these results. Furthermore, recent work suggests that TAD reorganization is not associated with extensive changes in gene expression. Here, we investigate the evolutionary conservation of TADs among 11 species of Drosophila. We use Hi-C data to identify TADs in each species and employ a comparative phylogenetic approach to derive empirical estimates of the rate of TAD evolution. Surprisingly, we find that TADs evolve rapidly. However, we also find that the rate of evolution depends on the chromatin state of the TAD, with TADs enriched for developmentally regulated chromatin evolving significantly slower than TADs enriched for broadly expressed, active chromatin. We also find that, after controlling for differences in chromatin state, highly conserved TADs do not exhibit higher levels of gene expression constraint. These results suggest that, in general, most TADs evolve rapidly and their divergence is not associated with widespread changes in gene expression. However, higher levels of evolutionary conservation and gene expression constraints in TADs enriched for developmentally regulated chromatin suggest that these TAD subtypes may be more important for regulating gene expression, likely due to the larger number of long-distance enhancer-promoter contacts associated with developmental genes.

Evolution of the ancestral mammalian karyotype and syntenic regions

Decrypting the rearrangements that drive mammalian chromosome evolution is critical to understanding the molecular bases of speciation, adaptation, and disease susceptibility. Using 8 scaffolded and 26 chromosome-scale genome assemblies representing 23/26 mammal orders, we computationally reconstructed ancestral karyotypes and syntenic relationships at 16 nodes along the mammalian phylogeny. Three different reference genomes (human, sloth, and cattle) representing phylogenetically distinct mammalian superorders were used to assess reference bias in the reconstructed ancestral karyotypes and to expand the number of clades with reconstructed genomes. The mammalian ancestor likely had 19 pairs of autosomes, with nine of the smallest chromosomes shared with the common ancestor of all amniotes (three still conserved in extant mammals), demonstrating a striking conservation of synteny for ∼320 My of vertebrate evolution. The numbers and types of chromosome rearrangements were classified for transitions between the ancestral mammalian karyotype, descendent ancestors, and extant species. For example, 94 inversions, 16 fissions, and 14 fusions that occurred over 53 My differentiated the therian from the descendent eutherian ancestor. The highest breakpoint rate was observed between the mammalian and therian ancestors (3.9 breakpoints/My). Reconstructed mammalian ancestor chromosomes were found to have distinct evolutionary histories reflected in their rates and types of rearrangements. The distributions of genes, repetitive elements, topologically associating domains, and actively transcribed regions in multispecies homologous synteny blocks and evolutionary breakpoint regions indicate that purifying selection acted over millions of years of vertebrate evolution to maintain syntenic relationships of developmentally important genes and regulatory landscapes of gene-dense chromosomes.

Distinct roles for CKM-Mediator in controlling Polycomb-dependent chromosomal interactions and priming genes for induction

Precise control of gene expression underpins normal development. This relies on mechanisms that enable communication between gene promoters and other regulatory elements. In embryonic stem cells (ESCs), the cyclin-dependent kinase module Mediator complex (CKM-Mediator) has been reported to physically link gene regulatory elements to enable gene expression and also prime genes for induction during differentiation. Here, we show that CKM-Mediator contributes little to three-dimensional genome organization in ESCs, but it has a specific and essential role in controlling interactions between inactive gene regulatory elements bound by Polycomb repressive complexes (PRCs). These interactions are established by the canonical PRC1 (cPRC1) complex but rely on CKM-Mediator, which facilitates binding of cPRC1 to its target sites. Importantly, through separation-of-function experiments, we reveal that this collaboration between CKM-Mediator and cPRC1 in creating long-range interactions does not function to prime genes for induction during differentiation. Instead, we discover that priming relies on an interaction-independent mechanism whereby the CKM supports core Mediator engagement with gene promoters during differentiation to enable gene activation.

Three-dimensional genome organization in immune cell fate and function

Immune cell development and activation demand the precise and coordinated control of transcriptional programmes. Three-dimensional (3D) organization of the genome has emerged as an important regulator of chromatin state, transcriptional activity and cell identity by facilitating or impeding long-range genomic interactions among regulatory elements and genes. Chromatin folding thus enables cell type-specific and stimulus-specific transcriptional responses to extracellular signals, which are essential for the control of immune cell fate, for inflammatory responses and for generating a diverse repertoire of antigen receptor specificities. Here, we review recent findings connecting 3D genome organization to the control of immune cell differentiation and function, and discuss how alterations in genome folding may lead to immune dysfunction and malignancy.

Causality in transcription and genome folding: Insights from X inactivation

The spatial organization of genomes is becoming increasingly understood. In mammals, where it is most investigated, this organization ties in with transcription, so an important research objective is to understand whether gene activity is a cause or a consequence of genome folding in space. In this regard, the phenomena of X-chromosome inactivation and reactivation open a unique window of investigation because of the singularities of the inactive X chromosome. Here we focus on the cause-consequence nexus between genome conformation and transcription and explain how recent results about the structural changes associated with inactivation and reactivation of the X chromosome shed light on this problem.

Non-coding genetic variation in regulatory elements determines thrombosis and hemostasis phenotypes

Since the early inception of genome-wide association studies (GWAS), it became clear that, in all diseases or traits studied, most genetic variants are likely to exert their effect on gene expression mainly by altering the function of regulatory elements. At the same time, the regulation of the gene expression field broadened its boundaries, from the univocal relationship between regulatory elements and genes to include genome organization, long-range DNA interactions, and epigenetics. Next-generation sequencing has introduced genome-wide approaches that have greatly improved our understanding of the general principles of gene expression. However, elucidating how these apply in every single genomic locus still requires painstaking experimental work, in which several independent lines of evidence are required, and often this is helped by rare genetic variants in individuals with rare diseases. This review will focus on the non-coding features of the genome involved in transcriptional regulation, that when altered, leads to known cases of inherited (familial) thrombotic and hemostatic phenotypes, emphasizing the role of enhancers and super-enhancers.

Smc3 acetylation, Pds5 and Scc2 control the translocase activity that establishes cohesin-dependent chromatin loops

Cohesin is a DNA translocase that is instrumental in the folding of the genome into chromatin loops, with functional consequences on DNA-related processes. Chromatin loop length and organization likely depend on cohesin processivity, translocation rate and stability on DNA. Here, we investigate and provide a comprehensive overview of the roles of various cohesin regulators in tuning chromatin loop expansion in budding yeast Saccharomyces cerevisiae. We demonstrate that Scc2, which stimulates cohesin ATPase activity, is also essential for cohesin translocation, driving loop expansion in vivo. Smc3 acetylation during the S phase counteracts this activity through the stabilization of Pds5, which finely tunes the size and stability of loops in G2.

Enhancer-gene specificity in development and disease

Enhancers control the establishment of spatiotemporal gene expression patterns throughout development. Over the past decade, the development of new technologies has improved our capacity to link enhancers with their target genes based on their colocalization within the same topological domains. However, the mechanisms that regulate how enhancers specifically activate some genes but not others within a given domain remain unclear. In this Review, we discuss recent insights into the factors controlling enhancer specificity, including the genetic composition of enhancers and promoters, the linear and 3D distance between enhancers and their target genes, and cell-type specific chromatin landscapes. We also discuss how elucidating the molecular principles of enhancer specificity might help us to better understand and predict the pathological consequences of human genetic, epigenetic and structural variants.

Large-Scale Chromatin Rearrangements in Cancer

Simple Summary

Cancers have many genetic mutations such as nucleotide changes, deletions, amplifications, and chromosome gains or losses. Some of these genetic alterations directly contribute to the initiation and progression of tumors. In parallel to these genetic changes, cancer cells acquire modifications to their chromatin landscape, i.e., to the marks that are carried by DNA and the histone proteins it is associated with. These “epimutations” have consequences for gene expression and genome stability, and also contribute to tumoral initiation and progression. Some of these chromatin changes are very local, affecting just one or a few genes. In contrast, some chromatin alterations observed in cancer are more widespread and affect a large part of the genome. In this review, we present different types of large-scale chromatin rearrangements in cancer, explain how they may occur, and why they are relevant for cancer diagnosis and treatment.

Abstract

Epigenetic abnormalities are extremely widespread in cancer. Some of them are mere consequences of transformation, but some actively contribute to cancer initiation and progression; they provide powerful new biological markers, as well as new targets for therapies. In this review, we examine the recent literature and focus on one particular aspect of epigenome deregulation: large-scale chromatin changes, causing global changes of DNA methylation or histone modifications. After a brief overview of the one-dimension (1D) and three-dimension (3D) epigenome in healthy cells and of its homeostasis mechanisms, we use selected examples to describe how many different events (mutations, changes in metabolism, and infections) can cause profound changes to the epigenome and fuel cancer. We then present the consequences for therapies and briefly discuss the role of single-cell approaches for the future progress of the field.

Nuclear receptor activation shapes spatial genome organization essential for gene expression control: lessons learned from the vitamin D receptor

Spatial genome organization is tightly controlled by several regulatory mechanisms and is essential for gene expression control. Nuclear receptors are ligand-activated transcription factors that modulate physiological and pathophysiological processes and are primary pharmacological targets. DNA binding of the important loop-forming insulator protein CCCTC-binding factor (CTCF) was modulated by 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3). We performed CTCF HiChIP assays to produce the first genome-wide dataset of CTCF long-range interactions in 1,25(OH)2D3-treated cells, and to determine whether dynamic changes of spatial chromatin interactions are essential for fine-tuning of nuclear receptor signaling. We detected changes in 3D chromatin organization upon vitamin D receptor (VDR) activation at 3.1% of all observed CTCF interactions. VDR binding was enriched at both differential loop anchors and within differential loops. Differential loops were observed in several putative functional roles including TAD border formation, promoter-enhancer looping, and establishment of VDR-responsive insulated neighborhoods. Vitamin D target genes were enriched in differential loops and at their anchors. Secondary vitamin D effects related to dynamic chromatin domain changes were linked to location of downstream transcription factors in differential loops. CRISPR interference and loop anchor deletion experiments confirmed the functional relevance of nuclear receptor ligand-induced adjustments of the chromatin 3D structure for gene expression regulation.

Conformational Changes in Three-Dimensional Chromatin Structure in Paulownia fortunei After Phytoplasma Infection

Higher-order chromatin structures play important roles in regulating multiple biological processes such as growth and development as well as biotic and abiotic stress response. However, little is known about three-dimensional chromatin structures in Paulownia or about whole-genome chromatin conformational changes that occur in response to Paulownia witches' broom (PaWB) disease. We used high-throughput chromosome conformation capture (Hi-C) to obtain genome-wide profiles of chromatin conformation in both healthy and phytoplasma-infected Paulownia fortunei genome. The heat map results indicated that the strongest interactions between chromosomes were in the telomeres. We confirmed that the main structural characteristics of A/B compartments, topologically associated domains, and chromatin loops were prominent in the Paulownia genome and were clearly altered in phytoplasma-infected plants. By combining chromatin immunoprecipitation sequencing, Hi-C signals, and RNA sequencing data, we inferred that the chromatin structure changed and the modification levels of three histones (H3K4me3/K9ac/K36me3) increased in phytoplasma-infected P. fortunei, which was associated with changes of transcriptional activity. We concluded that for epigenetic modifications, transcriptional activity might function in combination to shape chromatin packing in healthy and phytoplasm-infected Paulownia. Finally, 11 genes (e.g., RPN6, Sec61 subunit-α) that were commonly located at specific topologically associated domain boundaries, A/B compartment switching and specific loops, and had been associated with histone marks were identified and considered as closely related to PaWB stress. Our results provide new insights into the nexus between gene regulation and chromatin conformational alterations in nonmodel plants upon phytopathogen infection and plant disease resistance.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

A CTCF-Binding Element and Histone Deacetylation Cooperatively Maintain Chromatin Loops, Linking to Long-Range Gene Regulation in Cancer Genomes

Background

Genes spanning long chromosomal domains are coordinately regulated in human genome, which contribute to global gene dysregulation and carcinogenesis in cancer. It has been noticed that epigenetic modification and chromatin architecture may participate in the regulation process. However, the regulation patterns and functional elements of long-range gene regulation are unclear.

Methods

Based on the clinical transcriptome data from different tumor sets, a novel expressional correlation analysis pipeline was performed to classify the co-regulated regions and subsets of intercorrelated regions. The GLAM2 program was used to predict conserved DNA elements that enriched in regions. Two conserved elements were selected to delete in Ishikawa and HeLa cells by CRISPR-Cas9. SAHA treatment and HDAC knockdown were used to change the histone acetylation status. Using qPCR, MTT, and scratch healing assay, we evaluate the effect on gene expression and cancer cell phenotype. By DNA pull-down and ChIP, the element-binding proteins were testified. 3C and 3D-FISH were performed to depict the alteration in chromatin architecture.

Results

In multiple cancer genomes, we classified subsets of coordinately regulated regions (sub-CRRs) that possibly shared the same regulatory mechanisms and exhibited similar expression patterns. A new conserved DNA element (CRE30) was enriched in sub-CRRs and associated with cancer patient survival. CRE30 could restrict gene regulation in sub-CRRs and affect cancer cell phenotypes. DNA pull-down showed that multiple proteins including CTCF were recruited on the CRE30 locus, and ChIP assay confirmed the CTCF-binding signals. Subsequent results uncovered that as an essential element, CRE30 maintained chromatin loops and mediated a compact chromatin architecture. Moreover, we found that blocking global histone deacetylation induced chromatin loop disruption and CTCF dropping in the region containing CRE30, linked to promoted gene regulation. Additionally, similar effects were observed with CRE30 deletion in another locus of chromosome 8.

Conclusions

Our research clarified a new functional element that recruits CTCF and collaborates with histone deacetylation to maintain high-order chromatin organizations, linking to long-range gene regulation in cancer genomes. The findings highlight a close relationship among conserved DNA element, epigenetic modification, and chromatin architecture in long-range gene regulation process.

Contribution of 3D genome topological domains to genetic risk of cancers: a genome-wide computational study

Background

Genome-wide association studies have identified statistical associations between various diseases, including cancers, and a large number of single-nucleotide polymorphisms (SNPs). However, they provide no direct explanation of the mechanisms underlying the association. Based on the recent discovery that changes in three-dimensional genome organization may have functional consequences on gene regulation favoring diseases, we investigated systematically the genome-wide distribution of disease-associated SNPs with respect to a specific feature of 3D genome organization: topologically associating domains (TADs) and their borders.

Results

For each of 449 diseases, we tested whether the associated SNPs are present in TAD borders more often than observed by chance, where chance (i.e., the null model in statistical terms) corresponds to the same number of pointwise loci drawn at random either in the entire genome, or in the entire set of disease-associated SNPs listed in the GWAS catalog. Our analysis shows that a fraction of diseases displays such a preferential localization of their risk loci. Moreover, cancers are relatively more frequent among these diseases, and this predominance is generally enhanced when considering only intergenic SNPs. The structure of SNP-based diseasome networks confirms that localization of risk loci in TAD borders differs between cancers and non-cancer diseases. Furthermore, different TAD border enrichments are observed in embryonic stem cells and differentiated cells, consistent with changes in topological domains along embryogenesis and delineating their contribution to disease risk.

Conclusions

Our results suggest that, for certain diseases, part of the genetic risk lies in a local genetic variation affecting the genome partitioning in topologically insulated domains. Investigating this possible contribution to genetic risk is particularly relevant in cancers. This study thus opens a way of interpreting genome-wide association studies, by distinguishing two types of disease-associated SNPs: one with an effect on an individual gene, the other acting in interplay with 3D genome organization.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-022-00375-2.

Simultaneous Quantification of Spatial Genome Positioning and Transcriptomics in Single Cells with scDam&T-Seq

Spatial genome organization is considered to play an important role in mammalian cells, by guiding gene expression programs and supporting lineage specification. Yet it is still an outstanding question in the field what the direct impact of spatial genome organization on gene expression is. To elucidate this relationship further, we have recently developed scDam&T-seq, a method that simultaneously quantifies protein-DNA interactions and transcriptomes in single cells. This method efficiently combines two preexisting methods: DamID for measuring protein-DNA contacts and CEL-Seq2 for quantification of the transcriptome in single cells. scDam&T-seq has been successfully applied to measure DNA contacts with the nuclear lamina, while at the same time revealing the effect of these contacts on gene expression. This method is applicable to many different proteins of interest and can thereby aid in studying the relationship between protein-DNA interactions and gene expression in single cells.

A novel complex genomic rearrangement affecting the KCNJ2 regulatory region causes a variant of Cooks syndrome

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

EPIsHilbert: Prediction of Enhancer-Promoter Interactions via Hilbert Curve Encoding and Transfer Learning

Enhancer-promoter interactions (EPIs) play a significant role in the regulation of gene transcription. However, enhancers may not necessarily interact with the closest promoters, but with distant promoters via chromatin looping. Considering the spatial position relationship between enhancers and their target promoters is important for predicting EPIs. Most existing methods only consider sequence information regardless of spatial information. On the other hand, recent computational methods lack generalization capability across different cell line datasets. In this paper, we propose EPIsHilbert, which uses Hilbert curve encoding and two transfer learning approaches. Hilbert curve encoding can preserve the spatial position information between enhancers and promoters. Additionally, we use visualization techniques to explore important sequence fragments that have a high impact on EPIs and the spatial relationships between them. Transfer learning can improve prediction performance across cell lines. In order to further prove the effectiveness of transfer learning, we analyze the sequence coincidence of different cell lines. Experimental results demonstrate that EPIsHilbert is a state-of-the-art model that is superior to most of the existing methods both in specific cell lines and cross cell lines.

DeTOKI identifies and characterizes the dynamics of chromatin TAD-like domains in a single cell

Topologically associating domains (TAD) are a key structure of the 3D mammalian genomes. However, the prevalence and dynamics of TAD-like domains in single cells remain elusive. Here we develop a new algorithm, named deTOKI, to decode TAD-like domains with single-cell Hi-C data. By non-negative matrix factorization, deTOKI seeks regions that insulate the genome into blocks with minimal chance of clustering. deTOKI outperforms competing tools and reliably identifies TAD-like domains in single cells. Finally, we find that TAD-like domains are not only prevalent, but also subject to tight regulation in single cells.

How subtle changes in 3D structure can create large changes in transcription

Animal genomes are organized into topologically associated domains (TADs). TADs are thought to contribute to gene regulation by facilitating enhancer-promoter (E-P) contacts within a TAD and preventing these contacts across TAD borders. However, the absolute difference in contact frequency across TAD boundaries is usually less than 2-fold, even though disruptions of TAD borders can change gene expression by 10-fold. Existing models fail to explain this hypersensitive response. Here, we propose a futile cycle model of enhancer-mediated regulation that can exhibit hypersensitivity through bistability and hysteresis. Consistent with recent experiments, this regulation does not exhibit strong correlation between E-P contact and promoter activity, even though regulation occurs through contact. Through mathematical analysis and stochastic simulation, we show that this system can create an illusion of E-P biochemical specificity and explain the importance of weak TAD boundaries. It also offers a mechanism to reconcile apparently contradictory results from recent global TAD disruption with local TAD boundary deletion experiments. Together, these analyses advance our understanding of cis-regulatory contacts in controlling gene expression and suggest new experimental directions.

A description of the relationship in healthy longevity and aging-related disease: from gene to protein

Human longevity is a complex phenotype influenced by both genetic and environmental factors. It is also known to be associated with various types of age-related diseases, such as Alzheimer's disease (AD) and cardiovascular disease (CVD). The central dogma of molecular biology demonstrates the conversion of DNA to RNA to the encoded protein. These proteins interact to form complex cell signaling pathways, which perform various biological functions. With prolonged exposure to the environment, the in vivo homeostasis adapts to the changes, and finally, humans adopt the phenotype of longevity or aging-related diseases. In this review, we focus on two different states: longevity and aging-related diseases, including CVD and AD, to discuss the relationship between genetic characteristics, including gene variation, the level of gene expression, regulation of gene expression, the level of protein expression, both genetic and environmental influences and homeostasis based on these phenotypes shown in organisms.

What can go wrong in the non-coding genome and how to interpret whole genome sequencing data

Whole exome sequencing discovers causative mutations in less than 50 % of rare disease patients, suggesting the presence of additional mutations in the non-coding genome. So far, non-coding mutations have been identified in less than 0.2 % of individuals with genetic diseases listed in the ClinVar database and exhibit highly diverse molecular mechanisms. In contrast to our capability to sequence the whole genome, our ability to discover and functionally confirm such non-coding mutations is lagging behind severely. We discuss the problems and present examples of confirmed mutations in deep intronic sequences, non-coding triplet repeats, enhancers, and larger structural variants and highlight their proposed disease mechanisms. Finally, we discuss the type of data that would be required to establish non-coding mutation detection in routine diagnostics.

Altered chromatin conformation and transcriptional regulation in watermelon following genome doubling

Polyploidy has played a crucial role in plant evolution, development and function. Synthetic autopolyploid represents an ideal system to investigate the effects of polyploidization on transcriptional regulation. In this study, we deciphered the impact of genome duplication at phenotypic and molecular levels in watermelon. Overall, 88% of the genes in tetraploid watermelon followed a >1:1 dosage effect, and accordingly, differentially expressed genes were largely upregulated. In addition, a great number of hypomethylated regions (1688) were identified in an isogenic tetraploid watermelon. These differentially methylated regions were localized in promoters and intergenic regions and near transcriptional start sites of the identified upregulated genes, which enhances the importance of methylation in gene regulation. These changes were reflected in sophisticated higher-order chromatin structures. The genome doubling caused switching of 108 A and 626 B compartments that harbored genes associated with growth, development and stress responses.

Independence of chromatin conformation and gene regulation during Drosophila dorsoventral patterning

The relationship between chromatin organization and gene regulation remains unclear. While disruption of chromatin domains and domain boundaries can lead to misexpression of developmental genes, acute depletion of regulators of genome organization has a relatively small effect on gene expression. It is therefore uncertain whether gene expression and chromatin state drive chromatin organization or whether changes in chromatin organization facilitate cell-type-specific activation of gene expression. Here, using the dorsoventral patterning of the Drosophila melanogaster embryo as a model system, we provide evidence for the independence of chromatin organization and dorsoventral gene expression. We define tissue-specific enhancers and link them to expression patterns using single-cell RNA-seq. Surprisingly, despite tissue-specific chromatin states and gene expression, chromatin organization is largely maintained across tissues. Our results indicate that tissue-specific chromatin conformation is not necessary for tissue-specific gene expression but rather acts as a scaffold facilitating gene expression when enhancers become active.

Evolution of Genome-Organizing Long Non-coding RNAs in Metazoans

Long non-coding RNAs (lncRNAs) have important regulatory functions across eukarya. It is now clear that many of these functions are related to gene expression regulation through their capacity to recruit epigenetic modifiers and establish chromatin interactions. Several lncRNAs have been recently shown to participate in modulating chromatin within the spatial organization of the genome in the three-dimensional space of the nucleus. The identification of lncRNA candidates is challenging, as it is their functional characterization. Conservation signatures of lncRNAs are different from those of protein-coding genes, making identifying lncRNAs under selection a difficult task, and the homology between lncRNAs may not be readily apparent. Here, we review the evidence for these higher-order genome organization functions of lncRNAs in animals and the evolutionary signatures they display.

Genome architecture and expression 2019-2020: the transition phase

To set the stage for this collection, let us reflect on how the field of genome architecture and expression is undergoing current movement including, inter alia, the adoption of physico-chemical approaches of which the article's subtitle is a pun. This field had gained powerful new momentum 5-7 years ago and now one senses certain tipping points, as the foundational elements of the field are being refined and expanded.