PMCCPMCCPMCC

Search tips
Search criteria 

Advanced


 
Formats
Article sections
Authors
Related links
Logo of plosmedPLoS MedicineSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
 
PLoS Med. 2010 July; 7(7): e1000307.
Published online 2010 July 13. doi:  10.1371/journal.pmed.1000307
PMCID: PMC2903589
Copyright Stratford et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
A Six-Gene Signature Predicts Survival of Patients with Localized Pancreatic Ductal Adenocarcinoma
Jeran K. Stratford,1 David J. Bentrem,2 Judy M. Anderson,3 Cheng Fan,4 Keith A. Volmar,5 J. S. Marron,4,6 Elizabeth D. Routh,4 Laura S. Caskey,4 Jonathan C. Samuel,7 Channing J. Der,1,4 Leigh B. Thorne,5,6 Benjamin F. Calvo,4,7 Hong Jin Kim,4,7 Mark S. Talamonti,8 Christine A. Iacobuzio-Donahue,9 Michael A. Hollingsworth,3 Charles M. Perou,4,10 and Jen Jen Yeh1,4,7*
1Department of Pharmacology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
2Department of Surgery and Robert H. Lurie Comprehensive Cancer Center, Northwestern University, Feinberg School of Medicine, Chicago, Illinois, United States of America
3The Eppley Cancer Institute, University of Nebraska, Omaha, Nebraska, United States of America
4Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
5Department of Pathology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
6Department of Statistics and Operations Research, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
7Department of Surgery, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
8Department of Surgery, NorthShore University HealthSystem, Baltimore, Maryland, United States of America
9Department of Pathology, The Johns Hopkins Medical Institutions, Baltimore, Maryland, United States of America
10Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
Samir M. Hanash, Academic Editor
Fred Hutchinson Cancer Research Center, United States of America
* E-mail: jjyeh/at/med.unc.edu
ICMJE criteria for authorship read and met: JKS DJB JMA CF KAV JSM EDR LSC JCS CJD LBT BFC HJK MST CAID MAH CMP JJY. Agree with the manuscript's results and conclusions: JKS DJB JMA CF KAV JSM EDR LSC JCS CJD LBT BFC HJK MST CAID MAH CMP JJY. Designed the experiments/the study: JKS MAH JJY. Analyzed the data: JKS CF JSM CMP JJY. Collected data/did experiments for the study: JKS DJB JMA KAV EDR LSC JCS HJK JJY. Enrolled patients: DJB BFC HJK MST CAID JJY. Wrote the first draft of the paper: JJY. Contributed to the writing of the paper: JKS DJB CF JSM CJD HJK MST MAH CMP JJY. Created figures for the paper: JSM. Responsible for tissue procurement: LBT. Prepared tissue microarray: LBT.
Received January 4, 2010; Accepted June 3, 2010.
Small right arrow pointing to: This article has been cited by other articles in PMC.
Abstract
Background
Pancreatic ductal adenocarcinoma (PDAC) remains a lethal disease. For patients with localized PDAC, surgery is the best option, but with a median survival of less than 2 years and a difficult and prolonged postoperative course for most, there is an urgent need to better identify patients who have the most aggressive disease.
Methods and Findings
We analyzed the gene expression profiles of primary tumors from patients with localized compared to metastatic disease and identified a six-gene signature associated with metastatic disease. We evaluated the prognostic potential of this signature in a training set of 34 patients with localized and resected PDAC and selected a cut-point associated with outcome using X-tile. We then applied this cut-point to an independent test set of 67 patients with localized and resected PDAC and found that our signature was independently predictive of survival and superior to established clinical prognostic factors such as grade, tumor size, and nodal status, with a hazard ratio of 4.1 (95% confidence interval [CI] 1.7–10.0). Patients defined to be high-risk patients by the six-gene signature had a 1-year survival rate of 55% compared to 91% in the low-risk group.
Conclusions
Our six-gene signature may be used to better stage PDAC patients and assist in the difficult treatment decisions of surgery and to select patients whose tumor biology may benefit most from neoadjuvant therapy. The use of this six-gene signature should be investigated in prospective patient cohorts, and if confirmed, in future PDAC clinical trials, its potential as a biomarker should be investigated. Genes in this signature, or the pathways that they fall into, may represent new therapeutic targets.
Please see later in the article for the Editors' Summary
Editors' Summary
Background
Pancreatic cancer kills nearly a quarter of a million people every year. It begins when a cell in the pancreas (an organ lying behind the stomach that produces digestive enzymes and hormones such as insulin, which controls blood sugar levels) acquires genetic changes that allow it to grow uncontrollably and to spread around the body (metastasize). Nearly all pancreatic cancers are “pancreatic ductal adenocarcinomas” (PDACs)—tumors that start in the cells that line the tubes in the pancreas that take digestive juices to the gut. Because PDAC rarely causes any symptoms early in its development, it has already metastasized in about half of patients before it is diagnosed. Consequently, the average survival time after a diagnosis of PDAC is only 5–8 months. At present, the only chance for cure is surgical removal (resection) of the tumor, part of the pancreas, and other nearby digestive organs. The operation that is needed for the majority of patients—the Whipple procedure—is only possible in the fifth of patients whose tumor is found when it is small enough to be resectable but even with postoperative chemotherapy, these patients only live for 23 months after surgery on average, possibly because they have micrometastases at the time of their operation.
Why Was This Study Done?
Despite this poor overall outcome, about a quarter of patients with resectable PDAC survive for more than 5 years after surgery. Might some patients, therefore, have a less aggressive form of PDAC determined by the biology of the primary (original) tumor? If this is the case, it would be useful to be able to stratify patients according to the aggressiveness of their disease so that patients with very aggressive disease could be given chemotherapy before surgery (neoadjuvant therapy) to kill any micrometastases. At present neoadjuvant therapy is given to patients with locally advanced, unresectable tumors. In this study, the researchers compare gene expression patterns in primary tumor samples collected from patients with localized PDAC and from patients with metastatic PDAC between 1999 and 2007 to try to identify molecular markers that distinguish between more and less aggressive PDACs.
What Did the Researchers Do and Find?
The researchers identified a six-gene signature that was associated with metastatic disease using a molecular biology approach called microarray hybridization and a statistical method called significance analysis of microarrays to analyze gene expression patterns in primary tumor samples from 15 patients with localized PDAC and 15 patients with metastatic disease. Next, they used a training set of tumor samples from another 34 patients with localized and resected PDAC, microarray hybridization, and a graphical method called X-tile to select a combination of expression levels of the six genes that discriminated optimally between high-risk (aggressive) and low-risk (less aggressive) tumors on the basis of patient survival (a “cut-point”). When the researchers applied this cut-point to an independent set of 67 tumor samples from patients with localized and resected PDAC, they found that 42 patients had high-risk tumors. These patients had an average survival time of 15 months; 55% of them were alive a year after surgery. The remaining 25 patients, who had low-risk tumors, had an average survival time of 49 months and 91% of them were alive a year after resection.
What Do These Findings Mean?
These and other findings identify a six-gene signature that can predict outcomes in patients with localized, resectable PDAC better than, and independently of, established clinical markers of outcome. If the predictive ability of this signature can be confirmed in additional patients, it could be used to help patients make decisions about their treatment. For example, a patient wondering whether to risk the Whipple procedure (2%–6% of patients die during this operation and more than 50% have serious postoperative complications), the knowledge that their tumor was low risk might help them decide to have the operation. Conversely, a patient in poor health with a high-risk tumor might decide to spare themselves the trauma of major surgery. The six-gene signature might also help clinicians decide which patients would benefit most from neoadjuvant therapy. Finally, the genes in this signature, or the biological pathways in which they participate, might represent new therapeutic targets for the treatment of PDAC.
Additional Information
Please access these Web sites via the online version of this summary at http://dx.doi.org/10.1371/journal.pmed.1000307.
Pancreatic ductal adenocarcinoma (PDAC), comprising over 90% of all pancreatic cancers, remains a lethal disease with an estimated 232,000 new cases, 227,000 deaths per year worldwide, and a less than 5% 5-y survival rate [1],[2]. Currently the standard of care for the 20% of patients with localized disease is surgery followed by chemotherapy, and in some cases radiation. Unfortunately, despite the use of adjuvant therapy, median survival remains at best 23 mo [3]. It is important to note, however, that up to 27% of patients with resected PDAC can survive for 5 y [4]–[10]. However, in studies examining actual long-term survivors [4]–[10], only two have found that adjuvant therapy was associated with improved survival [9],[10]. In addition, randomized controlled trials of gemcitabine-based chemotherapy demonstrate an improvement in median survival of at best 3 mo [3],[11]. One possible conclusion from these studies is that tumor biology dictates outcome and that our current adjuvant therapy has only a modest impact on altering a patient's course.
Hypothesizing that the dismal outcome of patients with localized disease is due to the presence of micrometastatic disease, current clinical investigation has focused on preoperative or neoadjuvant therapy [12],[13]. This approach, in which patients who cannot tolerate the stress of therapy or who develop metastatic disease during treatment are spared surgery, has demonstrated an overall survival of 34 mo in this highly selected patient population [12],[13]. Therefore the ability to select patients who would most benefit from a neoadjuvant approach may be important. One way to select these individuals is to define a prognostic gene signature that can identify patients with more aggressive tumor biology upfront.
Expression profiling of PDAC has lead to further studies of additional molecular diagnostic and prognostic markers [14]–[19]. However, the search for genes of biological significance in these large datasets continues to be challenging. One approach to identify genes or pathways that are biologically relevant is to study those that are of prognostic significance [20]. Lowe and colleagues found differential gene expression changes associated with nodal status in primary PDAC [21], suggesting that molecular differences in primary PDAC do exist. We hypothesized that by comparing primary PDAC tumors at the extremes of disease, we would identify molecular changes reflective of differences in biology within primary PDAC tumors.
Patients
PDAC samples from 15 patients with resected primary PDAC from the University of North Carolina at Chapel Hill (UNC) and 15 patients with metastatic PDAC from the University of Nebraska Medical Center Rapid Autopsy Pancreatic Program (NEB) were used to derive differentially expressed genes associated with metastatic disease. For the NEB samples, human pancreatic tumors from decedents who had previously been diagnosed with PDAC, and who generously consented to post mortem examinations, were obtained from the institutional review board (IRB)-approved NEB Tissue Bank. To ensure minimal degradation of tissue, organs were harvested within 3 h post mortem and the specimens flash frozen in liquid nitrogen.
The training cohort included 34 patients with resected PDAC from Johns Hopkins Medical Institutions (JHMI). The testing or validation cohort included patients from two institutions: 48 from Northwestern Memorial Hospital (NW) and 19 from NorthShore University HealthSystem (NSU). All samples were collected between 1999 and 2007 at the time of operation and flash frozen in liquid nitrogen after approval by each individual IRB. The UNC IRB approved use of all de-identified samples for this study. All available samples were reviewed by a single pathologist (KAV). De-identified data including tumor, node, and metastasis (TNM), grade or differentiation, margin status, and survival were available for the majority of patients.
RNA Isolation and Microarray Hybridization
All RNA isolation and hybridization was performed on Agilent (Agilent Technologies) human whole genome 4×44 K DNA microarrays and at UNC. RNA was extracted from macrodissected snap-frozen tumor samples using Allprep Kits (Qiagen) and quantified using nanodrop spectrophotometry (ThermoScientific). RNA quality was assessed with the use of the Bioanalyzer 2100 (Agilent Technologies). RNA was selected for hybridization using RNA integrity number and by inspection of the 18S and 28S ribosomal RNA. Similar RNA quality was selected across samples. One microgram of RNA was used as a template for DNA preparations and hybridized to Agilent 4×44 K whole human genome arrays (Agilent Technologies). cDNA was labeled with Cy5-dUTP and a reference control (Stratagene) was labeled with Cy3-dUTP using the Agilent (Agilent Technologies) low RNA input linear amplification kit and hybridized overnight at 65°C to Agilent 4×44 K whole human genome arrays (Agilent Technologies). Arrays were washed and scanned using an Agilent scanner (Agilent Technologies). The data are publicly available in Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE21501).
Microarray and Statistical Analysis
All array data were normalized using Lowess normalization. Data were excluded for genes with poor spot quality or genes that did not have mean intensity greater than 10 for one of the two channels (green and red) in at least 70% of the experiments. The log2 ratio of the mean red intensity over mean green intensity was calculated for each gene and went through LOWESS normalization [22]. Missing data were imputed using the k-nearest neighbors imputation (KNN) with k = 10 [23]. A distance weighted discrimination (DWD) was used to detect the systematic biases between the different datasets and then global adjustments made to remove these biases [24]. Genes that were significantly up- or down-regulated were identified using significance analysis of microarrays (SAM) [25]. Two centroids were created using the mean gene expression profile of this significant gene list from the derivation set and used to develop a single sample predictor (SSP, nearest centroid algorithm) [26] for an objective classifier. After DWD, the SSP was applied to a 34-patient training set where any new sample was compared to the resected centroid and assigned by the SSP distance function to the resected centroid using (1 − Pearson correlation coefficient). The X-Tile software program, which assigns a two-population log-rank value to each sample and then determines the best cut-point, was used to determine the best threshold for classifying samples into high- and low-risk categories [27]. X-Tile predicted that the (1− Pearson correlation coefficient) distance of 1 would be the appropriate cut-point to stratify patients into a high- and low-risk group (p = 0.006). A second independent validation cohort was then used as a test set using this predetermined cut-point to evaluate outcome.
Survival analysis was performed using the statistical software programs R, the R-package “survival,” and SPSS (SPSS, Inc.). Overall survival (OS) was analyzed using the Kaplan-Meier product-limit method and the significance of our variables was measured by the log-rank test. The Fisher exact test was used to analyze associations between two variables, the Pearson Chi-square test was used to analyze association between more than two variables. Multivariable analysis and analysis of continuous and ordinal variables was performed using the Cox proportional hazards regression method.
Tissue Microarrays
Tissue microarrays (TMAs; UNC2) were prepared from formalin-fixed paraffin-embedded tissue sections using a 2-mm punch. The arrays contained triplicate cores of matched normal and tumor tissue as well as chronic pancreatitis when available, from each patient. We prepared 5-µm sections from each TMA block. Hematoxylin and eosin stained slides from each TMA block were reviewed by a pathologist (KAV) to ensure that tissues were cored accurately.
Immunohistochemistry
Slides with 5-µM sections from the paraffin-embedded specimens were deparaffinized and rehydrated. The slides were then subjected to alkaline heat antigen-retrieval using 1% Tris EDTA for 20 min in a steamer. All slides were incubated with 3% H2O2 for 5 min and washed with TBS. The slides were further treated with protein block solution (bovine serum albumin) for 20 min. The sections were incubated with primary KLF6 1[ratio]150 antibody (sc-7158, Santa Cruz Biotechnology) for 60 min at room temperature. Following a TBS wash, the slides were incubated with secondary labeled Polymer-HRP anti-rabbit (Dako K4002) for 30 min. This was followed by a 5-min incubation with the substrate-chromogen, 3,3′-diaminobenzidine (Vector SK-4100). The sections were counterstained with Harris hematoxylin. Positive KLF6 staining was defined as when more than 5% of cells expressed the marker and graded from 0 (no staining) to 4 (strong staining). The results of each protein marker were then expressed as intensity (I) and proportion (P) of positive epithelial cells and the score as the product of I and P [28],[29]. All stained slides were reviewed in a blinded fashion (JMA).
Patient and Tumor Characteristics
In order to study the extremes of PDAC tumor biology, we collected a diverse set of resected PDAC specimens from patients with and without metastases. As the tumor microenvironment is increasingly recognized to play a critical role in tumorigenesis [30]–[33], tissues were macrodissected in order to preserve the normal adjacent tissue and stroma of the tumors. The characteristics of the dataset used to derive the signature (derivation set) comprised 15 primary resected PDAC tumors (UNC1) and 15 primary tumors from patients with metastatic PDAC (NEB). The training set comprised 34 patients with primary PDAC and the independent validation test set comprised 67 patients with primary PDAC (Tables 1 and ​and2).2). There were no differences in RNA quality between the decedent and resected PDAC samples. Available treatment data of the patients in the training and test sets are also shown. One of 15 (7%) UNC1 patients received preoperative or neoadjuvant chemotherapy and 11/15 (73%) NEB patients received chemotherapy less than 6 mo prior to death. No patient in the 34-patient training set received neoadjuvant chemotherapy. Only 3% (2/67) of patients in the test set received neoadjuvant chemotherapy and 45% (30/67) of patients received postoperative or adjuvant chemotherapy.
Table 1
Table 1
Patient, tumor, and treatment characteristics in the derivation set.
Table 2
Table 2
Patient, tumor, and treatment characteristics in the training and testing sets.
Gene Expression Differences in Nonmetastatic And Metastatic Primary Tumors
We hypothesized that we could enrich for molecular differences in primary PDAC, which may be clinically and biologically relevant, through examining primary tumors representing opposite spectrums of PDAC: early (localized) and late (metastatic) stage. To accomplish this, we compared nonmetastatic (UNC1) with metastatic (NEB) primary PDAC tumors. As the methods of procurement for these tumors differed, we used DWD to identify systematic biases between the two datasets [24]. This method has been used previously to successfully combine three breast cancer datasets across three microarray platforms [26], across species [34], and across multiple datasets [35],[36]. We therefore used DWD to adjust for the systematic biases between the UNC1 and NEB datasets by taking advantage of the fact that each dataset also had 15 normal pancreas samples assayed. In short, we used DWD to adjust these 15 tumor-normal pairs from both datasets to have similar distributions in principal component (PC) 1×PC 2 space. After the DWD adjustment, we used SAM to identify differentially expressed genes [22],[25]. Using a false discovery rate of 5%, we identified six genes that were differentially overexpressed between nonmetastatic and metastatic primary tumors: FBJ murine osteosarcoma viral oncogene homolog B (Fos B), Kruppel-like factor 6 (KLF6), nuclear factor of kappa light polypeptide gene enhancer in B-cells inhibitor, zeta (NFKBIZ, IKBZ, MAIL), ATPase H+/K+ exchanging, alpha polypeptide (ATP4A), germ cell associated 1 (GSG1), and sialic acid binding Ig-like lectin 11 (SIGLEC11) (Figure 1A; Table S1).
Figure 1
Figure 1
Identification, development, and application of a six-gene signature for PDAC.
Development of a Classifier Using the Six-Gene Signature
We examined the relationship of our six-gene signature to outcome using a training set of 34 patients with localized and resected PDAC. After identifying and adjusting for systematic bias using DWD [24], a resected centroid-based predictor [26] was created using the 30 samples in the derivation dataset. The centroid was then applied to the DWD-adjusted training set of primary PDAC patients to determine the performance of the six-gene signature. X-tile [27] was used to determine the optimal distance function to the centroid cut-point for classifying this training set of patients into high-risk and low-risk groups on the basis of survival (Figure 1B and 1D). The optimal cut-point occurred at a Pearson correlation coefficient of zero (p = 0.006) with patients with Pearson correlation coefficients greater than zero in the low-risk and less than zero in the high-risk groups.
Application of the Six-Gene Signature to an Independent Validation Cohort of 67 Patients
In order to evaluate the performance of the cut-point determined by X-tile [27], we applied the cut-point to an independent validation test set of 67 patients with primary PDAC. Our predetermined Pearson correlation coefficient cut-point of zero distance to the centroid successfully stratified patients into high- (n = 42) and low-risk groups (n = 25) with a median overall survival (OS) of 15 versus 49 mo (p = 0.001) (Figure 1C and 1E). Patients in the high-risk group had 1-, 2-, and 3-y estimated survival rates of 55%, 34%, and 21%, compared to 91%, 64%, and 56% in the low-risk group.
Previous studies in PDAC have found that nodal status is the most predictive of outcome for patients with localized PDAC [37]. We compared our prognostic signature to current clinical prognostic benchmarks. We found that tumors that were node positive (p = 0.091) and grade 2 or 3 trended towards a shorter survival (p = 0.080). Neither T stage (p = 0.977) nor margin status (p = 0.223) were prognostic in this cohort. Treatment with adjuvant chemotherapy (p = 0.699) or with neoadjuvant chemotherapy (p = 0.409) was also not prognostic, although only two patients received neoadjuvant chemotherapy. We found no gene expression changes between the tumors of the two patients who received neoadjuvant chemotherapy and the tumors of patients who received no treatment prior to surgery.
An important feature of any prognostic signature is that it should be independent or additive to currently used clinicopathologic prognostic criteria. We therefore compared the prognostic importance of our molecular signature in the setting of grade (p = 0.417), nodal status (p = 0.381), T stage (p = 0.675), and margin status (p = 0.295). We found that our six-gene signature was the only independent predictor of survival in the 57 patients with complete data, with a hazard ratio of 4.1 (95% confidence interval 1.7–10.0) (Table 3).
Table 3
Table 3
Cox proportional hazards regression analysis of the six-gene signature.
We also looked at whether our six-gene signature was confounded by available clinicopathological variables. We found no association between our molecular signature, and tumor size, grade, margin status, nodal status, and neoadjuvant or adjuvant chemotherapy in our independent test set (Table 4).
Table 4
Table 4
Relationship between the six-gene signature and clinicopathological variables.
KLF6 Expression in Primary PDAC
In order to further validate the six-gene signature, we performed immunohistochemical analyses for KLF6, which showed a wide range of expression values between nonmetastatic versus metastatic samples (Figure 1A). To evaluate KLF6 protein expression, we obtained another independent dataset of 50 patients represented on a TMA with matched normal, chronic pancreatitis, and PDAC (UNC2, Table 2). First, using the median score of 1.5 as the cutoff, we found that KLF6 expression was much higher in tumors compared to normal pancreas (p<0.001) (Figure 2A and 2C). KLF6 expression was strong in normal islet cells in agreement with a previous study (Figure 2Ci, white arrowhead) [38]. Second, we found that KLF6 expression with a score greater than 1.5 (high) was associated with a shorter median survival of 11 mo compared to 24 mo for patients with KLF6 expression scores less than 1.5 (low) (p = 0.04) (Figure 2B).
Figure 2
Figure 2
Significance of KLF6 and Fos B expression in primary PDAC.
We profiled and compared nonmetastatic and metastatic primary PDAC tumors and identified a six-gene signature. Although this signature was not derived on the basis of outcomes, we show that it was prognostic in a true test set of resectable PDAC patients. Importantly, our six-gene signature was independently predictive of survival, stratifying patients with median survivals of 15 compared to 49 mo, outperforming current pathological staging criteria, suggesting that our signature will be a powerful prognostic tool for patients with localized PDAC.
PDAC continues to be a devastating disease with few long-term survivors. Surgery remains the standard therapy for patients diagnosed with resectable PDAC [39]. Yet with a median survival only of less than 2 y after surgery, the attendant postoperative mortality rate of 2%–6% [40],[41], and postoperative complication and hospital readmission rates of 59% [41],[42], the decision for surgery should be made cautiously. Therefore, improved patient selection for therapy is necessary. For the majority of patients who cannot undergo surgery, gemcitabine chemotherapy remains the best option, yet only 5%–10% of patients respond to the treatment [43],[44]. Given the current therapeutic limitations, additional prognostic tools are needed to help a patient decide whether to have surgery, and/or neoadjuvant chemotherapy, or when to consider participation in a clinical trial.
Our analysis identified a surprisingly small number of genes with differential expression between early compared to late stage primary PDAC (Table S1). This finding suggests that primary PDAC may be largely homogenous from a global gene expression standpoint. Nonetheless, the differences that we identified appear to be clinically and therefore biologically important. Our findings of molecular differences in resected primary PDAC tumors suggest that there is subtle biological variation in these tumors that influences outcome. A review of previous published studies did not identify differential expression of our six genes [15],[21],[45]–[56]. This finding is not surprising, as previous studies examined differential gene expression changes between either normal pancreas or chronic pancreatitis and PDAC [15],[45]–[56]. Only one study has looked at gene expression changes between PDAC of different stages [21]. Ours was the first, to our knowledge, to study molecular differences between nonmetastatic versus metastatic primary tumors and identify and validate a prognostic signature for PDAC.
Of the six genes identified in this study, most do not have an obvious role in carcinogenesis. Three of the six genes demonstrated significantly higher expression in the poor prognostic groups (SIGLEC11, KLF6, NFKBIZ; Table S2). ATP4A, GSG1, and SIGLEC-11 have not been studied in cancer. SIGLEC-11 is thought to be expressed by tissue macrophages and also the brain microglia [57]. Interestingly, a missense mutation of SIGLEC-11 (S465A) was identified in the mutation discovery screen of the recent genome-wide sequencing of PDAC [58]. NFKBIZ, also called IkappaB zeta, binds to the p50 subunit of nuclear factor (NF)-kappaB and is important for interleukin-6 (IL-6) induction and may be induced by IL-1 receptor and Toll-like receptors [57]. Given the prevalence of chronic pancreatitis and high degree of stromal fibrosis, it is possible that NFKBIZ may play a role in PDAC and inflammation.
KLF6 is a transcription factor and its full length transcript is thought to be a tumor suppressor gene involved in prostate, lung, and ovarian carcinogenesis [59]. However a splice variant KLF6-SV1 has been shown to have oncogenic properties. The oligonucleotide probes used in the Agilent whole human genome array and the antibody against KLF6 did not differentiate between the full-length and splice variant. In agreement with a previous study [38], we found that KLF6 protein expression was higher in tumors than normal pancreas. In addition we found that higher KLF6 expression was associated with worse survival. Hartel et al. further investigated KLF6-SV1 expression in their study using real-time PCR and demonstrated that the higher KLF6 expression seen in tissues was associated with a higher ratio of KLF6-SV1 compared to full-length KLF6. Therefore our findings that KLF6 expression is higher in tumors and is prognostic is likely in agreement with this study.
Only one patient in the UNC1 cohort was treated with neoadjuvant chemotherapy compared to 80% of NEB patients who were treated with palliative chemotherapy. Although there is a possibility that our signature may be reflective of gemcitabine treatment or perhaps resistance, as NEB patients died of metastatic disease despite gemcitabine treatment, the successful application of our six-gene signature on an independent test set of patients where only 3% of patients with localized PDAC were treated with neoadjuvant therapy suggests that it is a rigorous predictor of prognosis in previously untreated patients. We found no association between our six-gene signature and whether a patient received adjuvant chemotherapy. In addition, chemotherapy treatment in this cohort, either pre- or postoperative, did not demonstrate a survival advantage.
Another concern is the validity of our hypothesis that gene expression changes at different stages of primary PDAC development may occur and be important for prognosis. Our study is in agreement with Lowe and colleagues' findings that differential gene expression changes can be identified within primary PDAC [21]. However, they did not address the prognostic value of their findings. Several studies have also suggested that gene expression changes in metastasis may be found in primary tumors. In a study of molecular differences between primary tumors and metastases, Golub and colleagues identified a gene expression signature of metastasis present that could be identified in primary tumors [60]. In addition, studies in melanoma have suggested that metastatic cells may be found in the parent primary tumor [61]. Finally studies in breast cancer have demonstrated that gene expression changes found in breast cancer cells with metastatic potential may be prognostic and predictive of patients who will develop metastasis [62]–[64]. Our study is the first to demonstrate that molecular differences in metastatic PDAC can be identified at earlier stages, and that these differences are predictive of future behavior. Whether these molecular changes are biologically associated with metastatic potential will require further investigation.
We have applied our six-gene signature to an independent dataset of 67 patients, and have validated its prognostic value. In addition, we have validated the protein expression of KLF6 in a 50-patient TMA. Although not nearly as powerful a predictor of prognosis as our six-gene signature, we found that KLF6 expression was prognostic in our 50-patient TMA. Further validation studies will be needed to see if KLF6 alone may be a useful prognostic marker as others have shown [38]. Our findings suggest that the prognostic value of KLF6 is strengthened in evaluating the six genes in their entirety.
Studies of patients with resectable PDAC demonstrate median survivals of up to 22 mo, equivalent to the median survival of patients in our training and testing cohorts [3],[11],[65]. Our finding that our six-gene signature is able to stratify patients, with startling differences in survival, suggests that it may be used to select patients for therapies. For example, for patients who are at high operative risk, knowledge of a median survival of 49 compared to 15 mo, may be helpful in the operative decision-making process. Similarly, patients who have a poor prognosis based on the six-gene signature may be considered for neoadjuvant therapy. Currently, the minority of centers use neoadjuvant therapy as a standard of care, most instead reserve this for patients with locally advanced unresectable or borderline resectable tumors. Therefore the current decision-making process is based on anatomical considerations. Our prognostic signature may refine this paradigm such that neoadjuvant therapy is offered to patients on the basis of biological considerations, regardless of resectability, and may allow us to further study and maximize the benefits of neoadjuvant treatment. In addition, as new therapies are developed, it may help to determine whether patients may require more or less aggressive treatment. Finally, our findings that there are molecular differences associated with late-stage primary tumors, which translate into differences in prognosis, suggest that the six genes in this signature should be further studied for their potential as biomarkers, and some of these genes, or the pathways that they fall into, may represent new therapeutic targets.
Supporting Information
Table S1
SAM of metastatic compared to localized primary tumors.
(1.75 MB XLS)
Click here for additional data file.(1.6M, xls)
Table S2
Comparison of individual genes in the six-gene signature between the high- and low-risk groups.
(0.03 MB XLS)
Click here for additional data file.(26K, xls)
Acknowledgments
JJY thanks Ling Li, the UNC Genomics and Bioinformatics, Translational Pathology, and Tissue Procurement Core Laboratory Facilities for superb technical assistance and Joel Tepper and William Kim for critical review of the manuscript.
Abbreviations
DWDdistance weighted discrimination
NEBUniversity of Nebraska Medical Center Rapid Autopsy Pancreatic Program
PDACpancreatic ductal adenocarcinoma
SAMsignificance analysis of microarrays
TMAtissue microarray
TNMtumor, node, and metastasis
UNCUniversity of North Carolina at Chapel Hill

Footnotes
The authors have declared that no competing interests exist.
This work was supported by the Emerald Foundation (JJY), the Kimmel Foundation (JJY), the American College of Surgeons (JJY), the American Surgical Association Foundation (JJY), the North Carolina University Cancer Research Fund (JJY), and the NCI SPORE in Gastrointestinal Cancer P50 CA106991 (jjy, cjd), NCI SPORE in SPORE in Breast Cancer - Core 3 - Biostatistics 5-P50-CA58223-15 (JSM), NIH T32 Cancer Biology Training Grant (JKS), NIH CA106610 (CAI-D). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
References
1. Boyle P, Levin B, editors. Lyon: International Agency for Research on Cancer; 2008. World Cancer Report 2008.
2. Parkin DM, Bray F, Ferlay J, Pisani P. Global cancer statistics, 2002. CA Cancer J Clin. 2005;55:74–108. [PubMed]
3. Neuhaus P, Riess H, Post S, Gellert K, Ridwelski K, et al. CONKO-001: final results of the randomized, prospective, multicenter phase III trial of adjuvant chemotherapy with gemcitabine versus observation in patients with resected pancreatic cancer (PC). Journal of Clinical Oncology, 2008 ASCO Annual Meeting Proceedings. 2008;26:LBA4504.
4. Ferrone CR, Brennan MF, Gonen M, Coit DG, Fong Y, et al. Pancreatic adenocarcinoma: the actual 5-year survivors. J Gastrointest Surg. 2008;12:701–706. [PubMed]
5. Conlon KC, Klimstra DS, Brennan MF. Long-term survival after curative resection for pancreatic ductal adenocarcinoma. Clinicopathologic analysis of 5-year survivors. Ann Surg. 1996;223:273–279. [PubMed]
6. Schnelldorfer T, Ware AL, Sarr MG, Smyrk TC, Zhang L, et al. Long-term survival after pancreatoduodenectomy for pancreatic adenocarcinoma: is cure possible? Ann Surg. 2008;247:456–462. [PubMed]
7. Cleary SP, Gryfe R, Guindi M, Greig P, Smith L, et al. Prognostic factors in resected pancreatic adenocarcinoma: analysis of actual 5-year survivors. J Am Coll Surg. 2004;198:722–731. [PubMed]
8. Han SS, Jang JY, Kim SW, Kim WH, Lee KU, et al. Analysis of long-term survivors after surgical resection for pancreatic cancer. Pancreas. 2006;32:271–275. [PubMed]
9. Winter JM, Cameron JL, Campbell KA, Arnold MA, Chang DC, et al. 1423 pancreaticoduodenectomies for pancreatic cancer: a single-institution experience. J Gastrointest Surg. 2006;10:1199–1210; discussion 1210-1191. [PubMed]
10. Ahmad NA, Lewis JD, Ginsberg GG, Haller DG, Morris JB, et al. Long term survival after pancreatic resection for pancreatic adenocarcinoma. Am J Gastroenterol. 2001;96:2609–2615. [PubMed]
11. Ueno H, Kosuge T, Matsuyama Y, Yamamoto J, Nakao A, et al. A randomised phase III trial comparing gemcitabine with surgery-only in patients with resected pancreatic cancer: Japanese Study Group of Adjuvant Therapy for Pancreatic Cancer. Br J Cancer. 2009;101:908–915. [PMC free article] [PubMed]
12. Evans DB, Varadhachary GR, Crane CH, Sun CC, Lee JE, et al. Preoperative gemcitabine-based chemoradiation for patients with resectable adenocarcinoma of the pancreatic head. J Clin Oncol. 2008;26:3496–3502. [PubMed]
13. Varadhachary GR, Wolff RA, Crane CH, Sun CC, Lee JE, et al. Preoperative gemcitabine and cisplatin followed by gemcitabine-based chemoradiation for resectable adenocarcinoma of the pancreatic head. J Clin Oncol. 2008;26:3487–3495. [PubMed]
14. Kolbert CP, Chari S, Sreekumar R. Microarray technologies for gene transcript analysis in pancreatic cancer. Technol Cancer Res Treat. 2008;7:55–59. [PubMed]
15. Grutzmann R, Boriss H, Ammerpohl O, Luttges J, Kalthoff H, et al. Meta-analysis of microarray data on pancreatic cancer defines a set of commonly dysregulated genes. Oncogene. 2005;24:5079–5088. [PubMed]
16. Grutzmann R, Saeger HD, Luttges J, Schackert HK, Kalthoff H, et al. Microarray-based gene expression profiling in pancreatic ductal carcinoma: status quo and perspectives. Int J Colorectal Dis. 2004;19:401–413. [PubMed]
17. Goggins M. Identifying molecular markers for the early detection of pancreatic neoplasia. Semin Oncol. 2007;34:303–310. [PMC free article] [PubMed]
18. Grote T, Logsdon CD. Progress on molecular markers of pancreatic cancer. Curr Opin Gastroenterol. 2007;23:508–514. [PubMed]
19. Tonini G, Pantano F, Vincenzi B, Gabbrielli A, Coppola R, et al. Molecular prognostic factors in patients with pancreatic cancer. Expert Opin Ther Targets. 2007;11:1553–1569. [PubMed]
20. Yeh JJ. Prognostic signature for pancreatic cancer: are we close? Future Oncol. 2009;5:313–321. [PubMed]
21. Kim HN, Choi DW, Lee KT, Lee JK, Heo JS, et al. Gene expression profiling in lymph node-positive and lymph node-negative pancreatic cancer. Pancreas. 2007;34:325–334. [PubMed]
22. Yang YH, Dudoit S, Luu P, Lin DM, Peng V, et al. Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002;30:e15. [PMC free article] [PubMed]
23. Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, et al. Missing value estimation methods for DNA microarrays. Bioinformatics. 2001;17:520–525. [PubMed]
24. Benito M, Parker J, Du Q, Wu J, Xiang D, et al. Adjustment of systematic microarray data biases. Bioinformatics. 2004;20:105–114. [PubMed]
25. 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]
26. Hu Z, Fan C, Oh DS, Marron JS, He X, et al. The molecular portraits of breast tumors are conserved across microarray platforms. BMC Genomics. 2006;7:96. [PMC free article] [PubMed]
27. Camp RL, Dolled-Filhart M, Rimm DL. X-tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin Cancer Res. 2004;10:7252–7259. [PubMed]
28. Yeh JJ, Routh ED, Rubinas T, Peacock J, Martin TD, et al. KRAS/BRAF mutation status and ERK1/2 activation as biomarkers for MEK1/2 inhibitor therapy in colorectal cancer. Mol Cancer Ther. 2009;8:834–843. [PMC free article] [PubMed]
29. Hoos A, Cordon-Cardo C. Tissue microarray profiling of cancer specimens and cell lines: opportunities and limitations. Lab Invest. 2001;81:1331–1338. [PubMed]
30. Troester MA, Lee MH, Carter M, Fan C, Cowan DW, et al. Activation of host wound responses in breast cancer microenvironment. Clin Cancer Res. 2009;15:7020–7028. [PMC free article] [PubMed]
31. Allinen M, Beroukhim R, Cai L, Brennan C, Lahti-Domenici J, et al. Molecular characterization of the tumor microenvironment in breast cancer. Cancer Cell. 2004;6:17–32. [PubMed]
32. Mueller MM, Fusenig NE. Friends or foes - bipolar effects of the tumour stroma in cancer. Nat Rev Cancer. 2004;4:839–849. [PubMed]
33. Comoglio PM, Trusolino L. Cancer: the matrix is now in control. Nat Med. 2005;11:1156–1159. [PubMed]
34. Herschkowitz JI, Simin K, Weigman VJ, Mikaelian I, Usary J, 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]
35. Oh DS, Troester MA, Usary J, Hu Z, He X, et al. Estrogen-regulated genes predict survival in hormone receptor-positive breast cancers. J Clin Oncol. 2006;24:1656–1664. [PubMed]
36. Lu Y, Lemon W, Liu PY, Yi Y, Morrison C, et al. A gene expression signature predicts survival of patients with stage I non-small cell lung cancer. PLoS Med. 2006;3:e467. doi: 10.1371/journal.pmed.0030467. [PMC free article] [PubMed]
37. Sohn TA, Yeo CJ, Cameron JL, Koniaris L, Kaushal S, et al. Resected adenocarcinoma of the pancreas-616 patients: results, outcomes, and prognostic indicators. J Gastrointest Surg. 2000;4:567–579. [PubMed]
38. Hartel M, Narla G, Wente MN, Giese NA, Martignoni ME, et al. Increased alternative splicing of the KLF6 tumour suppressor gene correlates with prognosis and tumour grade in patients with pancreatic cancer. Eur J Cancer. 2008;44:1895–1903. [PubMed]
39. Yeo CJ, Cameron JL, Sohn TA, Lillemoe KD, Pitt HA, et al. Six hundred fifty consecutive pancreaticoduodenectomies in the 1990s: pathology, complications, and outcomes. Ann Surg. 1997;226:248–257; discussion 257-260. [PubMed]
40. Eppsteiner RW, Csikesz NG, McPhee JT, Tseng JF, Shah SA. Surgeon volume impacts hospital mortality for pancreatic resection. Ann Surg. 2009;249:635–640. [PubMed]
41. Yermilov I, Bentrem D, Sekeris E, Jain S, Maggard MA, et al. Readmissions following pancreaticoduodenectomy for pancreas cancer: a population-based appraisal. Ann Surg Oncol. 2009;16:554–561. [PubMed]
42. DeOliveira ML, Winter JM, Schafer M, Cunningham SC, Cameron JL, et al. Assessment of complications after pancreatic surgery: a novel grading system applied to 633 patients undergoing pancreaticoduodenectomy. Ann Surg. 2006;244:931–937; discussion 937-939. [PubMed]
43. Abou-Alfa GK, Letourneau R, Harker G, Modiano M, Hurwitz H, et al. Randomized phase III study of exatecan and gemcitabine compared with gemcitabine alone in untreated advanced pancreatic cancer. J Clin Oncol. 2006;24:4441–4447. [PubMed]
44. Van Cutsem E, Vervenne WL, Bennouna J, Humblet Y, Gill S, et al. Phase III trial of bevacizumab in combination with gemcitabine and erlotinib in patients with metastatic pancreatic cancer. J Clin Oncol. 2009;27:2231–2237. [PubMed]
45. Crnogorac-Jurcevic T, Efthimiou E, Capelli P, Blaveri E, Baron A, et al. Gene expression profiles of pancreatic cancer and stromal desmoplasia. Oncogene. 2001;20:7437–7446. [PubMed]
46. Crnogorac-Jurcevic T, Efthimiou E, Nielsen T, Loader J, Terris B, et al. Expression profiling of microdissected pancreatic adenocarcinomas. Oncogene. 2002;21:4587–4594. [PubMed]
47. Karanjawala ZE, Illei PB, Ashfaq R, Infante JR, Murphy K, et al. New markers of pancreatic cancer identified through differential gene expression analyses: claudin 18 and annexin A8. Am J Surg Pathol. 2008;32:188–196. [PMC free article] [PubMed]
48. Lowe AW, Olsen M, Hao Y, Lee SP, Taek Lee K, et al. Gene expression patterns in pancreatic tumors, cells and tissues. PLoS ONE. 2007;2:e323. doi: 10.1371/journal.pone.0000323. [PMC free article] [PubMed]
49. Ishikawa M, Yoshida K, Yamashita Y, Ota J, Takada S, et al. Experimental trial for diagnosis of pancreatic ductal carcinoma based on gene expression profiles of pancreatic ductal cells. Cancer Sci. 2005;96:387–393. [PubMed]
50. Segara D, Biankin AV, Kench JG, Langusch CC, Dawson AC, et al. Expression of HOXB2, a retinoic acid signaling target in pancreatic cancer and pancreatic intraepithelial neoplasia. Clin Cancer Res. 2005;11:3587–3596. [PubMed]
51. Logsdon CD, Simeone DM, Binkley C, Arumugam T, Greenson JK, et al. Molecular profiling of pancreatic adenocarcinoma and chronic pancreatitis identifies multiple genes differentially regulated in pancreatic cancer. Cancer Res. 2003;63:2649–2657. [PubMed]
52. Iacobuzio-Donahue CA, Ashfaq R, Maitra A, Adsay NV, Shen-Ong GL, et al. Highly expressed genes in pancreatic ductal adenocarcinomas: a comprehensive characterization and comparison of the transcription profiles obtained from three major technologies. Cancer Res. 2003;63:8614–8622. [PubMed]
53. Iacobuzio-Donahue CA, Maitra A, Olsen M, Lowe AW, van Heek NT, et al. Exploration of global gene expression patterns in pancreatic adenocarcinoma using cDNA microarrays. Am J Pathol. 2003;162:1151–1162. [PubMed]
54. Iacobuzio-Donahue CA, Maitra A, Shen-Ong GL, van Heek T, Ashfaq R, et al. Discovery of novel tumor markers of pancreatic cancer using global gene expression technology. Am J Pathol. 2002;160:1239–1249. [PubMed]
55. Crnogorac-Jurcevic T, Missiaglia E, Blaveri E, Gangeswaran R, Jones M, et al. Molecular alterations in pancreatic carcinoma: expression profiling shows that dysregulated expression of S100 genes is highly prevalent. J Pathol. 2003;201:63–74. [PubMed]
56. Grutzmann R, Pilarsky C, Ammerpohl O, Luttges J, Bohme A, et al. Gene expression profiling of microdissected pancreatic ductal carcinomas using high-density DNA microarrays. Neoplasia. 2004;6:611–622. [PMC free article] [PubMed]
57. Angata T, Kerr SC, Greaves DR, Varki NM, Crocker PR, et al. Cloning and characterization of human Siglec-11. A recently evolved signaling that can interact with SHP-1 and SHP-2 and is expressed by tissue macrophages, including brain microglia. J Biol Chem. 2002;277:24466–24474. [PubMed]
58. Jones S, Zhang X, Parsons DW, Lin JC, Leary RJ, et al. Core signaling pathways in human pancreatic cancers revealed by global genomic analyses. Science. 2008;321:1801–1806. [PMC free article] [PubMed]
59. DiFeo A, Martignetti JA, Narla G. The role of KLF6 and its splice variants in cancer therapy. Drug Resist Updat. 2009;12:1–7. [PubMed]
60. Ramaswamy S, Ross KN, Lander ES, Golub TR. A molecular signature of metastasis in primary solid tumors. Nat Genet. 2003;33:49–54. [PubMed]
61. Fidler IJ, Kripke ML. Metastasis results from preexisting variant cells within a malignant tumor. Science. 1977;197:893–895. [PubMed]
62. Minn AJ, Gupta GP, Padua D, Bos P, Nguyen DX, et al. Lung metastasis genes couple breast tumor size and metastatic spread. Proc Natl Acad Sci U S A. 2007;104:6740–6745. [PubMed]
63. Minn AJ, Gupta GP, Siegel PM, Bos PD, Shu W, et al. Genes that mediate breast cancer metastasis to lung. Nature. 2005;436:518–524. [PMC free article] [PubMed]
64. Minn AJ, Kang Y, Serganova I, Gupta GP, Giri DD, et al. Distinct organ-specific metastatic potential of individual breast cancer cells and primary tumors. J Clin Invest. 2005;115:44–55. [PMC free article] [PubMed]
65. Neoptolemos JP, Dunn JA, Stocken DD, Almond J, Link K, et al. Adjuvant chemoradiotherapy and chemotherapy in resectable pancreatic cancer: a randomised controlled trial. Lancet. 2001;358:1576–1585. [PubMed]
Articles from PLoS Medicine are provided here courtesy of
Public Library of Science