PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptNIH Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Nat Med. Author manuscript; available in PMC Aug 28, 2013.
Published in final edited form as:
Published online Apr 3, 2011. doi:  10.1038/nm.2344
PMCID: PMC3755490
NIHMSID: NIHMS278279
Subtypes of Pancreatic Ductal Adenocarcinoma and Their Differing Responses to Therapy
Eric A. Collisson,1,2* Anguraj Sadanandam,1,7* Peter Olson,3,8 William J. Gibb,1,9 Morgan Truitt,3 Shenda Gu,1 Janine Cooc,6 Jennifer Weinkle,1 Grace E. Kim,4 Lakshmi Jakkula,1 Heidi S. Feiler,1 Andrew H. Ko,2 Adam B. Olshen,5 Kathleen L. Danenberg,6 Margaret A. Tempero,2 Paul T. Spellman,1 Douglas Hanahan,3,7 and Joe W. Gray1,10
1Life Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
2Division of Hematology and Oncology, Diabetes Center, University of California, San Francisco, CA 94143, USA
3Department of Biochemistry and Biophysics, University of California, San Francisco, CA 94143, USA
4Department of Pathology, University of California, San Francisco, CA 94143, USA
5Department of Epidemiology and Biostatistics and Helen Diller Family Comprehensive Cancer Center, University of California, San Francisco, CA 94143, USA
6Response Genetics Inc., Los Angeles, CA 90033, USA
7Swiss Institute for Experimental Cancer Research (ISREC), Swiss Federal Institute of Technology Lausanne (EPFL), Lausanne CH-1015, Switzerland
10Biomedical Engineering, Oregon Health and Science University
Corresponding Author: Biomedical Engineering, 3303 SW Bond Ave., CH13B, Portland, OR 97239 Tel: 503-494-6500, Fax: 503-418-9311 ; grayjo/at/ohsu.edu
*Equal Contributors
8Current Address: Pfizer, 10724 Science Center Drive, La Jolla, CA 92121
9Current Address: Genomic Health, 301 Penobscot Drive, Redwood City, CA 94063
Pancreatic ductal adenocarcinoma (PDA) is a lethal disease. Overall survival is typically six months from diagnosis1. Numerous phase III trials of agents effective in other malignancies have failed to benefit unselected PDA populations, although patients do occasionally respond. Studies in other solid tumors have shown that heterogeneity in response is determined, in part, by molecular differences between tumors. Further, treatment outcomes are improved by targeting drugs to tumor subtypes in which they are selectively effective, with breast2 and lung3 cancers providing recent examples. Identification of PDA molecular subtypes has been frustrated by a paucity of tumor specimens available for study. We have overcome this problem by combined analysis of transcriptional profiles of primary PDA samples from several studies along with human and mouse PDA cell lines. We define three PDA subtypes: classical, quasi-mesenchymal, and exocrine-like, and present evidence for clinical outcome and therapeutic response differences between them. We further define gene signatures for these subtypes that may have utility in stratifying patients for treatment and present preclinical model systems that may be used to identify new subtype specific therapies.
Global gene expression analysis has proved useful for subtype identification in many human tumor types4. We approached PDA subtype identification by first identifying intrinsically variable (standard deviation > 0.8) genes in two gene expression microarray datasets from resected PDA. We generated one dataset using microdissected PDA material (UCSF tumors, n=27) from clinical samples for which information on survival duration was available and the second was previously published (Badea, et al.)5. To increase power, we merged these two datasets using the distance weighted discrimination (DWD) method6,7 and included intrinsically variable genes common to both studies. We then performed nonnegative matrix factorization (NMF) analysis with consensus clustering8 to identify subtypes of the disease. This analysis supported up to three subtypes (cophenetic coefficient >0.99; Supplementary Figs. 1, 2a and Supplementary Tables 13). We then developed a gene signature by using subtypes defined in NMF analysis of the merged clinical datasets to supervise significance analysis of microarrays (SAM) analysis9 with false discovery rate (FDR) less than 0.001. This resulted in a 62 gene signature, designated PDAssigner. The three PDA subtypes in the merged clinical dataset and their expression of PDAssigner genes are shown in Fig. 1a. We designated these subtypes as classical, quasi-mesenchymal (QM-PDA) and exocrine-like, based on our interpretation of subtype specific gene expression. The classical subtype had high expression of adhesion-associated and epithelial genes, the QM-PDA subtype showed high expression of mesenchyme associated genes. The exocrine-like subtype showed relatively high expression of tumor cell derived digestive enzyme genes, with immunohistochemical staining supporting this observation (Supplementary Fig. 3). Analysis of PDAssigner gene expression in three additional published PDA expression datasets of unique origin, platform or processing1012 also supported these three subtypes (Supplementary Fig. 4) demonstrating the robust nature of the subtype classification in early stage PDA.
Figure 1
Figure 1
Subtypes of PDA in tumors and cell lines and their prognostic significance
Survival after PDA resection has been associated with many factors including stage (tumor size and nodal involvement) and grade (degree of differentiation)13, but no one factor has been consistently prognostic14,15. We found that stratification by PDA transcriptional subtype provided useful prognostic information in one PDA dataset (UCSF) for which clinical annotation was available. Specifically, patients with classical subtype tumors fared better than patients with QM-PDA subtype tumors after resection (p=0.038, log rank, Fig. 1b). In this same data set, stage and grade were not significantly related (p>0.99), stage was not significantly associated with subtype (p=0.40), while grade was (p=0.041) (univariate analyses). In a multivariate Cox regression model including stage and subtype, subtype was an independent predictor of overall survival (p=0.024) indicating that PDA subtype independently contributed prognostic information to pathological staging in PDA. Associations of PDA subtype with other clinical variables are summarized in Supplementary Table 4. This analysis supports the use of subtypes (as defined using PDAssigner) as an independent prognostic indicator in resected PDA.
Further validation of PDA subtypes and preclinical identification of subtype specific therapeutic agents would be facilitated by the availability of laboratory models of the subtypes. Therefore, we asked if the PDA subtypes were represented in a collection of 19 human and 15 mouse PDA cell lines. The 19 human PDA cell lines were selected from publicly available PDA lines while the 15 mouse lines were derived by us from genetically engineered Tp53−/− and INK4A−/−16 models of PDA. The analyses of the human and mouse PDA cell lines using the 62 PDAssigner genes are shown in Fig. 1c,d, as well as in Supplementary Figs. 2b–e. Supplementary Tables 5 and 6. These cell line collections contain representatives of the classical and QM-PDA subtypes, but not the exocrine-like subtype. The absence of the exocrine-like subtype in the cell line collection raises the possibility that this subtype is an artifact of contaminating normal pancreas tissue adjacent to tumor. However, the UCSF samples were prepared by microdissection to enrich for PDA cancer cells thereby minimizing contaminating tumor-associated stroma and/or adjacent normal exocrine pancreas. This dataset includes the exocrine-like subtype, which argues that it is a bona fide PDA subtype. Thus, we conclude that the cell line collections model two of the PDA subtypes thereby enabling exploration of biological differences between these two PDA subtypes and may facilitate screening for classical and QM-PDA subtype-specific therapeutic agents and biological properties.
Two genes associated with PDA subtypes, GATA binding protein 6 (GATA6) and v-ki-ras2 kirsten rat sarcoma viral oncogene homolog (KRAS), two variable genes in our UCSF PDA dataset, Supplementary Table 1a), have been implicated in both aspects of normal development and cancer pathophysiology in published studies. GATA-family transcription factors are associated with tissue specific differentiation and have been demonstrated to be subtype specific markers in other cancers. For example, GATA binding protein 3 (GATA3) is required for luminal differentiation in the breast17, and high GATA3 characterizes luminal subtype breast cancers18. Likewise, GATA6 is essential for pancreatic development19 and has been implicated in PDA20,21. GATA6 is highly expressed in most classical subtype tumors and cell lines, and comparatively low in the QM-PDA cell lines and tumors. Additionally, a previously published gene signature associated with GATA6 overexpression20 is enriched in the classical subtype (Supplementary Fig. 5). Seeking to establish a possible functional role underlying the observed differences in GATA6 expression, we assessed the impact of GATA6-directed RNAi knockdown on colony formation in soft agar in the classical and QM-PDA cell lines. GATA6 depletion impaired anchorage-independent growth in classical PDA cell lines, but not in QM-PDA cell lines (Fig. 2), consistent with a functional, subtype-specific role for this transcription factor in the classical PDA subtype.
Figure 2
Figure 2
Classical PDA subtype and the GATA6 transcription factor
Recent work in the mouse demonstrates that PDA can arise from a variety of precursor cells by activating KRAS in distinct cellular compartments of the pancreas22. Others have found that cancer cell lines harboring mutant KRAS differ in their dependence on KRAS23. These findings imply plasticity in either reliance on KRAS signaling or a cell-type specific role for mutant KRAS in different cells of origin/lineages in PDA, or both. They further suggested to us that despite KRAS mutation in most PDAs, KRAS dependence might differ by PDA subtype. We found KRAS mRNA levels elevated in classical subtype PDAs relative to the other subtypes (Supplementary Fig. 6, p<0.05 in UCSF samples). We explored the relationship between KRAS dependence and subtype by using RNAi to probe KRAS-mutant human PDA cell lines for dependence on KRAS. Classical PDA lines proved to be relatively more dependent on KRAS than QM-PDA lines (Fig. 3). Further, a previously reported signature of KRAS-addiction23 is enriched in the classical subtype (Supplementary Fig. 7). These results suggest that KRAS-directed therapy, while not yet a clinical reality, might be best deployed in classical PDA. Mouse models capable of sequentially activating and then deleting mutant KRAS would further these observations to genetically define the respective roles mutant KRAS plays in both the initiation and ongoing maintenance of PDA.
Figure 3
Figure 3
Classical subtype cells depend on KRas
We tested the possibility that PDA subtypes with different biological characteristics might have subtype-specific drug responses by measuring responses to gemcitabine and erlotinib (the backbone of current treatment regimens24) in human PDA cell lines of known subtype. QM-PDA subtype lines were, on average, more sensitive to gemcitabine than the classical subtype (Fig. 4). Conversely, erlotinib was more effective in classical subtype cell lines. This suggests that KRAS mutation status is an imperfect predictor of sensitivity to EGFR-targeted therapy in PDA, an observation consistent with findings in nonsmall cell lung25 and colorectal cancers26, and implies that cancer cells dependent on mutant KRAS still employ the EGFR to some extent. These results further establish phenotypic differences between the classical and QM-PDA subtypes, and suggest that these and perhaps additional drugs will show subtype specificity in PDA, a distinction that could be exploited in clinical trial sensitivity enrichment schemes. More immediately, these results indicate that gemcitabine and erlotinib are preferentially active in different PDA subtypes, so that the current practice of combining them may increase toxicity without increasing efficacy for many patients. Combining agents with similar subtype specificity might be considered instead.
Figure 4
Figure 4
Drug Responses Differ by Subtype
In summary, our data support the existence of three intrinsic subtypes of PDA that progress at different rates clinically and may respond differently to selected therapies. The validity of these subtypes is supported by analysis of multiple primary clinical datasets as well as of PDA cell lines both from human tumors and from mouse models engineered with the hallmark genetic lesions of human PDA. Knowledge of these subtypes may motivate investigation of associations between clinico-pathologic variables and these subtypes by collaborative consortia examining the molecular diversity of PDA27. The markers that define these subtypes may have prognostic utility in risk-adapted surgical approaches28 or predictive utility in sensitivity enrichment schemes. The use of subtyped human and mouse PDA preclinical models promises to accelerate identification of subtype specific functional and regulatory processes that can be exploited to therapeutic benefit. In turn, such assay systems could be used to screen therapeutic approaches, empirically or based on mechanism, to identify those that are potent against PDA, either in a specific subtype that would then be used to personalize treatment29, or spanning the subtypes with possible therapeutic generality.
Clinical Samples
After institutional review board approval, we selected archival material from pancreatic ductal adenocarcinoma resections performed at the University of California, San Francisco between 1993 and 2006. G.E.K. reviewed all cases prior to inclusion in the study. Tissue processing is described in Supplementary Methods.
Merging of Microarray Datasets
After processing of microarrays (as described in Supplementary Information), we screened the UCSF and Badea et al.,5 PDA datasets separately by selecting probes with standard deviation (SD) > 0.8. We then merged SD-selected datasets using DWD7 as described6. We column (samples) normalized to N(0,1) and row (probe or gene) normalized each dataset by median centering. We merged the processed datasets using DWD and finally, we median centered the rows.
NMF and SAM Analysis
We analyzed the merged datasets by consensus clustering-based NMF8 to identify stable subtypes using R code from GenePattern30. See supplement for details regarding the interpretation of subtypes derived from consensus-based NMF clustering. We identified PDAssigner genes using three-class SAM analysis based on classes identified from NMF analysis using the Bioconductor31 package, Siggenes, and generated heatmaps of samples by PDAssigner genes using Cluster software32. For cell line classification, we merged the cell line datasets with core PDA clinical datasets after selection of the 62 PDAssigner genes from each dataset followed by DWD based merging. We visualized datasets using the Hierarchical Clustering Viewer (HCV) from GenePattern30.
Clinical Outcome Analysis
We calculated overall survival in days from the time of PDA resection until date of death as defined by the State of California Death Registry and clinical records. We employed the log-rank test for univariate associations with survival or the Cox proportional hazards model for multivariate modeling of survival. We applied Fisher’s exact test to investigate the relationships among subtype, stage and grade. We used the R language for all analyses. We drew the survival curve using web-based GenePattern30.
Drug Sensitivity
We plated 2.5x103 cells per well on day 0, treated with erlotinib or gemcitabine in nine, five-fold serial dilutions on day 1 and estimated cell number using Cell Titre Glow assay (CTG, Promega) on day 4. IC50 was calculated using the Calcusyn program (Biosoft).
RNAi
We obtained validated pLKO-based shRNA vectors shKRAS#533 from Dr. B.R. Stockwell (NYU) and shGATA6#5, and shLuc34 from Dr. R Adam, (Boston Children’s Hospital). We packaged lentiviruses, transduced cells and then selected in puromycin for 48 hours. For shKRAS proliferation experiments, we plated 2.5 x103 transduced cells per well on day 0 in 96 well plates, then counted one plate on day one and the other plate on day four. We confirmed protein knockdown by western blotting using the Odyssey system, with 10ug per lane of total protein and the c19 KRas antibody (Santa Cruz), normalized to total actin (I-19, Santa Cruz) and compared to pLKOshLuc -KRas levels. For GATA6 knockdown experiments, after puromycin selection cells we trypsinized and plated transduced cells in soft agar as described35. We assessed GATA6 transcript levels on the day of plating in soft agar as described34.
See Supplementary Information for detailed methods.
Supplementary Material
Acknowledgments
We are grateful to M. Lenburg, the Gray, Hanahan, and Speed labs for discussion. We thank L. Chin (Dana-Farber Cancer Institute), S. Batra (University of Nebraska Medical Center), M. McMahon (University of California San Francisco) and A. Singh (Massachusetts General Hospital) for cell lines, and B. Stockwell (New York University) and R. Adam (Children’s Hospital Boston) for shRNA. E.A.C was supported by a Young Investigator Award from the American Society of Clinical Oncology and US National Cancer Institute (NCI) K08 CA137153. A.S. was supported by Department of Defense (DOD) Postdoctoral Fellowship (BC087768). The research in the laboratory of D.H. was supported by a NCI Program Project Grant PO1 CA 117969; D.H. is an American Cancer Society Research Professor. This work was supported by the Director, Office of Science, Office of Biological & Environmental Research, of the United States Department of Energy under Contract No. DE-AC02-05CH11231, by the NIH/NCI grants P50 CA 58207, P50 CA 83639 and by the U54 CA 112970 to J.W.G.
Footnotes
Author Contributions:
E.A.C. and A.S. designed, conducted and interpreted experiments, and wrote the manuscript. P.O., W.J.G, M.T., S.G., J.C., J.W., L.J., and H.S.F. performed experiments. K.L.D., and P.T.S. provided support and interpreted experiments. G.E.K. and A.H.K. coordinated clinical sample acquisition. A.B.O. provided statistical expertise. M.A.T provided support, interpreted experiments and coordinated clinical sample acquisition. D.H. and J.W.G. designed and interpreted experiments, wrote the manuscript and supervised the project.
1. Stathis A, Moore MJ. Advanced pancreatic carcinoma: current treatment and future challenges. Nat Rev Clin Oncol. 2010;7:163–172. [PubMed]
2. Slamon DJ, et al. Use of chemotherapy plus a monoclonal antibody against HER2 for metastatic breast cancer that overexpresses HER2. N Engl J Med. 2001;344:783–792. [PubMed]
3. Lynch TJ, et al. Activating mutations in the epidermal growth factor receptor underlying responsiveness of non-small-cell lung cancer to gefitinib. N Engl J Med. 2004;350:2129–2139. [PubMed]
4. Alizadeh AA, et al. Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling. Nature. 2000;403:503–511. [PubMed]
5. Badea L, Herlea V, Dima SO, Dumitrascu T, Popescu I. Combined gene expression analysis of whole-tissue and microdissected pancreatic ductal adenocarcinoma identifies genes specifically overexpressed in tumor epithelia. Hepatogastroenterology. 2008;55:2016–2027. [PubMed]
6. Herschkowitz JI, et al. Identification of conserved gene expression features between murine mammary carcinoma models and human breast tumors. Genome Biol. 2007;8:R76. [PMC free article] [PubMed]
7. Benito M, et al. Adjustment of systematic microarray data biases. Bioinformatics. 2004;20:105–114. [PubMed]
8. Brunet JP, Tamayo P, Golub TR, Mesirov JP. Metagenes and molecular pattern discovery using matrix factorization. Proc Natl Acad Sci U S A. 2004;101:4164–4169. [PubMed]
9. Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98:5116–5121. [PubMed]
10. Balagurunathan Y, et al. Gene expression profiling-based identification of cell-surface targets for developing multimeric ligands in pancreatic cancer. Mol Cancer Ther. 2008;7:3071–3080. [PMC free article] [PubMed]
11. Pei H, et al. FKBP51 affects cancer cell response to chemotherapy by negatively regulating Akt. Cancer Cell. 2009;16:259–266. [PMC free article] [PubMed]
12. Grutzmann R, et al. Gene expression profiling of microdissected pancreatic ductal carcinomas using high-density DNA microarrays. Neoplasia. 2004;6:611–622. [PMC free article] [PubMed]
13. Brennan MF, Kattan MW, Klimstra D, Conlon K. Prognostic nomogram for patients undergoing resection for adenocarcinoma of the pancreas. Ann Surg. 2004;240:293–298. [PubMed]
14. Biankin AV, et al. Expression of S100A2 Calcium-Binding Protein Predicts Response to Pancreatectomy for Pancreatic Cancer. Gastroenterology. 2009 [PubMed]
15. Blackford A, et al. SMAD4 gene mutations are associated with poor prognosis in pancreatic cancer. Clin Cancer Res. 2009;15:4674–4679. [PMC free article] [PubMed]
16. Bardeesy N, et al. Both p16(Ink4a) and the p19(Arf)-p53 pathway constrain progression of pancreatic adenocarcinoma in the mouse. Proc Natl Acad Sci U S A. 2006;103:5947–5952. [PubMed]
17. Kouros-Mehr H, Slorach EM, Sternlicht MD, Werb Z. GATA-3 maintains the differentiation of the luminal cell fate in the mammary gland. Cell. 2006;127:1041–1055. [PMC free article] [PubMed]
18. Mehra R, et al. Identification of GATA3 as a breast cancer prognostic marker by global gene expression meta-analysis. Cancer Res. 2005;65:11259–11264. [PubMed]
19. Decker K, Goldman DC, Grasch CL, Sussel L. Gata6 is an important regulator of mouse pancreas development. Dev Biol. 2006;298:415–429. [PMC free article] [PubMed]
20. Kwei KA, et al. Genomic profiling identifies GATA6 as a candidate oncogene amplified in pancreatobiliary cancer. PLoS Genet. 2008;4:e1000081. [PMC free article] [PubMed]
21. Fu B, Luo M, Lakkur S, Lucito R, Iacobuzio-Donahue CA. Frequent genomic copy number gain and overexpression of GATA-6 in pancreatic carcinoma. Cancer Biol Ther. 2008;7:1593–1601. [PubMed]
22. Gidekel Friedlander SY, et al. Context-dependent transformation of adult pancreatic cells by oncogenic K-Ras. Cancer Cell. 2009;16:379–389. [PMC free article] [PubMed]
23. Singh A, et al. A gene expression signature associated with “K-Ras addiction” reveals regulators of EMT and tumor cell survival. Cancer Cell. 2009;15:489–500. [PMC free article] [PubMed]
24. Moore MJ, et al. Erlotinib plus gemcitabine compared with gemcitabine alone in patients with advanced pancreatic cancer: a phase III trial of the National Cancer Institute of Canada Clinical Trials Group. J Clin Oncol. 2007;25:1960–1966. [PubMed]
25. Yauch RL, et al. Epithelial versus mesenchymal phenotype determines in vitro sensitivity and predicts clinical activity of erlotinib in lung cancer patients. Clin Cancer Res. 2005;11:8686–8698. [PubMed]
26. De Roock W, et al. Association of KRAS p.G13D mutation with outcome in patients with chemotherapy-refractory metastatic colorectal cancer treated with cetuximab. Jama. 2010;304:1812–1820. [PubMed]
27. International network of cancer genome projects. Nature. 2007;464:993–998. [PMC free article] [PubMed]
28. Raut CP, Evans DB, Crane CH, Pisters PW, Wolff RA. Neoadjuvant therapy for resectable pancreatic cancer. Surg Oncol Clin N Am. 2004;13:639–661. ix. [PubMed]
29. McDermott U, et al. Identification of genotype-correlated sensitivity to selective kinase inhibitors by using high-throughput tumor cell line profiling. Proc Natl Acad Sci U S A. 2007;104:19936–19941. [PubMed]
30. Kuehn H, Liberzon A, Reich M, Mesirov JP. Using GenePattern for gene expression analysis. Curr Protoc Bioinformatics. 2008;Chapter 7(Unit 7):12. [PubMed]
31. Gentleman RC, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80. [PMC free article] [PubMed]
32. Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998;95:14863–14868. [PubMed]
33. Yagoda N, et al. RAS-RAF-MEK-dependent oxidative cell death involving voltage-dependent anion channels. Nature. 2007;447:864–868. [PMC free article] [PubMed]
34. Kanematsu A, Ramachandran A, Adam RM. GATA-6 mediates human bladder smooth muscle differentiation: involvement of a novel enhancer element in regulating alpha-smooth muscle actin gene expression. Am J Physiol Cell Physiol. 2007;293:C1093–1102. [PubMed]
35. Collisson EA, De A, Suzuki H, Gambhir SS, Kolodney MS. Treatment of metastatic melanoma with an orally available inhibitor of the Ras-Raf-MAPK cascade. Cancer Res. 2003;63:5669–5673. [PubMed]