Skip to main content

Association of leukocyte mitochondrial DNA copy number with longitudinal C-reactive protein levels and survival in older adults: a cohort study

Abstract

Background

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.

Methods

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.

Results

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.

Conclusions

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.

Introduction

Chronically high blood levels of C-reactive protein (CRP), an acute-phase protein induced by the systemic effects of proinflammatory cytokines [1], 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 [9]. 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. [19] showed that one standard deviation lower leukocytes mtDNACN was associated with greater risks of incident coronary heart disease and stroke. Tin et al. [20] 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 [1]. 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 [2], 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.

Methods

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 [27]. 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 [30]. 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).

Statistical analyses

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.

Fig. 1
figure 1

Directed acyclic graph depicting the assumed relationship among leukocyte mitochondrial DNA copy number, serum high-sensitive C-reactive protein, risk of death, and the potential confounders. Panel A displayed the relationships among leukocyte mitochondrial DNA copy number (exposure, represented by the green circle), serum high-sensitive C-reactive protein (outcome, represented by the blue circle), and covariates associated with both the leukocyte mitochondrial DNA copy number and serum high-sensitive C-reactive protein (ancestors of both the exposure and outcome, represented by the red circles). Panel B displayed the relationships among serum high-sensitive C-reactive protein (exposure, represented by the green circle), risk of death (outcome, represented by the blue circle), and the covariates associated with both serum high-sensitive C-reactive protein and risk of death (ancestors of both the exposure and outcome, represented by the red circles). The arrow from variable X to variable Y represents a direct effect of X on Y (an effect not mediated by other observed variables in the graph) that may have been present. The hypothesized causal pathway is presented in green. The biasing pathway is presented in red [34]

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 [35]. 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 [36]. 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).

Results

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).

Table 1 Baseline characteristics of study participantsa

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.

Table 2 Association of baseline leukocyte mitochondrial DNA copy number with high-sensitivity C-reactive protein trajectorya

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.

Table 3 Association between estimated high-sensitivity C-reactive protein trajectory and risk of mortalitya

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.

Discussion

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 [41]. 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 [42]. In a cohort of 9058 participants without chronic kidney disease [20], 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 [43], 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 [13]. Second, as the hub of multiple metabolic and signaling pathways, mitochondria in immune cells play pivotal roles in shaping immune and inflammatory responses [23]. 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 [24].

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.

Conclusions

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.

Abbreviations

BMI:

Body mass index

CRP:

C-reactive protein

DAG:

Directed acyclic graph

HALST:

Healthy Aging Longitudinal Study in Taiwan

HDL cholesterol:

High-density lipoprotein cholesterol

hs-CRP:

High-sensitivity CRP

LDL cholesterol:

Low-density lipoprotein cholesterol

mtDNACN :

Mitochondrial DNA copy number

References

  1. Pepys MB, Hirschfield GM. C-reactive protein: a critical update. J Clin Invest. 2003;111(12):1805–12.

    Article  CAS  Google Scholar 

  2. Singh T, Newman AB. Inflammatory markers in population studies of aging. Ageing Res Rev. 2011;10(3):319–29.

    Article  CAS  Google Scholar 

  3. Ferrucci L, Fabbri E. Inflammageing: chronic inflammation in ageing, cardiovascular disease, and frailty. Nat Rev Cardiol. 2018;15(9):505–22.

    Article  CAS  Google Scholar 

  4. 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.

    Article  CAS  Google Scholar 

  5. 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.

    Article  Google Scholar 

  6. Newman AB, Murabito JM. The epidemiology of longevity and exceptional survival. Epidemiol Rev. 2013;35(1):181–97.

    Article  Google Scholar 

  7. Franceschi C, Garagnani P, Vitale G, Capri M, Salvioli S. Inflammaging and “Garb-aging.” Trends Endocrinol Metab. 2017;28(3):199–212.

    Article  CAS  Google Scholar 

  8. Tchkonia T, Kirkland JL. Aging, cell senescence, and chronic disease: emerging therapeutic strategies. JAMA. 2018;320(13):1319–20.

    Article  Google Scholar 

  9. Lopez-Otin C, Blasco MA, Partridge L, Serrano M, Kroemer G. The hallmarks of aging. Cell. 2013;153(6):1194–217.

    Article  CAS  Google Scholar 

  10. 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.

    Article  CAS  Google Scholar 

  11. Malik AN, Czajka A. Is mitochondrial DNA content a potential biomarker of mitochondrial dysfunction? Mitochondrion. 2013;13(5):481–92.

    Article  CAS  Google Scholar 

  12. 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.

    Article  CAS  Google Scholar 

  13. 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.

    Article  Google Scholar 

  14. Khan SS, Singer BD, Vaughan DE. Molecular and physiological manifestations and measurement of aging in humans. Aging Cell. 2017;16(4):624–33.

    Article  CAS  Google Scholar 

  15. 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.

    Article  CAS  Google Scholar 

  16. 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.

    Article  CAS  Google Scholar 

  17. 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.

    Article  CAS  Google Scholar 

  18. 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.

    Article  Google Scholar 

  19. 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.

    Article  Google Scholar 

  20. 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.

    Article  CAS  Google Scholar 

  21. 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.

    Article  CAS  Google Scholar 

  22. 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.

    Article  CAS  Google Scholar 

  23. Mills EL, Kelly B, O’Neill LAJ. Mitochondria are the powerhouses of immunity. Nat Immunol. 2017;18(5):488–98.

    Article  CAS  Google Scholar 

  24. 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.

    Article  CAS  Google Scholar 

  25. 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.

    Article  CAS  Google Scholar 

  26. 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.

    Article  CAS  Google Scholar 

  27. 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.

    Article  Google Scholar 

  28. 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.

    Article  CAS  Google Scholar 

  29. 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.

    Google Scholar 

  30. 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.

    CAS  Google Scholar 

  31. 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.

    Book  Google Scholar 

  32. 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.

    Article  Google Scholar 

  33. Greenland S, Pearl J, Robins JM. Causal diagrams for epidemiologic research. Epidemiology. 1999;10(1):37–48.

    Article  CAS  Google Scholar 

  34. 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.

    Google Scholar 

  35. 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.

    Article  Google Scholar 

  36. VanderWeele TJ, Ding P. Sensitivity Analysis in Observational Research: Introducing the E-Value. Ann Intern Med. 2017;167(4):268–74.

    Article  Google Scholar 

  37. 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.

    Article  CAS  Google Scholar 

  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.

    Article  CAS  Google Scholar 

  39. Heinrich PC, Castell JV, Andus T. Interleukin-6 and the acute phase response. Biochem J. 1990;265(3):621–36.

    Article  CAS  Google Scholar 

  40. 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.

    Article  CAS  Google Scholar 

  41. 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.

    Article  CAS  Google Scholar 

  42. 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.

    Article  CAS  Google Scholar 

  43. 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.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We thank all participating men and women in the HALST study; and all members of the HALST study group.

Funding

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.

Author information

Authors and Affiliations

Authors

Contributions

I.C.W. had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. I.C.W. conceived and designed the study. I.C.W., C.S.L., W.L.C., T.T.L., H.L.C., R.C.W., C.A.H., C.C.H. contributed to the data acquisition. I.C.W., H.L.C., P.F.C., R.C.W., C.W.H. contributed to the data curation. I.C.W., P.F.C., C.W.H. contributed to the data analysis. I.C.W., C.S.L, C.A.H., C.C.H. contributed to data interpretation. I.C.W. drafted the manuscript. All authors contributed to critical revision of the manuscript for important intellectual content., C.A.H., C.C.H., I.C.W. obtained funding. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to I-Chien Wu.

Ethics declarations

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

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12979-022-00322-8

Keywords