Cellular heterogeneity during mouse pancreatic ductal adenocarcinoma progression at single-cell resolution

Epithelial to Mesenchymal Transition as Mechanism of Progression of Pancreatic Cancer: From Mice to Men

Owed to its aggressive yet subtle nature, pancreatic cancer remains unnoticed till an advanced stage so that in most cases the diagnosis is made when the cancer has already spread to other organs with deadly efficiency. The progression from primary tumor to metastasis involves an intricate cascade of events comprising the pleiotropic process of epithelial to mesenchymal transition (EMT) facilitating cancer spread. The elucidation of this pivotal phenotypic change in cancer cell morphology, initially heretic, moved from basic studies dissecting the progression of pancreatic cancer in animal models to move towards human disease, although no clinical translation of the concept emerged yet. Despite this transition, a full-blown mesenchymal phenotype may not be accomplished; rather, the plasticity of the program and its dependency on heterotopic signals implies a series of fluctuating modifications of cancer cells encompassing mesenchymal and epithelial features. Despite the evidence supporting the activation of EMT and MET during cancer progression, our understanding of the relationship between tumor microenvironment and EMT is not yet mature for a clinical application. In this review, we attempt to resume the knowledge on EMT and pancreatic cancer, aiming to include the EMT among the hallmarks of cancer that could potentially modify our clinical thinking with the purpose of filling the gap between the results pursued in basic research by animal models and those achieved in translational research by surrogate biomarkers, as well as their application for prognostic and predictive purposes.

Targeting TNF-α-producing macrophages activates antitumor immunity in pancreatic cancer via IL-33 signaling

Pancreatic ductal adenocarcinoma (PDA) remains resistant to immune therapies, largely owing to robustly fibrotic and immunosuppressive tumor microenvironments. It has been postulated that excessive accumulation of immunosuppressive myeloid cells influences immunotherapy resistance, and recent studies targeting macrophages in combination with checkpoint blockade have demonstrated promising preclinical results. Yet our understanding of tumor-associated macrophage (TAM) function, complexity, and diversity in PDA remains limited. Our analysis reveals significant macrophage heterogeneity, with bone marrow-derived monocytes serving as the primary source for immunosuppressive TAMs. These cells also serve as a primary source of TNF-α, which suppresses expression of the alarmin IL-33 in carcinoma cells. Deletion of Ccr2 in genetically engineered mice decreased monocyte recruitment, resulting in profoundly decreased TNF-α and increased IL-33 expression, decreased metastasis, and increased survival. Moreover, intervention studies targeting CCR2 with a new orthosteric inhibitor (CCX598) rendered PDA susceptible to checkpoint blockade, resulting in reduced metastatic burden and increased survival. Our data indicate that this shift in antitumor immunity is influenced by increased levels of IL-33, which increases dendritic cell and cytotoxic T cell activity. These data demonstrate that interventions to disrupt infiltration of immunosuppressive macrophages, or their signaling, have the potential to overcome barriers to effective immunotherapeutics for PDA.

Orchestration of mesenchymal plasticity and immune evasiveness via rewiring of the metabolic program in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is the most fatal cancer in humans, due to its difficulty of early detection and its high metastatic ability. The occurrence of epithelial to mesenchymal transition in preinvasive pancreatic lesions has been implicated in the early dissemination, drug resistance, and cancer stemness of PDAC. PDAC cells also have a reprogrammed metabolism, regulated by driver mutation-mediated pathways, a desmoplastic tumor microenvironment (TME), and interactions with stromal cells, including pancreatic stellate cells, fibroblasts, endothelial cells, and immune cells. Such metabolic reprogramming and its functional metabolites lead to enhanced mesenchymal plasticity, and creates an acidic and immunosuppressive TME, resulting in the augmentation of protumor immunity via cancer-associated inflammation. In this review, we summarize our recent understanding of how PDAC cells acquire and augment mesenchymal features via metabolic and immunological changes during tumor progression, and how mesenchymal malignancies induce metabolic network rewiring and facilitate an immune evasive TME. In addition, we also present our recent findings on the interesting relevance of the small G protein ADP-ribosylation factor 6-based signaling pathway driven by KRAS/TP53 mutations, inflammatory amplification signals mediated by the proinflammatory cytokine interleukin 6 and RNA-binding protein ARID5A on PDAC metabolic reprogramming and immune evasion, and finally discuss potential therapeutic strategies for the quasi-mesenchymal subtype of PDAC.

Differential integrated stress response and asparagine production drive symbiosis and therapy resistance of pancreatic adenocarcinoma cells

The pancreatic tumor microenvironment drives deregulated nutrient availability. Accordingly, pancreatic cancer cells require metabolic adaptations to survive and proliferate. Pancreatic cancer subtypes have been characterized by transcriptional and functional differences, with subtypes reported to exist within the same tumor. However, it remains unclear if this diversity extends to metabolic programming. Here, using metabolomic profiling and functional interrogation of metabolic dependencies, we identify two distinct metabolic subclasses among neoplastic populations within individual human and mouse tumors. Furthermore, these populations are poised for metabolic cross-talk, and in examining this, we find an unexpected role for asparagine supporting proliferation during limited respiration. Constitutive GCN2 activation permits ATF4 signaling in one subtype, driving excess asparagine production. Asparagine release provides resistance during impaired respiration, enabling symbiosis. Functionally, availability of exogenous asparagine during limited respiration indirectly supports maintenance of aspartate pools, a rate-limiting biosynthetic precursor. Conversely, depletion of extracellular asparagine with PEG-asparaginase sensitizes tumors to mitochondrial targeting with phenformin.

BRCA mutational status shapes the stromal microenvironment of pancreatic cancer linking clusterin expression in cancer associated fibroblasts with HSF1 signaling

Tumors initiate by mutations in cancer cells, and progress through interactions of the cancer cells with non-malignant cells of the tumor microenvironment. Major players in the tumor microenvironment are cancer-associated fibroblasts (CAFs), which support tumor malignancy, and comprise up to 90% of the tumor mass in pancreatic cancer. CAFs are transcriptionally rewired by cancer cells. Whether this rewiring is differentially affected by different mutations in cancer cells is largely unknown. Here we address this question by dissecting the stromal landscape of BRCA-mutated and BRCA Wild-type pancreatic ductal adenocarcinoma. We comprehensively analyze pancreatic cancer samples from 42 patients, revealing different CAF subtype compositions in germline BRCA-mutated vs. BRCA Wild-type tumors. In particular, we detect an increase in a subset of immune-regulatory clusterin-positive CAFs in BRCA-mutated tumors. Using cancer organoids and mouse models we show that this process is mediated through activation of heat-shock factor 1, the transcriptional regulator of clusterin. Our findings unravel a dimension of stromal heterogeneity influenced by germline mutations in cancer cells, with direct implications for clinical research.

Exploring the Biology of Cancer-Associated Fibroblasts in Pancreatic Cancer

Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy characterised by a stubbornly low 5-year survival which is essentially unchanged in the past 5 decades. Despite recent advances in chemotherapy and surgical outcomes, progress continues to lag behind that of other cancers. The PDAC microenvironment is characterised by a dense, fibrotic stroma of which cancer-associated fibroblasts (CAFs) are key players. CAFs and fibrosis were initially thought to be uniformly tumour-promoting, however this doctrine is now being challenged by a wealth of evidence demonstrating CAF phenotypic and functional heterogeneity. Recent technological advances have allowed for the molecular profiling of the PDAC tumour microenvironment at exceptional detail, and these technologies are being leveraged at pace to improve our understanding of this previously elusive cell population. In this review we discuss CAF heterogeneity and recent developments in CAF biology. We explore the complex relationship between CAFs and other cell types within the PDAC microenvironment. We discuss the potential for therapeutic targeting of CAFs, and we finally provide an overview of future directions for the field and the possibility of improving outcomes for patients with this devastating disease.

Systematic Comparison of Pancreatic Ductal Adenocarcinoma Models Identifies a Conserved Highly Plastic Basal Cell State

Intratumoral heterogeneity and cellular plasticity have emerged as hallmarks of cancer, including pancreatic ductal adenocarcinoma (PDAC). As PDAC portends a dire prognosis, a better understanding of the mechanisms underpinning cellular diversity in PDAC is crucial. Here, we investigated the cellular heterogeneity of PDAC cancer cells across a range of in vitro and in vivo growth conditions using single-cell genomics. Heterogeneity contracted significantly in two-dimensional and three-dimensional cell culture models but was restored upon orthotopic transplantation. Orthotopic transplants reproducibly acquired cell states identified in autochthonous PDAC tumors, including a basal state exhibiting coexpression and coaccessibility of epithelial and mesenchymal genes. Lineage tracing combined with single-cell transcriptomics revealed that basal cells display high plasticity in situ. This work defines the impact of cellular growth conditions on phenotypic diversity and uncovers a highly plastic cell state with the capacity to facilitate state transitions and promote intratumoral heterogeneity in PDAC.

SIGNIFICANCE

This work provides important insights into how different model systems of pancreatic ductal adenocarcinoma mold the phenotypic space of cancer cells, highlighting the power of in vivo models.

The roles of intratumour heterogeneity in the biology and treatment of pancreatic ductal adenocarcinoma

Intratumour heterogeneity (ITH) has become an important focus of cancer research in recent years. ITH describes the cellular variation that enables tumour evolution, including tumour progression, metastasis and resistance to treatment. The selection and expansion of genetically distinct treatment-resistant cancer cell clones provides one explanation for treatment failure. However, tumour cell variation need not be genetically encoded. In pancreatic ductal adenocarcinoma (PDAC) in particular, the complex tumour microenvironment as well as crosstalk between tumour and stromal cells result in exceptionally variable tumour cell phenotypes that are also highly adaptable. In this review we discuss four different types of phenotypic heterogeneity within PDAC, from morphological to metabolic heterogeneity. We suggest that these different types of ITH are not independent, but, rather, can inform one another. Lastly, we highlight recent findings that suggest how therapeutic efforts may halt PDAC progression by constraining cellular heterogeneity.

Single Cell RNA-Seq Identifies Immune-Related Prognostic Model and Key Signature-SPP1 in Pancreatic Ductal Adenocarcinoma

There are no reliable biomarkers for early diagnosis or prognosis evaluation in pancreatic ductal adenocarcinoma (PDAC). Multiple scRNA-seq datasets for PDAC were retrieved from online databases and combined with scRNA-seq results from our previous study. The malignant ductal cells were identified through calculating copy number variation (CNV) scores. The robust markers of malignant ductal cells in PDAC were found. Five immune-related signatures, including SPP1, LINC00683, SNHG10, LINC00237, and CASC19, were used to develop a risk score formula to predict the overall survival of PDAC patients. We also constructed an easy-to-use nomogram, combining risk score, N stage, and margin status. The expression level of SPP1 was related to the prognosis and immune regulators. We found that SPP1 was mainly expressed in ductal cells and macrophages in PDAC. In conclusion, we constructed a promising prognostic model based on immune-related signatures for PDAC using scRNA-seq and TCGA_PAAD datasets.

Single-cell sequencing and establishment of an 8-gene prognostic model for pancreatic cancer patients

Background

Single-cell sequencing (SCS) technologies enable analysis of gene structure and expression data at single-cell resolution. However, SCS analysis in pancreatic cancer remains largely unexplored.

Methods

We downloaded pancreatic cancer SCS data from different databases and applied appropriate dimensionality reduction algorithms. We identified 10 cell types and subsequently screened differentially expressed marker genes of these 10 cell types using FindAllMarkers analysis. Also, we evaluated the tumor immune microenvironment based on ESTIMATE and MCP-counter. Statistical enrichment was evaluated using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis. We used all candidate gene sets in KEGG database to perform gene set enrichment analysis. We used LASSO regression to reduce the number of genes in the pancreatic risk model by R package glmnet, followed by rtPCR to validate the expression of the signature genes in different pancreatic cancer cell lines.

Results

We identified 15 cell subpopulations by dimension reduction and data clustering. We divided the 15 subpopulations into 10 distinct cell types based on marker gene expression. Then, we performed functional enrichment analysis for the 352 marker genes in pancreatic cancer cells. Based on RNA expression data and prognostic information from TCGA and GEO datasets, we identified 42 prognosis-related genes, including 5 protective genes and 37 high-risk genes, which we used to identified two molecular subtypes. C1 subtype was associated with a better prognosis, whereas C2 subtype was associated with a worse prognosis. Moreover, chemokine and chemokine receptor genes were differentially expressed between C1 and C2 subtypes. Functional and pathway enrichment uncovered functional differences between C1 and C2 subtype. We identified eight genes that could serve as potential biomarkers for prognosis prediction in pancreatic cancer patients. These genes were used to establish an 8-gene pancreatic cancer prognostic model.

Conclusions

We established an 8-gene pancreatic cancer prognostic model. This model can meaningfully predict prognosis and treatment response in pancreatic cancer patients.

Spatially restricted drivers and transitional cell populations cooperate with the microenvironment in untreated and chemo-resistant pancreatic cancer

Pancreatic ductal adenocarcinoma is a lethal disease with limited treatment options and poor survival. We studied 83 spatial samples from 31 patients (11 treatment-naïve and 20 treated) using single-cell/nucleus RNA sequencing, bulk-proteogenomics, spatial transcriptomics and cellular imaging. Subpopulations of tumor cells exhibited signatures of proliferation, KRAS signaling, cell stress and epithelial-to-mesenchymal transition. Mapping mutations and copy number events distinguished tumor populations from normal and transitional cells, including acinar-to-ductal metaplasia and pancreatic intraepithelial neoplasia. Pathology-assisted deconvolution of spatial transcriptomic data identified tumor and transitional subpopulations with distinct histological features. We showed coordinated expression of TIGIT in exhausted and regulatory T cells and Nectin in tumor cells. Chemo-resistant samples contain a threefold enrichment of inflammatory cancer-associated fibroblasts that upregulate metallothioneins. Our study reveals a deeper understanding of the intricate substructure of pancreatic ductal adenocarcinoma tumors that could help improve therapy for patients with this disease.

DEMOC: a deep embedded multi-omics learning approach for clustering single-cell CITE-seq data

Advances in single-cell RNA sequencing (scRNA-seq) technologies has provided an unprecedent opportunity for cell-type identification. As clustering is an effective strategy towards cell-type identification, various computational approaches have been proposed for clustering scRNA-seq data. Recently, with the emergence of cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq), the cell surface expression of specific proteins and the RNA expression on the same cell can be captured, which provides more comprehensive information for cell analysis. However, existing single cell clustering algorithms are mainly designed for single-omic data, and have difficulties in handling multi-omics data with diverse characteristics efficiently. In this study, we propose a novel deep embedded multi-omics clustering with collaborative training (DEMOC) model to perform joint clustering on CITE-seq data. Our model can take into account the characteristics of transcriptomic and proteomic data, and make use of the consistent and complementary information provided by different data sources effectively. Experiment results on two real CITE-seq datasets demonstrate that our DEMOC model not only outperforms state-of-the-art single-omic clustering methods, but also achieves better and more stable performance than existing multi-omics clustering methods. We also apply our model on three scRNA-seq datasets to assess the performance of our model in rare cell-type identification, novel cell-subtype detection and cellular heterogeneity analysis. Experiment results illustrate the effectiveness of our model in discovering the underlying patterns of data.

Heterogeneity of Cancer-Associated Fibroblasts and the Tumor Immune Microenvironment in Pancreatic Cancer

Simple Summary

Stroma-targeting therapy in pancreatic ductal adenocarcinoma (PDAC) has been extensively investigated, but no candidates have shown efficacy at the clinical trial stage. Studies of cancer-associated fibroblast (CAF) depletion in a mouse model suggested that CAFs have not only tumor-promoting function but also tumor-suppressive activity. Recently, single-cell RNA sequencing (scRNA-seq) has revealed the complex tumor microenvironment within PDAC, and subpopulations of functionally distinct CAFs and their association with tumor immunity have been reported. However, the existence of tumor suppressive CAFs and CAFs involved in the maintenance of PDAC differentiation has also been reported. In the future, therapeutic strategies should be developed considering these CAF subpopulations, with the hope of improving the prognosis of PDAC.

Abstract

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers, with a 5-year survival rate of 9%. Cancer-associated fibroblasts (CAFs) have historically been considered tumor-promoting. However, multiple studies reporting that suppression of CAFs in PDAC mouse models resulted in more aggressive tumors and worse prognosis have suggested the existence of a tumor-suppressive population within CAFs, leading to further research on heterogeneity within CAFs. In recent years, the benefits of cancer immunotherapy have been reported in various carcinomas. Unfortunately, the efficacy of immunotherapies in PDAC has been limited, and the CAF-driven cancer immunosuppressive microenvironment has been suggested as the cause. Thus, clarification of heterogeneity within the tumor microenvironment, including CAFs and tumor immunity, is urgently needed to establish effective therapeutic strategies for PDAC. In this review, we report the latest findings on the heterogeneity of CAFs and the functions of each major CAF subtype, which have been revealed by single-cell RNA sequencing in recent years. We also describe reports of tumor-suppressive CAF subtypes and the existence of CAFs that maintain a differentiated PDAC phenotype and review the potential for targeted therapy.

Metabolic reprogramming and crosstalk of cancer-related fibroblasts and immune cells in the tumor microenvironment

It is notorious that cancer cells alter their metabolism to adjust to harsh environments of hypoxia and nutritional starvation. Metabolic reprogramming most often occurs in the tumor microenvironment (TME). TME is defined as the cellular environment in which the tumor resides. This includes surrounding blood vessels, fibroblasts, immune cells, signaling molecules and the extracellular matrix (ECM). It is increasingly recognized that cancer cells, fibroblasts and immune cells within TME can regulate tumor progression through metabolic reprogramming. As the most significant proportion of cells among all the stromal cells that constitute TME, cancer-associated fibroblasts (CAFs) are closely associated with tumorigenesis and progression. Multitudinous studies have shown that CAFs participate in and promote tumor metabolic reprogramming and exert regulatory effects via the dysregulation of metabolic pathways. Previous studies have demonstrated that curbing the substance exchange between CAFs and tumor cells can dramatically restrain tumor growth. Emerging studies suggest that CAFs within the TME have emerged as important determinants of metabolic reprogramming. Metabolic reprogramming also occurs in the metabolic pattern of immune cells. In the meanwhile, immune cell phenotype and functions are metabolically regulated. Notably, immune cell functions influenced by metabolic programs may ultimately lead to alterations in tumor immunity. Despite the fact that multiple previous researches have been devoted to studying the interplays between different cells in the tumor microenvironment, the complicated relationship between CAFs and immune cells and implications of metabolic reprogramming remains unknown and requires further investigation. In this review, we discuss our current comprehension of metabolic reprogramming of CAFs and immune cells (mainly glucose, amino acid, and lipid metabolism) and crosstalk between them that induces immune responses, and we also highlight their contributions to tumorigenesis and progression. Furthermore, we underscore potential therapeutic opportunities arising from metabolism dysregulation and metabolic crosstalk, focusing on strategies targeting CAFs and immune cell metabolic crosstalk in cancer immunotherapy.

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.

Metabolic requirement for GOT2 in pancreatic cancer depends on environmental context

Mitochondrial glutamate-oxaloacetate transaminase 2 (GOT2) is part of the malate-aspartate shuttle, a mechanism by which cells transfer reducing equivalents from the cytosol to the mitochondria. GOT2 is a key component of mutant KRAS (KRAS*)-mediated rewiring of glutamine metabolism in pancreatic ductal adenocarcinoma (PDA). Here, we demonstrate that the loss of GOT2 disturbs redox homeostasis and halts proliferation of PDA cells in vitro. GOT2 knockdown (KD) in PDA cell lines in vitro induced NADH accumulation, decreased Asp and α-ketoglutarate (αKG) production, stalled glycolysis, disrupted the TCA cycle, and impaired proliferation. Oxidizing NADH through chemical or genetic means resolved the redox imbalance induced by GOT2 KD, permitting sustained proliferation. Despite a strong in vitro inhibitory phenotype, loss of GOT2 had no effect on tumor growth in xenograft PDA or autochthonous mouse models. We show that cancer-associated fibroblasts (CAFs), a major component of the pancreatic tumor microenvironment (TME), release the redox active metabolite pyruvate, and culturing GOT2 KD cells in CAF conditioned media (CM) rescued proliferation in vitro. Furthermore, blocking pyruvate import or pyruvate-to-lactate reduction prevented rescue of GOT2 KD in vitro by exogenous pyruvate or CAF CM. However, these interventions failed to sensitize xenografts to GOT2 KD in vivo, demonstrating the remarkable plasticity and differential metabolism deployed by PDA cells in vitro and in vivo. This emphasizes how the environmental context of distinct pre-clinical models impacts both cell-intrinsic metabolic rewiring and metabolic crosstalk with the TME.

The Desmoplastic Stroma of Pancreatic Cancer: Multilayered Levels of Heterogeneity, Clinical Significance, and Therapeutic Opportunities

Simple Summary

Pancreatic cancer is a highly malignant disease with treatment resistance to standardized chemotherapies. In addition, only a small fraction of patients with pancreatic cancer has, to date, actionable genetic aberrations, leading to a narrow therapeutic window for molecularly targeted therapies or immunotherapies. A lot of preclinical and translational studies are ongoing to discover potential vulnerabilities to treat pancreatic cancer. Histologically, human pancreatic cancer is characterized by abundant cancer-associated fibrotic stroma, called “desmoplastic stroma”. Recent technological advances have revealed that desmoplastic stroma in pancreatic cancer is much more complicated than previously thought, playing pleiotropic roles in manipulating tumor cell fate and anti-tumor immunity. Moreover, real-world specimen-based analyses of pancreatic cancer stroma have also uncovered spatial heterogeneity and an intertumoral variety associated with molecular alterations, clinicopathological factors, and patient outcomes. This review describes an overview of the current efforts in the field of pancreatic cancer stromal biology and discusses treatment opportunities of stroma-modifying therapies against this hard-to-treat cancer.

Abstract

Pancreatic cancer remains one of the most lethal malignancies and is becoming a dramatically increasing cause of cancer-related mortality worldwide. Abundant desmoplastic stroma is a histological hallmark of pancreatic ductal adenocarcinoma. Emerging evidence suggests a promising therapeutic effect of several stroma-modifying therapies that target desmoplastic stromal elements in the pancreatic cancer microenvironment. The evidence also unveils multifaceted roles of cancer-associated fibroblasts (CAFs) in manipulating pancreatic cancer progression, immunity, and chemotherapeutic response. Current state-of-the-art technologies, including single-cell transcriptomics and multiplexed tissue imaging techniques, have provided a more profound knowledge of CAF heterogeneity in real-world specimens from pancreatic cancer patients, as well as in genetically engineered mouse models. In this review, we describe recent advances in the understanding of the molecular pathology bases of pancreatic cancer desmoplastic stroma at multilayered levels of heterogeneity, namely, (1) variations in cellular and non-cellular members, including CAF subtypes and extracellular matrix (ECM) proteins; (2) geographical heterogeneity in relation to cell–cell interactions and signaling pathways at niche levels and spatial heterogeneity at locoregional levels or organ levels; and (3) intertumoral stromal heterogeneity at individual levels. This review further discusses the clinicopathological significance of desmoplastic stroma and the potential opportunities for stroma-targeted therapies against this lethal malignancy.

Distinct mechanisms of innate and adaptive immune regulation underlie poor oncologic outcomes associated with KRAS-TP53 co-alteration in pancreatic cancer

Co-occurrent KRAS and TP53 mutations define a majority of patients with pancreatic ductal adenocarcinoma (PDAC) and define its pro-metastatic proclivity. Here, we demonstrate that KRAS-TP53 co-alteration is associated with worse survival compared with either KRAS-alone or TP53-alone altered PDAC in 245 patients with metastatic disease treated at a tertiary referral cancer center, and validate this observation in two independent molecularly annotated datasets. Compared with non-TP53 mutated KRAS-altered tumors, KRAS-TP53 co-alteration engenders disproportionately innate immune-enriched and CD8+ T-cell-excluded immune signatures. Leveraging in silico, in vitro, and in vivo models of human and murine PDAC, we discover a novel intersection between KRAS-TP53 co-altered transcriptomes, TP63-defined squamous trans-differentiation, and myeloid-cell migration into the tumor microenvironment. Comparison of single-cell transcriptomes between KRAS-TP53 co-altered and KRAS-altered/TP53WT tumors revealed cancer cell-autonomous transcriptional programs that orchestrate innate immune trafficking and function. Moreover, we uncover granulocyte-derived inflammasome activation and TNF signaling as putative paracrine mediators of innate immunoregulatory transcriptional programs in KRAS-TP53 co-altered PDAC. Immune subtyping of KRAS-TP53 co-altered PDAC reveals conflation of intratumor heterogeneity with progenitor-like stemness properties. Coalescing these distinct molecular characteristics into a KRAS-TP53 co-altered "immunoregulatory program" predicts chemoresistance in metastatic PDAC patients enrolled in the COMPASS trial, as well as worse overall survival.

Heterogeneity and plasticity of cancer-associated fibroblasts in the pancreatic tumor microenvironment

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive disease with a notably poor prognosis, in urgent need of improved treatment strategies. The desmoplastic PDAC tumor microenvironment (TME), marked by a high concentration of cancer-associated-fibroblasts (CAFs), is a dynamic part of PDAC pathophysiology which occasions a variety of effects throughout the course of pancreatic tumorigenesis and disease evolution. A better understanding of the desmoplastic TME and CAF biology in particular, should provide new opportunities for improving therapeutics. That CAFs have a tumor-supportive role in oncogenesis is well known, yet research evidence has shown that CAFs also have tumor-repressive functions. In this review, we seek to clarify the intriguing heterogeneity and plasticity of CAFs and their ambivalent role in PDAC tumorigenesis and progression. Additionally, we provide recommendations to advance the implementation of CAF-directed PDAC care. An improved understanding of CAFs' origins, spatial location, functional diversity, and marker determination, as well as CAF behavior during the course of PDAC progression and metastasis will provide essential knowledge for the future improvement of therapeutic strategies for patients suffering from PDAC.

Cancer-associated fibroblasts in the single-cell era

Cancer-associated fibroblasts (CAFs) are central players in the microenvironment of solid tumors, affecting cancer progression and metastasis. CAFs have diverse phenotypes, origins and functions and consist of distinct subpopulations. Recent progress in single-cell RNA-sequencing technologies has enabled detailed characterization of the complexity and heterogeneity of CAF subpopulations in multiple tumor types. In this Review, we discuss the current understanding of CAF subsets and functions as elucidated by single-cell technologies, their functional plasticity, and their emergent shared and organ-specific features that could potentially be harnessed to design better therapeutic strategies for cancer.

Tracing back to one of the origins of immune evasion in pancreatic cancer

In this issue of Cancer Cell, Huang et al. identify the origin and function of a subset of pancreatic cancer-associated fibroblasts. This work further highlights the power of lineage-tracing models for dissecting the tumor microenvironment in pancreatic cancer and provides new knowledge toward the development of novel therapeutic strategies.

Mesothelial cell-derived antigen-presenting cancer-associated fibroblasts induce expansion of regulatory T cells in pancreatic cancer

Recent studies have identified a unique cancer-associated fibroblast (CAF) population termed antigen-presenting CAFs (apCAFs), characterized by the expression of major histocompatibility complex class II molecules, suggesting a function in regulating tumor immunity. Here, by integrating multiple single-cell RNA-sequencing studies and performing robust lineage-tracing assays, we find that apCAFs are derived from mesothelial cells. During pancreatic cancer progression, mesothelial cells form apCAFs by downregulating mesothelial features and gaining fibroblastic features, a process induced by interleukin-1 and transforming growth factor β. apCAFs directly ligate and induce naive CD4+ T cells into regulatory T cells (Tregs) in an antigen-specific manner. Moreover, treatment with an antibody targeting the mesothelial cell marker mesothelin can effectively inhibit mesothelial cell to apCAF transition and Treg formation induced by apCAFs. Taken together, our study elucidates how mesothelial cells may contribute to immune evasion in pancreatic cancer and provides insight on strategies to enhance cancer immune therapy.

Eicosanoids in the pancreatic tumor microenvironment - a multicellular, multifaceted progression

BACKGROUND AND AIMS

Eicosanoids, oxidized fatty acids that serve as cell-signaling molecules, have been broadly implicated in tumorigenesis. Here, we aimed to identify eicosanoids associated with pancreatic tumorigenesis and the cell types responsible for their synthesis.

METHODS

We profiled normal pancreas and pancreatic ductal adenocarcinoma (PDAC) in mouse models and patient samples using mass spectrometry. We interrogated RNA sequencing datasets for eicosanoid synthase or receptor expression. Findings were confirmed by immunostaining.

RESULTS

In murine models, we identified elevated levels of PGD₂, prostacyclin, and thromboxanes in neoplasia while PGE₂, 12-HHTre, HETEs, and HDoHEs are elevated specifically in tumors. Analysis of scRNA-seq datasets suggests that PGE₂ and prostacyclins are derived from fibroblasts, PGD₂ and thromboxanes from myeloid cells, and PGD₂ and 5-HETE from tuft cells. In patient samples, we identified a transition from PGD₂ to PGE₂-producing enzymes in the epithelium during the transition to PDAC, fibroblast/tumor expression of PTGIS, and myeloid/tumor cell expression of TBXAS1.

CONCLUSIONS

Our analyses identify key changes in eicosanoid species during pancreatic tumorigenesis and the cell types that contribute to their synthesis. Thromboxane and prostacyclin expression is conserved between animal models and human disease and may represent new druggable targets.

Identification of Functional Heterogeneity of Carcinoma-Associated Fibroblasts with Distinct IL6-Mediated Therapy Resistance in Pancreatic Cancer

The tumor microenvironment in pancreatic ductal adenocarcinoma (PDAC) involves a significant accumulation of fibroblasts as part of the host response to cancer. Using single-cell RNA sequencing, multiplex immunostaining, and several genetic mouse models, we identify carcinoma-associated fibroblasts (CAF) with opposing functions in PDAC progression. Depletion of fibroblast activation protein (FAP)+ CAFs results in increased survival, in contrast to depletion of alpha smooth muscle actin (αSMA)+ CAFs, which leads to decreased survival. Tumor-promoting FAP+ CAFs (TP-CAF) and tumor-restraining αSMA+ CAFs (TR-CAF) differentially regulate cancer-associated pathways and accumulation of regulatory T cells. Improved efficacy of gemcitabine is observed when IL6 is deleted from αSMA+ CAFs but not from FAP+ CAFs using dual-recombinase genetic PDAC models. Improved gemcitabine efficacy due to lack of IL6 synergizes with anti-PD-1 immunotherapy to significantly improve survival of PDAC mice. Our study identifies functional heterogeneity of CAFs in PDAC progression and their different roles in therapy response.

SIGNIFICANCE

PDAC is associated with accumulation of dense stroma consisting of fibroblasts and extracellular matrix that regulate tumor progression. Here, we identify two distinct populations of fibroblasts with opposing roles in the progression and immune landscape of PDAC. Our findings demonstrate that fibroblasts are functionally diverse with therapeutic implications. This article is highlighted in the In This Issue feature, p. 1397.

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.

Current Pathology Model of Pancreatic Cancer

Pancreatic cancer (PC) is one of the most aggressive and lethal malignant neoplasms, ranking in seventh place in the world in terms of the incidence of death, with overall 5-year survival rates still below 10%. The knowledge about PC pathomechanisms is rapidly expanding. Daily reports reveal new aspects of tumor biology, including its molecular and morphological heterogeneity, explain complicated "cross-talk" that happens between the cancer cells and tumor stroma, or the nature of the PC-associated neural remodeling (PANR). Staying up-to-date is hard and crucial at the same time. In this review, we are focusing on a comprehensive summary of PC aspects that are important in pathologic reporting, impact patients' outcomes, and bring meaningful information for clinicians. Finally, we show promising new trends in diagnostic technologies that might bring a difference in PC early diagnosis.

Cytoglobin attenuates pancreatic cancer growth via scavenging reactive oxygen species

Pancreatic cancer is a highly challenging malignancy with extremely poor prognosis. Cytoglobin (CYGB), a hemeprotein involved in liver fibrosis and cancer development, is expressed in pericytes of all organs. Here, we examined the role of CYGB in the development of pancreatic cancer. CYGB expression appeared predominately in the area surrounding adenocarcinoma and negatively correlated with tumor size in patients with pancreatic cancer. Directly injecting 7, 12-dimethylbenz[a]anthracene into the pancreatic tail in wild-type mice resulted in time-dependent induction of severe pancreatitis, fibrosis, and oxidative damage, which was rescued by Cygb overexpression in transgenic mice. Pancreatic cancer incidence was 93% in wild-type mice but only 55% in transgenic mice. Enhanced CYGB expression in human pancreatic stellate cells in vitro reduced cellular collagen synthesis, inhibited cell activation, increased expression of antioxidant-related genes, and increased CYGB secretion into the medium. Cygb-overexpressing or recombinant human CYGB (rhCYGB) -treated MIA PaCa-2 cancer cells exhibited dose-dependent cell cycle arrest at the G1 phase, diminished cell migration, and reduction in colony formation. RNA sequencing in rhCYGB-treated MIA PaCa-2 cells revealed downregulation of cell cycle and oxidative phosphorylation pathways. An increase in MIA PaCa-2 cell proliferation and reactive oxygen species production by H₂O₂ challenge was blocked by rhCYGB treatment or Cygb overexpression. PANC-1, OCUP-A2, and BxPC-3 cancer cells showed similar responses to rhCYGB. Known antioxidants N-acetyl cysteine and glutathione also inhibited cancer cell growth. These results demonstrate that CYGB suppresses pancreatic stellate cell activation, pancreatic fibrosis, and tumor growth, suggesting its potential therapeutic application against pancreatic cancer.

MMP-7 marks severe pancreatic cancer and alters tumor cell signaling by proteolytic release of ectodomains

Pancreatic cancer incurs the worst survival rate of the major cancers. High levels of the protease matrix metalloproteinase-7 (MMP-7) in circulation correlate with poor prognosis and limited survival of patients. MMP-7 is required for a key path of pancreatic tumorigenesis in mice and is present throughout tumor progression. Enhancements to chemotherapies are needed for increasing the number of pancreatic tumors that can be removed and for preventing relapses after surgery. With these ends in mind, selective inhibition of MMP-7 may be worth investigation. An anti-MMP-7 monoclonal antibody was recently shown to increase the susceptibility of several pancreatic cancer cell lines to chemotherapeutics, increase their apoptosis, and decrease their migration. MMP-7 activities are most apparent at the surfaces of innate immune, epithelial, and tumor cells. Proteolytic shedding of multiple protein ectodomains by MMP-7 from such cell surfaces influence apoptosis, proliferation, migration, and invasion. These activities warrant targeting of MMP-7 selectively in pancreatic cancer and other tumors of mucosal epithelia. Competitive and non-competitive modes of MMP-7 inhibition are discussed.

Cancer-associated fibroblast-induced lncRNA UPK1A-AS1 confers platinum resistance in pancreatic cancer via efficient double-strand break repair

The tumor stroma of pancreatic ductal adenocarcinoma (PDAC) is characterized by an abundant and heterogeneous population of cancer-associated fibroblasts (CAFs), which are critically involved in chemoresistance. However, the underlying mechanism of CAFs in chemoresistance is unclear. Here, we show that CAFR, a CAF subset derived from platinum-resistant PDAC patients, assumes an iCAF phenotype and produces more IL8 than CAFS isolated from platinum-sensitive PDAC patients. CAFR-derived IL8 promotes oxaliplatin chemoresistance in PDAC. Based on long noncoding RNA (lncRNA) profiling in tumor cells incubated with CAF-CM, we found that UPK1A-AS1, whose expression is directly induced by IL8/NF-kappa B signaling, functions as a chemoresistance-promoting lncRNA and is critical for active IL8-induced oxaliplatin resistance. Impressively, blocking the activation of UPK1A-AS1 expression increases the oxaliplatin sensitivity of tumor cells in vivo. Mechanistically, UPK1A-AS1 strengthens the interaction between Ku70 and Ku80 to facilitate nonhomologous end joining (NHEJ), thereby enhancing DNA double-strand break (DSB) repair. Clinically, UPK1A-AS1 expression is positively correlated with IL8 expression, a poor chemotherapeutic response and a shorter progression-free survival (PFS) time in advanced PDAC patients. Collectively, our study reveals a lncRNA-mediated mechanism of CAF-derived paracrine IL8-dependent oxaliplatin resistance and highlights UPK1A-AS1 as a potential therapeutic target.

Tie2 Receptor in Tumor-Infiltrating Macrophages Is Dispensable for Tumor Angiogenesis and Tumor Relapse after Chemotherapy

Tumor relapse after chemotherapy relies on the reconstruction of damaged tumor vasculature. In this context, proangiogenic Tie2-expressing macrophages have been suggested to serve as crucial instructors of tumor revascularization by secreting angiogenic factors while being closely associated with the vessel wall. Although the proangiogenic nature of Tie2+ macrophages is well described, the functional contribution of macrophage Tie2 expression remains elusive. Here, we employed a Cre-loxP system to specifically delete Tie2 in macrophages. In multiple syngeneic solid tumor models and two distinct chemotherapeutic treatment regimens, macrophage-expressed Tie2 did not contribute to primary tumor growth, tumor revascularization after chemotherapy, tumor recurrence, or metastasis. Exposing cultured murine macrophage cell lines and bone marrow-derived macrophages to hypoxia or stimulating them with Ang2 did not induce expression of Tie2 at the RNA or protein level. Furthermore, a comprehensive meta-analysis of publicly available single cell RNA sequencing datasets of human and murine tumor-infiltrating CD11b+ myeloid cells did not reveal a transcriptionally distinct macrophage population marked by the expression of Tie2. Collectively, these data question the previously reported critical role of Tie2-expressing macrophages for tumor angiogenesis and tumor relapse after chemotherapy. Moreover, lack of Tie2 inducibility and absence of Tie2-positive macrophages in multiple recently published tumor studies refute a possible prognostic value of macrophage-expressed Tie2.

SIGNIFICANCE

Multiple preclinical tumor models, cell stimulation experiments, and meta-analysis of published tumor single cell RNA sequencing data challenge the reported role of Tie2-positive macrophages for tumor angiogenesis, metastasis, and relapse after chemotherapy. See related commentary by Zhang and Brekken, p. 1172.

Loss of Rnf43 Accelerates Kras-Mediated Neoplasia and Remodels the Tumor Immune Microenvironment in Pancreatic Adenocarcinoma

BACKGROUND & AIMS

RNF43 is an E3 ubiquitin ligase that is recurrently mutated in pancreatic ductal adenocarcinoma (PDAC) and precursor cystic neoplasms of the pancreas. The impact of RNF43 mutations on PDAC is poorly understood and autochthonous models have not been characterized sufficiently. In this study, we describe a genetically engineered mouse model (GEMM) of PDAC with conditional expression of oncogenic Kras and deletion of the catalytic domain of Rnf43 in exocrine cells.

METHODS

We generated Ptf1a-Cre;LSL-KrasG12D;Rnf43flox/flox (KRC) and Ptf1a-Cre; LSL-KrasG12D (KC) mice and animal survival was assessed. KRC mice were sacrificed at 2 months, 4 months, and at moribund status followed by analysis of pancreata by single-cell RNA sequencing. Comparative analyses between moribund KRC and a moribund Kras/Tp53-driven PDAC GEMM (KPC) was performed. Cell lines were isolated from KRC and KC tumors and interrogated by cytokine array analyses, ATAC sequencing, and in vitro drug assays. KRC GEMMs were also treated with an anti-CTLA4 neutralizing antibody with treatment response measured by magnetic response imaging.

RESULTS

We demonstrate that KRC mice display a marked increase in incidence of high-grade cystic lesions of the pancreas and PDAC compared with KC. Importantly, KRC mice have a significantly decreased survival compared with KC mice. Using single-cell RNA sequencing, we demonstrated that KRC tumor progression is accompanied by a decrease in macrophages, as well as an increase in T and B lymphocytes, with evidence of increased immune checkpoint molecule expression and affinity maturation, respectively. This was in stark contrast to the tumor immune microenvironment observed in the KPC PDAC GEMM. Furthermore, expression of the chemokine CXCL5 was found to be specifically decreased in KRC cancer cells by means of epigenetic regulation and emerged as a putative candidate for mediating the unique KRC immune landscape.

CONCLUSIONS

The KRC GEMM establishes RNF43 as a bona fide tumor suppressor gene in PDAC. This GEMM features a markedly different immune microenvironment compared with previously reported PDAC GEMMs and puts forth a rationale for an immunotherapy approach in this subset of PDAC cases.

Extrinsic KRAS Signaling Shapes the Pancreatic Microenvironment Through Fibroblast Reprogramming

BACKGROUND & AIMS

Oncogenic Kirsten Rat Sarcoma virus (KRAS) is the hallmark mutation of human pancreatic cancer and a driver of tumorigenesis in genetically engineered mouse models of the disease. Although the tumor cell-intrinsic effects of oncogenic Kras expression have been widely studied, its role in regulating the extensive pancreatic tumor microenvironment is less understood.

METHODS

Using a genetically engineered mouse model of inducible and reversible oncogenic Kras expression and a combination of approaches that include mass cytometry and single-cell RNA sequencing we studied the effect of oncogenic KRAS in the tumor microenvironment.

RESULTS

We have discovered that non-cell autonomous (ie, extrinsic) oncogenic KRAS signaling reprograms pancreatic fibroblasts, activating an inflammatory gene expression program. As a result, fibroblasts become a hub of extracellular signaling, and the main source of cytokines mediating the polarization of protumorigenic macrophages while also preventing tissue repair.

CONCLUSIONS

Our study provides fundamental knowledge on the mechanisms underlying the formation of the fibroinflammatory stroma in pancreatic cancer and highlights stromal pathways with the potential to be exploited therapeutically.

Effects of hyperinsulinemia on pancreatic cancer development and the immune microenvironment revealed through single-cell transcriptomics

BACKGROUND

Hyperinsulinemia is independently associated with increased risk and mortality of pancreatic cancer. We recently reported that genetically reduced insulin production resulted in ~ 50% suppression of pancreatic intraepithelial neoplasia (PanIN) precancerous lesions in mice. However, only female mice remained normoglycemic, and only the gene dosage of the rodent-specific Ins1 alleles was tested in our previous model. Moreover, we did not delve into the molecular and cellular mechanisms associated with modulating hyperinsulinemia.

METHODS

We studied how reduced Ins2 gene dosage affects PanIN lesion development in both male and female Ptf1a^CreER;Kras^LSL-G12D mice lacking the rodent-specific Ins1 gene (Ins1^-/-). We generated control mice having two alleles of the wild-type Ins2 gene (Ptf1a^CreER;Kras^LSL-G12D;Ins1^-/-;Ins2^+/+) and experimental mice having one allele of Ins2 gene (Ptf1a^CreER;Kras^LSL-G12D;Ins1^-/-;Ins2^+/-). We then performed thorough histopathological analyses and single-cell transcriptomics for both genotypes and sexes.

RESULTS

High-fat diet-induced hyperinsulinemia was transiently or modestly reduced in female and male mice, respectively, with only one allele of Ins2. This occurred without dramatically affecting glucose tolerance. Genetic reduction of insulin production resulted in mice with a tendency for less PanIN and acinar-to-ductal metaplasia (ADM) lesions. Using single-cell transcriptomics, we found hyperinsulinemia affected multiple cell types in the pancreas, with the most statistically significant effects on local immune cell types that were highly represented in our sampled cell population. Specifically, hyperinsulinemia modulated pathways associated with protein translation, MAPK-ERK signaling, and PI3K-AKT signaling, which were changed in epithelial cells and subsets of immune cells.

CONCLUSIONS

These data suggest a potential role for the immune microenvironment in hyperinsulinemia-driven PanIN development. Together with our previous work, we propose that mild suppression of insulin levels may be useful in preventing pancreatic cancer by acting on multiple cell types.

cPLA2 blockade attenuates S100A7-mediated breast tumorigenicity by inhibiting the immunosuppressive tumor microenvironment

BACKGROUND

Molecular mechanisms underlying inflammation-associated breast tumor growth are poorly studied. S100A7, a pro-inflammatory molecule has been shown to enhance breast cancer growth and metastasis. However, the S100A7-mediated molecular mechanisms in enhancing tumor growth and metastasis are unclear.

METHODS

Human breast cancer tissue and plasma samples were used to analyze the expression of S100A7, cPLA2, and PGE2. S100A7-overexpressing or downregulated human metastatic breast cancer cells were used to evaluate the S100A7-mediated downstream signaling mechanisms. Bi-transgenic mS100a7a15 overexpression, TNBC C3 (1)/Tag transgenic, and humanized patient-derived xenograft mouse models and cPLA2 inhibitor (AACOCF3) were used to investigate the role of S100A7/cPLA2/PGE2 signaling in tumor growth and metastasis. Additionally, CODEX, a highly advanced multiplexed imaging was employed to delineate the effects of S100A7/cPLA2 inhibition on the recruitment of various immune cells.

RESULTS

In this study, we found that S100A7 and cPLA2 are highly expressed and correlate with decreased overall survival in breast cancer patients. Further mechanistic studies revealed that S100A7/RAGE signaling promotes the expression of cPLA2 to mediate its oncogenic effects. Pharmacological inhibition of cPLA2 suppressed S100A7-mediated tumor growth and metastasis in multiple pre-clinical models including transgenic and humanized patient-derived xenograft (PDX) mouse models. The attenuation of cPLA2 signaling reduced S100A7-mediated recruitment of immune-suppressive myeloid cells in the tumor microenvironment (TME). Interestingly, we discovered that the S100A7/cPLA2 axis enhances the immunosuppressive microenvironment by increasing prostaglandin E2 (PGE2). Furthermore, CO-Detection by indEXing (CODEX) imaging-based analyses revealed that cPLA2 inhibition increased the infiltration of activated and proliferating CD4⁺ and CD8⁺ T cells in the TME. In addition, CD163⁺ tumor associated-macrophages were positively associated with S100A7 and cPLA2 expression in malignant breast cancer patients.

CONCLUSIONS

Our study provides new mechanistic insights on the cross-talk between S100A7/cPLA2 in enhancing breast tumor growth and metastasis by generating an immunosuppressive TME that inhibits the infiltration of cytotoxic T cells. Furthermore, our studies indicate that S100A7/cPLA2 could be used as novel prognostic marker and cPLA2 inhibitors as promising drugs against S100A7-overexpressing aggressive breast cancer.

Single-Cell Transcriptomics Reveals a Conserved Metaplasia Program in Pancreatic Injury

BACKGROUND & AIMS

Acinar to ductal metaplasia (ADM) occurs in the pancreas in response to tissue injury and is a potential precursor for adenocarcinoma. The goal of these studies was to define the populations arising from ADM, the associated transcriptional changes, and markers of disease progression.

METHODS

Acinar cells were lineage-traced with enhanced yellow fluorescent protein (EYFP) to follow their fate post-injury. Transcripts of more than 13,000 EYFP+ cells were determined using single-cell RNA sequencing (scRNA-seq). Developmental trajectories were generated. Data were compared with gastric metaplasia, KrasG12D-induced neoplasia, and human pancreatitis. Results were confirmed by immunostaining and electron microscopy. KrasG12D was expressed in injury-induced ADM using several inducible Cre drivers. Surgical specimens of chronic pancreatitis from 15 patients were evaluated by immunostaining.

RESULTS

scRNA-seq of ADM revealed emergence of a mucin/ductal population resembling gastric pyloric metaplasia. Lineage trajectories suggest that some pyloric metaplasia cells can generate tuft and enteroendocrine cells (EECs). Comparison with KrasG12D-induced ADM identifies populations associated with disease progression. Activation of KrasG12D expression in HNF1B+ or POU2F3+ ADM populations leads to neoplastic transformation and formation of MUC5AC+ gastric-pit-like cells. Human pancreatitis samples also harbor pyloric metaplasia with a similar transcriptional phenotype.

CONCLUSIONS

Under conditions of chronic injury, acinar cells undergo a pyloric-type metaplasia to mucinous progenitor-like populations, which seed disparate tuft cell and EEC lineages. ADM-derived EEC subtypes are diverse. KrasG12D expression is sufficient to drive neoplasia when targeted to injury-induced ADM populations and offers an alternative origin for tumorigenesis. This program is conserved in human pancreatitis, providing insight into early events in pancreas diseases.

Cancer-associated fibroblasts as accomplices to confer therapeutic resistance in cancer

The "seed and soil" concept has reformed paradigms for cancer treatment in the past decade. Accumulating evidence indicates that the intimate crosstalk between cancer cells and stromal cells plays a tremendous role in tumor progression. Cancer-associated fibroblasts (CAFs), the largest population of stroma cells, influence therapeutic effects through diverse mechanisms. Herein, we summarize the recent advances in the versatile functions of CAFs regarding their heterogeneity, and we mainly discuss the pro-tumorigenic functions of CAFs which promote tumorigenesis and confer therapeutic resistance to tumors. Targeting CAFs is emerging as one of the most appealing strategies in anticancer therapies. The endeavors to target or reprogram the specific subtypes of CAFs provide great cancer treatment opportunities, which may provide a better clinical benefit to cancer patients.

Estimation of tumor cell total mRNA expression in 15 cancer types predicts disease progression

Single-cell RNA sequencing studies have suggested that total mRNA content correlates with tumor phenotypes. Technical and analytical challenges, however, have so far impeded at-scale pan-cancer examination of total mRNA content. Here we present a method to quantify tumor-specific total mRNA expression (TmS) from bulk sequencing data, taking into account tumor transcript proportion, purity and ploidy, which are estimated through transcriptomic/genomic deconvolution. We estimate and validate TmS in 6,590 patient tumors across 15 cancer types, identifying significant inter-tumor variability. Across cancers, high TmS is associated with increased risk of disease progression and death. TmS is influenced by cancer-specific patterns of gene alteration and intra-tumor genetic heterogeneity as well as by pan-cancer trends in metabolic dysregulation. Taken together, our results indicate that measuring cell-type-specific total mRNA expression in tumor cells predicts tumor phenotypes and clinical outcomes.

Single-cell RNA transcriptome reveals the intra-tumoral heterogeneity and regulators underlying tumor progression in metastatic pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is the most frequent and aggressive pancreatic tumor characterized by high metastatic risk and special tumor microenvironment. To comprehensively delineate the complex intra-tumoral heterogeneity and the underlying mechanism during metastatic lesions malignant progression, single-cell RNA sequencing (scRNA-seq) was employed. PCA and TSNE were used for dimensionality reduction analysis and cell clustering. Find All Markers function was used to calculate differential genes in each cluster, and Do Heatmap function was used to plot the distribution of differential genes in each cluster. GSVA was employed to assign pathway activity estimates to individual cells. Lineage trajectory progression was inferred by monocle. CNV status was inferred to compare the heterogeneity among patients and subtypes by infercnv. Ligand-receptor interactions were identified by CellPhoneDB, and regulons network of cells was analyzed by SCENIC. Through RNA-sequencing of 6236 individual cells from 5 liver metastatic PDAC lesions, 10 major cell clusters are identified by using unbiased clustering analysis of expression profiling and well-known cell markers. Cells with high CNV level were considered as malignant cells and pathway analyses were carried out to highlight intratumor heterogeneity in PDAC. Pseudotime trajectory analysis revealed that components of multiple tumor-related pathways and transcription factors (TFs) were differentially expressed along PDAC progression. The complex cellular communication suggested potential immunotherapeutic targets in PDAC. Regulon network identified multiple candidates for promising cell-specific transcriptional factors. Finally, metastatic-related genes expression levels and signaling pathways were validated in bulk RNA Sequencing data. This study contributed a comprehensive single-cell transcriptome atlas and contributed into novel insight of intratumor heterogeneity and molecular mechanism in metastatic PDAC.

Elucidation of Tumor-Stromal Heterogeneity and the Ligand-Receptor Interactome by Single-Cell Transcriptomics in Real-world Pancreatic Cancer Biopsies

PURPOSE

Precision medicine approaches in pancreatic ductal adenocarcinoma (PDAC) are imperative for improving disease outcomes. With molecular subtypes of PDAC gaining relevance in the context of therapeutic stratification, the ability to characterize heterogeneity of cancer-specific gene expression patterns is of great interest. In addition, understanding patterns of immune evasion within PDAC is of importance as novel immunotherapeutic strategies are developed.

EXPERIMENTAL DESIGN

Single-cell RNA sequencing (scRNA-seq) is readily applicable to limited biopsies from human primary and metastatic PDAC and identifies most cancers as being an admixture of previously described epithelial transcriptomic subtypes.

RESULTS

Integrative analyses of our data provide an in-depth characterization of the heterogeneity within the tumor microenvironment, including cancer-associated fibroblast subclasses, and predicts for a multitude of ligand-receptor interactions, revealing potential targets for immunotherapy approaches.

CONCLUSIONS

Our analysis demonstrates that the use of de novo biopsies from patients with PDAC paired with scRNA-seq may facilitate therapeutic prediction from limited biopsy samples.

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.

Oncogenic KRAS-Induced Feedback Inflammatory Signaling in Pancreatic Cancer: An Overview and New Therapeutic Opportunities

Pancreatic ductal adenocarcinoma (PDAC) remains highly refractory to treatment. While the KRAS oncogene is present in almost all PDAC cases and accounts for many of the malignant feats of PDAC, targeting KRAS or its canonical, direct effector cascades remains unsuccessful in patients. The recalcitrant nature of PDAC is also heavily influenced by its highly fibro-inflammatory tumor microenvironment (TME), which comprises an acellular extracellular matrix and various types of non-neoplastic cells including fibroblasts, immune cells, and adipocytes, underscoring the critical need to delineate the bidirectional signaling interplay between PDAC cells and the TME in order to develop novel therapeutic strategies. The impact of tumor-cell KRAS signaling on various cell types in the TME has been well covered by several reviews. In this article, we critically reviewed evidence, including work from our group, on how the feedback inflammatory signals from the TME impact and synergize with oncogenic KRAS signaling in PDAC cells, ultimately augmenting their malignant behavior. We discussed past and ongoing clinical trials that target key inflammatory pathways in PDAC and highlight lessons to be learned from outcomes. Lastly, we provided our perspective on the future of developing therapeutic strategies for PDAC through understanding the breadth and complexity of KRAS and the inflammatory signaling network.

Cancer cell immune mimicry delineates onco-immunologic modulation

Immune transcripts are essential for depicting onco-immunologic interactions. However, whether cancer cells mimic immune transcripts to reprogram onco-immunologic interaction remains unclear. Here, single-cell transcriptomic analyses of 7,737 normal and 37,476 cancer cells reveal increased immune transcripts in cancer cells. Cells gradually acquire immune transcripts in malignant transformation. Notably, cancer cell-derived immune transcripts contribute to distinct prognoses of immune gene signatures. Optimized immune response signature (oIRS), obtained by excluding cancer-related immune genes from immune gene signatures, and offers a more reliable prognostic value. oIRS reveals that antigen presentation, NK cell killing and T cell signaling are associated with favorable prognosis. Patients with higher oIRS expression are associated with favorable responses to immunotherapy. Indeed, CD83⁺ cell infiltration, which indicates antigen presentation activity, predicts favorable prognosis in breast cancer. These findings unveil that immune mimicry is a distinct cancer hallmark, providing an example of cancer cell plasticity and a refined view of tumor microenvironment.

•

Single-cell transcriptome reveals immune mimicry as a distinct cancer hallmark

•

Mimicry transcripts contribute to distinct prognoses of immune gene signatures

•

oIRS refines prognostic evaluation and points to core favorable immune processes

•

oIRS signifies the role of antigen presentation and indicates immunotherapy response

Crosstalk between cancer-associated fibroblasts and immune cells in the tumor microenvironment: new findings and future perspectives

Cancer-associated fibroblasts (CAFs), a stromal cell population with cell-of-origin, phenotypic and functional heterogeneity, are the most essential components of the tumor microenvironment (TME). Through multiple pathways, activated CAFs can promote tumor growth, angiogenesis, invasion and metastasis, along with extracellular matrix (ECM) remodeling and even chemoresistance. Numerous previous studies have confirmed the critical role of the interaction between CAFs and tumor cells in tumorigenesis and development. However, recently, the mutual effects of CAFs and the tumor immune microenvironment (TIME) have been identified as another key factor in promoting tumor progression. The TIME mainly consists of distinct immune cell populations in tumor islets and is highly associated with the antitumor immunological state in the TME. CAFs interact with tumor-infiltrating immune cells as well as other immune components within the TIME via the secretion of various cytokines, growth factors, chemokines, exosomes and other effector molecules, consequently shaping an immunosuppressive TME that enables cancer cells to evade surveillance of the immune system. In-depth studies of CAFs and immune microenvironment interactions, particularly the complicated mechanisms connecting CAFs with immune cells, might provide novel strategies for subsequent targeted immunotherapies. Herein, we shed light on recent advances regarding the direct and indirect crosstalk between CAFs and infiltrating immune cells and further summarize the possible immunoinhibitory mechanisms induced by CAFs in the TME. In addition, we present current related CAF-targeting immunotherapies and briefly describe some future perspectives on CAF research in the end.

Heterogeneity in Pancreatic Cancer Fibroblasts-TGFβ as a Master Regulator?

Pancreatic ductal adenocarcinoma is an aggressive disease for which there are very few available therapies. It is notable for its high degree of tumour complexity, with the tumour microenvironment often accounting for the majority of the tumour volume. Until recently, the biology of the stroma was poorly understood, particularly in terms of heterogeneity. Recent research, however, has shed light on the intricacy of signalling within the stroma and particularly the molecular and functional heterogeneity of the cancer associated fibroblasts. In this review, we summarise the recent improvements in our understanding of the different fibroblast populations within PDAC, with a focus on the role TGFβ plays to dictate their formation and function. These studies have highlighted some of the reasons for the failure of trials targeting the tumour stroma, however, there are still considerable gaps in our knowledge, and more work is needed to make effective fibroblast targeting a reality in the clinic.

Mesenchymal Lineage Heterogeneity Underlies Non-Redundant Functions of Pancreatic Cancer-Associated Fibroblasts

Cancer-associated fibroblast (CAF) heterogeneity is increasingly appreciated, but the origins and functions of distinct CAF subtypes remain poorly understood. The abundant and transcriptionally diverse CAF population in pancreatic ductal adenocarcinoma (PDAC) is thought to arise from a common cell of origin, pancreatic stellate cells (PSCs), with diversification resulting from cytokine and growth factor gradients within the tumor microenvironment. Here we analyzed the differentiation and function of PSCs during tumor progression in vivo. Contrary to expectations, we found that PSCs give rise to a numerically minor subset of PDAC CAFs. Targeted ablation of PSC-derived CAFs within their host tissue revealed non-redundant functions for this defined CAF population in shaping the PDAC microenvironment, including production of specific extracellular matrix components and tissue stiffness regulation. Together, these findings link stromal evolution from distinct cells of origin to transcriptional heterogeneity among PDAC CAFs, and demonstrate unique functions for CAFs of a defined cellular origin.

Genomic Heterogeneity of Pancreatic Ductal Adenocarcinoma and Its Clinical Impact

Pancreatic ductal adenocarcinoma (PDAC) is one of the leading causes of cancer death due to limited advances in recent years in early diagnosis and personalized therapy capable of overcoming tumor resistance to chemotherapy. In the last decades, significant advances have been achieved in the identification of recurrent genetic and molecular alterations of PDAC including those involving the KRAS, CDKN2A, SMAD4, and TP53 driver genes. Despite these common genetic traits, PDAC are highly heterogeneous tumors at both the inter- and intra-tumoral genomic level, which might contribute to distinct tumor behavior and response to therapy, with variable patient outcomes. Despite this, genetic and genomic data on PDAC has had a limited impact on the clinical management of patients. Integration of genomic data for classification of PDAC into clinically defined entities-i.e., classical vs. squamous subtypes of PDAC-leading to different treatment approaches has the potential for significantly improving patient outcomes. In this review, we summarize current knowledge about the most relevant genomic subtypes of PDAC including the impact of distinct patterns of intra-tumoral genomic heterogeneity on the classification and clinical and therapeutic management of PDAC.

The Extracellular Matrix in Pancreatic Cancer: Description of a Complex Network and Promising Therapeutic Options

The stroma is a relevant player in driving and supporting the progression of pancreatic ductal adenocarcinoma (PDAC), and a large body of evidence highlights its role in hindering the efficacy of current therapies. In fact, the dense extracellular matrix (ECM) characterizing this tumor acts as a natural physical barrier, impairing drug penetration. Consequently, all of the approaches combining stroma-targeting and anticancer therapy constitute an appealing option for improving drug penetration. Several strategies have been adopted in order to target the PDAC stroma, such as the depletion of ECM components and the targeting of cancer-associated fibroblasts (CAFs), which are responsible for the increased matrix deposition in cancer. Additionally, the leaky and collapsing blood vessels characterizing the tumor might be normalized, thus restoring blood perfusion and allowing drug penetration. Even though many stroma-targeting strategies have reported disappointing results in clinical trials, the ECM offers a wide range of potential therapeutic targets that are now being investigated. The dense ECM might be bypassed by implementing nanoparticle-based systems or by using mesenchymal stem cells as drug carriers. The present review aims to provide an overview of the principal mechanisms involved in the ECM remodeling and of new promising therapeutic strategies for PDAC.

A Transcriptionally Distinct Subpopulation of Healthy Acinar Cells Exhibit Features of Pancreatic Progenitors and PDAC

Pancreatic ductal adenocarcinoma (PDAC) tumors can originate either from acinar or ductal cells in the adult pancreas. We re-analyze multiple pancreas and PDAC single-cell RNA-seq datasets and find a subset of nonmalignant acinar cells, which we refer to as acinar edge (AE) cells, whose transcriptomes highly diverge from a typical acinar cell in each dataset. Genes upregulated among AE cells are enriched for transcriptomic signatures of pancreatic progenitors, acinar dedifferentiation, and several oncogenic programs. AE-upregulated genes are upregulated in human PDAC tumors, and consistently, their promoters are hypomethylated. High expression of these genes is associated with poor patient survival. The fraction of AE-like cells increases with age in healthy pancreatic tissue, which is not explained by clonal mutations, thus pointing to a nongenetic source of variation. The fraction of AE-like cells is also significantly higher in human pancreatitis samples. Finally, we find edge-like states in lung, liver, prostate, and colon tissues, suggesting that subpopulations of healthy cells across tissues can exist in pre-neoplastic states. SIGNIFICANCE: These findings show "edge" epithelial cell states with oncogenic transcriptional activity in human organs without oncogenic mutations. In the pancreas, the fraction of acinar cells increases with age.

Myeloid Cell Mediated Immune Suppression in Pancreatic Cancer

Pancreatic ductal adenocarcinoma (PDA), the most common pancreatic cancer, is a nearly universally lethal malignancy. PDA is characterized by extensive infiltration of immunosuppressive myeloid cells, including tumor-associated macrophages and myeloid-derived suppressor cells. Myeloid cells in the tumor microenvironment inhibit cytotoxic T-cell responses promoting carcinogenesis. Immune checkpoint therapy has not been effective in PDA, most likely because of this robust immune suppression, making it critical to elucidate mechanisms behind this phenomenon. Here, we review myeloid cell infiltration and cellular crosstalk in PDA progression and highlight current therapeutic approaches to target myeloid cell-driven immune suppression.

Cancer-Associated Fibroblast (CAF) Heterogeneity and Targeting Therapy of CAFs in Pancreatic Cancer

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal disease that typically features a dramatic desmoplastic reaction, especially fibroblasts. The roles of cancer-associated fibroblasts (CAFs) in PDAC have received more attention in recent years. As increasing evidence suggests the heterogeneity of CAFs in PDAC, different CAF subtypes have been shown to support tumor growth, while others suppress cancer proliferation. Myofibrotic CAFs (myCAFs) show alpha-smooth muscle actin (α-SMA)^high interleukin-6 (IL-6)^low myofibroblastic features, are activated by direct contact with tumor cells, and are located in the periglandular region. Inflammatory CAFs (iCAFs) show α-SMA^low IL-6^high inflammatory features, are activated by paracrine factors secreted from tumor cells, and are located away from cancer cells. Antigen-presenting CAFs (apCAFs) show major histocompatibility complex II (MHC II) family genes that are highly expressed. CAFs have also been gradually explored as diagnostic and prognostic markers in pancreatic cancer. Targeted therapy of CAFs in PDAC has gradually attracted attention. With the deepening of related studies, some meaningful positive and negative results have surfaced, and CAFs may be the key to unlocking the door to pancreatic cancer treatment. Our review summarizes recent advances in the heterogeneity, function, and markers of CAFs in pancreatic cancer, as well as research and treatment targeting CAFs in pancreatic cancer.