Citations

Background Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal cancer with a 5-year survival rate of 12%. It has two major molecular subtypes: classical and basal, regulated by the master transcription factors (MTFs) GATA6 and ΔNp63, respectively. Objective This study sought to uncover the transcriptional regulatory mechanisms controlling PDAC subtype identity. Design We integrated primary tumour single-cell RNA-seq, patient-derived xenograft RNA-seq and multispectral imaging to identify MTF-dependent, subtype-specific markers. We created subtype-specific fluorescent reporter systems and conducted drug screenings to find actionable targets. We analysed chromatin accessibility (ATAC-seq), genome-wide occupancy (ChIP-seq) for epigenetic status (H3K27ac), MTFs (GATA6, ΔNp63), RNA polymerase II (Pol II), H3K4me3-anchored chromatin topology (HiChIP) and nascent RNA capture sequencing (PRO-seq). Additionally, we used nuclease-dead Cas9 (dCas9) to manipulate transcriptional regulatory mechanisms. Results Our approach identified glucocorticoid receptor (GR) agonists as agents that suppress the classical transcriptional programme by interacting with GATA6. GATA6 regulates classical-specific transcription through promoter-proximal pause release. Depletion of GATA6 increased Pol II occupancy at GATA6-bound enhancers and transcriptional start sites, stabilising enhancer–promoter interactions. Artificially inducing pausing at GATA6-bound enhancers with dCas9 abrogated target gene expression and induced pausing at both the enhancer and target gene promoter. Conversely, in basal PDAC ΔNp63 promotes Pol II recruitment and stabilises enhancer–promoter interactions. Conclusion This study provides new insights into the transcriptional control and role of GR agonists in controlling PDAC molecular subtype identity.

Abstract A major hurdle to the application of precision oncology in pancreatic cancer is the lack of molecular stratification approaches and targeted therapy for defined molecular subtypes. In this work, we sought to gain further insight and identify molecular and epigenetic signatures of the Basal-like A pancreatic ductal adenocarcinoma (PDAC) subgroup that can be applied to clinical samples for patient stratification and/or therapy monitoring. We generated and integrated global gene expression and epigenome mapping data from patient-derived xenograft models to identify subtype-specific enhancer regions that were validated in patient-derived samples. In addition, complementary nascent transcription and chromatin topology (HiChIP) analyses revealed a Basal-like A subtype-specific transcribed enhancer program in PDAC characterized by enhancer RNA (eRNA) production that is associated with more frequent chromatin interactions and subtype-specific gene activation. Importantly, we successfully confirmed the validity of eRNA detection as a possible histologic approach for PDAC patient stratification by performing RNA-ISH analyses for subtype-specific eRNAs on pathologic tissue samples. Thus, this study provides proof-of-concept that subtype-specific epigenetic changes relevant for PDAC progression can be detected at a single-cell level in complex, heterogeneous, primary tumor material. Implications: Subtype-specific enhancer activity analysis via detection of eRNAs on a single-cell level in patient material can be used as a potential tool for treatment stratification.

Objective Pancreatic ductal adenocarcinoma (PDAC) displays a remarkable propensity towards therapy resistance. However, molecular epigenetic and transcriptional mechanisms enabling this are poorly understood. In this study, we aimed to identify novel mechanistic approaches to overcome or prevent resistance in PDAC. Design We used in vitro and in vivo models of resistant PDAC and integrated epigenomic, transcriptomic, nascent RNA and chromatin topology data. We identified a JunD-driven subgroup of enhancers, called interactive hubs (iHUBs), which mediate transcriptional reprogramming and chemoresistance in PDAC. Results iHUBs display characteristics typical for active enhancers (H3K27ac enrichment) in both therapy sensitive and resistant states but exhibit increased interactions and production of enhancer RNA (eRNA) in the resistant state. Notably, deletion of individual iHUBs was sufficient to decrease transcription of target genes and sensitise resistant cells to chemotherapy. Overlapping motif analysis and transcriptional profiling identified the activator protein 1 (AP1) transcription factor JunD as a master transcription factor of these enhancers. JunD depletion decreased iHUB interaction frequency and transcription of target genes. Moreover, targeting either eRNA production or signaling pathways upstream of iHUB activation using clinically tested small molecule inhibitors decreased eRNA production and interaction frequency, and restored chemotherapy responsiveness in vitro and in vivo. Representative iHUB target genes were found to be more expressed in patients with poor response to chemotherapy compared with responsive patients. Conclusion Our findings identify an important role for a subgroup of highly connected enhancers (iHUBs) in regulating chemotherapy response and demonstrate targetability in sensitisation to chemotherapy.

Abstract Pancreatic ductal adenocarcinoma (PDAC) displays a dismal prognosis due to late diagnosis and high chemoresistance incidence. For advanced disease stages or patients with comorbidities, treatment options are limited to gemcitabine alone or in combination with other drugs. While gemcitabine resistance has been widely attributed to the levels of one of its targets, RRM1, the molecular consequences of gemcitabine resistance in PDAC remain largely elusive. Here we sought to identify genomic, epigenomic, and transcriptomic events associated with gemcitabine resistance in PDAC and their potential clinical relevance. We found that gemcitabine-resistant cells displayed a coamplification of the adjacent RRM1 and STIM1 genes. Interestingly, RRM1, but not STIM1, was required for gemcitabine resistance, while high STIM1 levels caused an increase in cytosolic calcium concentration. Higher STIM1-dependent calcium influx led to an impaired endoplasmic reticulum stress response and a heightened nuclear factor of activated T-cell activity. Importantly, these findings were confirmed in patient and patient-derived xenograft samples. Taken together, our study uncovers previously unknown biologically relevant molecular properties of gemcitabine-resistant tumors, revealing an undescribed function of STIM1 as a rheostat directing the effects of calcium signaling and controlling epigenetic cell fate determination. It further reveals the potential benefit of targeting STIM1-controlled calcium signaling and its downstream effectors in PDAC. Significance: Gemcitabine-resistant and some naïve tumors coamplify RRM1 and STIM1, which elicit gemcitabine resistance and induce a calcium signaling shift, promoting ER stress resistance and activation of NFAT signaling.

Abstract Recent studies have thoroughly described genome-wide expression patterns defining molecular subtypes of pancreatic ductal adenocarcinoma (PDAC), with different prognostic and predictive implications. Although the reversible nature of key regulatory transcription circuits defining the two extreme PDAC subtype lineages “classical” and “basal-like” suggests that subtype states are not permanently encoded but underlie a certain degree of plasticity, pharmacologically actionable drivers of PDAC subtype identity remain elusive. Here, we characterized the mechanistic and functional implications of the histone methyltransferase enhancer of zeste homolog 2 (EZH2) in controlling PDAC plasticity, dedifferentiation, and molecular subtype identity. Utilization of transgenic PDAC models and human PDAC samples linked EZH2 activity to PDAC dedifferentiation and tumor progression. Combined RNA- and chromatin immunoprecipitation sequencing studies identified EZH2 as a pivotal suppressor of differentiation programs in PDAC and revealed EZH2-dependent transcriptional repression of the classical subtype defining transcription factor Gata6 as a mechanistic basis for EZH2-dependent PDAC progression. Importantly, genetic or pharmacologic depletion of EZH2 sufficiently increased GATA6 expression, thus inducing a gene signature shift in favor of a less aggressive and more therapy-susceptible, classical PDAC subtype state. Consistently, abrogation of GATA6 expression in EZH2-deficient PDAC cells counteracted the acquisition of classical gene signatures and rescued their invasive capacities, suggesting that GATA6 derepression is critical to overcome PDAC progression in the context of EZH2 inhibition. Together, our findings link the EZH2-GATA6 axis to PDAC subtype identity and uncover EZH2 inhibition as an appealing strategy to induce subtype-switching in favor of a less aggressive PDAC phenotype. Significance: This study highlights the role of EZH2 in PDAC progression and molecular subtype identity and suggests EZH2 inhibition as a strategy to recalibrate GATA6 expression in favor of a less aggressive disease.

Epigenetic alterations play a central role in cancer development and progression. The acetylation of histone 3 at lysine 27 (H3K27ac) specifically marks active genes. While chromatin immunoprecipitation (ChIP) followed by next-generation sequencing (ChIP-seq) analyses are commonly performed in cell lines, only limited data are available from primary tumors. We therefore examined whether cancer-specific alterations in H3K27ac occupancy can be identified in primary rectal cancer. Tissue samples from primary rectal cancer and matched mucosa were obtained. ChIP-seq for H3K27ac was performed and differentially occupied regions were identified. The expression of selected genes displaying differential occupancy between tumor and mucosa were examined in gene expression data from an independent patient cohort. Differential expression of four proteins was further examined by immunohistochemistry. ChIP-seq for H3K27ac in primary rectal cancer and matched mucosa was successfully performed and revealed differential binding on 44 regions. This led to the identification of genes with increased H3K27ac, i.e., RIPK2, FOXQ1, KRT23, and EPHX4, which were also highly upregulated in primary rectal cancer in an independent dataset. The increased expression of these four proteins was confirmed by immunohistochemistry. This study demonstrates the feasibility of ChIP-seq-based epigenome mapping of primary rectal cancer and confirms the value of H3K27ac occupancy to predict gene expression differences.