Three-dimensional chromatin landscapes in T cell acute lymphoblastic leukemia

Oncogenic Enhancers in Leukemia

Although the study of leukemogenesis has traditionally focused on protein-coding genes, the role of enhancer dysregulation is becoming increasingly recognized. The advent of high-throughput sequencing, together with a better understanding of enhancer biology, has revealed how various genetic and epigenetic lesions produce oncogenic enhancers that drive transformation. These aberrations include translocations that lead to enhancer hijacking, point mutations that modulate enhancer activity, and copy number alterations that modify enhancer dosage. In this review, we describe these mechanisms in the context of leukemia and discuss potential therapeutic avenues to target these regulatory elements. Significance: Large-scale sequencing projects have uncovered recurrent gene mutations in leukemia, but the picture remains incomplete: some patients harbor no such aberrations, whereas others carry only a few that are insufficient to bring about transformation on their own. One of the missing pieces is enhancer dysfunction, which only recently has emerged as a critical driver of leukemogenesis. Knowledge of the various mechanisms of enhancer dysregulation is thus key for a complete understanding of leukemia and its causes, as well as the development of targeted therapies in the era of precision medicine.

Oncogenic dependency on SWI/SNF chromatin remodeling factors in T-cell acute lymphoblastic leukemia

T-cell acute lymphoblastic leukemia (T-ALL) is a hematological malignancy arising from immature thymocytes. Unlike well-known oncogenic transcription factors, such as NOTCH1 and MYC, the involvement of chromatin remodeling factors in T-ALL pathogenesis is poorly understood. Here, we provide compelling evidence on how SWI/SNF chromatin remodeling complex contributes to human T-ALL pathogenesis. Integrative analysis of transcriptomic and ATAC-Seq datasets revealed high expression of SMARCA4, one of the subunits of the SWI/SNF complex, in T-ALL patient samples and cell lines compared to normal T cells. Loss of SMARCA protein function resulted in apoptosis induction and growth inhibition in multiple T-ALL cell lines. ATAC-Seq analysis revealed a massive reduction in chromatin accessibility across the genome after the loss of SMARCA protein function. RUNX1 interacts with SMARCA4 protein and co-occupies the same genomic regions. Importantly, the NOTCH1-MYC pathway was primarily affected when SMARCA protein function was impaired, implicating SWI/SNF as a novel therapeutic target.

Clonal evolution of the 3D chromatin landscape in patients with relapsed pediatric B-cell acute lymphoblastic leukemia

Relapsed pediatric B-cell acute lymphoblastic leukemia (B-ALL) remains one of the leading causes of cancer mortality in children. We performed Hi-C, ATAC-seq, and RNA-seq on 12 matched diagnosis/relapse pediatric leukemia specimens to uncover dynamic structural variants (SVs) and 3D chromatin rewiring that may contribute to relapse. While translocations are assumed to occur early in leukemogenesis and be maintained throughout progression, we discovered novel, dynamic translocations and confirmed several fusion transcripts, suggesting functional and therapeutic relevance. Genome-wide chromatin remodeling was observed at all organizational levels: A/B compartments, TAD interactivity, and chromatin loops, including some loci shared by 25% of patients. Shared changes were found to drive the expression of genes/pathways previously implicated in resistance as well as novel therapeutic candidates, two of which (ATXN1 and MN1) we functionally validated. Overall, these results demonstrate chromatin reorganization under the selective pressure of therapy and offer the potential for discovery of novel therapeutic interventions.

Integrated analyses highlight interactions between the three-dimensional genome and DNA, RNA and epigenomic alterations in metastatic prostate cancer

The impact of variations in the three-dimensional structure of the genome has been recognized, but solid cancer tissue studies are limited. Here, we performed integrated deep Hi-C sequencing with matched whole-genome sequencing, whole-genome bisulfite sequencing, 5-hydroxymethylcytosine (5hmC) sequencing and RNA sequencing across a cohort of 80 biopsy samples from patients with metastatic castration-resistant prostate cancer. Dramatic differences were present in gene expression, 5-methylcytosine/5hmC methylation and in structural variation versus mutation rate between A and B (open and closed) chromatin compartments. A subset of tumors exhibited depleted regional chromatin contacts at the AR locus, linked to extrachromosomal circular DNA (ecDNA) and worse response to AR signaling inhibitors. We also identified topological subtypes associated with stark differences in methylation structure, gene expression and prognosis. Our data suggested that DNA interactions may predispose to structural variant formation, exemplified by the recurrent TMPRSS2-ERG fusion. This comprehensive integrated sequencing effort represents a unique clinical tumor resource.

Super-enhancer omics in stem cell

The hallmarks of stem cells, such as proliferation, self-renewal, development, differentiation, and regeneration, are critical to maintain stem cell identity which is sustained by genetic and epigenetic factors. Super-enhancers (SEs), which consist of clusters of active enhancers, play a central role in maintaining stemness hallmarks by specifically transcriptional model. The SE-navigated transcriptional complex, including SEs, non-coding RNAs, master transcriptional factors, Mediators and other co-activators, forms phase-separated condensates, which offers a toggle for directing diverse stem cell fate. With the burgeoning technologies of multiple-omics applied to examine different aspects of SE, we firstly raise the concept of "super-enhancer omics", inextricably linking to Pan-omics. In the review, we discuss the spatiotemporal organization and concepts of SEs, and describe links between SE-navigated transcriptional complex and stem cell features, such as stem cell identity, self-renewal, pluripotency, differentiation and development. We also elucidate the mechanism of stemness and oncogenic SEs modulating cancer stem cells via genomic and epigenetic alterations hijack in cancer stem cell. Additionally, we discuss the potential of targeting components of the SE complex using small molecule compounds, genome editing, and antisense oligonucleotides to treat SE-associated organ dysfunction and diseases, including cancer. This review also provides insights into the future of stem cell research through the paradigm of SEs.

CRISPR-inhibition screen for lncRNAs linked to melanoma growth and metastasis

Melanoma being one of the most common and deadliest skin cancers, has been rising since the past decade. Patients at advanced stages of the disease have very poor prognoses, as opposed to at the earlier stages. Nowadays the standard-of-care of advanced melanoma is resection followed by immune checkpoint inhibition based immunotherapy. However, a substantial proportion of patients either do not respond or develop resistances. This underscores a need for novel approaches and therapeutic targets as well as a better understanding of the mechanisms of melanoma pathogenesis. Long non-coding RNAs (lncRNAs) comprise a poorly characterized class of functional players and promising targets in promoting malignancy. Certain lncRNAs have been identified to play integral roles in melanoma progression and drug resistances, however systematic screens to uncover novel functional lncRNAs are scarce. Here, we profile differentially expressed lncRNAs in patient derived short-term metastatic cultures and BRAF-MEK-inhibition resistant cells. We conduct a focused growth-related CRISPR-inhibition screen of overexpressed lncRNAs, validate and functionally characterize lncRNA hits with respect to cellular growth, invasive capacities and apoptosis in vitro as well as the transcriptomic impact of our lead candidate the novel lncRNA XLOC_030781. In sum, we extend the current knowledge of ncRNAs and their potential relevance on melanoma.

Significance

Previously considered as transcriptional noise, lncRNAs have emerged as novel players in regulating many cellular aspects in health and disease including melanoma. However, the number and as well as the extent of functional significance of most lncRNAs remains elusive. We provide a comprehensive strategy to identify functionally relevant lncRNAs in melanoma by combining expression profiling with CRISPR-inhibition growths screens lowering the experimental effort. We also provide a larger resource of differentially expressed lncRNAs with potential implications in melanoma growth and invasion. Our results broaden the characterized of lncRNAs as potential targets for future therapeutic applications.

Connecting Chromatin Structures to Gene Regulation Using Dynamic Polymer Simulations

The transfer of regulatory information between distal loci on chromatin is thought to involve physical proximity, but key biophysical features of these contacts remain unclear. For instance, it is unknown how close and for how long two loci need to be in order to productively interact. The main challenge is that it is currently impossible to measure chromatin dynamics with high spatiotemporal resolution at scale. Polymer simulations provide an accessible and rigorous way to test biophysical models of chromatin regulation, yet there is a lack of simple and general methods for extracting the values of model parameters. Here we adapt the Nelder-Mead simplex optimization algorithm to select the best polymer model matching a given Hi-C dataset, using the MYC locus as an example. The model's biophysical parameters predict a compartmental rearrangement of the MYC locus in leukemia, which we validate with single-cell measurements. Leveraging trajectories predicted by the model, we find that loci with similar Hi-C contact frequencies can exhibit widely different contact dynamics. Interestingly, the frequency of productive interactions between loci exhibits a non-linear relationship with their Hi-C contact frequency when we enforce a specific capture radius and contact duration. These observations are consistent with recent experimental observations and suggest that the dynamic ensemble of chromatin configurations, rather than average contact matrices, is required to fully predict productive long-range chromatin interactions.

Enhancer-promoter hubs organize transcriptional networks promoting oncogenesis and drug resistance

Recent advances in high-resolution mapping of spatial interactions among regulatory elements support the existence of complex topological assemblies of enhancers and promoters known as enhancer-promoter hubs or cliques. Yet, organization principles of these multi-interacting enhancer-promoter hubs and their potential role in regulating gene expression in cancer remains unclear. Here, we systematically identified enhancer-promoter hubs in breast cancer, lymphoma, and leukemia. We found that highly interacting enhancer-promoter hubs form at key oncogenes and lineage-associated transcription factors potentially promoting oncogenesis of these diverse cancer types. Genomic and optical mapping of interactions among enhancer and promoter elements further showed that topological alterations in hubs coincide with transcriptional changes underlying acquired resistance to targeted therapy in T cell leukemia and B cell lymphoma. Together, our findings suggest that enhancer-promoter hubs are dynamic and heterogeneous topological assemblies with the potential to control gene expression circuits promoting oncogenesis and drug resistance.

3D Chromatin Alteration by Disrupting β-Catenin/CBP Interaction Is Enriched with Insulin Signaling in Pancreatic Cancer

The therapeutic potential of targeting the β-catenin/CBP interaction has been demonstrated in a variety of preclinical tumor models with a small molecule inhibitor, ICG-001, characterized as a β-catenin/CBP antagonist. Despite the high binding specificity of ICG-001 for the N-terminus of CBP, this β-catenin/CBP antagonist exhibits pleiotropic effects. Our recent studies found global changes in three-dimensional (3D) chromatin architecture in response to disruption of the β-catenin/CBP interaction in pancreatic cancer cells. However, an understanding of how the functional crosstalk between the antagonist and the β-catenin/CBP interaction affects changes in 3D chromatin architecture and, thereby, gene expression and downstream effects remains to be elucidated. Here, we perform Hi-C analyses on canonical and patient-derived pancreatic cancer cells before and after treatment with ICG-001. In addition to global alteration of 3D chromatin domains, we unexpectedly identify insulin signaling genes enriched in the altered chromatin domains. We further demonstrate that the chromatin loops associated with insulin signaling genes are significantly weakened after ICG-001 treatment. We finally elicit the deletion of a looping of IRS1-a key insulin signaling gene-significantly impeding pancreatic cancer cell growth, indicating that looping-mediated insulin signaling might act as an oncogenic pathway to promote pancreatic cancer progression. Our work shows that targeting aberrant insulin chromatin looping in pancreatic cancer might provide a therapeutic benefit.

Integrated enhancer regulatory network by enhancer-promoter looping in gastric cancer

Enhancer cis-regulatory elements play critical roles in gene regulation at many stages of cell growth. Enhancers in cancer cells also regulate the transcription of oncogenes. In this study, we performed a comprehensive analysis of long-range chromatin interactions, histone modifications, chromatin accessibility and expression in two gastric cancer (GC) cell lines compared to normal gastric epithelial cells. We found that GC-specific enhancers marked by histone modifications can activate a population of genes, including some oncogenes, by interacting with their proximal promoters. In addition, motif analysis of enhancer-promoter interacting enhancers showed that GC-specific transcription factors are enriched. Among them, we found that MYB is crucial for GC cell growth and activated by the enhancer with an enhancer-promoter loop and TCF7 upregulation. Clinical GC samples showed epigenetic activation of enhancers at the MYB locus and significant upregulation of TCF7 and MYB, regardless of molecular GC subtype and clinicopathological factors. Single-cell RNA sequencing of gastric mucosa with intestinal metaplasia showed high expression of TCF7 and MYB in intestinal stem cells. When we inactivated the loop-forming enhancer at the MYB locus using CRISPR interference (dCas9-KRAB), GC cell growth was significantly inhibited. In conclusion, we identified MYB as an oncogene activated by a loop-forming enhancer and contributing to GC cell growth.

Hijacked enhancer-promoter and silencer-promoter loops in cancer

Recent work has shown that besides inducing fusion genes, structural variations (SVs) can also contribute to oncogenesis by disrupting the three-dimensional genome organization and dysregulating gene expression. At the chromatin-loop level, SVs can relocate enhancers or silencers from their original genomic loci to activate oncogenes or repress tumor suppressor genes. On a larger scale, different types of alterations in topologically associating domains (TADs) have been reported in cancer, such as TAD expansion, shuffling, and SV-induced neo-TADs. Furthermore, the transformation from normal cells to cancerous cells is usually coupled with active or repressive compartmental switches, and cancer-specific compartments have been proposed. This review discusses the sites, and the other latest advances in studying how SVs disrupt higher-order genome structure in cancer, which in turn leads to oncogene dysregulation. We also highlight the clinical implications of these changes and the challenges ahead in this field.

Multiscale 3D genome organization underlies duck fatty liver with no adipose inflammation or serious injury

Non-alcoholic fatty liver disease (NAFLD) is the most common form of chronic liver disease. Little is known about how gene expression and chromatin structure are regulated in NAFLD due to lack of suitable model. Ducks naturally develop fatty liver similar to serious human non-alcoholic fatty liver (NAFL) without adipose inflammation and liver fibrosis, thus serves as a good model for investigating molecular mechanisms of adipose metabolism and anti-inflammation. Here, we constructed a NAFLD model without adipose inflammation and liver fibrosis in ducks. By performing dynamic pathological and transcriptomic analyses, we identified critical genes involving in regulation of the NF-κB and MHCII signaling, which usually lead to adipose inflammation and liver fibrosis. We further generated dynamic three-dimensional chromatin maps during liver fatty formation and recovery. This showed that ducks enlarged hepatocyte cell nuclei to reduce inter-chromosomal interaction, decompress chromatin structure, and alter strength of intra-TAD and loop interactions during fatty liver formation. These changes partially contributed to the tight control the NF-κB and the MHCII signaling. Our analysis uncovers duck chromatin reorganization might be advantageous to maintain liver regenerative capacity and reduce adipose inflammation. These findings shed light on new strategies for NAFLD control.

Probabilistic association of differentially expressed genes with cis-regulatory elements

Differential gene expression in response to perturbations is mediated at least in part by changes in binding of transcription factors (TFs) and other proteins at specific genomic regions. Association of these cis-regulatory elements (CREs) with their target genes is a challenging task that is essential to address many biological and mechanistic questions. Many current approaches rely on chromatin conformation capture techniques or single-cell correlational methods to establish CRE-to-gene associations. These methods can be effective but have limitations, including resolution, gaps in detectable association distances, and cost. As an alternative, we have developed DegCre, a nonparametric method that evaluates correlations between measurements of perturbation-induced differential gene expression and differential regulatory signal at CREs to score possible CRE-to-gene associations. It has several unique features, including the ability to use any type of CRE activity measurement, yield probabilistic scores for CRE-to-gene pairs, and assess CRE-to-gene pairings across a wide range of sequence distances. We apply DegCre to six data sets, each using different perturbations and containing a variety of regulatory signal measurements, including chromatin openness, histone modifications, and TF occupancy. To test their efficacy, we compare DegCre associations to Hi-C loop calls and CRISPR-validated CRE-to-gene associations, establishing good performance by DegCre that is comparable or superior to competing methods. DegCre is a novel approach to the association of CREs to genes from a perturbation-differential perspective, with strengths that are complementary to existing approaches and allow for new insights into gene regulation.

Pervasive structural heterogeneity rewires glioblastoma chromosomes to sustain patient-specific transcriptional programs

Glioblastoma multiforme (GBM) encompasses brain malignancies marked by phenotypic and transcriptional heterogeneity thought to render these tumors aggressive, resistant to therapy, and inevitably recurrent. However, little is known about how the spatial organization of GBM genomes underlies this heterogeneity and its effects. Here, we compile a cohort of 28 patient-derived glioblastoma stem cell-like lines (GSCs) known to reflect the properties of their tumor-of-origin; six of these were primary-relapse tumor pairs from the same patient. We generate and analyze 5 kbp-resolution chromosome conformation capture (Hi-C) data from all GSCs to systematically map thousands of standalone and complex structural variants (SVs) and the multitude of neoloops arising as a result. By combining Hi-C, histone modification, and gene expression data with chromatin folding simulations, we explain how the pervasive, uneven, and idiosyncratic occurrence of neoloops sustains tumor-specific transcriptional programs via the formation of new enhancer-promoter contacts. We also show how even moderately recurrent neoloops can relate to patient-specific vulnerabilities. Together, our data provide a resource for dissecting GBM biology and heterogeneity, as well as for informing therapeutic approaches.

CD-Loop: a chromatin loop detection method based on the diffusion model

Motivation

In recent years, there have been significant advances in various chromatin conformation capture techniques, and annotating the topological structure from Hi-C contact maps has become crucial for studying the three-dimensional structure of chromosomes. However, the structure and function of chromatin loops are highly dynamic and diverse, influenced by multiple factors. Therefore, obtaining the three-dimensional structure of the genome remains a challenging task. Among many chromatin loop prediction methods, it is difficult to fully extract features from the contact map and make accurate predictions at low sequencing depths.

Results

In this study, we put forward a deep learning framework based on the diffusion model called CD-Loop for predicting accurate chromatin loops. First, by pre-training the input data, we obtain prior probabilities for predicting the classification of the Hi-C contact map. Then, by combining the denoising process based on the diffusion model and the prior probability obtained by pre-training, candidate loops were predicted from the input Hi-C contact map. Finally, CD-Loop uses a density-based clustering algorithm to cluster the candidate chromatin loops and predict the final chromatin loops. We compared CD-Loop with the currently popular methods, such as Peakachu, Chromosight, and Mustache, and found that in different cell types, species, and sequencing depths, CD-Loop outperforms other methods in loop annotation. We conclude that CD-Loop can accurately predict chromatin loops and reveal cell-type specificity. The code is available at https://github.com/wangyang199897/CD-Loop.

Tet-mediated DNA methylation dynamics affect chromosome organization

DNA Methylation is a significant epigenetic modification that can modulate chromosome states, but its role in orchestrating chromosome organization has not been well elucidated. Here we systematically assessed the effects of DNA Methylation on chromosome organization with a multi-omics strategy to capture DNA Methylation and high-order chromosome interaction simultaneously on mouse embryonic stem cells with DNA methylation dioxygenase Tet triple knock-out (Tet-TKO). Globally, upon Tet-TKO, we observed weakened compartmentalization, corresponding to decreased methylation differences between CpG island (CGI) rich and poor domains. Tet-TKO could also induce hypermethylation for the CTCF binding peaks in TAD boundaries and chromatin loop anchors. Accordingly, CTCF peak generally weakened upon Tet-TKO, which results in weakened TAD structure and depletion of long-range chromatin loops. Genes that lost enhancer-promoter looping upon Tet-TKO showed DNA hypermethylation in their gene bodies, which may compensate for the disruption of gene expression. We also observed distinct effects of Tet1 and Tet2 on chromatin organization and increased DNA methylation correlation on spatially interacted fragments upon Tet inactivation. Our work showed the broad effects of Tet inactivation and DNA methylation dynamics on chromosome organization.

Hi-C profiling in tissues reveals 3D chromatin-regulated breast tumor heterogeneity and tumor-specific looping-mediated biological pathways

Current knowledge in three-dimensional (3D) chromatin regulation in normal and disease states was mostly accumulated through Hi-C profiling in in vitro cell culture system. The limitations include failing to recapitulate disease-specific physiological properties and often lacking clinically relevant disease microenvironment. In this study, we conduct tissue-specific Hi-C profiling in a pilot cohort of 12 breast tissues comprising of two normal tissues (NTs) and ten ER+ breast tumor tissues (TTs) including five primary tumors (PTs), and five tamoxifen-treated recurrent tumors (RTs). We find largely preserved compartments, highly heterogeneous topological associated domains (TADs) and intensively variable chromatin loops among breast tumors, demonstrating 3D chromatin-regulated breast tumor heterogeneity. Further cross-examination identifies RT-specific looping-mediated biological pathways and suggests CA2, an enhancer-promoter looping (EPL)-mediated target gene within the bicarbonate transport metabolism pathway, might play a role in driving the tamoxifen resistance. Remarkably, the inhibition of CA2 not only impedes tumor growth both in vitro and in vivo , but also reverses chromatin looping. Our study thus yields significant mechanistic insights into the role and clinical relevance of 3D chromatin architecture in breast cancer endocrine resistance.

A new chromosome-scale duck genome shows a major histocompatibility complex with several expanded multigene families

BACKGROUND

The duck (Anas platyrhynchos) is one of the principal natural hosts of influenza A virus (IAV), harbors almost all subtypes of IAVs and resists to many IAVs which cause extreme virulence in chicken and human. However, the response of duck's adaptive immune system to IAV infection is poorly characterized due to lack of a detailed gene map of the major histocompatibility complex (MHC).

RESULTS

We herein reported a chromosome-scale Beijing duck assembly by integrating Nanopore, Bionano, and Hi-C data. This new reference genome SKLA1.0 covers 40 chromosomes, improves the contig N50 of the previous duck assembly with highest contiguity (ZJU1.0) of more than a 5.79-fold, surpasses the chicken and zebra finch references in sequence contiguity and contains a complete genomic map of the MHC. Our 3D MHC genomic map demonstrated that gene family arrangement in this region was primordial; however, families such as AnplMHCI, AnplMHCIIβ, AnplDMB, NKRL (NK cell receptor-like genes) and BTN underwent gene expansion events making this area complex. These gene families are distributed in two TADs and genes sharing the same TAD may work in a co-regulated model.

CONCLUSIONS

These observations supported the hypothesis that duck's adaptive immunity had been optimized with expanded and diversified key immune genes which might help duck to combat influenza virus. This work provided a high-quality Beijing duck genome for biological research and shed light on new strategies for AIV control.

Topological reorganization and functional alteration of distinct genomic components in gallbladder cancer

Altered three-dimensional architecture of chromatin influences various genomic regulators and subsequent gene expression in human cancer. However, knowledge of the topological rearrangement of genomic hierarchical layers in cancer is largely limited. Here, by taking advantage of in situ Hi-C, RNA-sequencing, and chromatin immunoprecipitation sequencing (ChIP-seq), we investigated structural reorganization and functional changes in chromosomal compartments, topologically associated domains (TADs), and CCCTC binding factor (CTCF)-mediated loops in gallbladder cancer (GBC) tissues and cell lines. We observed that the chromosomal compartment A/B switch was correlated with CTCF binding levels and gene expression changes. Increased inter-TAD interactions with weaker TAD boundaries were identified in cancer cell lines relative to normal controls. Furthermore, the chromatin short loops and cancer unique loops associated with chromatin remodeling and epithelial-mesenchymal transition activation were enriched in cancer compared with their control counterparts. Cancer-specific enhancer-promoter loops, which contain multiple transcription factor binding motifs, acted as a central element to regulate aberrant gene expression. Depletion of individual enhancers in each loop anchor that connects with promoters led to the inhibition of their corresponding gene expressions. Collectively, our data offer the landscape of hierarchical layers of cancer genome and functional alterations that contribute to the development of GBC.

Epigenomic Characterization of Lymphoid Neoplasms

Lymphoid neoplasms represent a heterogeneous group of disease entities and subtypes with markedly different molecular and clinical features. Beyond genetic alterations, lymphoid tumors also show widespread epigenomic changes. These severely affect the levels and distribution of DNA methylation, histone modifications, chromatin accessibility, and three-dimensional genome interactions. DNA methylation stands out as a tracer of cell identity and memory, as B cell neoplasms show epigenetic imprints of their cellular origin and proliferative history, which can be quantified by an epigenetic mitotic clock. Chromatin-associated marks are informative to uncover altered regulatory regions and transcription factor networks contributing to the development of distinct lymphoid tumors. Tumor-intrinsic epigenetic and genetic aberrations cooperate and interact with microenvironmental cells to shape the transcriptome at different phases of lymphoma evolution, and intraclonal heterogeneity can now be characterized by single-cell profiling. Finally, epigenetics offers multiple clinical applications, including powerful diagnostic and prognostic biomarkers as well as therapeutic targets.

KnockTF 2.0: a comprehensive gene expression profile database with knockdown/knockout of transcription (co-)factors in multiple species

Transcription factors (TFs), transcription co-factors (TcoFs) and their target genes perform essential functions in diseases and biological processes. KnockTF 2.0 (http://www.licpathway.net/KnockTF/index.html) aims to provide comprehensive gene expression profile datasets before/after T(co)F knockdown/knockout across multiple tissue/cell types of different species. Compared with KnockTF 1.0, KnockTF 2.0 has the following improvements: (i) Newly added T(co)F knockdown/knockout datasets in mice, Arabidopsis thaliana and Zea mays and also an expanded scale of datasets in humans. Currently, KnockTF 2.0 stores 1468 manually curated RNA-seq and microarray datasets associated with 612 TFs and 172 TcoFs disrupted by different knockdown/knockout techniques, which are 2.5 times larger than those of KnockTF 1.0. (ii) Newly added (epi)genetic annotations for T(co)F target genes in humans and mice, such as super-enhancers, common SNPs, methylation sites and chromatin interactions. (iii) Newly embedded and updated search and analysis tools, including T(co)F Enrichment (GSEA), Pathway Downstream Analysis and Search by Target Gene (BLAST). KnockTF 2.0 is a comprehensive update of KnockTF 1.0, which provides more T(co)F knockdown/knockout datasets and (epi)genetic annotations across multiple species than KnockTF 1.0. KnockTF 2.0 facilitates not only the identification of functional T(co)Fs and target genes but also the investigation of their roles in the physiological and pathological processes.

Super-enhancers: Implications in gastric cancer

Gastric cancer (GC) is the fifth most prevalent malignancy and the third leading cause of cancer-related mortality globally. Despite intensive efforts to enhance the efficiencies of various therapeutics (chemotherapy, surgical interventions, molecular-targeted therapies, immunotherapies), the prognosis for patients with GC remains poor. This might be predominantly due to the limited understanding of the complicated etiology of GC. Importantly, epigenetic modifications and alterations are crucial during GC development. Super-enhancers (SEs) are a large cluster of adjacent enhancers that greatly activate transcription. SEs sustain cell-specific identity by enhancing the transcription of specific oncogenes. In this review, we systematically summarize how SEs are involved in GC development, including the SE landscape in GC, the SE target genes in GC, and the interventions related to SE functions for treating GC.

Super-enhancers and the super-enhancer reader BRD4: tumorigenic factors and therapeutic targets

Transcriptional super-enhancers and the BET bromodomain protein BRD4 are emerging as critical drivers of tumorigenesis and therapeutic targets. Characterized by substantial accumulation of histone H3 lysine 27 acetylation (H3K27ac) signals at the loci of cell identity genes and critical oncogenes, super-enhancers are recognized, bound and activated by BRD4, resulting in considerable oncogene over-expression, malignant transformation, cancer cell proliferation, survival, tumor initiation and progression. Small molecule compound BRD4 BD1 and BD2 bromodomain inhibitors block BRD4 binding to super-enhancers, suppress oncogene transcription and expression, reduce cancer cell proliferation and survival, and repress tumor progression in a variety of cancer types. Like other targeted therapy agents, BRD4 inhibitors show moderate anticancer effects on their own, and exert synergistic anticancer effects in vitro and in preclinical models, when combined with other anticancer agents including CDK7 inhibitors, CBP/p300 inhibitors and histone deacetylase inhibitors. More recently, BRD4 BD2 bromodomain selective inhibitors, proteolysis-targeting chimera (PROTAC) BRD4 protein degraders, and dual BRD4 and CBP/p300 bromodomain co-inhibitors have been developed and shown better anticancer efficacy and/or safety profile. Importantly, more than a dozen BRD4 inhibitors have entered clinical trials in patients with cancer of various organ origins. In summary, super-enhancers and their reader BRD4 are critical tumorigenic drivers, and BRD4 BD1 and BD2 bromodomain inhibitors, BRD4 BD2 bromodomain selective inhibitors, PROTAC BRD4 protein degraders, and dual BRD4 and CBP/p300 bromodomain co-inhibitors are promising novel anticancer agents for clinical translation.

Outward-oriented sites within clustered CTCF boundaries are key for intra-TAD chromatin interactions and gene regulation

CTCF plays an important role in 3D genome organization by adjusting the strength of chromatin insulation at TAD boundaries, where clustered CBS (CTCF-binding site) elements are often arranged in a tandem array with a complex divergent or convergent orientation. Here, using Pcdh and HOXD loci as a paradigm, we look into the clustered CTCF TAD boundaries and find that, counterintuitively, outward-oriented CBS elements are crucial for inward enhancer-promoter interactions as well as for gene regulation. Specifically, by combinatorial deletions of a series of putative enhancer elements in mice in vivo or CBS elements in cultured cells in vitro, in conjunction with chromosome conformation capture and RNA-seq analyses, we show that deletions of outward-oriented CBS elements weaken the strength of long-distance intra-TAD promoter-enhancer interactions and enhancer activation of target genes. Our data highlight the crucial role of outward-oriented CBS elements within the clustered CTCF TAD boundaries in developmental gene regulation and have interesting implications on the organization principles of clustered CTCF sites within TAD boundaries.

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.

Reprogramming of 3D chromatin domains by antagonizing the β-catenin/CBP interaction attenuates insulin signaling in pancreatic cancer

The therapeutic potential of targeting the β-catenin/CBP interaction has been demonstrated in a variety of preclinical tumor models with a small molecule inhibitor, ICG-001, characterized as a β-catenin/CBP antagonist. Despite the high binding specificity of ICG-001 for the N-terminus of CBP, this β-catenin/CBP antagonist exhibits pleiotropic effects. Our recent studies found global changes in three-dimensional (3D) chromatin architecture in response to disruption of the β-catenin/CBP interaction in pancreatic cancer cells. However, an understanding of the functional crosstalk between antagonizing the β-catenin/CBP interaction effect changes in 3D chromatin architecture and thereby gene expression and downstream effects remains to be elucidated. Here we perform Hi-C analyses on canonical and patient-derived pancreatic cancer cells before and after the treatment with ICG-001. In addition to global alteration of 3D chromatin domains, we unexpectedly identify insulin signaling genes enriched in the altered chromatin domains. We further demonstrate the chromatin loops associated with insulin signaling genes are significantly weakened after ICG-001 treatment. We finally elicit the deletion of a looping of IRS1, a key insulin signaling gene, significantly impede pancreatic cancer cell growth, indicating that looping-mediated insulin signaling might act as an oncogenic pathway to promote pancreatic cancer progression. Our work shows that targeting aberrant insulin chromatin looping in pancreatic cancer might provide a therapeutic benefit.

Enhancers dysfunction in the 3D genome of cancer cells

Eukaryotic genomes are spatially organized inside the cell nucleus, forming a threedimensional (3D) architecture that allows for spatial separation of nuclear processes and for controlled expression of genes required for cell identity specification and tissue homeostasis. Hence, it is of no surprise that mis-regulation of genome architecture through rearrangements of the linear genome sequence or epigenetic perturbations are often linked to aberrant gene expression programs in tumor cells. Increasing research efforts have shed light into the causes and consequences of alterations of 3D genome organization. In this review, we summarize the current knowledge on how 3D genome architecture is dysregulated in cancer, with a focus on enhancer highjacking events and their contribution to tumorigenesis. Studying the functional effects of genome architecture perturbations on gene expression in cancer offers a unique opportunity for a deeper understanding of tumor biology and sets the basis for the discovery of novel therapeutic targets.

Co-localization of clusters of TCR-regulated genes with TAD rearrangements

BACKGROUND

Gene expression has long been known to be influenced by the relative proximity of DNA regulatory elements. Topologically associating domains (TADs) are self-interacting genomic regions involved in regulating gene expression by controlling the proximity of these elements. Prior studies of TADs and their biological roles have revealed correlations between TAD changes and cellular differentiation. Here, we used Hi-C and RNA-seq data to correlate TCR-induced changes in TAD structure and gene expression in human CD4+ T cells.

RESULTS

We developed a pipeline, Differentially Expressed Gene Enrichment Finder (DEGEF), that identifies regions of differentially expressed gene enrichment. Using DEGEF, we found that TCR-regulated genes cluster non-uniformly across the genome and that these clusters preferentially localized in regions of TAD rearrangement. Interestingly, clusters of upregulated genes preferentially formed new Hi-C contacts compared to downregulated clusters, suggesting that TCR-activated CD4+ T cells may regulate genes by changing stimulatory contacts rather than inhibitory contacts.

CONCLUSIONS

Our observations support a significant relationship between TAD rearrangements and changes in local gene expression. These findings indicate potentially important roles for TAD rearrangements in shaping their local regulatory environments and thus driving differential expression of nearby genes during CD4+ T cell activation. Moreover, they provide new insights into global mechanisms that regulate gene expression.

Myeloid leukemia vulnerabilities embedded in long noncoding RNA locus MYNRL15

The noncoding genome presents a largely untapped source of new biological insights, including thousands of long noncoding RNA (lncRNA) loci. While lncRNA dysregulation has been reported in myeloid malignancies, their functional relevance remains to be systematically interrogated. We performed CRISPRi screens of lncRNA signatures from normal and malignant hematopoietic cells and identified MYNRL15 as a myeloid leukemia dependency. Functional dissection suggests an RNA-independent mechanism mediated by two regulatory elements embedded in the locus. Genetic perturbation of these elements triggered a long-range chromatin interaction and downregulation of leukemia dependency genes near the gained interaction sites, as well as overall suppression of cancer dependency pathways. Thus, this study describes a new noncoding myeloid leukemia vulnerability and mechanistic concept for myeloid leukemia. Importantly, MYNRL15 perturbation caused strong and selective impairment of leukemia cells of various genetic backgrounds over normal hematopoietic stem and progenitor cells in vitro, and depletion of patient-derived xenografts in vivo.

Perturbation of 3D nuclear architecture, epigenomic dysregulation and aging, and cannabinoid synaptopathy reconfigures conceptualization of cannabinoid pathophysiology: part 1-aging and epigenomics

Much recent attention has been directed toward the spatial organization of the cell nucleus and the manner in which three-dimensional topologically associated domains and transcription factories are epigenetically coordinated to precisely bring enhancers into close proximity with promoters to control gene expression. Twenty lines of evidence robustly implicate cannabinoid exposure with accelerated organismal and cellular aging. Aging has recently been shown to be caused by increased DNA breaks. These breaks rearrange and maldistribute the epigenomic machinery to weaken and reverse cellular differentiation, cause genome-wide DNA demethylation, reduce gene transcription, and lead to the inhibition of developmental pathways, which contribute to the progressive loss of function and chronic immune stimulation that characterize cellular aging. Both cell lineage-defining superenhancers and the superanchors that control them are weakened. Cannabis exposure phenocopies the elements of this process and reproduces DNA and chromatin breakages, reduces the DNA, RNA protein and histone synthesis, interferes with the epigenomic machinery controlling both DNA and histone modifications, induces general DNA hypomethylation, and epigenomically disrupts both the critical boundary elements and the cohesin motors that create chromatin loops. This pattern of widespread interference with developmental programs and relative cellular dedifferentiation (which is pro-oncogenic) is reinforced by cannabinoid impairment of intermediate metabolism (which locks in the stem cell-like hyper-replicative state) and cannabinoid immune stimulation (which perpetuates and increases aging and senescence programs, DNA damage, DNA hypomethylation, genomic instability, and oncogenesis), which together account for the diverse pattern of teratologic and carcinogenic outcomes reported in recent large epidemiologic studies in Europe, the USA, and elsewhere. It also accounts for the prominent aging phenotype observed clinically in long-term cannabis use disorder and the 20 characteristics of aging that it manifests. Increasing daily cannabis use, increasing use in pregnancy, and exponential dose-response effects heighten the epidemiologic and clinical urgency of these findings. Together, these findings indicate that cannabinoid genotoxicity and epigenotoxicity are prominent features of cannabis dependence and strongly indicate coordinated multiomics investigations of cannabinoid genome-epigenome-transcriptome-metabolome, chromatin conformation, and 3D nuclear architecture. Considering the well-established exponential dose-response relationships, the diversity of cannabinoids, and the multigenerational nature of the implications, great caution is warranted in community cannabinoid penetration.

In-depth understanding of higher-order genome architecture in orphan cancer

The human genome is intertwined, folded, condensed, and gradually constitutes the 3D architecture, thereby affecting transcription and widely involving in tumorigenesis. Incidence and mortality rates for orphan cancers increase due to poor early diagnosis and lack of effective medical treatments, which are now getting attention. In-depth understanding in tumorigenesis has fast-tracked over the last decade, however, the further role and mechanism of 3D genome organization in variant orphan tumorigenesis remains to be fully understood. We summarize for the first time that higher-order genome organization can provide novel insights into the occurrence mechanisms of orphan cancers, and discuss probable future research directions for drug development and anti-tumor therapies.

Multidimensional scaling methods can reconstruct genomic DNA loops using Hi-C data properties

This paper proposes multidimensional scaling (MDS) applied to high-throughput chromosome conformation capture (Hi-C) data on genomic interactions to visualize DNA loops. Currently, the mechanisms underlying the regulation of gene expression are poorly understood, and where and when DNA loops are formed remains undetermined. Previous studies have focused on reproducing the entire three-dimensional structure of chromatin; however, identifying DNA loops using these data is time-consuming and difficult. MDS is an unsupervised method for reconstructing the original coordinates from a distance matrix. Here, MDS was applied to high-throughput chromosome conformation capture (Hi-C) data on genomic interactions to visualize DNA loops. Hi-C data were converted to distances by taking the inverse to reproduce loops via MDS, and the missing values were set to zero. Using the converted data, MDS was applied to the log-transformed genomic coordinate distances and this process successfully reproduced the DNA loops in the given structure. Consequently, the reconstructed DNA loops revealed significantly more DNA-transcription factor interactions involved in DNA loop formation than those obtained from previously applied methods. Furthermore, the reconstructed DNA loops were significantly consistent with chromatin immunoprecipitation followed by sequencing (ChIP-seq) peak positions. In conclusion, the proposed method is an improvement over previous methods for identifying DNA loops.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.

Cell-type-specific prediction of 3D chromatin organization enables high-throughput in silico genetic screening

Investigating how chromatin organization determines cell-type-specific gene expression remains challenging. Experimental methods for measuring three-dimensional chromatin organization, such as Hi-C, are costly and have technical limitations, restricting their broad application particularly in high-throughput genetic perturbations. We present C.Origami, a multimodal deep neural network that performs de novo prediction of cell-type-specific chromatin organization using DNA sequence and two cell-type-specific genomic features-CTCF binding and chromatin accessibility. C.Origami enables in silico experiments to examine the impact of genetic changes on chromatin interactions. We further developed an in silico genetic screening approach to assess how individual DNA elements may contribute to chromatin organization and to identify putative cell-type-specific trans-acting regulators that collectively determine chromatin architecture. Applying this approach to leukemia cells and normal T cells, we demonstrate that cell-type-specific in silico genetic screening, enabled by C.Origami, can be used to systematically discover novel chromatin regulation circuits in both normal and disease-related biological systems.

Structural variation cooperates with permissive chromatin to control enhancer hijacking-mediated oncogenic transcription

Structural variants (SVs) involving enhancer hijacking can rewire chromatin topologies to cause oncogene activation in human cancers, including hematologic malignancies; however, because of the lack of tools to assess their effects on gene regulation and chromatin organization, the molecular determinants for the functional output of enhancer hijacking remain poorly understood. Here, we developed a multimodal approach to integrate genome sequencing, chromosome conformation, chromatin state, and transcriptomic alteration for quantitative analysis of transcriptional effects and structural reorganization imposed by SVs in leukemic genomes. We identified known and new pathogenic SVs, including recurrent t(5;14) translocations that cause the hijacking of BCL11B enhancers for the allele-specific activation of TLX3 in a subtype of pediatric leukemia. Epigenetic perturbation of SV-hijacked BCL11B enhancers impairs TLX3 transcription, which are required for the growth of t(5;14) leukemia cells. By CRISPR engineering of patient-derived t(5;14) in isogenic leukemia cells, we uncovered a new mechanism whereby the transcriptional output of SV-induced BCL11B enhancer hijacking is dependent on the loss of DNA hypermethylation at the TLX3 promoter. Our results highlight the importance of the cooperation between genetic alteration and permissive chromatin as a critical determinant of SV-mediated oncogene activation, with implications for understanding aberrant gene transcription after epigenetic therapies in patients with leukemia. Hence, leveraging the interdependency of genetic alteration on chromatin variation may provide new opportunities to reprogram gene regulation as targeted interventions in human disease.

Multi-omic analysis reveals dynamic changes of three-dimensional chromatin architecture during T cell differentiation

Cell differentiation results in widespread changes in transcriptional programs as well as multi-level remodeling of three-dimensional genome architecture. Nonetheless, few synthetically investigate the chromatin higher-order landscapes in different T helper (Th) cells. Using RNA-Seq, ATAC-Seq and Hi-C assays, we characterize dynamic changes in chromatin organization at different levels during Naive CD4⁺ T cells differentiation into T helper 17 (Th17) and T helper 1 (Th1) cells. Upon differentiation, we observe decreased short-range and increased extra-long-range chromatin interactions. Although there is no apparent global switch in the A/B compartments, Th cells display the weaker compartmentalization. A portion of topologically associated domains are rearranged. Furthermore, we identify cell-type specific enhancer-promoter loops, many of which are associated with functional genes in Th cells, such as Rorc facilitating Th17 differentiation and Hif1a responding to intracellular oxygen levels in Th1. Taken together, these results uncover the general patterns of chromatin reorganization and epigenetic landscapes of gene regulation during T helper cell differentiation.

Integrative CRISPR Activation and Small Molecule Inhibitor Screening for lncRNA Mediating BRAF Inhibitor Resistance in Melanoma

The incidence of melanoma, being one of the most commonly occurring cancers, has been rising since the past decade. Patients at advanced stages of the disease have very poor prognoses, as opposed to at the earlier stages. The conventional targeted therapy is well defined and effective for advanced-stage melanomas for patients not responding to the standard-of-care immunotherapy. However, targeted therapies do not prove to be as effective as patients inevitably develop V-Raf Murine Sarcoma Viral Oncogene Homolog B (BRAF)-inhibitor resistance to the respective drugs. Factors which are driving melanoma drug resistance mainly involve mutations in the mitogen-activated protein kinase (MAPK) pathway, e.g., BRAF splice variants, neuroblastoma RAS viral oncogene homolog (NRAS) amplification or parallel survival pathways. However, those mechanisms do not explain all cases of occurring resistances. Therefore, other factors accounting for BRAFi resistance must be better understood. Among them there are long non-coding RNAs (lncRNAs), but these remain functionally poorly understood. Here, we conduct a comprehensive, unbiased, and integrative study of lncRNA expression, coupled with a Clustered Regularly Interspaced Short Palindromic Repeats/Cas9-mediated activation (CRISPRa) and small molecule inhibitor screening for BRAF inhibitor resistance to expand the knowledge of potentially druggable lncRNAs, their function, and pave the way for eventual combinatorial treatment approaches targeting diverse pathways in melanoma.

Systematic mapping and modeling of 3D enhancer-promoter interactions in early mouse embryonic lineages reveal regulatory principles that determine the levels and cell-type specificity of gene expression

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages, the trophectoderm (TE), the epiblast (EPI) and the primitive endoderm (PrE). Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements via which transcriptional regulators enact these fates remain understudied. To address this gap, we have characterized, 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 observed extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although there are distinct groups of genes that 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 novel predictive model for transcriptional regulation (3D-HiChAT), which outperformed models that use only 1D promoter or proximal variables in predicting levels and cell-type specificity of gene expression. Using 3D-HiChAT, we performed genome-wide in silico perturbations to nominate candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments we validated several novel enhancers that control expression of one or more genes in their respective lineages. Our study comprehensively identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to understand lineage-specific transcriptional behaviors.

Cannabis- and Substance-Related Carcinogenesis in Europe: A Lagged Causal Inferential Panel Regression Study

Recent European data facilitate an epidemiological investigation of the controversial cannabis-cancer relationship. Of particular concern were prior findings associating high-dose cannabis use with reproductive problems and potential genetic impacts. Cancer incidence data age-standardised to the world population was obtained from the European Cancer Information System 2000-2020 and many European national cancer registries. Drug use data were obtained from the European Monitoring Centre for Drugs and Drug Addiction. Alcohol and tobacco consumption was sourced from the WHO. Median household income was taken from the World bank. Cancer rates in high-cannabis-use countries were significantly higher than elsewhere (β-estimate = 0.4165, p = 3.54 × 10^-115). Eighteen of forty-one cancers (42,675 individual rates) were significantly associated with cannabis exposure at bivariate analysis. Twenty-five cancers were linked in inverse-probability-weighted multivariate models. Temporal lagging in panel models intensified these effects. In multivariable models, cannabis was a more powerful correlate of cancer incidence than tobacco or alcohol. Reproductive toxicity was evidenced by the involvement of testis, ovary, prostate and breast cancers and because some of the myeloid and lymphoid leukaemias implicated occur in childhood, indicating inherited intergenerational genotoxicity. Cannabis is a more important carcinogen than tobacco and alcohol and fulfills epidemiological qualitative and quantitative criteria for causality for 25/41 cancers. Reproductive and transgenerational effects are prominent. These findings confirm the clinical and epidemiological salience of cannabis as a major multigenerational community carcinogen.

Spatial and clonality-resolved 3D cancer genome alterations reveal enhancer-hijacking as a potential prognostic marker for colorectal cancer

The regulatory effect of non-coding large-scale structural variations (SVs) on proto-oncogene activation remains unclear. This study investigated SV-mediated gene dysregulation by profiling 3D cancer genome maps from 40 patients with colorectal cancer (CRC). We developed a machine learning-based method for spatial characterization of the altered 3D cancer genome. This revealed a frequent establishment of "de novo chromatin contacts" that can span multiple topologically associating domains (TADs) in addition to the canonical TAD fusion/shuffle model. Using this information, we precisely identified super-enhancer (SE)-hijacking and its clonal characteristics. Clonal SE-hijacking genes, such as TOP2B, are recurrently associated with cell-cycle/DNA-processing functions, which can potentially be used as CRC prognostic markers. Oncogene activation and increased drug resistance due to SE-hijacking were validated by reconstructing the patient's SV using CRISPR-Cas9. Collectively, the spatial and clonality-resolved analysis of the 3D cancer genome reveals regulatory principles of large-scale SVs in oncogene activation and their clinical implications.

Comorbidity genetic risk and pathways impact SARS-CoV-2 infection outcomes

Understanding the genetic risk and mechanisms through which SARS-CoV-2 infection outcomes and comorbidities interact to impact acute and long-term sequelae is essential if we are to reduce the ongoing health burdens of the COVID-19 pandemic. Here we use a de novo protein diffusion network analysis coupled with tissue-specific gene regulatory networks, to examine putative mechanisms for associations between SARS-CoV-2 infection outcomes and comorbidities. Our approach identifies a shared genetic aetiology and molecular mechanisms for known and previously unknown comorbidities of SARS-CoV-2 infection outcomes. Additionally, genomic variants, genes and biological pathways that provide putative causal mechanisms connecting inherited risk factors for SARS-CoV-2 infection and coronary artery disease and Parkinson's disease are identified for the first time. Our findings provide an in depth understanding of genetic impacts on traits that collectively alter an individual's predisposition to acute and post-acute SARS-CoV-2 infection outcomes. The existence of complex inter-relationships between the comorbidities we identify raises the possibility of a much greater post-acute burden arising from SARS-CoV-2 infection if this genetic predisposition is realised.

Oncogenic super-enhancers in cancer: mechanisms and therapeutic targets

Activation of oncogenes to sustain proliferative signaling and initiate metastasis are important hallmarks of cancer. Oncogenes are amplified or overexpressed in cancer cells and overexpression is often controlled at the level of transcription. Gene expression is tightly controlled by many cis-regulatory elements and trans-acting factors. Large clusters of enhancers known as "super-enhancers" drive robust expression of cell-fate determining transcription factors in cell identity. Cancer cells can take advantage of super-enhancers and become transcriptionally addicted to them leading to tumorigenesis and metastasis. Additionally, the cis-regulatory landscape of cancer includes aberrant super-enhancers that are not present in normal cells. The landscape of super-enhancers in cancer is characterized by high levels of histone H3K27 acetylation and bromodomain-containing protein 4 (BRD4), and Mediator complex. These chromatin features facilitate the identification of cancer type-specific and cell-type-specific super-enhancers that control the expression of important oncogenes to stimulate their growth. Disruption of super-enhancers via inhibiting BRD4 or other epigenetic proteins is a potential therapeutic option. Here, we will describe the discovery of super-enhancers and their unique characteristics compared to typical enhancers. Then, we will highlight how super-enhancer-associated genes contribute to cancer progression in different solid tumor types. Lastly, we will cover therapeutic targets and their epigenetic modulators.

Elevated enhancer-oncogene contacts and higher oncogene expression levels by recurrent CTCF inactivating mutations in acute T cell leukemia

Monoallelic inactivation of CCCTC-binding factor (CTCF) in human cancer drives altered methylated genomic states, altered CTCF occupancy at promoter and enhancer regions, and deregulated global gene expression. In patients with T cell acute lymphoblastic leukemia (T-ALL), we find that acquired monoallelic CTCF-inactivating events drive subtle and local genomic effects in nearly half of t(5; 14) (q35; q32.2) rearranged patients, especially when CTCF-binding sites are preserved in between the BCL11B enhancer and the TLX3 oncogene. These solitary intervening sites insulate TLX3 from the enhancer by inducing competitive looping to multiple binding sites near the TLX3 promoter. Reduced CTCF levels or deletion of the intervening CTCF site abrogates enhancer insulation by weakening competitive looping while favoring TLX3 promoter to BCL11B enhancer looping, which elevates oncogene expression levels and leukemia burden.

3D genome mapping identifies subgroup-specific chromosome conformations and tumor-dependency genes in ependymoma

Ependymoma is a tumor of the brain or spinal cord. The two most common and aggressive molecular groups of ependymoma are the supratentorial ZFTA-fusion associated and the posterior fossa ependymoma group A. In both groups, tumors occur mainly in young children and frequently recur after treatment. Although molecular mechanisms underlying these diseases have recently been uncovered, they remain difficult to target and innovative therapeutic approaches are urgently needed. Here, we use genome-wide chromosome conformation capture (Hi-C), complemented with CTCF and H3K27ac ChIP-seq, as well as gene expression and DNA methylation analysis in primary and relapsed ependymoma tumors, to identify chromosomal conformations and regulatory mechanisms associated with aberrant gene expression. In particular, we observe the formation of new topologically associating domains ('neo-TADs') caused by structural variants, group-specific 3D chromatin loops, and the replacement of CTCF insulators by DNA hyper-methylation. Through inhibition experiments, we validate that genes implicated by these 3D genome conformations are essential for the survival of patient-derived ependymoma models in a group-specific manner. Thus, this study extends our ability to reveal tumor-dependency genes by 3D genome conformations even in tumors that lack targetable genetic alterations.

3D genome alterations and editing in pathology

The human genome is folded into a multi-level 3D structure that controls many nuclear functions including gene expression. Recently, alterations in 3D genome organization were associated with several genetic diseases and cancer. As a consequence, experimental approaches are now being developed to modify the global 3D genome organization and that of specific loci. Here, we discuss emerging experimental approaches of 3D genome editing that may prove useful in biomedicine.

NSD2 E1099K drives relapse in pediatric acute lymphoblastic leukemia by disrupting 3D chromatin organization

BACKGROUND

The NSD2 p.E1099K (EK) mutation is shown to be enriched in patients with relapsed acute lymphoblastic leukemia (ALL), indicating a role in clonal evolution and drug resistance.

RESULTS

To uncover 3D chromatin architecture-related mechanisms underlying drug resistance, we perform Hi-C on three B-ALL cell lines heterozygous for NSD2 EK. The NSD2 mutation leads to widespread remodeling of the 3D genome, most dramatically in terms of compartment changes with a strong bias towards A compartment shifts. Systematic integration of the Hi-C data with previously published ATAC-seq, RNA-seq, and ChIP-seq data show an expansion in H3K36me2 and a shrinkage in H3K27me3 within A compartments as well as increased gene expression and chromatin accessibility. These results suggest that NSD2 EK plays a prominent role in chromatin decompaction through enrichment of H3K36me2. In contrast, we identify few changes in intra-topologically associating domain activity. While compartment changes vary across cell lines, a common core of decompacting loci are shared, driving the expression of genes/pathways previously implicated in drug resistance. We further perform RNA sequencing on a cohort of matched diagnosis/relapse ALL patients harboring the relapse-specific NSD2 EK mutation. Changes in patient gene expression upon relapse significantly correlate with core compartment changes, further implicating the role of NSD2 EK in genome decompaction.

CONCLUSIONS

In spite of cell-context-dependent changes mediated by EK, there appears to be a shared transcriptional program dependent on compartment shifts which could explain phenotypic differences across EK cell lines. This core program is an attractive target for therapeutic intervention.

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.

MYC-driven synthesis of Siglec ligands is a glycoimmune checkpoint

The Siglecs (sialic acid-binding immunoglobulin-like lectins) are glycoimmune checkpoint receptors that suppress immune cell activation upon engagement of cognate sialoglycan ligands. The cellular drivers underlying Siglec ligand production on cancer cells are poorly understood. We find the MYC oncogene causally regulates Siglec ligand production to enable tumor immune evasion. A combination of glycomics and RNA-sequencing of mouse tumors revealed the MYC oncogene controls expression of the sialyltransferase St6galnac4 and induces a glycan known as disialyl-T. Using in vivo models and primary human leukemias, we find that disialyl-T functions as a "don't eat me" signal by engaging macrophage Siglec-E in mice or the human ortholog Siglec-7, thereby preventing cancer cell clearance. Combined high expression of MYC and ST6GALNAC4 identifies patients with high-risk cancers and reduced tumor myeloid infiltration. MYC therefore regulates glycosylation to enable tumor immune evasion. We conclude that disialyl-T is a glycoimmune checkpoint ligand. Thus, disialyl-T is a candidate for antibody-based checkpoint blockade, and the disialyl-T synthase ST6GALNAC4 is a potential enzyme target for small molecule-mediated immune therapy.

High-throughput Pore-C reveals the single-allele topology and cell type-specificity of 3D genome folding

Canonical three-dimensional (3D) genome structures represent the ensemble average of pairwise chromatin interactions but not the single-allele topologies in populations of cells. Recently developed Pore-C can capture multiway chromatin contacts that reflect regional topologies of single chromosomes. By carrying out high-throughput Pore-C, we reveal extensive but regionally restricted clusters of single-allele topologies that aggregate into canonical 3D genome structures in two human cell types. We show that fragments in multi-contact reads generally coexist in the same TAD. In contrast, a concurrent significant proportion of multi-contact reads span multiple compartments of the same chromatin type over megabase distances. Synergistic chromatin looping between multiple sites in multi-contact reads is rare compared to pairwise interactions. Interestingly, the single-allele topology clusters are cell type-specific even inside highly conserved TADs in different types of cells. In summary, HiPore-C enables global characterization of single-allele topologies at an unprecedented depth to reveal elusive genome folding principles.

Recruitment of TRIM33 to cell-context specific PML nuclear bodies regulates nodal signaling in mESCs

TRIM33 is a chromatin reader required for mammalian mesendoderm differentiation after activation of Nodal signaling, while its role in mESCs is still elusive. Here, we report that TRIM33 co-localizes with promyelocytic leukemia nuclear bodies (PML-NBs) specifically in mESCs, to mediate Nodal signaling-directed transcription of Lefty1/2. We show that TRIM33 puncta formation in mESCs depends on PML and on specific assembly of PML-NBs. Moreover, TRIM33 and PML co-regulate Lefty1/2 expression in mESCs, with both PML protein and formation of mESCs-specific PML-NBs being required for TRIM33 recruitment to these loci, and PML-NBs directly associating with the Lefty1/2 loci. Finally, a TurboID proximity-labeling experiment confirmed that TRIM33 is highly enriched only in mESCs-specific PML-NBs. Thus, our study supports a model in which TRIM33 condensates regulate Nodal signaling-directed transcription in mESCs and shows that PML-NBs can recruit distinct sets of client proteins in a cell-context-dependent manner.