Research Paper Volume 8, Issue 11 pp 2702—2712

Mutational load of the mitochondrial genome predicts pathological features and biochemical recurrence in prostate cancer

Anton M.F. Kalsbeek1,2, , Eva F.K. Chan1,2, , Judith Grogan3,4,5, , Desiree C. Petersen1,2, , Weerachai Jaratlerdsiri1, , Ruta Gupta3,4,5, , Ruth J. Lyons1, , Anne-Maree Haynes5, , Lisa G. Horvath5,6, , James G. Kench3,4,5, , Phillip D. Stricker7, , Vanessa M. Hayes1,2,4, ,

  • 1 Laboratory for Human Comparative and Prostate Cancer Genomics, Genomics and Epigenetics Division, Garvan Institute of Medical Research, Darlinghurst, NSW 2010, Australia
  • 2 School of Medical Sciences, University of New South Wales, Randwick, NSW 2031, Australia
  • 3 Department of Tissue Pathology and Diagnostic Oncology, Royal Prince Alfred Hospital, Camperdown, NSW 2050, Australia
  • 4 Central Clinical School, Sydney Medical School, University of Sydney, Camperdown, NSW 2050, Australia
  • 5 Cancer Research Division, The Kinghorn Cancer Centre/Garvan Institute of Medical Research, Darlinghurst, NSW 2010, Australia
  • 6 Chris O'Brien Lifehouse, Missenden Road, Camperdown, NSW 2050, Australia
  • 7 Department of Urology, St. Vincent's Hospital, Darlinghurst, NSW 2010, Australia

Received: June 21, 2016       Accepted: September 20, 2016       Published: October 5, 2016      

https://doi.org/10.18632/aging.101044
How to Cite

Abstract

Prostate cancer management is complicated by extreme disease heterogeneity, which is further limited by availability of prognostic biomarkers. Recognition of prostate cancer as a genetic disease has prompted a focus on the nuclear genome for biomarker discovery, with little attention given to the mitochondrial genome. While it is evident that mitochondrial DNA (mtDNA) mutations are acquired during prostate tumorigenesis, no study has evaluated the prognostic value of mtDNA variation. Here we used next-generation sequencing to interrogate the mitochondrial genomes from prostate tissue biopsies and matched blood of 115 men having undergone a radical prostatectomy for which there was a mean of 107 months clinical follow-up. We identified 74 unique prostate cancer specific somatic mtDNA variants in 50 patients, providing significant expansion to the growing catalog of prostate cancer mtDNA mutations. While no single variant or variant cluster showed recurrence across multiple patients, we observe a significant positive correlation between the total burden of acquired mtDNA variation and elevated Gleason Score at diagnosis and biochemical relapse. We add to accumulating evidence that total acquired genomic burden, rather than specific mtDNA mutations, has diagnostic value. This is the first study to demonstrate the prognostic potential of mtDNA mutational burden in prostate cancer.

Introduction

Prostate cancer is the most commonly diagnosed malignancy in men world-wide and a major cause of cancer death [1,2]. As the global population is aging, so is the total health burden of prostate cancer increasing. Clinical management is constraint by significant heterogeneity in disease course from indolence to mortality [3]. Currently, elevated serum prostate-specific antigen (PSA) and/or abnormal digital rectal examination (DRE) are used as the rationale for a needle biopsy and histopathological diagnosis and Gleason scoring, which is a grading system based on prostate gland histological architecture and used for prognosis [4,5]. Although the only biomarker in routine clinical use, PSA testing has led to over-diagnosis and overtreatment with no definitive mortality benefit [6,7]. Ultimately, the effective treatment of prostate cancer requires novel molecular biomarkers that provide early detection of lethal disease associated with metastasis, while appropriate detection of indolent disease is required to prevent overtreatment and, in turn, improve patient quality of life [8].

Genomic changes are significant contributors to prostate cancer from inherited risk [9] to somatic events initiating and driving tumor progression [10]. Although no definitive disease driver mutations have been identified, exome and whole genome sequencing efforts have identified several prostate cancer subgroups based on genome profiling [1113]. However, these studies have focused on the nuclear genome, largely ignoring the maternally inherited mitochondrial genome and its potential as a prostate cancer biomarker.

The significance of the mitochondrial genome can be directly attributed to mitochondrial function, in particular the generation of the majority of cellular energy in the form of adenosine triphosphate (ATP) through oxidative phosphorylation (OXPHOS) of glucose, as well as significant roles in the regulation of apoptosis and calcium homeostasis. Cells have multiple mitochondria, each containing several copies of the mitochondrial genome. The 16,569 base circular haploid genome encodes for 37 genes; 13 protein coding, 2 ribosomal rRNA and 22 transfer tRNAs. A 10-fold higher mutation rate than nuclear DNA, it is not surprising that mutations in the mtDNA have been associated with several cancers, including prostate cancer [1417]. These sporadic variants may become fixed or lost over time, creating a state of heteroplasmy. The assumption is that a variant need to reach a frequency threshold to have an impact on cell function [18,19].

Including 50 variants in the mtDNA variation database MITOMAP, as at 9th of April 2016, a total of 380 prostate cancer associated somatic mtDNA mutations have been reported [16,2023]. Most recently, McCrow et al., 2016 highlighted that prostate cancer presentation is more likely characterized by the total accumulation of mtDNA mutations rather than specific mtDNA mutations [22], while Ju et al., 2014, looking at mtDNA mutations across a range of cancer tissue types, suggested that functionally deleterious mtDNA mutations were more likely to be heteroplasmic as a result of negative selection [20]. These studies provide evidence for the pathogenic impact and diagnostic value of somatic mtDNA mutations, emphasize the need for further deep whole mitochondrial prostate cancer genome sequencing, and importantly highlight the need for associated clinical follow-up data to assess for prognostic value.

In this study, we compared normal-prostate tumor paired mitochondrial genomes from 115 men having undergone a mean follow-up of 107 months post radical prostatectomy. As well as assessing the biological impact of somatic variants on mitochondrial function and tumor formation, we aimed to correlate the mtDNA mutational burden with clinical impact including diagnosis, defined as pathological presentation, and patient outcomes, defined as biochemical and clinical relapse.

Results

Prostate cancer patient clinical and pathological characteristics

All 115 patients presented with prostate cancer and elected for surgical removal of the prostate gland, with a mean of 107 months (range 24-150 months) follow-up. Patient and pathological characteristics may be found in Table 1. In brief, the mean age at radical prostatectomy was 61 years (range 46 to 75 years), with pre-operative PSA levels ranging from 2.43 to18.6 ng/ml. Since our cohort is biased towards more aggressive disease presentation, defined by a Gleason score ≥7 (89/115, 77%), a substantial number of patients (41/115, 36%) had relapsed within three years (range 1 to 126 months) post-surgery. Radical prostatectomy was non-curative in eight patients presenting with advanced stage prostate cancer, two with reported mortality, while eight patients have advanced to metastatic disease. Pathological reassessment of the core biopsy from which DNA was extracted for this study corresponded with the total and primary Gleason score at surgical diagnosis in 46% (53/115) and 71% (82/115) of cases respectively, further detailed in Supplementary Table 1, while tumor purity estimations averaged 57% (range 5% to 85%).

Table 1. Patient and prostate tissue core characteristics.

Gleason 6Gleason 7Gleason 8-10ALL
Patient characteristics
Radical Prostatectomy (RP)
  Total patient numbers264049115
  Mean age in years (range)46 (58-67)63 (49-74)61 (50-75)61 (46-75)
  Mean PSA in ng/mL (range)6.35 (2.43-12)8.36 (2.5-16.4)8.61 (3.5-18.6)8 (2.43-18.6)
Outcomes
  Total follow up in months (range)105 (51-142)122 (63-146)83 (24-150)107 (24-150)
  Nil relapse124272374
  Relapse totals2132641
      Non curative RP0088
      BCR only2111023
      BCR & metastasis20268
      Death from PCa0022
  Months to recorded relapse (range)59 (36-82)54 (8-126)1 (21-78)33 (1-126)
Prostate tissue core characteristics
Pathology
  Total/Major Gleason score matching RP18/2622/2813/2853/82
  Mean est. % tumor purity (range)47 (5-75)60 (20-80)60 (25-85)57 (5-85)
mtDNA mutational spectrum
  Total no SNVs identified6274376
  Patients with SNVs34192750
  Mean SNV count per patient (range)0.23 (0-2)0.68 (0-3)0.88 (0-4)0.66 (1-4)
  Total non-synonymous SNV identified3101326
  Mean absolute SNV allele frequency0.47 (0.17-0.96)0.44 (0.15-0.89)0.49 (0.12-1)0.46 (0.12-1)
  Mean adjusted SNV allele frequency40.65 (0.28-1.37)0.77 (0.24-1.79)0.83 (0.17-1.43)0.75 (0.17-1.79)
  Mean cumulative SNV allele frequency1.35 (0.29-3.25)1.24 (0.28-2.81)1.94 (0.29-3.25)1.42 (0.28-3.25)
Abbreviations: RP, radical prostatectomy; BCR, biochemical recurrence; PCa, prostate cancer; SNVs, single nucleotide variants
1 Two recorded deaths not defined as a result of prostate cancer
2 Bone or visceral metastasis
3 One or more mtDNA SNVs observed
4 Adjusted for pathologically estimated tumor purity

Mitochondrial haplogroup distribution and prostate cancer risk association

Previous studies have suggested that inherited mitochondrial genome variation which defines population-specific haplogroups, may be associated with prostate cancer risk. Specifically, a study of 221 white North American men showed increased risk of haplogroup U for developing PCa [odds ratio 1.95] [24]. Sequencing complete mitochondrial genomes from blood, we found the haplogroup frequencies did not deviate significantly (p = 0.95 by two-sided Fischer’s exact test) from that expected for an Australian European-ancestral population [25] with haplogroup H most commonly represented at 43% and 44% for expected and observed frequencies respectively, followed by U (both 14%), T (both 9%), J (11% versus 9%) and K (8% versus 7%). Haplogroups were also not associated with specimen Gleason score (p = 0.81 by two-sided Fischer’s exact test). Patient specific haplogroups are depicted in Supplementary Table 1. Our study is in agreement with a more recent report that suggests that common European-derived mtDNA haplogroups are not correlated with increased prostate cancer risk [26,27] while the potential association between the earliest diverging L0 African-derived mtDNA haplogroups and aggressive prostate cancer [22] remains to be confirmed.

Spectrum and frequency of acquired prostate cancer specific mitochondrial genome variation

Sequencing mtDNA from matched tumors, we identified a total of 76 somatic single nucleotide variants (SNVs) in 50 patients (43.5%) (Table 1). A single variant recurred in three patients at position mt.16093, resulting in 74 unique variants (see Supplementary Table 2 for full list of variants). While sequencing was performed using the Ion Torrent PGM instrument, we used the Illumina HiSeq X Ten platform to validated 11/12 (91.7%) unique SNVs within 16 randomly selected tumor-normal sample pairs. The single SNV not identified in the HiSeq data represented a loss of a heteroplasmic low allele frequency (0.17) SNV and is most likely to be a discrepancy in variant calling between the platforms. We retained the SNV for our analysis. Compared with publicly available data, 91.9% (68/74) of the observed unique variants are novel to prostate cancer, 86.5% (64/74) novel to any cancer, and 77.0% (57/74) novel to any disease. Although distributed across the entire mitochondrial genome, we found a predominance of variation within the control region and ribosomal RNA genes (Figure 1A). When correcting for length, we found a greater, though not significant, frequency of variation within the non-protein coding than protein coding regions (p = 0.013 by two-sided Fischer’s exact test) (Figure 1B). This proved to be significant for the D-loop and the tRNA regions (p = 0.0044 by two-sided Fischer’s exact test), as previously reported [16]. Of the 41 SNVs mapping to protein coding regions, 26 directly alter an amino acid (non-synonymous mutation) of which 21 have biological potential defined by termination of translation (n=2) or by PolyPhen-2 computational prediction as either possibly (n=2) or probably damaging (n=17). The gene coding for protein Cytochrome b was affected most by damaging variants, both in absolute numbers and when corrected for gene length. Somatic SNVs further showed a predominance of pyrimidine nucleotide transition events, 54% C>T and 24% T>C in line with what was reported previously [20]. To determine the potential of SNVs to impact mitochondrial function (>70% variant frequency), we corrected for tumor purity estimates, resulting variant frequencies ranged 0.17-1.79, (average 0.75). Following this adjusted variant frequency the cumulative variant frequency (CVF) in each patient with variants ranged 0.28-3.25 (average 1.42). Investigating a previously reported ~3.4 kb deletion variant, long range PCR did not show this deletion to be detectable in any of our patients [28].

Number of somatic single nucleotide variants (SNVs) distributed over the gene, control and non-protein coding regions of the mitochondrial genome, defined as (A) absolute number of SNVs and (B) number of SNVs scaled by region length (SNVs per 1000 bases). Genes/regions are represented in order of appearance on the mitochondrial genome, with genes colored by functional prediction of SNV.

Figure 1. Number of somatic single nucleotide variants (SNVs) distributed over the gene, control and non-protein coding regions of the mitochondrial genome, defined as (A) absolute number of SNVs and (B) number of SNVs scaled by region length (SNVs per 1000 bases). Genes/regions are represented in order of appearance on the mitochondrial genome, with genes colored by functional prediction of SNV.

Correlating somatic mtDNA variation with clinical presentation and outcomes

While we found no single somatic mtDNA mutations common between patients, we observed a diagnostic and prognostic correlation with the total somatic variant burden and for the CVF per patient. Specifically, an increase in total number of somatic mtDNA SNVs and CVF was significantly associated with high-risk disease defined as an increased Gleason score at radical prostatectomy p = 0.0132 (Figure 2A) and p = 0.011 (Figure 2B), respectively, and within the sampled biopsy core p = 0.0013, (Figure 2C) and p = 0.0143 (Figure 2D), respectively, with p-values estimated from ANOVA models. Analyses for subgroups of variants showed a significant predictive value for number of non-synonymous variants (p = 0.034 by ANOVA) and synonymous and non-coding variants (p = 0.018 by ANOVA). We point out that, when adjusting for multiple testing using the conservative Bonferroni correction, all six tests drop below the significance threshold at the familywise alpha level of 5% (α = 0.05/7 = 0.0071). However, we note these tests are not independent, with two pairs of highly correlated variables, namely: number of SNV and CVF, and Gleason Scores at prostatectomy and of the biopsy cores.

Correlation of number of somatic single nucleotide variants (SNVs) with Gleason score. (A) Total number of somatic SNVs per tumor and (B) cumulative variant frequency (CVF) categorized by Gleason score at radical prostatectomy. (C) Total number of somatic SNVs per tumor and (D) CVF categorized by Gleason score of screened biopsy core. P-values represents linear model of SNV predicting pathology score. Diamonds and numbers represent the mean of each group.

Figure 2. Correlation of number of somatic single nucleotide variants (SNVs) with Gleason score. (A) Total number of somatic SNVs per tumor and (B) cumulative variant frequency (CVF) categorized by Gleason score at radical prostatectomy. (C) Total number of somatic SNVs per tumor and (D) CVF categorized by Gleason score of screened biopsy core. P-values represents linear model of SNV predicting pathology score. Diamonds and numbers represent the mean of each group.

When investigating the correlation between Gleason score and tumor purity we found a strongly significant correlation (p = 1.6e-5 by ANOVA for the linear model of: Gleason Score ~ Purity), consequently, the number of SNVs also correlates significantly with tumor purity (p = 0.003 by ANOVA for linear model of: SNV number ~ Purity). Additionally, patients with one or more somatic SNVs were significantly more likely than patients without any somatic mtDNA variants to present with a disease relapse (p = 0.012, hazard ratio (HR) 2.17 (confidence interval (ci) 1.17 – 4.04) by log rank test, Figure 3A), with risk for relapse increasing when considering 2 or more SNVs (p < 0.0001, HR 3.82 (ci 1.54 – 9.46) by log rank test (Figure 3B). As Gleason score is currently the most significant predictor for prostate cancer, we determined the power of total mtDNA variation burden to predict disease relapse in our study either alone and in combination with Gleason score, suggesting that number of mtDNA SNVs slightly improves relapse prediction in combination with Gleason score (p = 0.085 by DeLong test for area under curve (AUC)) (Figure 4).

Kaplan Meier curves for disease relapse after surgery, defined by total number of mtDNA single nucleotide variants (SNVs). (A) Patients divided by zero or any mtDNA SNVs. (B) Patients divided by 2 or more and less than 2 mtDNA SNVs. Time in months, with number of cases per mtNDA SNV category represented at time-points 0, 50 and 100 months.

Figure 3. Kaplan Meier curves for disease relapse after surgery, defined by total number of mtDNA single nucleotide variants (SNVs). (A) Patients divided by zero or any mtDNA SNVs. (B) Patients divided by 2 or more and less than 2 mtDNA SNVs. Time in months, with number of cases per mtNDA SNV category represented at time-points 0, 50 and 100 months.

ROC curves for predicting disease relapse. In orange disease relapse predicted by Gleason score of the radical prostatectomy alone (AUC = 0.704), green for number of single nucleotide variants (SNVs) in mtDNA alone (AUC = 0.659) and in blue the Gleason score combined with the number of mtDNA SNVs (AUC = 0.752).

Figure 4. ROC curves for predicting disease relapse. In orange disease relapse predicted by Gleason score of the radical prostatectomy alone (AUC = 0.704), green for number of single nucleotide variants (SNVs) in mtDNA alone (AUC = 0.659) and in blue the Gleason score combined with the number of mtDNA SNVs (AUC = 0.752).

Discussion

In this study, we contribute 68 novel prostate cancer associated mtDNA mutations to the growing somatic catalog, while confirming six others. Taken together, roughly 3% of the total mitochondrial genome has now been shown to have mutagenic potential during prostate tumourigenesis. Accounting for previous study biases towards the restricted analysis of the hypervariable regions (HVRs) within the control region (D-loop region) of mtDNA, the greatest overall impact appears to be within the smaller 1.2 kb non-coding portion of the genome, observing 14 SNVs per kb in our study. We confirm a previously reported bias towards mutational events within the tRNAs, reporting 8.7 SNVs per kb. We predict 32% (24/74) of our identified mutations to have biological potential, with roughly half (9/24) present in sufficient frequency (>70%) to alter mitochondrial function. Further biological studies are required to determine the overall impact of these somatic mutations on energy metabolism, specifically their potential to drive tumor respiration from one of oxidative phosphorylation to aerobic glycolysis [29], a common cancer phenomenon known as the Warburg effect [30].

The unique properties of mtDNA including its small size, high copy number, susceptibility to mutations and heteroplasmic populations has spearheaded investigations into the clinical significance of somatic mtDNA variation. A recent study suggested that, it is more likely the accumulation of mitochondrial mutations across the entire genome rather than specific mutations that contribute to prostate cancer progression [22]. The latter study showed an association between the mutational load, which is the total number and frequency of mtDNA somatic mutations, and higher Gleason score at diagnosis. We concur that the number and cumulative frequency of SNVs in the mitochondrial genome appears to be an independent predictor for prostate cancer presentation, with significant diagnostic potential for aggressive disease. While we confirm a diagnostic value to the number and frequency of somatic mtDNA SNVs based on a single core interrogation, we do not provide support for the previously reported 3.4 kb mtDNA deletion variant. It should be noted, however, that MitomicsTM proposes a minimum five-specimen requirement to achieve adequate sensitivity for early prostate cancer detection.

Historically studies focused on the identification of molecular biomarkers for prostate cancer have lacked prognostic potential. This is the first study, to our knowledge, to investigate the prognostic role of somatic mtDNA in prostate cancer, supported by the availability of extensive long-term clinical follow-up. We demonstrate that an increase in the number of mtDNA somatic SNVs, both the absolute and CVF is significantly correlated to disease relapse and that this prediction appears to improve stand alone Gleason score based risk prediction. The clinical significance of identifying patients prone to relapse has recently been highlighted by the reported benefit from salvage therapy at the earliest opportunity [31]. Associating an increase in number of mtDNA somatic SNVs with both aggressive disease at presentation and poor outcomes, questions whether we observe a limitation to the number of mutagenic events a cell can tolerate. Although we observed a single individual with four SNVs, the upper limit in our study appears to be three mutagenic events. Further cataloging of prostate cancer associated mitochondrial mutations is required not only to determine the spectrum and combination of oncogenic variants, as well as the upper limits for cellular viability, but importantly, to provide significant clinical potential.

Methods

Patient selection and sample preparation

Patients were recruited and consented and biospecimens biobanked according to St Vincent’s Hospital Human Research Ethics approval #HREC/12/SVH/323 at time of surgery. Study inclusion criteria was as follows; (i) a clinicopathological confirmed diagnosis of prostate cancer, (ii) no pre-operative treatment received (including hormone deprivation therapy, chemo or radio-therapy), (iii) a follow up of more than 24 months post radical prostatectomy, and (iv) the availability of a fresh frozen biopsy core with a matching blood sample. Biopsy cores were processed using cryostat-assisted frozen section sample preparation for pathology review (minimally two pathologists JG, RG, and/or JGK) to ensure tumor presence, evaluate the core pathology for percentage adenocarcinoma cells and Gleason score, and maximize for tumor purity during DNA extraction. DNA was extracted from matching blood and prostate tissue using standard methods (Qiagen) and made available under additional site-specific St Vincent’s Hospital Human Research Ethics approval #HREC/15/SVH/227 for whole mitochondrial genome sequencing. Based on these criteria a total of 115 patients were included in this study with clinic-pathological data related to age, PSA level and Gleason score at diagnosis, as well as follow-up with clinical and biochemical relapse data summarized in Table 1 and further elaborated in Supplementary Table 1.

Mitochondrial genome sequencing

To ensure mitochondrial genome specificity and thereby avoiding amplification of nuclear copies of mitochondrial-derived DNA (NuMTs), we performed a dual-amplicon long-range pre-amplification of the mitochondrial genome prior to sequencing on the Ion Torrent PGM instrument, as previously described [22]. In brief, after quantification, the 7.2 kilobase (kb) and 9.7 kb amplified products were combined in a 7:13 ratio, 200 bp libraries prepared using the Ion Xpress Plus Fragment Kit and Ion Xpress™ Barcode Adaptors (ThermoFisher), and six to eight matched normal-tumor pairs pooled and run on a 318v2 Chip. We generated per patient an average of 80 Megabytes (MB) of data, with an average coverage depth of more than 2000 X per mitochondrial genome.

Variant calling and annotation

Sequencing data were analyzed using the Ion Torrent suite v5.0.2.1. Specifically, sequencing reads were quality trimmed, aligned to the revised Cambridge Reference Sequence (rCRS, accession number: NC_012920) [32] and variants called using the Torrent Variant Caller v5.0.2.1 with customized parameters optimized for the mitochondrial genome. Somatic variants were identified by patient matched tumor-normal (blood) mitochondrial genome comparison, using an allele frequency difference greater than 10%. Variant annotation was performed using ANNOVAR [33] and functional impact was predicted using PolyPhen-2 [34]. In addition to single nucleotide variants (SNVs) and indels (small insertions and deletions), we used a previously published long-range 3.8 kb targeted amplification method [22] to screen for the 3379 bp mtDNA deletion suggested to be associated with prostate cancer in a north American cohort [28].

Variant validation

A total of 16 normal-tumor pairs (14% of the study) were randomly selected to undergo whole genome paired-end 2 x 150 bp sequencing using the Illumina HiSeq X Ten platform. A mean coverage of 44X and 71X was generated for the blood and tumor, respectively, resulting in mitochondrial genome coverage of roughly 3000 X and 14000 X, respectively. Mitochondrial sequences were trimmed and filtered and aligned to rCRS using bwa-mem v0.7.12 [35]. SNVs and indels were called using Strelka [36] and manually compared to sample matched SNVs/indels using the Ion Torrent PGM. Although several somatic indels were called using the Ion Torrent data, none were validated and no additional indels observed using the HiSeq X ten data. Conversely, 11/12 Ion Torrent called SNVs were HiSeq X ten validated and therefore only SNVs were used in subsequent mtDNA somatic mutation analyses.

Clinical correlation and statistical analysis

Clinically relevant data was extracted using CanSto ®, an in-house clinical database manager. Presentation at diagnosis was defined by three Gleason score categories (< 7, equal to 7, and > 7). Disease relapse was defined in this study as; non-curative surgery or a rise in PSA >0.2 ng/mL in two consecutive measurements (biochemical recurrence) either with or without the presence of distal visceral or bone metastasis, and/or associated death. Fisher’s exact test was used to test for difference in expected versus observed haplogroups, for difference in presentation at diagnosis and for testing expected versus observed variants in mtDNA regions. ANOVA linear regression models were used to test for clinical correlation, either Gleason score at radical prostatectomy or for the specific biopsy core, with mtDNA somatic mutation rate defined as SNV number or CVF. When more than two statistical tests were performed on subgroups, Bonferroni multiple comparison correction was performed using an alpha of α = 0.05. Statistical analysis and figures were generated in R (R Foundation, Vienna, Austria) using the ggplot2 package [37], with receiver operating characteristic (ROC) curves and area under the curve (AUC) displayed and analyzed using pROC [38]. Prism 7.0 was used to draw and test survival differences using Kaplan Meier Curves.

Acknowledgements

The authors thank patients who donated material for this study and the efforts of the St Vincent prostate cancer biobank for collecting and storing these samples and updating clinical data, specifically Sarita Tiwari and Ben Thompson for sample collection, Quoc Nguyen for management and overseeing ethical requirements.

Funding

This work was supported by the Australian Prostate Cancer Research Centre New South Wales (APCRC-NSW), the Petre Foundation and Sydney Medical School Foundation, the Movember Australia and the Prostate Cancer Foundation Australia (PCFA) Prostate Cancer Bone Metastasis (ProMis) Movember Revolutionary Team Award (MRTA) and a Perpetual IMPACT fund.

Conflicts of Interest

The authors have no conflict of interests to declare.

References

  • 1. Ferlay J, Parkin DM, Steliarova-Foucher E. Estimates of cancer incidence and mortality in Europe in 2008. Eur J Cancer. 2010; 46:765–81. https://doi.org/10.1016/j.ejca.2009.12.014 [PubMed]
  • 2. Canberra; AIHW; Australian Cancer Incidence and Mortality (ACIM) books. Prostate Cancer. 2016.
  • 3. Center MM, Jemal A, Lortet-Tieulent J, Ward E, Ferlay J, Brawley O, Bray F. International variation in prostate cancer incidence and mortality rates. Eur Urol. 2012; 61:1079–92. https://doi.org/10.1016/j.eururo.2012.02.054 [PubMed]
  • 4. Hayes JH, Barry MJ. Screening for prostate cancer with the prostate-specific antigen test: a review of current evidence. JAMA. 2014; 311:1143–49. https://doi.org/10.1001/jama.2014.2085 [PubMed]
  • 5. Miller DC, Hafez KS, Stewart A, Montie JE, Wei JT. Prostate carcinoma presentation, diagnosis, and staging: an update form the National Cancer Data Base. Cancer. 2003; 98:1169–78. https://doi.org/10.1002/cncr.11635 [PubMed]
  • 6. Auvinen A, Moss SM, Tammela TL, Taari K, Roobol MJ, Schröder FH, Bangma CH, Carlsson S, Aus G, Zappa M, Puliti D, Denis LJ, Nelen V, et al. Absolute Effect of Prostate Cancer Screening: Balance of Benefits and Harms by Center within the European Randomized Study of Prostate Cancer Screening. Clin Cancer Res. 2016; 22:243–49. https://doi.org/10.1158/1078-0432.CCR-15-0941 [PubMed]
  • 7. Heijnsdijk EA, der Kinderen A, Wever EM, Draisma G, Roobol MJ, de Koning HJ. Overdetection, overtreatment and costs in prostate-specific antigen screening for prostate cancer. Br J Cancer. 2009; 101:1833–38. https://doi.org/10.1038/sj.bjc.6605422 [PubMed]
  • 8. Heijnsdijk EA, Wever EM, Auvinen A, Hugosson J, Ciatto S, Nelen V, Kwiatkowski M, Villers A, Páez A, Moss SM, Zappa M, Tammela TL, Mäkinen T, et al. Quality-of-life effects of prostate-specific antigen screening. N Engl J Med. 2012; 367:595–605. https://doi.org/10.1056/NEJMoa1201637 [PubMed]
  • 9. Lichtenstein P, Holm NV, Verkasalo PK, Iliadou A, Kaprio J, Koskenvuo M, Pukkala E, Skytthe A, Hemminki K. Environmental and heritable factors in the causation of cancer--analyses of cohorts of twins from Sweden, Denmark, and Finland. N Engl J Med. 2000; 343:78–85. https://doi.org/10.1056/NEJM200007133430201 [PubMed]
  • 10. Baca SC, Prandi D, Lawrence MS, Mosquera JM, Romanel A, Drier Y, Park K, Kitabayashi N, MacDonald TY, Ghandi M, Van Allen E, Kryukov GV, Sboner A, et al. Punctuated evolution of prostate cancer genomes. Cell. 2013; 153:666–77. https://doi.org/10.1016/j.cell.2013.03.021 [PubMed]
  • 11. Lindberg J, Mills IG, Klevebring D, Liu W, Neiman M, Xu J, Wikström P, Wiklund P, Wiklund F, Egevad L, Grönberg H. The mitochondrial and autosomal mutation landscapes of prostate cancer. Eur Urol. 2013; 63:702–08. https://doi.org/10.1016/j.eururo.2012.11.053 [PubMed]
  • 12. Fraser M, Berlin A, Bristow RG, van der Kwast T. Genomic, pathological, and clinical heterogeneity as drivers of personalized medicine in prostate cancer. Urol Oncol. 2015; 33:85–94. https://doi.org/10.1016/j.urolonc.2013.10.020 [PubMed]
  • 13. Cancer Genome Atlas Research Network. The Molecular Taxonomy of Primary Prostate Cancer. Cell. 2015; 163:1011–25. https://doi.org/10.1016/j.cell.2015.10.025 [PubMed]
  • 14. Wallace DC. Mitochondria and cancer. Nat Rev Cancer. 2012; 12:685–98. https://doi.org/10.1038/nrc3365 [PubMed]
  • 15. Chen JZ, Gokden N, Greene GF, Mukunyadzi P, Kadlubar FF. Extensive somatic mitochondrial mutations in primary prostate cancer using laser capture microdissection. Cancer Res. 2002; 62:6470–74. [PubMed]
  • 16. Kloss-Brandstätter A, Schäfer G, Erhart G, Hüttenhofer A, Coassin S, Seifarth C, Summerer M, Bektic J, Klocker H, Kronenberg F. Somatic mutations throughout the entire mitochondrial genome are associated with elevated PSA levels in prostate cancer patients. Am J Hum Genet. 2010; 87:802–12. https://doi.org/10.1016/j.ajhg.2010.11.001 [PubMed]
  • 17. Mizumachi T, Muskhelishvili L, Naito A, Furusawa J, Fan CY, Siegel ER, Kadlubar FF, Kumar U, Higuchi M. Increased distributional variance of mitochondrial DNA content associated with prostate cancer cells as compared with normal prostate cells. Prostate. 2008; 68:408–17. https://doi.org/10.1002/pros.20697 [PubMed]
  • 18. King MP, Attardi G. Human cells lacking mtDNA: repopulation with exogenous mitochondria by complementation. Science. 1989; 246:500–03. https://doi.org/10.1126/science.2814477 [PubMed]
  • 19. Boulet L, Karpati G, Shoubridge EA. Distribution and threshold expression of the tRNA(Lys) mutation in skeletal muscle of patients with myoclonic epilepsy and ragged-red fibers (MERRF). Am J Hum Genet. 1992; 51:1187–200. [PubMed]
  • 20. Ju YS, Alexandrov LB, Gerstung M, Martincorena I, Nik-Zainal S, Ramakrishna M, Davies HR, Papaemmanuil E, Gundem G, Shlien A, Bolli N, Behjati S, Tarpey PS, et al, and ICGC Breast Cancer Group, and ICGC Chronic Myeloid Disorders Group, and ICGC Prostate Cancer Group. Origins and functional consequences of somatic mitochondrial DNA mutations in human cancer. eLife. 2014; 3:3. https://doi.org/10.7554/eLife.02935 [PubMed]
  • 21. Kogelnik AM, Lott MT, Brown MD, Navathe SB, Wallace DC. MITOMAP: a human mitochondrial genome database. Nucleic Acids Res. 1996; 24:177–79. https://doi.org/10.1093/nar/24.1.177 [PubMed]
  • 22. McCrow JP, Petersen DC, Louw M, Chan EK, Harmeyer K, Vecchiarelli S, Lyons RJ, Bornman MS, Hayes VM. Spectrum of mitochondrial genomic variation and associated clinical presentation of prostate cancer in South African men. Prostate. 2016; 76:349–58. https://doi.org/10.1002/pros.23126 [PubMed]
  • 23. Gómez-Zaera M, Abril J, González L, Aguiló F, Condom E, Nadal M, Nunes V. Identification of somatic and germline mitochondrial DNA sequence variants in prostate cancer patients. Mutat Res. 2006; 595:42–51. https://doi.org/10.1016/j.mrfmmm.2005.10.012 [PubMed]
  • 24. Booker LM, Habermacher GM, Jessie BC, Sun QC, Baumann AK, Amin M, Lim SD, Fernandez-Golarz C, Lyles RH, Brown MD, Marshall FF, Petros JA. North American white mitochondrial haplogroups in prostate and renal cancer. J Urol. 2006; 175:468–72. https://doi.org/10.1016/S0022-5347(05)00163-1 [PubMed]
  • 25. Manwaring N, Jones MM, Wang JJ, Rochtchina E, Mitchell P, Sue CM. Prevalence of mitochondrial DNA haplogroups in an Australian population. Intern Med J. 2006; 36:530–33. https://doi.org/10.1111/j.1445-5994.2006.01118.x [PubMed]
  • 26. Giorgi EE, Li Y, Caberto CP, Beckman KB, Lum-Jones A, Haiman CA, Le Marchand L, Stram DO, Saxena R, Cheng I. No association between the mitochondrial genome and prostate cancer risk: The Multiethnic Cohort. Cancer Epidemiol Biomarkers Prev. 2016; 25:1001–03. https://doi.org/10.1158/1055-9965.EPI-16-0111 [PubMed]
  • 27. Fachal L, Gómez-Caamaño A, Alvarez Iglesias V, Gómez Carballa A, Calvo P, Salas A, Vega A. No association between typical European mitochondrial variation and prostate cancer risk in a Spanish cohort. J Hum Genet. 2014; 59:411–14. https://doi.org/10.1038/jhg.2014.46 [PubMed]
  • 28. Maki J, Robinson K, Reguly B, Alexander J, Wittock R, Aguirre A, Diamandis EP, Escott N, Skehan A, Prowse O, Thayer RE, Froberg MK, Wilson MJ, et al. Mitochondrial genome deletion aids in the identification of false- and true-negative prostate needle core biopsy specimens. Am J Clin Pathol. 2008; 129:57–66. https://doi.org/10.1309/UJJTH4HFEPWAQ78Q [PubMed]
  • 29. Ngo DC, Ververis K, Tortorella SM, Karagiannis TC. Introduction to the molecular basis of cancer metabolism and the Warburg effect. Mol Biol Rep. 2015; 42:819–23. https://doi.org/10.1007/s11033-015-3857-y [PubMed]
  • 30. Kim A. Mitochondrial DNA somatic mutation in cancer. Toxicol Res. 2014; 30:235–42. https://doi.org/10.5487/TR.2014.30.4.235 [PubMed]
  • 31. Fossati N, Karnes RJ, Cozzarini C, Fiorino C, Gandaglia G, Joniau S, Boorjian SA, Goldner G, Hinkelbein W, Haustermans K, Tombal B, Shariat S, Karakiewicz PI, et al. Assessing the Optimal Timing for Early Salvage Radiation Therapy in Patients with Prostate-specific Antigen Rise After Radical Prostatectomy. Eur Urol. 2016; 69:728–33. https://doi.org/10.1016/j.eururo.2015.10.009 [PubMed]
  • 32. Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet. 1999; 23:147. https://doi.org/10.1038/13779 [PubMed]
  • 33. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010; 38:e164. https://doi.org/10.1093/nar/gkq603 [PubMed]
  • 34. 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:248–49. https://doi.org/10.1038/nmeth0410-248 [PubMed]
  • 35. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25:1754–60. https://doi.org/10.1093/bioinformatics/btp324 [PubMed]
  • 36. Saunders CT, Wong WS, Swamy S, Becq J, Murray LJ, Cheetham RK. Strelka: accurate somatic small-variant calling from sequenced tumor-normal sample pairs. Bioinformatics. 2012; 28:1811–17. https://doi.org/10.1093/bioinformatics/bts271 [PubMed]
  • 37. Wickham H. ggplot2. New York, NY: Springer New York, 2009.
  • 38. Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez J-C, Müller M. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011; 12:77. https://doi.org/10.1186/1471-2105-12-77 [PubMed]