REVIEW article

Front. Cell Dev. Biol., 05 November 2021

Sec. Molecular and Cellular Pathology

Volume 9 - 2021 | https://doi.org/10.3389/fcell.2021.743908

Clinical Impact of Molecular Subtyping of Pancreatic Cancer

    XZ

    Xu Zhou 1,2

    KH

    Kai Hu 1,2

    PB

    Peter Bailey 1,2,3

    CS

    Christoph Springfeld 4

    SR

    Susanne Roth 1

    RK

    Roma Kurilov 5

    BB

    Benedikt Brors 4,5

    TG

    Thomas Gress 6

    MB

    Malte Buchholz 6

    JA

    Jingyu An 1,2

    KW

    Kongyuan Wei 1,2

    TP

    Teresa Peccerella 1,2

    MW

    Markus W. Büchler 1

    TH

    Thilo Hackert 1

    JP

    John P. Neoptolemos 1,2*

  • 1. Department of General, Visceral and Transplantation Surgery, Heidelberg University Hospital, Heidelberg, Germany

  • 2. Section of Surgical Research, Heidelberg University Hospital, Heidelberg, Germany

  • 3. Institute of Cancer Sciences, University of Glasgow, Glasgow, United Kingdom

  • 4. Department of Medical Oncology, National Center for Tumor Diseases, Heidelberg University Hospital, Heidelberg, Germany

  • 5. Division of Applied Bioinformatics, German Cancer Research Center, Heidelberg, Germany

  • 6. Department of Gastroenterology and Endocrinology, Philipps University of Marburg, Marburg, Germany

Article metrics

View details

69

Citations

20,4k

Views

6,6k

Downloads

Abstract

Pancreatic ductal adenocarcinoma is a highly lethal malignancy, which has now become the seventh most common cause of cancer death in the world, with the highest mortality rates in Europe and North America. In the past 30 years, there has been some progress in 5-year survival (rates increasing from 2.5 to 10%), but this is still extremely poor compared to all other common cancer types. Targeted therapies for advanced pancreatic cancer based on actionable mutations have been disappointing, with only 3–5% showing even a short clinical benefit. There is, however, a molecular diversity beyond mutations in genes responsible for producing classical canonical signaling pathways. Pancreatic cancer is almost unique in promoting an excess production of other components of the stroma, resulting in a complex tumor microenvironment that contributes to tumor development, progression, and response to treatment. Various transcriptional subtypes have also been described. Most notably, there is a strong alignment between the Classical/Pancreatic progenitor and Quasi-mesenchymal/Basal-like/Squamous subtype signatures of Moffit, Collinson, Bailey, Puleo, and Chan-Seng-Yue, which have potential clinical impact. Sequencing of epithelial cell populations enriched by laser capture microscopy combined with single-cell RNA sequencing has revealed the potential genomic evolution of pancreatic cancer as being a consequence of a gene expression continuum from mixed Basal-like and Classical cell populations within the same tumor, linked to allelic imbalances in mutant KRAS, with metastatic tumors being more copy number-unstable compared to primary tumors. The Basal-like subtype appears more chemoresistant with reduced survival compared to the Classical subtype. Chemotherapy and/or chemoradiation will also enrich the Basal-like subtype. Squamous/Basal-like programs facilitate immune infiltration compared with the Classical-like programs. The immune infiltrates associated with Basal and Classical type cells are distinct, potentially opening the door to differential strategies. Single-cell and spatial transcriptomics will now allow single cell profiling of tumor and resident immune cell populations that may further advance subtyping. Multiple clinical trials have been launched based on transcriptomic response signatures and molecular subtyping including COMPASS, Precision Promise, ESPAC6/7, PREDICT-PACA, and PASS1. We review several approaches to explore the clinical relevance of molecular profiling to provide optimal bench-to-beside translation with clinical impact.

Introduction

Pancreatic ductal adenocarcinoma (PDAC), a distinct form of pancreatic cancer, remains a major oncological challenge (Kleeff et al., 2016). Globally the 5-year pancreatic cancer prevalence in 2020 was 4.87 per 105 per year (International Agency for Research on Cancer, 2021). The number of cases of pancreatic cancer worldwide in 2020 was 495,773 (world rank for all cancers = 13), with 466,003 deaths (world rank for all cancers = 7); incidence rates per 105 per year were 5.7 for men and 4.1 for women, and mortality rates of 4.9 and 4.5, respectively (International Agency for Research on Cancer, 2021). In Europe, there were 140,116 new cases with 132,134 deaths (International Agency for Research on Cancer, 2021). In North America, there were 62,643 new cases and 53,277 deaths, the fourth highest cancer mortality in both men and women (International Agency for Research on Cancer, 2021; Siegel et al., 2021). For Western Europe, the incidence rates per 105 per year were 9.9 for men and 7.4 for women, with mortality rates of 8.6 and 7.8, respectively (world rank first for pancreatic cancer) (International Agency for Research on Cancer, 2021). In 2017, in Germany, there were 18,687 new cases (with a rising incidence rate) and 18,005 deaths, but with a slight improvement in the 5-year survival rate from 8% in 2007–2008 to 9% in 2015–2016 (Zentrum für Krebsregisterdaten, 2021). In the United States, the 5-year survival rate for all stages has further improved to 10% (Siegel et al., 2021).

Standard Therapies for Pancreatic Cancer

Most patients present with metastatic disease, with only 10–20% being diagnosed with localized pancreatic cancer that can be surgically removed, while the remaining 20–30% have non-metastatic locally advanced disease that cannot be removed by standard surgical techniques (Kleeff et al., 2016). Systemic chemotherapy is the only conventional approach for improving survival in patients with advanced disease with the best achievable median survival rates being 8–12 months for metastatic disease and 12–15 months for locally advanced pancreatic cancer (Burris et al., 1997; Cunningham et al., 2009; Conroy et al., 2011; Von Hoff et al., 2013; Gill et al., 2016; Wang-Gillam et al., 2016; Springfeld et al., 2019). Although chemoradiotherapy is often used for locally advanced disease especially in the United States, there is increased toxicity without improvement in overall survival (Sultana et al., 2007; Chauffert et al., 2008; Hammel et al., 2016). In patients with locally resectable tumors but without metastatic disease, advances in surgical techniques and the use of adjuvant systematic chemotherapy have increased 5-year survival rates from 8% with resection alone to 30–50% in conjunction with adjuvant chemotherapy most notably using gemcitabine and capecitabine or modified folinic acid, 5-fluorouracil (5-FU), irinotecan, and oxaliplatin (mFOLFIRINOX) combinations (Neoptolemos et al., 2004; Oettle et al., 2007; Neoptolemos et al., 2010; Kleeff et al., 2016; Uesaka et al., 2016; Neoptolemos et al., 2017; Conroy et al., 2018; Strobel et al., 2019). Patients with borderline resectable disease may benefit from neoadjuvant chemotherapy regimens comprising gemcitabine with capecitabine as well as mFOLFIRINOX, while regimens with chemoradiotherapy are inferior to chemotherapy alone (Ghaneh et al., 2020; Katz et al., 2021). Neoadjuvant therapy may also increase resectabilty with improved survival in patients with otherwise unresectable local disease due to major vessel encasement using comprising mFOLFIRINOX or gemcitabine-based regimens with either capecitabine or nab-paclitaxel (Hackert et al., 2016; Diener et al., 2021; Kunzmann et al., 2021). An argument has been made to extend the use of neoadjuvant chemotherapy to patients with resectable disease, but this appears to be inferior for overall survival compared to upfront surgery and adjuvant treatment (Sohal et al., 2021).

The evolution of molecular targeted therapies, aimed at advancing tumor control and cell killing of pancreatic cancer, has so far met with only very limited progress (Davis et al., 2019; Golan et al., 2019; Mosele et al., 2020; O’Kane et al., 2020; Pishvaian et al., 2020; Cobain et al., 2021). While systemic chemotherapy is the mainstay of treatment when added to surgery, its impact is limited by a wide variation in responsiveness that is related to intrinsic and acquired mechanisms of sensitivity and resistance by both the cancer cells themselves and the stromal environment (Greenhalf et al., 2014; Martinez-Balibrea et al., 2015; Noll et al., 2016; Amrutkar and Gladhaug, 2017; Geller et al., 2017; Martinelli et al., 2017; Schlitter et al., 2017; Neoptolemos et al., 2018; Tiriac et al., 2018; Dominguez et al., 2020; Kalimuthu et al., 2020; Sahai et al., 2020). Going one step further, the integration of molecular subtypes derived from global genomic and transcriptomic analyses into clinical trials is enabling translational insights into how we might better refine existing and evolving therapy modalities to improve pancreatic cancer treatment. Pancreatic cancer is almost unique in promoting an excess production of other components of the admixture of general tissue (stroma), resulting in a complex tumor microenvironment that contributes to tumor development, progression, and response to treatment (Kleeff et al., 2016).

Single Gene Alterations in Pancreatic Cancer

Mutations in genes responsible for producing classical canonical signaling pathways including driver oncogenes and dysfunction of tumor suppressor genes include KRAS, TP53, CDKN2A, and SMAD4 in most cases, and ARID1A, KDM6A, MLL3, TGFBR2, RBM10, BCORL1, and ROBO2 in 5–10% of tumors (Jones et al., 2008; Waddell et al., 2015; Bailey et al., 2016; Chan-Seng-Yue et al., 2020; Escobar-Hoyos et al., 2020). Genetic alterations occur in each of a core set of 12 cellular signaling pathways in 67–100% of the tumors, with representative genes listed below (Jones et al., 2008).

  • Apoptosis: CASP10, VCP, CAD, HIP1

  • DNA damage control: ERCC4, ERCC6, EP300, RANBP2, TP53, BRCA1/2, PALB2, ATM, ATR, MLH1, MSH2, MSH6, RPA1, STK11, FANCA, FANCC

  • Regulation of G1/S phase transition: CDKN2A, FBXW7, CHD1, APC2

  • Hedgehog signaling: TBX5, SOX3, LRP2, GLI1, GLI3, BMPR2, CREBBP

  • Homophilic cell adhesion: CDH1, FAT, PCDH15, PCDHB16, PCDHGA1

  • Integrin signaling: ITGA4, LAMA1, LAMA4, LAMA5, FN1, ILK

  • c-Jun N-terminal kinase signaling: MAP4K3, TNF, ATF2, NFATC3

  • KRAS signaling: KRAS, MAP2K4, RASGRP3

  • Regulation of invasion: ADAM11, DPP6, MEP1A, PCSK6, APG4A

  • Small GTPase–dependent signaling: AGHGEF7, ARHGEF9, CDC42BPA

  • TGF-β signaling: TGFBR2, BMPR2, SMAD4, SMAD3

  • Wnt/Notch signaling: MYC, PPP2R3A, WNT9A, MAP2, TSC2, GATA6.

Single genetic alterations occur in <5% of tumors, notably BRCA1/2 mutations, BRAF gene fusions/mutations, ERBB2 amplifications/mutations, RNF43, TGFBR2, MAP2K4, MLL3, PIK3CA, RBM10, SMARCA4, PBRM1, SLIT2, KDM6A, GATA6, BRAF, ATM, and mismatch repair (MMR) gene mutations (Waddell et al., 2015). MMR genes (MLH1, MSH2, MSH6, and PMS2) normally recognize mistakes in insertion, deletion, or mismatched incorporation of nucleotides arising from errors by DNA polymerases and then replacing them with the correct nucleotides. As well as gene mutations, loss of MMR protein function may arise through promoter methylation especially in the case of MLH1. The consequence is an accumulation of errors in DNA microsatellites (short repetitive sequences in DNA) causing high microsatellite instability (MSI-H) (Table 1; Waddell et al., 2015; Mosele et al., 2020).

TABLE 1

GeneAlterationPrevalenceESCAT Level
BRCA1/2Germline mutations1–4%I
BRCA1/2Somatic mutations3%III
MSH1, PMS2, MLH1, and MSH6MSI-H1–3%I
NTRKFusions<1%I
KRASG12CMutation1–2%II
PIK3CAHotspot mutations3%III
BRAFV600EMutations3%III
MDM2Amplifications2%III
ERBB2Amplifications/mutations1–2%III
NRG1Fusions1%III
ALKFusions<1%III
RETFusions<1%III
ROS1Fusions<1%III

List of actionable single gene alterations to in advanced pancreatic ductal adenocarcinoma in accordance with ESMO Scale for Clinical Actionability of Molecular Targets (ESCAT) levels I–III (Mateo et al., 2018; Mosele et al., 2020).

Clinical applicability of genetic biomarkers has been classified by the European Society for Medical Oncology (ESMO) Translational Research and Precision Medicine Working Group into the ESMO Scale of Clinical Actionability for molecular Targets (ESCAT) (Mateo et al., 2018). There are four main levels defined as follows: I = the match of an alteration and a drug has been validated in clinical trials, and should drive treatment decision in daily practice; II = a drug that matches the alteration has been associated with responses in phase I/II or in retrospective analyses of randomized trials; III = alterations that are validated in another cancer, but not in the disease-to-treat; IV = hypothetically targetable alterations based on preclinical data (Mateo et al., 2018; Mosele et al., 2020). So far, the clinical utility of targeting drugs to specific molecular alterations is rather limited (Hong et al., 2020; Mosele et al., 2020):

  • mFOLFIRINOX –preferred for known germline BRCA1/2 or PALB2 mutations.

  • Olaparib – a PARP inhibitor as maintenance therapy in patients who have a germline BRCA1 or BRCA2 mutation and with metastatic pancreatic cancer that had not progressed during first-line platinum-based chemotherapy, resulting in improved progression-free survival.

  • Entrectinib – an inhibitor of tropomyosin receptor kinases (TRKs) of tumors with NTRK or ROS-1 gene fusions.

  • Laroctrenib – an inhibitor of tropomyosin receptor kinases (TRKs) of tumors with NTRK gene fusions.

  • Afatinib – an EGFR tyrosine kinase inhibitor in KRAS wild-type tumors with NRG1 gene fusions.

  • Sotorasib – a small molecule that targets the KRAS p.G12C mutation that is present in 1–2% of PDAC patients (Hong et al., 2020).

Also, erlotinib, a multiple tyrosine kinase inhibitor (including EGFR) used with gemcitabine, produces an improved survival in metastatic pancreatic cancer, but this benefit is only marginal with increased toxicity.

Actionable Genomic Subtypes

Structural variations amongst the 25,000 defined human genomes include deletions, amplifications, duplications, and translocations (International Cancer Genome Consortium et al., 2010; Waddell et al., 2015). The Waddell signature based on whole-genome sequencing and copy number variation identified four subtypes based on patterns of chromosomal structural variation with potential clinical utility (Table 2; Waddell et al., 2015).

TABLE 2

Study nameYearSample source and numberMethodologyTumor cellularitySource of the DNA/RNASubtyping methodMolecular subtypesPromising molecular biomarkersClinical relevance
Waddell et al., 20152015Primary resected tumors (n = 100)WGSCryo bulk tissueStructural rearrangementsStable, locally rearranged, scattered, unstableBRCA mutation (frequent in unstable tumors)Unstable genome and/or BRCA mutation tumors responded well to platinum-based therapy.
Collisson et al., 20112011Primary resected PDAC (n = 27); PDXs (n = 19); Mouse cell lines (n = 15)Gene expression microarrayHighCryo samples underwent microdissectionTumor transcriptional profilesClassical, QM-PDA, exocrine-likeGATA6 (high in classical and low in QM-PDA)QM-PDA subtype was correlated with poor outcome and was sensitive to gemcitabine; classical subtype was correlated with good outcome and was sensitive to erlotinib.
Moffitt et al., 20152015Primary (n = 145) and metastatic (n = 61) PDAC; Cell lines (n = 17); Pancreas (n = 46) and distant adjacent normal samples (n = 88); PDXs (n = 37); CAF (n = 6)Whole-genome DNA microarrays, RNA-seq, virtual microdissectionCryo bulk tissue, cell lines, PDXs, CAFsTumor and stroma transcriptional profilesTumor: classical, basal-like Stroma: normal, activatedSMAD4 and GATA6 (high in classical and low in basal-like)Classical tumors with normal stroma had best outcome and basal-like tumors with activated stroma had worst outcome; basal-like tumors responded to adjuvant therapy better than classical tumors.
Bailey et al., 20162016Primary resected tumors (n = 456); Patient-derived cell lines (n = 41); Mouse cell linesWGS, deep-exome sequencing, RNA-seq>40% for WGS, 12–40% for deep-exome sequencingCryo bulk tissueTumor transcriptional profilesSquamous, pancreatic progenitor, immunogenic, ADEXTP53 and KDM6A mutations (frequent in squamous tumors); FOXA2/3, PDX1 and MNX1 (high in progenitor tumors)Squamous tumors had worst outcome.
Raphael and Cancer Genome Atlas Research Network, 20172017Primary resected PDAC (n = 150)WES, custom targeted gene panel sequencing, RNA-seqCryo bulk tissueTumor transcriptional profilesValidation of former classifications (Moffitt’s classification was independent of tumor purity, while Collisson’s and Bailey’s classifications were correlated with tumor purity)Low EMT and apoptosis pathway activity, high TSC-mTOR and RTK activity in better survival groups
Wartenberg et al., 20182018Primary resected PDAC (n = 110)ngTMAFFPE samplesTumor transcriptional profilesImmune escape, immune rich and immune exhaustedImmune escape tumors had poor outcome, and immune rich tumors had better outcome
Puleo et al., 20182018Primary resected PDAC (n = 381)RNA microarray, immunohistochemistry, DNA panel sequencingFFPE samplesTumor transcriptional profilesPure classical, immune classical, desmoplastic, stroma activated, pure basal-likePure basal-like tumors had the worst outcome.
Maurer et al., 20192019Primary resected PDAC (n = 60)RNA-seqHighCryo samples underwent LCMTumor and stroma transcriptional profilesTumor: classical, basal-like Stroma: immune-rich, ECM-richBasal-like/ECM-rich tumors had worse outcome than classical/immune-rich tumors
Dijk et al., 20202020Primary resected PDAC (n = 90)RNA-seqCryo bulk tissue, PDXsTumor transcriptional profilesSecretory, epithelial, compound pancreatic and mesenchymalKRAS (highest transcript level in mesenchymal subtype)Mesenchymal and secretory tumors had worse outcome than epithelial and compound pancreatic tumors.
Chan-Seng-Yue et al., 20202020Primary resected PDAC (n = 206) and advanced PDAC (n = 111)WGS, RNA-seq, single cell analysisHighCryo samples underwent LCM, single cellTumor transcriptional profilesBasal-like-A, basal-like-B, hybrid, classical-A, classical-BSMAD4 (high in basal-like-A tumors and low in classical-A tumors), CDKN2A and TP53 (completely loss in basal-like-A/B tumors), GATA6 (high in Classical-A/B tumors)Classical-A/B tumors were more frequent in early stage, while basal-like tumors in late stage; in resectable disease, basal-like-B and hybrid tumors identified two prognostic subgroups considered to be uniformly aggressive before; in advanced disease, basal-like-A was highly chemo-resistant and trended toward worse survival.

Main molecular subtypes of pancreatic cancer.

  • Stable (20%), with <50 structural variations per genome, with widespread aneuploidy.

  • Locally rearranged (30%), with >200 structural variants clustered on 1–2 chromosomes. Of these, about a 35% had focal amplifications in KRAS, SOX9 and GATA6, as well as ERBB2, MET, CDK6, PIK3CA, and PIK3R3 but were only present in 1–2% of patients. The remaining local rearrangements involved complex genomic events such as breakage–fusion–bridge or chromothripsis (thousands of clustered chromosomal rearrangements occurring in a single event in localized and confined genomic regions in one or two chromosomes).

  • Scattered (36%), with non-random chromosomal damage in 50–200 structural variants per genome.

  • Unstable (14%), with >200 structural variants distributed across the genome indicating defects in DNA maintenance (BRCA1/2, and PALB2 gene defects) and a mutational DNA damage repair (DDR) deficiency, with potential sensitivity to DNA-damaging agents. The unstable structural variation subtype is responsive to platinum therapy and BRCA1/2 germline carriers also sensitive to both platinum and PARP inhibitors.

It is estimated that 24% of all pancreatic cancers may be sensitive to platinum therapy based on an unstable genomic structural variation subtype, and/or somatic and germline mutations in BRCA genes, and/or a BRCA-type mutational signature (Waddell et al., 2015). MSI-H occurs in 1–3% of pancreatic cancers, which is commonly associated with mutations in the MSH2 and MLH1 MMR genes, and can be detected by immunohistochemistry (MSH1, PMS2, MLH1, and MSH6 expression) or sequencing (single gene mutations and MMR mutational signature) (Waddell et al., 2015; Connor et al., 2017).

MSI-H tumors express a large number of neoantigens, potentially rendering them more susceptible to immunotherapy in comparison to those tumors with relatively few mutations. DNA replication stress producing single-stranded DNA will induce DDR of which the DNA damage checkpoint kinase ATR [Ataxia-Telangiectasia Mutated (ATM) and Rad3-related protein kinase] is a critical component. Genes encoding subunits of SWI/SNF (BAF) chromatin remodeling complexes, each composed of approximately 15 protein subunits, include ARID1A, ARID1B, ARID2, PBRM1, SMARCA4, and SMARCB1. ARID1A deficiency will impair cells to recruit topoisomerase 2A to chromatin causing cell cycle defects. The consequence is increased reliance on ATR checkpoint activity and thereby increased sensitivity to ATR inhibitor therapy (Williamson et al., 2016).

Transcriptomic Subtypes

Various transcriptional pancreatic cancer subtypes have also been described, most notably the Moffit, Collinson, Bailey, Puleo, and Chan-Seng-Yue signatures amongst others, which have potential clinical impact (Table 2; Collisson et al., 2011; Moffitt et al., 2015; Bailey et al., 2016; Cancer Genome Atlas Research Network, 2017; Puleo et al., 2018; Wartenberg et al., 2018; Maurer et al., 2019; Chan-Seng-Yue et al., 2020; Dijk et al., 2020). Each study has used a different approach to deal with the low cellularity and stromal contribution, leading to some debate regards the value of those subtypes. Two dominant transcriptional subtypes have emerged: a Classical subtype that tends to be more responsive to chemotherapy and a very aggressive poorly differentiated Squamous/Basal-like subtype.

Collisson et al. (2011) used micro-dissected tumor samples from resected primary PDAC from two different clinical series to define three specific gene expression subtypes.

  • Exocrine-like: characterized by relatively high expression of tumor cell derived digestive enzyme genes.

  • Classical: demonstrating high expression of adhesion-associated and epithelial genes, and epithelial cell terminal differentiation genes, notably GATA6; KRAS mRNA levels elevated relative to the other subtypes; Classical subtype cell lines are more sensitive to erlotinib.

  • Quasi-mesenchymal: has high expression of mesenchyme associated genes; a relatively high proportion of high-grade tumors and poor patient outcomes; low GATA6 expression; QM-PDA subtype cell lines are relatively more sensitive to gemcitabine than those with the Classical subtype.

Moffitt et al. (2015) used a diverse collection of pancreatic gene expression microarray data, including normal pancreata samples as well as primary and metastatic cancer samples, to identify two tumor-specific subtypes as well as additional stromal Normal and Activated subtypes which were independently prognostic. To develop their two tumor-specific subtypes, Moffitt et al. excluded transcripts thought to be specifically enriched in either the normal pancreas or the tumor microenvironment. The two tumor-specific subtypes were referred to as Classical and Basal-Like (Moffitt et al., 2015).

  • Classical: characterized by overlapping signature with the genes described in the Collisson classification including GATA6, and overall a better prognosis.

  • Basal-like: associated with a worse prognosis than the Classical subtype but may have a better response to adjuvant therapy.

Bailey et al. (2016) described four subtypes using samples with >40% cellularity from resectable primary pancreatic cancer, based on differential transcription factor expression and downstream targets responsible for lineage specification and differentiation during development and regeneration.

  • Squamous: is characterized by enrichment for TP53 and KDM6A mutations; upregulation of the TP63ΔN transcriptional network; hypermethylation of pancreatic endodermal cell-fate determining genes and is associated with a poor clinical prognosis.

  • Pancreatic progenitor: is defined by preferential expression of genes involved in early pancreatic development notably FOXA2/3, PDX1, and MNX1 and also by gene programs involved in metabolism.

  • Immunogenic: is classed by the enrichment of genes associated with specific immune cell populations, including T-cells and B-cells.

  • Aberrantly differentiated endocrine exocrine (ADEX): is featured by upregulation of genes that regulate networks involved in KRAS activation, and exocrine (NR5A2 and RBPJL) and endocrine differentiation (NEUROD1 and NKX2-2).

The Squamous subtype overlaps with the Quasi-mesenchymal subtype of Collisson but has notable pan-squamous features, including a significant association with adenosquamous PDAC histology (Bailey et al., 2016). There is a marked epigenetic shift, with changes in DNA methylation down-regulating key transcription factors controlling pancreatic cell fate determination (PDX1, MNX1, GATA6, HNF1B), and the activation of subtype-driver multigene programs regulated by ΔNTP63 and c-MYC, leading to a loss of endodermal identity (Bailey et al., 2016). In addition, the Squamous subtype was also found to be enriched for mutations in KDM6A, MLL2, and MLL3 chromatin modifying enzymes that belong to the COMPASS complex (COMplex of Proteins Associated with Set1-like) (Bailey et al., 2016).

The Pancreatic progenitor subtype has four key characteristics.

  • (i) Transcriptional networks containing transcription factors PDX1, MNX1, HNF4G, HNF4A, HNF1B, HNF1A, FOXA2, FOXA3, and HES1, which are pivotal for pancreatic endoderm cell-fate determination toward a pancreatic lineage and are linked to maturity onset diabetes of the young

  • (ii) Gene programs regulating metabolism notably fatty acid oxidation, steroid hormone biosynthesis, and drug metabolism

  • (iii) O-linked glycosylation of mucins, notably apomucins MUC5AC and MUC1, but not MUC2 or MUC6, that define the IPMN pancreatobiliary subtype with PDAC-associated IPMN clustering

  • (iv) TGFBR2 inactivating mutations.

The ADEX subtype was defined by both exocrine and endocrine lineage features in later stages of pancreatic development and differentiation (rather than one or the other as is in normal pancreas development), and could be considered a subclass of Pancreatic progenitor tumors. There are two main transcriptional networks.

  • (i) Acinar cell differentiation and pancreatitis/regeneration, transcription factors NR5A2, MIST1 (BHLHA15A) and RBPJL and their downstream targets.

  • (ii) Endocrine differentiation and maturity onset diabetes of the young, including INS, NEUROD1, NKX2-2, and MAFA.

Puleo et al. (2018) proposed two classifications, one specifically for the transformed neoplastic tumor cells and the other for the complete tumor entity, including the stroma: pure basal-like, stroma activated, desmoplastic, pure classical, and immune classical:

  • Pure basal-like tumors are composed of poorly differentiated tumors with predominant Gly12Asp and Gly12Val KRAS mutations; they have a low stromal signal.

  • Stroma Activated tumors are moderately differentiated, specifically enriched in the activated stroma component defined by high a-SMA, SPARC, and FAP.

  • Desmoplastic tumors are also moderately differentiated with a predominant basal association, characterized by a low tumoral component and a large stromal transcriptomic signal, including immune and inflammatory stroma components and, particularly, a high expression of structural and vascularized stroma components.

  • Pure-classical and Immune classical tumors are histologically well differentiated with fewer CDKN2A and TP53 mutations than basal-like tumors, and are also enriched with the Gly12Arg KRAS mutation, and associated with hENT1 expression; predicted to be Moffitt–Classical, and Bailey–Progenitor subtypes.

Maurer et al. used laser capture microdissection (LCM) epithelial cell enriched samples for mRNA sequencing to profile the expression of 60 matched pairs of human PDAC malignant epithelial and stroma samples (Maurer et al., 2019). They developed a computational model that could infer tissue composition and generate virtual compartment-specific expression profiles from bulk gene expression cohorts (Maurer et al., 2019). This study was able to provide a clearer understanding on the previous molecular gene signatures built from bulk tumor tissue samples with the following conclusions.

  • (1)

    Genes used to define the Collisson–Classical, Moffitt–Classical, Moffitt–Basal-like, and Bailey–Progenitor subtypes predominantly provide information about the malignant compartment regardless of the amount of stromal cell infiltration.

  • (2)

    Genes used to define the Moffitt–Activated, Moffitt–Normal, and Bailey–Immunogenic subtypes report on stromal expression that is largely independent of the malignant compartment (but see below).

  • (3)

    Gene sets in the Collisson–Quasi-Mesenchymal and Bailey–Squamous subtypes represent a mixture of epithelial and stromal identity, indicative of a more poorly differentiated state.

  • (4)

    Most genes that define the Collisson–Exocrine and Bailey–ADEX subtypes are largely derived from bulk tumor tissue samples are arguably mostly absent from LCM samples (but see below).

Chan-Seng-Yue et al. (2020) used LCM-purified pancreatic cancers for whole-genome sequencing in tumors from 314 patients, and whole-transcriptome sequencing of tumors from 248 patients, accompanied by single-cell RNA sequencing on 13 resectable and two metastatic tumors. For this classification, tumors with homologous recombination defects and MMR deficiency were excluded due to their unique mutational signatures (Chan-Seng-Yue et al., 2020).

  • Basal-like A: these tumors were associated with the epithelial mesenchymal transition (EMT) program; TP53 gene and TGF-β signaling enriched; 5% of stage I/II (resectable); and 24% of stage IV (metastatic) tumors.

  • Basal-like B: these tumors were associated with the EMT program; TP53 and TGF-β signaling enriched; 9% of stage I/II (resectable), 7% of stage III (locally advanced), and 12% of stage IV (metastatic) tumors.

  • Hybrid: this subtype was found in 24% of stage I/II (resectable), 43% of stage III (locally advanced), and 18% of stage IV (metastatic) tumors.

  • Classical A: this subtype was found in 44% of stage I/II (resectable), 43% of stage III (locally advanced), and 36% of stage IV (metastatic) tumors.

  • Classical B: this subtype was found in 8% of stage I/II (resectable), 7% of stage III (locally advanced), and 10% of stage IV (metastatic) tumors.

This classification split each of the previously defined Basal-like and Classical subtypes into two disease subtypes, while the Hybrid subtype was inconsistently classified by previous systems arising from multiple expression profiles (Chan-Seng-Yue et al., 2020). Single-cell RNA sequencing revealed that both Basal-like and Classical clusters were present in the same tumor found in 13 out of 15 patients (Chan-Seng-Yue et al., 2020). The EMT program was positively correlated with Basal-like signatures and negatively correlated with Classical signatures (Chan-Seng-Yue et al., 2020). Moreover, they found that a major imbalance of allelic states of KRAS (KRASMa) favoring the mutant allele over the wild-type allele occurred in only 4% of primary tumors compared to 29% in metastatic disease (Chan-Seng-Yue et al., 2020). Basal-like A/B tumors were enriched for the major imbalance KRASMa allelic states (44%) compared to metastatic Classical A/B tumors (14%), and KRASMa tumors were also more chemoresistant (Chan-Seng-Yue et al., 2020).

They proposed a possible model for the genomic evolution of pancreatic cancer as being a consequence of a gene expression continuum from (a) both Basal-like and Classical cell populations, and (b) linked to allelic imbalances in mutant (mt) KRAS, with metastatic tumors being more copy number-unstable compared to primary tumors (Chan-Seng-Yue et al., 2020). In primary tumors, the Basal-like phenotype is linked to minor mtKRAS allelic imbalances, whist in metastatic tumors, it is linked to major mtKRAS allelic imbalances (Chan-Seng-Yue et al., 2020).

Potential Influence of Tumor Cellularity on Transcriptomic Subtypes

Raphael and Cancer Genome Atlas Research Network (2017) performed genomic, transcriptomic, and proteomic profiling of 150 PDAC specimens, including samples with low neoplastic cellularity, provided by the Cancer Genome Atlas Research Network. They applied clustering techniques to reproduce the four-group classification of Bailey et al. (2016; Squamous, Immunogenic, Pancreatic Progenitor, and ADEX), the three-group classification (Classical, Quasi-mesenchymal, and Exocrine-like) of Collisson et al. (2011) and the two-group classification (Basal-like or Classical) of Moffitt et al. (2015). They found that while the Basal-like and Classical subtypes were independent of cancer cell purity, the Collisson Exocrine-like and Quasi-Mesenchymal subtypes, and the Bailey ADEX and Immunogenic subtypes were all associated with lower tumor purity (Raphael and Cancer Genome Atlas Research Network, 2017). Raphael and Cancer Genome Atlas Research Network (2017) also found that, among low purity tumors, a higher estimated leukocyte fraction was associated with the Immunogenic subtype and that the ADEX subtype was a subset of the Collisson Exocrine-like subtype.

Puleo et al. (2018) using formalin-fixed and paraffin-embedded tissues also concluded that the ADEX tumor subtype largely resulted from contamination with pancreatic acinar cells. As few as 39 of the most highly expressed genes of normal acinar cells from healthy pancreas single-cell transcriptomes can alone constitute 50% of the total number of expressed transcripts such that even a low level of normal pancreas contamination can materially affect any otherwise presumed subtype (Puleo et al., 2018).

Maurer et al. (2019) also suggested that the Collisson–Exocrine and Bailey–ADEX subtypes might be a function of the degree of tumor cellularity rather than being a distinct subtype as most of the subtype defining genes are largely derived from bulk tumor tissue samples and are mostly absent from LCM epithelial cell enriched samples.

Nevertheless, the assertion that the Collisson–Exocrine and Bailey ADEX and Immunogenic subtypes were all associated with lower tumor purity cannot be entirely true, since the same gene expression signatures seen in patient clinical PDAC tumors are identified in derived cell lines—cell lines and xenografts from these same tumors, and specifically the Classical, Quasi-mesenchymal and Exocrine-like gene expression profiles (Jones et al., 2008; Knudsen et al., 2018). Moreover, most clinical PDAC tumors have low cellularity, so these too should be included to avoid observer bias. The Bailey Immunogenic subtype as well as containing gene expression profiles derived from tumor stroma immune infiltration predominantly related to B and T cells also contains an underlying Pancreatic progenitor-like gene expression character (Bailey et al., 2016). Both cytotoxic (CD8+) and regulatory T cells (CD4+CD25+FOXP3+ Tregs) are predominant (Bailey et al., 2016). It has been suggested that a distinct Immunogenic subtype does not exist as distinct since immune infiltrates are enriched across all tumor-intrinsic subtypes, and their prevalence is primarily driven by tumor cellularity of the sequenced samples (International Cancer Genome Consortium et al., 2010; Puleo et al., 2018; Maurer et al., 2019). Nevertheless, by allowing for different degrees of cellularity, the strong Pancreatic Progenitor-like signals can still be split into an immune high signature (Immunogenic subtype) and an immune low signature (Progenitor subtype) indicating that the signals from the underlying epithelium drive and formulate tumor cell immunogenicity (Figure 1; Bailey et al., 2016; Bear et al., 2020; Sahai et al., 2020). Hwang et al. using single-nucleus RNA sequencing found an association between Basal-like programs and higher immune infiltration with increased lymphocytic content, whereas Classical-like programs were associated with sparser macrophage-predominant microniches (Hwang et al., 2020).

FIGURE 1

FIGURE 1

Epithelial and stromal cell interactions in pancreatic cancer (Jones et al., 2008; Waddell et al., 2015; Bailey et al., 2016; Bear et al., 2020; Dominguez et al., 2020; Escobar-Hoyos et al., 2020; Sahai et al., 2020). PDAC immune resistance is driven by complex genetic background. Expression of tumor-intrinsic GM-CSF and CXCL1 is increased by oncogenic KRAS to mediate T cell exclusion and MDSC infiltration. Downstream signaling initiated by mutant KRAS (mKRAS) mediates innate and adaptive immune escape through enhancing autophagy to downregulate MHC-1 expression and upregulate the expression of PD-L1 and CD47. In addition to increased IL-6-mediated systemic dysregulation of conventional type 1 dendritic cell (DC), activation of WNT/β-catenin mediated by mKRAS signaling further downregulates CCL4 expression to inhibit DC recruitment. Tumor group 2 innate lymphoid cells (TILC2s) infiltrate the tumor microenvironment and are activated by IL-33 through binding to the ST2 receptor, further leading to an enhancement of anti-tumor immunity by expressing the inhibitory checkpoint receptor PD-1 and recruiting DCs potentially through CCL5 production. Furthermore, mKRAS signaling enhances chronic inflammation signaling such as Sonic Hedgehog, COX2, and pSTAT3 signaling, and promotes multiple inflammation-associated factors such as IL-1, IL-6, tumor necrosis factor (TNF), and matrix metalloproteinase 7 (MMP) to activate cancer-associated fibroblasts (CAF). Additional factors leading to activation of CAFs include TGF-β, extracellular matrix (ECM) stiffness and composition, RTK ligands such as PDGF and FGF, DNA damage caused by chemotherapy and radiotherapy, physiological stress, and contact signals such as Notch and Eph-ephrins. Activated CAFs further regulate macrophage and endothelial functions by factors such as VEGF, HGF, and GAS6 and participate in immune crosstalk through TGF-β activation, IL-6, CXCL12, and CCL2 production. Deficiency of p53 mediates transition of TAM toward an immunosuppressive M2 phenotype. Mutant p53 (such as R175H) increases expression of the splicing regulator hnRNPK to promote inclusion of cytosine-rich exons (+polyC exons) within GTPase-activating proteins (GAPs), particularly GAP17, leading to enhanced KRAS activity. CCL2/4/5, CC-chemokine ligand 2/4/5; CCR, CC-chemokine receptor; COX2, cyclooxygenase 2; CSF-1, colony-stimulating factor 1; CSF-1R, colony stimulating factor 1 receptor; CTLA-4, cytotoxic T-lymphocyte-associated protein 4; CXCL1/12, CXC-chemokine ligand 1/12; CXCR4, CXC-chemokine receptor type 4; DC, conventional type 1 dendritic cell; FGF, fibroblast growth factor; Flt3L, Fms related receptor tyrosine kinase 3 ligand; GAS6, growth arrest-specific protein 6; GM-CSF, granulocyte-macrophage colony-stimulating factor; hnRNPK, heterogeneous nuclear ribonucleoprotein K; HA, hyaluronic acid; HGF, hepatocyte growth factor; IL-1/-6/-33, interleukin-1/-6/-33; MHC-1, major histocompatibility complex 1; PDGF, platelet-derived growth factor; PD-1, programmed cell death protein 1; PD-L1, programmed death-ligand 1; SHH, sonic hedgehog; ST2, suppression of tumorigenicity 2; STAT3, signal transducer and activator of transcription 3; TAM, tumor-associated macrophage; TGF-β, transforming growth factor-β; VEGF, vascular endothelial growth factor; VISTA, V-domain Ig suppressor of T cell activation.

At the present time, it is not entirely clear that the Collisson Exocrine-like and Quasi-Mesenchymal subtypes, and the Bailey ADEX and Immunogenic subtypes should be discarded, as there is a considerable variation in the way samples have been retrieved, stored, analyzed for mRNA expression, and assessed for epithelial cell purity by direct and indirect methodologies (Table 2; International Cancer Genome Consortium et al., 2010; Collisson et al., 2011; Carter et al., 2012; Moffitt et al., 2015; Bailey et al., 2016; Williamson et al., 2016; Cancer Genome Atlas Research Network, 2017; Connor et al., 2017; Raphael and Cancer Genome Atlas Research Network, 2017; Knudsen et al., 2018; Mateo et al., 2018; Puleo et al., 2018; Wartenberg et al., 2018; Maurer et al., 2019; Bear et al., 2020; Chan-Seng-Yue et al., 2020; Dijk et al., 2020; Escobar-Hoyos et al., 2020; Hong et al., 2020; Hwang et al., 2020). The validity of the Collisson Exocrine-like and Quasi-mesenchymal subtypes, and the Bailey ADEX and immunogenic subtypes requires further investigation.

Commonality of Transcriptomic Signatures

There is a strong alignment between the Classical/Pancreatic Progenitor and Quasi-mesenchymal/Basal-like/Squamous subtypes signatures of Moffit, Collinson, Bailey, Puleo and Chan-Seng-Yue (Figure 2; Moffitt et al., 2015; Bailey et al., 2016; Cancer Genome Atlas Research Network, 2017; Puleo et al., 2018; Chan-Seng-Yue et al., 2020).

FIGURE 2

FIGURE 2

Comparison of different transcriptional classifications of PDAC. Comparison of previously published transcriptional classifications of PDAC, two major consensus subtypes have been identified (Chan-Seng-Yue et al., 2020): (A) Consensus Classical, which is named as “Classical” in the classifications of Collisson et al. (2011) and Moffitt et al. (2015), “Progenitor” by Bailey et al. (2016), “Pure Classical” by Puleo et al. (2018), and “Classical-A/-B” by Chan-Seng-Yue et al. (2020); (B) Consensus Basal, which is named as “Basal-like” by Moffitt et al. (2015), “Quasi-Mesenchymal” by Collisson et al. (2011), “Squamous” by Bailey et al. (2016), “Pure Basal-like” by Puleo et al. (2018), and “Basal-like A/B” by Chan-Seng-Yue et al. (2020).

Stromal Immune Cell and Cancer-Associated Fibroblast Infiltrate

Signals from the stroma play an important role in disease progression (Bailey et al., 2016; Ligorio et al., 2019; Bear et al., 2020; Sahai et al., 2020). PDAC is characterized by a complex and dense microenvironment with an extensive desmoplastic stromal reaction. Typically, around 5–30% of cells in pancreatic tumors are epithelial cancer cells. Activation of pancreatic stellate cells and cancer-associated fibroblasts (CAF), along with inflammatory and immune cell accumulation, occurs during early pancreatic tumorigenesis, creating an immunosuppressive microenvironment that restricts immune surveillance and supports tumor growth and invasiveness (Figure 1; Bear et al., 2020; Sahai et al., 2020). Oncogenic driver mutations promote immunosuppression from the earliest stages of tumor inception that accompanies oncogenesis. Beyond immunogenic prognostic subtypes, patient-specific immune changes should be considered in combination immune-modulatory therapies targeting roadblocks in antitumor immunity. An immune-signature-based stratification may guide personalized therapy of PDAC patients and enable the design of novel combinatorial treatments with improved clinical efficacy (Kandimalla et al., 2020).

Pancreatic cancer has relatively few coding mutations and thus only few neo-antigenic targets, and is embedded in an immunosuppressive cold tumor microenvironment, which impedes intratumoral CD8+ T cell infiltration and activation (Bear et al., 2020). Therefore, endogenous PDAC-reactive T cells are limited in quantity and quality and single agent immunotherapies with immune checkpoint inhibitors, which unleash pre-existing T cell immunity, are mostly ineffective in pancreatic cancer (Bear et al., 2020). Yet, exceptionally high neoantigen numbers, with robust antitumor CD8+ T cell responses have been associated with long-term survival in pancreatic cancer patients, and immune checkpoint blockade has shown clinical responses in patients with hypermutated MMR-deficient tumors (Balachandran et al., 2017; Le et al., 2017). To induce specific antitumor adaptive immune responses, tumor-derived antigens must not only be taken up by innate immune cells; they must also be efficiently processed and cross-presented to CD8+ T cells in the presence of a costimulatory signal. This mechanism is impeded in the vast majority of PDAC tumors by immunosuppressive mechanisms. CD8+ lymphocytes are trapped in peritumoral compartments and mostly display an exhausted gene expression profile (Steele et al., 2020). Th1-polarized CD4 + T cells are less frequent at the tumor site compared to Th2-polarized CD4+ T cells (De Monte et al., 2011). Tumor-associated macrophages (TAM) and myeloid-derived suppressor cells (MDSCs), which thwart the generation of cytotoxic T cell responses, are predominant over dendritic cells (DCs) that are largely dysfunctional, and immune checkpoint ligands are upregulated on myeloid cells (Hegde et al., 2020; Steele et al., 2020). In order to overcome the web of immune resistance and achieve durable antitumor effects, immunotherapeutic regimens need to target different steps in the cancer-immunity cycle, combining ideal antigen presenting cell (APC) activation that mediates priming of tumor-specific T cells, with strategies that enhance T cell effector function, and disrupt immunosuppressive myeloid cell programs. Thus, immune-modulatory strategies must be multi-modal aiming to (1) enhance endogenous T-cell function, (2) adoptively transfer tumor-specific T-cell immunity, and (3) attempt to devise an immunologically hot tumor microenvironment (Bear et al., 2020).

Moral et al. have shown that group 2 innate lymphoid cells (ILC2s) infiltrate PDACs to activate tissue-specific tumor immunity, inferring another novel immunoregulatory target (Moral et al., 2020). Enhanced anti-tumor immunity ensued blockade of the T cell checkpoint receptor programmed death (PD) receptor-1, which released ILC2 cell-intrinsic inhibition to expand and activate the tumor ILC2s to produce CCL5, thereby resulting in CD103+ dendritic cell expansion and then CD8 + T-cell activation (Moral et al., 2020). Tumor infiltrating ILC2s which express the programmed cell death protein (PD-1) receptor were enriched in long-term survivors with an immunologically hot tumor microenvironment containing abundant activated CD8+ T-cells, and containing higher bulk tumor RNA expression of the ILC2-activating cytokine IL33 (Moral et al., 2020). Pre-clinical studies in KPC mouse models have also suggested that specific targeting of macrophages and neutrophils using small molecule inhibitors, specific for either macrophage receptor CSF1R or the neutrophil receptor CXCR2, might facilitate better outcomes by enhancing endogenous T-cell cancer killing functions and the reprogramming of tumor cell intrinsic phenotypes (Steele et al., 2016; Candido et al., 2018). The PD-1 antibody pembrolizumab has FDA approval to for the treatment of MSI-H solid tumors, although this is present in only 1–3% of pancreatic cancers (Marcus et al., 2019; Marabelle et al., 2020). Inhibition of the CXCR4–CXCL12 pathway in pancreatic cancer also enhances tumor sensitivity to anti-PD-1 ligand-1 treatment. In the two-cohort phase IIa, COMBAT study (NCT02826486) pembrolizumab was combined with BL-8040 (a CXCR4 antagonist) in metastatic pancreatic cancer with promising responses and survival rates (Bockorny et al., 2020).

Taking Molecular Subtyping Into Clinical Trials

Targeted therapies for advance pancreatic cancer based on next generation sequencing has been disappointing, with only 3–5% showing any clinical benefit in terms of actionable mutations, and limited to only a few months of additional survival (Pishvaian et al., 2020; Cobain et al., 2021). In reality the greatest sensitivity of pancreatic cancer to systemic therapies is chemotherapy, with increasing interest being shown in developing treatment response transcriptomic signatures to different agents (Tiriac et al., 2018). Moreover, it is important to distinguish treatment response signatures from the two main molecular subtypes as the Basal-like subtype appears to be more chemoresistant compared with reduced patient survival to the Classical subtype (Bailey et al., 2016; Aung et al., 2018; Chan-Seng-Yue et al., 2020; O’Kane et al., 2020). Several groups have now established informatic approaches that proport to accurately stratify patients based on PDAC subtype for clinical use. These include PuRIST, a single sample classifier that can stratify patients into two tumor-cell intrinsic subtypes based on the Moffitt classification scheme (Rashid et al., 2020). While PuRIST and other similar approaches promise better patient selection for chemotherapy, they have not been assessed in clinical trials.

Noll et al. identified hepatocyte nuclear factor (HNF)-1A and KRT81 that enabled stratification of tumors into different molecular subtypes by using immunohistochemistry (Noll et al., 2016). The two-marker combination identified the QM-PDA (KRT81+/HNF1A) subtype, which was associated with the shortest survival; the Exocrine-like (KRT81/HNF1A+) subtype which was associated with the longest survival; and the Classical (KRT81/HNF1A) subtype, which was associated with intermediate survival (Noll et al., 2016). Exocrine-like subtype tumors were resistant to tyrosine kinase inhibitors (erlotinib and dasatinib) and paclitaxel, which induced cytochrome P450 (CYP) 3A5 (CYP3A5) in the tumors, leading to the metabolism of these compounds (Noll et al., 2016). CYP3A5 expression was correlated positively with HNF1A+ and negatively with KRT81, and also contributed to acquired resistance in the QM-PDA and Classical subtypes (Noll et al., 2016).

Hwang et al. (2020) performed single-cell RNA sequencing on 26 flash-frozen pancreatic cancers from patients who underwent surgical resection, with upfront surgery in 11 and in 15 after neoadjuvant chemoradiotherapy. Following chemoradiation, there was a relative increase in Basal-like cells (including the master transcription factor ΔTP63 for the Squamous subtype), and a decrease in Classical-like cells (including the hallmark transcription factor GATA6) (Hwang et al., 2020). Thus, there appears to be a commonality with the effects of chemotherapy such as FOLFIRINOX which will also enrich for the Basal-like subtype (Aung et al., 2018; Porter et al., 2019). Following chemoradiotherapy, there was enhanced expression of genes needed to maintain the Wnt/β-catenin niche, which is critical for treatment resistance and can be mediated by autocrine signaling of the epithelial cells and/or paracrine interactions with CAFs (Hwang et al., 2020). Squamous/Basal-like programs facilitate immune infiltration compared with the Classical-like programs (Hwang et al., 2020; Somerville et al., 2020). Importantly, the immune infiltrates associated with Basal-like and Classical-like malignant cells are distinct, pointing to differential strategies choosing checkpoint inhibitors for the Basal-like subtype, and for the Classical-like subtype choosing myeloid directed therapies such as CD40 agonists and TGF-β modulators (Hwang et al., 2020). The study by Hwang et al. (2020) has yet to be published following review and the findings and conclusions will need further evaluation. Other approaches to subtyping may be required to understand more fully the extent of interpatient heterogeneity such as differential DNA methylation, associated with interferon (IFN) signaling (Espinet et al., 2021).

The encouragement from COMBAT along with other immunotherapy approaches currently being tested will expand the armamentarium against pancreatic cancer. A detailed understanding of the individual patient’s response to the different forms of treatment will be necessary to further improve the prognosis of pancreatic cancer patients. The Bailey Immunogenic subtype is associated with immune gene programs involving B-cell signaling pathways, antigen presentation, CD4+ T-cell, CD8+ T-cell, and Toll-like receptor signaling pathways (Bailey et al., 2016). Acquired tumor immune suppression pathways through upregulation of the T cell checkpoint receptor PD-1 and cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) in this Immunogenic subtype may offer therapeutic opportunities (Bailey et al., 2016). Puleo et al. (2018) found that the expression of CTLA-4 was higher in the Immune classical and Desmoplastic subtypes and, to a lesser extent, in the Pure basal-like subtypes, making these subtypes potentially sensitive to anti-CTLA-4 therapy such as ipilumumab. Also other promising therapeutic targets identified were the inhibitory checkpoint membrane receptors CD276 (B7-H3) and HAVCR2 (TIM3), both of which were highly expressed in the Desmoplastic, Stroma activated, and Pure basal-like subtypes (Wartenberg et al., 2018). Immune classical and Desmoplastic subtypes also showed high expression of the T-cell checkpoint inhibitor receptor CTLA-4, the costimulatory T-cell receptor CD27, and the tumor inhibitory T4+ cell subset CD26 marker protein; Basal-like tumors were enriched in CD276 (B7-H3) and HAVCR2 (TIM3); and PDL-2 (PDCD1LG2) was expressed in all but the Pure classical subtypes, which overall do not up-regulate any of the immune checkpoints (Puleo et al., 2018).

In an effort to improve outcomes in pancreatic cancer through the use of more effective therapeutics, large-scale efforts are required with multiple centers and cooperating disciplines. Multiple clinical trials have been launched, pursuing better treatment schemes and ideal medication regimens based on molecular profiling and subtyping including COMPASS, PREDICT-PACA, Precision Promise, Know Your Tumor, ESPAC6/7, PANCuRx, and PASS1 (Table 3).

TABLE 3

Trial nameYearType of studyCountryTumor stageMethodologyTreatmentAims/Results
COMPASS (Aung et al., 2018; O’Kane et al., 2020) NCT027506572018Cohort studyCanadaLocally advanced or metastatic PDACRNA-seq156 Classical vs. 39 Basal-like subtypes, 91 mFFX vs. 63 GnPORR: 33% vs. 10%; Classical mFFX 29.6%, Basal-like mFFX 10%, Classical GA 39%, Basal-like GA 11%. mOS: 9.3 vs. 5.9 mths; Classical mFFX 10.6 mths, Basal-like mFFX 6.5 mths, Classical GA 8.19 mths, Basal-like GA 8.12 mths.
ESPAC-6 EudraCT: 2020-004906-792021Phase IIIEuropeResectable PDACLCM, organoid generation, RNA-seq, WGS/WESAdjuvant chemotherapy, mFFX vs. Gemcitabine with capecitabinePatients randomized to allocation of oxaliplatin- or gemcitabine-based chemotherapy by standard clinical criteria or by a transcriptomic treatment specific stratification signature after surgery.
(1) Primary outcome: disease free survival, time from randomization to disease recurrence or death from any cause
(2) Secondary outcomes: overall survival; survival based on targeted signatures in test versus control arms; response to targeted therapies
ESPAC-72022Randomized phase IIEuropeNon-metastatic locally advanced PDACOrganoid generation, RNA-seq, WGS/WESNeoadjuvant chemotherapy, mFFX vs. GnPPatients randomized to allocation of oxaliplatin- or gemcitabine-based chemotherapy by standard clinical criteria or by a transcriptomic treatment specific stratification signature before surgery.
(1) Primary outcome: resection rate following neoadjuvant therapy Disease free survival, time from randomization to disease recurrence or death from any cause
(2) Secondary outcomes: resection rates and disease-free survival based on targeted signatures in test versus control arms; overall survival; response to targeted therapies
PASS-01 NCT044695562019Randomized phase IIUnited States, CanadaMetastatic PDACOrganoid generation, RNA-seqmFFX vs. GA(1) To determine the PFS benefit, ORR, DOR, and OS of mFFX compared to GA in patients with metastatic PDAC.
(2) To evaluate chemotherapy sensitive signatures.
(3) To determine if organoid transcriptomic profiles parallel patient transcriptomic profiles.
(4) To explore biomarkers as responders to mFFX/GA (including GATA6 validation).
Know Your Tumor (Pishvaian et al., 2020)2020Registry studyUSAPancreatic cancerGenomic testing and tailored treatment recommendation based on the molecular profiling.Molecularly matched therapy group vs. unmatched therapy group vs. no molecular marker groupmOS 2.58 vs. 1.51 vs. 1.32 y; patients with 2/more lines of therapy: mOS 1.81 vs. 0.85 vs. 0.73 y, mPFS 10.93 vs. 4.53 vs. 5.37 mths;
Precision Promise NCT042290042020Phase II/III trialsUnited StatesMetastatic PDACGenomic, transcriptomic and immune analysisSM-88 vs. SOC (GA/mFFX)(1) To compare each investigational arm in OS in 1st/2nd line metastatic pancreatic cancer patients and to determine which patients benefit from which regimen.
(2) To determine short-/long-term safety signals of each investigational arm.
(3) To determine clinical benefit (PFS, OS, CR, and PR) and the duration of it.
PREDICT-PACA DRKS000224292021Co-clinical validation cohort studyGermanyMetastatic PDACOrganoid generation/chemosensitivity testing, PDX analyses, NGS, transcriptomics,Gemcitabine monotherapy, GnP, 5-FU with oxaliplatin(1) To validate molecular signatures predicting response to specific chemotherapy regimens and establish point-of-care diagnostics.
(2) To establish PDO-based chemosensitivity testing for individual tumors.
(3) To evaluate the influence of the tumor microenvironment’s composition on chemotherapy response.
PANCuRx2021Cohort studiesCanadaPDACWGS, RNA-seqTranslational Research Initiative(1) To seek solutions to the high fatality rate of PDAC by generating genetic and biologic knowledge of the cancer, the mechanisms of how tumors grow and tailored treatment options.
(2) To expand the COMPASS trial into the other cancer centers across Canada.

Translational clinical trials with molecular subtyping.

COMPASS, Study of Changes and Characteristics of Genes in Patients with Pancreatic Cancer for Better Treatment Selection; ESPAC, European Study Group for Pancreatic Cancer; PASS-01, Pancreatic Adenocarcinoma Signature Stratification for Treatment; PREDICT-PACA, Integrated Biopsy-Based Approach to Predict Response to Chemotherapy for Patients with Stage IV Pancreatic Cancer; NGS, next-generation sequencing; mFFX, modified FOLFIRINOX; GnP, gemcitabine plus nab-paclitaxel; GA, gemcitabine-abraxane; SOC, standard of care; mths, months; y, year; ORR, overall response rate; (m)OS, (median) overall survival; PFS, progression-free survival; DOR, duration of response; CR, complete response; PR, partial response; PDO, patient-derived organoid.

The EPPIC (Enhanced Pancreatic Cancer Profiling For Individualized Care Study) study based in Canada aims to sequence metastatic pancreatic tumors of 400 patients through two clinical trials (COMPASS and PanGen), both of which are generating molecular and phenotypic signatures of individual tumors in a clinically relevant timeframe and related to chemotherapy responses1 [Aung et al., 2018; O’Kane et al., 2020].

Precision Promise was created by the Pancreatic Cancer Action Network and 15 USA clinical academic sites, in cooperation with the FDA and pharmaceutical partners. Precision Promise is an active adaptive Phase II/III clinical trial platform (NCT04229004) that allows rapid evaluation of novel therapeutic options in patients with metastatic pancreatic cancer. The protocol utilizes an adaptive randomization design and includes several trial designs and statistical innovations. All patients undergo pre- and on-treatment biopsies with state-of-the-art genomic, transcriptomic, and immune analysis, along with collection of blood research samples throughout the study. Focused on both first- and second-line treatment of metastatic PDAC, 30% of patients are randomized to one of the two common control arms (gemcitabine plus nab-paclitaxel and FOLFIRINOX), while 70% of patients are randomized to an experimental treatment arm. The platform currently has one experimental arm open (SM-88, Tyme), with two additional experimental arms to be added in 2021. Compared to traditional trial designs, Precision Promise has several advantages: multiple investigational treatments can be evaluated simultaneously using common controls; only 175 patients per experimental arm are required to initiate a regulatory registration, and it is expected that this platform will significantly accelerate the time to evaluate a new therapy, with an anticipated cost savings of 30–50%. With its unique design and novel method of data sharing, Precision Promise serves as a new clinical trial ecosystem to accelerate drug development for PDAC.

The ESPAC-6 adjuvant trial in resectable patients and the ESPAC-7 neoadjuvant trial in locally advanced patients are evaluating oxaliplatin- or gemcitabine-based chemotherapy response of PDAC patients that will be randomized according to standard clinical criteria (control arms) or by transcriptomic stratification signatures (experimental arms).

The PREDICT-PACA co-clinical trial funded by the German Cancer Aid is being conducted as a biopsy-based approach to predict response to chemotherapy for patients with metastasized pancreatic cancer. The consortium has established robust and highly predictive transcriptomic signatures linked to specific chemotherapy response profiles. Using microfluidic card-based qRT-PCR marker panels, clinical utility of the signatures is validated in a prospective collection of core biopsies from metastases of stage IV PDAC patients that receive one of the current chemotherapies upon clinical decision. In parallel, potential alternative therapies for third-line treatment are identified by high-throughput drug screening using patient-derived organoids and by next-generation sequencing-based detection of actionable mutations.

A major area of discussion is whether to discard the tumor microenvironment for molecular classification and replace this with a tumor-cell intrinsic classification system, but this may seem short-sighted as the stroma is a uniquely powerful biological phenomenon in pancreatic cancer. The relationship between molecular subtypes, most notably the immunogenic subtype as described by Bailey et al. (2016), and the responsiveness to evolving multimodality and immunotherapeutic strategies needs further investigation. Single-cell and spatial transcriptomic approaches now allow single-cell profiling of tumor and immune cell populations resident in patient tumors (Kuboki et al., 2019; Hwang et al., 2020; Moncada et al., 2020). These approaches are providing unparalleled insights into immune-tumor interactions and offer new opportunities for targeted immunotherapeutic intervention.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Author contributions

XZ, KH, PB, CS, SR, TG, MB, and JN prepared the initial drafts of the manuscript which were then critically reviewed by RK, BB, JA, KW, TP, MWB, and TH. All of the authors approved the final version of the manuscript.

Acknowledgments

We thank Diane Simeone (Department of Surgery, New York University School of Medicine, New York, NY, United States) for providing an update on the Precision Promise.

Conflict of interest

XZ is in receipt of funding from the Chinese Scholarship Council; PB declares grants from Horizon Europe 2020 (Marie Skłodowska Curie Innovative Training Network); SR declares grants from German Cancer Aid grants (70112720 and 70113167); JN declares grants from the Dietmar Hopp Stiftung GmbH, the Stiftung Deutsche Krebshilfe, and the Heidelberger Stiftung Chirurgie. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    AmrutkarM.GladhaugI. P. (2017). Pancreatic cancer chemoresistance to gemcitabine.Cancers9:157. 10.3390/cancers9110157

  • 2

    AungK. L.FischerS. E.DenrocheR. E.JangG. H.DoddA.CreightonS.et al (2018). Genomics-driven precision medicine for advanced pancreatic cancer: early results from the COMPASS trial.Clin. Cancer Res.2413441354. 10.1158/1078-0432.CCR-17-2994

  • 3

    BaileyP.ChangD. K.NonesK.JohnsA. L.PatchA. M.GingrasM. C.et al (2016). Genomic analyses identify molecular subtypes of pancreatic cancer.Nature5314752. 10.1038/nature16965

  • 4

    BalachandranV. P.LukszaM.ZhaoJ. N.MakarovV.MoralJ. A.RemarkR.et al (2017). Identification of unique neoantigen qualities in long-term survivors of pancreatic cancer.Nature551512516. 10.1038/nature24462

  • 5

    BearA. S.VonderheideR. H.O’HaraM. H. (2020). Challenges and opportunities for pancreatic cancer immunotherapy.Cancer Cell38788802. 10.1016/j.ccell.2020.08.004

  • 6

    BockornyB.SemenistyV.MacarullaT.BorazanciE.WolpinB. M.StemmerS. M.et al (2020). BL-8040, a CXCR4 antagonist, in combination with pembrolizumab and chemotherapy for pancreatic cancer: the COMBAT trial.Nat. Med.26878885. 10.1038/s41591-020-0880-x

  • 7

    BurrisH. A.IIIMooreM. J.AndersenJ.GreenM. R.RothenbergM. L.ModianoM. R.et al (1997). Improvements in survival and clinical benefit with gemcitabine as first-line therapy for patients with advanced pancreas cancer: a randomized trial.J. Clin. Oncol.1524032413. 10.1200/JCO.1997.15.6.2403

  • 8

    Cancer Genome Atlas Research Network (2017). Integrated genomic characterization of pancreatic ductal adenocarcinoma.Cancer Cell32185203.e13. 10.1016/j.ccell.2017.07.007

  • 9

    CandidoJ. B.MortonJ. P.BaileyP.CampbellA. D.KarimS. A.JamiesonT.et al (2018). CSF1R+ macrophages sustain pancreatic tumor growth through T cell suppression and maintenance of key gene programs that define the squamous subtype.Cell Rep.2314481460. 10.1016/j.celrep.2018.03.131

  • 10

    CarterS. L.CibulskisK.HelmanE.McKennaA.ShenH.ZackT.et al (2012). Absolute quantification of somatic DNA alterations in human cancer.Nat. Biotechnol.30413421. 10.1038/nbt.2203

  • 11

    Chan-Seng-YueM.KimJ. C.WilsonG. W.NgK.FigueroaE. F.O’KaneG. M.et al (2020). Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution.Nat. Genet.52231240. 10.1038/s41588-019-0566-9

  • 12

    ChauffertB.MornexF.BonnetainF.RougierP.MarietteC.BouchéO.et al (2008). Phase III trial comparing intensive induction chemoradiotherapy (60 Gy, infusional 5-FU and intermittent cisplatin) followed by maintenance gemcitabine with gemcitabine alone for locally advanced unresectable pancreatic cancer. Definitive results of the 2000-01 FFCD/SFRO study.Ann. Oncol.1915921599. 10.1093/annonc/mdn281

  • 13

    CobainE. F.WuY. M.VatsP.ChughR.WordenF.SmithD. C.et al (2021). Assessment of clinical benefit of integrative genomic profiling in advanced solid tumors.JAMA Oncol.7525533. 10.1001/jamaoncol.2020.7987

  • 14

    CollissonE. A.SadanandamA.OlsonP.GibbW. J.TruittM.GuS.et al (2011). Subtypes of pancreatic ductal adenocarcinoma and their differing responses to therapy.Nat. Med.17500503. 10.1038/nm.2344

  • 15

    ConnorA. A.DenrocheR. E.JangG. H.TimmsL.KalimuthuS. N.SelanderI.et al (2017). Association of distinct mutational signatures with correlates of increased immune activity in pancreatic ductal adenocarcinoma.JAMA Oncol.3774783. 10.1001/jamaoncol.2016.3916

  • 16

    ConroyT.DesseigneF.YchouM.BoucheO.GuimbaudR.BecouarnY.et al (2011). FOLFIRINOX versus gemcitabine for metastatic pancreatic cancer.N. Engl. J. Med.36418171825. 10.1056/NEJMoa1011923

  • 17

    ConroyT.HammelP.HebbarM.Ben AbdelghaniM.WeiA. C.RaoulJ. L.et al (2018). FOLFIRINOX or gemcitabine as adjuvant therapy for pancreatic cancer.N. Engl. J. Med.37923952406. 10.1056/NEJMoa1809775

  • 18

    CunninghamD.ChauI.StockenD. D.ValleJ. W.SmithD.StewardW.et al (2009). Phase III randomized comparison of gemcitabine versus gemcitabine plus capecitabine in patients with advanced pancreatic cancer.J. Clin. Oncol.2755135518. 10.1200/JCO.2009.24.2446

  • 19

    DavisW.MakarG.MehtaP.ZhuG. G.SomerR.MorrisonJ.et al (2019). Next-generation sequencing in 305 consecutive patients: clinical outcomes and management changes.J. Oncol. Pract.15e1028e1034. 10.1200/JOP.19.00269

  • 20

    De MonteL.ReniM.TassiE.ClavennaD.PapaI.RecaldeH.et al (2011). Intratumor T helper type 2 cell infiltrate correlates with cancer-associated fibroblast thymic stromal lymphopoietin production and reduced survival in pancreatic cancer.J. Exp. Med.208469478.

  • 21

    DienerM. K.MihaljevicA. L.StrobelO.LoosM.SchmidtT.SchneiderM.et al (2021). Periarterial divestment in pancreatic cancer surgery.Surgery16910191025. 10.1016/j.surg.2020.08.030

  • 22

    DijkF.VeenstraV. L.SoerE. C.DingsM. P. G.ZhaoL.HalfwerkJ. B.et al (2020). Unsupervised class discovery in pancreatic ductal adenocarcinoma reveals cell-intrinsic mesenchymal features and high concordance between existing classification systems.Sci. Rep.10:337. 10.1038/s41598-019-56826-9

  • 23

    DominguezC. X.MüllerS.KeerthivasanS.KoeppenH.HungJ.GierkeS.et al (2020). Single-cell RNA sequencing reveals stromal evolution into LRRC15+ myofibroblasts as a determinant of patient response to cancer immunotherapy.Cancer Discov.10232253. 10.1158/2159-8290.CD-19-0644

  • 24

    Escobar-HoyosL. F.PensonA.KannanR.ChoH.PanC. H.SinghR. K.et al (2020). Altered RNA splicing by mutant p53 activates oncogenic RAS ignaling in pancreatic cancer.Cancer Cell38198211.e8. 10.1016/j.ccell.2020.05.010

  • 25

    EspinetE.GuZ.ImbuschC. D.GieseN. A.BüscherM.SafaviM.et al (2021). Aggressive PDACs show hypomethylation of repetitive elements and the execution of an intrinsic IFN program linked to a ductal cell of origin.Cancer Discov.11638659. 10.1158/2159-8290.CD-20-1202

  • 26

    GellerL. T.Barzily-RokniM.DaninoT.JonasO. H.ShentalN.NejmanD.et al (2017). Potential role of intratumor bacteria in mediating tumor resistance to the chemotherapeutic drug gemcitabine.Science35711561160. 10.1126/science.aah5043

  • 27

    GhanehP.PalmerD. H.CicconiS.HalloranC.PsarelliE. E.RawcliffeC. L.et al (2020). ESPAC-5F: Four-arm, prospective, multicenter, international randomized phase II trial of immediate surgery compared with neoadjuvant gemcitabine plus capecitabine (GEMCAP) or FOLFIRINOX or chemoradiotherapy (CRT) in patients with borderline resectable pancreatic cancer.J. Clin. Oncol.38 (Suppl. 15):4505. 10.1200/JCO.2020.38.15_suppl.4505

  • 28

    GillS.KoY. J.CrippsC.BeaudoinA.Dhesy-ThindS.ZulfiqarM.et al (2016). PANCREOX: a randomized phase III study of fluorouracil/leucovorin with or without oxaliplatin for second-line advanced pancreatic cancer in patients who have received gemcitabine-based chemotherapy.J. Clin. Oncol.3439143920. 10.1200/JCO.2016.68.5776

  • 29

    GolanT.HammelP.ReniM.Van CutsemE.MacarullaT.HallM. J.et al (2019). Maintenance olaparib for germline BRCA-mutated metastatic pancreatic cancer.N. Engl. J. Med.381317327. 10.1056/NEJMoa1903387

  • 30

    GreenhalfW.GhanehP.NeoptolemosJ. P.PalmerD. H.CoxT. F.LambR. F.et al (2014). Pancreatic cancer hENT1 expression and survival from gemcitabine in patients from the ESPAC-3 trial.J. Natl. Cancer Inst.106:djt347. 10.1093/jnci/djt347

  • 31

    HackertT.SachsenmaierM.HinzU.SchneiderL.MichalskiC. W.SpringfeldC.et al (2016). Locally advanced pancreatic cancer: neoadjuvant therapy with Folfirinox results in resectability in 60% of the patients.Ann. Surg.264457463. 10.1097/SLA.0000000000001850

  • 32

    HammelP.HuguetF.van LaethemJ. L.GoldsteinD.GlimeliusB.ArtruP.et al (2016). Effect of chemoradiotherapy vs chemotherapy on survival in patients with locally advanced pancreatic cancer controlled after 4 months of gemcitabine with or without erlotinib: the LAP07 randomized clinical trial.JAMA31518441853. 10.1001/jama.2016.4324

  • 33

    HegdeS.KrisnawanV. E.HerzogB. H.ZuoC.BredenM. A.KnolhoffB. L.et al (2020). Dendritic cell paucity leads to dysfunctional immune surveillance in pancreatic cancer.Cancer Cell37289307.e9. 10.1016/j.ccell.2020.02.008

  • 34

    HongD. S.FakihM. G.StricklerJ. H.DesaiJ.DurmG. A.ShapiroG. I.et al (2020). KRASG12C inhibition with sotorasib in advanced solid tumors.N. Engl. J. Med.38312071217. 10.1056/NEJMoa1917239

  • 35

    HwangW. L.JagadeeshK. A.GuoJ. A.HoffmanH. I.YadollahpourP.MohanR.et al (2020). Single-nucleus and spatial transcriptomics of archival pancreatic cancer reveals multi-compartment reprogramming after neoadjuvant treatment.bioRxiv [Preprint]. 10.1101/2020.08.25.267336

  • 36

    International Agency for Research on Cancer (2021). GLOBOCAN: Cancer Incidence and Mortality Worldwide: International Agency for Research on Cancer [online]. Available online at: http://globocan.iarc.fr(accessed May 21, 2021).

  • 37

    International Cancer Genome ConsortiumHudsonT. J.AndersonW.ArtezA.BarkerA. D.BellC.et al (2010). International network of cancer genome projects.Nature464993998. 10.1038/nature08987

  • 38

    JonesS.ZhangX.ParsonsD. W.LinJ. C.LearyR. J.AngenendtP.et al (2008). Core signalling pathways in human pancreatic cancers revealed by global genomic analyses.Science32118011806. 10.1126/science.1164368

  • 39

    KalimuthuS. N.WilsonG. W.GrantR. C.SetoM.O’KaneG.VajpeyiR.et al (2020). Morphological classification of pancreatic ductal adenocarcinoma that predicts molecular subtypes and correlates with clinical outcome.Gut69317328. 10.1136/gutjnl-2019-318217

  • 40

    KandimallaR.TomiharaH.BanwaitJ. K.YamamuraK.SinghG.BabaH.et al (2020). A 15-gene immune, stromal, and proliferation gene signature that significantly associates with poor survival in patients with pancreatic ductal adenocarcinoma.Clin. Cancer Res.2636413648. 10.1158/1078-0432

  • 41

    KatzM. H. G.ShiQ.MeyersJ. P.HermanJ. M.ChoungM.WolpinB. M.et al (2021). Alliance A021501: preoperative mFOLFIRINOX or mFOLFIRINOX plus hypofractionated radiation therapy (RT) for borderline resectable (BR) adenocarcinoma of the pancreas.J. Clin. Oncol.39 (Suppl. 3):377. 10.1200/JCO.2021.39.3_suppl.377

  • 42

    KleeffJ.KorcM.ApteM.La VecchiaC.JohnsonC. D.BiankinA. V.et al (2016). Pancreatic cancer.Nat. Rev. Dis. Primers2:16022. 10.1038/nrdp.2016.22

  • 43

    KnudsenE. S.BalajiU.MannakeeB.VailP.EslingerC.MoxomC.et al (2018). Pancreatic cancer cell lines as patient-derived avatars: genetic characterisation and functional utility.Gut67508520. 10.1136/gutjnl-2016-313133

  • 44

    KubokiY.FischerC. G.Beleva GuthrieV.HuangW.YuJ.ChianchianoP.et al (2019). Single-cell sequencing defines genetic heterogeneity in pancreatic cancer precursor lesions.J. Pathol.247347356. 10.1002/path.5194

  • 45

    KunzmannV.SivekeJ. T.AlgülH.GoekkurtE.SieglerG.MartensU.et al (2021). Nab-paclitaxel plus gemcitabine versus nab-paclitaxel plus gemcitabine followed by FOLFIRINOX induction chemotherapy in locally advanced pancreatic cancer (NEOLAP-AIO-PAK-0113): a multicentre, randomised, phase 2 trial.Lancet Gastroenterol. Hepatol.6128138. 10.1016/S2468-1253(20)30330-7

  • 46

    LeD. T.DurhamJ. N.SmithK. N.WangH.BartlettB. R.AulakhL. K.et al (2017). Mismatch repair deficiency predicts response of solid tumors to PD-1 blockade.Science357409413. 10.1126/science.aan6733

  • 47

    LigorioM.SilS.Malagon-LopezJ.NiemanL. T.MisaleS.Di PilatoM.et al (2019). Stromal microenvironment shapes the intratumoral architecture of pancreatic cancer.Cell178160175.e27. 10.1016/j.cell.2019.05.012

  • 48

    LomberkG.BlumY.NicolleR.NairA.GaonkarK. S.MarisaL.et al (2018). Distinct epigenetic landscapes underlie the pathobiology of pancreatic cancer subtypes.Nat. Commun.9:1978. 10.1038/s41467-018-04383-6

  • 49

    MarabelleA.LeD. T.AsciertoP. A.Di GiacomoA. M.De Jesus-AcostaA.DelordJ. P.et al (2020). Efficacy of pembrolizumab in patients with noncolorectal high microsatellite instability/mismatch repair-deficient cancer: results from the phase II KEYNOTE-158 study.J. Clin. Oncol.38110. 10.1200/JCO.19.02105

  • 50

    MarcusL.LemeryS. J.KeeganP.PazdurR. (2019). FDA approval summary: pembrolizumab for the treatment of microsatellite instability-High solid tumors.Clin. Cancer Res.2537533758. 10.1158/1078-0432

  • 51

    MartinelliP.Carrillo-de Santa PauE.CoxT.SainzB.Jr.DusettiN.GreenhalfW.et al (2017). GATA6 regulates EMT and tumour dissemination, and is a marker of response to adjuvant chemotherapy in pancreatic cancer.Gut6616651676. 10.1136/gutjnl-2015-311256

  • 52

    Martinez-BalibreaE.Martínez-CardúsA.GinésA.Ruiz de PorrasV.MoutinhoC.LayosL.et al (2015). Tumor-related molecular mechanisms of oxaliplatin resistance.Mol. Cancer Ther.1417671776. 10.1158/1535-7163

  • 53

    MateoJ.ChakravartyD.DienstmannR.JezdicS.Gonzalez-PerezA.Lopez-BigasN.et al (2018). A framework to rank genomic alterations as targets for cancer precision medicine: the ESMO Scale for Clinical Actionability of Molecular Targets (ESCAT).Ann. Oncol.2918951902. 10.1093/annonc/mdy263

  • 54

    MaurerC.HolmstromS. R.HeJ.LaiseP.SuT.AhmedA.et al (2019). Experimental microdissection enables functional harmonisation of pancreatic cancer subtypes.Gut6810341043. 10.1136/gutjnl-2018-317706

  • 55

    MoffittR. A.MarayatiR.FlateE. L.VolmarK. E.LoezaS. G.HoadleyK. A.et al (2015). Virtual microdissection identifies distinct tumor- and stroma-specific subtypes of pancreatic ductal adenocarcinoma.Nat. Genet.4711681178. 10.1038/ng.3398

  • 56

    MoncadaR.BarkleyD.WagnerF.ChiodinM.DevlinJ. C.BaronM.et al (2020). Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas.Nat. Biotechnol.38333342.

  • 57

    MoralJ. A.LeungJ.RojasL. A.RuanJ.ZhaoJ.SethnaZ.et al (2020). ILC2s amplify PD-1 blockade by activating tissue-specific cancer immunity.Nature579130135. 10.1038/s41586-020-2015-4

  • 58

    MoseleF.RemonJ.MateoJ.WestphalenC. B.BarlesiF.LolkemaM. P.et al (2020). Recommendations for the use of next-generation sequencing (NGS) for patients with metastatic cancers: a report from the ESMO Precision Medicine Working Group.Ann. Oncol.3114911505. 10.1016/j.annonc.2020.07.014

  • 59

    NeoptolemosJ. P.KleeffJ.MichlP.CostelloE.GreenhalfW.PalmerD. H. (2018). Therapeutic developments in pancreatic cancer: current and future perspectives.Nat. Rev. Gastroenterol. Hepatol.15333348. 10.1038/s41575-018-0005-x

  • 60

    NeoptolemosJ. P.PalmerD. H.GhanehP.PsarelliE. E.ValleJ. W.HalloranC. M.et al (2017). Comparison of adjuvant gemcitabine and capecitabine with gemcitabine monotherapy in patients with resected pancreatic cancer (ESPAC-4): a multicentre, open-label, randomised, phase 3 trial.Lancet38910111024. 10.1016/S0140-6736(16)32409-6

  • 61

    NeoptolemosJ. P.StockenD. D.BassiC.GhanehP.CunninghamD.GoldsteinD.et al (2010). Adjuvant chemotherapy with fluorouracil plus folinic acid vs gemcitabine following pancreatic cancer resection: a randomized controlled trial.JAMA30410731081. 10.1001/jama.2010.1275

  • 62

    NeoptolemosJ. P.StockenD. D.FriessH.BassiC.DunnJ. A.HickeyH.et al (2004). A randomized trial of chemoradiotherapy and chemotherapy after resection of pancreatic cancer.N. Engl. J. Med.35012001210. 10.1056/NEJMoa032295

  • 63

    NollE. M.EisenC.StenzingerA.EspinetE.MuckenhuberA.KleinC.et al (2016). CYP3A5 mediates basal and acquired therapy resistance in different subtypes of pancreatic ductal adenocarcinoma.Nat. Med.22278287. 10.1038/nm.4038

  • 64

    O’KaneG. M.GrünwaldB. T.JangG. H.MasoomianM.PicardoS.GrantR. C.et al (2020). GATA6 expression distinguishes classical and basal-like subtypes in advanced pancreatic cancer.Clin. Cancer Res.2649014910. 10.1158/1078-0432.CCR-19-3724

  • 65

    OettleH.PostS.NeuhausP.GellertK.LangrehrJ.RidwelskiK.et al (2007). Adjuvant chemotherapy with gemcitabine vs observation in patients undergoing curative-intent resection of pancreatic cancer: a randomized controlled trial.JAMA297267277. 10.1001/jama.297.3.267

  • 66

    PishvaianM. J.BlaisE. M.BrodyJ. R.LyonsE.DeArbeloaP.HendifarA.et al (2020). Overall survival in patients with pancreatic cancer receiving matched therapies following molecular profiling: a retrospective analysis of the know your tumor registry trial.Lancet Oncol.21508518. 10.1016/S1470-2045(20)30074-7

  • 67

    PorterR. L.MagnusN. K. C.ThaparV.MorrisR.SzabolcsA.NeyazA.et al (2019). Epithelial to mesenchymal plasticity and differential response to therapies in pancreatic ductal adenocarcinoma.Proc. Natl. Acad. Sci. U.S.A.1162683526845. 10.1073/pnas.1914915116

  • 68

    PuleoF.NicolleR.BlumY.CrosJ.MarisaL.DemetterP.et al (2018). Stratification of pancreatic ductal adenocarcinomas based on tumor and microenvironment features.Gastroenterology15519992013.e3. 10.1053/j.gastro.2018.08.033

  • 69

    RaphaelB. J.Cancer Genome Atlas Research Network (2017). Integrated genomic characterization of pancreatic ductal adenocarcinoma.Cancer Cell32185203.e13. 10.1016/j.ccell.2017.07.007Electronic address: ; Cancer Genome Atlas Research Network

  • 70

    RashidN. U.PengX. L.JinC.MoffittR. A.VolmarK. E.BeltB. A.et al (2020). Purity independent subtyping of tumors (PurIST), a clinically robust, single-sample classifier for tumor subtyping in pancreatic cancer.Clin. Cancer Res.268292. 10.1158/1078-0432.CCR-19-1467

  • 71

    SahaiE.AstsaturovI.CukiermanE.DeNardoD. G.EgebladM.EvansR. M.et al (2020). A framework for advancing our understanding of cancer-associated fibroblasts.Nat. Rev. Cancer20174186. 10.1038/s41568-019-0238-1

  • 72

    SchlitterA. M.SeglerA.SteigerK.MichalskiC. W.JägerC.KonukiewitzB.et al (2017). Molecular, morphological and survival analysis of 177 resected pancreatic ductal adenocarcinomas (PDACs): identification of prognostic subtypes.Sci. Rep.7:41064. 10.1038/srep41064

  • 73

    SiegelR. L.MillerK. D.FuchsH. E.JemalA. (2021). Cancer statistics, 2021.CA Cancer J. Clin.71733. 10.3322/caac.21654

  • 74

    SohalD. P. S.DuongM.AhmadS. A.GandhiN. S.BegM. S.Wang-GillamA.et al (2021). Efficacy of perioperative chemotherapy for resectable pancreatic adenocarcinoma: a phase 2 randomized clinical trial.JAMA Oncol.7421427. 10.1001/jamaoncol.2020.7328

  • 75

    SomervilleT. D.BiffiG.Daßler-PlenkerJ.HurS. K.HeX. Y.VanceK. E.et al (2020). Squamous trans-differentiation of pancreatic cancer cells promotes stromal inflammation.Elife9:e53381. 10.7554/eLife.53381

  • 76

    SpringfeldC.JägerD.BüchlerM. W.StrobelO.HackertT.PalmerD. H.et al (2019). Chemotherapy for pancreatic cancer.Presse Med.48 (3 Pt 2) e159e174. 10.1016/j.lpm.2019.02.025

  • 77

    SteeleC. W.KarimS. A.LeachJ. D. G.BaileyP.Upstill-GoddardR.RishiL.et al (2016). CXCR2 inhibition profoundly suppresses metastases and augments immunotherapy in pancreatic uctal adenocarcinoma.Cancer Cell29832845. 10.1016/j.ccell.2016.04.014

  • 78

    SteeleN. G.CarpenterE. S.KempS. B.SirihorachaiV. R.TheS.DelrosarioL.et al (2020). Multimodal mapping of the tumor and peripheral blood immune landscape in human pancreatic cancer.Nat. Cancer110971112. 10.1038/s43018-020-00121-4

  • 79

    StrobelO.NeoptolemosJ.JägerD.BüchlerM. W. (2019). Optimizing the outcomes of pancreatic cancer surgery.Nat. Rev. Clin. Oncol.161126. 10.1038/s41571-018-0112-1

  • 80

    SultanaA.Tudur SmithC.CunninghamD.StarlingN.TaitD.NeoptolemosJ. P.et al (2007). Systematic review, including meta-analyses, on the management of locally advanced pancreatic cancer using radiation/combined modality therapy.Br. J. Cancer9611831190. 10.1038/sj.bjc.6603719

  • 81

    TiriacH.BelleauP.EngleD. D.PlenkerD.DeschênesA.SomervilleT. D. D.et al (2018). Organoid profiling identifies common responders to chemotherapy in pancreatic cancer.Cancer Discov.811121129. 10.1158/2159-8290

  • 82

    UesakaK.BokuN.FukutomiA.OkamuraY.KonishiM.MatsumotoI.et al (2016). Adjuvant chemotherapy of S-1 versus gemcitabine for resected pancreatic cancer: a phase 3, open-label, randomised, non-inferiority trial (JASPAC 01).Lancet388248257. 10.1016/S0140-6736(16)30583-9

  • 83

    Von HoffD. D.ErvinT.ArenaF. P.ChioreanE. G.InfanteJ.MooreM.et al (2013). Increased survival in pancreatic cancer with nab-paclitaxel plus gemcitabine.N. Engl. J. Med.36916911703. 10.1056/NEJMoa1304369

  • 84

    WaddellN.PajicM.PatchA. M.ChangD. K.KassahnK. S.BaileyP.et al (2015). Whole genomes redefine the mutational landscape of pancreatic cancer.Nature518495501. 10.1038/nature14169

  • 85

    Wang-GillamA.LiC. P.BodokyG.DeanA.ShanY. S.JamesonG.et al (2016). Nanoliposomal irinotecan with fluorouracil and folinic acid in metastatic pancreatic cancer after previous gemcitabine-based therapy (NAPOLI-1): a global, randomised, open-label, phase 3 trial.Lancet387545557. 10.1016/S0140-6736(15)00986-1

  • 86

    WartenbergM.CibinS.ZlobecI.VassellaE.Eppenberger-CastoriS.TerraccianoL.et al (2018). Integrated genomic and immunophenotypic classification of pancreatic cancer reveals three distinct subtypes with prognostic/predictive significance.Clin. Cancer Res.2444444454. 10.1158/1078-0432.CCR-17-3401

  • 87

    WilliamsonC. T.MillerR.PembertonH. N.JonesS. E.CampbellJ.KondeA.et al (2016). ATR inhibitors as a synthetic lethal therapy for tumours deficient in ARID1A.Nat. Commun.7:13837. 10.1038/ncomms13837

  • 88

    Zentrum für Krebsregisterdaten (2021). Krebsregister in Deutschland e. V. (GEKID) und des Zentrums für Krebsregisterdaten (ZfKD) im Robert Koch-Institut herausgegeben. Available online at: https://www.krebsdaten.de/(accessed May 21, 2021).

Summary

Keywords

molecular subtypes, transcriptomes, structural variants, precision medicine, next generation sequencing, clinical trials, ESPAC

Citation

Zhou X, Hu K, Bailey P, Springfeld C, Roth S, Kurilov R, Brors B, Gress T, Buchholz M, An J, Wei K, Peccerella T, Büchler MW, Hackert T and Neoptolemos JP (2021) Clinical Impact of Molecular Subtyping of Pancreatic Cancer. Front. Cell Dev. Biol. 9:743908. doi: 10.3389/fcell.2021.743908

Received

19 July 2021

Accepted

17 September 2021

Published

05 November 2021

Volume

9 - 2021

Edited by

Zhi-Gang Zhang, Shanghai Jiao Tong University, China

Reviewed by

Mark Dings, Amsterdam University Medical Center, Netherlands; Andrew Liss, Massachusetts General Hospital and Harvard Medical School, United States

Updates

Copyright

*Correspondence: John P. Neoptolemos,

This article was submitted to Molecular and Cellular Pathology, a section of the journal Frontiers in Cell and Developmental Biology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics