Dissection of transcriptional and cis-regulatory control of differentiation in human pancreatic cancer

Insights into the mechanisms driven by H3K4 KMTs in pancreatic cancer

Pancreatic cancer is a malignancy arising from the endocrine or exocrine compartment of this organ. Tumors from exocrine origin comprise over 90% of all pancreatic cancers diagnosed. Of these, pancreatic ductal adenocarcinoma (PDAC) is the most common histological subtype. The five-year survival rate for PDAC ranged between 5 and 9% for over four decades, and only recently saw a modest increase to ∼12-13%, making this a severe and lethal disease. Like other cancers, PDAC initiation stems from genetic changes. However, therapeutic targeting of PDAC genetic drivers has remained relatively unsuccessful, thus the focus in recent years has expanded to the non-genetic factors underlying the disease pathogenesis. Specifically, it has been proposed that dynamic changes in the epigenetic landscape promote tumor growth and metastasis. Emphasis has been given to the re-organization of enhancers, essential regulatory elements controlling oncogenic gene expression, commonly marked my histone 3 lysine 4 monomethylation (H3K4me1). H3K4me1 is typically deposited by histone lysine methyltransferases (KMTs). While well characterized as oncogenes in other cancer types, recent work has expanded the role of KMTs as tumor suppressor in pancreatic cancer. Here, we review the role and translational significance for PDAC development and therapeutics of KMTs.

The xenobiotic transporter ABCC4/MRP4 promotes epithelial mesenchymal transition in pancreatic cancer

The xenobiotic transporter ABCC4/MRP4 is highly expressed in pancreatic ductal adenocarcinoma (PDAC) and correlates with a more aggressive phenotype and metastatic propensity. Here, we show that ABCC4 promotes epithelial-mesenchymal transition (EMT) in PDAC, a hallmark process involving the acquisition of mesenchymal traits by epithelial cells, enhanced cell motility, and chemoresistance. Modulation of ABCC4 levels in PANC-1 and BxPC-3 cell lines resulted in the dysregulation of genes present in the EMT signature. Bioinformatic analysis on several cohorts including tumor samples, primary patient-derived cultured cells, patient-derived xenografts, and cell lines, revealed a positive correlation between ABCC4 expression and EMT markers. We also characterized the ABCC4 cistrome and identified four candidate clusters in the distal promoter and intron one that showed differential binding of pro-epithelial FOXA1 and pro-mesenchymal GATA2 transcription factors in low ABCC4-expressing HPAF-II and high ABCC4-expressing PANC-1 xenografts. HPAF-II xenografts showed exclusive binding of FOXA1, and PANC-1 xenografts exclusive binding of GATA2, at ABCC4 clusters, consistent with their low and high EMT phenotype respectively. Our results underscore ABCC4/MRP4 as a valuable prognostic marker and a potential therapeutic target to treat PDAC subtypes with prominent EMT features, such as the basal-like/squamous subtype, characterized by worse prognosis and no effective therapies.

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.

A mucus production programme promotes classical pancreatic ductal adenocarcinoma

OBJECTIVE

The optimal therapeutic response in cancer patients is highly dependent upon the differentiation state of their tumours. Pancreatic ductal adenocarcinoma (PDA) is a lethal cancer that harbours distinct phenotypic subtypes with preferential sensitivities to standard therapies. This study aimed to investigate intratumour heterogeneity and plasticity of cancer cell states in PDA in order to reveal cell state-specific regulators.

DESIGN

We analysed single-cell expression profiling of mouse PDAs, revealing intratumour heterogeneity and cell plasticity and identified pathways activated in the different cell states. We performed comparative analysis of murine and human expression states and confirmed their phenotypic diversity in specimens by immunolabeling. We assessed the function of phenotypic regulators using mouse models of PDA, organoids, cell lines and orthotopically grafted tumour models.

RESULTS

Our expression analysis and immunolabeling analysis show that a mucus production programme regulated by the transcription factor SPDEF is highly active in precancerous lesions and the classical subtype of PDA - the most common differentiation state. SPDEF maintains the classical differentiation and supports PDA transformation in vivo. The SPDEF tumour-promoting function is mediated by its target genes AGR2 and ERN2/IRE1β that regulate mucus production, and inactivation of the SPDEF programme impairs tumour growth and facilitates subtype interconversion from classical towards basal-like differentiation.

CONCLUSIONS

Our findings expand our understanding of the transcriptional programmes active in precancerous lesions and PDAs of classical differentiation, determine the regulators of mucus production as specific vulnerabilities in these cell states and reveal phenotype switching as a response mechanism to inactivation of differentiation states determinants.

Abrogation of KLF5 sensitizes BRCA1-proficient pancreatic cancer to PARP inhibition

Poly ADP-ribose polymerase (PARP) inhibitor monotherapies are selectively effective in patients with pancreatic, breast, prostate, and ovarian cancers with BRCA1 mutations. Cancer patients with more frequent wild-type BRCA show poor responses to PARP inhibitors. Moreover, patients who are initially sensitive to these inhibitors eventually respond poorly to drugs. In the present study, we discover that abrogation of Kruppel-like factor 5 (KLF5) significantly inhibits homologous recombination, which is the main mechanism for DNA double-stranded repair. Furthermore, the downregulation of KLF5 expression promotes the DNA damage induced by olaparib and significantly reduces the IC 50 of the RARP inhibitor in pancreatic cancer cells. Overexpression of BRCA1 reverses the above effects caused by silencing of KLF5. Olaparib combined with a KLF5 inhibitor has an enhanced cytotoxic effect. Mechanistically, we identify BRCA1 as a KLF5 target gene. BRCA1 is positively correlated with KLF5 in PDAC tissue. Our results indicate that inhibition of KLF5 may induce BRCAness in a larger pancreatic cancer subset with proficient BRCA. The combination of KLF5 inhibitors and PARP inhibitors provides a novel treatment strategy to enhance the sensitivity of BRCA1-proficient pancreatic cancer to PARP inhibitors.

Mapping functional to morphological variation reveals the basis of regional extracellular matrix subversion and nerve invasion in pancreatic cancer

Intratumor morphological heterogeneity of pancreatic ductal adenocarcinoma (PDAC) predicts clinical outcomes but is only partially understood at the molecular level. To elucidate the gene expression programs underpinning intratumor morphological variation in PDAC, we investigated and deconvoluted at single cell level the molecular profiles of histologically distinct clusters of PDAC cells. We identified three major morphological and functional variants that co-exist in varying proportions in all PDACs, display limited genetic diversity, and are associated with a distinct organization of the extracellular matrix: a glandular variant with classical ductal features; a transitional variant displaying abortive ductal structures and mixed endodermal and myofibroblast-like gene expression; and a poorly differentiated variant lacking ductal features and basement membrane, and showing neuronal lineage priming. Ex vivo and in vitro evidence supports the occurrence of dynamic transitions among these variants in part influenced by extracellular matrix composition and stiffness and associated with local, specifically neural, invasion.

Altered histone acetylation patterns in pancreatic cancer cell lines induce subtype‑specific transcriptomic and phenotypical changes

Pancreatic ductal adenocarcinoma (PDAC) is often diagnosed at advanced tumor stages with chemotherapy as the only treatment option. Transcriptomic analysis has defined a classical and basal‑like PDAC subtype, which are regulated by epigenetic modification. The present study aimed to determine if drug‑induced epigenetic reprogramming of pancreatic cancer cells affects PDAC subtype identity and chemosensitivity. Classical and basal‑like PDAC cell lines PaTu‑S, Capan‑1, Capan‑2, Colo357, PaTu‑T, PANC‑1 and MIAPaCa‑2, were treated for a short (up to 96 h) and long (up to 30 weeks) period with histone acetyltransferase (HAT) and histone deacetylase (HDAC) inhibitors. The cells were analyzed using gene expression approaches, immunoblot analysis, and various cell assays to assess cell characteristics, such as proliferation, colony formation, cell migration and sensitivity to chemotherapeutic drugs. Classical and basal‑like PDAC cell lines showed pronounced epigenetic regulation of subtype‑specific genes through acetylation of lysine 27 on Histone H3 (H3K27ac). Moreover, classical cell lines revealed a significantly decreased expression of HDAC2 and increased total levels of H3K27ac in comparison with the basal‑like cell lines. Following HAT inhibitor treatment, classical cell lines exhibited a loss of epithelial marker gene expression, decreased chemotherapy response gene score and increased cell migration in vitro, indicating a tumor‑promoting phenotype. HDAC inhibitor treatment, however, exerted minimal reprogramming effects in both subtypes. Epigenetic reprogramming of classical and basal‑like tumor cells did not have a major impact on gemcitabine response, although the gemcitabine transporter gene SLC29A1 (solute carrier family 29 member 1) was epigenetically regulated.

An alternative splicing signature defines the basal-like phenotype and predicts worse clinical outcome in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is characterized by extremely poor prognosis. PDAC presents with molecularly distinct subtypes, with the basal-like one being associated with enhanced chemoresistance. Splicing dysregulation contributes to PDAC; however, its involvement in subtype specification remains elusive. Herein, we uncover a subtype-specific splicing signature associated with prognosis in PDAC and the splicing factor Quaking (QKI) as a determinant of the basal-like signature. Single-cell sequencing analyses highlight QKI as a marker of the basal-like phenotype. QKI represses splicing events associated with the classical subtype while promoting basal-like events associated with shorter survival. QKI favors a plastic, quasi-mesenchymal phenotype that supports migration and chemoresistance in PDAC organoids and cell lines, and its expression is elevated in high-grade primary tumors and metastatic lesions. These studies identify a splicing signature that defines PDAC subtypes and indicate that QKI promotes an undifferentiated, plastic phenotype, which renders PDAC cells chemoresistant and adaptable to environmental changes.

Cancer Cell Biomechanical Properties Accompany Tspan8-Dependent Cutaneous Melanoma Invasion

The intrinsic biomechanical properties of cancer cells remain poorly understood. To decipher whether cell stiffness modulation could increase melanoma cells' invasive capacity, we performed both in vitro and in vivo experiments exploring cell stiffness by atomic force microscopy (AFM). We correlated stiffness properties with cell morphology adaptation and the molecular mechanisms underlying epithelial-to-mesenchymal (EMT)-like phenotype switching. We found that melanoma cell stiffness reduction was systematically associated with the acquisition of invasive properties in cutaneous melanoma cell lines, human skin reconstructs, and Medaka fish developing spontaneous MAP-kinase-induced melanomas. We observed a systematic correlation of stiffness modulation with cell morphological changes towards mesenchymal characteristic gains. We accordingly found that inducing melanoma EMT switching by overexpressing the ZEB1 transcription factor, a major regulator of melanoma cell plasticity, was sufficient to decrease cell stiffness and transcriptionally induce tetraspanin-8-mediated dermal invasion. Moreover, ZEB1 expression correlated with Tspan8 expression in patient melanoma lesions. Our data suggest that intrinsic cell stiffness could be a highly relevant marker for human cutaneous melanoma development.

HOXC6 drives a therapeutically targetable pancreatic cancer growth and metastasis pathway by regulating MSK1 and PPP2R2B

Pancreatic ductal adenocarcinoma (PDAC) is one of the deadliest cancers, which lacks effective therapies. Here, we demonstrate that the transcription factor, homeobox C6 (HOXC6), is overexpressed in most PDACs, and its inhibition blocks PDAC tumor growth and metastasis. HOXC6 transcriptionally activates tumor-promoting kinase MSK1 and suppresses tumor-inhibitory protein PPP2R2B in PDAC. HOXC6-induced PPP2R2B suppression causes mammalian target of rapamycin (mTOR) pathway activation, which facilitates PDAC growth. Also, MSK1 upregulation by HOXC6 is necessary for PDAC growth because of its ability to suppress apoptosis via its substrate DDX17. Combinatorial pharmacological inhibition of MSK1 and mTOR potently suppressed PDAC tumor growth and metastasis in PDAC mouse models. PDAC cells with acquired resistance to MSK1/mTOR-inhibitors displayed activated insulin-like growth factor 1 receptor (IGF1R) signaling and were successfully eradicated by IGF1R inhibitor. Furthermore, MEK inhibitor trametinib enhanced the efficacy of dual MSK1 and mTOR inhibition. Collectively, these results identify therapeutic vulnerabilities of PDAC and an approach to overcome acquired drug resistance to prolong therapeutic benefit.

Krüppel-like Factor 5 Plays an Important Role in the Pathogenesis of Chronic Pancreatitis

Chronic pancreatitis results in the formation of pancreatic intraepithelial neoplasia (PanIN) and poses a risk of developing pancreatic cancer. Our previous study demonstrated that Krüppel-like factor 5 (KLF5) is necessary for forming acinar-to-ductal metaplasia (ADM) in acute pancreatitis. Here, we investigated the role of KLF5 in response to chronic injury in the pancreas. Human tissues originating from chronic pancreatitis patients showed increased levels of epithelial KLF5. An inducible genetic model combining the deletion of Klf5 and the activation of KrasG12D mutant expression in pancreatic acinar cells together with chemically induced chronic pancreatitis was used. The chronic injury resulted in increased levels of KLF5 in both control and KrasG12D mutant mice. Furthermore, it led to numerous ADM and PanIN lesions and extensive fibrosis in the KRAS mutant mice. In contrast, pancreata with Klf5 loss (with or without KrasG12D) failed to develop ADM, PanIN, or significant fibrosis. Furthermore, the deletion of Klf5 reduced the expression level of cytokines and fibrotic components such as Il1b, Il6, Tnf, Tgfb1, Timp1, and Mmp9. Notably, using ChIP-PCR, we showed that KLF5 binds directly to the promoters of Il1b, Il6, and Tgfb1 genes. In summary, the inactivation of Klf5 inhibits ADM and PanIN formation and the development of pancreatic fibrosis.

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.

Revealing and harnessing CD39 for the treatment of colorectal cancer and liver metastases by engineered T cells

OBJECTIVE

Colorectal tumours are often densely infiltrated by immune cells that have a role in surveillance and modulation of tumour progression but are burdened by immunosuppressive signals, which might vary from primary to metastatic stages. Here, we deployed a multidimensional approach to unravel the T-cell functional landscape in primary colorectal cancers (CRC) and liver metastases, and genome editing tools to develop CRC-specific engineered T cells.

DESIGN

We paired high-dimensional flow cytometry, RNA sequencing and immunohistochemistry to describe the functional phenotype of T cells from healthy and neoplastic tissue of patients with primary and metastatic CRC and we applied lentiviral vectors (LV) and CRISPR/Cas9 genome editing technologies to develop CRC-specific cellular products.

RESULTS

We found that T cells are mainly localised at the front edge and that tumor-infiltrating T cells co-express multiple inhibitory receptors, which largely differ from primary to metastatic sites. Our data highlighted CD39 as the major driver of exhaustion in both primary and metastatic colorectal tumours. We thus simultaneously redirected T-cell specificity employing a novel T-cell receptor targeting HER-2 and disrupted the endogenous TCR genes (TCR editing (TCRED)) and the CD39 encoding gene (ENTPD1), thus generating TCREDENTPD1KOHER-2-redirected lymphocytes. We showed that the absence of CD39 confers to HER-2-specific T cells a functional advantage in eliminating HER-2+ patient-derived organoids in vitro and in vivo.

CONCLUSION

HER-2-specific CD39 disrupted engineered T cells are promising advanced medicinal products for primary and metastatic CRC.

Monoclonal antibodies binding to different epitopes of CD20 differentially sensitize DLBCL to different classes of chemotherapy

Introduction

Rituximab (R), an anti-CD20 monoclonal antibody (mAb) and the world's first approved antibody for oncology patients, was combined with the CHOP chemotherapy regimen and markedly improved the prognosis of all B- cell-derived lymphomas, the most common hematological malignancy worldwide. However, there is a 35% disease recurrence with no advancement in the first-line treatment since R was combined with the archetypal CHOP chemotherapy regimen nearly 30 years ago. There is evidence that R synergizes with chemotherapy, but the pharmacological interactions between R and CHOP or between newer anti-CD20 mAbs and CHOP remain largely unexplored.

Methods

We used in vitro models to score pharmacological interactions between R and CHOP across various lymphoma cell lines. We compared these pharmacological interactions to ofatumumab, a second-generation anti-CD20 mAb, and CHOP. Lastly, we used RNA-sequencing to characterize the transcriptional profiles induced by these two antibodies and potential molecular pathways that mediate their different effects.

Results

We discovered vast heterogeneity in the pharmacological interactions between R and CHOP in a way not predicted by the current clinical classification. We then discovered that R and ofatumumab differentially synergize with the cytotoxic and cytostatic capabilities of CHOP in separate distinct subsets of B-cell lymphoma cell lines, thereby expanding favorable immunochemotherapy interactions across a greater range of cell lines beyond those induced by R-CHOP. Lastly, we discovered these two mAbs differentially modulate genes enriched in the JNK and p38 MAPK family, which regulates apoptosis and proliferation.

Discussion

Our findings were completely unexpected because these mAbs were long considered to be biological and clinical equivalents but, in practice, may perform better than the other in a patient-specific manner. This finding may have immediate clinical significance because both immunochemotherapy combinations are already FDA-approved with no difference in toxicity across phase I, II, and III clinical trials. Therefore, this finding could inform a new precision medicine strategy to provide additional therapeutic benefit to patients with B-cell lymphoma using immunochemotherapy combinations that already meet the clinical standard of care.

A TEAD2-Driven Endothelial-Like Program Shapes Basal-Like Differentiation and Metastasis of Pancreatic Cancer

BACKGROUND & AIMS

Pancreatic ductal adenocarcinoma (PDA), with its highly metastatic propensity, is one of the most lethal subtypes of pancreatic cancer. Although recent large-scale transcriptomic studies have demonstrated that heterogeneous gene expressions play an essential role in determining molecular phenotypes of PDA, biological cues for and consequences of distinct transcriptional programs remain unclear.

METHODS

We developed an experimental model that enforces the transition of PDA cells toward a basal-like subtype. We combined epigenome and transcriptome analyses with extensive in vitro and in vivo evaluations of tumorigenicity to demonstrate the validity of basal-like subtype differentiation in association with endothelial-like enhancer landscapes via TEA domain transcription factor 2 (TEAD2). Finally, we used loss-of-function experiments to investigate the importance of TEAD2 in regulating reprogrammed enhancer landscape and metastasis in basal-like PDA cells.

RESULTS

Aggressive characteristics of the basal-like subtype are faithfully recapitulated in vitro and in vivo, demonstrating the physiological relevance of our model. Further, we showed that basal-like subtype PDA cells acquire a TEAD2-dependent proangiogenic enhancer landscape. Genetic and pharmacologic inhibitions of TEAD2 in basal-like subtype PDA cells impair their proangiogenic phenotypes in vitro and cancer progression in vivo. Last, we identify CD109 as a critical TEAD2 downstream mediator that maintains constitutively activated JAK-STAT signaling in basal-like PDA cells and tumors.

CONCLUSIONS

Our findings implicate a TEAD2-CD109-JAK/STAT axis in the basal-like differentiated pancreatic cancer cells and as a potential therapeutic vulnerability.

Defining the fine structure of promoter activity on a genome-wide scale with CISSECTOR

Classic promoter mutagenesis strategies can be used to study how proximal promoter regions regulate the expression of particular genes of interest. This is a laborious process, in which the smallest sub-region of the promoter still capable of recapitulating expression in an ectopic setting is first identified, followed by targeted mutation of putative transcription factor binding sites. Massively parallel reporter assays such as survey of regulatory elements (SuRE) provide an alternative way to study millions of promoter fragments in parallel. Here we show how a generalized linear model (GLM) can be used to transform genome-scale SuRE data into a high-resolution genomic track that quantifies the contribution of local sequence to promoter activity. This coefficient track helps identify regulatory elements and can be used to predict promoter activity of any sub-region in the genome. It thus allows in silico dissection of any promoter in the human genome to be performed. We developed a web application, available at cissector.nki.nl, that lets researchers easily perform this analysis as a starting point for their research into any promoter of interest.

The Role of Krüppel-like Factors in Pancreatic Physiology and Pathophysiology

Krüppel-like factors (KLFs) belong to the family of transcription factors with three highly conserved zinc finger domains in the C-terminus. They regulate homeostasis, development, and disease progression in many tissues. It has been shown that KLFs play an essential role in the endocrine and exocrine compartments of the pancreas. They are necessary to maintain glucose homeostasis and have been implicated in the development of diabetes. Furthermore, they can be a vital tool in enabling pancreas regeneration and disease modeling. Finally, the KLF family contains proteins that act as tumor suppressors and oncogenes. A subset of members has a biphasic function, being upregulated in the early stages of oncogenesis and stimulating its progression and downregulated in the late stages to allow for tumor dissemination. Here, we describe KLFs' function in pancreatic physiology and pathophysiology.

Sirtuin 6 is required for the integrated stress response and resistance to inhibition of transcriptional cyclin-dependent kinases

Pancreatic ductal adenocarcinoma (PDAC) is classified into two key subtypes, classical and basal, with basal PDAC predicting worse survival. Using in vitro drug assays, genetic manipulation experiments, and in vivo drug studies in human patient-derived xenografts (PDXs) of PDAC, we found that basal PDACs were uniquely sensitive to transcriptional inhibition by targeting cyclin-dependent kinase 7 (CDK7) and CDK9, and this sensitivity was recapitulated in the basal subtype of breast cancer. We showed in cell lines, PDXs, and publicly available patient datasets that basal PDAC was characterized by inactivation of the integrated stress response (ISR), which leads to a higher rate of global mRNA translation. Moreover, we identified the histone deacetylase sirtuin 6 (SIRT6) as a critical regulator of a constitutively active ISR. Using expression analysis, polysome sequencing, immunofluorescence, and cycloheximide chase experiments, we found that SIRT6 regulated protein stability by binding activating transcription factor 4 (ATF4) in nuclear speckles and protecting it from proteasomal degradation. In human PDAC cell lines and organoids as well as in murine PDAC genetically engineered mouse models where SIRT6 was deleted or down-regulated, we demonstrated that SIRT6 loss both defined the basal PDAC subtype and led to reduced ATF4 protein stability and a nonfunctional ISR, causing a marked vulnerability to CDK7 and CDK9 inhibitors. Thus, we have uncovered an important mechanism regulating a stress-induced transcriptional program that may be exploited with targeted therapies in particularly aggressive PDAC.

A Kaleidoscope of Keratin Gene Expression and the Mosaic of Its Regulatory Mechanisms

Keratins are a family of intermediate filament-forming proteins highly specific to epithelial cells. A combination of expressed keratin genes is a defining property of the epithelium belonging to a certain type, organ/tissue, cell differentiation potential, and at normal or pathological conditions. In a variety of processes such as differentiation and maturation, as well as during acute or chronic injury and malignant transformation, keratin expression undergoes switching: an initial keratin profile changes accordingly to changed cell functions and location within a tissue as well as other parameters of cellular phenotype and physiology. Tight control of keratin expression implies the presence of complex regulatory landscapes within the keratin gene loci. Here, we highlight patterns of keratin expression in different biological conditions and summarize disparate data on mechanisms controlling keratin expression at the level of genomic regulatory elements, transcription factors (TFs), and chromatin spatial structure.

ESE1/AGR2 axis antagonizes TGF-β-induced epithelial-mesenchymal transition in low-grade pancreatic cancer

Epithelium-specific ETS transcription factor 1 (ESE1) has been implicated in epithelial homeostasis, inflammation, as well as tumorigenesis, and cancer progression. However, numerous studies have reported contradictory roles-as an oncogene or a tumor suppressor of ESE1 in different cancers, and its function in the development and progression of pancreatic ductal adenocarcinoma (PDAC) has remained largely unexplored. Herein, we report that ESE1 was found upregulated in primary PDAC compared to normal pancreatic tissue, but high expression of ESE1 correlated to better relapse-free survival in patients with PDAC. Interestingly, ESE1 was found to exhibit dual roles in regulation of malignant properties of PDAC cells in that its overexpression promoted cell proliferation, whereas its downregulation enhanced epithelial-mesenchymal transition (EMT) phenotype. In the context of TGF-β-induced EMT, ESE1 is markedly downregulated at post-transcriptional level, and reconstituted ESE1 expression partially reversed TGF-β-induced EMT marker expression. Furthermore, we identify AGR2 as a novel transcriptional target of ESE1 that participates in TGF-β-induced EMT in PDAC. Collectively, our findings reveal an ESE1/AGR2 axis that interacts with TGF-β signaling to modulate EMT phenotype in PDAC.

Regulatory enhancer profiling of mesenchymal-type gastric cancer reveals subtype-specific epigenomic landscapes and targetable vulnerabilities

OBJECTIVE

Gastric cancer (GC) comprises multiple molecular subtypes. Recent studies have highlighted mesenchymal-subtype GC (Mes-GC) as a clinically aggressive subtype with few treatment options. Combining multiple studies, we derived and applied a consensus Mes-GC classifier to define the Mes-GC enhancer landscape revealing disease vulnerabilities.

DESIGN

Transcriptomic profiles of ~1000 primary GCs and cell lines were analysed to derive a consensus Mes-GC classifier. Clinical and genomic associations were performed across >1200 patients with GC. Genome-wide epigenomic profiles (H3K27ac, H3K4me1 and assay for transposase-accessible chromatin with sequencing (ATAC-seq)) of 49 primary GCs and GC cell lines were generated to identify Mes-GC-specific enhancer landscapes. Upstream regulators and downstream targets of Mes-GC enhancers were interrogated using chromatin immunoprecipitation followed by sequencing (ChIP-seq), RNA sequencing, CRISPR/Cas9 editing, functional assays and pharmacological inhibition.

RESULTS

We identified and validated a 993-gene cancer-cell intrinsic Mes-GC classifier applicable to retrospective cohorts or prospective single samples. Multicohort analysis of Mes-GCs confirmed associations with poor patient survival, therapy resistance and few targetable genomic alterations. Analysis of enhancer profiles revealed a distinctive Mes-GC epigenomic landscape, with TEAD1 as a master regulator of Mes-GC enhancers and Mes-GCs exhibiting preferential sensitivity to TEAD1 pharmacological inhibition. Analysis of Mes-GC super-enhancers also highlighted NUAK1 kinase as a downstream target, with synergistic effects observed between NUAK1 inhibition and cisplatin treatment.

CONCLUSION

Our results establish a consensus Mes-GC classifier applicable to multiple transcriptomic scenarios. Mes-GCs exhibit a distinct epigenomic landscape, and TEAD1 inhibition and combinatorial NUAK1 inhibition/cisplatin may represent potential targetable options.

Smarcd3 is an epigenetic modulator of the metabolic landscape in pancreatic ductal adenocarcinoma

Pancreatic cancer is characterized by extensive resistance to conventional therapies, making clinical management a challenge. Here we map the epigenetic dependencies of cancer stem cells, cells that preferentially evade therapy and drive progression, and identify SWI/SNF complex member SMARCD3 as a regulator of pancreatic cancer cells. Although SWI/SNF subunits often act as tumor suppressors, we show that SMARCD3 is amplified in cancer, enriched in pancreatic cancer stem cells and upregulated in the human disease. Diverse genetic mouse models of pancreatic cancer and stage-specific Smarcd3 deletion reveal that Smarcd3 loss preferentially impacts established tumors, improving survival especially in context of chemotherapy. Mechanistically, SMARCD3 acts with FOXA1 to control lipid and fatty acid metabolism, programs associated with therapy resistance and poor prognosis in cancer. These data identify SMARCD3 as an epigenetic modulator responsible for establishing the metabolic landscape in aggressive pancreatic cancer cells and a potential target for new therapies.

Role of epigenetics in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive cancers, associated with poor survival outcomes. Lack of early diagnosis, resistance to conventional therapeutic treatments (including immunotherapy) and recurrence are some of the major hurdles in PDAC and contribute to its poor survival rate. While the risk of genetic predisposition to cancers is widely acknowledged and understood, recent advances in whole-genome and next-generation sequencing techniques have led to the acknowledgment of the role played by epigenetics, especially in PDAC. Epigenetic changes are heritable genetic modifications that influence gene expression without altering the DNA sequence. Epigenetic mechanisms (e.g., DNA methylation, post-translational modification of histone complexes and ncRNA) that result in reversible changes in gene expression are increasingly understood to be responsible for tumor initiation, development and even escape from immune surveillance. Our review seeks to highlight the various components of the epigenetic machinery that are known to be implicated in PDAC initiation and development and the feasibility of targeting these components to identify novel pharmacological strategies that could potentially lead to breakthroughs in PDAC treatment.

Interferon regulatory factor 1 (IRF1) controls the metabolic programmes of low-grade pancreatic cancer cells

OBJECTIVE

Pancreatic ductal adenocarcinomas (PDACs) include heterogeneous mixtures of low-grade cells forming pseudoglandular structures and compact nests of high-grade cells organised in non-glandular patterns. We previously reported that low-grade PDAC cells display high expression of interferon regulatory factor 1 (IRF1), a pivotal transcription factor of the interferon (IFN) system, suggesting grade-specific, cell-intrinsic activation of IFN responses. Here, we set out to determine the molecular bases and the functional impact of the activation of IFN-regulated responses in human PDACs.

DESIGN

We first confirmed the correlation between glandular differentiation and molecular subtypes of PDAC on the one hand, and the expression of IRF1 and IFN-stimulated genes on the other. We next used unbiased omics approaches to systematically analyse basal and IFN-regulated responses in low-grade and high-grade PDAC cells, as well as the impact of IRF1 on gene expression programmes and metabolic profiles of PDAC cells.

RESULTS

High-level expression of IRF1 in low-grade PDAC cells was controlled by endodermal lineage-determining transcription factors. IRF1-regulated gene expression equipped low-grade PDAC cells with distinctive properties related to antigen presentation and processing as well as responsiveness to IFN stimulation. Notably, IRF1 also controlled the characteristic metabolic profile of low-grade PDAC cells, suppressing both mitochondrial respiration and fatty acid synthesis, which may in part explain its growth-inhibiting activity.

CONCLUSION

IRF1 links endodermal differentiation to the expression of genes controlling antigen presentation and processing as well as to the specification of the metabolic profile characteristic of classical PDAC cells.

Dysregulation of HNF1B/Clusterin axis enhances disease progression in a highly aggressive subset of pancreatic cancer patients

Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy and is largely refractory to available treatments. Identifying key pathways associated with disease aggressiveness and therapeutic resistance may characterize candidate targets to improve patient outcomes. We used a strategy of examining the tumors from a subset of PDAC patient cohorts with the worst survival to understand the underlying mechanisms of aggressive disease progression and to identify candidate molecular targets with potential therapeutic significance. Non-negative matrix factorization (NMF) clustering, using gene expression profile, revealed three patient subsets. A 142-gene signature specific to the subset with the worst patient survival, predicted prognosis and stratified patients with significantly different survival in the test and validation cohorts. Gene-network and pathway analysis of the 142-gene signature revealed dysregulation of Clusterin (CLU) in the most aggressive patient subset in our patient cohort. Hepatocyte nuclear factor 1 b (HNF1B) positively regulated CLU, and a lower expression of HNF1B and CLU was associated with poor patient survival. Mechanistic and functional analyses revealed that CLU inhibits proliferation, 3D spheroid growth, invasiveness and epithelial-to-mesenchymal transition (EMT) in pancreatic cancer cell lines. Mechanistically, CLU enhanced proteasomal degradation of EMT-regulator, ZEB1. In addition, orthotopic transplant of CLU-expressing pancreatic cancer cells reduced tumor growth in mice. Furthermore, CLU enhanced sensitivity of pancreatic cancer cells representing aggressive patient subset, to the chemotherapeutic drug gemcitabine. Taken together, HNF1B/CLU axis negatively regulates pancreatic cancer progression and may potentially be useful in designing novel strategies to attenuate disease progression in PDAC patients.

HNF1B-driven three-dimensional chromatin structure for molecular classification in pancreatic cancers

The molecular subtypes of pancreatic cancer (PC), either classical/progenitor-like or basal/squamous-like, are currently a major topic of research because of their direct association with clinical outcomes. Some transcription factors (TFs) have been reported to be associated with these subtypes. However, the mechanisms by which these molecular signatures of PCs are established remain unknown. Epigenetic regulatory processes, supported by dynamic changes in the chromatin structure, are essential for transcriptional profiles. Previously, we reported the importance of open chromatin profiles in the biological features and transcriptional status of PCs. Here, we aimed to analyze the relationships between three-dimensional (3D) genome structures and the molecular subtypes of human PCs using Hi-C analysis. We observed a correlation of the specific elements of 3D genome modules, including compartments, topologically associating domains, and enhancer-promoter loops, with the expression of related genes. We focused on HNF1B, a TF that is implicated in the progenitor subtype. Forced expression of HNF1B in squamous-type PC organoids induced the upregulation and downregulation of genes associated with progenitor and squamous subtypes, respectively. Long-range genomic interactions induced by HNF1B were accompanied by compartment modulation and H3K27ac redistribution. We also found that these HNF1B-induced changes in subtype-related gene expression required an intrinsically disordered region, suggesting a possible involvement of phase separation in compartment modulation. Thus, mapping of 3D structural changes induced by TFs, such as HNF1B, may become a useful resource for further understanding the molecular features of PCs.

ACAP1 Deficiency Predicts Inferior Immunotherapy Response in Solid Tumors

BACKGROUND

ACAP1 plays a key role in endocytic recycling, which is essential for the normal function of lymphocytes. However, the expression and function of ACAP1 in lymphocytes have rarely been studied.

METHODS

Large-scale genomic data, including multiple bulk RNA-sequencing datasets, single-cell sequencing datasets, and immunotherapy cohorts, were exploited to comprehensively characterize ACAP1 expression, regulation, and function. Gene set enrichment analysis (GSEA) was used to uncover the pathways associated with ACAP1 expression. Eight algorithms, including TIMER, CIBERSORT, CIBERSORT-ABS, QUANTISEQ, xCELL, MCPCOUNTER, EPIC, and TIDE, were applied to estimate the infiltrating level of immune cells. Western blotting, qPCR, and ChIP-PCR were used to validate the findings from bioinformatic analyses. A T-cell co-culture killing assay was used to investigate the function of ACAP1 in lymphocytes.

RESULTS

ACAP1 was highly expressed in immune-related tissues and cells and minimally in other tissues. Moreover, single-cell sequencing analysis in tumor samples revealed that ACAP1 is expressed primarily in tumor-infiltrating lymphocytes (TILs), including T, B, and NK cells. ACAP1 expression is negatively regulated by promoter DNA methylation, with its promoter hypo-methylated in immune cells but hyper-methylated in other cells. Furthermore, SPI1 binds to the ACAP1 promoter and positively regulates its expression in immune cells. ACAP1 levels positively correlate with the infiltrating levels of TILs, especially CD8+ T cells, across a broad range of solid cancer types. ACAP1 deficiency is associated with poor prognosis and immunotherapeutic response in multiple cancer types treated with checkpoint blockade therapy (ICT). Functionally, the depletion of ACAP1 by RNA interference significantly impairs the T cell-mediated killing of tumor cells.

CONCLUSIONS

Our study demonstrates that ACAP1 is essential for the normal function of TILs, and its deficiency indicates an immunologically "cold" status of tumors that are resistant to ICT.

Structural variants drive context-dependent oncogene activation in cancer

Higher-order chromatin structure is important for the regulation of genes by distal regulatory sequences^1,2. Structural variants (SVs) that alter three-dimensional (3D) genome organization can lead to enhancer-promoter rewiring and human disease, particularly in the context of cancer³. However, only a small minority of SVs are associated with altered gene expression^4,5, and it remains unclear why certain SVs lead to changes in distal gene expression and others do not. To address these questions, we used a combination of genomic profiling and genome engineering to identify sites of recurrent changes in 3D genome structure in cancer and determine the effects of specific rearrangements on oncogene activation. By analysing Hi-C data from 92 cancer cell lines and patient samples, we identified loci affected by recurrent alterations to 3D genome structure, including oncogenes such as MYC, TERT and CCND1. By using CRISPR-Cas9 genome engineering to generate de novo SVs, we show that oncogene activity can be predicted by using 'activity-by-contact' models that consider partner region chromatin contacts and enhancer activity. However, activity-by-contact models are only predictive of specific subsets of genes in the genome, suggesting that different classes of genes engage in distinct modes of regulation by distal regulatory elements. These results indicate that SVs that alter 3D genome organization are widespread in cancer genomes and begin to illustrate predictive rules for the consequences of SVs on oncogene activation.

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.

Efficient Selection of Enhancers and Promoters from MIA PaCa-2 Pancreatic Cancer Cells by ChIP-lentiMPRA

A library of active genome regulatory elements (putative promoters and enhancers) from MIA PaCa-2 pancreatic adenocarcinoma cells was constructed using a specially designed lentiviral vector and a massive parallel reporter assay (ChIP-lentiMPRA). Chromatin immunoprecipitation of the cell genomic DNA by H3K27ac antibodies was used for primary enrichment of the library for regulatory elements. Totally, 11,264 unique genome regions, many of which are capable of enhancing the expression of the CopGFP reporter gene from the minimal CMV promoter, were identified. The regions tend to be located near promoters. Based on the proximity assay, we found an enrichment of highly expressed genes among those associated with three or more mapped distal regions (2 kb distant from the 5'-ends of genes). It was shown significant enrichment of genes related to carcinogenesis or Mia PaCa-2 cell identity genes in this group. In contrast, genes associated with 1-2 distal regions or only with proximal regions (within 2 kbp of the 5'-ends of genes) are more often related to housekeeping functions. Thus, ChIP-lentiMPRA is a useful strategy for creating libraries of regulatory elements for the study of tumor-specific gene transcription.

Extensive germline-somatic interplay contributes to prostate cancer progression through HNF1B co-option of TMPRSS2-ERG

Genome-wide association studies have identified 270 loci conferring risk for prostate cancer (PCa), yet the underlying biology and clinical impact remain to be investigated. Here we observe an enrichment of transcription factor genes including HNF1B within PCa risk-associated regions. While focused on the 17q12/HNF1B locus, we find a strong eQTL for HNF1B and multiple potential causal variants involved in the regulation of HNF1B expression in PCa. An unbiased genome-wide co-expression analysis reveals PCa-specific somatic TMPRSS2-ERG fusion as a transcriptional mediator of this locus and the HNF1B eQTL signal is ERG fusion status dependent. We investigate the role of HNF1B and find its involvement in several pathways related to cell cycle progression and PCa severity. Furthermore, HNF1B interacts with TMPRSS2-ERG to co-occupy large proportion of genomic regions with a remarkable enrichment of additional PCa risk alleles. We finally show that HNF1B co-opts ERG fusion to mediate mechanistic and biological effects of the PCa risk-associated locus 17p13.3/VPS53/FAM57A/GEMIN4. Taken together, we report an extensive germline-somatic interaction between TMPRSS2-ERG fusion and genetic variations underpinning PCa risk association and progression.

Axon guidance receptor ROBO3 modulates subtype identity and prognosis via AXL-associated inflammatory network in pancreatic cancer

Metastatic pancreatic cancer (PDAC) has a poor clinical outcome with a 5-year survival rate below 3%. Recent transcriptome profiling of PDAC biopsies has identified 2 clinically distinct subtypes - the "basal-like" (BL) subtype with poor prognosis and therapy resistance compared with the less aggressive and drug-susceptible "classical" (CLA) subtype. However, the mechanistic events and environmental factors that promote the BL subtype identity are not very clear. Using preclinical models, patient-derived xenografts, and FACS-sorted PDAC patient biopsies, we report here that the axon guidance receptor, roundabout guidance receptor 3 (ROBO3), promotes the BL metastatic program via a potentially unique AXL/IL-6/phosphorylated STAT3 (p-STAT3) regulatory axis. RNA-Seq identified a ROBO3-mediated BL-specific gene program, while tyrosine kinase profiling revealed AXL as the key mediator of the p-STAT3 activation. CRISPR/dCas9-based ROBO3 silencing disrupted the AXL/p-STAT3 signaling axis, thereby halting metastasis and enhancing therapy sensitivity. Transcriptome analysis of resected patient tumors revealed that AXLhi neoplastic cells associated with the inflammatory stromal program. Combining AXL inhibitor and chemotherapy substantially restored a CLA phenotypic state and reduced disease aggressiveness. Thus, we conclude that a ROBO3-driven hierarchical network determines the inflammatory and prometastatic programs in a specific PDAC subtype.

ID1 marks the tumorigenesis of pancreatic ductal adenocarcinoma in mouse and human

Pancreatic Ductal Adenocarcinoma (PDAC) is a deadly disease that has an increasing death rate but no effective treatment to now. Although biological and immunological hallmarks of PDAC have been frequently reported recently, early detection and the particularly aggressive biological features are the major challenges remaining unclear. In the current study, we retrieved multiple scRNA-seq datasets and illustrated the genetic programs of PDAC development in genetically modified mouse models. Notably, the transcription levels of Id1 were elevated specifically along with the PDAC development. Pseudotime trajectory analysis revealed that Id1 was closely correlated with the malignancy of PDAC. The gene expression patterns of human PDAC cells were determined by the comparative analysis of the scRNA-seq data on human PDAC and normal pancreas tissues. ID1 levels in human PDAC cancer cells were dramatically increased compared to normal epithelial cells. ID1 deficiency in vitro significantly blunt the invasive tumor-formation related phenotypes. IPA analysis on the differentially expressed genes suggested that EIF2 signaling was the core pathway regulating the development of PDAC. Blocking EFI2 signaling remarkably decreased the expression of ID1 and attenuated the tumor-formation related phenotypes. These observations confirmed that ID1 was regulated by EIF2 signaling and was the critical determinator of PDAC development and progression. This study suggests that ID1 is a potential malignant biomarker of PDAC in both mouse models and human and detecting and targeting ID1 may be a promising strategy to treat or even rescue PDAC.

MNX1-HNF1B Axis Is Indispensable for Intraductal Papillary Mucinous Neoplasm Lineages

BACKGROUND & AIMS

Chromatin architecture governs cell lineages by regulating the specific gene expression; however, its role in the diversity of cancer development remains unknown. Among pancreatic cancers, pancreatic ductal adenocarcinoma (PDAC) and intraductal papillary mucinous neoplasms (IPMN) with an associated invasive carcinoma (IPMNinv) arise from 2 distinct precursors, and their fundamental differences remain obscure. Here, we aimed to assess the difference of chromatin architecture regulating the transcriptional signatures or biological features in pancreatic cancers.

METHODS

We established 28 human organoids from distinct subtypes of pancreatic tumors, including IPMN, IPMNinv, and PDAC. We performed exome sequencing (seq), RNA-seq, assay for transposase-accessible chromatin-seq, chromatin immunoprecipitation-seq, high-throughput chromosome conformation capture, and phenotypic analyses with short hairpin RNA or clustered regularly interspaced short palindromic repeats interference.

RESULTS

Established organoids successfully reproduced the histology of primary tumors. IPMN and IPMNinv organoids harbored GNAS, RNF43, or KLF4 mutations and showed the distinct expression profiles compared with PDAC. Chromatin accessibility profiles revealed the gain of stomach-specific open regions in IPMN and the pattern of diverse gastrointestinal tissues in IPMNinv. In contrast, PDAC presented an impressive loss of accessible regions compared with normal pancreatic ducts. Transcription factor footprint analysis and functional assays identified that MNX1 and HNF1B were biologically indispensable for IPMN lineages. The upregulation of MNX1 was specifically marked in the human IPMN lineage tissues. The MNX1-HNF1B axis governed a set of genes, including MYC, SOX9, and OLFM4, which are known to be essential for gastrointestinal stem cells. High-throughput chromosome conformation capture analysis suggested the HNF1B target genes to be 3-dimensionally connected in the genome of IPMNinv.

CONCLUSIONS

Our organoid analyses identified the MNX1-HNF1B axis to be biologically significant in IPMN lineages.

Zinc Finger E-Box Binding Homeobox 2 as a Prognostic Biomarker in Various Cancers and Its Correlation with Infiltrating Immune Cells in Ovarian Cancer

This study investigated the expression of zinc finger E-box binding homeobox 2 (ZEB2), its prognostic significance in various cancers, and the correlation between ZEB2 and infiltrating immune cells and ZEB2-related proteins in ovarian cancer (OV). The Gene Expression Profiling Interactive Analysis tool was used to analyze RNA sequencing data and cancer survival rates, based on normal and tumor tissue data available in The Cancer Genome Atlas (TCGA) database. The Kaplan–Meier plotter and PrognoScan databases were used to analyze the prognostic value of ZEB2 in OV (n = 1144). The Tumor Immune Estimation Resource was used to investigate the correlation between ZEB2 and infiltrating immune cells in various cancers, including OV. High ZEB2 expression was associated with a poorer prognosis in OV. In OV, ZEB2 is positively correlated with CD8+T cells, neutrophils, macrophages, and dendritic cell invasion; and ZEB2 is negatively correlated with tumor-infiltrating B cells. The STRING database was used to investigate the correlations with ZEB2-related proteins. The results reveal that ZEB2 was positively correlated with SMAD1 and SMAD2 in OV. Our findings may serve as a potential prognostic biomarker, and provide novel insights into the tumor immunology in OV. Thus, ZEB2 may be a potential diagnostic and therapeutic target in OV.

Identification of the gene expression changes and gene regulatory aspects in ELF3 mutant bladder cancer

BACKGROUND

Recent genome-wide studies revealed the molecular subtypes and mutational landscape of bladder cancer, which is the 10th most common cancer causing many deaths. ELF3 is one of the frequently mutated genes in bladder cancer with 14% alteration rate. It mainly functions as an epithelial transcription factor and its proper function is critical for the urothelium development. However, the impact of ELF3 mutations in bladder cancer is currently unknown.

METHODS AND RESULTS

In this study, we analysed the gene expression data available for primary bladder cancer and bladder cancer cell lines according to the mutation status of ELF3. Our results show that de-regulated genes common in cell lines and primary tissue are primarily involved in ameboidal type cell migration and cell-cell junction organization. Additionally, we identify that ELF3-mutant cases in primary samples significantly overexpress PIK3C2B and ELF3 and PIK3C2B and ELF3 are significantly co-mutated in many cancer types. Our integrative analysis with existing Hi-C data further revealed the genes proximally located to ELF3, including PIK3C2B to be upregulated in ELF3 mutant cases, potentially as a result of truncated ELF3 protein product and subsequent changes in regulatory interactions.

CONCLUSIONS

Our results provide important insights about how ELF3 mutation contributes to bladder tumorigenesis and uncover previously unknown dependencies.

HDAC2 Facilitates Pancreatic Cancer Metastasis

The mortality of patients with pancreatic ductal adenocarcinoma (PDAC) is strongly associated with metastasis, a multistep process that is incompletely understood in this disease. Although genetic drivers of PDAC metastasis have not been defined, transcriptional and epigenetic rewiring can contribute to the metastatic process. The epigenetic eraser histone deacetylase 2 (HDAC2) has been connected to less differentiated PDAC, but the function of HDAC2 in PDAC has not been comprehensively evaluated. Using genetically defined models, we show that HDAC2 is a cellular fitness factor that controls cell cycle in vitro and metastasis in vivo, particularly in undifferentiated, mesenchymal PDAC cells. Unbiased expression profiling detected a core set of HDAC2-regulated genes. HDAC2 controlled expression of several prosurvival receptor tyrosine kinases connected to mesenchymal PDAC, including PDGFRα, PDGFRβ, and EGFR. The HDAC2-maintained program disabled the tumor-suppressive arm of the TGFβ pathway, explaining impaired metastasis formation of HDAC2-deficient PDAC. These data identify HDAC2 as a tractable player in the PDAC metastatic cascade. The complexity of the function of epigenetic regulators like HDAC2 implicates that an increased understanding of these proteins is needed for implementation of effective epigenetic therapies.

SIGNIFICANCE

HDAC2 has a context-specific role in undifferentiated PDAC and the capacity to disseminate systemically, implicating HDAC2 as targetable protein to prevent metastasis.

Loss of FGFR4 promotes the malignant phenotype of PDAC

Transcriptomic analyses of pancreatic ductal adenocarcinoma (PDAC) have identified two major epithelial subtypes with distinct biology and clinical behaviours. Here, we aimed to clarify the role of FGFR1 and FGFR4 in the definition of aggressive PDAC phenotypes. We found that the expression of FGFR4 is exclusively detected in epithelial cells, significantly elevated in the classical PDAC subtype, and associates with better outcomes. In highly aggressive basal-like/squamous PDAC, reduced FGFR4 expression aligns with hypermethylation of the gene and lower levels of histone marks associated with active transcription in its regulatory regions. Conversely, FGFR1 has more promiscuous expression in both normal and malignant pancreatic tissues and is strongly associated with the EMT phenotype but not with the basal-like cell lineage. Regardless of the genetic background, the increased proliferation of FGFR4-depleted PDAC cells correlates with hyperactivation of the mTORC1 pathway both in vitro and in vivo. Downregulation of FGFR4 in classical cell lines invariably leads to the enrichment of basal-like/squamous gene programs and is associated with either partial or full switch of phenotype. In sum, we show that endogenous levels of FGFR4 limit the malignant phenotype of PDAC cells. Finally, we propose FGFR4 as a valuable marker for the stratification of PDAC patients.

Deregulation of Transcription Factor Networks Driving Cell Plasticity and Metastasis in Pancreatic Cancer

Pancreatic cancer is a very aggressive disease with 5-year survival rates of less than 10%. The constantly increasing incidence and stagnant patient outcomes despite changes in treatment regimens emphasize the requirement of a better understanding of the disease mechanisms. Challenges in treating pancreatic cancer include diagnosis at already progressed disease states due to the lack of early detection methods, rapid acquisition of therapy resistance, and high metastatic competence. Pancreatic ductal adenocarcinoma, the most prevalent type of pancreatic cancer, frequently shows dominant-active mutations in KRAS and TP53 as well as inactivation of genes involved in differentiation and cell-cycle regulation (e.g. SMAD4 and CDKN2A). Besides somatic mutations, deregulated transcription factor activities strongly contribute to disease progression. Specifically, transcriptional regulatory networks essential for proper lineage specification and differentiation during pancreas development are reactivated or become deregulated in the context of cancer and exacerbate progression towards an aggressive phenotype. This review summarizes the recent literature on transcription factor networks and epigenetic gene regulation that play a crucial role during tumorigenesis.

Chromatin-Independent Interplay of NFATc1 and EZH2 in Pancreatic Cancer

Background: The Nuclear Factor of Activated T-cells 1 (NFATc1) transcription factor and the methyltransferase Enhancer of Zeste Homolog 2 (EZH2) significantly contribute to the aggressive phenotype of pancreatic ductal adenocarcinoma (PDAC). Herein, we aimed at dissecting the mechanistic background of their interplay in PDAC progression. Methods: NFATc1 and EZH2 mRNA and protein expression and complex formation were determined in transgenic PDAC models and human PDAC specimens. NFATc1 binding on the Ezh2 gene and the consequences of perturbed NFATc1 expression on Ezh2 transcription were explored by Chromatin Immunoprecipitation (ChIP) and upon transgenic or siRNA-mediated interference with NFATc1 expression, respectively. Integrative analyses of RNA- and ChIP-seq data was performed to explore NFATc1-/EZH2-dependent gene signatures. Results: NFATc1 targets the Ezh2 gene for transcriptional activation and biochemically interacts with the methyltransferase in murine and human PDAC. Surprisingly, our genome-wide binding and expression analyses do not link the protein complex to joint gene regulation. In contrast, our findings provide evidence for chromatin-independent functions of the NFATc1:EZH2 complex and reveal posttranslational EZH2 phosphorylation at serine 21 as a prerequisite for robust complex formation. Conclusion: Our findings disclose a previously unknown NFATc1-EZH2 axis operational in the pancreas and provide mechanistic insights into the conditions fostering NFATc1:EZH2 complex formation in PDAC.

HSP90 Inhibition Synergizes with Cisplatin to Eliminate Basal-like Pancreatic Ductal Adenocarcinoma Cells

Simple Summary

Pancreatic cancer is currently difficult to treat, but the drug cisplatin represents one of the most important therapeutic options. We find that cells derived from this cancer fall into two classes regarding their sensitivity towards cisplatin, and we observe that cells with high expression levels of GATA6 and microRNA 200 are mostly sensitive. However, those cells that respond poorly to cisplatin can be sensitized by drugs that inhibit HSP90, a protein that helps other proteins to fold properly. This was also found in a mouse model of pancreatic cancer. Our results suggest that the combination of cisplatin with HSP90-inhibitory drugs might improve the treatment of pancreatic cancer.

Abstract

To improve the treatment of pancreatic ductal adenocarcinoma (PDAC), a promising strategy consists of personalized chemotherapy based on gene expression profiles. Investigating a panel of PDAC-derived human cell lines, we found that their sensitivities towards cisplatin fall in two distinct classes. The platinum-sensitive class is characterized by the expression of GATA6, miRNA-200a, and miRNA-200b, which might be developable as predictive biomarkers. In the case of resistant PDAC cells, we identified a synergism of cisplatin with HSP90 inhibitors. Mechanistic explanations of this synergy include the degradation of Fanconi anemia pathway factors upon HSP90 inhibition. Treatment with the drug combination resulted in increased DNA damage and chromosome fragmentation, as we have reported previously for ovarian cancer cells. On top of this, HSP90 inhibition also enhanced the accumulation of DNA-bound platinum. We next investigated an orthotopic syngeneic animal model consisting of tumors arising from KPC cells (LSL-KrasG12D/+; LSL-Trp53R172H/+; Pdx-1-Cre, C57/BL6 genetic background). Here again, when treating established tumors, the combination of cisplatin with the HSP90 inhibitor onalespib was highly effective and almost completely prevented further tumor growth. We propose that the combination of platinum drugs and HSP90 inhibitors might be worth testing in the clinics for the treatment of cisplatin-resistant PDACs.

TNF-α-producing macrophages determine subtype identity and prognosis via AP1 enhancer reprogramming in pancreatic cancer

Large-scale genomic profiling of pancreatic cancer (PDAC) has revealed two distinct subtypes: 'classical' and 'basal-like'. Their variable coexistence within the stromal immune microenvironment is linked to differential prognosis; however, the extent to which these neoplastic subtypes shape the stromal immune landscape and impact clinical outcome remains unclear. By combining preclinical models, patient-derived xenografts, as well as FACS-sorted PDAC patient biopsies, we show that the basal-like neoplastic state is sustained via BRD4-mediated cJUN/AP1 expression, which induces CCL2 to recruit tumor necrosis factor (TNF)-α-secreting macrophages. TNF-α+ macrophages force classical neoplastic cells into an aggressive phenotypic state via lineage reprogramming. Integration of ATAC-, ChIP- and RNA-seq data revealed distinct JUNB/AP1 (classical) and cJUN/AP1 (basal-like)-driven regulation of PDAC subtype identity. Pharmacological inhibition of BRD4 led to suppression of the BRD4-cJUN-CCL2-TNF-α axis, restoration of classical subtype identity and a favorable prognosis. Hence, patient-tailored therapy for a cJUNhigh/TNF-αhigh subtype is paramount in overcoming highly inflamed and aggressive PDAC states.

Hepatocyte nuclear factor HNF1A is a potential regulator in shaping the super-enhancer landscape in colorectal cancer liver metastasis

Super-enhancers (SEs) play essential roles in colorectal cancer (CRC) progression. However, how the SE landscape is orchestrated by transcriptional regulators and evolves is not clear. Using de novo motif analysis, we show that the hepatocyte nuclear factor 1 (HNF1)-binding motif is enriched in SEs in cell lines derived from liver metastases, but not in those from primary tumors. This finding was further validated by extending the method to pancreatic cancer and a pair of isogenic CRC lines. Next, we revealed HNF1-alpha (HNF1A) was majorly expressed and upregulated in CRC liver metastatic cell lines. Clinically, HNF1A was remarkably upregulated in synchronous liver metastases as compared to localized tumors. Collectively, our study implicates HNF1A as a key regulator in shaping the SE landscape in CRC liver metastasis.

Epithelial memory of inflammation limits tissue damage while promoting pancreatic tumorigenesis

Inflammation is a major risk factor for pancreatic ductal adenocarcinoma (PDAC). When occurring in the context of pancreatitis, KRAS mutations accelerate tumor development in mouse models. We report that long after its complete resolution, a transient inflammatory event primes pancreatic epithelial cells to subsequent transformation by oncogenic KRAS. Upon recovery from acute inflammation, pancreatic epithelial cells display an enduring adaptive response associated with sustained transcriptional and epigenetic reprogramming. Such adaptation enables the reactivation of acinar-to-ductal metaplasia (ADM) upon subsequent inflammatory events, thereby limiting tissue damage through a rapid decrease of zymogen production. We propose that because activating mutations of KRAS maintain an irreversible ADM, they may be beneficial and under strong positive selection in the context of recurrent pancreatitis.

Krüppel-like factor 5 regulates wound repair and the innate immune response in human airway epithelial cells

A complex network of transcription factors regulates genes involved in establishing and maintaining key biological properties of the human airway epithelium. However, detailed knowledge of the contributing factors is incomplete. Here we characterize the role of Krüppel-like factor 5 (KLF5), in controlling essential pathways of epithelial cell identity and function in the human lung. RNA-seq following siRNA-mediated depletion of KLF5 in the Calu-3 lung epithelial cell line identified significant enrichment of genes encoding chemokines and cytokines involved in the proinflammatory response and also components of the junctional complexes mediating cell adhesion. To determine direct gene targets of KLF5, we defined the cistrome of KLF5 using ChIP-seq in both Calu-3 and 16HBE14o- lung epithelial cell lines. Occupancy site concordance analysis revealed that KLF5 colocalized with the active histone modification H3K27ac and also with binding sites for the transcription factor CCAAT enhancer-binding protein beta (C/EBPβ). Depletion of KLF5 increased both the expression and secretion of cytokines including IL-1β, a response that was enhanced following exposure to Pseudomonas aeruginosa lipopolysaccharide. Calu-3 cells exhibited faster rates of repair after KLF5 depletion compared with negative controls in wound scratch assays. Similarly, CRISPR-mediated KLF5-null 16HBE14o- cells also showed enhanced wound closure. These data reveal a pivotal role for KLF5 in coordinating epithelial functions relevant to human lung disease.

Tumor cell heterogeneity and its transcriptional bases in pancreatic cancer: a tale of two cell types and their many variants

Pancreatic ductal adenocarcinoma (PDAC), one of the most highly lethal tumors, is characterized by complex histology, with a massive fibrotic stroma in which both pseudo-glandular structures and compact nests of abnormally differentiated tumor cells are embedded, in different proportions and with different mutual relationships in space. This complexity and the heterogeneity of the tumor component have hindered the development of a broadly accepted, clinically actionable classification of PDACs, either on a morphological or a molecular basis. Here, we discuss evidence suggesting that such heterogeneity can to a large extent, albeit not exclusively, be traced back to two main classes of PDAC cells that commonly coexist in the same tumor: cells that maintained their ability to differentiate toward endodermal, mucin-producing epithelia and epithelial cells unable to form glandular structures and instead characterized by various levels of squamous differentiation and the expression of mesenchymal lineage genes. The underlying gene regulatory networks and how they are controlled by distinct transcription factors, as well as the practical implications of these two different populations of tumor cells, are discussed.

Drivers of Gene Expression Dysregulation in Pancreatic Cancer

Pancreatic ductal adenocarcinoma (PDAC) remains a devastating disease with a poor prognosis. The functional consequences of common genetic aberrations and their roles in treatment strategies have been extensively reviewed. In addition to these genomic aberrations, consideration of non-genetic drivers of altered oncogene expression is essential to account for the diversity in PDAC phenotypes. In this review we seek to assess our current understanding of mechanisms of gene expression dysregulation. We focus on four drivers of gene expression dysregulation, including mutations, transcription factors, epigenetic regulators, and RNA stability/isoform regulation, in the context of PDAC pathogenesis. Recent studies provide much-needed insight into the role of gene expression dysregulation in dissecting tumor heterogeneity and stratifying patients for the development of personalized treatment strategies.

Interpreting type 1 diabetes risk with genetics and single-cell epigenomics

Genetic risk variants that have been identified in genome-wide association studies of complex diseases are primarily non-coding1. Translating these risk variants into mechanistic insights requires detailed maps of gene regulation in disease-relevant cell types2. Here we combined two approaches: a genome-wide association study of type 1 diabetes (T1D) using 520,580 samples, and the identification of candidate cis-regulatory elements (cCREs) in pancreas and peripheral blood mononuclear cells using single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) of 131,554 nuclei. Risk variants for T1D were enriched in cCREs that were active in T cells and other cell types, including acinar and ductal cells of the exocrine pancreas. Risk variants at multiple T1D signals overlapped with exocrine-specific cCREs that were linked to genes with exocrine-specific expression. At the CFTR locus, the T1D risk variant rs7795896 mapped to a ductal-specific cCRE that regulated CFTR; the risk allele reduced transcription factor binding, enhancer activity and CFTR expression in ductal cells. These findings support a role for the exocrine pancreas in the pathogenesis of T1D and highlight the power of large-scale genome-wide association studies and single-cell epigenomics for understanding the cellular origins of complex disease.

C646 inhibits G2/M cell cycle-related proteins and potentiates anti-tumor effects in pancreatic cancer

The activity of histone acetyltransferases (HATs) plays a central role in an epigenetic modification in cooperation with HDACs (histone deacetyl transferases). It is likely that malfunction of this enzymatic machinery controlling epigenetic modification is relevant to carcinogenesis and tumor progression. However, in pancreatic cancer, the clinical relevance of HAT activity and histone acetylation has remained unclear. We identified that H3 acetylation was expressed in all pancreatic cancer patients, indicating that H3 acetylation may be essential in pancreatic cancer cells. We also found that the HAT inhibitor C646 augmented anti-tumor effects in vitro by inhibiting cell proliferation and cell cycle progression concomitantly with suppression of acetylated H3K9 and H3K27 expression. C646 or p300 and CBP (CREB-binding protein)-specific siRNA treatment inhibited the transcription of the G2/M cell cycle regulatory proteins cyclin B1 and CDK1 (cyclin-dependent kinase 1). C646 treatment also inhibited tumor growth in vivo in a xenograft mouse model. C646 could be an effective therapeutic agent for pancreatic cancer. The epigenetic status of pancreatic cancers based on their level of histone H3 acetylation may influence patient survival. Epigenetic stratification according to H3K27 acetylation could be useful for predicting disease prognosis as well as the therapeutic efficacy of C646 in pancreatic cancer.