Three-dimensional analysis reveals altered chromatin interaction by enhancer inhibitors harbors TCF7L2-regulated cancer gene signature

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.

BRD4-binding enhancer promotes CRC progression by interacting with YY1 to activate the Wnt pathway through upregulation of TCF7L2

Colorectal carcinoma (CRC), one of the most life-threatening cancer types, is associated with aberrant expression of epigenetic modifiers and activation of the Wnt pathway. However, the role of epigenetic regulators in driving cancer cell proliferation and their potential as therapeutic targets affecting the Wnt pathway remain unclear. In this study, BRD4 was found to promote the progression of CRC both in vitro and in vivo. The expression of BRD4 correlated with shortened CRC patient survival. In addition, BRD4 function was strongly correlated with the Wnt pathway, but rather through regulation of TCF7L2 at transcriptional levels. BRD4 and H3K27ac have overlapping occupancies in the cis-regulatory elements of TCF7L2, suggesting enhancer-based epigenetic regulation. Numerous YY1 binding sites were found in the abovementioned region. YY1 recruited BRD4 to bind to cis-regulatory elements of TCF7L2, thereby regulating the expression of TCF7L2. Altogether, this study validates that BRD4 performs a canonical epigenetic regulatory function in CRC and can be used in the treatment of Wnt pathway-dependent CRC or other malignancies with clinically available bromodomain inhibitors.

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.

Interplay Between the Histone Variant H2A.Z and the Epigenome in Pancreatic Cancer

BACKGROUND

The oncogenic process is orchestrated by a complex network of chromatin remodeling elements that shape the cancer epigenome. Histone variant H2A.Z regulates DNA control elements such as promoters and enhancers in different types of cancer; however, the interplay between H2A.Z and the pancreatic cancer epigenome is unknown.

OBJECTIVE

This study analyzed the role of H2A.Z in different DNA regulatory elements.

METHODS

We performed Chromatin Immunoprecipitation Sequencing assays (ChiP-seq) with total H2A.Z and acetylated H2A.Z (acH2A.Z) antibodies and analyzed published data from ChIP-seq, RNA-seq, bromouridine labeling-UV and sequencing (BruUV-seq), Hi-C and ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) in the pancreatic cancer cell line PANC-1.

RESULTS

The results indicate that total H2A.Z facilitates the recruitment of RNA polymerase II and transcription factors at promoters and enhancers allowing the expression of pro-oncogenic genes. Interestingly, we demonstrated that H2A.Z is enriched in super-enhancers (SEs) contributing to the transcriptional activation of key genes implicated in tumor development. Importantly, we established that H2A.Z contributes to the three-dimensional (3D) genome organization of pancreatic cancer and that it is a component of the Topological Associated Domains (TADs) boundaries in PANC-1 and that total H2A.Z and acH2A.Z are associated with A and B compartments, respectively.

CONCLUSIONS

H2A.Z participates in the biology and development of pancreatic cancer by generating a pro-oncogenic transcriptome through its posttranslational modifications, interactions with different partners, and regulatory elements, contributing to the oncogenic 3D genome organization. These data allow us to understand the molecular mechanisms that promote an oncogenic transcriptome in pancreatic cancer mediated by H2A.Z.

Analysis of immune related gene expression profiles and immune cell components in patients with Barrett esophagus

Barrett's esophagus (BE) is a well-known precancerous condition of esophageal adenocarcinoma. However, the immune cells and immune related genes involved in BE development and progression are not fully understood. Therefore, our study attempted to investigate the roles of immune cells and immune related genes in BE patients. The raw gene expression data were downloaded from the GEO database. The limma package in R was used to screen differentially expressed genes (DEGs). Then we performed the least absolute shrinkage and selection operator (LASSO) and random forest (RF) analyses to screen key genes. The proportion of infiltrated immune cells was evaluated using the CIBERSORT algorithm between BE and normal esophagus (NE) samples. The spearman index was used to show the correlations of immune genes and immune cells. Receiver operating characteristic (ROC) curves were used to assess the diagnostic value of key genes in BE. A total of 103 differentially expressed immune-related genes were identified between BE samples and normal samples. Then, 7 genes (CD1A, LTF, FABP4, PGC, TCF7L2, INSR,SEMA3C) were obtained after Lasso analysis and RF modeling. CIBERSORT analysis revealed that resting CD4 T memory cells and gamma delta T cells were present at significantly lower levels in BE samples. Moreover, plasma cell and regulatory T cells were present at significantly higher levels in BE samples than in NE samples. INSR had the highest AUC values in ROC analysis. We identified 7 immune related genes and 4 different immune cells in our study, that may play vital roles in the occurrence and development of BE. Our findings improve the understanding of the molecular mechanisms of BE.

Chromatin interaction maps identify Wnt responsive cis-regulatory elements coordinating Paupar-Pax6 expression in neuronal cells

Central nervous system-expressed long non-coding RNAs (lncRNAs) are often located in the genome close to protein coding genes involved in transcriptional control. Such lncRNA-protein coding gene pairs are frequently temporally and spatially co-expressed in the nervous system and are predicted to act together to regulate neuronal development and function. Although some of these lncRNAs also bind and modulate the activity of the encoded transcription factors, the regulatory mechanisms controlling co-expression of neighbouring lncRNA-protein coding genes remain unclear. Here, we used high resolution NG Capture-C to map the cis-regulatory interaction landscape of the key neuro-developmental Paupar-Pax6 lncRNA-mRNA locus. The results define chromatin architecture changes associated with high Paupar-Pax6 expression in neurons and identify both promoter selective as well as shared cis-regulatory-promoter interactions involved in regulating Paupar-Pax6 co-expression. We discovered that the TCF7L2 transcription factor, a regulator of chromatin architecture and major effector of the Wnt signalling pathway, binds to a subset of these candidate cis-regulatory elements to coordinate Paupar and Pax6 co-expression. We describe distinct roles for Paupar in Pax6 expression control and show that the Paupar DNA locus contains a TCF7L2 bound transcriptional silencer whilst the Paupar transcript can act as an activator of Pax6. Our work provides important insights into the chromatin interactions, signalling pathways and transcription factors controlling co-expression of adjacent lncRNAs and protein coding genes in the brain.

High-resolution Hi-C maps highlight multiscale 3D epigenome reprogramming during pancreatic cancer metastasis

Background

Pancreatic cancer’s poor prognosis is caused by distal metastasis, which is associated with epigenetic changes. However, the role of the 3D epigenome in pancreatic cancer biology, especially its metastasis, remains unclear.

Methods

Here, we developed high-resolution 3D epigenomic maps of cells derived from normal pancreatic epithelium, primary and metastatic pancreatic cancer by in situ Hi-C, ChIP-seq, ATAC-seq, and RNA-seq to identify key genes involved in pancreatic cancer metastasis

Results

We found that A/B compartments, contact domains, and chromatin loops changed significantly in metastatic pancreatic cancer cells, which are associated with epigenetic state alterations. Moreover, we found that upregulated genes, which were located in switched compartments, changed contact domains, and metastasis-specific enhancer-promoter loops, were related to cancer metastasis and poor prognosis of patients with pancreatic cancer. We also found that transcription factors in specific enhancer-promoter loop formation were also associated with metastasis. Finally we demonstrated that LIPC, looped to metastasis-specific enhancers, could promote pancreatic cancer metastasis.

Conclusions

These results highlight the multiscale 3D epigenome reprogramming during pancreatic cancer metastasis and expand our knowledge of mechanisms of gene regulation during pancreatic cancer metastasis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13045-021-01131-0.

The 3D genomic landscape of differential response to EGFR/HER2 inhibition in endocrine-resistant breast cancer cells

BACKGROUND

Recent studies suggested that crosstalk between ERα and EGFR/HER2 pathways plays a critical role in mediating endocrine therapy resistance. Several inhibitors targeting EGFR/HER2 signaling, including FDA-approved lapatinib and gefitinib as well as a novel dual tyrosine kinase inhibitor (TKI) sapitinib, showed greater therapeutic efficacies. However, how 3D chromatin landscape responds to the inhibition of EGFR/HER2 pathway remains to be elucidated.

METHODS

In this study, we conducted in situ Hi-C and RNA-seq in two ERα+ breast cancer cell systems, 1) parental MCF7 cells and its associated tamoxifen-resistant MCF7TR cells; and 2) parental T47D cells and its associated tamoxifen-resistant T47DTR cells, before and after the treatment of sapitinib.

RESULTS

We identified differential responses in topologically associated domains (TADs), looping genes and expressed genes. Interestingly, we found that many differential TADs and looping genes are reversible after sapitinib treatment, indicating that EGFR/HER2 signaling may play a role in reshaping and rewiring the high order genome organization. We further examined and recapitulated the reversible looping genes in 3D spheroids of breast cancer cells, demonstrating that 3D cell culture spheroid of breast cancer cells could be a potential preclinical breast cancer model for studying 3D chromatin regulation.

CONCLUSIONS

Our study has provided significant insights into our understanding of 3D genomic landscape changes in response to EGFR/HER2 Inhibition in endocrine-resistant breast cancer cells. Our data provides a rich resource for further evaluating chromatin structural responses to EGFR/HER2 targeted therapies in endocrine-resistant breast cancer cells. Our analyses suggest that these alterations of chromatin structures and transcriptional programs may provide new avenues for intervention or designing of patient selection for targeted endocrine treatment.

Disruption of Broad Epigenetic Domains in PDAC Cells by HAT Inhibitors

The spreading of epigenetic domains has emerged as a distinguishing epigenomic phenotype for diverse cell types. In particular, clusters of H3K27ac- and H3K4me3-marked elements, referred to as super-enhancers, and broad H3K4me3 domains, respectively, have been linked to cell identity and disease states. Here, we characterized the broad domains from different pancreatic ductal adenocarcinoma (PDAC) cell lines that represent distinct histological grades. Our integrative genomic analysis found that human derived cell line models for distinct PDAC grades exhibit characteristic broad epigenetic features associated with gene expression patterns that are predictive of patient prognosis and provide insight into pancreatic cancer cell identity. In particular, we find that genes marked by overlapping Low-Grade broad domains correspond to an epithelial phenotype and hold potential as markers for patient stratification. We further utilize ChIP-seq to compare the effects of histone acetyltransferase (HAT) inhibitors to detect global changes in histone acetylation and methylation levels. We found that HAT inhibitors impact certain broad domains of pancreatic cancer cells. Overall, our results reveal potential roles for broad domains in cells from distinct PDAC grades and demonstrate the plasticity of particular broad epigenomic domains to epigenetic inhibitors.