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.

3D chromatin hubs as regulatory units of identity and survival in human acute leukemia

Cancer progression involves genetic and epigenetic changes that disrupt chromatin 3D organization, affecting enhancer-promoter interactions and promoting growth. Here, we provide an integrative approach, combining chromatin conformation, accessibility, and transcription analysis, validated by in silico and CRISPR-interference screens, to identify relevant 3D topologies in pediatric T cell leukemia (T-ALL and ETP-ALL). We characterize 3D hubs as regulatory centers for oncogenes and disease markers, linking them to biological processes like cell division, inflammation, and stress response. Single-cell mapping reveals heterogeneous gene activation in discrete epigenetic clones, aiding in patient stratification for relapse risk after chemotherapy. Finally, we identify MYB as a 3D hub regulator in leukemia cells and show that the targeting of key regulators leads to hub dissolution, thereby providing a novel and effective anti-leukemic strategy. Overall, our work demonstrates the relevance of studying oncogenic 3D hubs to better understand cancer biology and tumor heterogeneity and to propose novel therapeutic strategies.

Epigenetic regulation in cancer

Epigenetic modifications are defined as heritable changes in gene activity that do not involve changes in the underlying DNA sequence. The oncogenic process is driven by the accumulation of alterations that impact genome's structure and function. Genetic mutations, which directly disrupt the DNA sequence, are complemented by epigenetic modifications that modulate gene expression, thereby facilitating the acquisition of malignant characteristics. Principals among these epigenetic changes are shifts in DNA methylation and histone mark patterns, which promote tumor development and metastasis. Notably, the reversible nature of epigenetic alterations, as opposed to the permanence of genetic changes, positions the epigenetic machinery as a prime target in the discovery of novel therapeutics. Our review delves into the complexities of epigenetic regulation, exploring its profound effects on tumor initiation, metastatic behavior, metabolic pathways, and the tumor microenvironment. We place a particular emphasis on the dysregulation at each level of epigenetic modulation, including but not limited to, the aberrations in enzymes responsible for DNA methylation and histone modification, subunit loss or fusions in chromatin remodeling complexes, and the disturbances in higher-order chromatin structure. Finally, we also evaluate therapeutic approaches that leverage the growing understanding of chromatin dysregulation, offering new avenues for cancer treatment.

Involvement of CCCTC-binding factor in epigenetic regulation of cancer

A major global health burden continues to be borne by the complex and multifaceted disease of cancer. Epigenetic changes, which are essential for the emergence and spread of cancer, have drawn a huge amount of attention recently. The CCCTC-binding factor (CTCF), which takes part in a wide range of cellular processes including genomic imprinting, X chromosome inactivation, 3D chromatin architecture, local modifications of histone, and RNA polymerase II-mediated gene transcription, stands out among the diverse array of epigenetic regulators. CTCF not only functions as an architectural protein but also modulates DNA methylation and histone modifications. Epigenetic regulation of cancer has already been the focus of plenty of studies. Understanding the role of CTCF in the cancer epigenetic landscape may lead to the development of novel targeted therapeutic strategies for cancer. CTCF has already earned its status as a tumor suppressor gene by acting like a homeostatic regulator of genome integrity and function. Moreover, CTCF has a direct effect on many important transcriptional regulators that control the cell cycle, apoptosis, senescence, and differentiation. As we learn more about CTCF-mediated epigenetic modifications and transcriptional regulations, the possibility of utilizing CTCF as a diagnostic marker and therapeutic target for cancer will also increase. Thus, the current review intends to promote personalized and precision-based therapeutics for cancer patients by shedding light on the complex interplay between CTCF and epigenetic processes.

Highly connected 3D chromatin networks established by an oncogenic fusion protein shape tumor cell identity

Cell fate transitions observed in embryonic development involve changes in three-dimensional genomic organization that provide proper lineage specification. Whether similar events occur within tumor cells and contribute to cancer evolution remains largely unexplored. We modeled this process in the pediatric cancer Ewing sarcoma and investigated high-resolution looping and large-scale nuclear conformation changes associated with the oncogenic fusion protein EWS-FLI1. We show that chromatin interactions in tumor cells are dominated by highly connected looping hubs centered on EWS-FLI1 binding sites, which directly control the activity of linked enhancers and promoters to establish oncogenic expression programs. Conversely, EWS-FLI1 depletion led to the disassembly of these looping networks and a widespread nuclear reorganization through the establishment of new looping patterns and large-scale compartment configuration matching those observed in mesenchymal stem cells, a candidate Ewing sarcoma progenitor. Our data demonstrate that major architectural features of nuclear organization in cancer cells can depend on single oncogenes and are readily reversed to reestablish latent differentiation programs.

Structural variations in cancer and the 3D genome

Structural variations (SVs) affect more of the cancer genome than any other type of somatic genetic alteration but difficulties in detecting and interpreting them have limited our understanding. Clinical cancer sequencing also increasingly aims to detect SVs, leading to a widespread necessity to interpret their biological and clinical relevance. Recently, analyses of large whole-genome sequencing data sets revealed features that impact rates of SVs across the genome in different cancers. A striking feature has been the extent to which, in both their generation and their influence on the selective fitness of cancer cells, SVs are more specific to individual cancer types than other genetic alterations such as single-nucleotide variants. This Perspective discusses how the folding of the 3D genome, and differences in its folding across cell types, affect observed SV rates in different cancer types as well as how SVs can impact cancer cell fitness.

Proteomic Analysis of Lung Cancer Types—A Pilot Study

Lung cancer is the leading cause of tumor-related mortality, therefore significant effort is directed towards understanding molecular alterations occurring at the origin of the disease to improve current treatment options. The aim of our pilot-scale study was to carry out a detailed proteomic analysis of formalin-fixed paraffin-embedded tissue sections from patients with small cell or non-small cell lung cancer (adenocarcinoma, squamous cell carcinoma, and large cell carcinoma). Tissue surface digestion was performed on relatively small cancerous and tumor-adjacent normal regions and differentially expressed proteins were identified using label-free quantitative mass spectrometry and subsequent statistical analysis. Principal component analysis clearly distinguished cancerous and cancer adjacent normal samples, while the four lung cancer types investigated had distinct molecular profiles and gene set enrichment analysis revealed specific dysregulated biological processes as well. Furthermore, proteins with altered expression unique to a specific lung cancer type were identified and could be the targets of future studies.

STAG2 regulates interferon signaling in melanoma via enhancer loop reprogramming

The cohesin complex participates in the organization of 3D genome through generating and maintaining DNA loops. Stromal antigen 2 (STAG2), a core subunit of the cohesin complex, is frequently mutated in various cancers. However, the impact of STAG2 inactivation on 3D genome organization, especially the long-range enhancer-promoter contacts and subsequent gene expression control in cancer, remains poorly understood. Here we show that depletion of STAG2 in melanoma cells leads to expansion of topologically associating domains (TADs) and enhances the formation of acetylated histone H3 lysine 27 (H3K27ac)-associated DNA loops at sites where binding of STAG2 is switched to its paralog STAG1. We further identify Interferon Regulatory Factor 9 (IRF9) as a major direct target of STAG2 in melanoma cells via integrated RNA-seq, STAG2 ChIP-seq and H3K27ac HiChIP analyses. We demonstrate that loss of STAG2 activates IRF9 through modulating the 3D genome organization, which in turn enhances type I interferon signaling and increases the expression of PD-L1. Our findings not only establish a previously unknown role of the STAG2 to STAG1 switch in 3D genome organization, but also reveal a functional link between STAG2 and interferon signaling in cancer cells, which may enhance the immune evasion potential in STAG2-mutant cancer.

Profiling of 3D Genome Organization in Nasopharyngeal Cancer Needle Biopsy Patient Samples by a Modified Hi-C Approach

Nasopharyngeal cancer (NPC), a cancer derived from epithelial cells in the nasopharynx, is a cancer common in China, Southeast Asia, and Africa. The three-dimensional (3D) genome organization of nasopharyngeal cancer is poorly understood. A major challenge in understanding the 3D genome organization of cancer samples is the lack of a method for the characterization of chromatin interactions in solid cancer needle biopsy samples. Here, we developed Biop-C, a modified in situ Hi-C method using solid cancer needle biopsy samples. We applied Biop-C to characterize three nasopharyngeal cancer solid cancer needle biopsy patient samples. We identified topologically associated domains (TADs), chromatin interaction loops, and frequently interacting regions (FIREs) at key oncogenes in nasopharyngeal cancer from the Biop-C heatmaps. We observed that the genomic features are shared at some important oncogenes, but the patients also display extensive heterogeneity at certain genomic loci. On analyzing the super enhancer landscape in nasopharyngeal cancer cell lines, we found that the super enhancers are associated with FIREs and can be linked to distal genes via chromatin loops in NPC. Taken together, our results demonstrate the utility of our Biop-C method in investigating 3D genome organization in solid cancers.

3D chromatin architecture and epigenetic regulation in cancer stem cells

Dedifferentiation of cell identity to a progenitor-like or stem cell-like state with increased cellular plasticity is frequently observed in cancer formation. During this process, a subpopulation of cells in tumours acquires a stem cell-like state partially resembling to naturally occurring pluripotent stem cells that are temporarily present during early embryogenesis. Such characteristics allow these cancer stem cells (CSCs) to give rise to the whole tumour with its entire cellular heterogeneity and thereby support metastases formation while being resistant to current cancer therapeutics. Cancer development and progression are demarcated by transcriptional dysregulation. In this article, we explore the epigenetic mechanisms shaping gene expression during tumorigenesis and cancer stem cell formation, with an emphasis on 3D chromatin architecture. Comparing the pluripotent stem cell state and epigenetic reprogramming to dedifferentiation in cellular transformation provides intriguing insight to chromatin dynamics. We suggest that the 3D chromatin architecture could be used as a target for re-sensitizing cancer stem cells to therapeutics.

Scanning window analysis of non-coding regions within normal-tumor whole-genome sequence samples

Genomics has benefited from an explosion in affordable high-throughput technology for whole-genome sequencing. The regulatory and functional aspects in non-coding regions may be an important contributor to oncogenesis. Whole-genome tumor-normal paired alignments were used to examine the non-coding regions in five cancer types and two races. Both a sliding window and a binning strategy were introduced to uncover areas of higher than expected variation for additional study. We show that the majority of cancer associated mutations in 154 whole-genome sequences covering breast invasive carcinoma, colon adenocarcinoma, kidney renal papillary cell carcinoma, lung adenocarcinoma and uterine corpus endometrial carcinoma cancers and two races are found outside of the coding region (4 432 885 in non-gene regions versus 1 412 731 in gene regions). A pan-cancer analysis found significantly mutated windows (292 to 3881 in count) demonstrating that there are significant numbers of large mutated regions in the non-coding genome. The 59 significantly mutated windows were found in all studied races and cancers. These offer 16 regions ripe for additional study within 12 different chromosomes-2, 4, 5, 7, 10, 11, 16, 18, 20, 21 and X. Many of these regions were found in centromeric locations. The X chromosome had the largest set of universal windows that cluster almost exclusively in Xq11.1-an area linked to chromosomal instability and oncogenesis. Large consecutive clusters (super windows) were found (19 to 114 in count) providing further evidence that large mutated regions in the genome are influencing cancer development. We show remarkable similarity in highly mutated non-coding regions across both cancer and race.

Enhancer rewiring in tumors: an opportunity for therapeutic intervention

Enhancers are cis-regulatory sequences that fine-tune expression of their target genes in a spatiotemporal manner. They are recognized by sequence-specific transcription factors, which in turn recruit transcriptional coactivators that facilitate transcription by promoting assembly and activation of the basal transcriptional machinery. Their functional importance is underscored by the fact that they are often the target of genetic and nongenetic events in human disease that disrupt their sequence, interactome, activation potential, and/or chromatin environment. Dysregulation of transcription and addiction to transcriptional effectors that interact with and modulate enhancer activity are common features of cancer cells and are amenable to therapeutic intervention. Here, we discuss the current knowledge on enhancer biology, the broad spectrum of mechanisms that lead to their malfunction in tumor cells, and recent progress in developing drugs that efficaciously target their dependencies.

Revisiting 3D chromatin architecture in cancer development and progression

Cancer development and progression are demarcated by transcriptional dysregulation, which is largely attributed to aberrant chromatin architecture. Recent transformative technologies have enabled researchers to examine the genome organization at an unprecedented dimension and precision. In particular, increasing evidence supports the essential roles of 3D chromatin architecture in transcriptional homeostasis and proposes its alterations as prominent causes of human cancer. In this article, we will discuss the recent findings on enhancers, enhancer-promoter interaction, chromatin topology, phase separation and explore their potential mechanisms in shaping transcriptional dysregulation in cancer progression. In addition, we will propose our views on how to employ state-of-the-art technologies to decode the unanswered questions in this field. Overall, this article motivates the study of 3D chromatin architecture in cancer, which allows for a better understanding of its pathogenesis and develop novel approaches for diagnosis and treatment of cancer.

Toward understanding the dynamic state of 3D genome

The three-dimensional (3D) genome organization and its role in biological activities have been investigated for over a decade in the field of cell biology. Recent studies using live-imaging and polymer simulation have suggested that the higher-order chromatin structures are dynamic; the stochastic fluctuations of nucleosomes and genomic loci cannot be captured by bulk-based chromosome conformation capture techniques (Hi-C). In this review, we focus on the physical nature of the 3D genome architecture. We first describe how to decode bulk Hi-C data with polymer modeling. We then introduce our recently developed PHi-C method, a computational tool for modeling the fluctuations of the 3D genome organization in the presence of stochastic thermal noise. We also present another new method that analyzes the dynamic rheology property (represented as microrheology spectra) as a measure of the flexibility and rigidity of genomic regions over time. By applying these methods to real Hi-C data, we highlighted a temporal hierarchy embedded in the 3D genome organization; chromatin interaction boundaries are more rigid than the boundary interior, while functional domains emerge as dynamic fluctuations within a particular time interval. Our methods may bridge the gap between live-cell imaging and Hi-C data and elucidate the nature of the dynamic 3D genome organization.

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

The causal relationship between 3D chromatin domains and gene regulation has been of considerable debate in recent years. Initial Hi-C studies profiling the 3D chromatin structure of the genome described evolutionarily conserved Topologically Associating Domains (TADs) that correlated with gene expression. Subsequent evidence from mouse models and human disease directly linked TADs to gene regulation. However, a number of focused genetic and genome-wide studies questioned the relevance of 3D chromatin domains for orchestrating gene expression, ultimately yielding a more multi-layered view of 3D chromatin structure and gene regulation. We review the evidence for and against the importance of 3D chromatin structure for gene regulation and argue for a more comprehensive classification of regulatory chromatin domains that integrates 3D chromatin structure with genomic, functional, and evolutionary conservation.

Heat Shock Factor 1 Drives Intergenic Association of Its Target Gene Loci upon Heat Shock

Transcriptional induction of heat shock protein (HSP) genes is accompanied by dynamic changes in their 3D structure and spatial organization, yet the molecular basis for these phenomena remains unknown. Using chromosome conformation capture and single-cell imaging, we show that genes transcriptionally activated by Hsf1 specifically interact across chromosomes and coalesce into diffraction-limited intranuclear foci. Genes activated by the alternative stress regulators Msn2/Msn4, in contrast, do not interact among themselves nor with Hsf1 targets. Likewise, constitutively expressed genes, even those interposed between HSP genes, show no detectable interaction. Hsf1 forms discrete subnuclear puncta when stress activated, and these puncta dissolve in concert with transcriptional attenuation, paralleling the kinetics of HSP gene coalescence and dissolution. Nuclear Hsf1 and RNA Pol II are both necessary for intergenic HSP gene interactions, while DNA-bound Hsf1 is necessary and sufficient to drive heterologous gene coalescence. Our findings demonstrate that Hsf1 can dynamically restructure the yeast genome.

While gene repositioning is thought to be a general feature of transcription, Chowdhary et al. provide evidence that argues against this concept. The authors demonstrate that Hsf1-regulated genes in Saccharomyces cerevisiae distinctively coalesce into intranuclear foci upon their transcriptional activation, while those activated by alternative transcription factors do not.

Chemical and Light Inducible Epigenome Editing

The epigenome defines the unique gene expression patterns and resulting cellular behaviors in different cell types. Epigenome dysregulation has been directly linked to various human diseases. Epigenome editing enabling genome locus-specific targeting of epigenome modifiers to directly alter specific local epigenome modifications offers a revolutionary tool for mechanistic studies in epigenome regulation as well as the development of novel epigenome therapies. Inducible and reversible epigenome editing provides unique temporal control critical for understanding the dynamics and kinetics of epigenome regulation. This review summarizes the progress in the development of spatiotemporal-specific tools using small molecules or light as inducers to achieve the conditional control of epigenome editing and their applications in epigenetic research.

Distinct Classes of Chromatin Loops Revealed by Deletion of an RNA-Binding Region in CTCF

Mammalian genomes are folded into topologically associating domains (TADs), consisting of chromatin loops anchored by CTCF and cohesin. Some loops are cell-type specific. Here we asked whether CTCF loops are established by a universal or locus-specific mechanism. Investigating the molecular determinants of CTCF clustering, we found that CTCF self-association in vitro is RNase sensitive and that an internal RNA-binding region (RBRi) mediates CTCF clustering and RNA interaction in vivo. Strikingly, deleting the RBRi impairs about half of all chromatin loops in mESCs and causes deregulation of gene expression. Disrupted loop formation correlates with diminished clustering and chromatin binding of RBRi mutant CTCF, which in turn results in a failure to halt cohesin-mediated extrusion. Thus, CTCF loops fall into at least two classes: RBRi-independent and RBRi-dependent loops. We speculate that evidence for RBRi-dependent loops may provide a molecular mechanism for establishing cell-specific CTCF loops, potentially regulated by RNA(s) or other RBRi-interacting partners.

Determining cellular CTCF and cohesin abundances to constrain 3D genome models

Achieving a quantitative and predictive understanding of 3D genome architecture remains a major challenge, as it requires quantitative measurements of the key proteins involved. Here, we report the quantification of CTCF and cohesin, two causal regulators of topologically associating domains (TADs) in mammalian cells. Extending our previous imaging studies (Hansen et al., 2017), we estimate bounds on the density of putatively DNA loop-extruding cohesin complexes and CTCF binding site occupancy. Furthermore, co-immunoprecipitation studies of an endogenously tagged subunit (Rad21) suggest the presence of cohesin dimers and/or oligomers. Finally, based on our cell lines with accurately measured protein abundances, we report a method to conveniently determine the number of molecules of any Halo-tagged protein in the cell. We anticipate that our results and the established tool for measuring cellular protein abundances will advance a more quantitative understanding of 3D genome organization, and facilitate protein quantification, key to comprehend diverse biological processes.

Modeling three-dimensional genomic organization in evolution and pathogenesis

The regulation of gene expression is mediated via the complex three-dimensional (3D) conformation of the genetic material and its interactions with various intracellular factors. Various experimental and computational approaches have been developed in recent years for understating the relation between the 3D conformation of the genome and the phenotypes of cells in normal condition and diseases. In this review, we will discuss novel approaches for analyzing and modeling the 3D genomic conformation, focusing on deciphering disease-causing mutations that affect gene expression. We conclude that as this is a very challenging mission, an important direction should involve the comparative analysis of various 3D models from various organisms or cells.

Transcriptional Dysregulation of MYC Reveals Common Enhancer-Docking Mechanism

Transcriptional dysregulation of the MYC oncogene is among the most frequent events in aggressive tumor cells, and this is generally accomplished by acquisition of a super-enhancer somewhere within the 2.8 Mb TAD where MYC resides. We find that these diverse cancer-specific super-enhancers, differing in size and location, interact with the MYC gene through a common and conserved CTCF binding site located 2 kb upstream of the MYC promoter. Genetic perturbation of this enhancer-docking site in tumor cells reduces CTCF binding, super-enhancer interaction, MYC gene expression, and cell proliferation. CTCF binding is highly sensitive to DNA methylation, and this enhancer-docking site, which is hypomethylated in diverse cancers, can be inactivated through epigenetic editing with dCas9-DNMT. Similar enhancer-docking sites occur at other genes, including genes with prominent roles in multiple cancers, suggesting a mechanism by which tumor cell oncogenes can generally hijack enhancers. These results provide insights into mechanisms that allow a single target gene to be regulated by diverse enhancer elements in different cell types.