Complex multi-enhancer contacts captured by Genome Architecture Mapping (GAM)

Biochemical properties of chromatin domains define genome compartmentalization

Chromatin three-dimensional (3D) organization inside the cell nucleus determines the separation of euchromatin and heterochromatin domains. Their segregation results in the definition of active and inactive chromatin compartments, whereby the local concentration of associated proteins, RNA and DNA results in the formation of distinct subnuclear structures. Thus, chromatin domains spatially confined in a specific 3D nuclear compartment are expected to share similar epigenetic features and biochemical properties, in terms of accessibility and solubility. Based on this rationale, we developed the 4f-SAMMY-seq to map euchromatin and heterochromatin based on their accessibility and solubility, starting from as little as 10 000 cells. Adopting a tailored bioinformatic data analysis approach we reconstruct also their 3D segregation in active and inactive chromatin compartments and sub-compartments, thus recapitulating the characteristic properties of distinct chromatin states. A key novelty of the new method is the capability to map both the linear segmentation of open and closed chromatin domains, as well as their compartmentalization in one single experiment.

Enhancing Hi-C contact matrices for loop detection with Capricorn: a multiview diffusion model

MOTIVATION

High-resolution Hi-C contact matrices reveal the detailed three-dimensional architecture of the genome, but high-coverage experimental Hi-C data are expensive to generate. Simultaneously, chromatin structure analyses struggle with extremely sparse contact matrices. To address this problem, computational methods to enhance low-coverage contact matrices have been developed, but existing methods are largely based on resolution enhancement methods for natural images and hence often employ models that do not distinguish between biologically meaningful contacts, such as loops and other stochastic contacts.

RESULTS

We present Capricorn, a machine learning model for Hi-C resolution enhancement that incorporates small-scale chromatin features as additional views of the input Hi-C contact matrix and leverages a diffusion probability model backbone to generate a high-coverage matrix. We show that Capricorn outperforms the state of the art in a cross-cell-line setting, improving on existing methods by 17% in mean squared error and 26% in F1 score for chromatin loop identification from the generated high-coverage data. We also demonstrate that Capricorn performs well in the cross-chromosome setting and cross-chromosome, cross-cell-line setting, improving the downstream loop F1 score by 14% relative to existing methods. We further show that our multiview idea can also be used to improve several existing methods, HiCARN and HiCNN, indicating the wide applicability of this approach. Finally, we use DNA sequence to validate discovered loops and find that the fraction of CTCF-supported loops from Capricorn is similar to those identified from the high-coverage data. Capricorn is a powerful Hi-C resolution enhancement method that enables scientists to find chromatin features that cannot be identified in the low-coverage contact matrix.

AVAILABILITY AND IMPLEMENTATION

Implementation of Capricorn and source code for reproducing all figures in this paper are available at https://github.com/CHNFTQ/Capricorn.

Regulating epithelial-mesenchymal plasticity from 3D genome organization

Epithelial-mesenchymal transition (EMT) is a dynamic process enabling polarized epithelial cells to acquire mesenchymal features implicated in development and carcinoma progression. As our understanding evolves, it is clear the reversible execution of EMT arises from complex epigenomic regulation involving histone modifications and 3-dimensional (3D) genome structural changes, leading to a cascade of transcriptional events. This review summarizes current knowledge on chromatin organization in EMT, with a focus on hierarchical structures of the 3D genome and chromatin accessibility changes.

FISHnet: Detecting chromatin domains in single-cell sequential Oligopaints imaging data

Sequential Oligopaints DNA FISH is an imaging technique that measures higher-order genome folding at single-allele resolution via multiplexed, probe-based tracing. Currently there is a paucity of algorithms to identify 3D genome features in sequential Oligopaints data. Here, we present FISHnet, a graph theory method based on optimization of network modularity to detect chromatin domains and boundaries in pairwise distance matrices. FISHnet uncovers cell type-specific domain-like folding patterns on single alleles, thus enabling future studies aiming to elucidate the role for single-cell folding variation on genome function.

Single-cell mapping of cell-type specific chromatin architecture in the central nervous system

Determining how chromatin is structured in the nucleus is critical to studying its role in gene regulation. Recent advances in the analysis of single-cell chromatin architecture have considerably improved our understanding of cell-type-specific chromosome conformation and nuclear architecture. In this review, we discuss the methods used for analysis of 3D chromatin conformation, including sequencing-based methods, imaging-based techniques, and computational approaches. We further review the application of these methods in the study of the role of chromatin topology in neural development and disorders.

Chromatin organization in myelodysplastic syndrome

Disordered chromatin organization has emerged as a new aspect of the pathogenesis of myelodysplastic syndrome (MDS). Characterized by lineage dysplasia and a high transformation rate to acute myeloid leukemia (AML), the genetic determinant of MDS is thought to be the main driver of the disease's progression. Among the recurrently mutated pathways, alterations in chromatin organization, such as the cohesin complex, have a profound impact on hematopoietic stem cell (HSC) function and lineage commitment. The cohesin complex is a ring-like structure comprised of structural maintenance of chromosomes (SMC), RAD21, and STAG proteins that involve three-dimensional (3D) genome organization via loop extrusion in mammalian cells. The partial loss of the functional cohesin ring leads to altered chromatin accessibility specific to key hematopoietic transcription factors, which is thought to be the molecular mechanism of cohesin dysfunction. Currently, there are no specific targeting agents for cohesin mutant MDS/AML. Potential therapeutic strategies have been proposed based on the current understanding of cohesin mutant leukemogenesis. Here, we will review the recent advances in investigation and targeting approaches against cohesin mutant MDS/AML.

Attraction and disruption: how loop extrusion and compartmentalisation shape the nuclear genome

Chromatin loops, which bring two distal loci of the same chromosome into close physical proximity, are the ubiquitous units of the three-dimensional genome. Recent advances in understanding the spatial organisation of chromatin suggest that several distinct mechanisms control chromatin interactions, such as loop extrusion by cohesin complexes, compartmentalisation by phase separation, direct protein-protein interactions and others. Here, we review different types of chromatin loops and highlight the factors and processes involved in their regulation. We discuss how loop extrusion and compartmentalisation shape chromatin interactions and how these two processes can either positively or negatively influence each other.

Identification of drug responsive enhancers by predicting chromatin accessibility change from perturbed gene expression profiles

Individual may response to drug treatment differently due to their genetic variants located in enhancers. These variants can alter transcription factor's (TF) binding strength, affect enhancer's chromatin activity or interaction, and eventually change expression level of downstream gene. Here, we propose a computational framework, PERD, to Predict the Enhancers Responsive to Drug. A machine learning model was trained to predict the genome-wide chromatin accessibility from transcriptome data using the paired expression and chromatin accessibility data collected from ENCODE and ROADMAP. Then the model was applied to the perturbed gene expression data from Connectivity Map (CMAP) and Cancer Drug-induced gene expression Signature DataBase (CDS-DB) and identify drug responsive enhancers with significantly altered chromatin accessibility. Furthermore, the drug responsive enhancers were related to the pharmacogenomics genome-wide association studies (PGx GWAS). Stepping on the traditional drug-associated gene signatures, PERD holds the promise to enhance the causality of drug perturbation by providing candidate regulatory element of those drug associated genes.

Statistical modeling and significance estimation of multi-way chromatin contacts with HyperloopFinder

Recent advances in chromatin conformation capture technologies, such as SPRITE and Pore-C, have enabled the detection of simultaneous contacts among multiple chromatin loci. This has made it possible to investigate the cooperative transcriptional regulation involving multiple genes and regulatory elements at the resolution of a single molecule. However, these technologies are unavoidably subject to the random polymer looping effect and technical biases, making it challenging to distinguish genuine regulatory relationships directly from random polymer interactions. Here, we present HyperloopFinder, a method for identifying regulatory multi-way chromatin contacts (hyperloops) by jointly modeling the random polymer looping effect and technical biases to estimate the statistical significance of multi-way contacts. The results show that our model can accurately estimate the expected interaction frequency of multi-way contacts based on the distance distribution of pairwise contacts, revealing that most multi-way contacts can be formed by randomly linking the pairwise contacts adjacent to each other. Moreover, we observed the spatial colocalization of the interaction sites of hyperloops from image-based data. Our results also revealed that hyperloops can function as scaffolds for the cooperation among multiple genes and regulatory elements. In summary, our work contributes novel insights into higher-order chromatin structures and functions and has the potential to enhance our understanding of transcriptional regulation and other cellular processes.

A single dose of cocaine rewires the 3D genome structure of midbrain dopamine neurons

Midbrain dopamine neurons (DNs) respond to a first exposure to addictive drugs and play key roles in chronic drug usage1-3. As the synaptic and transcriptional changes that follow an acute cocaine exposure are mostly resolved within a few days4,5, the molecular changes that encode the long-term cellular memory of the exposure within DNs remain unknown. To investigate whether a single cocaine exposure induces long-term changes in the 3D genome structure of DNs, we applied Genome Architecture Mapping and single nucleus transcriptomic analyses in the mouse midbrain. We found extensive rewiring of 3D genome architecture at 24 hours past exposure which remains or worsens by 14 days, outlasting transcriptional responses. The cocaine-induced chromatin rewiring occurs at all genomic scales and affects genes with major roles in cocaine-induced synaptic changes. A single cocaine exposure triggers extensive long-lasting changes in chromatin condensation in post-synaptic and post-transcriptional regulatory genes, for example the unfolding of Rbfox1 which becomes most prominent 14 days post exposure. Finally, structurally remodeled genes are most expressed in a specific DN sub-type characterized by low expression of the dopamine auto-receptor Drd2, a key feature of highly cocaine-sensitive cells. These results reveal an important role for long-lasting 3D genome remodelling in the cellular memory of a single cocaine exposure, providing new hypotheses for understanding the inception of drug addiction and 3D genome plasticity.

Extensive folding variability between homologous chromosomes in mammalian cells

Genetic variation and 3D chromatin structure have major roles in gene regulation. Due to challenges in mapping chromatin conformation with haplotype-specific resolution, the effects of genetic sequence variation on 3D genome structure and gene expression imbalance remain understudied. Here, we applied Genome Architecture Mapping (GAM) to a hybrid mouse embryonic stem cell (mESC) line with high density of single nucleotide polymorphisms (SNPs). GAM resolved haplotype-specific 3D genome structures with high sensitivity, revealing extensive allelic differences in chromatin compartments, topologically associating domains (TADs), long-range enhancer-promoter contacts, and CTCF loops. Architectural differences often coincide with allele-specific differences in gene expression, mediated by Polycomb repression. We show that histone genes are expressed with allelic imbalance in mESCs, are involved in haplotype-specific chromatin contact marked by H3K27me3, and are targets of Polycomb repression through conditional knockouts of Ezh2 or Ring1b. Our work reveals highly distinct 3D folding structures between homologous chromosomes, and highlights their intricate connections with allelic gene expression.

Super-enhancer interactomes from single cells link clustering and transcription

Regulation of gene expression hinges on the interplay between enhancers and promoters, traditionally explored through pairwise analyses. Recent advancements in mapping genome folding, like GAM, SPRITE, and multi-contact Hi-C, have uncovered multi-way interactions among super-enhancers (SEs), spanning megabases, yet have not measured their frequency in single cells or the relationship between clustering and transcription. To close this gap, here we used multiplexed imaging to map the 3D positions of 376 SEs across thousands of mammalian nuclei. Notably, our single-cell images reveal that while SE-SE contacts are rare, SEs often form looser associations we termed "communities". These communities, averaging 4-5 SEs, assemble cooperatively under the combined effects of genomic tethers, Pol2 clustering, and nuclear compartmentalization. Larger communities are associated with more frequent and larger transcriptional bursts. Our work provides insights about the SE interactome in single cells that challenge existing hypotheses on SE clustering in the context of transcriptional regulation.

Potential roles of inter-chromosomal interactions in cell fate determination

Mammalian genomic DNA is packed in a small nucleus, and its folding and organization in the nucleus are critical for gene regulation and cell fate determination. In interphase, chromosomes are compartmentalized into certain nuclear spaces and territories that are considered incompatible with each other. The regulation of gene expression is influenced by the epigenetic characteristics of topologically associated domains and A/B compartments within chromosomes (intrachromosomal). Previously, interactions among chromosomes detected via chromosome conformation capture-based methods were considered noise or artificial errors. However, recent studies based on newly developed ligation-independent methods have shown that inter-chromosomal interactions play important roles in gene regulation. This review summarizes the recent understanding of spatial genomic organization in mammalian interphase nuclei and discusses the potential mechanisms that determine cell identity. In addition, this review highlights the potential role of inter-chromosomal interactions in early mouse development.

SRRM2 splicing factor modulates cell fate in early development

Embryo development is an orchestrated process that relies on tight regulation of gene expression to guide cell differentiation and fate decisions. The Srrm2 splicing factor has recently been implicated in developmental disorders and diseases, but its role in early mammalian development remains unexplored. Here, we show that Srrm2 dosage is critical for maintaining embryonic stem cell pluripotency and cell identity. Srrm2 heterozygosity promotes loss of stemness, characterised by the coexistence of cells expressing naive and formative pluripotency markers, together with extensive changes in gene expression, including genes regulated by serum-response transcription factor (SRF) and differentiation-related genes. Depletion of Srrm2 by RNA interference in embryonic stem cells shows that the earliest effects of Srrm2 heterozygosity are specific alternative splicing events on a small number of genes, followed by expression changes in metabolism and differentiation-related genes. Our findings unveil molecular and cellular roles of Srrm2 in stemness and lineage commitment, shedding light on the roles of splicing regulators in early embryogenesis, developmental diseases and tumorigenesis.

Predicting chromatin conformation contact maps

Over the past 15 years, a variety of next-generation sequencing assays have been developed for measuring the 3D conformation of DNA in the nucleus. Each of these assays gives, for a particular cell or tissue type, a distinct picture of 3D chromatin architecture. Accordingly, making sense of the relationship between genome structure and function requires teasing apart two closely related questions: how does chromatin 3D structure change from one cell type to the next, and how do different measurements of that structure differ from one another, even when the two assays are carried out in the same cell type? In this work, we assemble a collection of chromatin 3D datasets-each represented as a 2D contact map- spanning multiple assay types and cell types. We then build a machine learning model that predicts missing contact maps in this collection. We use the model to systematically explore how genome 3D architecture changes, at the level of compartments, domains, and loops, between cell type and between assay types.

The scales, mechanisms, and dynamics of the genome architecture

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

Footprints of loop extrusion in statistics of intra-chromosomal distances: An analytically solvable model

Active loop extrusion-the process of formation of dynamically growing chromatin loops due to the motor activity of DNA-binding protein complexes-is a firmly established mechanism responsible for chromatin spatial organization at different stages of a cell cycle in eukaryotes and bacteria. The theoretical insight into the effect of loop extrusion on the experimentally measured statistics of chromatin conformation can be gained with an appropriately chosen polymer model. Here, we consider the simplest analytically solvable model of an interphase chromosome, which is treated as an ideal chain with disorder of sufficiently sparse random loops whose conformations are sampled from the equilibrium ensemble. This framework allows us to arrive at the closed-form analytical expression for the mean-squared distance between pairs of genomic loci, which is valid beyond the one-loop approximation in diagrammatic representation. In addition, we analyze the loop-induced deviation of chain conformations from the Gaussian statistics by calculating kurtosis of probability density of the pairwise separation vector. The presented results suggest the possible ways of estimating the characteristics of the loop extrusion process based on the experimental data on the scale-dependent statistics of intra-chromosomal pair-wise distances.

Multifaceted roles of cohesin in regulating transcriptional loops

Cohesin is required for chromatin loop formation. However, its precise role in regulating gene transcription remains largely unknown. We investigated the relationship between cohesin and RNA Polymerase II (RNAPII) using single-molecule mapping and live-cell imaging methods in human cells. Cohesin-mediated transcriptional loops were highly correlated with those of RNAPII and followed the direction of gene transcription. Depleting RAD21, a subunit of cohesin, resulted in the loss of long-range (>100 kb) loops between distal (super-)enhancers and promoters of cell-type-specific genes. By contrast, the short-range (<50 kb) loops were insensitive to RAD21 depletion and connected genes that are mostly housekeeping. This result explains why only a small fraction of genes are affected by the loss of long-range chromatin interactions due to cohesin depletion. Remarkably, RAD21 depletion appeared to up-regulate genes located in early initiation zones (EIZ) of DNA replication, and the EIZ signals were amplified drastically without RAD21. Our results revealed new mechanistic insights of cohesin's multifaceted roles in establishing transcriptional loops, preserving long-range chromatin interactions for cell-specific genes, and maintaining timely order of DNA replication.

Oncogenic enhancers prime quiescent metastatic cells to escape NK immune surveillance by eliciting transcriptional memory

Metastasis arises from disseminated tumour cells (DTCs) that are characterized by intrinsic phenotypic plasticity and the capability of seeding to secondary organs. DTCs can remain latent for years before giving rise to symptomatic overt metastasis. In this context, DTCs fluctuate between a quiescent and proliferative state in response to systemic and microenvironmental signals including immune-mediated surveillance. Despite its relevance, how intrinsic mechanisms sustain DTCs plasticity has not been addressed. By interrogating the epigenetic state of metastatic cells, we find that tumour progression is coupled with the activation of oncogenic enhancers that are organized in variable interconnected chromatin domains. This spatial chromatin context leads to the activation of a robust transcriptional response upon repeated exposure to retinoic acid (RA). We show that this adaptive mechanism sustains the quiescence of DTCs through the activation of the master regulator SOX9. Finally, we determine that RA-stimulated transcriptional memory increases the fitness of metastatic cells by supporting the escape of quiescent DTCs from NK-mediated immune surveillance. Overall, these findings highlight the contribution of oncogenic enhancers in establishing transcriptional memories as an adaptive mechanism to reinforce cancer dormancy and immune escape, thus amenable for therapeutic intervention.

Machine and deep learning methods for predicting 3D genome organization

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

Chromosome compartmentalization: causes, changes, consequences, and conundrums

The spatial segregation of the genome into compartments is a major feature of 3D genome organization. New data on mammalian chromosome organization across different conditions reveal important information about how and why these compartments form and change. A combination of epigenetic state, nuclear body tethering, physical forces, gene expression, and replication timing (RT) can all influence the establishment and alteration of chromosome compartments. We review the causes and implications of genomic regions undergoing a 'compartment switch' that changes their physical associations and spatial location in the nucleus. About 20-30% of genomic regions change compartment during cell differentiation or cancer progression, whereas alterations in response to a stimulus within a cell type are usually much more limited. However, even a change in 1-2% of genomic bins may have biologically relevant implications. Finally, we review the effects of compartment changes on gene regulation, DNA damage repair, replication, and the physical state of the cell.

Giant pandas in captivity undergo short-term adaptation in nerve-related pathways

BACKGROUND

Behaviors in captive animals, including changes in appetite, activity level, and social interaction, are often seen as adaptive responses. However, these behaviors may become progressively maladaptive, leading to stress, anxiety, depression, and other negative reactions in animals.

RESULTS

In this study, we investigated the whole-genome sequencing data of 39 giant panda individuals, including 11 in captivity and 28 in the wild. To eliminate the mountain range effect and focus on the factor of captivity only, we first performed a principal component analysis. We then enumerated the 21,474,180 combinations of wild giant pandas (11 chosen from 28) and calculated their distances from the 11 captive individuals. The 11 wild individuals with the closest distances were used for the subsequent analysis. The linkage disequilibrium (LD) patterns demonstrated that the population was almost eliminated. We identified 505 robust selected genomic regions harboring at least one SNP, and the absolute frequency difference was greater than 0.6 between the two populations. GO analysis revealed that genes in these regions were mainly involved in nerve-related pathways. Furthermore, we identified 22 GO terms for which the selection strength significantly differed between the two populations, and there were 10 nerve-related pathways among them. Genes in the differentially abundant regions were involved in nerve-related pathways, indicating that giant pandas in captivity underwent minor genomic selection. Additionally, we investigated the relationship between genetic variation and chromatin conformation structures. We found that nucleotide diversity (θπ) in the captive population was correlated with chromatin conformation structures, which included A/B compartments, topologically associated domains (TADs) and TAD-cliques. For each GO term, we then compared the expression level of genes regulated by the above four factors (AB index, TAD intactness, TAD clique and PEI) with the corresponding genomic background. The retained 10 GO terms were all coordinately regulated by the four factors, and three of them were associated with nerve-related pathways.

CONCLUSIONS

This study revealed that giant pandas in captivity undergo short-term adaptation in nerve-related pathways. Furthermore, it provides new insights into the molecular mechanism of gene expression regulation under short-term adaptation to environmental change.

Computational methods for analysing multiscale 3D genome organization

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

Epigenetic regulatory layers in the 3D nucleus

Nearly 7 decades have elapsed since Francis Crick introduced the central dogma of molecular biology, as part of his ideas on protein synthesis, setting the fundamental rules of sequence information transfer from DNA to RNAs and proteins. We have since learned that gene expression is finely tuned in time and space, due to the activities of RNAs and proteins on regulatory DNA elements, and through cell-type-specific three-dimensional conformations of the genome. Here, we review major advances in genome biology and discuss a set of ideas on gene regulation and highlight how various biomolecular assemblies lead to the formation of structural and regulatory features within the nucleus, with roles in transcriptional control. We conclude by suggesting further developments that will help capture the complex, dynamic, and often spatially restricted events that govern gene expression in mammalian cells.

3D Enhancer-promoter networks provide predictive features for gene expression and coregulation in early embryonic lineages

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages: the trophectoderm, the epiblast and the primitive endoderm. Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements through which transcriptional regulators enact these fates remain understudied. Here, we characterize, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observe extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although distinct groups of genes are irresponsive to topological changes. In each lineage, a high degree of connectivity, or 'hubness', positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a predictive model for transcriptional regulation (3D-HiChAT) that outperforms models using only 1D promoter or proximal variables to predict levels and cell-type specificity of gene expression. Using 3D-HiChAT, we identify, in silico, candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments, we validate several enhancers that control gene expression in their respective lineages. Our study identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to comprehensively understand lineage-specific transcriptional behaviors.

Revisiting chromatin packaging in mouse sperm

Mammalian sperm show an unusual and heavily compacted genomic packaging state. In addition to its role in organizing the compact and hydrodynamic sperm head, it has been proposed that sperm chromatin architecture helps to program gene expression in the early embryo. Scores of genome-wide surveys in sperm have reported patterns of chromatin accessibility, nucleosome localization, histone modification, and chromosome folding. Here, we revisit these studies in light of recent reports that sperm obtained from the mouse epididymis are contaminated with low levels of cell-free chromatin. In the absence of proper sperm lysis, we readily recapitulate multiple prominent genome-wide surveys of sperm chromatin, suggesting that these profiles primarily reflect contaminating cell-free chromatin. Removal of cell-free DNA, and appropriate lysis conditions, are together required to reveal a sperm chromatin state distinct from most previous reports. Using ATAC-seq to explore relatively accessible genomic loci, we identify a landscape of open loci associated with early development and transcriptional control. Histone modification and chromosome folding profiles also strongly support the hypothesis that prior studies suffer from contamination, but technical challenges associated with reliably preserving the architecture of the compacted sperm head prevent us from confidently assaying true localization patterns for these epigenetic marks. Together, our studies show that our knowledge of chromosome packaging in mammalian sperm remains largely incomplete, and motivate future efforts to more accurately characterize genome organization in mature sperm.

Super-enhancers include classical enhancers and facilitators to fully activate gene expression

Super-enhancers are compound regulatory elements that control expression of key cell identity genes. They recruit high levels of tissue-specific transcription factors and co-activators such as the Mediator complex and contact target gene promoters with high frequency. Most super-enhancers contain multiple constituent regulatory elements, but it is unclear whether these elements have distinct roles in activating target gene expression. Here, by rebuilding the endogenous multipartite α-globin super-enhancer, we show that it contains bioinformatically equivalent but functionally distinct element types: classical enhancers and facilitator elements. Facilitators have no intrinsic enhancer activity, yet in their absence, classical enhancers are unable to fully upregulate their target genes. Without facilitators, classical enhancers exhibit reduced Mediator recruitment, enhancer RNA transcription, and enhancer-promoter interactions. Facilitators are interchangeable but display functional hierarchy based on their position within a multipartite enhancer. Facilitators thus play an important role in potentiating the activity of classical enhancers and ensuring robust activation of target genes.

Spatial promoter-enhancer hubs in cancer: organization, regulation, and function

Transcriptional dysregulation is a hallmark of cancer and can be driven by altered enhancer landscapes. Recent studies in genome organization have revealed that multiple enhancers and promoters can spatially coalesce to form dynamic topological assemblies, known as promoter-enhancer hubs, which strongly correlate with elevated gene expression. In this review, we discuss the structure and complexity of promoter-enhancer hubs recently identified in multiple cancer types. We further discuss underlying mechanisms driving dysregulation of promoter-enhancer hubs and speculate on their functional role in pathogenesis. Understanding the role of promoter-enhancer hubs in transcriptional dysregulation can provide insight into new therapeutic approaches to target these complex features of genome organization.

Enhancers within the Ig V Gene Region Orchestrate Chromatin Topology and Regulate V Gene Rearrangement Frequency to Shape the B Cell Receptor Repertoire Specificities

Effective Ab-mediated responses depend on a highly diverse Ab repertoire with the ability to bind a wide range of epitopes in disease-causing agents. The generation of this repertoire depends on the somatic recombination of the variable (V), diversity (D), and joining (J) genes in the Ig loci of developing B cells. It has been known for some time that individual V, D, and J gene segments rearrange at different frequencies, but the mechanisms behind this unequal V gene usage have not been well understood. However, recent work has revealed that newly described enhancers scattered throughout the V gene-containing portion of the Ig loci regulate the V gene recombination frequency in a regional manner. Deletion of three of these enhancers revealed that these elements exert many layers of control during V(D)J recombination, including long-range chromatin interactions, epigenetic milieu, chromatin accessibility, and compartmentalization.

Epigenome editing in cancer: Advances and challenges for potential therapeutic options

Cancers are diseases caused by genetic and non-genetic environmental factors. Epigenetic alterations, some attributed to non-genetic factors, can lead to cancer development. Epigenetic changes can occur in tumor suppressors or oncogenes, or they may contribute to global cell state changes, making cells abnormal. Recent advances in gene editing technology show potential for cancer treatment. Herein, we will discuss our current knowledge of epigenetic alterations occurring in cancer and epigenetic editing technologies that can be applied to developing therapeutic options.

Be-1DCNN: a neural network model for chromatin loop prediction based on bagging ensemble learning

The chromatin loops in the three-dimensional (3D) structure of chromosomes are essential for the regulation of gene expression. Despite the fact that high-throughput chromatin capture techniques can identify the 3D structure of chromosomes, chromatin loop detection utilizing biological experiments is arduous and time-consuming. Therefore, a computational method is required to detect chromatin loops. Deep neural networks can form complex representations of Hi-C data and provide the possibility of processing biological datasets. Therefore, we propose a bagging ensemble one-dimensional convolutional neural network (Be-1DCNN) to detect chromatin loops from genome-wide Hi-C maps. First, to obtain accurate and reliable chromatin loops in genome-wide contact maps, the bagging ensemble learning method is utilized to synthesize the prediction results of multiple 1DCNN models. Second, each 1DCNN model consists of three 1D convolutional layers for extracting high-dimensional features from input samples and one dense layer for producing the prediction results. Finally, the prediction results of Be-1DCNN are compared to those of the existing models. The experimental results indicate that Be-1DCNN predicts high-quality chromatin loops and outperforms the state-of-the-art methods using the same evaluation metrics. The source code of Be-1DCNN is available for free at https://github.com/HaoWuLab-Bioinformatics/Be1DCNN.

Emerging insights into enhancer biology and function

Cell type-specific gene expression is coordinated by DNA-encoded enhancers and the transcription factors (TFs) that bind to them in a sequence-specific manner. As such, these enhancers and TFs are critical mediators of normal development and altered enhancer or TF function is associated with the development of diseases such as cancer. While initially defined by their ability to activate gene transcription in reporter assays, putative enhancer elements are now frequently defined by their unique chromatin features including DNase hypersensitivity and transposase accessibility, bidirectional enhancer RNA (eRNA) transcription, CpG hypomethylation, high H3K27ac and H3K4me1, sequence-specific transcription factor binding, and co-factor recruitment. Identification of these chromatin features through sequencing-based assays has revolutionized our ability to identify enhancer elements on a genome-wide scale, and genome-wide functional assays are now capitalizing on this information to greatly expand our understanding of how enhancers function to provide spatiotemporal coordination of gene expression programs. Here, we highlight recent technological advances that are providing new insights into the molecular mechanisms by which these critical cis-regulatory elements function in gene control. We pay particular attention to advances in our understanding of enhancer transcription, enhancer-promoter syntax, 3D organization and biomolecular condensates, transcription factor and co-factor dependencies, and the development of genome-wide functional enhancer screens.

Technological advances in probing 4D genome organization

The last two decades of work on chromosome conformation in eukaryotic nuclei have revealed a complex and highly regulated hierarchy of architectural features, from self-associating domains and compartmental interactions to locus-specific loops. Recent findings have shown that these structures are dynamic and heterogeneous, with emerging insights into the factors that shape them and implications for the control of transcription and other nuclear processes. Here, we review the latest advances in the DNA sequencing- and microscopy-based technologies for probing these features in space and time (4D) and discuss how they have been combined with complementary approaches such as genetic perturbations, protein and RNA measurements, and modeling to gain mechanistic insights about genome regulation across space and time.

scNanoHi-C: a single-cell long-read concatemer sequencing method to reveal high-order chromatin structures within individual cells

The high-order three-dimensional (3D) organization of regulatory genomic elements provides a topological basis for gene regulation, but it remains unclear how multiple regulatory elements across the mammalian genome interact within an individual cell. To address this, herein, we developed scNanoHi-C, which applies Nanopore long-read sequencing to explore genome-wide proximal high-order chromatin contacts within individual cells. We show that scNanoHi-C can reliably and effectively profile 3D chromatin structures and distinguish structure subtypes among individual cells. This method could also be used to detect genomic variations, including copy-number variations and structural variations, as well as to scaffold the de novo assembly of single-cell genomes. Notably, our results suggest that extensive high-order chromatin structures exist in active chromatin regions across the genome, and multiway interactions between enhancers and their target promoters were systematically identified within individual cells. Altogether, scNanoHi-C offers new opportunities to investigate high-order 3D genome structures at the single-cell level.

Unravelling animal-microbiota evolution on a chip

Whether and how microorganisms have shaped the evolution of their animal hosts is a major question in biology. Although many animal evolutionary processes appear to correlate with changes in their associated microbial communities, the mechanistic processes leading to these patterns and their causal relationships are still far from being resolved. Gut-on-a-chip models provide an innovative approach that expands beyond the potential of conventional microbiome profiling to study how different animals sense and react to microbes by comparing responses of animal intestinal tissue models to different microbial stimuli. This complementary knowledge can contribute to our understanding of how host genetic features facilitate or prevent different microbiomes from being assembled, and in doing so elucidate the role of host-microbiota interactions in animal evolution.

Three-way contact analysis characterizes the higher order organization of the Tcra locus

The generation of highly diverse antigen receptors in T and B lymphocytes relies on V(D)J recombination. The enhancer Eα has been implicated in regulating the accessibility of Vα and Jα genes through long-range interactions during rearrangements of the T-cell antigen receptor gene Tcra. However, direct evidence for Eα physically mediating the interaction of Vα and Jα genes is still lacking. In this study, we utilized the 3C-HTGTS assay, a chromatin interaction technique based on 3C, to analyze the higher order chromatin structure of the Tcra locus. Our analysis revealed the presence of sufficient information in the 3C-HTGTS data to detect multiway contacts. Three-way contact analysis of the Tcra locus demonstrated the co-occurrence of the proximal Jα genes, Vα genes and Eα in CD4+CD8+ double-positive thymocytes. Notably, the INT2-TEAp loop emerged as a prominent structure likely to be responsible for bringing the proximal Jα genes and the Vα genes into proximity. Moreover, the enhancer Eα utilizes this loop to establish physical proximity with the proximal Vα gene region. This study provides insights into the higher order chromatin structure of the Tcra locus, shedding light on the spatial organization of chromatin and its impact on V(D)J recombination.

Systematic identification of inter-chromosomal interaction networks supports the existence of RNA factories

Most studies of genome organization have focused on intra-chromosomal (cis) contacts because they harbor key features such as DNA loops and topologically associating domains. Inter-chromosomal (trans) contacts have received much less attention, and tools for interrogating potential biologically relevant trans structures are lacking. Here, we develop a computational framework to identify sets of loci that jointly interact in trans from Hi-C data. This method, trans-C, initiates probabilistic random walks with restarts from a set of seed loci to traverse an input Hi-C contact network, thereby identifying sets of trans-contacting loci. We validate trans-C in three increasingly complex models of established trans contacts: the Plasmodium falciparum var genes, the mouse olfactory receptor "Greek islands", and the human RBM20 cardiac splicing factory. We then apply trans-C to systematically test the hypothesis that genes co-regulated by the same trans-acting element (i.e., a transcription or splicing factor) co-localize in three dimensions to form "RNA factories" that maximize the efficiency and accuracy of RNA biogenesis. We find that many loci with multiple binding sites of the same transcription factor interact with one another in trans, especially those bound by transcription factors with intrinsically disordered domains. Similarly, clustered binding of a subset of RNA binding proteins correlates with trans interaction of the encoding loci. These findings support the existence of trans interacting chromatin domains (TIDs) driven by RNA biogenesis. Trans-C provides an efficient computational framework for studying these and other types of trans interactions, empowering studies of a poorly understood aspect of genome architecture.

SnapHiC-D: a computational pipeline to identify differential chromatin contacts from single-cell Hi-C data

Single-cell high-throughput chromatin conformation capture technologies (scHi-C) has been used to map chromatin spatial organization in complex tissues. However, computational tools to detect differential chromatin contacts (DCCs) from scHi-C datasets in development and through disease pathogenesis are still lacking. Here, we present SnapHiC-D, a computational pipeline to identify DCCs between two scHi-C datasets. Compared to methods designed for bulk Hi-C data, SnapHiC-D detects DCCs with high sensitivity and accuracy. We used SnapHiC-D to identify cell-type-specific chromatin contacts at 10 Kb resolution in mouse hippocampal and human prefrontal cortical tissues, demonstrating that DCCs detected in the hippocampal and cortical cell types are generally associated with cell-type-specific gene expression patterns and epigenomic features. SnapHiC-D is freely available at https://github.com/HuMingLab/SnapHiC-D.

iSMOD: an integrative browser for image-based single-cell multi-omics data

Genomic and transcriptomic image data, represented by DNA and RNA fluorescence in situ hybridization (FISH), respectively, together with proteomic data, particularly that related to nuclear proteins, can help elucidate gene regulation in relation to the spatial positions of chromatins, messenger RNAs, and key proteins. However, methods for image-based multi-omics data collection and analysis are lacking. To this end, we aimed to develop the first integrative browser called iSMOD (image-based Single-cell Multi-omics Database) to collect and browse comprehensive FISH and nucleus proteomics data based on the title, abstract, and related experimental figures, which integrates multi-omics studies focusing on the key players in the cell nucleus from 20 000+ (still growing) published papers. We have also provided several exemplar demonstrations to show iSMOD's wide applications-profiling multi-omics research to reveal the molecular target for diseases; exploring the working mechanism behind biological phenomena using multi-omics interactions, and integrating the 3D multi-omics data in a virtual cell nucleus. iSMOD is a cornerstone for delineating a global view of relevant research to enable the integration of scattered data and thus provides new insights regarding the missing components of molecular pathway mechanisms and facilitates improved and efficient scientific research.

EpiMCI: Predicting Multi-Way Chromatin Interactions from Epigenomic Signals

The recently emerging high-throughput Pore-C (HiPore-C) can identify whole-genome high-order chromatin multi-way interactions with an ultra-high output, contributing to deciphering three-dimensional (3D) genome organization. However, it also brings new challenges to relevant data analysis. To alleviate this problem, we proposed the EpiMCI, a model for multi-way chromatin interaction prediction based on a hypergraph neural network with epigenomic signals as the input. The EpiMCI integrated separate hyperedge representations with coupling hyperedge information and obtained AUCs of 0.981 and 0.984 in the GM12878 and K562 datasets, respectively, which outperformed the current available method. Moreover, the EpiMCI can be applied to denoise the HiPore-C data and improve the data quality efficiently. Furthermore, the vertex embeddings extracted from the EpiMCI reflected the global chromatin architecture accurately. The principal component analysis suggested that it was well aligned with the activities of genomic regions at the chromatin compartment level. Taken together, the EpiMCI can accurately predict multi-way chromatin interactions and can be applied to studies relying on chromatin architecture.

Chromosome-level organization of the regulatory genome in the Drosophila nervous system

Previous studies have identified topologically associating domains (TADs) as basic units of genome organization. We present evidence of a previously unreported level of genome folding, where distant TAD pairs, megabases apart, interact to form meta-domains. Within meta-domains, gene promoters and structural intergenic elements present in distant TADs are specifically paired. The associated genes encode neuronal determinants, including those engaged in axonal guidance and adhesion. These long-range associations occur in a large fraction of neurons but support transcription in only a subset of neurons. Meta-domains are formed by diverse transcription factors that are able to pair over long and flexible distances. We present evidence that two such factors, GAF and CTCF, play direct roles in this process. The relative simplicity of higher-order meta-domain interactions in Drosophila, compared with those previously described in mammals, allowed the demonstration that genomes can fold into highly specialized cell-type-specific scaffolds that enable megabase-scale regulatory associations.

High-sensitive nascent transcript sequencing reveals BRD4-specific control of widespread enhancer and target gene transcription

Gene transcription by RNA polymerase II (Pol II) is under control of promoters and distal regulatory elements known as enhancers. Enhancers are themselves transcribed by Pol II correlating with their activity. How enhancer transcription is regulated and coordinated with transcription at target genes has remained unclear. Here, we developed a high-sensitive native elongating transcript sequencing approach, called HiS-NET-seq, to provide an extended high-resolution view on transcription, especially at lowly transcribed regions such as enhancers. HiS-NET-seq uncovers new transcribed enhancers in human cells. A multi-omics analysis shows that genome-wide enhancer transcription depends on the BET family protein BRD4. Specifically, BRD4 co-localizes to enhancer and promoter-proximal gene regions, and is required for elongation activation at enhancers and their genes. BRD4 keeps a set of enhancers and genes in proximity through long-range contacts. From these studies BRD4 emerges as a general regulator of enhancer transcription that may link transcription at enhancers and genes.

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.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

Evaluating the role of the nuclear microenvironment in gene function by population-based modeling

The nuclear folding of chromosomes relative to nuclear bodies is an integral part of gene function. Here, we demonstrate that population-based modeling-from ensemble Hi-C data-provides a detailed description of the nuclear microenvironment of genes and its role in gene function. We define the microenvironment by the subnuclear positions of genomic regions with respect to nuclear bodies, local chromatin compaction, and preferences in chromatin compartmentalization. These structural descriptors are determined in single-cell models, thereby revealing the structural variability between cells. We demonstrate that the microenvironment of a genomic region is linked to its functional potential in gene transcription, replication, and chromatin compartmentalization. Some chromatin regions feature a strong preference for a single microenvironment, due to association with specific nuclear bodies in most cells. Other chromatin shows high structural variability, which is a strong indicator of functional heterogeneity. Moreover, we identify specialized nuclear microenvironments, which distinguish chromatin in different functional states and reveal a key role of nuclear speckles in chromosome organization. We demonstrate that our method produces highly predictive three-dimensional genome structures, which accurately reproduce data from a variety of orthogonal experiments, thus considerably expanding the range of Hi-C data analysis.

Cell-type specific transcriptional adaptations of nucleus accumbens interneurons to amphetamine

Parvalbumin-expressing (PV+) interneurons of the nucleus accumbens (NAc) play an essential role in the addictive-like behaviors induced by psychostimulant exposure. To identify molecular mechanisms of PV+ neuron plasticity, we isolated interneuron nuclei from the NAc of male and female mice following acute or repeated exposure to amphetamine (AMPH) and sequenced for cell type-specific RNA expression and chromatin accessibility. AMPH regulated the transcription of hundreds of genes in PV+ interneurons, and this program was largely distinct from that regulated in other NAc GABAergic neurons. Chromatin accessibility at enhancers predicted cell-type specific gene regulation, identifying transcriptional mechanisms of differential AMPH responses. Finally, we assessed expression of PV-enriched, AMPH-regulated genes in an Mecp2 mutant mouse strain that shows heightened behavioral sensitivity to psychostimulants to explore the functional importance of this transcriptional program. Together these data provide novel insight into the cell-type specific programs of transcriptional plasticity in NAc neurons that underlie addictive-like behaviors.

A role for SETD2 loss in tumorigenesis through DNA methylation dysregulation

SETD2-dependent H3 Lysine-36 trimethylation (H3K36me3) has been recently linked to the deposition of de-novo DNA methylation. SETD2 is frequently mutated in cancer, however, the functional impact of SETD2 loss and depletion on DNA methylation across cancer types and tumorigenesis is currently unknown. Here, we perform a pan-cancer analysis and show that both SETD2 mutation and reduced expression are associated with DNA methylation dysregulation across 21 out of the 24 cancer types tested. In renal cancer, these DNA methylation changes are associated with altered gene expression of oncogenes, tumour suppressors, and genes involved in neoplasm invasiveness, including TP53, FOXO1, and CDK4. This suggests a new role for SETD2 loss in tumorigenesis and cancer aggressiveness through DNA methylation dysregulation. Moreover, using a robust machine learning methodology, we develop and validate a 3-CpG methylation signature which is sufficient to predict SETD2 mutation status with high accuracy and correlates with patient prognosis.

RNA polymerase II dynamics shape enhancer-promoter interactions

How enhancers control target gene expression over long genomic distances remains an important unsolved problem. Here we investigated enhancer-promoter communication by integrating data from nucleosome-resolution genomic contact maps, nascent transcription and perturbations affecting either RNA polymerase II (Pol II) dynamics or the activity of thousands of candidate enhancers. Integration of new Micro-C experiments with published CRISPRi data demonstrated that enhancers spend more time in close proximity to their target promoters in functional enhancer-promoter pairs compared to nonfunctional pairs, which can be attributed in part to factors unrelated to genomic position. Manipulation of the transcription cycle demonstrated a key role for Pol II in enhancer-promoter interactions. Notably, promoter-proximal paused Pol II itself partially stabilized interactions. We propose an updated model in which elements of transcriptional dynamics shape the duration or frequency of interactions to facilitate enhancer-promoter communication.

Ultra-long-range interactions between active regulatory elements

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

Activation, decommissioning, and dememorization: enhancers in a life cycle

Spatiotemporal regulation of cell type-specific gene expression is essential to convert a zygote into a complex organism that contains hundreds of distinct cell types. A class of cis-regulatory elements called enhancers, which have the potential to enhance target gene transcription, are crucial for precise gene expression programs during development. Following decades of research, many enhancers have been discovered and how enhancers become activated has been extensively studied. However, the mechanisms underlying enhancer silencing are less well understood. We review current understanding of enhancer decommissioning and dememorization, both of which enable enhancer silencing. We highlight recent progress from genome-wide perspectives that have revealed the life cycle of enhancers and how its dynamic regulation underlies cell fate transition, development, cell regeneration, and epigenetic reprogramming.