- Research article
- Open Access
Mitochondrial DNA copy number and incident atrial fibrillation
BMC Medicine volume 18, Article number: 246 (2020)
Mechanistic studies suggest that mitochondria DNA (mtDNA) dysfunction may be associated with increased risk of atrial fibrillation (AF). The association between mtDNA copy number (mtDNA-CN) and incident AF in the general population, however, remains unknown.
We conducted prospective analyses of 19,709 participants from the Atherosclerosis Risk in Communities Study (ARIC), the Multi-Ethnic Study of Atherosclerosis (MESA), and the Cardiovascular Health Study (CHS). mtDNA-CN from the peripheral blood was calculated from probe intensities on the Affymetrix Genome-Wide Human single nucleotide polymorphisms (SNP) Array 6.0 in ARIC and MESA and from multiplexed real-time quantitative polymerase chain reaction (qPCR) in CHS. Incident AF cases were identified through electrocardiograms, review of hospital discharge codes, Medicare claims, and death certificates.
The median follow-up time was 21.4 years in ARIC, 12.9 years in MESA, and 11.0 years in CHS, during which 4021 participants developed incident atrial fibrillation (1761 in ARIC, 790 in MESA, and 1470 in CHS). In fully adjusted models, participants with the lowest quintile of mitochondria DNA copy number had an overall 13% increased risk (95% CI 1 to 27%) of incident atrial fibrillation compared to those with the highest quintile. Dose-response spline analysis also showed an inverse association between mitochondria DNA copy number and hazard for atrial fibrillation for all three cohorts. These associations were consistent across subgroups.
Mitochondria DNA copy number was inversely associated with the risk of AF independent of traditional cardiovascular risk factors. These findings implicate mitochondria DNA copy number as a novel risk factor for atrial fibrillation. Further research is warranted to understand the underlying mechanisms and to evaluate the role of mitochondria DNA copy number in the management of atrial fibrillation risk.
Atrial fibrillation (AF) is the most common form of clinical cardiac arrhythmia, with rising prevalence and incidence worldwide. The lifetime risk of developing AF ranges from 20 to 37% in Whites and Blacks [1,2,3], and it is estimated that the number of adults with AF will double in the USA by the year 2050, affecting more than 8 million people . AF imposes considerable mortality and morbidity risks related to cardiovascular events and thromboembolism and is associated with tremendous healthcare costs [5, 6]. The high lifetime risk and adverse consequences of AF highlight the critical need for identifying novel risk markers that may provide insights into AF prevention and treatment.
Mitochondria generate energy for the cell through converting nutrients and oxygen into adenosine triphosphate (ATP) . Unlike other organelles, mitochondria have their own circular DNA (mtDNA), which encodes essential genes for oxidative phosphorylation. Each cell contains on average 103 to 104 copies of mtDNA, with variations by cell type and development phase . Mitochondrial DNA copy number (mtDNA-CN) is proportional to the transcription of mitochondrial genes and is a marker of mitochondrial dysfunction . Indeed, reduced mtDNA-CN from the peripheral blood is associated with adverse cardiovascular disease (CVD) events including heart failure, all-cause mortality, sudden cardiac death, and atherosclerotic CVD [10,11,12,13], as well as with CVD risk factors, including hypertension, diabetes, atherosclerosis, and chronic kidney disease [14,15,16].
Emerging evidence from mechanistic studies suggests that mtDNA dysfunction may be associated with an increased risk of AF through reduced ATP production and elevated reactive oxygen species [17, 18]. The association between mtDNA-CN and incident AF in the general population, however, remains unknown. In the present study, we examined the prospective association between baseline mtDNA-CN and the risk of incident AF among participants from 3 community-based prospective cohort studies: the Atherosclerosis Risk in Communities (ARIC) study, the Multi-Ethnic Study of Atherosclerosis (MESA), and the Cardiovascular Health Study (CHS).
ARIC is a prospective cohort study of 15,792 men and women 45–64 years of age at baseline (1987–1989) . Participants were randomly selected from 4 communities in the USA: Forsyth County, NC; Jackson, MS; Minneapolis suburbs, MN; and Washington County, MD. We excluded participants who reported race other than Black or White (n = 48), Blacks from the Minnesota and Maryland centers because the numbers are too small for adequate adjustment or within-community comparisons (n = 55), and participants with prevalent CHD (n = 667) or prevalent AF at the time of mtDNA-CN measurement (n = 292). We further excluded participants whose mtDNA-CN was not measured due to sample or measurement availability (n = 4259) and those missing other covariates (n = 322). The final sample included 10,149 participants (Fig. 1).
MESA is a prospective cohort study of 6814 men and women aged 45–84 years of age who were free of clinical cardiovascular disease (CVD) at the baseline visit (2000–2002) . We excluded participants whose mtDNA-CN was not measured due to sample or technique availability (n = 898), who had no information on AF during follow-up (n = 82), or who were missing other covariates (n = 53). The final study sample included 5781 participants.
CHS is a prospective cohort study of 5888 men and women aged 65 years and older . The original cohort of 5201 participants was recruited in 1989–1990 from random samples of Medicare eligibility lists, and an additional predominantly African-American cohort of 687 participants was recruited in 1992–1993. We excluded participants who reported race other than Black or White (n = 39), who lacked genetic consent (n = 273), and who had prevalent CHD (n = 1106) or prevalent AF at the time of mtDNA-CN measurement (n = 123). After further excluding participants whose mtDNA-CN was not measured due to sample or technique availability (n = 471), and those missing other covariates (n = 97), the final study sample included 3779 participants.
mtDNA copy number
In ARIC, DNA for mtDNA-CN analysis was collected at visit 1 (1987–1989) for 389 participants (3.8%), at visit 2 (1990–1992) for 8221 participants (81.0%), at visit 3 (1993–1995) for 1485 participants (14.6%), and at visit 4 for (1996–1998) for 54 participants (0.5%). In MESA and CHS, DNA for mtDNA-CN analysis was collected at visit 1. The visit of DNA collection was considered as the baseline visit for each participant in the present analysis.
In ARIC and MESA, DNA samples were isolated from buffy coat and genotyped using Affymetrix Genome-Wide Human SNP Arrays 6.0 (the Genvisis software package [www.genvisis.org]) [11, 12, 22]. Mitochondrial SNPs were collected across all samples and were signaled with high-quality mitochondrial probes. Unadjusted mtDNA-CN was determined as the median of normalized probe intensity differences across all mitochondrial SNPs. To correct for technical artifacts, batch effects, DNA quality, and starting DNA quantity, we applied surrogate variable analysis (ARIC) and principal component analysis (MESA) to probe intensities of 43,316 Affymetrix autosomal SNPs . We then calculated the residuals of a linear regression model with unadjusted mtDNA-CN as the dependent variable and age, sex, enrollment center, technical covariates, and white blood cell (WBC) count as independent variables in ARIC, and age, sex, collection center, race and principal components as independent variables in MESA. In ARIC, WBC count was missing in 14.9% of participants, and we imputed missing WBC as the study mean. In MESA, the rank-based inverse normal transformation was also performed to reduce the impact of outliers.
In CHS, mtDNA-CN was measured using multiplexed real-time quantitative polymerase chain reaction (qPCR) utilizing ABI TaqMan chemistry (Applied Biosystems) . We calculated residuals using a linear mixed effect model stratified by race, with unadjusted mtDNA-CN as the dependent variable, and age, sex, collection site (fixed effects), and qPCR plate (random effect) as independent variables. The mtDNA-CN residuals were standardized within each study (mean of 0 and standard deviation of 1), and this measure was used as our estimate of mtDNA-CN. The normal distribution of mtDNA-CN residuals depended on the goodness of fit of the model. The residuals were standardized within each study (mean of 0 and standard deviation of 1), and this measure was used as our estimate of mtDNA-CN.
In ARIC, AF cases were identified through December 31, 2014, from three sources: electrocardiograms (ECGs) performed during study visits, review of hospital discharge codes, and death certificates . At each study exam, a supine 12-lead resting ECG was performed and transmitted to the ARIC ECG Reading Center (Epidemiological Cardiology Research Center, Wake Forest School of Medicine, Winston Salem, NC) for automatic coding with E Marquette 12-SL program (GE Marquette, Milwaukee, WI). AF or atrial flutter was detected automatically by a computer and confirmed by a cardiologist. Hospitalization information during follow-up was obtained through annual follow-up phone calls and surveillance of local hospitals. Trained abstractors collected hospital discharge codes. AF cases detected in the same hospitalization with open cardiac surgery were excluded. The validity of identifying AF from hospital discharge codes has been established in epidemiological studies [24, 25]. The presence of AF was identified if ICD-9-CM codes 427.31 (AF) or 427.32 (atrial flutter) were listed. Finally, AF was identified from death certificates if ICD-9 427.3 or ICD-10 I48 codes were listed as a cause of death. The AF date was determined as the date of the first ECG with AF (4%), the time of first hospital discharge with AF codes (96%), or when AF was listed as a cause of death (0.1%), whichever occurred first.
In MESA, AF cases were identified through December 31, 2014, from three sources: hospital discharge diagnosis codes, Medicare claims data, and study ECGs . Hospitalization information during follow-up was obtained through phone calls every 9–12 months, and medical records and discharge diagnoses were obtained subsequently. Additionally, for participants enrolled in fee-for-service Medicare, AF diagnoses were identified from inpatient, outpatient, and physician claims. Study ECGs from visit 5 (2010–2012) were also used to identify incident AF.
In CHS, AF cases were identified through December 31, 2012, from three sources: annual ECGs at each study visit through 1999, discharge diagnoses for all hospitalizations (ICD-9-CM code 427.31 or 427.32), and for those enrolled in fee-for-service Medicare, from inpatient, outpatient, or physician claims in Medicare data . The date of AF diagnosis was based on the date of the first ECG indicating AF, the time of the first hospital discharge with AF codes, or the time of the first qualifying outpatient or physician claim, whichever occurred first.
The measurement of other covariates in the three cohorts has been described previously [19,20,21]. Age, sex, race/ethnicity, alcohol intake, smoking status, physical activity, and medication use were self-reported. Alcohol consumption was categorized into never, former, and current for ARIC and MESA and into non-current and current for CHS. Body mass index was calculated as weight (kg) divided by height (m) squared. Hypertension was defined as systolic blood pressure ≥ 140 mmHg, diastolic blood pressure ≥ 90 mmHg, or current use of anti-hypertension medication. Physical activity was assessed via a modified Baecke Questionnaire in ARIC (scales 1–5), a modified Minnesota Leisure Time Physical Activity Questionnaire in CHS (scales 1–4), and as the total amount of intentional moderate or vigorous exercise performed in a usual week in MESA (MET-min/week). Prevalent heart failure was defined by hospital records, physician diagnosis, or self-reported history of treatment.
Plasma total cholesterol, HDL cholesterol, fasting glucose, and creatinine were measured in each study are previously described [19,20,21]. N-terminal pro-brain natriuretic peptide (NT-proBNP) was measured at visits 2 and 4 using an electrochemiluminescent immunoassay in ARIC and at baseline using the Elecsys 2010 analyzer in MESA and CHS. Diabetes was defined as fasting glucose ≥ 126 mg/dL, non-fasting glucose ≥ 200 mg/dL, or use of glycemic control medication. Estimated glomerular filtration rate (eGFR) was calculated using the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation ARIC and MESA , and Modification of Diet in Renal Disease (CKD-MDRD) equation in CHS . Kidney disease was defined as eGFR < 60 mL/min/1.73 m2.
Follow-up started from the baseline visit and continued until the development of AF, death, dropout, or through December 31, 2014, in ARIC and MESA or December 31, 2012, in CHS, whichever occurred first. mtDNA-CN was categorized into cohort-specific quintiles. We used a Cox proportional hazards model to estimate hazard ratios (HR) and 95% confidence intervals (CI) for the association between mtDNA-CN and incident AF in each cohort. HRs compared quintiles 1st to 4th with the 5th quintile (reference). Linear trends across quintiles were tested by including a variable containing the median mtDNA-CN level of each quintile in the models. We also modeled mtDNA-CN as a continuous variable and estimated the HR comparing the 10th to the 90th percentile of mtDNA-CN. In addition, to evaluate non-linear dose-response relationships between mtDNA-CN and incident AF, we modeled mtDNA-CN as restricted cubic splines with knots at the 5th, 35th, 65th, and 95th percentiles of its distribution. Finally, we tested for potential interactions by age, sex, race, smoking, alcohol intake, BMI, hypertension, diabetes, and kidney disease.
All analyses were conducted separately in each cohort, and cohort-specific HRs were combined using a fixed-effects meta-analysis approach. In each cohort, we used 4 multivariate models with progressive degrees of adjustment. Model 1 was adjusted for age, sex, and race enrollment center groups. Model 2 was further adjusted for body mass index, height, smoking, alcohol intake, and physical activity. Model 3 was further adjusted for total and HDL cholesterol, cholesterol medication, hypertension, diabetes, and prevalent HF. Model 4 was further adjusted for NT-proBNP. As additional analyses, we computed missing mtDNA-CN in ARIC using multiple imputation with chained equation (MICE) based on age, sex, race/center groups, BMI, height, physical activity, total cholesterol, HDL cholesterol, hyperlipidemia medication, eGFR, BNP, smoking, alcohol consumption, diabetes, hypertension, and HF. We used a series of 20 imputations to derive the additional estimates. We also conducted a post hoc analysis using the inverse probability of attrition weighting (IPAW) to account for potential informative death censoring [28,29,30]. Time-varying predicted probabilities of the remaining alive at a 6-month binned follow-up interval for each participant were calculated using the Cox proportional hazards model with death as the outcome and all variables in model 3 as independent variables. Attrition weights were calculated as the inverse of the predicted probability.
Statistical analyses were performed using Stata version 15 for ARIC and MESA study and Stata version 12 for the CHS study (StataCorp LP, College Station, TX). All p values were 2-sided, and statistical significance was declared at p < 0.05.
The study included 19,709 participants in total (10,149 from ARIC, 5781 from MESA, and 3779 from CHS). The mean (standard deviation) age of study participants was 57.2 (5.9), 62.3 (10.3), and 72.2 (5.3) for ARIC, MESA, and CHS, respectively (Table 1); 43.5% of the participants were men (ARIC 42.4%, MESA 47.8%, and CHS 39.9%), and 69.4% participants were Whites (ARIC 79.0%, MESA 42.6%, and CHS 84.7%). Participants with lower mtDNA-CN were more likely to be current smokers, to have higher NT-proBNP, and to have diabetes and prevalent CKD (Additional files: Tables S1-S3). Age, male gender, BMI, current smoking, diabetes, and NT-proBNP are the most strongly associated with the risk of incident AF (Table 2).
The median follow-up time was 21.4 (IQR 14.8–23.3) years in ARIC, 12.9 (9.6–13.6) years in MESA, and 11.0 (5.9–17.1) years in CHS. During follow-up, 4021 participants developed incident AF (1761 in ARIC, 790 in MESA, and 1470 in CHS). The overall hazard ratios for incident AF for the 1st to 4th quintiles of mtDNA-CN compared to the 5th quintile are shown in Table 3. In the meta-analysis across the 3 cohorts using the fully adjusted model, participants in the 1st quintile of mitochondria DNA-CN had a 13% increased risk of incident AF compared to those in the 5th quintile (overall hazard ratio 1.13 [1.01, 1.27]) (Fig. 2). In models with mtDNA-CN as a continuous variable, participants at the 10th percentile of mtDNA-CN had a 13% increased risk of incident AF compared to those at the 90th percentile (overall hazard ratio 1.13 [1.04, 1.24]). Dose-response spline analysis also showed an inverse association between mtDNA-CN and AF for all three cohorts, with an approximately linear trend (the p values for non-linear spline terms were 0.20, 0.31, and 0.11 in ARIC, MESA, and CHS, respectively; Fig. 3).
In the stratified analysis, there was no evidence of interaction for the associations between mtDNA-CN and AF across all subgroups evaluated (Additional Files: Figure S1), except that they were stronger in hypertensive than in non-hypertensive participants from ARIC (p interaction = 0.03). Sensitivity analysis excluding non-Black and non-White participants in MESA, as well as additional analyses using MICE for missing mtDNA-CN in ARIC and IPAW to account for informative censoring in all cohorts, showed similar results (Additional Files: Table S4-S6).
In three large population-based prospective cohort studies, mtDNA-CN was inversely associated with the risk of incident AF, independent of traditional risk factors. The association was not statistically different across race and sex groups. This novel association indicates a potential role of mitochondrial dysfunction in atrial arrhythmias and adds to the pathophysiological evidence from basic science studies supporting the role of mitochondrial mechanisms in the genesis of AF.
Animal models and molecular studies suggest that mitochondrial dysfunction is associated with adverse CVD outcomes and subclinical atherosclerosis [31, 32], and some of these associations have been confirmed in population-based studies. However, previous studies on mtDNA-CN and other CVD outcomes reported their results using different metrics, and the magnitude of reported estimates was incomparable. In ARIC, MESA, and CHS, low levels of mtDNA-CN were associated with increased risk of incident CVD (HR for a 1-SD decrease in mtDNA-CN, 1.29 [95% CI, 1.24–1.33]), coronary heart disease (HR for a 1-SD decrease in mtDNA-CN, 1.29 [1.24–1.33]), sudden cardiac death (HR comparing the 1st to the 5th quintiles of mtDNA-CN, 2.24 [1.58–3.19]), and all-cause mortality (HR comparing the 1st to the 5th quintiles of mtDNA-CN, 1.47 [95% CI, 1.33–1.62]) [10,11,12]. Low levels of mtDNA-CN were also associated with CVD risk factors, including hypertension, diabetes, and chronic kidney disease [14,15,16]. Since clinical CVD events and CVD risk factors are also risk factors for AF, the inverse associations between mtDNA-CN and AF could also be mediated through the traditional CVD pathways. As participants in our analysis were free of CHD at baseline and we adjusted for CVD risk factors, we identified an independent and statistically significant association between mtDNA-CN and incident AF after accounting for CVD and traditional risk factors. Furthermore, since low mtDNA-CN levels preceded the development of AF and other CVD events in ARIC, CHS, and MESA, maintaining CVD health may require adequate cell energy mechanisms and mitochondrial function for preserved cardiac contractility, electrical activity, and endothelial function.
AF is triggered by structural and electrophysiological remodeling changes in the atrial myocardium, which in turn adversely affect cardiac function and increase the risk of mortality, stroke, and peripheral embolism . Risk factors of AF include older age, male sex, and the presence of clinical cardiomyopathy, coronary artery disease, and traditional CVD risk factors (hypertension, diabetes, obesity, and smoking) [25, 34, 35]. Our results suggest that mtDNA-CN may be a novel risk factor for AF. As low levels of mtDNA-CN are a marker for mitochondrial dysfunction and abnormal ATP production, our findings suggest that myocyte electrical activity could be compromised by mitochondrial dysfunction and insufficient energy supply.
The mechanisms underlying the association between mtDNA-CN and AF are unknown, but previous mechanistic research and in vitro studies provide some leads. Most of the energy for cardiomyocyte electrical activity and cardiac muscle contraction is supplied by the mitochondria through the oxidative phosphorylation pathway . The generation of cellular reactive oxygen species (ROS) associated with mitochondrial dysfunction could increase arrhythmia susceptibility by affecting energy-dissipating ion channels and transporters, including the proteins involved in excitation-contraction coupling. Indeed, disrupted intracellular Ca2+ homeostasis has been shown to contribute to the pathogenesis of AF [17, 18, 36]. Increased ROS can also impair gap junction regulation and affect voltage-gated sodium-potassium channels, which increases electrical heterogeneity and causes early afterdepolarizations .
Apart from ROS, mitochondrial dysfunction impairs ATP synthesis, affecting cardiomyocyte energy metabolism, sarcolemmal and sarcoplasmic ion channel function, intracellular cation homeostasis, and membrane excitability, which are essential in maintaining the electrical activity of cardiac cells . The reduced energy production promotes the opening of sarcolemmal KATP channels, which has arrhythmogenic effects by conferring shortened action potential duration and excitation wavelengths .
Some limitations of this study need to be considered. Differences in measurement methods in mtDNA-CN, AF, and covariates may contribute to the heterogeneity of the results across the 3 cohorts. Reassuringly, the direction of the association of mtDNA-CN and incident AF was consistent in ARIC, CHS, and MESA, adding weight to the validity of our findings. mtDNA-CN was collected from the peripheral blood and was thus not a direct measurement of mitochondrial function in atrial myocytes. However, mtDNA from the peripheral blood has been shown to be strongly correlated with mtDNA from cardiomyocyte (coefficient of correlation > 0.5) [13, 38]. Therefore, the mtDNA-CN from the peripheral blood could be used as a marker for myocardial mitochondrial function. Furthermore, mtDNA-CN was measured only at a single time point, and variability due to changes over time was not captured. Finally, we could not evaluate the association between mtDNA-CN with various subtypes of AF. Since AF was obtained primarily based on hospital discharge codes in all 3 studies, the diagnoses were mostly persistent or permanent forms of AF, whereas asymptomatic paroxysmal AF, the most common type of AF, is often undetected. Future studies with repeated measures of mtDNA markers may provide a better evaluation of longitudinal changes in mtDNA function and its impact on AF risk.
The strengths of this study include the prospective design, the long follow-up, the rigorous quality control procedures of the individual cohorts, the large sample size, and the heterogeneous composition of the study population, including men and women from multiple race/ethnicity groups and a wide age range from middle-aged adults to elderly participants.
In three prospective community-based cohorts, mtDNA-CN levels in the peripheral blood were inversely associated with the risk of AF independent of traditional risk factors. These findings implicate mtDNA-CN as a novel risk factor for AF. Further research is warranted to better understand the underlying mechanisms, to better characterize the dose-response shape of the association, and to evaluate the role of mtDNA-CN in the prevention and management of AF risk.
Availability of data and materials
The data that support the findings of this study are available from the ARIC, MESA, and CHS repository but are not publicly available. Restrictions apply to the availability of these data, which were used under license for the current study.
Atherosclerosis Risk in Communities Study
Cardiovascular Health Study
Multi-Ethnic Study of Atherosclerosis
Mitochondria deoxyribonucleic acid copy number
Reactive oxygen species
Single nucleotide polymorphisms
Lloyd-Jones DM, Wang TJ, Leip EP, Larson MG, Levy D, Vasan RS, D’Agostino RB, Massaro JM, Beiser A, Wolf PA, et al. Lifetime risk for development of atrial fibrillation: the Framingham Heart Study. Circulation. 2004;110(9):1042–6.
Weng LC, Preis SR, Hulme OL, Larson MG, Choi SH, Wang B, Trinquart L, McManus DD, Staerk L, Lin H, et al. Genetic predisposition, clinical risk factor burden, and lifetime risk of atrial fibrillation. Circulation. 2018;137(10):1027–38.
Mou L, Norby FL, Chen LY, O’Neal WT, Lewis TT, Loehr LR, Soliman EZ, Alonso A. Lifetime risk of atrial fibrillation by race and socioeconomic status: ARIC Study (Atherosclerosis Risk in Communities). Circ Arrhythm Electrophysiol. 2018;11(7):e006350.
Colilla S, Crow A, Petkun W, Singer DE, Simon T, Liu X. Estimates of current and future incidence and prevalence of atrial fibrillation in the U.S. adult population. Am J Cardiol. 2013;112(8):1142–7.
Chugh SS, Havmoeller R, Narayanan K, Singh D, Rienstra M, Benjamin EJ, Gillum RF, Kim YH, McAnulty JH Jr, Zheng ZJ, et al. Worldwide epidemiology of atrial fibrillation: a Global Burden of Disease 2010 Study. Circulation. 2014;129(8):837–47.
Fuster V, Ryden LE, Cannom DS, Crijns HJ, Curtis AB, Ellenbogen KA, Halperin JL, Kay GN, Le Huezey JY, Lowe JE, et al. 2011 ACCF/AHA/HRS focused updates incorporated into the ACC/AHA/ESC 2006 guidelines for the management of patients with atrial fibrillation: a report of the American College of Cardiology Foundation/American Heart Association Task Force on practice guidelines. Circulation. 2011;123(10):e269–367.
Friedman JR, Nunnari J. Mitochondrial form and function. Nature. 2014;505(7483):335–43.
Rooney JP, Ryde IT, Sanders LH, Howlett EH, Colton MD, Germ KE, Mayer GD, Greenamyre JT, Meyer JN. PCR based determination of mitochondrial DNA copy number in multiple species. Methods Mol Biol. 2015;1241:23–38.
Malik AN, Czajka A. Is mitochondrial DNA content a potential biomarker of mitochondrial dysfunction? Mitochondrion. 2013;13(5):481–92.
Zhang Y, Guallar E, Ashar FN, Longchamps RJ, Castellani CA, Lane J, Grove ML, Coresh J, Sotoodehnia N, Ilkhanoff L, et al. Association between mitochondrial DNA copy number and sudden cardiac death: findings from the Atherosclerosis Risk in Communities study (ARIC). Eur Heart J. 2017;38(46):3443–8.
Ashar FN, Zhang Y, Longchamps RJ, Lane J, Moes A, Grove ML, Mychaleckyj JC, Taylor KD, Coresh J, Rotter JI, et al. Association of mitochondrial DNA copy number with cardiovascular disease. JAMA Cardiol. 2017;2(11):1247–55.
Ashar FN, Moes A, Moore AZ, Grove ML, Chaves PH, Coresh J, Newman AB, Matteini AM, Bandeen-Roche K, Boerwinkle E, et al. Association of mitochondrial DNA levels with frailty and all-cause mortality. J Mol Med (Berl). 2015;93(2):177–86.
Huang J, Tan L, Shen R, Zhang L, Zuo H, Wang DW. Decreased peripheral mitochondrial DNA copy number is associated with the risk of heart failure and long-term outcomes. Medicine (Baltimore). 2016;95(15):e3323.
Tin A, Grams ME, Ashar FN, Lane JA, Rosenberg AZ, Grove ML, Boerwinkle E, Selvin E, Coresh J, Pankratz N, et al. Association between mitochondrial DNA copy number in peripheral blood and incident CKD in the Atherosclerosis Risk in Communities Study. J Am Soc Nephrol. 2016;27(8):2467–73.
Lee HK, Song JH, Shin CS, Park DJ, Park KS, Lee KU, Koh CS. Decreased mitochondrial DNA content in peripheral blood precedes the development of non-insulin-dependent diabetes mellitus. Diabetes Res Clin Pract. 1998;42(3):161–7.
Tang X, Luo YX, Chen HZ, Liu DP. Mitochondria, endothelial cell function, and vascular diseases. Front Physiol. 2014;5:175.
Xie W, Santulli G, Reiken SR, Yuan Q, Osborne BW, Chen BX, Marks AR. Mitochondrial oxidative stress promotes atrial fibrillation. Sci Rep. 2015;5:11427.
Yang KC, Bonini MG, Dudley SC Jr. Mitochondria and arrhythmias. Free Radic Biol Med. 2014;71:351–61.
The Atherosclerosis Risk in Communities (ARIC) Study: design and objectives. The ARIC investigators. Am J Epidemiol 1989, 129(4):687–702.
Bild DE, Bluemke DA, Burke GL, Detrano R, Diez Roux AV, Folsom AR, Greenland P, Jacob DR Jr, Kronmal R, Liu K, et al. Multi-ethnic study of atherosclerosis: objectives and design. Am J Epidemiol. 2002;156(9):871–81.
Fried LP, Borhani NO, Enright P, Furberg CD, Gardin JM, Kronmal RA, Kuller LH, Manolio TA, Mittelmark MB, Newman A, et al. The Cardiovascular Health Study: design and rationale. Ann Epidemiol. 1991;1(3):263–76.
Longchamps R, Castellani C, Newcomb C, Sumpter J, Lane J, Grove M, Guallar E, Pankratz N, Taylor K, Rotter J et al: Evaluation of mitochondrial DNA copy number estimation techniques. 2019.
Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28(6):882–3.
Alonso A, Agarwal SK, Soliman EZ, Ambrose M, Chamberlain AM, Prineas RJ, Folsom AR. Incidence of atrial fibrillation in whites and African-Americans: the Atherosclerosis Risk in Communities (ARIC) study. Am Heart J. 2009;158(1):111–7.
Psaty BM, Manolio TA, Kuller LH, Kronmal RA, Cushman M, Fried LP, White R, Furberg CD, Rautaharju PM. Incidence of and risk factors for atrial fibrillation in older adults. Circulation. 1997;96(7):2455–61.
Heckbert SR, Wiggins KL, Blackshear C, Yang Y, Ding J, Liu J, McKnight B, Alonso A, Austin TR, Benjamin EJ, et al. Pericardial fat volume and incident atrial fibrillation in the Multi-Ethnic Study of Atherosclerosis and Jackson Heart Study. Obesity (Silver Spring). 2017;25(6):1115–21.
Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF 3rd, Feldman HI, Kusek JW, Eggers P, Van Lente F, Greene T, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–12.
JM R, DM F. Correcting for noncompliance and dependent censoring in an AIDS clinical trial with inverse probability of censoring weighted (IPCW) log-rank tests. Biometrics. 2000;56(3):779–88.
SR S, IR W. Review of inverse probability weighting for dealing with missing data. Stat Methods Med Res. 2013;22(3):278–95.
S W, A S, MS vN, M F. Correcting for dependent censoring in routine outcome monitoring data by applying the inverse probability censoring weighted estimator. Stat Methods Med Res. 2018;27(2):323–35.
Stockigt F, Beiert T, Knappe V, Baris OR, Wiesner RJ, Clemen CS, Nickenig G, Andrie RP, Schrickel JW. Aging-related mitochondrial dysfunction facilitates the occurrence of serious arrhythmia after myocardial infarction. Biochem Biophys Res Commun. 2017;493(1):604–10.
Ballinger SW. Mitochondrial dysfunction in cardiovascular disease. Free Radic Biol Med. 2005;38(10):1278–95.
Chugh SS, Blackshear JL, Shen WK, Hammill SC, Gersh BJ. Epidemiology and natural history of atrial fibrillation: clinical implications. J Am Coll Cardiol. 2001;37(2):371–8.
Benjamin EJ, Levy D, Vaziri SM, D’Agostino RB, Belanger AJ, Wolf PA. Independent risk factors for atrial fibrillation in a population-based cohort. The Framingham Heart Study. JAMA. 1994;271(11):840–4.
Magnani JW, Rienstra M, Lin H, Sinner MF, Lubitz SA, McManus DD, Dupuis J, Ellinor PT, Benjamin EJ. Atrial fibrillation: current knowledge and future directions in epidemiology and genomics. Circulation. 2011;124(18):1982–93.
Montaigne D, Marechal X, Lefebvre P, Modine T, Fayad G, Dehondt H, Hurt C, Coisne A, Koussa M, Remy-Jouet I, et al. Mitochondrial dysfunction as an arrhythmogenic substrate: a translational proof-of-concept study in patients with metabolic syndrome in whom post-operative atrial fibrillation develops. J Am Coll Cardiol. 2013;62(16):1466–73.
Brown DA, O’Rourke B. Cardiac mitochondria and arrhythmias. Cardiovasc Res. 2010;88(2):241–9.
Knez J, Lakota K, Bozic N, Okrajsek R, Cauwenberghs N, Thijs L, Knezevic I, Vrtovec B, Tomsic M, Cucnik S, et al. Correlation between mitochondrial DNA content measured in myocardium and peripheral blood of patients with non-ischemic heart failure. Genet Test Mol Biomarkers. 2017;21(12):736–41.
The authors thank the other investigators, the staff, and the participants of the ARIC, MESA, and CHS studies.
This research was supported by grant R01HL131573 from the US National Institutes of Health (Drs. Zhao, Guallar, Longchamps, Castellani, and Arking) and by grant P30AG021334 from the Johns Hopkins University Claude D. Pepper Older Americans Independence Center National Institute on Aging (Dr. Arking). Dr. Sotoodehnia is supported by grant R01HL116747, R01HL111089, and R01HL141989. The Atherosclerosis Risk in Communities study has been funded in whole or in part with Federal funds from the National Heart, Lung, and Blood Institute, National Institutes of Health, Department of Health and Human Services, under Contract nos. HHSN268201700001I, HHSN268201700002I, HHSN268201700003I, HHSN268201700005I, and HHSN268201700004I. The Multi-Ethnic Study of Atherosclerosis is supported by contracts HHSN268201500003I, N01-HC-95159, N01-HC-95160, N01-HC-95161, N01-HC-95162, N01-HC-95163, N01-HC-95164, N01-HC-95165, N01-HC-95166, N01-HC-95167, N01-HC-95168, and N01-HC-95169 from the National Heart, Lung, and Blood Institute and by grants UL1-TR-000040 and UL1-TR-001079 from the National Center for Research Resources. Funding for SHARe genotyping was provided by the National Heart, Lung, and Blood Institute contract N02-HL-64278. Also supported in part by the National Center for Advancing Translational Sciences, CTSI grant UL1TR001881, and the National Institute of Diabetes and Digestive and Kidney Disease Diabetes Research Center (DRC) grant DK063491 to the Southern California Diabetes Endocrinology Research Center. The Cardiovascular Health Study is supported by contracts HHSN268201200036C, HHSN268200800007C, HHSN268201800001C, N01HC55222, N01HC85079, N01HC85080, N01HC85081, N01HC85082, N01HC85083, and N01HC85086 and grants U01HL080295 and U01HL130114 from the National Heart, Lung, and Blood Institute, with additional contribution from the National Institute of Neurological Disorders and Stroke. Additional support was provided by R01AG023629 from the National Institute on Aging and R01HL127659 from the National Heart, Lung, and Blood Institute. A full list of principal CHS investigators and institutions can be found at CHS-NHLBI.org. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. Infrastructure for the CHARGE Consortium is supported in part by the National Heart, Lung, and Blood Institute grant R01HL105756.
Ethics approval and consent to participate
Each study obtained approval from their institutional review boards at all centers (ARIC study: Johns Hopkins University, University of Mississippi, University of Minnesota, University of North Carolina at Chapel Hill, University of South Florida, Wake Forrest University; MESA study: University of California Los Angeles, University of Minnesota, Northwestern University, Columbia University, Wake Forrest University, Johns Hopkins University; CHS study: University of California Davis, Johns Hopkins University, Wake Forest University, University of Pittsburgh), and all participants provided written informed consent.
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.
. Hazard ratios for incident atrial fibrillation by levels of mtDNA copy number in pre-specified subgroups. Age groups were ≤65 and >65 for ARIC and MESA, ≤70 and >70 for CHS. Hazard ratios are for comparing the 10th to the 90th percentile of mtDNA copy number. Models were adjusted for age, sex, race/enrollment center, body mass index, height, smoking, alcohol intake, physical activity, total and HDL cholesterol, cholesterol medication, hypertension, diabetes, prevalent CHD, prevalent heart failure, eGFR and log-transformed NT-proBNP at baseline. Table S1. Baseline characteristics of study participants by mtDNA-CN quintiles in ARIC. Table S2. Baseline characteristics of study participants by mtDNA-CN quintiles in MESA. Table S3. Baseline characteristics of study participants by mtDNA-CN quintiles in CHS. Table S4. Hazard ratios for incident AF by quintiles of mtDNA copy number, among Black and White MESA participants. Table S5. Hazard ratios for incident AF by quintiles of mtDNA copy number in ARIC, using multiple imputation by chained equations (MICE) to impute for missing mtDNA-CN. Table S6. Hazard ratios for incident AF by quintiles of mtDNA copy number, using inverse probability weighting to account for death as dependent censoring.
About this article
Cite this article
Zhao, D., Bartz, T.M., Sotoodehnia, N. et al. Mitochondrial DNA copy number and incident atrial fibrillation. BMC Med 18, 246 (2020). https://doi.org/10.1186/s12916-020-01715-6
- Mitochondria DNA copy number
- Atrial fibrillation