Skip to main content

Polymorphisms in the mitochondrial oxidative phosphorylation chain genes as prognostic markers for colorectal cancer

Abstract

Background

Currently, the TNM classification of malignant tumours based on clinicopathological staging remains the standard for colorectal cancer (CRC) prognostication. Recently, we identified the mitochondrial oxidative phosphorylation chain as a consistently overrepresented category in the published gene expression profiling (GEP) studies on CRC prognosis.

Methods

We evaluated associations of putative regulatory single nucleotide polymorphisms (SNPs) in genes from the oxidative phosphorylation chain with survival and disease prognosis in 613 CRC patients from Northern Germany (PopGen cohort).

Results

Two SNPs in the 3′ untranslated region of UQCRB (complex III), rs7836698 and rs10504961, were associated with overall survival (HR = 0.52, 95% CI 0.32–0.85 and HR = 0.64, 95% CI 0.42–0.99, for TT carriers). These associations were restricted to the group of patients with cancer located in the colon (HR = 0.42, 95% CI 0.22–0.82 and HR = 0.46, 95% CI 0.25–0.83). Multivariate analysis indicated that both markers might act as independent prognostic markers. Additionally, the TT carriers were ~2 times more likely to develop tumours in the colon than in the rectum. Two SNPs in COX6B1 (complex IV) were associated with lymph node metastasis in a dominant model (rs6510502, OR = 1.75, 95% CI 1.20–2.57; rs10420252, OR = 1.68, 95% CI 1.11–2.53); rs6510502 was associated also with distant metastasis (OR = 1.67, 95% CI 1.09–2.56 in a dominant model).

Conclusions

This is the first report suggesting that markers in genes from the mitochondrial oxidative chain might be prognostic factors for CRC. Additional studies replicating the presented findings are needed.

Peer Review reports

Background

Colorectal cancer (CRC) is the third most common cancer and the fourth-leading cause of cancer death worldwide [1]. Currently, the TNM classification of malignant tumours based on clinicopathological staging remains the standard for CRC prognostication [2]. Many gene expression profiling (GEP) studies on prognosis of CRC have been performed in the last decade using microarray technology, with the aim to identify a gene expression profile to discriminate more aggressive from less aggressive CRC. However, due to their low overlap in the identified genes, no reliable signature useful in the clinical practice has been found. Recently, we published a systematic pathway-based enrichment analysis of 23 independent GEP studies on prognosis of CRC [3]. This analysis indicated the mitochondrial oxidative phosphorylation (OXPHOS) chain as a significantly and consistently overrepresented prognostic category for CRC.

Already several decades ago, it was suggested that impaired oxidative metabolism may cause malignant growth [4]. In cancer cells, there is an enhanced glucose use, slowing the rate of the tricarboxylic acid cycle and oxidative phosphorylation and increasing glycolysis, as the way to generate energy in form of adenosine triphosphate (ATP), despite aerobic conditions [5]. This assumption, known as Warburg’s hypothesis, has been rediscovered by a broad range of experimental approaches showing interaction of mitochondrial metabolism and tumour growth [68]. Additionally, germline mutations in the mitochondrial succinate dehydrogenase (complex II of the OXPHOS chain) subunits SDHD, SDHC, and SDHB are a frequent cause of paragangliomas of the head and neck and of phaeochromocytomas [9].

Thus, the OXPHOS chain appears as a promising candidate pathway for CRC progression. The OXPHOS system, whose proteins are encoded by both nuclear and mitochondrial DNA (mtDNA), consists of five major protein complexes named complex I, II, III, IV, and V, localized on the inner mitochondrial membrane. Its main function is the generation of ATP [10]. We therefore selected putatively functional single nucleotide polymorphisms (SNPs) in the genes from the pathway, and investigated their possible role in the progression of CRC and the prognosis of patients in a group of 613 German CRC patients.

Methods

Study population

The 613 patients included in the study belonged to the population-based PopGen project in Schleswig-Holstein (Germany) [11]. All individuals were of German origin, i.e. with both parents born in Germany. Patients fulfilling either the clinical Amsterdam or Bethesda criteria for hereditary non-polyposis colorectal cancer (HNPCC) were excluded from the study, as were patients with a history of malignant disease or inflammatory bowel disease. For a full description of the characteristics of the 613 patients see Table 1 and Castro et al. [12]. The work has been approved by the corresponding ethical committees and subjects gave informed consent to the work.

Table 1 Characteristics of the 613 CRC patients at the time of diagnosis

Gene and SNP selection

Candidate genes in the OXPHOS pathway (map00190 in the KEGG database [13]) were selected based on a recently published meta-analysis of GEP studies on CRC prognosis [3]. Three out of the seven genes reported in two GEP studies had SNPs in the putative regulatory regions (ATP5C1 from complex V, COX5B and COX6B1 from complex IV). Two additional genes with SNPs in putative regulatory regions, each reported in one GEP study, were also included, namely NDUFS2 from complex I and UQCRB from complex III. With these selection criteria, the four complexes from the OXPHOS chain which allow flowing of protons between the mitochondrial matrix and the intermembrane space, resulting into energy generation in form of ATP, were included in the study. Finally, SNPs in the genes GAPDH and HSPD1 were also selected, since both genes might define the cellular bioenergetic activity of the mitochondria [6].

SNPs within the set of candidate genes (Table 2 and Additional file 1: Table S1) were selected according to the following criteria: minor allele frequency (MAF) > 0.05 in HapMap CEPH population (Utah residents with ancestry from northern and western Europe - CEU), and location in potentially regulatory gene regions (5′ and 3′-UTR, promoter or non-synonymous coding SNPs). In regions with several SNPs in high linkage disequilibrium (LD), tagSNPs were selected for genotyping. LD between the SNPs was calculated using Haploview [14]. The promoter SNP rs2802460 in ATP5C1 was reported to regulate expression of the gene in lymphoblastoid cell lines (eSNPs) on a genome-wide association study of global gene expression [15]. Predictions of the possible impact of amino acid substitutions on the structure and function of the protein were performed with PolyPhen2 [16].

Table 2 SNPs evaluated in the study

Genotyping

The 613 DNAs used in the study were extracted from peripheral venous blood using the FlexiGene chemistry (Qiagen, Hilden, Germany) according to the manufacturer’s protocols. From them, 272 (44.4%) samples were subjected to genome-wide amplification using the illustra GenomiPhi V2 DNA Amplification Kit (GE Healthcare, Freiburg, Germany) according to manufacturer’s instructions. Genotyping was performed using 5 ng of either genomic or whole-genome amplified DNA in 384-well plate format, using KASPar assays on demand (KBiosciences, Hoddesdon, UK) and following the manufacturer’s protocols. Average genotyping call rate was 98.8%. All polymorphisms were in Hardy-Weinberg equilibrium.

Statistical analysis

SAS software version 9.2 (SAS Institute) was used in all analysis. Unadjusted association between the genotyped markers and patient characteristics were evaluated by chi-square test. Odds ratios (ORs) with 95% confidence intervals (CIs) were estimated using logistic regression. Effect of the different genotypes on CRC survival was evaluated using Kaplan-Meier method and was compared using log-rank testing. Analysis of different parameters for prognostic significance was done by univariate and multivariate Cox proportional hazard models. P values < 0.05 were considered statistically significant. Follow-up time was calculated from the date of CRC diagnosis to the CRC specific death or death by any cause or to the end of follow-up (date of last contact with the treating physician). Taking into account the possible biological effects of the selected genes/variants on cancer progression and a possible correlation between the variants, correction for multiple comparisons was excluded.

Results

SNP selection

A total of 29 SNPs were identified in the regions of interest (5′- and 3′-UTR, promoter and non-synonymous coding SNPs) of the seven genes NDUFS2, UQCRB, COX5B, COX6B1, ATP5C1, GAPDH and HSPD1 ( Additional file 1: Table S1). Out of them, 18 SNPs were selected for the association studies, covering all gene regions with a putative regulatory function.

Association with survival

The two polymorphisms in the 3′-UTR region of UQCRB (complex III) showed association with CRC survival (Table 3 and Additional file 1: Table S2). The marker rs7836698 was statistically significantly associated with a decreased risk of death due to CRC (HR = 0.53, 95% CI 0.31–0.91 for TT), or due to any cause (HR = 0.52, 95% CI 0.32–0.85). For rs10504961, only the association with overall survival was statistically significant (HR = 0.64, 95% CI 0.42–0.99 for TT). Interestingly, for both UQCRB polymorphisms, the TT carriers were more likely to develop tumours in the colon than in the rectum, with ORs of 2.02 (95% CI 1.24–3.28) for rs7836698 (TT vs. CC) and 1.74 (95% CI 1.08–2.78) for rs10504961 (TT vs. CC). When survival was analysed according to the cancer site, both SNPs showed association exclusively in the group of patients with cancer located in the colon (HR for overall survival 0.42, 95% CI 0.22–0.82 for rs7836698 TT carriers and 0.46, 95% CI 0.25–0.83 for rs10504961 TT carriers), but not in the rectum (Table 3). The Kaplan-Meier survival curves representing the survival rates of the patients according to their genotypes, as well as the P values for the log-rank test are presented in Figure 1.

Table 3 Association of two polymorphisms in UQCRB (complex III) with survival
Figure 1
figure 1

Kaplan-Meier estimates of CRC specific survival (months) according to UQCRB genotypes, (a) rs7836698, all CRC patients (n = 606), log rank p value = 0.07, (b) rs7836698, all patients diagnosed with colon cancer (n = 324), log rank p value = 0.05, (c) rs10504961, all CRC patients (n = 606), log rank p value = 0.55, and (d) rs10504961, all patients diagnosed with colon cancer (n = 324), log rank p value = 0.06.

We examined also whether inclusion of other variables associated with CRC survival affected the parameter estimates for the UQCRB polymorphisms. The risk of death for carriers of the TT genotype at rs7836698 remained statistically significantly decreased after adjustment for age at diagnosis and TNM stage. Among patients with tumour located in the colon, the hazard ratios were similar both in the uni- and multivariate analyses for rs7836698 and rs10504961, although the statistical significance was sometimes lost.

Association with CRC progression

Three of the 16 SNPs investigated showed an association with one of the different components that form the UICC stage (T, N and M) in a dominant model. Two SNPs in COX6B1 (complex IV) showed an association with the presence of affected lymph nodes (rs6510502, OR = 1.75, 95% CI 1.20–2.57; rs10420252, OR = 1.68, 95% CI 1.11–2.53). The polymorphism rs6510502 was also associated with distant metastasis (OR = 1.67, 95% CI 1.09–2.56). The marker rs7971637 in GAPDH was associated with tumour stage T3/T4 (OR = 1.60, 95% CI 1.11–2.30).

Discussion

In a homogeneous Northern German population of 613 CRC patients, we observed several associations of polymorphisms in the OXPHOS chain with survival and progression of the disease. Mainly, two SNPs in the 3′-UTR region of UQCRB (complex III), rs7836698 and rs10504961, were associated with survival. The two polymorphisms were in relatively strong linkage disequilibrium (D′ = 0.98, r2 = 0.65 in the 613 samples), explaining the association observed with the two SNPs. Added to that, the two studied polymorphisms in COX6B1 (complex IV), rs6510502 in the promoter and rs10420252 in the 5′-UTR, were associated with lymph node metastasis, and rs6510502 also with distant metastasis. The two markers were in high LD (D′ = 0.97, r2 = 0.71 in the 613 samples).

Patients carrying the TT genotype for any of the two SNPs in the UQCRB gene had an increased CRC survival. Interestingly, the association appeared to be restricted to the group of patients with tumour located in the colon. Multivariate analyses confirmed that the association was independent of other variables associated with CRC survival, indicating that both markers might act as independent prognostic markers. In the group of patients with tumour located in the colon the hazard ratios remained low, even though the significance sometimes disappeared when other variables were included. Our results may partially explain the previously reported higher familial risk for colon than for rectal cancer, which was shown to decrease gradually from proximal colon to distal colon and to rectum [17]. Also for survival, differences have been observed between colon and rectal cancer, with rectal cancer having somewhat better survival compared to colon cancer [18].

UQCRB is one of the ten nucleus-encoded subunits of the complex III in the mitochondrial respiratory chain, the ubiquinol-cytochrome c oxidoreductase. In a genome-wide analysis of pancreatic cancer with CGH (comparative genomic hybridization) arrays, a DNA copy number variation was identified at the UQCRB locus, and this gene was suggested as a potential diagnostic marker or a therapeutic target [19]. Recently, it has been demonstrated that UQCRB plays a crucial role in the oxygen sensing mechanism that regulates responses to hypoxia. Using the small molecule terpestacin to inhibit UQCRB, tumour angiogenesis was blocked in vivo [20]. COX6B1 is one of the ten nucleus-encoded subunits of the complex IV in the mitochondrial respiratory chain, the cytochrome c oxidase. Mutations in this gene have been linked to severe infantile encephalomyopathy, a phenotype associated with cytochrome c oxidase deficiency [21]. Our results suggest that genetic variation in the regulatory regions of genes in the OXPHOS chain may associate with CRC progression. Whether such variation affects the gene expression pattern, impairs the function of mitochondria and leads to cancer development and progression remains to be elucidated.

One of the polymorphisms we investigated, rs2802460, located on the promoter of ATP5C1 (complex V), was reported to regulate expression of the gene in lymphoblastoid cell lines in a genome-wide association study of global gene expression [15]. However, this SNP did not show any association with CRC survival or disease progression. A recent report has evaluated the role of 376 tagSNPs (one included in our study, rs11576415 in NDUFS2) in 78 nuclear-encoded mitochondrial genes (including ATP5C1, NDUFS2, and UQCR) in prostate cancer risk, without clear evidence of association [22]. Another recent large association study included variants from 90 genes involved in oxidative phosphorylation (including ATP5C1, COX5B, COX6B1, GAPDH, and NDUFS2) and tested their association with epithelial ovarian cancer risk; no indication of association was found [23]. To our knowledge, this is the first association study of variants in the respiratory chain with CRC survival and progression.

Conclusions

Our findings suggest a possible influence of genes in the OXPHOS chain in CRC survival and progression, and support the newly discovered Warburg’s hypothesis. This is the first report of a possible role of genes in the respiratory chain in the prognosis of CRC patients, but further independent studies are clearly needed to elucidate the validity of the results.

References

  1. Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM: Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int J Cancer. 2010, 127 (12): 2893-2917. 10.1002/ijc.25516.

    Article  CAS  PubMed  Google Scholar 

  2. Walther A, Johnstone E, Swanton C, Midgley R, Tomlinson I, Kerr D: Genetic prognostic and predictive markers in colorectal cancer. Nat Rev Cancer. 2009, 9 (7): 489-499. 10.1038/nrc2645.

    Article  CAS  PubMed  Google Scholar 

  3. Lascorz J, Chen B, Hemminki K, Forsti A: Consensus pathways implicated in prognosis of colorectal cancer identified through systematic enrichment analysis of gene expression profiling studies. PLoS One. 2011, 6 (4): e18867-10.1371/journal.pone.0018867.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Warburg O: On respiratory impairment in cancer cells. Science. 1956, 124 (3215): 269-270.

    CAS  PubMed  Google Scholar 

  5. Levine AJ, Puzio-Kuter AM: The control of the metabolic switch in cancers by oncogenes and tumor suppressor genes. Science. 2010, 330 (6009): 1340-1344. 10.1126/science.1193494.

    Article  CAS  PubMed  Google Scholar 

  6. Cuezva JM, Sanchez-Arago M, Sala S, Blanco-Rivero A, Ortega AD: A message emerging from development: the repression of mitochondrial beta-F1-ATPase expression in cancer. J Bioenerg Biomembr. 2007, 39 (3): 259-265. 10.1007/s10863-007-9087-9.

    Article  CAS  PubMed  Google Scholar 

  7. Ristow M: Oxidative metabolism in cancer growth. Curr Opin Clin Nutr Metab Care. 2006, 9 (4): 339-345. 10.1097/01.mco.0000232892.43921.98.

    Article  CAS  PubMed  Google Scholar 

  8. Koppenol WH, Bounds PL, Dang CV: Otto Warburg’s contributions to current concepts of cancer metabolism. Nat Rev Cancer. 2011, 11 (5): 325-337. 10.1038/nrc3038.

    Article  CAS  PubMed  Google Scholar 

  9. Eng C, Kiuru M, Fernandez MJ, Aaltonen LA: A role for mitochondrial enzymes in inherited neoplasia and beyond. Nat Rev Cancer. 2003, 3 (3): 193-202. 10.1038/nrc1013.

    Article  CAS  PubMed  Google Scholar 

  10. Chandra D, Singh KK: Genetic insights into OXPHOS defect and its role in cancer. Biochim Biophys Acta. 2010, 1807 (6): 620-625.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Krawczak M, Nikolaus S, von Eberstein H, Croucher PJ, El Mokhtari NE, Schreiber S: PopGen: population-based recruitment of patients and controls for the analysis of complex genotype-phenotype relationships. Community Genet. 2006, 9 (1): 55-61. 10.1159/000090694.

    PubMed  Google Scholar 

  12. Castro FA, Forsti A, Buch S, Kalthoff H, Krauss C, Bauer M, Egberts J, Schniewind B, Broering DC, Schreiber S, et al: TLR-3 polymorphism is an independent prognostic marker for stage II colorectal cancer. Eur J Cancer. 2011, 47 (8): 1203-1210. 10.1016/j.ejca.2010.12.011.

    Article  CAS  PubMed  Google Scholar 

  13. Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004, 32: D277-D280. 10.1093/nar/gkh063. (Database issue)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.

    Article  CAS  PubMed  Google Scholar 

  15. Dixon AL, Liang L, Moffatt MF, Chen W, Heath S, Wong KC, Taylor J, Burnett E, Gut I, Farrall M, et al: A genome-wide association study of global gene expression. Nat Genet. 2007, 39 (10): 1202-1207. 10.1038/ng2109.

    Article  CAS  PubMed  Google Scholar 

  16. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR: A method and server for predicting damaging missense mutations. Nat Methods. 2010, 7 (4): 248-249. 10.1038/nmeth0410-248.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Hemminki K, Chen B: Familial risks for colorectal cancer show evidence on recessive inheritance. Int J Cancer. 2005, 115 (5): 835-838. 10.1002/ijc.20938.

    Article  CAS  PubMed  Google Scholar 

  18. Storm HH, Engholm G, Hakulinen T, Tryggvadottir L, Klint A, Gislum M, Kejs AM, Bray F: Survival of patients diagnosed with cancer in the Nordic countries up to 1999–2003 followed to the end of 2006. A critical overview of the results. Acta Oncol. 2010, 49 (5): 532-544. 10.3109/02841861003801148.

    Article  PubMed  Google Scholar 

  19. Harada T, Chelala C, Crnogorac-Jurcevic T, Lemoine NR: Genome-wide analysis of pancreatic cancer using microarray-based techniques. Pancreatology. 2009, 9 (1–2): 13-24.

    Article  CAS  PubMed  Google Scholar 

  20. Jung HJ, Shim JS, Lee J, Song YM, Park KC, Choi SH, Kim ND, Yoon JH, Mungai PT, Schumacker PT, et al: Terpestacin inhibits tumor angiogenesis by targeting UQCRB of mitochondrial complex III and suppressing hypoxia-induced reactive oxygen species production and cellular oxygen sensing. J Biol Chem. 2010, 285 (15): 11584-11595. 10.1074/jbc.M109.087809.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Massa V, Fernandez-Vizarra E, Alshahwan S, Bakhsh E, Goffrini P, Ferrero I, Mereghetti P, D’Adamo P, Gasparini P, Zeviani M: Severe infantile encephalomyopathy caused by a mutation in COX6B1, a nucleus-encoded subunit of cytochrome c oxidase. Am J Hum Genet. 2008, 82 (6): 1281-1289. 10.1016/j.ajhg.2008.05.002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Wang L, McDonnell SK, Hebbring SJ, Cunningham JM, St Sauver J, Cerhan JR, Isaya G, Schaid DJ, Thibodeau SN: Polymorphisms in mitochondrial genes and prostate cancer risk. Cancer Epidemiol Biomarkers Prev. 2008, 17 (12): 3558-3566. 10.1158/1055-9965.EPI-08-0434.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Permuth-Wey J, Chen YA, Tsai YY, Chen Z, Qu X, Lancaster JM, Stockwell H, Dagne G, Iversen E, Risch H, et al: Inherited variants in mitochondrial biogenesis genes may influence epithelial ovarian cancer risk. Cancer Epidemiol Biomarkers Prev. 2011, 20 (6): 1131-1145. 10.1158/1055-9965.EPI-10-1224.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Pre-publication history

Download references

Acknowledgement

This work was supported by the German National Genome Research Network (NGFN-Plus) [01GS08181]; the Deutsche Krebshilfe (German Cancer AID) [107318]; and European Union (EU) [HEALTH-F4-2007-200767].

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Asta Försti.

Additional information

Competing interests

The authors declare that they have no competing interests.

Author’s contributions

JL, KH, and AF conceived and designed the study. JL carried out the molecular genetic studies. JL wrote the initial manuscript. MB performed the statistical analyses. KH and AF provided oversight and conceptual guidance to the project. MB, WS, HK, SB, JH, KH, AF and CS contributed to the final manuscript. WS, HK, HA, JB, JE, SB, TB, SS, JH, and CS participated in the clinical survey and sample collection. All authors read and approved the final manuscript.

Electronic supplementary material

12881_2011_945_MOESM1_ESM.pdf

Additional file 1: Table S1. SNPs in the regions of interest (5′- and 3′-UTR, promoter region and non-synonymous coding SNPs) with MAF > 0 in European population. Table S2. Association of the 16 investigated polymorphisms with survival. (PDF 59 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Lascorz, J., Bevier, M., Schönfels, W.V. et al. Polymorphisms in the mitochondrial oxidative phosphorylation chain genes as prognostic markers for colorectal cancer. BMC Med Genet 13, 31 (2012). https://doi.org/10.1186/1471-2350-13-31

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2350-13-31

Keywords