Mitochondrial DNA copy number variation, leukocyte telomere length, and breast cancer risk in the European Prospective Investigation into Cancer and Nutrition (EPIC) study
Breast Cancer Research volume 20, Article number: 29 (2018)
Leukocyte telomere length (LTL) and mitochondrial genome (mtDNA) copy number and deletions have been proposed as risk markers for various cancer types, including breast cancer (BC).
To gain a more comprehensive picture on how these markers can modulate BC risk, alone or in conjunction, we performed simultaneous measurements of LTL and mtDNA copy number in up to 570 BC cases and 538 controls from the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort. As a first step, we measured LTL and mtDNA copy number in 96 individuals for which a blood sample had been collected twice with an interval of 15 years.
According to the intraclass correlation (ICC), we found very good stability over the time period for both measurements, with ICCs of 0.63 for LTL and 0.60 for mtDNA copy number. In the analysis of the entire study sample, we observed that longer LTL was strongly associated with increased risk of BC (OR 2.71, 95% CI 1.58–4.65, p = 3.07 × 10− 4 for highest vs. lowest quartile; OR 3.20, 95% CI 1.57–6.55, p = 1.41 × 10− 3 as a continuous variable). We did not find any association between mtDNA copy number and BC risk; however, when considering only the functional copies, we observed an increased risk of developing estrogen receptor-positive BC (OR 2.47, 95% CI 1.05–5.80, p = 0.04 for highest vs. lowest quartile).
We observed a very good correlation between the markers over a period of 15 years. We confirm a role of LTL in BC carcinogenesis and suggest an effect of mtDNA copy number on BC risk.
Mitochondria are responsible for essential functions in the eukaryotic cell, such as energy metabolism, calcium homeostasis, and apoptosis [1, 2]. Mitochondria are also responsible for the generation of ROS, which is a by-product of energy metabolism. They possess their own genome (mitochondrial DNA [mtDNA]), and in each eukaryotic cell, there can be hundreds or even thousands of copies of their genomes . mtDNA copy number reflects the number of mitochondria within each cell, and this number is maintained within a constant range in order to sustain the energetic needs of the cell . Although mtDNA copy number differs by cell type, it has been observed that there is a good correlation between the number of copies of mtDNA across different cell types within the same individual . Thus, measuring mtDNA copy number in blood has been proposed as a reliable and noninvasive method for estimating the average number of mitochondria in .
Over the last decade or so, several reports have suggested that mtDNA content measured in blood may be a risk factor for various diseases, including metabolic diseases, diabetes, obesity, neurodegenerative diseases, and cancer [5,6,7,8,9,10,11,12,13]. The majority of these studies suggest that high mtDNA copy numbers are associated with increased disease risk. One possible explanation for the association of mtDNA copy number with chronic disease risk is oxidative stress. The mitochondrial genome is particularly prone to DNA damage and mutations caused by oxidative stress, and mtDNA mutations are associated with increased mtDNA copy number as a possible compensatory mechanism to cope with mitochondrial dysfunction.
The study of mtDNA copy number variation in breast cancer (BC) is of particular interest, given the role of oxidative stress in the disease’s etiology [14,15,16,17]. Only four studies have investigated the role of mtDNA copy number in BC, and these suggested an association between high copy number and increased risk of developing the disease [5, 9, 13, 18]. All of these studies were focused entirely or predominantly on estrogen receptor (ER)-positive BC.
mtDNA deletions are a heterogeneous assortment of mtDNA mutations in terms of length and position and represent the major contributor to mitochondrial dysfunction . Along the mtDNA circle, there are hotspot positions for deletions, and 90% of these include the loss of the nicotinamide adenine dinucleotide dehydrogenase 4 (ND4) sequence . The most frequent deletion is the “common” deletion of 4977 bp, including the ND4 sequence, which causes primary mitochondrial diseases but has also been associated with age-related disease susceptibility, such as cancer, muscle atrophy, and neurodegeneration [21, 22]. Moreover, in a recent study, Nie and colleagues reported that the common deletion in the mitochondrial genome could increase the risk of developing BC .
Besides oxidative stress, mitochondrial dysfunction (measured as decreased mitochondrial mass and energy production) recently has been associated with telomere attrition and shortening, and mtDNA copy number has been associated with leukocyte telomere length (LTL) . A possible mechanistic link between telomere length and mitochondrial function has been proposed through the peroxisome proliferator-activated receptor-γ pathway . To date, a number of studies have examined LTL in surrogate tissues in relation to cancer risks. However, the results have been inconsistent, showing positive, inverse, or null associations between telomere length and cancer risk, with the majority reporting that shorter telomere length increases the risk (reviewed in [26,27,28,29]). In relation to BC risk, studies have shown conflicting results [26, 29]; however, none of these studies were focused specifically on ER− BC risk.
To date, LTL and mtDNA copy number and deletion levels have been examined mostly as independent contributors to cancer risk. Recent studies have indicated that LTL and mtDNA copy number are positively correlated in healthy individuals, in pregnant women, and in patients with psychological disorders [24, 30, 31]. Moreover, very recently, it has been demonstrated that LTL and mtDNA copy number are correlated in intestinal gastric cancer . Hence, there is emerging evidence that the markers are linked biologically and that their joint measurement may increase their predictive value for cancer risk. Therefore, to gain a more comprehensive picture of how the markers mentioned above can modulate BC risk, alone or in conjunction, we measured mtDNA copy number, LTL, and mtDNA deletions in the context of the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort. In addition, given the fact that ER− BC has been greatly understudied for these markers, we collected a fairly large number of them (n = 251).
Study population: the EPIC cohort
The EPIC cohort has been described in full detail elsewhere . Briefly, EPIC consists of about 520,000 healthy volunteers, aged 35–69 years, who were recruited between 1992 and 2005 in 10 European countries. All EPIC study subjects provided anthropometric measurements (height, weight, and waist and hip circumferences) and extensive standardized questionnaire information about medical history, diet, physical activity, smoking, and other lifestyle factors. The women also answered questions about menstrual and reproductive history, hysterectomy, ovariectomy, and use of exogenous hormones for contraception or treatment of menopausal symptoms. About 260,000 women and 140,000 men provided a blood sample, which was split into aliquots of plasma, serum, buffy coat, and erythrocytes and stored frozen for later laboratory analyses.
Cases of cancer occurring after recruitment into the cohort and blood donation are identified through local and national cancer registries or by a combination of contacts with national health insurance and/or active follow-up through the study subjects or their next of kin. Cancer incidence data are classified according to the International Classification of Diseases, Tenth Revision (ICD-10), system. Incident cases of BC were identified as first occurrence of primary invasive tumors, ICD-10 code C50, occurring among women who had no previous diagnosis of cancer. Since 2001, an expanding series of nested case-control studies have been conducted on hormonal, metabolic, and other blood-based risk factors for BC [34,35,36,37], as well as on genetic determinants [14, 38,39,40,41,42,43,44,45]. For these latter studies, the cases and controls were not individually matched; however, care was taken to select the controls randomly from the cohort to avoid selection bias. Because there were insufficient amounts of biological material (extracted buffy coat DNA) for some of the selected controls, we retained slightly uneven numbers of invasive BC cases. We had measurements on telomere length for 570 controls and 533 cases, measurements on mtDNA copy number on 548 controls and 522 cases, and measurements on mtDNA deletions on 539 controls and 519 cases. This study was approved by the ethical review board of the International Agency for Research on Cancer.
Sample preparation and DNA extraction
All blood specimens were kept frozen in liquid nitrogen at the International Agency for Research on Cancer. DNA was extracted from blood samples on an Autopure instrument (Qiagen, Hilden, Germany) with Puregene chemistry (Qiagen). All of the samples were extracted using the same method and in batches of 96 to avoid any possible bias due to sample handling. The order of DNA from cases and controls was randomized on PCR plates, and each batch cases and controls was chosen to be from the same study center.
LTL and mtDNA copy number/deletion measurement
LTL was measured by real-time qPCR as the number of copies of telomeric repeats compared with a single copy gene (RPLP0) used as a quantitative control. For each sample plate, a standard curve was prepared by using a serially diluted reference DNA (range 0.468–30 ng/μl). The relative length of telomeres (T) compared with the control gene (S) was calculated as the ratio T/S. All the samples were analyzed in triplicate. mtDNA copy number and deletions were also analyzed by qPCR. The mtDNA/nuclear DNA content was assessed using specific primers designed for the mitochondrial ND1 gene and the nuclear β-actin gene. The ratio between these two genes (ND1/β-actin) identifies the relative mtDNA copy number . To quantify mtDNA deletions, primers for the mitochondrial gene ND4, located in the major arch of the mtDNA where deletion frequently occurs, were used in combination with those of the mitochondrial gene ND1, which is rarely deleted. The mtDNA deletion level is reported as the ratio between ND4 and ND1  and thus represented by the quantification of the nondeleted fraction of mtDNA (ND4) vs. the total amount of mtDNA (ND1). As for the LTL measurement, for the analysis of mtDNA, a standard curve was also used in each plate. Four samples were included in all plates for all the assays for quality control.
Longitudinal assessment of LTL and mtDNA copy number
To assess the reliability of the measure of mtDNA content and LTL over time, we compared the measurements of these two markers obtained from samples collected 15 years apart from a set of 96 samples of healthy individuals randomly selected from among the EPIC Heidelberg cohort. We calculated intraclass correlations (ICCs) by dividing the between-person variance by the total variance (sum of between- and within-person variances) . ICCs ≥ 0.75 are considered an indicator of excellent reliability; ICCs between 0.51 and 0.74, good reliability; ICCs between 0.40 and 0.50, fair reliability; and ICCs ≤ 0.40, poor reliability. This threshold has been selected according to a similar study done on 100 EPIC subjects in which researchers evaluated the reliability of serum metabolite measurements .
Quantitative measurements of LTL and mtDNA copy number were log-transformed to obtain variables with an approximately normal distribution. The relationships of age, smoking, and body mass index (BMI) on the biomarkers (LTL and mtDNA copy number/deletion) were examined by linear regression, adjusting for plate number, center, and case-control status. Correlation between the biomarkers (LTL vs. mtDNA copy number, LTL vs. mtDNA, mtDNA copy number vs. deletion) was examined by calculating Pearson’s correlation coefficients and corresponding p values.
To examine the relationships of LTL and mtDNA copy number in relation to BC risk, we used unconditional logistic regression to obtain ORs and 95% CIs adjusted for age, center, smoking status, BMI, and plate. We also performed stratified analyses by ER status. LTL and mtDNA copy number were analyzed as continuous variables, and as ordinal variables using quartile cutpoints for the variable distributions among the controls. To investigate the aggregated effects of mtDNA copy number together with LTL for each individual, the quartiles of the log-transformed values of the markers were added together, creating a score with values ranging from 2 to 8. The score was investigated as a continuous variable accounting for age, center, and plate. Finally, we investigated the ratio between mtDNA deletion level and mtDNA copy number in relation to BC risk. This measure represents the proportion of intact molecules relative to the total number of mtDNA (calculated as ND1/nuclear DNA content) molecules per subject.
Considering the number of statistical tests carried out in the present set of analyses, the threshold for statistical significance had to be corrected to avoid false-positive findings. Although it is not reasonable to assume that the three markers are independent, we preferred a conservative approach and thus considered the p value threshold of (0.05/3) < 0.017.
In the first step, we measured LTL and mtDNA copy number in 96 individuals belonging to the EPIC Heidelberg cohort, from whom a blood sample was collected at 2 time points approximately 15 years apart. Both measurements showed good stability over this 14- to 15-year period, with ICCs of 0.63 for LTL and 0.60 for mtDNA. Given the good reliability of the estimates, we then measured LTL and mtDNA copy number in 570 invasive BC cases and 538 control subjects recruited in the context of the prospective EPIC study. Four samples were included in all plates and analyzed in all assays for quality control. The highest observed coefficient of variation (CV) among these samples and across assays was 5.52%, and the average CV over all samples and all assays was 3.06%.
We also calculated CVs on the seven dilutions of the sample used as standard across genotyping plates and assays. The highest observed CV was 3.46%, and the average CV over all samples and all assays was 2.02%.
Regarding the distribution of LTL and mtDNA measures across different strata of risk covariates (Table 1), age was, as expected, inversely related to LTL (p = 1.66 × 10− 11) but was not associated with the mtDNA genome variables. Smoking did not show any association with the three markers measured, whereas BMI was significantly associated with mtDNA copy number, though with a modest effect size (correlation coefficient 0.009, p = 0.005). Similarly to a previous study , we also observed a modest but statistically significant, correlation between telomere length and mtDNA copy number (correlation coefficient 0.16, p = 7.3 × 10− 4), whereas there was no correlation between the other two markers (Additional file 1: Table S1).
We observed that longer LTL was strongly associated with increased risk of BC when analyzing LTL as a categorical variable (OR 2.71, 95% CI 1.58–4.65, p = 3.07 × 10− 4 for highest vs. lowest quartile) and also as a continuous variable (OR 3.20, 95% CI 1.57–6.55, p = 1.41 × 10− 3). When we stratified the analysis by ER status of cancer, we observed that longer telomeres were associated with increased risk of both ER− BC (OR 2.12, 95% CI 1.04–4.31, p = 3.84 × 10− 2 for highest vs. lowest quartile) and ER+ BC (OR 2.55, 95% CI 1.10–5.86, p = 2.82 × 10− 2 for highest vs. lowest quartile). All results from the fully adjusted analyses on LTL and BC risk are shown in Table 2. (Minimally adjusted models are shown in Additional file 2: Table S2.)
For mtDNA copy numbers, we did not observe any differential association between cases and controls when considering quartiles. Also, the stratified analysis (ER+/ER−) did not show any statistically significant differences between cases and controls (Table 3). The results of the minimally adjusted models are shown in Additional file 3: Table S3. However, we included in the analysis the ratio between mtDNA deletion level and mtDNA copy (ND4/ND1 divided by ND1/nuclear genome). When comparing this measure across cases and controls, we observed an increase in risk with increased number of copies. This difference was statistically significant only for ER+ BC. When analyzed as a categorical variable, unconditional logistic regression models showed a significant association of ER+ BC with mtDNA integrity (OR 2.47, 95% CI 1.05–5.80, for highest vs. lowest quartile) (Table 4). The results of the minimally adjusted models are shown in Additional file 4: Table S4. Finally when analyzing the score (computing the aggregate effect of mtDNA copy number and LTL) as a continuous variable, we observed an increased risk for individuals with a high score compared with individuals with a low score (OR 1.20, 95% CI 1.05–1.37, p = 0.009). This trend was present in all the strata but was significant only in the analysis considering ER+ and ER− together (Table 5).
Given their central role in several essential functions of the eukaryotic cell, mitochondria—and in particular mitochondrial mutations—have been extensively investigated in relation to human cancer. A growing number of somatic mitochondrial mutations have been implicated in several cancer types, including BC [50, 51]. In addition, the common germline variability in the mitochondrial genome (deletions, mitochondrial single-nucleotide polymorphisms, and resulting mitochondrial haplogroups) has been associated with risk of various cancers, such as of the lung, pancreas, and breast [50, 52,53,54]. Several studies have suggested that, to counterbalance mitochondrial dysfunction caused by genetic mutations, the eukaryotic cell increases their number. Mitochondrial copy number has been proposed as a risk marker for various cancer types [5,6,7, 9, 10, 12, 13]. In a small exploratory study on BC, Shen and colleagues reported an association between high mtDNA copy number and increased risk of the disease . This association was subsequently replicated in two independent studies focused primarily on ER+ BC [9, 13] and again by Shen and colleagues .
Telomeres are specialized structures that adorn the telomeric end of the eukaryotic chromosomes and are essential for the correct segregation of chromosomes to daughter cells  as well as for the control of chromosomal stability and regulation of cell growth [56,57,58]. In the past several years, there has been an overwhelming number of studies aimed at understanding the role of telomere length, measured in blood, in relation to cancer risk and progression. The results, however, are often inconclusive, with some articles reporting an association between increased risk of developing various cancer types and longer telomeres and others reporting the opposite, no association at all, or even a nonlinear dose-response association with increased risks only at the two extremes of the LTL distribution (very short and very long LTL) [26,27,28,29].
LTL and mtDNA copy number have largely been examined as independent contributors to cancer risk; yet, there is emerging evidence that the markers are linked biologically or at least that measuring them together increases the efficiency of the estimation of their effect. For example, Bao and colleagues  analyzed LTL and mtDNA as progression markers for hepatocellular carcinoma and observed that the combination of the two was a better indicator than each of the two separately. Moreover, it has been demonstrated that LTL and mtDNA copy number are correlated in intestinal gastric cancer . Recently, an association between mitochondrial dysfunction and telomere length attrition has been suggested, and a correlation between LTL and mtDNA has been observed by several authors [24, 31, 60]. Both telomere length and mitochondrial function have been proposed as markers of aging (although the evidence is much stronger for telomeres). Mice null for either the telomerase reverse transcriptase (Tert) or telomerase RNA component (Terc) gene had marked repression of peroxisome proliferator-activated receptor-γ, coactivator 1 alpha and beta, and the downstream network. Consistent with peroxisome proliferator-activated receptors as master regulators of mitochondrial physiology and metabolism, telomere dysfunction, including shortening, was associated with impaired mitochondrial biogenesis and function, decreased gluconeogenesis, cardiomyopathy, and increased ROS .
Given the mounting evidence of possible combined associations of cancer risk with LTL and mtDNA, we sought to investigate, for the first time in the same study, the two markers with the addition of mtDNA deletion in order to gain a more comprehensive picture of their involvement in BC risk. Given the large number of caveats linked to these measurements as a first preliminary step, we wanted to determine how stable the measurements were. To do so, we isolated DNA from blood samples taken 15 years apart from 96 individuals. We observed a very good stability over time for both LTL (ICC 0.63) and mtDNA copy number (ICC 0.60). These values probably underestimate the real correlations because the starting material was not exactly the same: The DNA collected at baseline was isolated from buffy coat, whereas for the samples collected after 15 years, the extraction was done starting from whole blood.
In the cross-sectional analysis, our results suggest a clear association between longer telomeres and increased risk of BC independent of ER status. The association between longer telomeres and increased BC risk is biologically plausible because longer LTL may be a marker of an actively reproducing cell and therefore a cell that keeps dividing and that is exposed to an increased risk of acquiring tumor-causing mutations. The association between longer or shorter LTL and cancer risk could also be tumor-specific and may reflect the complex interaction between the genetic determination of telomere length and the environmental risk factors specific for each cancer type.
The analyses of our data did not show any strong association between mtDNA, when analyzed alone, and BC risk, contrary to what has been observed by others [5, 9, 13, 18]. Authors of two recent meta-analyses investigated the role of mtDNA copy number in a large number of studies across different tumors and highlighted the great heterogeneity of the mtDNA effect in different cancer types and a null overall effect [61, 62]. Interestingly, also the four studies in which researchers reported the association with BC showed a high degree of heterogeneity, possibly explaining the difference between those studies and our findings. However, when normalizing for mitochondrial deletions and using the mitochondrial integrity variable, we did observe an increase in risk with a higher number of functional copies. This association is biologically plausible because an increase in mtDNA copy number may reflect increased oxidative stress that in turn may increase damage levels, fuel inflammation , and finally increase the risk of developing BC. One of the novelties of our report is the introduction of a measure of mitochondrial genome integrity and the possibility of observing the effect of functional mitochondrial copy number in cancer risk instead of just measuring the effect mtDNA copy number regardless of their status. Given that when considering mtDNA copy number alone we did not see any association with BC risk, we hypothesize that this “correction,” at least in part, could explain the differences and heterogeneity observed across previous reports. Analyzing the score of mtDNA and LTL, we observed an association with increased risk for the individuals with a high score; however, this probably reflects the association between LTL and BC risk.
This study has several clear strengths, one being its prospective design, which, to a certain extent, may reduce the possibility of reverse causation, which is a critical point when looking for the relationship between LTL, mtDNA, and cancer risk. A second strength is the possibility of avoiding possible confounders by adjusting for all known potential confounding factors on which information was available, such as smoking, BMI, and alcohol consumption. In addition, our data suggest that our calculations of the LTL and mtDNA ratios are unlikely to be heavily affected by sample handling or technical errors, given that even at 15 years apart the measurements remained comparable. A possible limitation of the present study is the relatively small sample size, which is partially counterbalanced, however, by the rather large number of ER− BC cases.
We observed a very good correlation between the markers across a period of 15 years, highlighting how sample handling is crucial for this kind of analysis. Moreover, we also observed a strong and consistent association between longer telomere length and increased BC risk, and additionally, we found a novel association between the mitochondrial genome stability and ER+ BC risk.
Body mass index
Coefficient of variation
European Prospective Investigation into Cancer and Nutrition
International Classification of Diseases, Tenth Revision
Leukocyte telomere length
Nicotinamide adenine dinucleotide dehydrogenase
Wallace DC. A mitochondrial paradigm of metabolic and degenerative diseases, aging, and cancer: a dawn for evolutionary medicine. Annu Rev Genet. 2005;39:359–407.
Wallace DC. Bioenergetics and the epigenome: interface between the environment and genes in common diseases. Dev Disabil Res Rev. 2010;16(2):114–9.
Clay Montier LL, Deng JJ, Bai Y. Number matters: control of mammalian mitochondrial DNA copy number. J Genet Genomics. 2009;36(3):125–31.
Veltri KL, Espiritu M, Singh G. Distinct genomic copy number in mitochondria of different mammalian organs. J Cell Physiol. 1990;143(1):160–4.
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.
Bonner MR, Shen M, Liu CS, Divita M, He X, Lan Q. Mitochondrial DNA content and lung cancer risk in Xuan Wei, China. Lung Cancer. 2009;63(3):331–4.
Hosnijeh FS, Lan Q, Rothman N, San Liu C, Cheng WL, Nieters A, Guldberg P, Tjonneland A, Campa D, Martino A, et al. Mitochondrial DNA copy number and future risk of B-cell lymphoma in a nested case-control study in the prospective EPIC cohort. Blood. 2014;124(4):530–5.
Lan Q, Lim U, Liu CS, Weinstein SJ, Chanock S, Bonner MR, Virtamo J, Albanes D, Rothman N. A prospective study of mitochondrial DNA copy number and risk of non-Hodgkin lymphoma. Blood. 2008;112(10):4247–9.
Lemnrau A, Brook MN, Fletcher O, Coulson P, Tomczyk K, Jones M, Ashworth A, Swerdlow A, Orr N, Garcia-Closas M. Mitochondrial DNA copy number in peripheral blood cells and risk of developing breast cancer. Cancer Res. 2015;75(14):2844–50.
Lynch SM, Weinstein SJ, Virtamo J, Lan Q, Liu CS, Cheng WL, Rothman N, Albanes D, Stolzenberg-Solomon RZ. Mitochondrial DNA copy number and pancreatic cancer in the Alpha-Tocopherol Beta Carotene Cancer Prevention. Cancer Prev Res (Phila). 2011;4(11):1912–9.
Qu F, Liu X, Zhou F, Yang H, Bao G, He X, Xing J. Association between mitochondrial DNA content in leukocytes and colorectal cancer risk: a case-control analysis. Cancer. 2011;117(14):3148–55.
Yu M. Generation, function and diagnostic value of mitochondrial DNA copy number alterations in human cancers. Life Sci. 2011;89(3–4):65–71.
Thyagarajan B, Wang R, Nelson H, Barcelo H, Koh WP, Yuan JM. Mitochondrial DNA copy number is associated with breast cancer risk. PLoS One. 2013;8(6):e65968.
Blein S, Berndt S, Joshi AD, Campa D, Ziegler RG, Riboli E, Cox DG. Factors associated with oxidative stress and cancer risk in the Breast and Prostate Cancer Cohort Consortium. Free Radic Res. 2014;48(3):380–6.
Cox DG, Dostal L, Hunter DJ, Le Marchand L, Hoover R, Ziegler RG, Thun MJ. N-acetyltransferase 2 polymorphisms, tobacco smoking, and breast cancer risk in the breast and prostate cancer cohort consortium. Am J Epidemiol. 2011;174(11):1316–22.
Gonenc A, Erten D, Aslan S, Akinci M, Simsek B, Torun M. Lipid peroxidation and antioxidant status in blood and tissue of malignant breast tumor and benign breast disease. Cell Biol Int. 2006;30(4):376–80.
Sener DE, Gonenc A, Akinci M, Torun M. Lipid peroxidation and total antioxidant status in patients with breast cancer. Cell Biochem Funct. 2007;25(4):377–82.
Shen J, Wan J, Song R, Zhao H. Peripheral blood mitochondrial DNA copy number, length heteroplasmy and breast cancer risk: a replication study. Carcinogenesis. 2015;36(11):1307–13.
Taylor RW, Turnbull DM. Mitochondrial DNA mutations in human disease. Nat Rev. 2005;6(5):389–402.
Damas J, Samuels DC, Carneiro J, Amorim A, Pereira F. Mitochondrial DNA rearrangements in health and disease—a comprehensive study. Hum Mutat. 2014;35(1):1–14.
Brierley EJ, Johnson MA, Lightowlers RN, James OF, Turnbull DM. Role of mitochondrial DNA mutations in human aging: implications for the central nervous system and muscle. Ann Neurol. 1998;43(2):217–23.
Chinnery PF, Samuels DC, Elson J, Turnbull DM. Accumulation of mitochondrial DNA mutations in ageing, cancer, and mitochondrial disease: is there a common mechanism? Lancet. 2002;360(9342):1323–5.
Nie H, Chen G, He J, Zhang F, Li M, Wang Q, Zhou H, Lyu J, Bai Y. Mitochondrial common deletion is elevated in blood of breast cancer patients mediated by oxidative stress. Mitochondrion. 2016;26:104–12.
Kim JH, Kim HK, Ko JH, Bang H, Lee DC. The relationship between leukocyte mitochondrial DNA copy number and telomere length in community-dwelling elderly women. PLoS One. 2013;8(6):e67227.
Sahin E, Colla S, Liesa M, Moslehi J, Muller FL, Guo M, Cooper M, Kotton D, Fabian AJ, Walkey C, et al. Telomere dysfunction induces metabolic and mitochondrial compromise. Nature. 2011;470(7334):359–65.
Hou L, Zhang X, Gawron AJ, Liu J. Surrogate tissue telomere length and cancer risk: shorter or longer? Cancer Lett. 2012;319(2):130–5.
Campa D, Martino A, Varkonyi J, Lesueur F, Jamroziak K, Landi S, Jurczyszyn A, Marques H, Andersen V, Jurado M, et al. Risk of multiple myeloma is associated with polymorphisms within telomerase genes and telomere length. Int J Cancer. 2015;136(5):E351–8.
Campa D, Mergarten B, De Vivo I, Boutron-Ruault MC, Racine A, Severi G, Nieters A, Katzke VA, Trichopoulou A, Yiannakouris N, et al. Leukocyte telomere length in relation to pancreatic cancer risk: a prospective study. Cancer Epidemiol Biomarkers Prev. 2014;23(11):2447–54.
Weischer M, Nordestgaard BG, Cawthon RM, Freiberg JJ, Tybjærg-Hansen A, Bojesen SE. Short telomere length, cancer survival, and cancer risk in 47102 individuals. J Natl Cancer Inst. 2013;105(7):459–68.
Qiu C, Enquobahrie DA, Gelaye B, Hevner K, Williams MA. The association between leukocyte telomere length and mitochondrial DNA copy number in pregnant women: a pilot study. Clin Lab. 2015;61(3–4):363–9.
Tyrka AR, Carpenter LL, Kao HT, Porton B, Philip NS, Ridout SJ, Ridout KK, Price LH. Association of telomere length and mitochondrial DNA copy number in a community sample of healthy adults. Exp Gerontol. 2015;66:17–20.
Jung SJ, Cho JH, Park WJ, Heo YR, Lee JH. Telomere length is correlated with mitochondrial DNA copy number in intestinal, but not diffuse, gastric cancer. Oncol Lett. 2017;14(1):925–9.
Riboli E, Hunt KJ, Slimani N, Ferrari P, Norat T, Fahey M, Charrondiere UR, Hemon B, Casagrande C, Vignat J, et al. European Prospective Investigation into Cancer and Nutrition (EPIC): study populations and data collection. Public Health Nutr. 2002;5(6b):1113–24.
Kaaks R, Berrino F, Key T, Rinaldi S, Dossus L, Biessy C, Secreto G, Amiano P, Bingham S, Boeing H, et al. Serum sex steroids in premenopausal women and breast cancer risk within the European Prospective Investigation into Cancer and Nutrition (EPIC). J Natl Cancer Inst. 2005;97(10):755–65.
Kaaks R, Johnson T, Tikk K, Sookthai D, Tjonneland A, Roswall N, Overvad K, Clavel-Chapelon F, Boutron-Ruault MC, Dossus L, et al. Insulin-like growth factor I and risk of breast cancer by age and hormone receptor status—a prospective study within the EPIC cohort. Int J Cancer. 2014;134(11):2683–90.
Kaaks R, Rinaldi S, Key TJ, Berrino F, Peeters PH, Biessy C, Dossus L, Lukanova A, Bingham S, Khaw KT, et al. Postmenopausal serum androgens, oestrogens and breast cancer risk: the European Prospective Investigation into Cancer and Nutrition. Endocr Relat Cancer. 2005;12(4):1071–82.
Kaaks R, Tikk K, Sookthai D, Schock H, Johnson T, Tjønneland A, Olsen A, Overvad K, Clavel-Chapelon F, Dossus L, et al. Premenopausal serum sex hormone levels in relation to breast cancer risk, overall and by hormone receptor status—results from the EPIC cohort. Int J Cancer. 2014;134(8):1947–57.
Barrdahl M, Canzian F, Joshi AD, Travis RC, Chang-Claude J, Auer PL, Gapstur SM, Gaudet M, Diver WR, Henderson BE, et al. Post-GWAS gene-environment interplay in breast cancer: results from the Breast and Prostate Cancer Cohort Consortium and a meta-analysis on 79,000 women. Hum Mol Genet. 2014;23(19):5260–70.
Barrdahl M, Canzian F, Lindstrom S, Shui I, Black A, Hoover RN, Ziegler RG, Buring JE, Chanock SJ, Diver WR, et al. Association of breast cancer risk loci with breast cancer survival. Int J Cancer. 2015;137(12):2837–45.
Campa D, Barrdahl M, Gaudet MM, Black A, Chanock SJ, Diver WR, Gapstur SM, Haiman C, Hankinson S, Hazra A, et al. Genetic risk variants associated with in situ breast cancer. Breast Cancer Res. 2015;17:82.
Campa D, Claus R, Dostal L, Stein A, Chang-Claude J, Meidtner K, Boeing H, Olsen A, Tjønneland A, Overvad K, et al. Variation in genes coding for AMP-activated protein kinase (AMPK) and breast cancer risk in the European Prospective Investigation on Cancer (EPIC). Breast Cancer Res Treat. 2011;127(3):761–7.
Campa D, Husing A, McKay JD, Sinilnikova O, Vogel U, Tjonneland A, Overvad K, Stegger J, Clavel-Chapelon F, Chabbert-Buffet N, et al. The INSIG2 rs7566605 polymorphism is not associated with body mass index and breast cancer risk. BMC Cancer. 2010;10:563.
Campa D, McKay J, Sinilnikova O, Husing A, Vogel U, Hansen RD, Overvad K, Witt PM, Clavel-Chapelon F, Boutron-Ruault MC, et al. Genetic variation in genes of the fatty acid synthesis pathway and breast cancer risk. Breast Cancer Res Treat. 2009;118(3):565–74.
Garcia-Closas M, Couch FJ, Lindstrom S, Michailidou K, Schmidt MK, Brook MN, Orr N, Rhie SK, Riboli E, Feigelson HS, et al. Genome-wide association studies identify four ER negative-specific breast cancer risk loci. Nat Genet. 2013;45(4):392–8.e2.
Siddiq A, Couch FJ, Chen GK, Lindstrom S, Eccles D, Millikan RC, Michailidou K, Stram DO, Beckmann L, Rhie SK, et al. A meta-analysis of genome-wide association studies of breast cancer identifies two novel susceptibility loci at 6q14 and 20q11. Hum Mol Genet. 2012;21(24):5373–84.
Altilia S, Santoro A, Malagoli D, Lanzarini C, Ballesteros Alvarez JA, Galazzo G, Porter DC, Crocco P, Rose G, Passarino G, et al. TP53 codon 72 polymorphism affects accumulation of mtDNA damage in human cells. Aging. 2012;4(1):28–39.
He L, Chinnery PF, Durham SE, Blakely EL, Wardell TM, Borthwick GM, Taylor RW, Turnbull DM. Detection and quantification of mitochondrial DNA deletions in individual cells by real-time PCR. Nucleic Acids Res. 2002;30(14):e68.
Rosner B. The intraclass correlation coefficient. In: Fundamentals of biostatistics. 6th ed. Belmont: Thomson-Brooks/Cole; 2006. p. 609–13.
Floegel A, Drogan D, Wang-Sattler R, Prehn C, Illig T, Adamski J, Joost HG, Boeing H, Pischon T. Reliability of serum metabolite concentrations over a 4-month period using a targeted metabolomic approach. PLoS One. 2012;6(6):e21103.
Brandon M, Baldi P, Wallace DC. Mitochondrial mutations in cancer. Oncogene. 2006;25(34):4647–62.
Yadav N, Chandra D. Mitochondrial DNA mutations and breast tumorigenesis. Biochim Biophys Acta. 2013;1836(2):336–44.
Ebner S, Lang R, Mueller EE, Eder W, Oeller M, Moser A, Koller J, Paulweber B, Mayr JA, Sperl W, et al. Mitochondrial haplogroups, control region polymorphisms and malignant melanoma: a study in middle European Caucasians. PLoS One. 2011;6(12):e27192.
Lam ET, Bracci PM, Holly EA, Chu C, Poon A, Wan E, White K, Kwok PY, Pawlikowska L, Tranah GJ. Mitochondrial DNA sequence variation and risk of pancreatic cancer. Cancer Res. 2012;72(3):686–95.
Skonieczna K, Malyarchuk BA, Grzybowski T. The landscape of mitochondrial DNA variation in human colorectal cancer on the background of phylogenetic knowledge. Biochim Biophys Acta. 2012;1825(2):153–9.
McEachern MJ, Krauskopf A, Blackburn EH. Telomeres and their control. Annu Rev Genet. 2000;34:331–58.
de Lange T. Shelterin: the protein complex that shapes and safeguards human telomeres. Genes Dev. 2005;19:2100–10.
de Lange T. How telomeres solve the end-protection problem. Science. 2009;326:948–52.
Martinez P, Blasco MA. Role of shelterin in cancer and aging. Aging Cell. 2010;9:653–66.
Bao D, Ba Y, Zhou F, Zhao J, Yang Q, Ge N, Guo X, Wu Z, Zhang H, Yang H, et al. Alterations of telomere length and mtDNA copy number are associated with overall survival in hepatocellular carcinoma patients treated with transarterial chemoembolization. Cancer Chemother Pharmacol. 2016;78(4):791–9.
Alegria-Torres JA, Velazquez-Villafana M, Lopez-Gutierrez JM, Chagoyan-Martinez MM, Rocha-Amador DO, Costilla-Salazar R, Garcia-Torres L. Association of leukocyte telomere length and mitochondrial DNA copy number in children from Salamanca, Mexico. Genet Test Mol Biomarkers. 2016;20(11):654–9.
Hu L, Yao X, Shen Y. Altered mitochondrial DNA copy number contributes to human cancer risk: evidence from an updated meta-analysis. Sci Rep. 2016;6:35859.
Mi J, Tian G, Liu S, Li X, Ni T, Zhang L, Wang B. The relationship between altered mitochondrial DNA copy number and cancer risk: a meta-analysis. Sci Rep. 2015;5:10039.
Kamp DW, Shacter E, Weitzman SA. Chronic inflammation and cancer: the role of the mitochondria. Oncology (Williston Park). 2011;25(5):400–10. 413
We thank the participants and the investigators of EPIC study.
The coordination of the EPIC study is financially supported by the European Commission (DG-SANCO) and the International Agency for Research on Cancer. The national cohorts are supported by the following organizations: Ligue Contre le Cancer, Institut Gustave Roussy, Mutuelle Générale de l’Education Nationale, Institut National de la Santé et de la Recherche Médicale (INSERM) (France); German Cancer Aid, German Cancer Research Center (DKFZ), Federal Ministry of Education and Research (BMBF) (Germany); the Hellenic Health Foundation, the Stavros Niarchos Foundation (Greece); Associazione Italiana per la Ricerca sul Cancro (AIRC) and National Research Council (Italy); Dutch Ministry of Public Health, Welfare and Sports (VWS), Netherlands Cancer Registry (NKR), LK Research Funds, Dutch Prevention Funds, Dutch ZON (Zorg Onderzoek Nederland), World Cancer Research Fund (WCRF), Statistics Netherlands (The Netherlands); Health Research Fund (FIS) (PI13/00061 to Granada, PI13/01162 to EPIC-Murcia), the regional Governments of Andalucía, Asturias, Basque Country, Murcia, and Navarra, and Instituto de Salud Carlos III (ISCIII) Redes temáticas de investigación cooperativa en salud (RETICS) (RD06/0020) (Spain); and Cancer Research UK (14136 to EPIC-Norfolk; C570/A16491 and C8221/A19170 to EPIC-Oxford) and the Medical Research Council (1000143 to EPIC-Norfolk, MR/M012190/1 to EPIC-Oxford) (United Kingdom).
Availability of data and materials
The datasets on which the present analyses are based are not publicly available, owing to incompatibility with the regulations of the different ethics committees in the contributing study centers. However, these datasets are available from the corresponding author upon reasonable request.
Ethics approval and consent to participate
All individuals signed a written consent form, and the study was reviewed by ethics committees in the different study centers. The following are the ethics committees from the individual EPIC centers: University of Athens Medical School, Norwich District Ethics Committee, the National Committee on Health Research Ethics, Comité de Protection des Personnes, Ethics Committee of the Heidelberg University Medical School, Imperial College Research Ethics Committee, Comitato Etico Indipendente, Fondazione IRCCS Istituto Nazionale dei Tumori, Comitato Etico Locale Azienda Sanitaria di Firenze, Ethics Committee of Lund University, the Medical Ethics Committee of the University Medical Center Utrecht, regional ethics committee for northern Norway and the Norwegian Data Inspectorate, Scotland A Research Ethics Committee, Ethikkommission der Landesärztekammer Brandenburg Cottbus, Comité de Ética de Investigación Clínica (CEIC), Human Genetics Foundation Torino Ethics Committee, Umea Regional Ethical Review Board, and IARC Ethics Committee.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Correlation between the biomarkers in control subjects. (DOCX 13 kb)
Table S2. Associations between LTL and BC risk. (DOCX 15 kb)
Table S3. Associations between mtDNA copy number and BC risk. (DOCX 15 kb)
Table S4. Associations between mtDNA integrity (ND4/ND1 divided by ND1/nuclear genome) and BC risk. (DOCX 15 kb)
About this article
Cite this article
Campa, D., Barrdahl, M., Santoro, A. et al. Mitochondrial DNA copy number variation, leukocyte telomere length, and breast cancer risk in the European Prospective Investigation into Cancer and Nutrition (EPIC) study. Breast Cancer Res 20, 29 (2018). https://doi.org/10.1186/s13058-018-0955-5
- Mitochondrial copy number
- Telomere length
- Breast cancer
- Cancer epidemiology