It has become broadly acknowledged that the redox status of cancer cells could be manipulated to achieve cancer-specific killing since high levels of intracellular ROS activate death signal pathways (reviewed in Ref. [27]). Given that mitochondria are a primary intracellular site of ROS production via OXPHOS during ATP generation [26], it is not surprising that mitochondria have attracted increasing research interest as a promising target for anticancer therapy. It has been reported that mtDNA mutations could influence intracellular ROS levels [45, 46]. In particular, since complexes I and III are the main ROS generation sites, variations in their mtDNA genes could have a significant impact on the overall intracellular ROS level (reviewed in Ref. [47]). However, how specific mtDNA variations influence the efficacy of mitochondria-targeting/ROS-stimulating therapy remains unclear. Our previous data have shown that elevated baseline intracellular ROS levels among different cancer cell lines correlate positively with increased resistance towards cisplatin and dequalinium, both being ROS-stimulating drugs [36]. We hypothesised that certain mitochondrial genetic abnormalities, including variations and copy number (mtDNAcn) changes, could influence intracellular ROS levels, and therefore be utilised as biomarkers to predict cancer cells’ response to ROS-stimulating agents. The present study characterised the interrelationship between mitochondrial genetic parameters (variation and mtDNAcn), ROS levels and ROS-mediated drug response with the aim of generating future novel biomarkers that have a significant association with drug response of cancer cells. Such biomarkers have the potential to be integrated into future clinical practice in order to facilitate personalised medicine.
As a first step to investigating the possible link between mitochondrial genetic abnormalities and the cells’ response to ROS-stimulating agents, a range of molecular, biochemical and bioinformatics methods were employed in the present study to analyse ROS levels, cytotoxicity of cisplatin and dequalinium, as well as mtDNA variations and relative copy numbers in 4 cancer and 1 non-cancerous cell lines (two of the cell lines, cancerous PC-3 and non-cancerous PNT-2, originated from the same tissue type, i.e. prostate).
We found that cancer cells generally carry more variations and a higher mtDNAcn in their mitochondrial genome than non-cancerous cells and have higher baseline ROS levels. This is consistent with the observations of others [47–49] and likely reflects the cancer cells’ adaption to oxidative stress via enhanced DNA repair mechanisms, upregulated anti-apoptotic pathways and consequently increased drug resistance/cell survival (reviewed in Ref. [28]).
We also confirmed our previous observation that increasing baseline ROS level correlates positively with increasing resistance to the ROS-stimulating agents in cancer cell lines. We further demonstrated that increased ROS and drug resistance levels also correlate positively with greater mtDNAcn, higher number of non-synonymous complex I and III variations, and higher number of complex I and III functional variations predicted by the 3D structural modelling (i.e. A10398G, T11120C, C12084T, A13681G, G13708A, C13802T, A13966G and T14798C) (Table S4). The positive correlation between the number of predicted complex I/III functional variations and ROS levels appears more significant for mitochondrial ROS than for overall intracellular ROS. Since complexes I and III are the main contributors to mitochondrial ROS generation, it is not surprising that cancer cells carrying larger numbers of functional variations located in these two complexes would produce higher levels of mitochondrial ROS, as shown in this study.
Intriguingly, the positive correlation between the number of predicted complex I/III functional variations and drug IC50 seems more marked for cisplatin than for dequalinium. This suggests that cisplatin, a conventional nDNA-targeting compound, may rely more on elevated mitochondrial ROS to promote drug resistance in cancer cells compared to dequalinium, a mitochondria-targeting compound. Our observation was echoed by a recent study in which cisplatin-resistant lung cancer cells demonstrated increased mitochondrial mass through upregulation of PGC-1α (the predominant mitochondrial biogenesis promoter) and consequently enhanced intracellular ROS production upon cisplatin treatment [50]. Our data support the previously proposed theory that cancer cells may develop resistance towards cisplatin by upregulating mitochondrial ROS production in order to stimulate anti-apoptotic pathways and promote cell survival [50]. We speculate that this mechanism could be triggered by certain nucleus-mitochondria crosstalk in response to cisplatin treatment. On the contrary, unlike cisplatin, dequalinium is a mitochondria-targeting compound and therefore may not initiate such crosstalk, at least not to the same degree, upon entering the cells. This may explain the much lower IC50s of dequalinium observed in the cancer cells compared to the IC50s of cisplatin in the present study. Furthermore, our data and the previously published data support our hypothesis that mitochondria-targeting therapy can be more effective than conventional therapy since the latter may trigger nucleus-mitochondria crosstalk to promote cell survival.
Using 3D structural modelling, we predicted 8 functional OXPHOS variations: A10398G, T11120C, C12084T, A13681G, G13708A, C13802T, and A13966G are likely to destabilise complex I, impacting on corresponding enzyme activity, which in turn may affect ROS levels and drug response. In support of this view, A10398G was reported in various types of cancer including brain, breast and cervical cancer, and strongly associated with elevated levels of ROS [51–54]; T11120C and C13802T was found in the PC3 cells (but not the PNT-2 cells) previously and shown to inhibit OXPHOS and increase ROS production, and consequently promote tumour growth in vivo when those PC-3 cells were injected into nude mice [55]; C12084T and A13966G variations found in the MDA-MB-231 cells were previously reported to be responsible for reduced complex I activity and enhanced metastatic potential via ROS overproduction and the resultant overexpression of Hif-1α [56, 57]. While G13708A in complex I has been reported in various types of cancer such as breast and colorectal cancer [58, 59], A13681G has no previously reported disease associations. There are no previous reports documenting G13708A or A13681G and ROS levels, although our predictions suggest this would be worth investigating. On the other hand, we predict that the final functional OXPHOS mutation: T14798C, despite frequently occurring in the healthy population (like some of the other variations; Table S3), is likely to affect the association/dissociation of ubiquinone at the Qi-site, impacting on the activity of complex III, and this in turn may affect the levels of ROS produced by the complex and consequently drug response. In support of this view, a study conducted by Keatley et al. demonstrated that glioblastoma cell lines that carried T14798C had elevated complex III activity, oxidative stress, and different drug sensitivity levels compared to cells that did not contain the mutation. Further, glioblastoma patients with T14798C had worse prognosis than non-carriers [10]. This is in keeping with an earlier study which demonstrated that the activity of yeast mutant complex III genetically manipulated to contain the equivalent variation of T14798C was significantly impaired and the mutant showed more sensitivity to mitochondrial-targeting compounds [60].
In contrast to increased ROS and drug resistance levels correlating positively with greater mtDNAcn, larger number of non-synonymous complex I and III variations, and larger number of predicted complex I and III functional variations, none of these parameters correlated with the number of D-loop variations. The mtDNAcn result is surprising given that the D-loop region is a major control site for mtDNA replication and transcription, and higher mtDNAcn has been reported in various types of cancer carrying D-loop variations [61, 62]. One possibility is that specific variations located within the light and heavy strand promoters in the D-loop region rather than the total number of D-loop variations are more important. Such variations could affect the binding affinities of the initiators and modulators of mtDNA transcription, thus, disturbing the rate of transcription, RNA primer synthesis as well as mtDNA replication at the H-strand origin of replication (OH) [63]. In other words, quality is more important than quantity in terms of functional variations in the D-loop region. In the present study, only one variation appeared to be located in a key regulatory region. Indeed, 310InsC, present in Ishikawa, MDA-MB-231 and PC-3 cells, is located in CSB2, one of the three conserved sequence blocks within the H-strand origin of replication (OH) that harbours critical functional motifs for the initiation and regulation of H-strand replication. Given that insertion variations have a real fundamental impact on DNA function, as they completely change the sequence of the region and result in the shift of the original functional motifs, it is likely that the 310InsC variation may interfere with DNA-primer interaction and have a negative impact on the replication of the H-strand in those cell lines. This observation might explain why the above cell lines had lower mtDNA copy numbers compared to Caco-2 that did not harbour this insertion.
The apparent relationship between mtDNAcn, ROS levels and drug IC50s could also represent cellular adaptation to the level of OXPHOS activity that in turn might be related to the number of functional OXPHOS variations. One example of such adaptation includes the retrograde signalling pathway, where dysfunctional mitochondria (caused by functional OXPHOS variations) can lead to differences in ROS levels, which in turn lead to the activation of various nuclear responses, consequently promoting multiple pathways that regulate energy homeostasis, oxidative stress, mitophagy, fission, fusion and other functions to facilitate cellular adaption strategies and hence regulate the transcription and translation of genes responsible for mitochondrial biogenesis (e.g. PGC-1α), mtDNA replication and OXPHOS functions [64]. Taking these observations together, we propose that cell line-specific functional OXPHOS variation profiles could mechanistically contribute to the mitochondrial activity, baseline ROS levels, mtDNAcn and response to ROS-stimulating agents observed in cancer (Fig. 8), and the best illustration of this from our study is the comparison between the cancerous PC-3 cells and non-cancerous PNT-2 cells, as they are both derived from prostate.
A recent publication by Cocetta and co-authors provides a sophisticated review on mitochondrial involvement in cisplatin resistance [65]. As pointed out by the authors, so far, existing knowledge of the retrograde signalling is still very limited. Here we present a hypothesis based on our data that specific functional variations in the OXPHOS complexes I and III are responsible for enhanced ROS production and the subsequent retrograde communication, rendering resistance to cisplatin. Our data not only support previous findings that redox imbalance and nucleus-mitochondria crosstalk play important roles in cisplatin resistance, but also provide a potential mechanistic link to specific mtDNA variations. Furthermore, our studies were not restricted to cisplatin only, hence our findings would provide insights into the broader area of ROS-stimulating/mitochondria-targeting therapy as well as the related biomarker studies.
Measuring mitochondrial activity and the expression level of genes controlling the mitochondrial biogenesis and mtDNA replication, such as PGC-1α, DNA polymerase γ (POLγ) and mitochondrial transcriptional factor A (TFAM), in the future would obviously help corroborate this hypothesis.
Now that we know the characteristics of the cell lines in detail, there are a number of approaches that could be employed in the future to further validate our hypothesis that functional OXPHOS variation or a group of variations contribute to the intracellular ROS/drug response phenotypes observed. One way would be to create cybrids which involves the transfer the cancer cell mitochondria/mtDNA containing the functional OXPHOS variations onto a common wild type nuclear background, and the observation whether the cellular phenotype is transferred. However, the process of generating cybrids is not without its problems and can create irreversible epigenetic changes [66], as well as nuclear genome instability and variation [67], potentially confounding the ability to identify one single mtDNA variation (or a group of mtDNA variations) as the contributing factor(s), and so the insights generated in this paper, especially the comparison between PC-3 and PNT-2, are still useful. Another possibility would be to introduce the functional OXPHOS variations into Saccharomyces cerevisiae (Baker’s yeast) cells, as unlike mammalian mtDNA, its mitochondrial genome is amenable to genetic manipulation. This has proved useful for the functional validation of T14798C [60], but the lack of respiratory complex I in S. cerevisiae precludes a similar validation of functional complex I variations identified in this study. Recent technological advances suggest that such precise manipulation of mammalian mtDNA (and thus the functional validation of complex I variations) may be possible in the future, and may in fact provide a better alternative for functionally validating mtDNA variations than both cybrids and the yeast model system [68].
It is worth mentioning that two cell lines of distinctly opposite characteristics were identified among the tested cancer cell lines in our studies that was either the most sensitive (i.e. Ishikawa) or the most resistant (i.e. Caco-2) to ROS-stimulating therapy. The profiles of Ishikawa and Caco-2 based on the ROS/mtDNAcn/functional complexes I/III mutation parameters could be further validated as benchmarks to assess cancer cells’ response towards ROS-stimulating therapy. In addition, these two cell lines would make for interesting candidates for future downstream cybrid/mtDNA gene editing studies to help validate that the differences observed in the present study are due to mtDNA pattern and not due to differences in their nuclear proteomes.