TP63-Mediated Enhancer Reprogramming Drives the Squamous Subtype of Pancreatic Ductal Adenocarcinoma

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.

Patient-derived organoids of pancreatic ductal adenocarcinoma for subtype determination and clinical outcome prediction

BACKGROUND

Recently, two molecular subtypes of pancreatic ductal adenocarcinoma (PDAC) have been proposed: the "Classical" and "Basal-like" subtypes, with the former showing better clinical outcomes than the latter. However, the "molecular" classification has not been applied in real-world clinical practice. This study aimed to establish patient-derived organoids (PDOs) for PDAC and evaluate their application in subtype classification and clinical outcome prediction.

METHODS

We utilized tumor samples acquired through endoscopic ultrasound-guided fine-needle biopsy and established a PDO library for subsequent use in morphological assessments, RNA-seq analyses, and in vitro drug response assays. We also conducted a prospective clinical study to evaluate whether analysis using PDOs can predict treatment response and prognosis.

RESULTS

PDOs of PDAC were established at a high efficiency (> 70%) with at least 100,000 live cells. Morphologically, PDOs were classified as gland-like structures (GL type) and densely proliferating inside (DP type) less than 2 weeks after tissue sampling. RNA-seq analysis revealed that the "morphological" subtype (GL vs. DP) corresponded to the "molecular" subtype ("Classical" vs. "Basal-like"). The "morphological" classification predicted the clinical treatment response and prognosis; the median overall survival of patients with GL type was significantly longer than that with DP type (P < 0.005). The GL type showed a better response to gemcitabine than the DP type in vitro, whereas the drug response of the DP type was improved by the combination of ERK inhibitor and chloroquine.

CONCLUSIONS

PDAC PDOs help in subtype determination and clinical outcome prediction, thereby facilitating the bench-to-bedside precision medicine for PDAC.

Whole-genome bisulfite sequencing identifies stage- and subtype-specific DNA methylation signatures in pancreatic cancer

In pancreatic ductal adenocarcinoma (PDAC), no recurrent metastasis-specific mutation has been found, suggesting that epigenetic mechanisms, such as DNA methylation, are the major contributors of late-stage disease progression. Here, we performed the first whole-genome bisulfite sequencing (WGBS) on mouse and human PDAC organoid models to identify stage-specific and molecular subtype-specific DNA methylation signatures. With this approach, we identified thousands of differentially methylated regions (DMRs) that can distinguish between the stages and molecular subtypes of PDAC. Stage-specific DMRs are associated with genes related to nervous system development and cell-cell adhesions, and are enriched in promoters and bivalent enhancers. Subtype-specific DMRs showed hypermethylation of GATA6 foregut endoderm transcriptional networks in the squamous subtype and hypermethylation of EMT transcriptional networks in the progenitor subtype. These results indicate that aberrant DNA methylation contributes to both PDAC progression and subtype differentiation, resulting in significant and reoccurring DNA methylation patterns with diagnostic and prognostic potential.

Crosstalk between metabolic remodeling and epigenetic reprogramming: A new perspective on pancreatic cancer

Pancreatic cancer is a highly malignant solid tumor with a poor prognosis and a high mortality rate. Thus, exploring the mechanisms underlying the development and progression of pancreatic cancer is critical for identifying targets for diagnosis and treatment. Two important hallmarks of cancer-metabolic remodeling and epigenetic reprogramming-are interconnected and closely linked to regulate one another, creating a complex interaction landscape that is implicated in tumorigenesis, invasive metastasis, and immune escape. For example, metabolites can be involved in the regulation of epigenetic enzymes as substrates or cofactors, and alterations in epigenetic modifications can in turn regulate the expression of metabolic enzymes. The crosstalk between metabolic remodeling and epigenetic reprogramming in pancreatic cancer has gained considerable attention. Here, we review the emerging data with a focus on the reciprocal regulation of metabolic remodeling and epigenetic reprogramming. We aim to highlight how these mechanisms could be applied to develop better therapeutic strategies.

TP63-TRIM29 axis regulates enhancer methylation and chromosomal instability in prostate cancer

BACKGROUND

Prostate adenocarcinoma (PRAD) is the second leading cause of cancer-related deaths in men. High variability in DNA methylation and a high rate of large genomic rearrangements are often observed in PRAD.

RESULTS

To investigate the reasons for such high variance, we integrated DNA methylation, RNA-seq, and copy number alterations datasets from The Cancer Genome Atlas (TCGA), focusing on PRAD, and employed weighted gene co-expression network analysis (WGCNA). Our results show that only single cluster of co-expressed genes is associated with genomic and epigenomic instability. Within this cluster, TP63 and TRIM29 are key transcription regulators and are downregulated in PRAD. We discovered that TP63 regulates the level of enhancer methylation in prostate basal epithelial cells. TRIM29 forms a complex with TP63 and together regulates the expression of genes specific to the prostate basal epithelium. In addition, TRIM29 binds DNA repair proteins and prevents the formation of the TMPRSS2:ERG gene fusion typically observed in PRAD.

CONCLUSION

Our study demonstrates that TRIM29 and TP63 are important regulators in maintaining the identity of the basal epithelium under physiological conditions. Furthermore, we uncover the role of TRIM29 in PRAD development.

Engrailed-1 Promotes Pancreatic Cancer Metastasis

Engrailed-1 (EN1) is a critical homeodomain transcription factor (TF) required for neuronal survival, and EN1 expression has been shown to promote aggressive forms of triple negative breast cancer. Here, it is reported that EN1 is aberrantly expressed in a subset of pancreatic ductal adenocarcinoma (PDA) patients with poor outcomes. EN1 predominantly repressed its target genes through direct binding to gene enhancers and promoters, implicating roles in the activation of MAPK pathways and the acquisition of mesenchymal cell properties. Gain- and loss-of-function experiments demonstrated that EN1 promoted PDA transformation and metastasis in vitro and in vivo. The findings nominate the targeting of EN1 and downstream pathways in aggressive PDA.

Integrative analysis of spatial and single-cell transcriptome data from human pancreatic cancer reveals an intermediate cancer cell population associated with poor prognosis

BACKGROUND

Recent studies using single-cell transcriptomic analysis have reported several distinct clusters of neoplastic epithelial cells and cancer-associated fibroblasts in the pancreatic cancer tumor microenvironment. However, their molecular characteristics and biological significance have not been clearly elucidated due to intra- and inter-tumoral heterogeneity.

METHODS

We performed single-cell RNA sequencing using enriched non-immune cell populations from 17 pancreatic tumor tissues (16 pancreatic cancer and one high-grade dysplasia) and generated paired spatial transcriptomic data from seven patient samples.

RESULTS

We identified five distinct functional subclusters of pancreatic cancer cells and six distinct cancer-associated fibroblast subclusters. We deeply profiled their characteristics, and we found that these subclusters successfully deconvoluted most of the features suggested in bulk transcriptome analysis of pancreatic cancer. Among those subclusters, we identified a novel cancer cell subcluster, Ep_VGLL1, showing intermediate characteristics between the extremities of basal-like and classical dichotomy, despite its prognostic value. Molecular features of Ep_VGLL1 suggest its transitional properties between basal-like and classical subtypes, which is supported by spatial transcriptomic data.

CONCLUSIONS

This integrative analysis not only provides a comprehensive landscape of pancreatic cancer and fibroblast population, but also suggests a novel insight to the dynamic states of pancreatic cancer cells and unveils potential therapeutic targets.

Epigenetic therapeutic strategies in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal solid malignancies, characterized by its aggressiveness and metastatic potential, with a 5-year survival rate of only 8-11%. Despite significant improvements in PDAC treatment and management, therapeutic alternatives are still limited. One of the main reasons is its high degree of intra- and inter-individual tumor heterogeneity which is established and maintained through a complex network of transcription factors and epigenetic regulators. Epigenetic drugs, have shown promising preclinical results in PDAC and are currently being evaluated in clinical trials both for their ability to sensitize cancer cells to cytotoxic drugs and to counteract the immunosuppressive characteristic of PDAC tumor microenvironment. In this review, we discuss the current status of epigenetic treatment strategies to overcome molecular and cellular PDAC heterogeneity in order to improve response to therapy.

High-resolution and quantitative spatial analysis reveal intra-ductal phenotypic and functional diversification in pancreatic cancer

A 'classical' and a 'basal-like' subtype of pancreatic cancer have been reported, with differential expression of GATA6 and different dosages of mutant KRAS. We established in situ detection of KRAS point mutations and mRNA panels for the consensus subtypes aiming to project these findings to paraffin-embedded clinical tumour samples for spatial quantitative analysis. We unveiled that, next to inter-patient and intra-patient inter-ductal heterogeneity, intraductal spatial phenotypes exist with anti-correlating expression levels of GATA6 and KRASG12D . The basal-like mRNA panel better captured the basal-like cell states than widely used protein markers. The panels corroborated the co-existence of the classical and basal-like cell states in a single tumour duct with functional diversification, i.e. proliferation and epithelial-to-mesenchymal transition respectively. Mutant KRASG12D detection ascertained an epithelial origin of vimentin-positive cells in the tumour. Uneven spatial distribution of cancer-associated fibroblasts could recreate similar intra-organoid diversification. This extensive heterogeneity with functional cooperation of plastic tumour cells poses extra challenges to therapeutic approaches. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

The urothelial gene regulatory network: understanding biology to improve bladder cancer management

The urothelium is a stratified epithelium composed of basal cells, one or more layers of intermediate cells, and an upper layer of differentiated umbrella cells. Most bladder cancers (BLCA) are urothelial carcinomas. Loss of urothelial lineage fidelity results in altered differentiation, highlighted by the taxonomic classification into basal and luminal tumors. There is a need to better understand the urothelial transcriptional networks. To systematically identify transcription factors (TFs) relevant for urothelial identity, we defined highly expressed TFs in normal human bladder using RNA-Seq data and inferred their genomic binding using ATAC-Seq data. To focus on epithelial TFs, we analyzed RNA-Seq data from patient-derived organoids recapitulating features of basal/luminal tumors. We classified TFs as "luminal-enriched", "basal-enriched" or "common" according to expression in organoids. We validated our classification by differential gene expression analysis in Luminal Papillary vs. Basal/Squamous tumors. Genomic analyses revealed well-known TFs associated with luminal (e.g., PPARG, GATA3, FOXA1) and basal (e.g., TP63, TFAP2) phenotypes and novel candidates to play a role in urothelial differentiation or BLCA (e.g., MECOM, TBX3). We also identified TF families (e.g., KLFs, AP1, circadian clock, sex hormone receptors) for which there is suggestive evidence of their involvement in urothelial differentiation and/or BLCA. Genomic alterations in these TFs are associated with BLCA. We uncover a TF network involved in urothelial cell identity and BLCA. We identify novel candidate TFs involved in differentiation and cancer that provide opportunities for a better understanding of the underlying biology and therapeutic intervention.

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

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

Cell-in-Cell-Mediated Entosis Reveals a Progressive Mechanism in Pancreatic Cancer

BACKGROUND & AIMS

Pancreatic ductal adenocarcinoma (PDAC) is a deadly malignancy with high intratumoral heterogeneity. There is a lack of effective therapeutics for PDAC. Entosis, a form of nonapoptotic regulated cell death mediated by cell-in-cell structures (CICs), has been reported in multiple cancers. However, the role of entosis in PDAC progression remains unclear.

METHODS

CICs were evaluated using immunohistochemistry and immunofluorescence staining. The formation of CICs was induced by suspension culture. Through fluorescence-activated cell sorting and single-cell RNA sequencing, entosis-forming cells were collected and their differential gene expression was analyzed. Cell functional assays and mouse models were used to investigate malignant phenotypes. Clinical correlations between entosis and PDAC were established by retrospective analysis.

RESULTS

Entosis was associated with an unfavorable prognosis for patients with PDAC and was more prevalent in liver metastases than in primary tumors. The single-cell RNA sequencing results revealed that several oncogenes were up-regulated in entosis-forming cells compared with parental cells. These highly entotic cells demonstrated higher oncogenic characteristics in vitro and in vivo. NET1, neuroepithelial cell transforming gene 1, is an entosis-related gene that plays a pivotal role in PDAC progression and is correlated with poor outcomes.

CONCLUSIONS

Entosis is correlated with PDAC progression, especially in liver metastasis. NET1 is a newly validated entosis-related gene and a molecular marker of poor outcomes. PDAC cells generate a highly aggressive subpopulation marked by up-regulated NET1 via entosis, which may drive PDAC progression.

Epigenetic control of pancreatic cancer metastasis

Surgical resection, when combined with chemotherapy, has been shown to significantly improve the survival rate of patients with pancreatic ductal adenocarcinoma (PDAC). However, this treatment option is only feasible for a fraction of patients, as more than 50% of cases are diagnosed with metastasis. The multifaceted process of metastasis is still not fully understood, but recent data suggest that transcriptional and epigenetic plasticity play significant roles. Interfering with epigenetic reprogramming can potentially control the adaptive processes responsible for metastatic progression and therapy resistance, thereby enhancing treatment responses and preventing recurrence. This review will focus on the relevance of histone-modifying enzymes in pancreatic cancer, specifically on their impact on the metastatic cascade. Additionally, it will also provide a brief update on the current clinical developments in epigenetic therapies.

Dual RNA-seq to catalogue host and parasite gene expression changes associated with virulence of T. annulata-transformed bovine leukocytes: towards identification of attenuation biomarkers

The apicomplexan parasite Theileria annulata is transmitted by Hyalomma ticks and causes an acute lymphoproliferative disease that is invariably lethal in exotic cattle breeds. The unique ability of the schizont stage of T. annulata to transform infected leukocytes to a cancer-like phenotype and the simplicity of culturing and passaging T. annulata-transformed cells in vitro have been explored for live vaccine development by attenuating the transformed cells using lengthy serial propagation in vitro. The empirical in vivo evaluation of attenuation required for each batch of long-term cultured cells is a major constraint since it is resource intensive and raises ethical issues regarding animal welfare. As yet, the molecular mechanisms underlying attenuation are not well understood. Characteristic changes in gene expression brought about by attenuation are likely to aid in the identification of novel biomarkers for attenuation. We set out to undertake a comparative transcriptome analysis of attenuated (passage 296) and virulent (passage 26) bovine leukocytes infected with a Tunisian strain of T. annulata termed Beja. RNA-seq was used to analyse gene expression profiles and the relative expression levels of selected genes were verified by real-time quantitative PCR (RT-qPCR) analysis. Among the 3538 T. annulata genes analysed, 214 were significantly differentially expressed, of which 149 genes were up-regulated and 65 down-regulated. Functional annotation of differentially expressed T. annulata genes revealed four broad categories of metabolic pathways: carbon metabolism, oxidative phosphorylation, protein processing in the endoplasmic reticulum and biosynthesis of secondary metabolites. It is interesting to note that of the top 40 genes that showed altered expression, 13 were predicted to contain a signal peptide and/or at least one transmembrane domain, suggesting possible involvement in host-parasite interaction. Of the 16,514 bovine transcripts, 284 and 277 showed up-regulated and down-regulated expression, respectively. These were assigned to functional categories relevant to cell surface, tissue morphogenesis and regulation of cell adhesion, regulation of leucocyte, lymphocyte and cell activation. The genetic alterations acquired during attenuation that we have catalogued herein, as well as the accompanying in silico functional characterization, do not only improve understanding of the attenuation process, but can also be exploited by studies aimed at identifying attenuation biomarkers across different cell lines focusing on some host and parasite genes that have been highlighted in this study, such as bovine genes (CD69, ZNF618, LPAR3, and APOL3) and parasite genes such as TA03875.

Identification of a ΔNp63-Dependent Basal-Like A Subtype-Specific Transcribed Enhancer Program (B-STEP) in Aggressive Pancreatic Ductal Adenocarcinoma

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.

An integrated genomic approach identifies follistatin as a target of the p63-epidermal growth factor receptor oncogenic network in head and neck squamous cell carcinoma

Although numerous putative oncogenes have been associated with the etiology of head and neck squamous cell carcinoma (HNSCC), the mechanisms by which these oncogenes and their downstream targets mediate tumor progression have not been fully elucidated. We performed an integrative analysis to identify a crucial set of targets of the oncogenic transcription factor p63 that are common across multiple transcriptomic datasets obtained from HNSCC patients, and representative cell line models. Notably, our analysis revealed FST which encodes follistatin, a secreted glycoprotein that inhibits the transforming growth factor TGFβ/activin signaling pathways, to be a direct transcriptional target of p63. In addition, we found that FST expression is also driven by epidermal growth factor receptor EGFR signaling, thus mediating a functional link between the TGF-β and EGFR pathways. We show through loss- and gain-of-function studies that FST predominantly imparts a tumor-growth and migratory phenotype in HNSCC cells. Furthermore, analysis of single-cell RNA sequencing data from HNSCC patients unveiled cancer cells as the dominant source of FST within the tumor microenvironment and exposed a correlation between the expression of FST and its regulators with immune infiltrates. We propose FST as a prognostic biomarker for patient survival and a compelling candidate mediating the broad effects of p63 on the tumor and its associated microenvironment.

Post-operative mortality and recurrence patterns in pancreatic cancer according to KRAS mutation and CDKN2A, p53, and SMAD4 expression

Alterations in KRAS, CDKN2A (p16), TP53, and SMAD4 genes have been major drivers of pancreatic carcinogenesis. The clinical course of patients with pancreatic cancer in relation to these driver alterations has not been fully characterised in large populations. We hypothesised that pancreatic carcinomas with different combinations of KRAS mutation and aberrant expression of CDKN2A, p53, and SMAD4 might show distinctive recurrence patterns and post-operative survival outcomes. To test this hypothesis, we utilised a multi-institutional cohort of 1,146 resected pancreatic carcinomas and assessed KRAS mutations by droplet digital polymerase chain reaction and CDKN2A, p53, and SMAD4 expression by immunohistochemistry. Multivariable hazard ratios (HRs) and 95% confidence intervals (CIs) for disease-free survival (DFS) and overall survival (OS) were computed according to each molecular alteration and the number of altered genes using the Cox regression models. Multivariable competing risks regression analyses were conducted to assess the associations of the number of altered genes with specific patterns of recurrence. Loss of SMAD4 expression was associated with short DFS (multivariable HR, 1.24; 95% CI, 1.09-1.43) and OS times (multivariable HR, 1.27; 95% CI, 1.10-1.46). Compared to cases with 0-2 altered genes, cases with three and four altered genes had multivariable HRs for OS of 1.28 (95% CI, 1.09-1.51) and 1.47 (95% CI, 1.22-1.78), respectively (ptrend  < 0.001). Patients with an increasing number of altered genes were more likely to have short DFS time (ptrend  = 0.003) and to develop liver metastasis (ptrend  = 0.006) rather than recurrence at local or other distant sites. In conclusion, loss of SMAD4 expression and an increasing number of altered genes were associated with unfavourable outcomes in pancreatic cancer patients. This study suggests that the accumulation of the four major driver alterations can confer a high metastatic potential to the liver, thereby impairing post-operative survival among patients with pancreatic cancer.

The Transcriptional and Epigenetic Landscape of Cancer Cell Lineage Plasticity

Lineage plasticity, a process whereby cells change their phenotype to take on a different molecular and/or histologic identity, is a key driver of cancer progression and therapy resistance. Although underlying genetic changes within the tumor can enhance lineage plasticity, it is predominantly a dynamic process controlled by transcriptional and epigenetic dysregulation. This review explores the transcriptional and epigenetic regulators of lineage plasticity and their interplay with other features of malignancy, such as dysregulated metabolism, the tumor microenvironment, and immune evasion. We also discuss strategies for the detection and treatment of highly plastic tumors.

SIGNIFICANCE

Lineage plasticity is a hallmark of cancer and a critical facilitator of other oncogenic features such as metastasis, therapy resistance, dysregulated metabolism, and immune evasion. It is essential that the molecular mechanisms of lineage plasticity are elucidated to enable the development of strategies to effectively target this phenomenon. In this review, we describe key transcriptional and epigenetic regulators of cancer cell plasticity, in the process highlighting therapeutic approaches that may be harnessed for patient benefit.

Overcoming therapeutic resistance in pancreatic cancer: Emerging opportunities by targeting BRCAs and p53

Pancreatic cancer (PC) is characterized by (epi)genetic and microenvironmental alterations that negatively impact the treatment outcomes. New targeted therapies have been pursued to counteract the therapeutic resistance in PC. Aiming to seek for new therapeutic options for PC, several attempts have been undertaken to exploit BRCA1/2 and TP53 deficiencies as promising actionable targets. The elucidation of the pathogenesis of PC highlighted the high prevalence of p53 mutations and their connection with the aggressiveness and therapeutic resistance of PC. Additionally, PC is associated with dysfunctions in several DNA repair-related genes, including BRCA1/2, which sensitize tumours to DNA-damaging agents. In this context, poly(ADP-ribose) polymerase (PARP) inhibitors (PARPi) were approved for mutant BRCA1/2 PC patients. However, acquired drug resistance has become a major drawback of PARPi. This review emphasizes the importance of targeting defective BRCAs and p53 pathways for advancing personalized PC therapy, with particular focus on how this approach may provide an opportunity to tackle PC resistance.

Pancreatic cancer epigenetics: adaptive metabolism reprograms starving primary tumors for widespread metastatic outgrowth

Pancreatic cancer is a paradigm for adaptation to extreme stress. That is because genetic drivers are selected during tissue injury with epigenetic imprints encoding wound healing responses. Ironically, epigenetic memories of trauma that facilitate neoplasia can also recreate past stresses to restrain malignant progression through symbiotic tumor:stroma crosstalk. This is best exemplified by positive feedback between neoplastic chromatin outputs and fibroinflammatory stromal cues that encase malignant glands within a nutrient-deprived desmoplastic stroma. Because epigenetic imprints are chemically encoded by nutrient-derived metabolites bonded to chromatin, primary tumor metabolism adapts to preserve malignant epigenetic fidelity during starvation. Despite these adaptations, stromal stresses inevitably awaken primordial drives to seek more hospitable climates. The invasive migrations that ensue facilitate entry into the metastatic cascade. Metastatic routes present nutrient-replete reservoirs that accelerate malignant progression through adaptive metaboloepigenetics. This is best exemplified by positive feedback between biosynthetic enzymes and nutrient transporters that saturate malignant chromatin with pro-metastatic metabolite byproducts. Here we present a contemporary view of pancreatic cancer epigenetics: selection of neoplastic chromatin under fibroinflammatory pressures, preservation of malignant chromatin during starvation stresses, and saturation of metastatic chromatin by nutritional excesses that fuel lethal metastasis.

Interactive enhancer hubs (iHUBs) mediate transcriptional reprogramming and adaptive resistance in pancreatic cancer

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.

Integrated analysis revealed hypoxia signatures and LDHA related to tumor cell dedifferentiation and unfavorable prognosis in pancreatic adenocarcinoma: Hypoxia in PDAC

Pancreatic ductal adenocarcinoma (PDAC) is a highly heterogeneous cancer with limited understanding of its classification and tumor microenvironment. Here, by analyzing single-nucleus RNA sequencing of 43, 817 tumor cells from 15 PDAC tumors and non-tumor, we find that hypoxia signatures were heterogeneous across samples and were potential regulators for tumor progression and more aggressive phenotype. Hypoxia-high PDAC tends to present a basal/squamous-like phenotype and has significantly increased outgoing signaling, which enhances tumor cell stemness and promotes metastasis, angiogenesis, and fibroblast differentiation in PDAC. Hypoxia is related to an extracellular matrix enriched microenvironment, and increased possibility of TP53 mutation in PDAC. TP63 is a specific marker of squamous-like phenotype, and presents elevated transcriptome levels in most hypoxia PDAC tumors. In summary, our research highlights the potential linkage of hypoxia, tumor progression and genome alteration in PDAC, leading to further understand of the formation of inter-tumoral and intra-tumoral heterogenous in PDAC. Our study extends the understanding of the diversity and transition of tumor cells in PDAC, which provides insight into future PDAC management.

Convergence of YAP/TAZ, TEAD and TP63 activity is associated with bronchial premalignant severity and progression

BACKGROUND

Bronchial premalignant lesions (PMLs) are composed primarily of cells resembling basal epithelial cells of the airways, which through poorly understood mechanisms have the potential to progress to lung squamous cell carcinoma (LUSC). Despite ongoing efforts that have mapped gene expression and cell diversity across bronchial PML pathologies, signaling and transcriptional events driving malignancy are poorly understood. Evidence has suggested key roles for the Hippo pathway effectors YAP and TAZ and associated TEAD and TP63 transcription factor families in bronchial basal cell biology and LUSC. In this study we examine the functional association of YAP/TAZ, TEADs and TP63 in bronchial epithelial cells and PMLs.

METHODS

We performed RNA-seq in primary human bronchial epithelial cells following small interfering RNA (siRNA)-mediated depletion of YAP/TAZ, TEADs or TP63, and combined these data with ChIP-seq analysis of these factors. Directly activated or repressed genes were identified and overlapping genes were profiled across gene expression data obtained from progressive or regressive human PMLs and across lung single cell RNA-seq data sets.

RESULTS

Analysis of genes regulated by YAP/TAZ, TEADs, and TP63 in human bronchial epithelial cells revealed a converged transcriptional network that is strongly associated with the pathological progression of bronchial PMLs. Our observations suggest that YAP/TAZ-TEAD-TP63 associate to cooperatively promote basal epithelial cell proliferation and repress signals associated with interferon responses and immune cell communication. Directly repressed targets we identified include the MHC Class II transactivator CIITA, which is repressed in progressive PMLs and associates with adaptive immune responses in the lung. Our findings provide molecular insight into the control of gene expression events driving PML progression, including those contributing to immune evasion, offering potential new avenues for lung cancer interception.

CONCLUSIONS

Our study identifies important gene regulatory functions for YAP/TAZ-TEAD-TP63 in the early stages of lung cancer development, which notably includes immune-suppressive roles, and suggest that an assessment of the activity of this transcriptional complex may offer a means to identify immune evasive bronchial PMLs and serve as a potential therapeutic target.

KMT2D links TGF-β signaling to noncanonical activin pathway and regulates pancreatic cancer cell plasticity

Although KMT2D, also known as MLL2, is known to play an essential role in development, differentiation, and tumor suppression, its role in pancreatic cancer development is not well understood. Here, we discovered a novel signaling axis mediated by KMT2D, which links TGF-β to the activin A pathway. We found that TGF-β upregulates a microRNA, miR-147b, which in turn leads to post-transcriptional silencing of KMT2D. Loss of KMT2D induces the expression and secretion of activin A, which activates a noncanonical p38 MAPK-mediated pathway to modulate cancer cell plasticity, promote a mesenchymal phenotype, and enhance tumor invasion and metastasis in mice. We observed a decreased KMT2D expression in human primary and metastatic pancreatic cancer. Furthermore, inhibition or knockdown of activin A reversed the protumoral role of KMT2D loss. These findings support a tumor-suppressive role of KMT2D in pancreatic cancer and identify miR-147b and activin A as novel therapeutic targets.

p73 activates transcriptional signatures of basal lineage identity and ciliogenesis in pancreatic ductal adenocarcinoma

During the progression of pancreatic ductal adenocarcinoma (PDAC), tumor cells are known to acquire transcriptional and morphological properties of the basal (also known as squamous) epithelial lineage, which leads to more aggressive disease characteristics. Here, we show that a subset of basal-like PDAC tumors aberrantly express p73 (TA isoform), which is a known transcriptional activator of basal lineage identity, ciliogenesis, and tumor suppression in normal tissue development. Using gain- and loss- of function experiments, we show that p73 is necessary and sufficient to activate genes related to basal identity (e.g. KRT5), ciliogenesis (e.g. FOXJ1), and p53-like tumor suppression (e.g. CDKN1A) in human PDAC models. Owing to the paradoxical combination of oncogenic and tumor suppressive outputs of this transcription factor, we propose that PDAC cells express a low level of p73 that is optimal for promoting lineage plasticity without severe impairment of cell proliferation. Collectively, our study reinforces how PDAC cells exploit master regulators of the basal epithelial lineage during disease progression.

Pancreatic cancer: Advances and challenges

Pancreatic ductal adenocarcinoma (PDAC) remains one of the deadliest cancers. Significant efforts have largely defined major genetic factors driving PDAC pathogenesis and progression. Pancreatic tumors are characterized by a complex microenvironment that orchestrates metabolic alterations and supports a milieu of interactions among various cell types within this niche. In this review, we highlight the foundational studies that have driven our understanding of these processes. We further discuss the recent technological advances that continue to expand our understanding of PDAC complexity. We posit that the clinical translation of these research endeavors will enhance the currently dismal survival rate of this recalcitrant disease.

Ras-dependent activation of BMAL2 regulates hypoxic metabolism in pancreatic cancer

To identify novel drivers of malignancy in pancreatic ductal adenocarcinoma (PDAC), we employed regulatory network analysis, which calculates the activity of transcription factors and other regulatory proteins based on the integrated expression of their positive and negative target genes. We generated a regulatory network for the malignant epithelial cells of human PDAC using gene expression data from a set of 197 laser capture microdissected human PDAC samples and 45 low-grade precursors, for which we had matched histopathological, clinical, and epidemiological annotation. We then identified the most highly activated and repressed regulatory proteins (e.g. master regulators or MRs) associated with four malignancy phenotypes: precursors vs. PDAC (initiation), low-grade vs. high grade histopathology (progression), survival post resection, and association with KRAS activity. Integrating across these phenotypes, the top MR of PDAC malignancy was found to be BMAL2, a member of the PAS family of bHLH transcription factors. Although the canonical function of BMAL2 is linked to the circadian rhythm protein CLOCK, annotation of BMAL2 target genes highlighted a potential role in hypoxia response. We previously demonstrated that PDAC is hypovascularized and hypoperfused, and here show that PDAC from the genetically engineered KPC model exists in a state of extreme hypoxia, with a partial oxygen pressure of <1mmHg. Given the close homology of BMAL2 to HIF1β (ARNT) and its potential to heterodimerize with HIF1A and HIF2A, we investigated whether BMAL2 plays a role in the hypoxic response of PDAC. Indeed, BMAL2 controlled numerous hypoxia response genes and could be inhibited following treatment with multiple RAF, MEK, and ERK inhibitors, validating its association with RAS activity. Knockout of BMAL2 in four human PDAC cell lines led to defects in growth and invasion in the setting of hypoxia. Strikingly, BMAL2 null cells failed to induce glycolysis upon exposure to severe hypoxia and this was associated with a loss of expression of the glycolytic enzyme LDHA. Moreover, HIF1A was no longer stabilized under hypoxia in BMAL2 knockout cells. By contrast, HIF2A was hyper-stabilized under hypoxia, indicating a dysregulation of hypoxia metabolism in response to BMAL2 loss. We conclude that BMAL2 is a master regulator of hypoxic metabolism in PDAC, serving as a molecular switch between the disparate metabolic roles of HIF1A- and HIF2A-dependent hypoxia responses.

Computational evaluation of bioactive compounds from Viscum album (mistletoe) as inhibitors of p63 for pancreatic cancer treatment

Pancreatic ductal adenocarcinoma is an aggressive malignancy usually detectable at the advanced stage, with a 5-year survival rate of less than 8%. It has been reported that a gene called tumor-protein 63 (TP63) is expressed in an aggressive form of pancreatic cancer with a squamous signature. Thus, inhibiting the activity of p63 can be a means of treating and managing PDA. Different studies have shown that plant constituents are rich and can be a promising source for discovering drug candidates. The extract from mistletoe (Viscum album) is known to contain anticancer compounds; however, the specific molecular mechanism of the bioactive compounds is unknown. This study examines the pancreatic cancer therapeutic potential of the bioactive compounds in the flavonoid and phenolic acid constituents of mistletoe by adopting structural bioinformatics and advanced theoretical chemistry techniques via molecular docking, molecular dynamics simulation, molecular mechanics/generalized Born surface area (MM/GBSA) calculations, pharmacokinetic analysis, and density functional theory analysis. The six best compounds from the flavonoid constituent with the highest binding affinity ranging from -6.8 kcal/mol to -6.7 kcal/mol were selected with the control gemcitabine (-5.5 kcal/mol) for further computational analysis after molecular docking. Furthermore, MM/GBSA calculation showed the highest binding energy for the selected docked compounds, which validates their inhibitory potential. Hence, the molecular dynamics simulation, post-simulation analysis, pharmacokinetics model, and DFT results showed that mistletoe compounds are reliable due to their stable interaction with the target protein and drug-likeness properties.Communicated by Ramaswamy H. Sarma.

Spontaneously evolved progenitor niches escape Yap oncogene addiction in advanced pancreatic ductal adenocarcinomas

Lineage plasticity has been proposed as a major source of intratumoral heterogeneity and therapeutic resistance. Here, by employing an inducible genetic engineered mouse model, we illustrate that lineage plasticity enables advanced Pancreatic Ductal Adenocarcinoma (PDAC) tumors to develop spontaneous relapse following elimination of the central oncogenic driver - Yap. Transcriptomic and immunohistochemistry analysis of a large panel of PDAC tumors reveals that within high-grade tumors, small niches of PDAC cells gradually evolve to re-activate pluripotent transcription factors (PTFs), which lessen their dependency on Yap. Comprehensive Cut&Tag analysis demonstrate that although acquisition of PTF expression is coupled with the process of epithelial-to-mesenchymal transition (EMT), PTFs form a core transcriptional regulatory circuitry (CRC) with Jun to overcome Yap dependency, which is distinct from the classic TGFb-induced EMT-TF network. A chemical-genetic screen and follow-up functional studies establish Brd4 as an epigenetic gatekeeper for the PTF-Jun CRC, and strong synergy between BET and Yap inhibitors in blocking PDAC growth.

Poor-prognosis molecular subtypes in adenocarcinomas of pancreato-biliary and gynecological origin: A systematic review

Pancreato-biliary and gynecological adenocarcinomas need better tools to predict clinical outcome. Potential prognostic mesenchymal(-like) transcriptome-based subtypes have been identified in these cancers. In this systematic review, we include studies into molecular subtyping and summarize biological and clinical features of the subtypes within and across sites of origin, searching for suggestions to improve classification and prognostication. PubMed and Embase were searched for original research articles describing potential mesenchymal(-like) mRNA-based subtypes in pancreato-biliary or gynecological adenocarcinomas. Studies limited to supervised clustering were excluded. Fourty-four studies, discussing cholangiocarcinomas, gallbladder, ampullary, pancreatic, ovarian, and endometrial adenocarcinomas were included. There was overlap in molecular and clinical features in mesenchymal(-like) subtypes across all adenocarcinomas. Approaches including microdissection were more likely to identify prognosis-associating subtypes. To conclude, molecular subtypes in pancreato-biliary and gynecological adenocarcinomas share biological and clinical characteristics. Furthermore, separation of stromal and epithelial signals should be applied in future studies into biliary and gynecological adenocarcinomas.

MUC16 and TP53 family co-regulate tumor-stromal heterogeneity in pancreatic adenocarcinoma

MUC16/CA125 is one of the few oldest cancer biomarkers still used in current clinical practice. As mesothelium is an abundant source of MUC16 and a major contributor to stromal heterogeneity in PDAC, we investigated the regulation of MUC16 in tumor and stromal compartments individually. The trajectories constructed using the single-cell transcriptomes of stromal cells from KPC tumors demonstrated continuity in the trajectory path between MUC16-expressing mesothelial cells and other CAF subsets. Further, the tumor tissues of MUC16 whole-body knockout (KPCM) showed dysregulation in the markers of actomyosin assembly and fibroblast differentiation (iCAF and myCAF), indicating that MUC16 has an extra-tumoral role in controlling CAF differentiation. Although we found mesothelium-derivative stromal cells to be bystanders in normal pancreas, the proportion of these cells was higher in invasive PDAC, particularly in TP53 deficient tumors. Moreover, we also detail the regulation of MUC16, KRAS, and SOX9 by TP53 family members (TP53 and TP63) using multi-omics data from knockout models, PDAC cell lines, and human PDAC tissues.

Research advances and treatment perspectives of pancreatic adenosquamous carcinoma

BACKGROUND

As a malignant tumor, pancreatic cancer has an extremely low overall 5-year survival rate. Pancreatic adenosquamous carcinoma (PASC), a rare pancreatic malignancy, owns clinical presentation similar to pancreatic ductal adenocarcinoma (PDAC), which is the most prevalent pancreatic cancer subtype. PASC is generally defined as a pancreatic tumor consisting mainly of adenocarcinoma tissue and squamous carcinoma tissue. Compared with PDAC, PASC has a higher metastatic potential and worse prognosis, and lacks of effective treatment options to date. However, the pathogenesis and treatment of PASC are not yet clear and are accompanied with difficulties.

CONCLUSION

The present paper systematically summarizes the possible pathogenesis, diagnosis methods, and further suggests potential new treatment directions through reviewing research results of PASC, including the clinical manifestations, pathological manifestation, the original hypothesis of squamous carcinoma and the potential regulatory mechanism. In short, the present paper provides a systematic review of the research progress and new ideas for the development mechanism and treatment of PASC.

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.

A Systemic and Integrated Analysis of p63-Driven Regulatory Networks in Mouse Oral Squamous Cell Carcinoma

Oral squamous cell carcinoma (OSCC) is the most common malignancy of the oral cavity and is linked to tobacco exposure, alcohol consumption, and human papillomavirus infection. Despite therapeutic advances, a lack of molecular understanding of disease etiology, and delayed diagnoses continue to negatively affect survival. The identification of oncogenic drivers and prognostic biomarkers by leveraging bulk and single-cell RNA-sequencing datasets of OSCC can lead to more targeted therapies and improved patient outcomes. However, the generation, analysis, and continued utilization of additional genetic and genomic tools are warranted. Tobacco-induced OSCC can be modeled in mice via 4-nitroquinoline 1-oxide (4NQO), which generates a spectrum of neoplastic lesions mimicking human OSCC and upregulates the oncogenic master transcription factor p63. Here, we molecularly characterized established mouse 4NQO treatment-derived OSCC cell lines and utilized RNA and chromatin immunoprecipitation-sequencing to uncover the global p63 gene regulatory and signaling network. We integrated our p63 datasets with published bulk and single-cell RNA-sequencing of mouse 4NQO-treated tongue and esophageal tumors, respectively, to generate a p63-driven gene signature that sheds new light on the role of p63 in murine OSCC. Our analyses reveal known and novel players, such as COTL1, that are regulated by p63 and influence various oncogenic processes, including metastasis. The identification of new sets of potential biomarkers and pathways, some of which are functionally conserved in human OSCC and can prognosticate patient survival, offers new avenues for future mechanistic studies.

CXCR1/2 dual-inhibitor ladarixin reduces tumour burden and promotes immunotherapy response in pancreatic cancer

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal malignancy with few therapeutic options available. Despite immunotherapy has revolutionised cancer treatment, the results obtained in PDAC are still disappointing. Emerging evidence suggests that chemokines/CXCRs-axis plays a pivotal role in immune tumour microenvironment modulation, which may influence immunotherapy responsiveness. Here, we evaluated the effectiveness of CXCR1/2 inhibitor ladarixin, alone or in combination with anti-PD-1, against immunosuppression in PDAC.

METHODS

A set of preclinical models was obtained by engrafting mouse PDAC-derived cells into syngeneic immune-competent mice, as well as by orthotopically transplanting patient-derived PDAC tumour into human immune-system-reconstituted (HIR) mice (HuCD34-NSG-mice). Tumour-bearing mice were randomly assigned to receive vehicles, ladarixin, anti-PD-1 or drugs combination.

RESULTS

CXCR1/2 inhibition by ladarixin reverted in vitro tumour-mediated M2 macrophages polarisation and migration. Ladarixin as single agent reduced tumour burden in cancer-derived graft (CDG) models with high-immunogenic potential and increased the efficacy of ICI in non-immunogenic CDG-resistant models. In a HIR mouse model bearing the immunogenic subtype of human PDAC, ladarixin showed high efficacy increasing the antitumor effect of anti-PD-1.

CONCLUSION

Ladarixin in combination with anti-PD-1 might represent an extremely effective approach for the treatment of immunotherapy refractory PDAC, allowing pro-tumoral to immune-permissive microenvironment conversion.

Metastatic phenotype and immunosuppressive tumour microenvironment in pancreatic ductal adenocarcinoma: Key role of the urokinase plasminogen activator (PLAU)

Background

Previous studies have revealed the role of dysregulated urokinase plasminogen activator (encoded by PLAU) expression and activity in several pathways associated with cancer progression. However, systematic investigation into the association of PLAU expression with factors that modulate PDAC (pancreatic ductal adenocarcinoma) progression is lacking, such as those affecting stromal (pancreatic stellate cell, PSC)-cancer cell interactions, tumour immunity, PDAC subtypes and clinical outcomes from potential PLAU inhibition.

Methods

This study used an integrated bioinformatics approach to identify prognostic markers correlated with PLAU expression using different transcriptomics, proteomics, and clinical data sets. We then determined the association of dysregulated PLAU and correlated signatures with oncogenic pathways, metastatic phenotypes, stroma, immunosuppressive tumour microenvironment (TME) and clinical outcome. Finally, using an in vivo orthotopic model of pancreatic cancer, we confirmed the predicted effect of inhibiting PLAU on tumour growth and metastasis.

Results

Our analyses revealed that PLAU upregulation is not only associated with numerous other prognostic markers but also associated with the activation of various oncogenic signalling pathways, aggressive phenotypes relevant to PDAC growth and metastasis, such as proliferation, epithelial-mesenchymal transition (EMT), stemness, hypoxia, extracellular cell matrix (ECM) degradation, upregulation of stromal signatures, and immune suppression in the tumour microenvironment (TME). Moreover, the upregulation of PLAU was directly connected with signalling pathways known to mediate PSC-cancer cell interactions. Furthermore, PLAU upregulation was associated with the aggressive basal/squamous phenotype of PDAC and significantly reduced overall survival, indicating that this subset of patients may benefit from therapeutic interventions to inhibit PLAU activity. Our studies with a clinically relevant orthotopic pancreatic model showed that even short-term PLAU inhibition is sufficient to significantly halt tumour growth and, importantly, eliminate visible metastasis.

Conclusion

Elevated PLAU correlates with increased aggressive phenotypes, stromal score, and immune suppression in PDAC. PLAU upregulation is also closely associated with the basal subtype type of PDAC; patients with this subtype are at high risk of mortality from the disease and may benefit from therapeutic targeting of PLAU.

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.

Computational Modeling of TP63-TP53 Interaction and Rational Design of Inhibitors: Implications for Therapeutics

Tumor protein p63 (TP63) is a member of the TP53 protein family that are important for development and in tumor suppression. Unlike TP53, TP63 is rarely mutated in cancer, but instead different TP63 isoforms regulate its activity. TA isoforms (TAp63) act as tumor suppressors, whereas ΔN isoforms are strong drivers of squamous or squamous-like cancers. Many of these tumors become addicted to ΔN isoforms and removal of ΔN isoforms result in cancer cell death. Furthermore, some TP53 conformational mutants (TP53CM) gain the ability to interact with TAp63 isoforms and inhibit their antitumorigenic function, while indirectly promoting tumorigenic function of ΔN isoforms, but the exact mechanism of TP63-TP53CM interaction is unclear. The changes in the balance of TP63 isoform activity are crucial to understanding the transition between normal and tumor cells. Here, we modeled TP63-TP53CM complex using computational approaches. We then used our models to design peptides to disrupt the TP63-TP53CM interaction and restore antitumorigenic TAp63 function. In addition, we studied ΔN isoform oligomerization and designed peptides to inhibit its oligomerization and reduce their tumorigenic activity. We show that some of our peptides promoted cell death in a TP63 highly expressed cancer cell line, but not in a TP63 lowly expressed cancer cell line. Furthermore, we performed kinetic-binding assays to validate binding of our peptides to their targets. Our computational and experimental analyses present a detailed model for the TP63-TP53CM interaction and provide a framework for potential therapeutic peptides for the elimination of TP53CM cancer cells.

A functional analysis of 180 cancer cell lines reveals conserved intrinsic metabolic programs

Cancer cells reprogram their metabolism to support growth and invasion. While previous work has highlighted how single altered reactions and pathways can drive tumorigenesis, it remains unclear how individual changes propagate at the network level and eventually determine global metabolic activity. To characterize the metabolic lifestyle of cancer cells across pathways and genotypes, we profiled the intracellular metabolome of 180 pan-cancer cell lines grown in identical conditions. For each cell line, we estimated activity for 49 pathways spanning the entirety of the metabolic network. Upon clustering, we discovered a convergence into only two major metabolic types. These were functionally confirmed by ¹³ C-flux analysis, lipidomics, and analysis of sensitivity to perturbations. They revealed that the major differences in cancers are associated with lipid, TCA cycle, and carbohydrate metabolism. Thorough integration of these types with multiomics highlighted little association with genetic alterations but a strong association with markers of epithelial-mesenchymal transition. Our analysis indicates that in absence of variations imposed by the microenvironment, cancer cells adopt distinct metabolic programs which serve as vulnerabilities for therapy.

Signaling pathways and targeted therapies in lung squamous cell carcinoma: mechanisms and clinical trials

Lung cancer is the leading cause of cancer-related death across the world. Unlike lung adenocarcinoma, patients with lung squamous cell carcinoma (LSCC) have not benefitted from targeted therapies. Although immunotherapy has significantly improved cancer patients' outcomes, the relatively low response rate and severe adverse events hinder the clinical application of this promising treatment in LSCC. Therefore, it is of vital importance to have a better understanding of the mechanisms underlying the pathogenesis of LSCC as well as the inner connection among different signaling pathways, which will surely provide opportunities for more effective therapeutic interventions for LSCC. In this review, new insights were given about classical signaling pathways which have been proved in other cancer types but not in LSCC, including PI3K signaling pathway, VEGF/VEGFR signaling, and CDK4/6 pathway. Other signaling pathways which may have therapeutic potentials in LSCC were also discussed, including the FGFR1 pathway, EGFR pathway, and KEAP1/NRF2 pathway. Next, chromosome 3q, which harbors two key squamous differentiation markers SOX2 and TP63 is discussed as well as its related potential therapeutic targets. We also provided some progress of LSCC in epigenetic therapies and immune checkpoints blockade (ICB) therapies. Subsequently, we outlined some combination strategies of ICB therapies and other targeted therapies. Finally, prospects and challenges were given related to the exploration and application of novel therapeutic strategies for LSCC.

Regulatory chromatin rewiring promotes metabolic switching during adaptation to oncogenic receptor tyrosine kinase inhibition

Oesophageal adenocarcinoma (OAC) patients show poor survival rates and there are few targeted molecular therapies available. However, components of the receptor tyrosine kinase (RTK) driven pathways are commonly mutated in OAC, typified by high frequency amplifications of the RTK ERBB2. ERBB2 can be therapeutically targeted, but this has limited clinical benefit due to the acquisition of drug resistance. Here we examined how OAC cells adapt to ERBB2 inhibition as they transition to a drug resistant state. ERBB2 inhibition triggers widespread remodelling of the accessible chromatin landscape and the underlying gene regulatory networks. The transcriptional regulators HNF4A and PPARGC1A play a key role in this network rewiring. Initially, inhibition of cell cycle associated gene expression programmes is observed, with compensatory increases in the programmes driving changes in metabolic activity. Both PPARGC1A and HNF4A are required for the acquisition of resistance to ERBB2 inhibition and PPARGC1A is instrumental in promoting a switch to dependency on oxidative phosphorylation. Our work therefore reveals the molecular pathways that support the acquisition of a resistant state and points to potential new therapeutic strategies to combat cellular adaptation and ensuing drug resistance.

The Role of the Microbiome in Pancreatic Cancer

Simple Summary

Pancreatic cancer is deadly cancer characterized by dense stroma creating an immunosuppressive tumor microenvironment. Accumulating evidences indicate that the microbiome plays an important role in pancreatic cancer development and progression via the local and systemic inflammation and immune responses. The alteration of the microbiome modulates the tumor microenvironment and immune system in pancreatic cancer, which affects the efficacy of chemotherapies including immune-targeted therapies. Understanding the role of microbiome and underlying mechanisms may lead to novel biomarkers and therapeutic strategies for pancreatic cancer. This review summarizes the current evidence on the role of the microbiome in pancreatic cancer.

Abstract

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignancies, with little improvement in outcomes in recent decades, although the molecular and phenotypic characterization of PDAC has contributed to advances in tailored therapies. PDAC is characterized by dense stroma surrounding tumor cells, which limits the efficacy of treatment due to the creation of a physical barrier and immunosuppressive environment. Emerging evidence regarding the microbiome in PDAC implies its potential role in the initiation and progression of PDAC. However, the underlying mechanisms of how the microbiome affects the local tumor microenvironment (TME) as well as the systemic immune system have not been elucidated in PDAC. In addition, therapeutic strategies based on the microbiome have not been established. In this review, we summarize the current evidence regarding the role of the microbiome in the development of PDAC and discuss a possible role for the microbiome in the early detection of PDAC in relation to premalignant pancreatic diseases, such as chronic pancreatitis and intraductal papillary mucinous neoplasm (IPMN). In addition, we discuss the potential role of the microbiome in the treatment of PDAC, especially in immunotherapy, although the biomarkers used to predict the efficacy of immunotherapy in PDAC are still unknown. A comprehensive understanding of tumor-associated immune responses, including those involving the microbiome, holds promise for new treatments in PDAC.

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.

Multi-omics data integration and modeling unravels new mechanisms for pancreatic cancer and improves prognostic prediction

Pancreatic ductal adenocarcinoma (PDAC), has recently been found to be a heterogeneous disease, although the extension of its diversity remains to be fully understood. Here, we harmonize transcriptomic profiles derived from both PDAC epithelial and microenvironment cells to develop a Master Regulators (MR)-Gradient model that allows important inferences on transcriptional networks, epigenomic states, and metabolomics pathways that underlies this disease heterogeneity. This gradient model was generated by applying a blind source separation based on independent components analysis and robust principal component analyses (RPCA), following regulatory network inference. The result of these analyses reveals that PDAC prognosis strongly associates with the tumor epithelial cell phenotype and the immunological component. These studies were complemented by integration of methylome and metabolome datasets generated from patient-derived xenograft (PDX), together experimental measurements of metabolites, immunofluorescence microscopy, and western blot. At the metabolic level, PDAC favorable phenotype showed a positive correlation with enzymes implicated in complex lipid biosynthesis. In contrast, the unfavorable phenotype displayed an augmented OXPHOS independent metabolism centered on the Warburg effect and glutaminolysis. Epigenetically, we find that a global hypermethylation profile associates with the worst prognosis. Lastly, we report that, two antagonistic histone code writers, SUV39H1/SUV39H2 (H3K9Me3) and KAT2B (H3K9Ac) were identified key deregulated pathways in PDAC. Our analysis suggests that the PDAC phenotype, as it relates to prognosis, is determined by a complex interaction of transcriptomic, epigenomic, and metabolic features. Furthermore, we demonstrated that PDAC prognosis could be modulated through epigenetics.

EZH2 regulates a SETDB1/ΔNp63α axis via RUNX3 to drive a cancer stem cell phenotype in squamous cell carcinoma

Enhancer of zeste homolog 2 (EZH2) and SET domain bifurcated 1 (SETDB1, also known as ESET) are oncogenic methyltransferases implicated in a number of human cancers. These enzymes typically function as epigenetic repressors of target genes by methylating histone H3 K27 and H3-K9 residues, respectively. Here, we show that EZH2 and SETDB1 are essential to proliferation in 3 SCC cell lines, HSC-5, FaDu, and Cal33. Additionally, we find both of these proteins highly expressed in an aggressive stem-like SCC sub-population. Depletion of either EZH2 or SETDB1 disrupts these stem-like cells and their associated phenotypes of spheroid formation, invasion, and tumor growth. We show that SETDB1 regulates this SCC stem cell phenotype through cooperation with ΔNp63α, an oncogenic isoform of the p53-related transcription factor p63. Furthermore, EZH2 is upstream of both SETDB1 and ΔNp63α, activating these targets via repression of the tumor suppressor RUNX3. We show that targeting this pathway with inhibitors of EZH2 results in activation of RUNX3 and repression of both SETDB1 and ΔNp63α, antagonizing the SCC cancer stem cell phenotype. This work highlights a novel pathway that drives an aggressive cancer stem cell phenotype and demonstrates a means of pharmacological intervention.

A p300/GATA6 axis determines differentiation and Wnt dependency in pancreatic cancer models

Wnt signaling regulates the balance between stemness and differentiation in multiple tissues and in cancer. RNF43-mutant pancreatic cancers are dependent on Wnt production, and pharmacologic blockade of the pathway, e.g., by PORCN inhibitors, leads to tumor differentiation. However, primary resistance to these inhibitors has been observed. To elucidate potential mechanisms, we performed in vivo CRISPR screens in PORCN inhibitor-sensitive RNF43-mutant pancreatic cancer xenografts. As expected, genes in the Wnt pathway whose loss conferred drug resistance were identified, including APC, AXIN1, and CTNNBIP1. Unexpectedly, the screen also identified the histone acetyltransferase EP300 (p300), but not its paralog, CREBBP (CBP). We found that EP300 is silenced due to genetic alterations in all the existing RNF43-mutant pancreatic cancer cell lines that are resistant to PORCN inhibitors. Mechanistically, loss of EP300 directly downregulated GATA6 expression, thereby silencing the GATA6-regulated differentiation program and leading to a phenotypic transition from the classical subtype to the dedifferentiated basal-like/squamous subtype of pancreatic cancer. EP300 mutation and loss of GATA6 function bypassed the antidifferentiation activity of Wnt signaling, rendering these cancer cells resistant to Wnt inhibition.

Understanding Tricky Cellular and Molecular Interactions in Pancreatic Tumor Microenvironment: New Food for Thought

Pancreatic ductal adenocarcinoma (PDAC) represents 90% of all pancreatic cancer cases and shows a high mortality rate among all solid tumors. PDAC is often associated with poor prognosis, due to the late diagnosis that leads to metastasis development, and limited efficacy of available treatments. The tumor microenvironment (TME) represents a reliable source of novel targets for therapy, and even if many of the biological interactions among stromal, immune, and cancer cells that populate the TME have been studied, much more needs to be clarified. The great limitation in the efficacy of current standard chemoterapy is due to both the dense fibrotic inaccessible TME barrier surrounding cancer cells and the immunological evolution from a tumor-suppressor to an immunosuppressive environment. Nevertheless, combinatorial therapies may prove more effective at overcoming resistance mechanisms and achieving tumor cell killing. To achieve this result, a deeper understanding of the pathological mechanisms driving tumor progression and immune escape is required in order to design rationale-based therapeutic strategies. This review aims to summarize the present knowledge about cellular interactions in the TME, with much attention on immunosuppressive functioning and a specific focus on extracellular matrix (ECM) contribution.

A GATA6-centred gene regulatory network involving HNFs and ΔNp63 controls plasticity and immune escape in pancreatic cancer

OBJECTIVE

Molecular taxonomy of tumours is the foundation of personalised medicine and is becoming of paramount importance for therapeutic purposes. Four transcriptomics-based classification systems of pancreatic ductal adenocarcinoma (PDAC) exist, which consistently identified a subtype of highly aggressive PDACs with basal-like features, including ΔNp63 expression and loss of the epithelial master regulator GATA6. We investigated the precise molecular events driving PDAC progression and the emergence of the basal programme.

DESIGN

We combined the analysis of patient-derived transcriptomics datasets and tissue samples with mechanistic experiments using a novel dual-recombinase mouse model for Gata6 deletion at late stages of KRasG12D-driven pancreatic tumorigenesis (Gata6LateKO).

RESULTS

This comprehensive human-to-mouse approach showed that GATA6 loss is necessary, but not sufficient, for the expression of ΔNp63 and the basal programme in patients and in mice. The concomitant loss of HNF1A and HNF4A, likely through epigenetic silencing, is required for the full phenotype switch. Moreover, Gata6 deletion in mice dramatically increased the metastatic rate, with a propensity for lung metastases. Through RNA-Seq analysis of primary cells isolated from mouse tumours, we show that Gata6 inhibits tumour cell plasticity and immune evasion, consistent with patient-derived data, suggesting that GATA6 works as a barrier for acquiring the fully developed basal and metastatic phenotype.

CONCLUSIONS

Our work provides both a mechanistic molecular link between the basal phenotype and metastasis and a valuable preclinical tool to investigate the most aggressive subtype of PDAC. These data, therefore, are important for understanding the pathobiological features underlying the heterogeneity of pancreatic cancer in both mice and human.

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.

The biology of pancreatic cancer morphology

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal of all human malignancies. PDAC precursor lesions, invasive primary PDAC, and metastatic PDAC each display distinct morphologies that reflect unique biology. This 'biomorphology' is determined by a complex neoplastic history of clonal phylogenetic relationships, geographic locations, external environmental exposures, intrinsic metabolic demands, and tissue migration patterns. Understanding the biomorphological evolution of PDAC progression is not only of academic interest but also of great practical value. Applying this knowledge to surgical pathology practice facilitates the correct diagnosis on routine H&E stains without additional ancillary studies in most cases. Here I provide a concise overview of the entire biomorphological spectrum of PDAC progression beginning with initial neoplastic transformation and ending in terminal distant metastasis. Most biopsy and resection specimens are currently obtained prior to treatment. As such, our understanding of untreated PDAC biomorphology is mature. The biomorphology of treated PDAC is less defined but will assume greater importance as the frequency of neoadjuvant therapy increases. Although this overview is slanted towards pathology, it is written so that pathologists, clinicians, and scientists alike might find it instructive for their respective disciplines.