Skip to main content

High leukocyte mitochondrial DNA copy number contributes to poor prognosis in breast cancer patients

Abstract

Background

Compelling evidence has indicated a significant association between leukocyte mitochondrial DNA copy number (mtDNAcn) and prognosis of several malignancies in a cancer-specific manner. However, whether leukocyte mtDNAcn can predict the clinical outcome of breast cancer (BC) patients has not been well investigated.

Methods

The mtDNA copy number of peripheral blood leukocytes from 661 BC patients was measured using a Multiplex AccuCopy™Kit based on a multiplex fluorescence competitive PCR principle. Kaplan–Meier curves and Cox proportional hazards regression model were applied to investigate the association of mtDNAcn with invasive disease-free survival (iDFS), distant disease-free survival (DDFS), breast cancer special survival (BCSS), and overall survival (OS) of patients. The possible mtDNAcn-environment interactions were also evaluated by the Cox proportional hazard regression models.

Results

BC patients with higher leukocyte mtDNA-CN exhibited a significantly worse iDFS than those with lower leukocyte mtDNAcn (5-year iDFS: fully-adjusted model: HR = 1.433[95%CI 1.038–1.978], P = 0.028). Interaction analyses showed that mtDNAcn was significantly associated with hormone receptor status (adjusted p for interaction: 5-year BCSS: 0.028, 5-year OS: 0.022), so further analysis was mainly in the HR subgroup. Multivariate Cox regression analysis demonstrated that mtDNAcn was an independent prognostic factor for both BCSS and OS in HR-positive patients (HR+: 5-year BCSS: adjusted HR (aHR) = 2.340[95% CI 1.163–4.708], P = 0.017 and 5-year OS: aHR = 2.446 [95% CI 1.218–4.913], P = 0.011).

Conclusions

For the first time, our study demonstrated that leukocyte mtDNA copy number might influence the outcome of early-stage breast cancer patients depending on intrinsic tumor subtypes in Chinese women.

Peer Review reports

Introduction

Breast cancer (BC) is the first common malignancy and the fifth leading cause of cancer death among women in China [1]. BC is highly heterogeneous at molecular and clinical levels [2]. Generally, some prognostic factors for early-stage breast cancer (EBC) survival include tumor size, lymph node involvement, tumor grade, and hormone receptor (HR) status. The survival rate has increased due to improved management strategies in the past few decades [3]. However, EBC patients with the same traditional stage often show different clinical outcomes because of tumor heterogeneity.

Consequently, biomarkers are urgently needed to complement the traditional system for a more accurate prognostic prediction of EBC and thus benefit individualized treatment. Some theories indicate that genetic background is an important determinant of metastatic potential. The theories also suggest that it might be possible to define metastasis susceptibility based on gene variation in readily accessible tissues (for example, blood) rather than tumors [4]. The study of mtDNA copy number variation in breast cancer is of particular interest, given the role of oxidative stress in the disease’s etiology.

Mitochondria are essential cellular organelles in eukaryotic cells and perform multiple cellular functions, notably in energy metabolism, reactive oxygen species (ROS) generation, and apoptosis [5]. Mitochondria contain their DNA (mtDNA), which is susceptible to damage caused by high reactive oxygen species due to the lack of protective histones and its diminished DNA repair capacity [6]. There can be several 100 or even more than 1000 copies of mitochondrial genomes in each eukaryotic cell. Although mtDNA copy number is significantly varied among cells from different tissue origins [7], it has been observed that there is a good correlation between the number of mtDNA in various cell types within the same individual [8]. The mtDNA copy number may change significantly under different internal or external microenvironments, damaging the OXPHOS system and the enhanced ROS production. This scenario has been believed to contribute to the initiation and development of tumors [9].

Quantitative changes in mtDNA content have been observed in a variety of malignant tumors, such as BC [10], gastric cancer [11], head and neck cancer [12] and colorectal cancer [13]. Previous studies have also shown that changes in mtDNA content in tumor tissues were related to tumor staging, prognosis, and treatment response, again in a cancer type-specific manner [14, 15]. Yu et al. have found that reduced mitochondrial DNA copy number in tumor tissues was correlated with tumor progression and prognosis in Chinese BC patients [10]. There have been a series of reports on the relationship between mtDNA content in peripheral blood leukocytes and cancer susceptibility [16,17,18,19]. However, only a few studies have shown that higher mtDNA content in peripheral blood leukocytes is associated with poor prognosis in colorectal cancer [20], hepatocellular carcinoma [21], and glioma [22]. In addition, Xia et al. have found that leukocyte mtDNA content is associated with the T stage of BC, indicating that high mtDNA content may facilitate cancer progression, but did not explore the relationship with the prognosis of BC [23]. Thus, the effect of mtDNA copy number in peripheral blood leukocytes on BC patient prognosis has not been explored.

This study measured the leukocyte mtDNA copy number in peripheral blood from EBC patients and analyzed the association between leukocyte mtDNA copy number and the survival of patients in Southeastern China with EBC. To the best of our knowledge, this is the first study to assess mtDNA content’s value in predicting the prognosis of EBC patients.

Materials and methods

Study population

Our hospital-based study recruited 681 BC cases from Fujian Medical University Union Hospital between July 2000 and October 2014. Patient eligibility criteria in this study were as follows: (i) histopathologically confirmed with invasive BC; (ii) subsequently treated with curative surgical resection and systemic therapy; (iii) availability of complete clinical and follow-up data, (iv) no history of other malignancy and (v) alive at least 1 month after surgery. According to this screening criteria, we ruled out 20 patients, containing 6 patients with breast carcinoma in situ and 14 patients with incomplete clinical information. Finally, 661 patients with resected BC were included in the present study for predictive analysis.

Clinicopathological and demographic information was collected from the hospital records, and survival data were obtained from the followed-up database renewed annually. The patients were staged according to the 7th version of the American Joint Commission on Cancer (AJCC) tumor-node-metastasis (TNM) staging system [24]. Estrogen receptor (ER)/progesterone receptor (PR) positivity was determined by IHC analysis of the number of positively stained nuclei (≥ 10%), and HR positivity was defined as being either ER + and/or PR+. Tumors were considered human epidermal growth factor-2 (HER2) positive when cells exhibited strong membrane staining (3+). Expressions of 2 + would require further in situ hybridization testing for HER2 gene amplification, while expressions of 0 or 1 + were regarded as negative. The subtypes were categorized [25]: luminal A (ER+, PR+ >20%, HER2-, Ki67<14% or grade I when Ki67 was unavailable), luminal B (HR+, HER2-, Ki67>14% or grade II/III when Ki67 was unavailable or HR+, HER2+); HER2 enriched (HR-, HER2+) and triple-negative (HR- and HER2-). The Institutional Ethics Committee approved the study, and all participants consented to the testing at the time of their participation and contributed data.

DNA extraction

Blood samples were collected in EDTA anticoagulant tubes and were centrifuged within 30 min. Genomic DNA was extracted from blood samples by using the Whole-Blood DNA Extraction Kit (Bioteke, Beijing, China), according to the manufacturer’s protocol, and then the genomic DNA was aliquoted and stored at -80 °C for future analysis. The purity of the DNA samples was assessed before mtDNA copy number assessment, and samples without sufficient DNA yield were not included in mtDNA copy number analysis.

MtDNA copy number assessment

Relative mtDNA content was measured by a custom-by-design Multiplex AccuCopyâ„¢Kit (Genesky Biotechnologies Inc., Shanghai, China) based on a multiplex fluorescence competitive PCR principle as previously described [26] which can interrogate copy number variation (CNV) status at multiple genomic loci in the same assay reaction.

The methods below briefly describe the manufacturer’s process. A total of 2 target genomic segments within the ND1 gene and 6 reference segments (2, 10, 16, 18, 19, and 20 p) were chosen for the AccuCopy assay. The primers for both target and reference segments were synthesized. The forward primers were fluorescent-labeled at Genesky Biotechnologies (Shanghai, China). The competitive DNAs for the ND1 and 6 reference segments were synthesized in double-strand and provided a mixture from Genesky Biotechnologies (Shanghai, China). These competitive DNAs are almost identical to their homologies in the human reference genome except for 1–2 base pairs deleted. The primers of target segments, reference segments, and probes information were shown in Table S1.

The PCR reaction was prepared in 20 µl for each sample, containing a mixture of 2 µl target genomic DNA (5 ng/µl) with 2 µl reference segment DNA, 1 µl Multiplex PCR Fluorescence Primer Mix (AccuCopy™), 10 µl 2 × PCR Master Mix (Genesky Biotechnologies), and 5 µl ddH2O. 10 µl 2 × PCR Master Mix (Genesky Biotechnologies), and 5 µl ddH2O. The program used was an initial denaturation step of 95 °C for 10 min followed by the first 11 cycles of denaturation at 94 °C for 20 s, annealing at 60 °C for 40 s (the annealing temperature was decreased by 0.5 °C in each consecutive cycle), and elongation at 72 °C for 1.5 min, followed by the second 24 cycles of denaturation at 94 °C for 20 s, annealing at 59 °C for 30 s, and elongation at 72 °C for 1.5 min, and a final extension step at 60 °C for 60 min and then at 4 °C forever. The PCR product was diluted in a 1:5 ratio, and 1 µl diluted products were mixed with 0.5 µl 500 (Liz) size standard and 8.5 µl Hi-Di formamide (both from Applied Biosystems, Foster City, CA). The mixture was subjected to a denaturation step of 95 °C for 5 min and electrophoresed in a 3730XL genetic analyzer (ABI, Carlsbad, CA, USA). Raw data were analyzed by GeneMapper 4.0 (ABI). The height and area data for all specific peaks were exported into a Microsoft Excel file.

The sample/competitive (S/C) peak ratio was calculated for the 2 target and 6 reference segments. The S/C ratio for each target fragment was first normalized based on four reference segments, respectively. The 6 normalized S/C ratios were further normalized to the median value in all samples for each reference segment, respectively, and then averaged. If one of the six normalized S/C ratios deviated >25% from the average of the other five, it was excluded from further analysis. The copy number of each target segment was determined by the average S/C ratio times two, given that the copy numbers of four reference segments are two in the diploid genome. Twenty samples were randomly selected as blinded duplicates for quality assessment purposes, and excellent concordance was obtained.

Statistical analyses

Overall survival (OS) and BC-specific survival (BCSS) were our primary endpoints and defined as the time from cancer diagnosis to mortality for all causes and BC, respectively. Invasive disease-free survival (iDFS) and distant disease-free survival (DDFS) were our secondary endpoints and calculated separately from the date of diagnosis to the date of any recurrence and distant recurrence to follow-up cut-off time [27]. Survival data were analyzed using the Kaplan–Meier method with the log-rank test. Survival analysis was discontinued at 5 years of follow-up because 5-year survival is an important prognostic indicator for BC and the survival data after five years is not yet mature. To adjust for possible confounders and evaluate the form of the relationship, Cox proportional hazards models were fitted using different adjustment parameters: model 1 adjusted for age at diagnosis; model 2 additionally adjusted for hormone receptor status and HER2 status; model 3 was the same as model 2 plus an adjustment for tumor size, Lymph node involvement, and grade. The Cox proportional hazards model results showed the non-linearity in the effects of mitochondrial copy number on cancer progression.

Moreover, nonlinear p-splines with knots were created to evaluate the functional form of the relationship between mtDNA copy number and BC patients’ survival as previously described [28, 29]. The splines’ knots were at the quintiles of the mtDNA copy number distribution. The results of the splines and the estimates for quintiles 2, 3, 4, and 5 compared with quintile 1 revealed that a categorical analysis of the data by combining the 2nd, 3rd, 4th and 5th quintile of mtDNA copy number and set quintile 1 as a reference fits the data better than modeling the data on a continuous scale. The Chi-square test was used to examine differences in categorical variables between subgroups. The student’s t-test was used to analyze the difference of normally distributed continuous variables between two groups. The Cox regression model calculated the hazard ratios (HRs) and 95% confidence interval (CI) for each factor in multivariate analyses. The main analysis started with a full Cox regression model (including copy number and all prognostic factors). A global likelihood ratio test was carried out to test whether all factors together significantly influenced BC survival. In case of a nonsignificant result, no further analyses were carried out, to avoid false-positive results. If, however, the P-value was significant, mtDNA copy number and all prognostic factors will be further analyzed. A stratified multivariate analysis of clinical parameters was performed to evaluate the impact of mtDNA copy number on BC prognosis in specific patients. The possible mtDNAcn-environment interactions were also evaluated by the Cox proportional hazard regression models. All tests were 2-sided, and P-values < 0.05 were considered statistically significant. SAS 9.4 (SAS Institute Inc., Cary, NC) was used for all statistical analyses.

Results

Patient characteristics and clinical features

Participants in the 661 early BC cohort were female, and their mean age was 46.73 ± 10.28 years old at BC diagnosis. The clinical characteristics and survival of the 661 BC patients are summarized in Table 1. During a follow-up time of 5 years, 279 cases experienced recurrence (61 locoregional and 254 distant), and 180 died (176 died of BC and 4 died of another disease). No significant difference in BC-iDFS, BC-DDFS, BCSS, and OS was shown in the subgroup of age at diagnosis 40 (P = 0.136, 0.331, 0.754 and 0.934). Nevertheless, patients with a tumor size > 2 cm, lymph node-positive, grade III, clinical-stage II + III, or HER2 positive had significantly shorter survival times, whereas HR positivity remarkably improved the survival of BC patients (log-rank P < 0.05, Table 1). Furthermore, our intrinsic molecular subtypes (luminal A, luminal B, HER2-enriched, and triple-negative) were also associated with significantly different survival (log-rank P < 0.05, Table 1). Treatment characteristics of patients showed in Table S3.

Table 1 Patients’ clinicopathological characteristics and clinical outcome

Effects of mtDNA copy number on survival of BC

We performed a Cox regression analysis based on quintile categorization of mtDNA copy numbers applying different adjustments for confounders. Since the estimates for quintiles 2, 3, 4, and 5 pointed in the same direction with no clear indication of a linear association, we set the first quintile as the reference category to illustrate the quintiles 2–5 values’ influence. Details for each quintile and adjustment model can be found in Supplementary Table S2. In line with this observation, nonlinear splines showed an increased risk for patients within the 2–5 quintiles of mtDNA copies, compared with quintile 1 for all 4 outcome variables (Figure S1). For reasons of statistical power, we collapsed quintiles 2–5 and compared them to quintile 1, which was a reference in all subsequent analyses. The detailed characteristics comparing combined quintiles 2–5 against quintile 1 showed in Table 2. Patients with higher copy numbers (quintiles 2–5) have higher grades (P = 0.041). The Kaplan-Meier curves and log-rank tests comparing combined quintiles 2–5 against quintile 1 (Fig. 1) showed significant differences concerning 5-year iDFS, DDFS, and OS, except for BCSS (P = 0.018 for iDFS, P = 0.042 for DDFS, P = 0.061 for BCSS, and P = 0.043 for OS).

Table 2 Associations between mtDNA copy number and clinicopathological characteristics in BC patients
Fig. 1
figure 1

Kaplan-Meier curves for iDFS (a), DDFS (b), BCSS (c), and OS (d) due to BC by Quintile 2–5 AND Quintile 1 of mitochondrial DNA copy number. Abbreviations: iDFS, invasive disease-free survival; DDFS, distant disease-free survival; BCSS, BC special survival; OS, overall survival

The global likelihood ratio test for the full Cox regression model (including copy number and all prognostic factors) showed that all factors together significantly affect breast cancer prognosis(P<0.0001). Cox regression analyses model 1 in Table 3 shows that a higher mtDNA copy number (quintile 2–5) was significantly associated with a 1.476-fold higher risk of 5-year BC-specific recurrence (HR = 1.476[95% CI 1.073–2.031], P = 0.016), 1.418-fold higher risk of 5-year distant recurrence (HR = 1.412[95% CI 1.015–1.981], P = 0.040), and 1.522-fold higher risk of 5-year overall death (HR = 1.522[95% CI 1.011–2.291], P = 0.044). However, after adjusting for age at BC diagnosis, tumor size, lymph node involvement, grade, hormone receptor status, and HER2 status, only association with iDFS remained significant. In model 3, quintile 2–5 was significantly associated with shorter iDFS (5-year iDFS: adjusted HR = 1.433[95% CI 1.038–1.978], P = 0.028; Table 3).

Table 3 Results of Cox model on breast cancer risk for the quintile 2 to 5 of mtDNA copy number versus quintile 1 (as reference category)

Stratification and interaction analysis

The associations between mtDNA copy number and BC survival were then evaluated by stratified analysis of age at diagnosis, tumor size, lymph node involvement, grade, hormone-receptor status, and HER2 status, while interaction analysis was also performed between these covariates and mtDNAcn on BC survival. Statistically significant multiplicative interactions on BC survival were found between mtDNA copy number and hormone-receptor status (adjusted p for interaction: 5-year BCSS: 0.028, 5-year OS: 0.022).

Prognostic analysis in HR + patients

As shown in Table 4, we found that higher mtDNA copy numbers (quintile 2–5) were associated with shorter survival of the patients who were positive hormone-receptor status (HR+: 5-year BCSS: aHR = 2.340[95% CI 1.163–4.708], P = 0.017 and 5-year OS: aHR = 2.446[95% CI 1.218–4.913], P = 0.011). However, compared with HR-positive patients, the performance of mitochondrial copy number in HR-negative patients is the opposite and has no significance (Table 4). We further conducted a stepwise Cox proportional hazard analysis to evaluate the effects of clinical variables and mtDNA copy number on BCSS and OS in HR-positive patients. Three variables (node status, grade, mtDNAcn) were selected for the final regression model of BC-specific survival, while four variables (node status, mtDNAcn, grade, HER2) were selected for the final regression model of overall survival (Table 5). In the HR-positive subgroup final model, the mtDNAcn was an independent risk factor for 5-year BCSS with a 2.240-fold (95% CI, 1.204–4.863) increased risk for 5-year OS with a 2.464-fold (95% CI, 1.227–4.948) increased risk.

Table 4 Hazard ratio of tumor survival for all 661 patients with breast cancer according to the mtDNA copy number in different subgroups stratified by clinical parameters
Table 5 Results of stepwise Cox regression analysis on breast cancer-specific survival and overall survival in HR-positive breast cancer patients

Discussion

This study evaluated the possible relation between mtDNA copy number in peripheral blood and BC clinical outcome in a cohort of 661 EBC cases. We found that patients with high mtDNA copy numbers showed significantly shorter iDFS than those with low mtDNA (5-year iDFS: adjusted HR = 1.433[95% CI 1.038–1.978], P = 0.028). The most significant and novel result of this study is that, in the HR + subgroup, patients with high copy number have a worse prognosis than those with low mtDNA (HR+: 5-year BCSS: aHR = 2.340[95% CI 1.163–4.708], P = 0.017 and 5-year OS: aHR = 2.446[95% CI 1.218–4.913], P = 0.011). In short, the high mtDNA copy number in peripheral blood led to a worse prognosis for EBC patients, especially in the HR-positive subgroup patients.

Previous studies have proposed that high leukocyte mtDNA copy number in patients may impair the immune defense against tumor cells and lead to poor prognosis. MtDNA content is regulated by a wide range of environmental and genetic factors. For example, raised mtDNA content has altered oxidative stress, aging, immune response activation, and environmental exposures [30,31,32]. Several studies have demonstrated that high leukocyte mtDNA content is associated with a worse prognosis for patients with glioma [22], HCC [21], and colorectal cancer [20]. These studies have reported that cancer patients with high leukocyte mtDNAcn may be in an immunosuppressive state. Furthermore, mitochondrial ROS, generally increased with mtDNA copy number, can inhibit the expansion and functions of helper T cells and killer cells in peripheral blood [33, 34]. In line with these findings, our study demonstrated that a higher leukocyte mtDNA copy number was associated with poor prognosis in HR-positive EBC patients.

However, several reports showed that lower mtDNAcn in breast tumor tissues was associated with a worse prognosis [10, 35, 36]. Only one previous study evaluated the progression role of mtDNA copy number in peripheral blood leucocytes in BC patients. Xia et al. reported that leukocyte mtDNA content is associated with the T stage of BC, indicating that high mtDNA content may facilitate cancer progression [23]. Nevertheless, the relationship between tumor size and prognosis is unstable, and Xia et al. did not have any follow-up data. We are the first and only study to explore the relationship between the mitochondrial DNA copy number of peripheral blood cells and the prognosis of BC patients. We mainly found that higher leukocyte mtDNA content was associated with poor BCSS and OS in HR-positive EBC patients.

In contrast to mtDNAcn reduction in breast tumor tissues, the studies of leukocyte mtDNAcn of BC patients and its association with prognosis have been heterogeneous. Rai et al. showed an inverse relation in blood vs. tissue samples in BC patients [37]. Our results were consistent with a meta-analysis that summarized the results of 18 studies involving 3961 cancer patients and concluded that elevated mtDNAcns in peripheral blood leukocytes and tumor tissues predict the opposite outcome of cancer [38]. Biologically, mtDNAcn is specific to different tissue types, reflecting differing energy requirements [39, 40]. Collectively, these findings support the point of view that changes in mtDNA content are involved in the progression of BC.

High mtDNA content was significantly associated with poor prognosis, mainly in the HR-positive subgroup but not in the HR-negative subgroup, indicating host characteristics’ modulating effects through unknown mechanisms. HR positivity was either ER + and/or PR+. After binding to estrogen or progesterone, HR can interact with the mtDNA to maintain mitochondrial function, inhibit cell apoptosis, and overcome oxidative damage [41, 42]. The different performance of mtDNA copy numbers in the HR + and HR-subgroups may be due to this interaction. Gao et al. found that estrogen indirectly regulated mitochondrial DNA replication and repair functions [43, 44]. Previous studies have identified that the instability of mtDNA due to exposure to estrogen induces cancer cell metastasis [45]. As a result, we proposed that in HR-positive patients, the oxidative damage in the tumor caused by estrogen or progesterone and the immunosuppression of the high copy number in the peripheral blood cells together lead to a significantly worse prognosis.

As the first exploration of a possible relationship between leukocyte mtDNAcn and survival in EBC patients, our study has several strengths, including long-term follow-up, more cases than previous studies, and comprehensive clinical data. There are several limitations to our study. First, as all patients were of Chinese origin, it is unclear whether our findings are Chinese Han population-specific or common in other populations. Second, we only performed association analyses between mtDNA content and the prognosis of BC patients. The underlying mechanisms through which mtDNA content affects immune functions need further investigation. Lastly, future correlation analysis of mtDNAcn in different tissues (e.g., leukocytes vs. tumor tissues) would help us better understand the relationship between mitochondrial biology and carcinogenesis.

Conclusions

For the first time, we proposed that the prognostic value of leukocyte mtDNA content in EBC patients may depend on intrinsic tumor subtypes. Our study reminds us that we should pay special attention to HR-positive BC patients with high peripheral blood mitochondrial copy numbers when treating them because the prognosis of HR-positive BC patients may be worse. Once confirmed, leukocyte mtDNA content can serve as a useful biomarker to improve the prognosis prediction for BC patients.

Data Availability

The data that support the findings of this study are available from Fujian Medical University Union Hospital, but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. However, data are available from the authors upon reasonable request and with permission of Fujian Medical University Union Hospital. Fangmeng Fu, the corresponding author, should be contacted if someone wants to request the data from this study.

Abbreviations

mtDNAcn:

mitochondrial DNA copy number

BC:

breast cancer

EBC:

early-stage breast cancer

iDFS:

invasive disease-free survival

DDFS:

distant disease-free survival

BCSS:

BC special survival

OS:

overall survival

HR:

hormone receptor

ROS:

reactive oxygen species

AJCC:

American Joint Commission on Cancer

TNM:

tumor-node-metastasis

ER:

estrogen receptor

PR:

progesterone receptor

HER2:

Human epidermal growth factor-2

HRs:

hazard ratios

CI:

confidence interval

aHR:

adjusted hazard ratio

References

  1. Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, Jemal A, Yu XQ, He J. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66(2):115–32.

    Article  PubMed  Google Scholar 

  2. D’Aiuto F, Callari M, Dugo M, Merlino G, Musella V, Miodini P, Paolini B, Cappelletti V, Daidone MG. miR-30e* is an independent subtype-specific prognostic marker in breast cancer. Br J Cancer. 2015;113(2):290–8.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Loibl S, Poortmans P, Morrow M, Denkert C, Curigliano G. Breast cancer. Lancet. 2021;397(10286):1750–69.

    Article  CAS  PubMed  Google Scholar 

  4. Hunter K, Welch DR, Liu ET. Genetic background is an important determinant of metastatic potential. Nat Genet. 2003;34(1):23–4. author reply 25.

    Article  CAS  PubMed  Google Scholar 

  5. Greaves LC, Reeve AK, Taylor RW, Turnbull DM. Mitochondrial DNA and disease. J Pathol. 2012;226(2):274–86.

    Article  CAS  PubMed  Google Scholar 

  6. Van Houten B, Woshner V, Santos JH. Role of mitochondrial DNA in toxic responses to oxidative stress. DNA Repair (Amst). 2006;5(2):145–52.

    Article  PubMed  Google Scholar 

  7. Veltri KL, Espiritu M, Singh G. Distinct genomic copy number in mitochondria of different mammalian organs. J Cell Physiol. 1990;143(1):160–4.

    Article  CAS  PubMed  Google Scholar 

  8. Shen J, Platek M, Mahasneh A, Ambrosone CB, Zhao H. Mitochondrial copy number and risk of breast cancer: a pilot study. Mitochondrion. 2010;10(1):62–8.

    Article  CAS  PubMed  Google Scholar 

  9. Lee HC, Wei YH. Mitochondrial DNA instability and metabolic shift in human cancers. Int J Mol Sci. 2009;10(2):674–701.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Yu M, Zhou Y, Shi Y, Ning L, Yang Y, Wei X, Zhang N, Hao X, Niu R. Reduced mitochondrial DNA copy number is correlated with tumor progression and prognosis in chinese breast cancer patients. IUBMB Life. 2007;59(7):450–7.

    Article  CAS  PubMed  Google Scholar 

  11. Wu CW, Yin PH, Hung WY, Li AF, Li SH, Chi CW, Wei YH, Lee HC. Mitochondrial DNA mutations and mitochondrial DNA depletion in gastric cancer. Genes Chromosomes Cancer. 2005;44(1):19–28.

    Article  CAS  PubMed  Google Scholar 

  12. Kim MM, Clinger JD, Masayesva BG, Ha PK, Zahurak ML, Westra WH, Califano JA. Mitochondrial DNA quantity increases with histopathologic grade in premalignant and malignant head and neck lesions. Clin Cancer Res. 2004;10(24):8512–5.

    Article  CAS  PubMed  Google Scholar 

  13. van Osch FH, Voets AM, Schouten LJ, Gottschalk RW, Simons CC, van Engeland M, Lentjes MH, van den Brandt PA, Smeets HJ, Weijenberg MP. Mitochondrial DNA copy number in colorectal cancer: between tissue comparisons, clinicopathological characteristics and survival. Carcinogenesis. 2015;36(12):1502–10.

    PubMed  Google Scholar 

  14. Mizumachi T, Suzuki S, Naito A, Carcel-Trullols J, Evans TT, Spring PM, Oridate N, Furuta Y, Fukuda S, Higuchi M. Increased mitochondrial DNA induces acquired docetaxel resistance in head and neck cancer cells. Oncogene. 2008;27(6):831–8.

    Article  CAS  PubMed  Google Scholar 

  15. Hsu CW, Yin PH, Lee HC, Chi CW, Tseng LM. Mitochondrial DNA content as a potential marker to predict response to anthracycline in breast cancer patients. Breast J. 2010;16(3):264–70.

    Article  CAS  PubMed  Google Scholar 

  16. Gentiluomo M, Katzke VA, Kaaks R, Tjonneland A, Severi G, Perduca V, Boutron-Ruault MC, Weiderpass E, Ferrari P, Johnson T, et al. Mitochondrial DNA Copy-Number variation and pancreatic Cancer risk in the prospective EPIC cohort. Cancer Epidemiol Biomarkers Prev. 2020;29(3):681–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Xu J, Chang WS, Tsai CW, Bau DT, Davis JW, Thompson TC, Logothetis CJ, Gu J. Mitochondrial DNA copy number in peripheral blood leukocytes is associated with biochemical recurrence in prostate cancer patients in African Americans. Carcinogenesis. 2020;41(3):267–73.

    Article  CAS  PubMed  Google Scholar 

  18. Yang K, Li X, Forman MR, Monahan PO, Graham BH, Joshi A, Song M, Hang D, Ogino S, Giovannucci EL, et al. Pre-diagnostic leukocyte mitochondrial DNA copy number and colorectal cancer risk. Carcinogenesis. 2019;40(12):1462–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Xu E, Sun W, Gu J, Chow WH, Ajani JA, Wu X. Association of mitochondrial DNA copy number in peripheral blood leukocytes with risk of esophageal adenocarcinoma. Carcinogenesis. 2013;34(11):2521–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Qu F, Chen Y, Wang X, He X, Ren T, Huang Q, Zhang J, Liu X, Guo X, Gu J, et al. Leukocyte mitochondrial DNA content: a novel biomarker associated with prognosis and therapeutic outcome in colorectal cancer. Carcinogenesis. 2015;36(5):543–52.

    Article  CAS  PubMed  Google Scholar 

  21. He X, Qu F, Zhou F, Zhou X, Chen Y, Guo X, Li J, Huang Q, Yang Y, Lyu Z, et al. High leukocyte mtDNA content contributes to poor prognosis through ROS-mediated immunosuppression in hepatocellular carcinoma patients. Oncotarget. 2016;7(16):22834–45.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Chen Y, Zhang J, Huang X, Zhang J, Zhou X, Hu J, Li G, He S, Xing J. High leukocyte mitochondrial DNA content contributes to poor prognosis in glioma patients through its immunosuppressive effect. Br J Cancer. 2015;113(1):99–106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Xia P, An HX, Dang CX, Radpour R, Kohler C, Fokas E, Engenhart-Cabillic R, Holzgreve W, Zhong XY. Decreased mitochondrial DNA content in blood samples of patients with stage I breast cancer. BMC Cancer. 2009;9:454.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Edge SB, Compton CC. The american Joint Committee on Cancer: the 7th edition of the AJCC cancer staging manual and the future of TNM. Ann Surg Oncol. 2010;17(6):1471–4.

    Article  PubMed  Google Scholar 

  25. Goldhirsch A, Winer EP, Coates AS, Gelber RD, Piccart-Gebhart M, Thurlimann B, Senn HJ. Panel m: personalizing the treatment of women with early breast cancer: highlights of the St Gallen International Expert Consensus on the primary therapy of early breast Cancer 2013. Ann Oncol. 2013;24(9):2206–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Du R, Lu C, Jiang Z, Li S, Ma R, An H, Xu M, An Y, Xia Y, Jin L, et al. Efficient typing of copy number variations in a segmental duplication-mediated rearrangement hotspot using multiplex competitive amplification. J Hum Genet. 2012;57(8):545–51.

    Article  CAS  PubMed  Google Scholar 

  27. Gourgou-Bourgade S, Cameron D, Poortmans P, Asselain B, Azria D, Cardoso F, A’Hern R, Bliss J, Bogaerts J, Bonnefoi H, et al. Guidelines for time-to-event end point definitions in breast cancer trials: results of the DATECAN initiative (definition for the Assessment of Time-to-event endpoints in CANcer trials)dagger. Ann Oncol. 2015;26(5):873–9.

    Article  CAS  PubMed  Google Scholar 

  28. Fazzini F, Lamina C, Fendt L, Schultheiss UT, Kotsis F, Hicks AA, Meiselbach H, Weissensteiner H, Forer L, Krane V, et al. Mitochondrial DNA copy number is associated with mortality and infections in a large cohort of patients with chronic kidney disease. Kidney Int. 2019;96(2):480–8.

    Article  CAS  PubMed  Google Scholar 

  29. Koller A, Fazzini F, Lamina C, Rantner B, Kollerits B, Stadler M, Klein-Weigel P, Fraedrich G, Kronenberg F. Mitochondrial DNA copy number is associated with all-cause mortality and cardiovascular events in patients with peripheral arterial disease. J Intern Med. 2020;287(5):569–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Lee HC, Yin PH, Lu CY, Chi CW, Wei YH. Increase of mitochondria and mitochondrial DNA in response to oxidative stress in human cells. Biochem J 2000, 348 Pt 2:425–432.

  31. Shen M, Zhang L, Bonner MR, Liu CS, Li G, Vermeulen R, Dosemeci M, Yin S, Lan Q. Association between mitochondrial DNA copy number, blood cell counts, and occupational benzene exposure. Environ Mol Mutagen. 2008;49(6):453–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Liu CS, Tsai CS, Kuo CL, Chen HW, Lii CK, Ma YS, Wei YH. Oxidative stress-related alteration of the copy number of mitochondrial DNA in human leukocytes. Free Radic Res. 2003;37(12):1307–17.

    Article  CAS  PubMed  Google Scholar 

  33. Jain M, Rivera S, Monclus EA, Synenki L, Zirk A, Eisenbart J, Feghali-Bostwick C, Mutlu GM, Budinger GR, Chandel NS. Mitochondrial reactive oxygen species regulate transforming growth factor-beta signaling. J Biol Chem. 2013;288(2):770–7.

    Article  CAS  PubMed  Google Scholar 

  34. Giannelli G, Villa E, Lahn M. Transforming growth factor-beta as a therapeutic target in hepatocellular carcinoma. Cancer Res. 2014;74(7):1890–4.

    Article  CAS  PubMed  Google Scholar 

  35. Weerts MJ, Sieuwerts AM, Smid M, Look MP, Foekens JA, Sleijfer S, Martens JW. Mitochondrial DNA content in breast cancer: impact on in vitro and in vivo phenotype and patient prognosis. Oncotarget. 2016;7(20):29166–76.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Bai RK, Chang J, Yeh KT, Lou MA, Lu JF, Tan DJ, Liu H, Wong LJ. Mitochondrial DNA content varies with pathological characteristics of breast cancer. J Oncol. 2011;2011:496189.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Rai NK, Panjwani G, Ghosh AK, Haque R, Sharma LK. Analysis of mitochondrial DNA copy number variation in blood and tissue samples of metastatic breast cancer patients (a pilot study). Biochem Biophys Rep. 2021;26:100931.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Chen N, Wen S, Sun X, Fang Q, Huang L, Liu S, Li W, Qiu M. Elevated mitochondrial DNA Copy Number in Peripheral Blood and tissue predict the Opposite Outcome of Cancer: a Meta-analysis. Sci Rep. 2016;6:37404.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. D’Erchia AM, Atlante A, Gadaleta G, Pavesi G, Chiara M, De Virgilio C, Manzari C, Mastropasqua F, Prazzoli GM, Picardi E, et al. Tissue-specific mtDNA abundance from exome data and its correlation with mitochondrial transcription, mass and respiratory activity. Mitochondrion. 2015;20:13–21.

    Article  PubMed  Google Scholar 

  40. Fernandez-Vizarra E, Enriquez JA, Perez-Martos A, Montoya J, Fernandez-Silva P. Tissue-specific differences in mitochondrial activity and biogenesis. Mitochondrion. 2011;11(1):207–13.

    Article  CAS  PubMed  Google Scholar 

  41. Pedram A, Razandi M, Wallace DC, Levin ER. Functional estrogen receptors in the mitochondria of breast cancer cells. Mol Biol Cell. 2006;17(5):2125–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Behera MA, Dai Q, Garde R, Saner C, Jungheim E, Price TM. Progesterone stimulates mitochondrial activity with subsequent inhibition of apoptosis in MCF-10A benign breast epithelial cells. Am J Physiol Endocrinol Metab. 2009;297(5):E1089–1096.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Gao W, Wu MH, Wang N, Ying MZ, Zhang YY, Hua J, Chuan L, Wang YJ. Mitochondrial transcription factor A contributes to cisplatin resistance in patients with estrogen receptorpositive breast cancer. Mol Med Rep. 2016;14(6):5304–10.

    Article  CAS  PubMed  Google Scholar 

  44. Cam H, Balciunaite E, Blais A, Spektor A, Scarpulla RC, Young R, Kluger Y, Dynlacht BD. A common set of gene regulatory networks links metabolism and growth inhibition. Mol Cell. 2004;16(3):399–411.

    Article  CAS  PubMed  Google Scholar 

  45. Ishikawa K, Takenaga K, Akimoto M, Koshikawa N, Yamaguchi A, Imanishi H, Nakada K, Honma Y, Hayashi J. ROS-generating mitochondrial DNA mutations can regulate tumor cell metastasis. Science. 2008;320(5876):661–4.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by Joint Funds for the Innovation of Science and Technology, Fujian Province [grant numbers 2018Y9055, 2019Y9103, 2020Y9053]; Joint Key Funds for the Health and Education of Fujian Province [grant number 2019-WJ-23]; and Joint Key Funds for the Natural Science Foundation of Fujian Province [grant number 2021J01737].

Author information

Authors and Affiliations

Authors

Contributions

Chuan Wang and Fangmeng Fu conceived and designed the study. Yan Li, Shengmei Li, Xuan Jin, Weifeng Cai, Kun Zhang, Qian Nie, Hanxi Chen, Jing Li, Peng He, Qindong Cai, and Yibin Qiu collected the data. Liuwen Yu and Jialin Hou analyzed the data. Xiaobin Chen, Lili Chen, Minyan Chen, Wenhui Guo, and Yuxiang Lin offer technical assistance. Wenzhe Zhang, Songping Lin, and Bangwei Zeng wrote the manuscript. All authors contributed toward data analysis, drafting, and critically revising the paper, gave final approval of the version to be published, and agreed to be accountable for all aspects of the work.

Corresponding authors

Correspondence to Chuan Wang or Fangmeng Fu.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the ethics committee of Fujian Medical University Union Hospital (2019KJTZD003), and written informed consent was obtained from all subjects before their participation. All the procedures performed in studies involving human participants adhere to the World Medical Association (Declaration of Helsinki).

Consent for publication

We have obtained consent to publish from the participants.

Competing interests

The authors declare that they have no conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, W., Lin, S., Zeng, B. et al. High leukocyte mitochondrial DNA copy number contributes to poor prognosis in breast cancer patients. BMC Cancer 23, 377 (2023). https://doi.org/10.1186/s12885-023-10838-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-10838-x

Keywords