- Open Access
Association of leukocyte mitochondrial DNA copy number with longitudinal C-reactive protein levels and survival in older adults: a cohort study
Immunity & Ageing volume 19, Article number: 62 (2022)
Systemic chronic inflammation occurs with age. The association of the leukocyte mitochondrial DNA copy number, a measure of mitochondrial function in aging, with the temporal profile of serum high-sensitivity C-reactive protein and mortality risk remains uncertain. The objectives of this study were to examine the association of the leukocyte mitochondrial DNA copy number with longitudinal high-sensitivity C-reactive protein levels and the association of the longitudinal high-sensitivity C-reactive protein levels with mortality risk.
This prospective cohort study included 3928 adults aged ≥ 55 years without systemic inflammation in the baseline examination of the Healthy Aging Longitudinal Study in Taiwan, which started in 2009. Each participant received leukocyte mitochondrial DNA copy number measurement using a fluorescence-based quantitative polymerase chain reaction at baseline, serum high-sensitivity C-reactive protein measurements at baseline and the follow-up examination five years later, and the ascertainment of all-cause death (until November 30, 2021). The relationships among the leukocyte mitochondrial DNA copy number, longitudinal serum high-sensitivity C-reactive protein levels, and time to all-cause mortality were examined using the joint longitudinal and survival modeling analysis.
Of the 3928 participants (mean age: 69 years; 2060 [52%] were women), 837 (21%) died during follow-up. In the adjusted analysis, one standard deviation lower natural log-transformed baseline leukocyte mitochondrial DNA copy number was associated with an increase of 0.05 (95% confidence interval [CI], 0.02 to 0.08) standard deviation in serum high-sensitivity C-reactive protein in subsequent years. An increase of 1 standard deviation in instantaneous high-sensitivity C-reactive protein levels was associated with a hazard ratio (HR) for all-cause mortality of 1.22 (95% CI, 1.14 to 1.30). Similar results were obtained after further adjusting for baseline high-sensitivity C-reactive protein levels (HR [95% CI], 1.27 [1.16 to 1.38]) and after excluding those with serum high-sensitivity C-reactive protein above 10 mg/L (HR [95% CI], 1.21[1.11 to 1.31]) or 3 mg/L (HR [95% CI], 1.19 [1.06 to 1.31]) during follow-up.
A lower leukocyte mitochondrial DNA copy number was associated with persistently higher high-sensitivity C-reactive protein levels. Moreover, these higher time-varying high-sensitivity C-reactive protein levels were instantaneously associated with a higher risk of death.
Chronically high blood levels of C-reactive protein (CRP), an acute-phase protein induced by the systemic effects of proinflammatory cytokines , and some proinflammatory cytokines themselves are frequently observed in older individuals, even those who appear healthy [2, 3]. Of note, multiple lines of evidence have indicated that simply having a higher blood level of CRP is associated with greater risks of age-related chronic diseases, chronic health conditions, and multiple adverse health outcomes [2,3,4,5,6]. As the population ages, the prevalence of these inflammaging-related health problems, which is already sizable, is likely to increase.
However, the systemic chronic inflammation that tends to develop with age remains incompletely understood [2, 3]. The systemic chronic inflammation has been theorized as a key feature of the aging process [3, 7,8,9]. Aging is associated with alterations in mitochondria . The mitochondrion has its own genome, which encodes components of the oxidative phosphorylation system and exists in multiple copies [10, 11]. The mitochondrial DNA content of peripheral blood cells, usually measured as the quantity of mitochondrial DNA relative to that of nuclear DNA (mitochondrial DNA copy number [mtDNACN]) [9, 11,12,13,14], often decreases with age and is a potential biomarker of mitochondrial function in aging [15,16,17,18]. Studies have demonstrated that the lower the mtDNACN of blood cells is, the greater the risk of specific age-related diseases is [12, 19,20,21]. These diseases include, but not limited to, coronary heart disease, stroke, kidney function decline, and neurodegenerative disorders. For instance, Ashar et al.  showed that one standard deviation lower leukocytes mtDNACN was associated with greater risks of incident coronary heart disease and stroke. Tin et al.  demonstrated that blood cells mtDNACN at the highest quartile was associated with a lower risk of incident chronic kidney disease. These findings suggest that blood cells mtDNACN is inversely related to the risk of these diseases. Notably, chronic inflammation is a common feature of many of these diseases, in particular coronary heart disease, ischemic stroke, and kidney function decline [5, 22]. Therefore, an inference from these observations is that a lower blood cell mtDNACN may be related to systemic chronic inflammation. In line with this, accumulating experimental evidence suggests that immune cell mitochondria are active players in immunity and inflammation and may play etiological roles in the inflammatory pathogenic process [23, 24].
The blood levels of high-sensitivity CRP (hs-CRP) are primarily determined by its synthesis rate in hepatocytes, which in turn is determined by the levels of circulating proinflammatory cytokines raised by a distal inflammatory pathogenic process . As such, serum hs-CRP is a commonly used clinical laboratory test, and its high level typically indicates an active inflammatory pathogenic process. Most studies on hs-CRP have focused on levels measured at a single time point. Although its blood levels at a single time point carry rich prognostic information , more recent studies suggested that its temporal profile may reflect the activity of the corresponding inflammatory pathogenic process in a more timely and proportional manner and that it may additionally provide more detailed information about the activity of the process [25, 26].
The aim of this study was twofold (Additional file 1: Supplementary Fig. 1). The first aim was to prospectively evaluate the relationship of the peripheral blood leukocyte mtDNACN with longitudinal hs-CRP levels in a large cohort of older persons without systemic inflammation. The second aim was to assess the relationship between the longitudinal hs-CRP levels and mortality risk. We hypothesized that a lower peripheral blood leukocyte mtDNACN would be associated with persistently higher time-varying hs-CRP levels and that these higher time-varying hs-CRP levels would be associated with a higher risk of death after adjusting for the effects of multiple potential confounders.
Study design and participants
We analyzed prospectively collected data from the Healthy Aging Longitudinal Study in Taiwan (HALST). The details of the HALST were published previously . In brief, the HALST is an ongoing prospective population-based cohort study in Taiwan that started in 2009 with 5663 community-dwelling adults aged 55 years and older enrolled. After being recruited, these participants received a baseline (first wave) examination and a follow-up (second wave) examination five years later. Sociodemographic status, lifestyles, diseases, and health status were assessed during these examinations. In addition, venous blood samples were collected, immediately processed, and maintained at -80 °C until assayed. This study was approved by the institutional review boards of the National Health Research Institutes and participating hospitals.
This study examined HALST participants without systemic inflammation, defined as a serum hs-CRP level of 2 mg per liter or higher [25, 28, 29], at baseline examination, leaving 3928 HALST participants entering this study. Supplementary Table 1 (of Additional file 1) showed the characteristics of the excluded participants. Each participant received measurement of leukocyte mtDNACN at baseline and measurements of serum hs-CRP at baseline and the follow-up examination.
Measurement of peripheral blood leukocyte mitochondrial DNA copy number
The leukocyte mtDNACN was measured in accordance with protocols described in previous studies . In brief, we amplified a mitochondrial gene (ND1 gene) and a nuclear gene (β-globin gene) in the leukocyte total cellular DNA extracted from venous blood samples using a fluorescence-based quantitative polymerase chain reaction. The relative mtDNACN was estimated on the basis of these two genes’ threshold cycle numbers.
Measurement of serum high-sensitivity C-reactive protein
Baseline and follow-up serum hs-CRP levels were measured using a latex-enhanced immunoturbidimetric assay (ADVIA 1800 Chemistry system, Siemens AG, Munich, Germany). The intra-assay coefficient of variation was 5.61%, and the lowest detectable concentration was 0.12 mg/L.
Assessment of mortality
By linking to data from the Bureau of National Health Insurance of Taiwan death registry, any death events that occurred during the follow-up period were identified systematically, and the date of each death was ascertained. We followed each participant from the date of baseline examination (index date) until death or November 30, 2021, whichever occurred first.
Measurement of covariate variables
Other variables were measured during the baseline examinations. These variables were age (years), sex, education level (illiteracy, elementary school, junior high school, ≥ high school), smoking status (current, former, never), body mass index (BMI), dyslipidemia (serum triglycerides, serum high-density lipoprotein [HDL] cholesterol, serum low-density lipoprotein [LDL] cholesterol), and comorbidities (hypertension, cardiometabolic diseases [diabetes mellitus, stroke, cardiovascular disease], chronic kidney disease, lung disease). Further details are presented in the Supplementary Materials (Additional file 1: Supplementary Methods).
The baseline characteristics of the participants were summarized using descriptive statistics (mean [standard deviation] for normally distributed continuous variables, median [interquartile range] for nonnormally distributed continuous variables, and numbers [percentage] for categorical variables). Because both hs-CRP and the mtDNACN have a skewed distribution, both were analyzed as natural log-transformed values throughout the remainder of the analysis.
We jointly modeled the association of the leukocyte mtDNACN with longitudinal serum hs-CRP levels and time to all-cause mortality using maximum likelihood [31, 32]. The modeling details are presented in the Supplementary Materials (Additional file 1: Supplementary Methods). In the longitudinal part of the joint analysis, we modeled serum hs-CRP levels over time using linear mixed models (random-intercept-and-random-slope models). Each model incorporated fixed effects for the leukocyte mtDNACN, time (in years) since baseline examination, time’s interactions with the leukocyte mtDNACN, and random effects for intercept and slope with an unstructured variance–covariance structure. Departure from linearity was detected using restricted cubic splines with three knots. The correlations of longitudinal data for each individual were accommodated by the random effects in these models. We examined the directions and magnitudes of associations of the baseline leukocyte mtDNACN with hs-CRP based on parameter estimates.
In the survival part of the joint analysis, we modeled the association of the instantaneous hs-CRP level at time t (current value parameterization), as estimated by the longitudinal sub-models, with the risk of death at that time by using a Weibull parametric survival model. Based on the parameter estimates in the model, the hazard ratio (HR) per 1 standard deviation (SD) increase in log-transformed hs-CRP was estimated and presented.
Adjustments were made for covariates in the longitudinal (fixed effects for covariates and their interactions with time) and survival sub-models of joint analysis. As suggested in the literature, we identified a priori covariates that could act as potential confounders and a minimally sufficient set of covariates for adjustment by using a directed acyclic graph (DAG; Fig. 1) [33, 34]. In the longitudinal sub-model, the set of covariates was age, sex, educational level, smoking status, obesity (BMI), dyslipidemia, hypertension, and chronic kidney disease. In the survival sub-model, the set of covariates was age, sex, educational level, smoking status, obesity (BMI), dyslipidemia, hypertension, cardiometabolic diseases, lung disease, and chronic kidney disease. The joint models’ assumptions were examined using residual plots. The proportional hazards assumption of the survival part was evaluated using log–log survival curves and observed-versus-predicted survival curves.
We conducted several sensitivity analyses. First, we added the baseline value of hs-CRP in the survival part of the joint analysis. Second, in the survival part of the joint analysis, we assumed that the risk of death at time t depends on the instantaneous hs-CRP level at that time. However, the risk of death at time t may depend on the levels of hs-CRP prior to time t. To examine whether the assumption might affect the results, we additionally modeled the association of the area under the hs-CRP trajectory curve (in mg/L × years) up to time t (cumulative effects parameterization), which captured the effects of the whole inflammation exposure history up to time t, with the risk of death at time t. Third, to minimize the effects of infectious processes, we excluded all participants with serum hs-CRP above 10 mg/L during follow-up (n = 82) and then repeated the analysis . We then repeated the analysis again after excluding those with serum hs-CRP above 3 mg/L during follow-up (n = 319). Fourth, we conducted a sensitivity analysis to evaluate the robustness of our observations against potential unmeasured confounders by calculating the E-value . E-value quantifies the minimal strength (on the risk-ratio scale) of the association of unmeasured confounders with both exposure and outcomes, conditional on the measured confounders, that can fully explain the observed exposure-outcome association.
Statistical significance was defined as a two-sided P value lower than 0.05. Analyses were performed using STATA, version 16.1 (StataCorp LP., College Station, TX, USA) and R, versions 4.1.2 (R Foundation for Statistical Computing).
The characteristics of the study participants at baseline are presented in Table 1. The mean age was 69 years, and 52% of the participants were women. The median leukocyte mtDNACN at baseline was 72.00. Of the 3928 participants, 837 died during a median follow-up period of 10.6 years. Participants who died during the follow-up period had a significantly lower leukocyte mtDNACN at baseline (68.00) than those who survived (73.00).
Relationship of peripheral blood leukocyte MtDNACN with longitudinal hs-CRP levels
The longitudinal part of the joint analysis revealed that among the participants without systemic inflammation at baseline, those with a lower baseline leukocyte mtDNACN were more likely to have higher serum hs-CRP levels in subsequent years (Table 2). After adjustment, the baseline leukocyte mtDNACN remained inversely associated with serum hs-CRP. A 1-SD lower natural log-transformed baseline leukocyte mtDNACN was associated with an increase of 0.05-SD in serum hs-CRP levels throughout the subsequent years. The participants with a baseline leukocyte mtDNACN within the lowest tertile (< 63.5) had 0.11-SD (95% CI, 0.04 to 0.18) increased serum hs-CRP levels compared with those with a baseline leukocyte mtDNACN within the highest tertile (> 81). We did not observe a non-linear relationship.
We observed similar results when we excluded the participants with serum hs-CRP above 10 mg/L during follow-up (Additional file 1: Supplementary Table 2). A 1-SD lower natural log-transformed baseline leukocyte mtDNACN was associated with an increase of 0.05-SD in serum hs-CRP levels throughout the subsequent years. After the further exclusion of all participants with serum hs-CRP above 3 mg/L during follow-up, a 1-SD lower natural log-transformed baseline leukocyte mtDNACN was associated with an increase of 0.04-SD in serum hs-CRP levels (Additional file 1: Supplementary Table 3).
Relationship between the longitudinal hs-CRP levels and mortality risk
Next, we examined whether the trajectory of serum hs-CRP during follow-up was associated with mortality risk (Table 3). The survival part of the joint analysis revealed that a 1-SD increase in instantaneous hs-CRP levels was associated with an HR for all-cause mortality of 1.36 (95% CI, 1.27 to 1.46) in the unadjusted model and 1.22 (95% CI, 1.14 to 1.30) in the fully adjusted model.
When the baseline value of hs-CRP was added in the survival part of the joint analysis, we observed an HR for all-cause mortality of 1.27 (95% CI, 1.16 to 1.38) associated with a 1-SD increase in instantaneous hs-CRP levels. The term for baseline hs-CRP did not reach statistical significance. We likewise observed an association between the entire hs-CRP exposure during the follow-up period (mg/L × years) and the risk of death (Additional file 1: Supplementary Table 4). One standard deviation increase in hs-CRP levels × years was associated with an HR for all-cause mortality of 1.05 (95% CI, 1.04 to 1.07) in the unadjusted model and 1.03 (95% CI, 1.02 to 1.04) in the fully adjusted models. We also examined whether the estimated effects of hs-CRP on mortality risk were explained by hs-CRP levels exceeding clinically relevant thresholds. In general, the analysis yielded similar results after we excluded the participants with serum hs-CRP above 10 mg/L during follow-up and even after the further exclusion of those with serum hs-CRP above 3 mg/L during follow-up (Table 3) (Additional file 1: Supplementary Table 4). The E-values on the risk ratio scales for the unadjusted and adjusted HR estimates and their lower confidence limit are presented in Supplementary Tables 5, 6, and 7 of additional file 1.
In this longitudinal study of community-dwelling older adults, we identified the relationship of peripheral blood leukocyte mtDNACN with longitudinal serum hs-CRP levels and the relationship of the longitudinal serum hs-CRP levels with the risk of death. This study was the first to demonstrate that even without accompanying signs of systemic inflammation, a lower peripheral blood leukocyte mtDNACN at baseline was associated with persistently higher hs-CRP levels. Moreover, our analysis indicated that these higher longitudinal hs-CRP levels were associated with a higher risk of death. These associations were independent of baseline hs-CRP levels, infectious processes, high-risk cardiovascular inflammation and potential confounders.
Blood levels of hs-CRP are a well-recognized and extensively used biomarker of inflammation [1, 25, 26, 35, 37, 38]. As an acute-phase protein, CRP is rapidly synthesized in hepatocytes, stimulated mainly by circulating cytokine interleukin-6, and released into circulation in response to a range of inflammatory pathologies at distant sites in a dose–response manner [1, 39, 40]. With a constant half-life in the circulation, the level of hs-CRP in an individual’s blood is mainly determined by its rate of synthesis, which in turn is determined by the activity of the underlying cytokines-producing inflammatory pathologies . Therefore, an examination of an individual’s temporal profile of serum hs-CRP by measuring the levels of hs-CRP in their blood serially—rather than at a single time point, as in most studies—offers a unique opportunity to reliably track the activity of the underlying inflammatory pathogenic process [1, 25, 26, 37, 38].
To the best of our knowledge, this study is the first to examine the relationship between the peripheral blood leukocyte mtDNACN and the temporal profile of serum hs-CRP. We demonstrated that among older adults free of clinically relevant systemic inflammation, a lower peripheral blood leukocyte mtDNACN was associated with a persistently higher level of hs-CRP in the subsequent years. This observation was consistent across the multiple analyses that took into account the potential effects of confounders and intercurrent hs-CRP-raising disease processes. These results suggest that a lower leukocyte mtDNACN may be associated with an active inflammatory pathogenic process throughout the years. In general, the observation of an inverse relationship between the peripheral blood leukocyte mtDNACN and serum hs-CRP is consistent with the results of other studies, although most analyses have been cross-sectional in nature and thus have precluded the determination of the chronicity of hs-CRP elevation. For instance, in another study, we revealed that the peripheral blood leukocyte mtDNACN was negatively correlated with hs-CRP in a cohort of 1990 older adults . In a cohort of 9058 participants without chronic kidney disease , Tin et al. noted that those with lower peripheral blood cells mtDNACN had a higher serum hs-CRP level. In another cohort of 4812 patients with chronic kidney disease , those with a lower whole blood mtDNACN were also observed to have a higher hs-CRP.
In this study, by conducting the joint longitudinal and survival analysis, we further demonstrated that a higher level of hs-CRP during follow-up was instantaneously associated with a higher risk of death independent of confounders. Most importantly, this effect of follow-up hs-CRP levels on death risk was independent of baseline hs-CRP levels. On the other hand, the baseline hs-CRP levels were not associated with the death risk after adjusting for the longitudinal hs-CRP levels. Interestingly, after excluding individuals with follow-up hs-CRP levels exceeding thresholds that may indicate an intercurrent infection or high-risk cardiovascular inflammation, the study sample consisted of older adults likely at average or lower risk of diseases (e.g., cardiovascular diseases). Nevertheless, the effects were still observed in this sample of older adults. In sensitivity analysis, we examined whether a different analytic assumption about how the death risk depends on the longitudinal hs-CRP levels might affect the results. Assuming that the risk of death at time t depends on the whole hs-CRP exposure history up to time t, we still observed significant, albeit small, effects of hs-CRP level on the risk of death. The overall results of the survival analysis agree with the argued values of serially determining hs-CRP levels and imply that the longitudinal hs-CRP levels may carry additional information about the near-term health event hazard. Even below clinically relevant thresholds, the chronic low-grade inflammatory process associated with a lower leukocyte mtDNACN might uniquely confer an immediately higher mortality risk in older adults, including low-risk individuals. However, we could not exclude the possibility that the chronic low-grade inflammatory process caused by factors other than lower leukocyte mtDNACN confers a higher mortality risk.
The findings of this longitudinal study are consistent with recent results of experimental investigations that have suggested a causal role of blood cell mitochondria in the development of pathogenic systemic chronic inflammation. First, emerging evidence from human studies has indicated that a lower blood cells mtDNACN may reflect a lower expression of many mitochondrial- and nuclear-encoded genes involved in mitochondrial DNA replication and function in blood cells . Second, as the hub of multiple metabolic and signaling pathways, mitochondria in immune cells play pivotal roles in shaping immune and inflammatory responses . Third, experiments in mice have demonstrated that decreased activity in immune cells' mitochondria metabolic pathways, including mitochondrial respiration, drives systemic chronic inflammation in aging, thereby, the development of age-related diseases .
The major strengths of this study include its prospective design with serial measurements of hs-CRP and the use of an advanced analytical method suitable for studying events with a mechanism of significant chronicity (e.g., aging). In addition, this study analyzed prospectively collected data from a well-established cohort with comprehensive baseline information obtained through a standardized procedure. Notably, the measurements of the mtDNACN and hs-CRP in this study were centralized and standardized. All these characteristics helped this study minimize the likelihood of observing spurious associations. Nevertheless, our study has some limitations. First, as an observational study, inferences regarding causality and mechanism should be drawn with caution. Second, although we identified potential confounders a priori and, guided by DAGs, attempted to adjust for confounding bias during model building, residual confounding by unmeasured or uncontrolled covariates still could not be excluded. We could not adjust for the potential confounding effects of time-varying covariates on the association between higher longitudinal hs-CRP levels and a higher risk of death. To assess how likely the observed association could be due to uncontrolled confounding, we calculated the E-value in the sensitivity analysis. We found that an uncontrolled covariate would need to be associated with both the time-varying hs-CRP levels and mortality risk with a risk ratio of at least 1.56 (lower limit: 1.42) to explain the observed association fully. Baseline stroke was the covariate most strongly associated with mortality, with an HR of 1.74. After excluding stroke from the model, we still obtained an E-value of 1.56 with a lower confidence limit of 1.42. We reasoned that the observed association was unlikely to be due to uncontrolled confounding. Third, the sample of this observational study consisted of adults aged 55 years and older at study entry. We could not completely exclude the possibility of survival bias and the possibility of underestimating the effects of lower leukocyte mtDNACN. Fourth, we excluded those with higher baseline hs-CRP from this study to examine the association of leukocyte mtDNACN with hs-CRP trajectory among those who had not yet developed systemic inflammation at the time of mtDNACN measurement. As such, our findings may not be generalizable to other populations of older adults who might be more ill. Fifth, further research would be worthwhile to determine the clinical relevance of the sub-threshold hs-CRP levels in older adults. Sixth, we showed that the inflammatory process might have a cumulative effect on the risk of death, as demonstrated by the sensitivity analysis. The clinical relevance of such findings is unclear and is needed to be further determined. Finally, both high and low blood cells mtDNACN have been associated with higher risks of a variety of cancers. The roles of inflammation in this context would also be worth investigating.
In conclusion, among older persons without systemic inflammation at baseline, a lower peripheral blood leukocyte mtDNACN was associated with persistently higher hs-CRP levels. Moreover, the higher time-varying hs-CRP levels were immediately associated with higher risks of death. Further studies are warranted to clarify the underlying mechanisms.
Availability of data and materials
The dataset supporting the conclusions of this article is not publicly available according to the Personal Information Protection Act of Taiwan and the regulations of the National Health Research Institutes.
Body mass index
Directed acyclic graph
Healthy Aging Longitudinal Study in Taiwan
- HDL cholesterol:
High-density lipoprotein cholesterol
- LDL cholesterol:
Low-density lipoprotein cholesterol
- mtDNACN :
Mitochondrial DNA copy number
Pepys MB, Hirschfield GM. C-reactive protein: a critical update. J Clin Invest. 2003;111(12):1805–12.
Singh T, Newman AB. Inflammatory markers in population studies of aging. Ageing Res Rev. 2011;10(3):319–29.
Ferrucci L, Fabbri E. Inflammageing: chronic inflammation in ageing, cardiovascular disease, and frailty. Nat Rev Cardiol. 2018;15(9):505–22.
Wang X, Bao W, Liu J, Ouyang YY, Wang D, Rong S, et al. Inflammatory markers and risk of type 2 diabetes: a systematic review and meta-analysis. Diabetes Care. 2013;36(1):166–75.
Kaptoge S, Di Angelantonio E, Lowe G, Pepys MB, Thompson SG, Collins R, et al. C-reactive protein concentration and risk of coronary heart disease, stroke, and mortality: an individual participant meta-analysis. Lancet. 2010;375(9709):132–40.
Newman AB, Murabito JM. The epidemiology of longevity and exceptional survival. Epidemiol Rev. 2013;35(1):181–97.
Franceschi C, Garagnani P, Vitale G, Capri M, Salvioli S. Inflammaging and “Garb-aging.” Trends Endocrinol Metab. 2017;28(3):199–212.
Tchkonia T, Kirkland JL. Aging, cell senescence, and chronic disease: emerging therapeutic strategies. JAMA. 2018;320(13):1319–20.
Lopez-Otin C, Blasco MA, Partridge L, Serrano M, Kroemer G. The hallmarks of aging. Cell. 2013;153(6):1194–217.
Anderson S, Bankier AT, Barrell BG, de Bruijn MH, Coulson AR, Drouin J, et al. Sequence and organization of the human mitochondrial genome. Nature. 1981;290(5806):457–65.
Malik AN, Czajka A. Is mitochondrial DNA content a potential biomarker of mitochondrial dysfunction? Mitochondrion. 2013;13(5):481–92.
Chong M, Mohammadi-Shemirani P, Perrot N, Nelson W, Morton R, Narula S, et al. GWAS and ExWAS of blood Mitochondrial DNA copy number identifies 71 loci and highlights a potential causal role in dementia. Elife. 2022;11:e70382.
Yang SY, Castellani CA, Longchamps RJ, Pillalamarri VK, O’Rourke B, Guallar E, et al. Blood-derived mitochondrial DNA copy number is associated with gene expression across multiple tissues and is predictive for incident neurodegenerative disease. Genome Res. 2021;31(3):349–58.
Khan SS, Singer BD, Vaughan DE. Molecular and physiological manifestations and measurement of aging in humans. Aging Cell. 2017;16(4):624–33.
Mengel-From J, Thinggaard M, Dalgard C, Kyvik KO, Christensen K, Christiansen L. Mitochondrial DNA copy number in peripheral blood cells declines with age and is associated with general health among elderly. Hum Genet. 2014;133(9):1149–59.
Ashar FN, Moes A, Moore AZ, Grove ML, Chaves PH, Coresh J, et al. Association of mitochondrial DNA levels with frailty and all-cause mortality. J Mol Med (Berl). 2015;93(2):177–86.
Liu CS, Tsai CS, Kuo CL, Chen HW, Lii CK, Ma YS, et al. Oxidative stress-related alteration of the copy number of mitochondrial DNA in human leukocytes. Free Radic Res. 2003;37(12):1307–17.
Zhang R, Wang Y, Ye K, Picard M, Gu Z. Independent impacts of aging on mitochondrial DNA quantity and quality in humans. BMC Genomics. 2017;18(1):890.
Ashar FN, Zhang Y, Longchamps RJ, Lane J, Moes A, Grove ML, et al. Association of Mitochondrial DNA Copy Number With Cardiovascular Disease. JAMA cardiology. 2017;2(11):1247–55.
Tin A, Grams ME, Ashar FN, Lane JA, Rosenberg AZ, Grove ML, 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.
Castellani CA, Longchamps RJ, Sun J, Guallar E, Arking DE. Thinking outside the nucleus: Mitochondrial DNA copy number in health and disease. Mitochondrion. 2020;53:214–23.
Mc Causland FR, Claggett B, Burdmann EA, Eckardt KU, Kewalramani R, Levey AS, et al. C-Reactive Protein and Risk of ESRD: Results From the Trial to Reduce Cardiovascular Events With Aranesp Therapy (TREAT). Am J Kidney Dis. 2016;68(6):873–81.
Mills EL, Kelly B, O’Neill LAJ. Mitochondria are the powerhouses of immunity. Nat Immunol. 2017;18(5):488–98.
Minhas PS, Latif-Hernandez A, McReynolds MR, Durairaj AS, Wang Q, Rubin A, et al. Restoring metabolism of myeloid cells reverses cognitive decline in ageing. Nature. 2021;590(7844):122–8.
Ridker PM, Everett BM, Thuren T, MacFadyen JG, Chang WH, Ballantyne C, et al. Antiinflammatory therapy with canakinumab for atherosclerotic disease. N Engl J Med. 2017;377(12):1119–31.
Ridker PM, MacFadyen JG, Everett BM, Libby P, Thuren T, Glynn RJ. Relationship of C-reactive protein reduction to cardiovascular event reduction following treatment with canakinumab: a secondary analysis from the CANTOS randomised controlled trial. Lancet. 2018;391(10118):319–28.
Hsu CC, Chang HY, Wu IC, Chen CC, Tsai HJ, Chiu YF, et al. Cohort Profile: The Healthy Aging Longitudinal Study in Taiwan (HALST). Int J Epidemiol. 2017;46(4):1106–1106j.
Ridker PM, Thuren T, Zalewski A, Libby P. Interleukin-1β inhibition and the prevention of recurrent cardiovascular events: rationale and design of the Canakinumab Anti-inflammatory Thrombosis Outcomes Study (CANTOS). Am Heart J. 2011;162(4):597–605.
Arnett DK, Blumenthal RS, Albert MA, Buroker AB, Goldberger ZD, Hahn EJ, et al. 2019 ACC/AHA Guideline on the Primary Prevention of Cardiovascular Disease: A Report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines. Circulation. 2019;140(11):e596–646.
Liu CS, Cheng WL, Lee CF, Ma YS, Lin CY, Huang CC, et al. Alteration in the copy number of mitochondrial DNA in leukocytes of patients with mitochondrial encephalomyopathies. Acta Neurol Scand. 2006;113(5):334–41.
Rizopoulos D. Joint Models for Longitudinal and Time-to-Event Data: With Applications in R. 1st ed. Boca Raton, FL: Chapman and Hall/CRC; 2012.
Asar Ö, Ritchie J, Kalra PA, Diggle PJ. Joint modelling of repeated measurement and time-to-event data: an introductory tutorial. Int J Epidemiol. 2015;44(1):334–44.
Greenland S, Pearl J, Robins JM. Causal diagrams for epidemiologic research. Epidemiology. 1999;10(1):37–48.
Textor J, van der Zander B, Gilthorpe MS, Liskiewicz M, Ellison GT. Robust causal inference using directed acyclic graphs: the R package “dagitty.” Int J Epidemiol. 2016;45(6):1887–94.
Pearson TA, Mensah GA, Alexander RW, Anderson JL, Cannon RO 3rd, Criqui M, et al. Markers of inflammation and cardiovascular disease: application to clinical and public health practice: a statement for healthcare professionals from the Centers for Disease Control and Prevention and the American Heart Association. Circulation. 2003;107(3):499–511.
VanderWeele TJ, Ding P. Sensitivity Analysis in Observational Research: Introducing the E-Value. Ann Intern Med. 2017;167(4):268–74.
Nissen SE, Tuzcu EM, Schoenhagen P, Crowe T, Sasiela WJ, Tsai J, et al. Statin therapy, LDL cholesterol, C-reactive protein, and coronary artery disease. N Engl J Med. 2005;352(1):29–38.
Ridker PM, Cannon CP, Morrow D, Rifai N, Rose LM, McCabe CH, et al. C-reactive protein levels and outcomes after statin therapy. N Engl J Med. 2005;352(1):20–8.
Heinrich PC, Castell JV, Andus T. Interleukin-6 and the acute phase response. Biochem J. 1990;265(3):621–36.
Pepys MB, Baltz ML. Acute phase proteins with special reference to C-reactive protein and related proteins (pentaxins) and serum amyloid A protein. Adv Immunol. 1983;34:141–212.
Vigushin DM, Pepys MB, Hawkins PN. Metabolic and scintigraphic studies of radioiodinated human C-reactive protein in health and disease. J Clin Invest. 1993;91(4):1351–7.
Wu IC, Lin CC, Liu CS, Hsu CC, Chen CY, Hsiung CA. Interrelations between mitochondrial DNA copy number and inflammation in older adults. J Gerontol A Biol Sci Med Sci. 2017;72(7):937–44.
Fazzini F, Lamina C, Fendt L, Schultheiss UT, Kotsis F, Hicks AA, 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.
We thank all participating men and women in the HALST study; and all members of the HALST study group.
This work was supported by the National Health Research Institutes (Project nos. BS-097-SP-04, PH-098-SP-02, PH-099-SP-01, PH-100-SP-01, PH-101-SP-01, PH-102-SP-01, PH-103-SP-01, PH-104-SP-01, PH-105-SP-01, PH-106-SP-01, PH-107-SP-01, PH-107-PP-22, PH-108-SP-01, PH-108-PP-22, PH-109-SP-01, PH-109- PP-22, PH-110-SP-01, PH-110-PP-17, PH-111-SP-01, PH-111-PP-15) and Ministry of Science and Technology (Project no. MOST 109–2314-B-400–026) in Taiwan. The sponsors had no roles in the design and conduct of the study; collections and analysis of the data; preparation and approval of the manuscript.
Ethics approval and consent to participate
This study was approved by the institutional review boards of the National Health Research Institutes and participating hospitals. All participants provided written informed consent according to the Declaration of Helsinki.
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.
Additional file 1:
Supplementary Figure 1. The assumed relationship among leukocyte mitochondrial DNA copy number, serum high-sensitive C-reactive protein, risk of death. Supplementary Methods. Measurement of Covariate Variables, Statistical Analyses, References. Supplementary Table 1. Characteristics of HALST Participants Excluded from This Study. Supplementary Table 2. Association of Baseline Leukocyte Mitochondrial DNA Copy Number with High-Sensitivity C-Reactive Protein Trajectory After Excluding Participants with High-Sensitivity C-Reactive Protein >10 mg/L During Follow-Up. Supplementary Table 3. Association of Baseline Leukocyte Mitochondrial DNA Copy Number with High-Sensitivity C-Reactive Protein Trajectory After Excluding Participants with High-Sensitivity C-Reactive Protein >3 mg/L During Follow-Up. Supplementary Table 4. Association between Estimated High-Sensitivity C-Reactive Protein Trajectory and Risk of Mortality. Supplementary Table 5. E-Values for the Association of Features of Estimated High-Sensitivity C-Reactive Protein Trajectory Over Time With Risk of Mortality in All Participants. Supplementary Table 6. E-Values for the Association of Features of Estimated High-Sensitivity C-Reactive Protein Trajectory Over Time With Risk of Mortality after the Exclusion of Participants With High-Sensitivity C-Reactive Protein > 10 mg/L During Follow-Up. Supplementary Table 7. E-Values for the Association of Features of Estimated High-Sensitivity C-Reactive Protein Trajectory Over Time With Risk of Mortality after the Exclusion of Participants With High-Sensitivity C-Reactive Protein > 3 mg/L During Follow-Up.
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.
About this article
Cite this article
Wu, IC., Liu, CS., Cheng, WL. et al. Association of leukocyte mitochondrial DNA copy number with longitudinal C-reactive protein levels and survival in older adults: a cohort study. Immun Ageing 19, 62 (2022). https://doi.org/10.1186/s12979-022-00322-8