Research article | Open | Open Peer Review | Published:
Metabolic markers GAPDH, PKM2, ATP5B and BEC-index in advanced serous ovarian cancer
BMC Clinical Pathologyvolume 13, Article number: 30 (2013)
A deregulated energy metabolism is a hallmark of malignant disease that offers possible future targets for treatment. We investigated the prognostic value of the glycolytic enzymes glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and pyruvate kinase type M2 (PKM2), mitochondrial β-F1-ATPase (ATP5B) and the bioenergetic cellular (BEC) index in advanced ovarian cancer.
Fresh tumor samples were prospectively collected from 123 patients undergoing primary surgery for suspected advanced ovarian cancer. Of these, 57 met the eligibility criteria; stage IIC-IV, serous or endometrioid subtype, specimens containing ≥ 50% tumor cells and patients receiving platinum-based chemotherapy. An adequate amount of mRNA could be extracted in all but one case, with a resultant study population of 56 patients. Eighty-six percent of cases had serous tumors, and 93% were grade 2–3. GAPDH, PKM2 and ATP5B mRNA- and protein expression was assessed by real-time PCR and immunohistochemistry. We estimated the association with platinum-free interval (PFI) and overall survival (OS) by Cox proportional hazards models. Median follow-up was 60 months.
High GAPDH mRNA levels (HR 2.1, 95% CI 1.0-4.5) and low BEC-index (HR 0.47, 95% CI 0.23-0.95) were both independently associated with shorter PFI. Median PFI for patients with high GAPDH mRNA was 5.0 months compared to 10.1 months for low expression cases (p = 0.031). Similarly, median PFI for patients with low BEC-index based on mRNA was 5.3 months compared to 9.8 months for high BEC-index cases (p = 0.028).
High GAPDH or low BEC-index mRNA expression indicate early disease progression in advanced serous ovarian cancer.
Epithelial ovarian cancer (EOC) is a heterogeneous disease comprising several subtypes with very different response to treatment and prognosis. The majority of EOC-patients are detected at advanced stage and the overall 5-year survival rate is 50% or less . Treatment for primary EOC consists of cytoreductive surgery and platinum-taxane combination chemotherapy. Despite emerging knowledge of subtype-specific variations in tumor biology and chemo-responsiveness, EOC is still treated as one entity. The high-grade serous carcinomas account for 60–70% of diagnosed tumors. This subtype is, although often initially platinum-sensitive, associated with a poor prognosis. Identification of EOC subtype-specific prognostic markers is important for development of new targeted therapies.
A reprogrammed energy metabolism has recently been designated a “new” hallmark of cancer . These metabolic changes, including increased aerobic glycolysis, decreased oxidative phosphorylation and increased glutaminolysis and fatty acid-metabolism, are known features of cancer, which often correlate with an overall poor prognosis [3, 4]. Several glycolytic enzymes, for example glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and pyruvate kinase M2 (PKM2), have been reported overexpressed in cancer [5–8]. First introduced by Cuezva and co-workers in 2002 , a bioenergetic cellular index (BEC-index), based on the protein levels of mitochondrial β-F1-ATPase (ATP5B, the catalytic subunit of the rate-limiting enzyme of oxidative phosphorylation), heat-shock protein 60 (HSP60, a mitochondrial chaperone), and glycolytic GAPDH, has been shown to be of prognostic value in several cancer types such as colon, lung and breast carcinomas [3, 9, 10]. The expression of above metabolic markers has not yet been studied in EOC.
We have previously shown, in cell lines and ovarian cancer cells from ascites, that inhibitors of glycolysis can potentiate the effects of different chemotherapeutics, including platinum compounds [11, 12]. In the present study we investigated the mRNA and protein expression of GAPDH, PKM2 and ATP5B in 56 advanced serous or endometrioid ovarian cancers, to test their ability to predict length of platinum-free interval and survival. We also assessed, at the mRNA-level, the prognostic significance of the BEC-index.
We prospectively collected fresh tumor specimens from patients undergoing surgery for epithelial ovarian, fallopian tube or primary peritoneal cancer. Patients were operated at Karolinska University Hospital, Stockholm, Sweden, between April 2003 and July 2008. They all gave informed consent and the Ethics Committee at Karolinska Institutet, Stockholm, approved the study.
Tumor and, if present, ascites samples were obtained from 123 patients of which 62 fulfilled the study’s inclusion criteria, i.e. age ≥ 18 years, FIGO stage IIC-IV serous or endometrioid ovarian, fallopian tube or primary peritoneal cancer (Figure 1). Also, patients had to be scheduled for platinum-based chemotherapy. Tumors were classified and graded according to WHO standards . Patient baseline characteristics, treatment and relapse data were prospectively collected in Case Report Forms (CRFs). All but four patients underwent primary debulking surgery aiming at cytoreduction to microscopic disease. In the remaining cases only biopsy sampling was done. At the time of study, the staging procedure for ovarian cancer did not include lymphadenectomy.
Tumor sample collection, processing and RNA extraction
The collection and processing of fresh tumor biopsies has been previously described in detail . For solid tumor specimens, imprints were prepared for pathological review. Only samples with ≥ 50% tumor cells were included in the study, and five of 62 cases were thus excluded. We did smear preparations to assess tumor cell yield from ascites, with the same cut-off of 50% tumor cells as for biopsies. Ascites samples were processed using the Lymphoprep™ (Axis-Shield, Oslo, Norway) density gradient separation method  for enrichment of tumor cells. We used the RNeasy Midi Kit™ (Qiagen, Hilden, Germany) for RNA extraction. RNA quality was assessed using Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA). An adequate quantity (≥1.5 μg) and quality (RIN 8–10) RNA could be extracted in all but one case. Thus, the present study population consists of 56 patients (Figure 1). In five cases the RNA yield from the ascites sample was better than that from the corresponding solid tumor, and therefore used in the analysis.
RNA was reverse transcribed by Superscript III (Invitrogene AB, Stockholm, Sweden). We used the ABI 7500 computerized system (Applied Biosystems, Foster City, CA) and SYBR Green™ for amplification reactions. Relative quantity expression values were calculated by the ΔΔCt method using GenEx software (MultiD, Gothenburg, Sweden). Primer sequences are listed in Additional file 1. HPRT1 (hypoxanthine guanine phosphoribosyltransferase 1) and Β2-microglobulin were used as reference genes. Universal Human Reference RNA (Stratagene, Santa Clara, CA) served as positive control. Due to technical failure in one case, PCR analyses were performed in 55 of 56 tumor samples.
For analyses of GAPDH, PKM2 and ATP5B real-time PCR results, we divided cases into three equally sized groups according to their relative quantity of mRNA expression, and the group with lowest expression (one third of patients, n = 18) was then compared to cases with higher expression (two thirds of patients, n = 37). A cut-off at one third was chosen because it approximates the proportion of patients still alive at median follow-up. The BEC-index was calculated, at the mRNA-level, by dividing the ATP5B:HSP60 ratio with the GAPDH expression value. For this, we used recently published real-time PCR data for HSP60 mRNA . The third of cases with the highest BEC-index value (n = 19) was compared to cases with lower values (n = 36).
Formalin-fixed, paraffin-embedded tumor blocks could be obtained in 54 of 56 cases. Immunohistochemistry was performed as previously described . In brief, following heat-induced antigen retrieval, primary antibodies (Additional file 2) were incubated at 8° overnight. Reaction products were visualized using the Vectastain Elite ABC kit (Vector Laboratories, Burlingame, CA). Negative controls were produced by omission of the primary antibody. A breast carcinoma specimen served as positive control for GAPDH and PKM2, and for ATP5B a colon carcinoma sample was used.
Three observers (EH, SE and JC), blinded for clinical data, independently evaluated all slides by assessing the whole tumor area. The maximum staining intensity of tumor cells was scored 0-3+, and the percentage of cancer cells thus stained estimated. Cases were dichotomized into low and high expression groups. Cut-offs were set arbitrarily to make groups equally sized; for GAPDH: low expression, 1+ or < 50% of tumor cells staining 2+, and high expression, 3+ or ≥ 50% staining 2+, for PKM2: low, 1-2+ or < 20% of tumor cells staining 3+, and high, ≥ 20% staining 3+, and for ATP5B: low, staining intensity of tumor cells 1-2+, and high, staining intensity 3+ (regardless of proportion). Discrepancies between observers were found for 13% of examined GAPDH-slides, 6% of PKM2-slides and 15% of ATP5B-slides, in which cases consensus was reached on further review.
Statistical analyses were performed using STATA 11.2 software. The platinum-free interval (PFI) was defined as the time interval from last course of primary platinum-based chemotherapy (end of treatment, EOT) to documented disease progression, death, or last follow-up. Overall survival (OS) was defined as the time from diagnosis to death or last follow-up. We estimated PFI and OS by the Kaplan-Meier method  and curves were compared using Wald’s test. Cox proportional hazards regression models , adjusted for standard confounding risk factors age, stage, grade and postoperative residual tumor, were used for multivariate analyses. P-values below 0.05 were considered statistically significant.
Patients were monitored until 31 March 2011 with a median follow-up of 59.9 (30.2-82.8) months. Clinical and pathological characteristics are presented in Table 1. Over 85% of patients had poorly differentiated (grade 3) carcinomas, and 89% had FIGO stage IIIC-IV disease. In 70% of cases the postoperative tumor residuals measured at least 10 mm. All patients received platinum-based chemotherapy. Of 49 evaluable patients (with measurable disease at start of chemotherapy), 82% responded to treatment. Relapse within 6 months from EOT was detected in 25 of 56 cases (45%). Median PFI was 7.2 months, and median OS 34.6 months. At study closure, 19 patients (34%) were alive, with no evidence of disease in 9 cases (16%).
Eighteen patients had a tumor GAPDH mRNA expression of less than 0.25, and 37 cases expressed at least 0.25 relative to reference genes and positive control. Kaplan-Meier curves for PFI and OS, based on GAPDH mRNA expression, are shown in Figure 2A and B. The median PFI was significantly shorter in the group with high expression compared to the low expression group (5.0 and 10.1 months, respectively, p = 0.031). Table 2 shows uni- and multivariate analyses of PFI and OS. In univariate analysis, high GAPDH mRNA expression was associated with shorter PFI (HR 2.1, 95% CI 1.1-4.0, p = 0.031) and OS (HR 2.8, 95% CI 1.2-6.5, p = 0.015). In multivariate analysis, high GAPDH mRNA remained significant for shorter PFI (HR 2.1, 95% CI 1.0-4.5, p = 0.043) but not for OS (HR 2.3, 95% CI 1.0-5.6, p = 0.063).
GAPDH protein expression by IHC was low in 26 and high in 28 of 54 cases (Table 1, Figure 3). GAPDH-reactivity was predominantly cytoplasmatic, but in 22 cases (41%) we also found nuclear localization of the enzyme. There were no statistically significant differences in PFI or survival between groups with high- or low GAPDH protein expression (data not shown).
Eighteen patients had a tumor PKM2 mRNA expression of less than 0.70, and 37 cases had higher mRNA expression. There were no statistically significant differences in PFI or survival between groups with high or low mRNA expression (Table 2).
Twenty-six cases had low PKM2 protein expression and in 28 cases the expression was high (Table 1, Figure 3). The assessed PKM2-reactivity was cytoplasmatic, but in 11 cases (20%) nuclear localization of the protein was observed in mitotic cells. The immunoreaction does not provide any information whether the enzyme present is in its dimeric or more glycolytically active tetrameric form. There were no statistically significant differences in PFI or survival between high and low PKM2 protein expressing groups.
Eighteen patients had a tumor ATP5B mRNA expression of less than 0.31, and in 37 cases the mRNA expression was higher. Uni- and multivariate analyses of PFI and OS are shown in Table 2. In univariate analysis, high ATP5B mRNA expression predicted poor OS (HR 2.6, 95% CI 1.1-5.9, p = 0.025), but it did not remain significant in multivariate analysis (HR 2.3, 95% CI 1.0-5.3, p = 0.062).
Twenty-six cases had low ATP5B protein expression and in 28 cases the expression was high (Table 1, Figure 3). ATP5B-reactivity was cytoplasmatic and often granular. In univariate analysis, a high ATP5B protein expression predicted short PFI (HR 1.9, 95% CI 1.0-3.4, p = 0.039), but it did not reach significance in the multivariate model (data not shown).
Thirty-six patients had an mRNA-based tumor BEC-index of ≤ 2.65 and 19 cases had higher expression. Kaplan-Meier curves for PFI and OS, based on BEC-index mRNA expression, are shown in Figure 2C and D. The median PFI was significantly longer in the group with high compared to the group with low BEC-index (9.8 and 5.3 months, respectively, p = 0.028). Table 2 shows uni- and multivariate analyses for PFI and OS. Multivariate analysis showed a high BEC-index mRNA to independently predict longer PFI (HR 0.47, 95% CI 0.23-0.95, p = 0.035). Also, univariate analysis indicated high BEC-index to predict longer survival (HR 0.42, 95% CI 0.20-0.94, p = 0.033), but it did not remain significant in multivariate calculations.
Our data indicate that high GAPDH as well as low BEC-index based on mRNA are associated with early disease progression in patients with advanced high-grade serous adenocarcinomas of the ovary, fallopian tube or peritoneum. To our knowledge, no other study has investigated the prognostic impact of these markers in ovarian cancer.
GAPDH, a tetramer of four identical 37-kDa subunits, is a key glycolytic enzyme present in all tissues. It has commonly been used as a house-keeping reference gene/protein for expression analyses. However, this practice has recently been challenged due to findings of GAPDH playing a role in cancer pathogenesis. In a recent review GAPDH expression was reported deregulated in e.g. cancers of the lung, breast, colon, pancreas, liver, kidney and prostate .
In agreement with our results, Révillion and co-workers reported reduced OS and relapse-free survival among breast cancer patients with enhanced GAPDH mRNA expression as assessed by real-time PCR . The study included 404 unselected breast cancer patients with a median follow-up of 82 months. In multivariate analyses, GAPDH was not an independent prognostic factor, but its expression was inversely correlated to oestradiol and progesterone receptor concentrations, young age and grading. Thus, GAPDH mRNA expression in breast cancer seems to reflect tumor aggressiveness. In the present study, GAPDH mRNA expression was an independent marker for early tumor progression in a cohort of ovarian cancer patients with poor prognosis. Also, our results demonstrate that GAPDH should not be used as reference RNA in ovarian cancer.
In contrast to a recently published study in breast cancer we found no association between GAPDH protein expression and platinum-free interval or OS. Descotes et al. used 2D-electrophoresis and Western blotting in 41 node-negative breast cancers of which 20 had metastatic relapse and 21 had no metastatic recurrence . The former group had significantly higher GAPDH expression. Both our studies are based on small cohorts and should be considered as exploratory. In addition to tumor type, the methodological approaches also differed. There are several possible biological explanations to the apparent discrepancy between our GAPDH mRNA and protein data, e. g. posttranslational regulation. A plausible reason can also be the limited sensitivity of immunohistochemistry, which requires major differences in expression for visual detection.
Taking both glycolytic and oxidative capacity into account, a low BEC-index has been shown to predict shorter survival in colon, lung and breast carcinomas [3, 9, 10]. Although we used mRNA instead of protein expression data, the results of the present study in ovarian cancer are in accordance with these reports. Apart from GAPDH and ATP5B, the BEC-index is also influenced by the expression of the mitochondrial multifunctional chaperone HSP60. We have earlier shown HSP60 overexpression to predict poor outcome in serous ovarian cancer , and the BEC-index data from the present study probably reflects both the tumor material’s GAPDH and HSP60 expression.
ATP5B-expression per se has also been investigated, but the results seem contradictory. At the protein level, high ATP5B has been shown to predict longer survival in colon cancer [3, 20]. In contrast, it was found overexpressed  and associated with poor survival  in breast cancer. In the present study, high ATP5B mRNA expression in ovarian cancer was associated with worse OS, but was not found to have independent prognostic importance in multivariate analysis (p = 0.062).
PKM2 has been linked to poor survival in signet ring cell cancer and to advanced stage in colorectal cancer [7, 8]. However, the results from our work indicate that PKM2 might not provide prognostic information in ovarian cancer.
The mechanisms by which the deregulated metabolic enzymes in cancer mediate the observed poor prognosis are complex. Rapidly proliferating cells have a high demand for energy and glycolytic intermediates for anabolic processes. Some glycolytic enzymes also directly participate in pathways of proliferation and apoptosis [22–24]. In addition, PKM2 has been implicated in platinum-resistance in ovarian cancer cell lines , and down-regulation of ATP5B protein expression has in vitro been shown to induce 5-FU-resistance in colon cancer . However, targeting breast cancer cell lines with ATP5B-inhibitor Aurovertin B decreases proliferation .
We present a homogeneous clinical material of advanced, mainly high-grade serous cancers, in which all patients had peritoneal spread. In this uniform cohort with dismal prognosis, GAPDH and BEC-index mRNA could distinguish groups with different time to relapse. In line with this, other studies have supported the notion that there are different subpopulations of high-grade serous carcinomas [14, 27].
This study is limited by its small size, so our results cannot be generalized to all patients with advanced serous carcinomas. Its strengths are the prospective study design and application of CRFs for collection of data. We also used fresh tumor samples to minimize analytical errors and our statistical analyses were adjusted for known risk factors.
Our findings indicate that GAPDH and BEC-index may identify groups of advanced high-grade serous carcinomas with different prognosis. However, before transferring these results into clinical practice, validation in large independent cohorts is needed. Targeting metabolic enzymes could be a future strategy to improve survival in this group of patients with poor prognosis.
Heintz APM, Odicino F, Maisonneuve P, Quinn MA, Benedet JL, Creasman WT, Ngan HYS, Pecorelli S, Beller U: Carcinoma of the Ovary. Int J Gynecol Obstet. 2006, 95 (Supplement 1): S161-S192.
Hanahan D, Weinberg Robert A: Hallmarks of Cancer: The Next Generation. Cell. 2011, 144 (5): 646-674.
Cuezva JM, Krajewska M, de Heredia ML, Krajewski S, Santamaria G, Kim H, Zapata JM, Marusawa H, Chamorro M, Reed JC: The Bioenergetic Signature of Cancer: A Marker of Tumor Progression. Cancer Res. 2002, 62 (22): 6674-6681.
DeBerardinis RJ, Lum JJ, Hatzivassiliou G, Thompson CB: The Biology of Cancer: Metabolic Reprogramming Fuels Cell Growth and Proliferation. Cell Metab. 2008, 7 (1): 11-20.
Altenberg B, Greulich KO: Genes of glycolysis are ubiquitously overexpressed in 24 cancer classes. Genomics. 2004, 84 (6): 1014-1020.
Descotes F, Jézéquel P, Spyratos F, Campion L, Grenot C, Lerebours F, Campone M, Guérin-Charbonnel C, Lanoe D, Adams M, et al.: Identification of potential prognostic biomarkers for node-negative breast tumours by proteomic analysis: A multicentric 2004 national PHRC study. Int J Oncol. 2012, 41 (1): 92-104.
Lim JY, Yoon SO, Seol SY, Hong SW, Kim JW, Choi SH, Cho JY: Overexpression of the M2 isoform of pyruvate kinase is an adverse prognostic factor for signet ring cell gastric cancer. World J Gastroenterol. 2012, 18 (30): 4037-4043.
Zhou C-F, Li X-B, Sun H, Zhang B, Han Y-S, Jiang Y, Zhuang Q-L, Fang J, Wu G-H: Pyruvate kinase type M2 is upregulated in colorectal cancer and promotes proliferation and migration of colon cancer cells. IUBMB Life. 2012, 64 (9): 775-782.
Cuezva JM, Chen G, Alonso AM, Isidoro A, Misek DE, Hanash SM, Beer DG: The bioenergetic signature of lung adenocarcinomas is a molecular marker of cancer diagnosis and prognosis. Carcinogenesis. 2004, 25 (7): 1157-1163.
Isidoro A, Casado E, Redondo A, Acebo P, Espinosa E, Alonso AM, Cejas P, Hardisson D, Fresno Vara JA, Belda-Iniesta C, et al.: Breast carcinomas fulfill the Warburg hypothesis and provide metabolic markers of cancer prognosis. Carcinogenesis. 2005, 26 (12): 2095-2104.
Hernlund E, Hjerpe E, Åvall-Lundqvist E, Shoshan M: Ovarian carcinoma cells with low levels of β-F1-ATPase are sensitive to combined platinum and 2-deoxy-d-glucose treatment. Mol Cancer Ther. 2009, 8 (7): 1916-1923.
Hernlund E, Ihrlund LS, Khan O, Ates YO, Linder S, Panaretakis T, Shoshan MC: Potentiation of chemotherapeutic drugs by energy metabolism inhibitors 2-deoxyglucose and etomoxir. Int J Cancer. 2008, 123 (2): 476-483.
Tavassoéli F, Devilee P: World Health Organization Classification of Tumours. 2003, Lyon: IARC Press
Hjerpe E, Egyhazi S, Carlson J, Frostvik Stolt M, Shedvins K, Johansson H, Shoshan M, Åvall-Lundqvist E: HSP60 predicts survival in advanced serous ovarian cancer. Int J Gynecol Cancer. 2013, 23 (3): 448-455.
Bjorge L, Hakulinen J, Vintermyr OK, Jarva H, Jensen TS, Iversen OE, Meri S: Ascitic complement system in ovarian cancer. Br J Cancer. 2005, 92 (5): 895-905.
Kaplan EL, Meier P: Nonparametric estimation from incomplete observations. J Am Stat Assoc. 1958, 53: 457-481.
Cox DR: Regression models and life tables. J R Stat Soc B. 1972, 34: 187-220.
Guo C, Liu S, Sun M-Z: Novel insight into the role of GAPDH playing in tumor. Clin Transl Oncol. 2013, 15 (3): 167-172.
Révillion F, Pawlowski V, Hornez L, Peyrat JP: Glyceraldehyde-3-phosphate dehydrogenase gene expression in human breast cancer. Eur J Cancer. 2000, 36 (8): 1038-1042.
Lin P-C, Lin J-K, Yang S-H, Wang H-S, Li A-Y, Chang S-C: Expression of β-F1-ATPase and mitochondrial transcription factor A and the change in mitochondrial DNA content in colorectal cancer: clinical data analysis and evidence from an in vitro study. Int J Colorectal Dis. 2008, 23 (12): 1223-1232.
Huang T-C, Chang H-Y, Hsu C-H, Kuo W-H, Chang K-J, Juan H-F: Targeting Therapy for Breast Carcinoma by ATP Synthase Inhibitor Aurovertin B. J Proteome Res. 2008, 7 (4): 1433-1444.
Suh DH, Kim M-K, No JH, Chung HH, Song YS: Metabolic approaches to overcoming chemoresistance in ovarian cancer. Ann N Y Acad Sci. 2011, 1229 (1): 53-60.
Tamada M, Suematsu M, Saya H: Pyruvate Kinase M2: Multiple Faces for Conferring Benefits on Cancer Cells. Clin Cancer Res. 2012, 18 (20): 5554-5561.
Tristan C, Shahani N, Sedlak TW, Sawa A: The diverse functions of GAPDH: Views from different subcellular compartments. Cell Signal. 2011, 23 (2): 317-323.
Li S-L, Ye F, Cai W-J, Hu H-D, Hu P, Ren H, Zhu F-F, Zhang D-Z: Quantitative proteome analysis of multidrug resistance in human ovarian cancer cell line. J Cell Biochem. 2010, 109 (4): 625-633.
Shin Y-K, Yoo BC, Chang HJ, Jeon E, Hong S-H, Jung M-S, Lim S-J, Park J-G: Down-regulation of Mitochondrial F1F0-ATP Synthase in Human Colon Cancer Cells with Induced 5-Fluorouracil Resistance. Cancer Res. 2005, 65 (8): 3162-3170.
Espinosa I, Catasus L, Canet B, D’Angelo E, Munoz J, Prat J: Gene expression analysis identifies two groups of ovarian high-grade serous carcinomas with different prognosis. Mod Pathol. 2011, 24 (6): 846-854.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6890/13/30/prepub
We would like to thank cytologist Lambert Skoog, MD, PhD, for evaluating imprints, Liss Garberg and Inger Bodin for help with the immunohistochemistry and Lawrence Lundgren for linguistic advice.
This study was funded by The Swedish Cancer Society, The Cancer Research Funds of Radiumhemmet and The Swedish State under the ALF agreement, Stockholm.
The authors declare that they have no competing interests.
EH conceived the study, participated in the PCR and IHC analyses, analyzed data and wrote the manuscript. SE helped design the study and analyzed data. Pathologist JC evaluated IHC. MFS performed the PCR. KS performed surgery and provided tumor samples. HJ did statistical analyses. MS helped in study design and provided conceptual input. EÅL designed the study, analyzed data and revised the manuscript. All authors read and approved the final manuscript.