Abstract

Cellular aging markers, including telomere length and mitochondrial function, as well as oxidative stress and inflammation markers influence each other and form a complex network, which is affected in diabetes. However, it remains unknown whether these markers could independently predict future diabetes after adjustment for their mutual effects. We conducted a 3-year longitudinal study in a Chinese cohort that comprised 108 nondiabetic individuals at baseline. The 2-hour 75 g oral glucose tolerance tests were performed at baseline and at 3-year follow-up. At baseline, leukocyte telomere length (LTL) and mitochondrial DNA copy number (mtDNAcn) in leukocytes were determined using the polymerase chain reaction method. Tumor necrosis factor (TNF-α), interleukin-6, 8-hydroxy-2-deoxyguanosine levels, and superoxide dismutase (SOD) activity were measured by the enzyme-linked immunosorbent assay. Participants who developed diabetes at the 3-year follow-up () had shorter LTL and higher levels of TNF-α and SOD activity at baseline. Baseline LTL was found to be independently associated with the development of diabetes at the 3-year follow-up after the adjustment for mtDNAcn, markers of oxidative stress and inflammation, and conventional diabetes risk factors. Our findings suggest that LTL is an independent predictor for 3-year diabetes risk, which might inform timely prevention and treatment of diabetes. Telomere shortening might be involved in the pathogenesis of diabetes independently of conventional diabetes risk factors, mtDNAcn, or oxidative stress and inflammation pathways.

1. Introduction

Over the last several decades, numerous studies have revealed that telomere length, mitochondrial function, oxidative stress, and inflammation are cellular aging markers intertwined in the complex relationships that determine their roles in cellular senescence [13]. These cellular aging markers have also been reported to be associated with various age-related diseases [4, 5], including diabetes [6, 7]. The relationships between these cellular aging markers and diabetes are particularly complex. The potential associations between leukocyte telomere length (LTL), mitochondrial DNA copy number (mtDNAcn), markers of oxidative stress/antioxidants, and proinflammatory cytokines on the one hand and diabetes on the other hand have been intensively scrutinized. Unfortunately, the cross-sectional design of the majority of the previous population-based studies failed to uncover causal relationships between these cellular aging markers and diabetes, with many findings being inconsistent [814]. Prospective studies are scarce, with only three published reports [1517] investigating the prospective association of LTL with the development of diabetes. Furthermore, the results of those studies were discrepant and confusing. In a multiracial cohort of postmenopausal women, just a weak association was found between baseline LTL and future diabetes risk, which was attenuated after adjusting for conventional diabetes risk factors [15]. In contrast, in an American Indian cohort, baseline LTL was identified as a predictive marker independent of conventional diabetes risk factors [16]. Similar results were obtained in a European cohort [17]. The telomere length is characterized by strong heritability, and it varies substantially among different races [18], which could partly explain the discrepancy among these studies. Notably, none of these prospective studies included East Asian or, more specifically, Chinese subjects. Additionally, studies exploring the predictive role of mtDNAcn as well as markers of oxidative stress and inflammation in diabetes development are rare in general. Furthermore, considering the complex relationships between these biomarkers, it should be established whether they could affect each other’s predictive role in diabetes risk. Taking into account these considerations, we conducted a 3-year longitudinal study in a Chinese population to investigate the prospective associations of baseline LTL and mtDNAcn, as well as oxidative stress and inflammatory markers with the development of diabetes 3 years later.

2. Methods

2.1. Study Participants

The current study was conducted in a Chinese cohort from the ChangPing suburb of Beijing between March 2014 and May 2017. The study protocols and consent procedures had been approved by the Ethics Committee of the Peking Union Medical College Hospital. All participants have signed the informed consent forms. A total of 142 individuals who had completed data in 2014 and 2017 were recruited into this study. According to the World Health Organization criteria [19], participants with and 2-hour postload during the 75 g oral glucose tolerance test (OGTT) were identified as nondiabetic individuals, including those with normal glucose tolerance (NGT) and prediabetes (impaired fasting glucose ( and ) and/or impaired glucose tolerance ( and )). There were 34 individuals who had diabetes at baseline and were thus excluded from the study. Hence, a total of 108 nondiabetic individuals at baseline were selected as the study population, including those with NGT () and prediabetes () at baseline. Demographic and anthropometric data were collected at baseline and at 3-year follow-up, including age, sex, height, weight, waist circumference (WC), and hip circumference (HC).

2.2. Biochemical Analysis

All the participants underwent a 75 g 2-hour OGTT, in which blood samples were drawn at four time points at baseline (0 min, 30 min, 60 min, and 120 min), whereas blood samples were obtained at two time points at the 3-year follow-up (0 min and 120 min). The levels of PG were determined by the glucose oxidase assay. Serum insulin (INS) and C-peptide (CP) were determined by the chemiluminescence immunoassay (Insulin IRI #02230141 (128434) and C-peptide #03649928 (129026), both from Siemens Medical Solutions Diagnostics, Tarrytown, NY, USA). Glycosylated hemoglobin (HbA1c) was measured by high-performance liquid chromatography, whereas total cholesterol (TC), total triglyceride (TG), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), uric acid (UA), creatinine (Cr), and urea were measured using an automatic analyzer.

2.3. Calculations of Insulin Sensitivity and β-Cell Function Indices

Based on the results of the 75 g 2-hour OGTT, we evaluated the sensitivity to insulin and β-cell function of the participants by using the following equations [2022]:

2.4. LTL Assays

The LTL assays were performed based on the blood samples collected at baseline. The details of LTL measurement have been described in our previous publication [23]. LTL was determined as the relative ratio of the telomere repeat copy number to the single copy number (T/S ratio) using the monochrome multiplex quantitative polymerase chain reaction protocol [24].

2.5. Measurement of mtDNAcn in Peripheral Blood

mtDNAcn was measured based on the blood samples collected at baseline. As described previously [25], genomic DNA was extracted from leukocytes using a QIAamp DNA Blood Midi Kit (Qiagen, Hilden, Germany) and subsequently purified, diluted, and quantified using a NanoDrop 1000 Spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA). The relative mtDNAcn value was determined using the real-time polymerase chain reaction and adjusted by simultaneously measuring the amount of nuclear DNA.

2.6. Assessment of Oxidative Stress and Inflammatory Markers

The measurements of oxidative stress and inflammatory markers were performed based on the blood samples collected at baseline. The measurements of oxidative stress and inflammatory markers have also been described in detail previously [23]. Superoxide dismutase (SOD) activity, as well as the levels of 8-hydroxy-2-deoxyguanosine (8-OHdG), tumor necrosis factor (TNF-α), and interleukine-6 (IL-6) were determined according to the manufacturer’s instructions (Cloud-Clone Corp., Houston, USA). Absorbance kinetics were measured using an enzyme-linked immunosorbent assay reader.

2.7. Statistical Analysis

Statistical analysis was conducted using the SPSS 25.0 software package (IBM), R Studio (version 1.2.1335, https://www.r-project.org/), and GraphPad Prism 8.0 (https://www.graphpad.com). Continuous variables are presented as the , whereas categorical variables are presented as percentages. To examine if the data were normally distributed, the Shapiro-Wilk and Kolmogorov-Smirnov tests were used. Normality transformations were performed on the variables that did not meet the normality assumptions using the appropriate formula when necessary. Comparisons of the continuous variables between groups were performed by Student’s -test or nonparametric Mann-Whitney’s -test, where appropriate. Pearson’s or Spearman’s correlation analysis was performed to determine the bivariate correlations, where appropriate. To explore whether mtDNAcn mediated the effect of TNF-α on 30-minute postload plasma glucose (PG30min), PROCESS macro Version 3.3 was used to generate simple mediation models with ordinary least squares. Mediation hypotheses were tested via a bias-corrected bootstrap method with 5,000 samples to calculate 95% confidence intervals (95% CI). If the 95% CI did not encompass zero, this represented statistical significance of the mediating effect. To identify the independent predictors for diabetes risk, multivariate logistic regression analysis was performed. In all comparisons, the statistical significance was set at .

3. Results

3.1. Relationships between LTL, TNF-α Level, and SOD Activity and the Risk of Diabetes Development

According to whether the participants developed diabetes at the 3-year follow-up, they were divided into progressors () and nonprogressors (). Tables 1 and 2 present the differences in characteristics between progressors and nonprogressors at baseline and at the 3-year follow-up, respectively. The proportions of NGT subjects in progressors and nonprogressors at baseline were similar (51.7% versus 50.6%, ). At baseline, the progressors had a significantly higher level of PG30min () and PG60min () than nonprogressors, with slightly higher HbA1c (). However, the differences of other conventional risk factors of diabetes at baseline, including age, BMI, WC, lipid profile, UA, insulin sensitivity indices, and insulin secretion indices between the two groups did not reach statistical significance (). Baseline age-adjusted LTL of progressors was significantly shorter than that of nonprogressors (). TNF-α levels () and SOD activity () of progressors were higher than those of nonprogressors, whereas we failed to observe statistically significant differences in the values of mtDNAcn (), IL-6 (), and 8-OHdG () between groups. As expected, at the 3-year follow-up, HbA1c, FPG, PG120min, fasting insulin (FINS), CP120min, TG, and HOMA-IR values of progressors were significantly higher than those of nonprogressors.

3.2. Correlation of Cellular Aging Markers with Clinical Characteristics at Baseline

As indicated in Table 3, we found a modest positive correlation between baseline LTL and the Matsuda index (, ), which proved to be a favorable OGTT-stimulated insulin sensitivity index. Nevertheless, we did not find an association between baseline LTL and any other conventional diabetes risk factor at baseline. mtDNAcn negatively correlated with PG30min, TG, and UA. The level of TNF-α positively correlated with HbA1c, PG30min, and UA, and negatively correlated with the Matsuda index. 8-OHdG level correlated with age, HDL-C, HOMA-IR, and HOMA-β. As for the correlations among the cellular aging markers per se, both LTL and mtDNAcn correlated with TNF-α level and SOD activity. To explore whether mtDNAcn mediated the effect of TNF-α on PG30min, mediation model analysis was performed. As shown in Figure 1, mtDNAcn fully mediated the effect of TNF-α on PG30min (direct effect ; 95% CI -0.0260–0.0571; indirect effect ; 95% CI 0.0077–0.0635).

3.3. Prospective Correlations of Baseline LTL, TNF-α, 8-OHdG, and SOD Activity with PG120min at the 3-Year Follow-Up

Subsequently, we performed Spearman’s correlation analysis to explore the associations between parameters of cellular aging markers at baseline and clinical characteristics at the 3-year follow-up. As indicated in Table 4, baseline LTL negatively correlated whereas TNF-α level and SOD activity positively correlated with PG120min at the 3-year follow-up (, 0.324, and 0.272, respectively, ). Baseline 8-OHdG level inversely correlated with PG120min and HOMA-β at the 3-year follow-up (, , and , , respectively). Additionally, both baseline LTL and TNF-α correlated with TC at the 3-year follow-up ( and 0.203, respectively, ).

3.4. Independent and Prospective Association between LTL and Future Diabetes Risk

Multivariate logistic regression analysis was performed to explore the independent predictors for the 3-year diabetes risk. As shown in Table 5, age and sex-adjusted LTL at baseline were significantly associated with the development of diabetes at the 3-year follow-up (, 95% CI 0.231–0.776, adjusted , Model 1). Following further and progressive adjustments for conventional diabetes risk factors (Model 2), as well as oxidative stress and inflammatory markers at baseline (Model 3), the association remained significant (Model 2: , 95% CI 0.187–0.919, adjusted ; Model 3: , 95% CI 0.139–0.839, adjusted ). As shown in Table S1, multivariate logistic regression analysis according to LTL quintiles demonstrated that individuals with LTL values in the second quintile were more than twelvefold likely to develop diabetes as those with the longest LTL values in the highest quintile (, 95% CI 1.240–124.252, adjusted ). However, the elevation of diabetes risk in individuals with LTL values in the first (, 95% CI 0.646–115.257, adjusted ), third (, 95% CI 0.155–16.684, adjusted ), and fourth quintiles (, 95% CI 0.315–33.411, adjusted ) compared to those in individuals with the longest LTL values in the highest quintile did not reach statistical significance, suggesting a nonlinear association between LTL at baseline and diabetes risk at the 3-year follow-up.

4. Discussion

Our longitudinal findings indicate that short LTL at baseline may be predictive of increased future diabetes risk. This relationship was independent of the effects of conventional diabetes risk factors, mtDNAcn, or oxidative stress and inflammatory markers. Therefore, telomere shortening is possibly involved in the pathogenesis of diabetes via the pathways other than those mediating the effects of conventional diabetes risk factors, mitochondrial dysfunction, oxidative stress, and inflammation.

Our observation that LTL predicts diabetes risk independently of conventional diabetes risk factors (e.g., age, BMI, WC, FPG, and lipid profile) was consistent with findings of previous prospective studies [16, 17]. We considerably extended the conclusions of those reports by showing that the correlation between LTL for diabetes risk is also independent of mtDNAcn and markers of oxidative stress and inflammation. Thus far, the role of the latter biomarkers in the development of diabetes remained confusing. Previous data indicated that telomere shortening and mitochondrial dysfunction likely affect each other and thus form a self-amplifying cycle [26]. The impaired mitochondria produce excessive reactive oxygen species that subsequently cause oxidative stress and inflammation [27]. In turn, oxidative stress and inflammation have also been demonstrated to damage telomere [2830] and mtDNA [31, 32], which results in a vicious cycle. These changes may contribute to the development of diabetes via their possible involvement in insulin resistance [31, 33, 34] and β-cell dysfunction [35, 36]. Our finding of the positive correlation between LTL and the insulin sensitivity index, the Matsuda index, at baseline indicates an association between telomere shortening and resistance to insulin. The prospective inverse correlation between baseline 8-OHdG and HOMA-β at the 3-year follow-up suggests that oxidative DNA damage might participate in β-cell dysfunction, which is consistent with a previous finding of pancreatic β-cells being particularly vulnerable to oxidative damage [37]. The levels of 8-OHdG, a well-known biomarker of DNA oxidation [38], were demonstrated to be elevated in prediabetes and in patients with T2D in numerous previous retrospective studies [39, 40]. Unexpectedly, we found a modest inverse correlation between baseline 8-OHdG and PG120min at the 3-year follow-up, suggesting that elevated 8-OHdG might predict or contribute to decreased postload plasma glucose, which seemingly contradicts the previous findings. Nevertheless, considering the retrospective design of the previous studies and the lack of studies exploring the prospective association between 8-OHdG and future hyperglycemia, further large-scale prospective investigations are needed to validate our observations. Despite the important role of mitochondrial dysfunction, oxidative stress, and inflammation in the development of diabetes, the prospective association between LTL and the development of diabetes, which in our study was independent of mtDNAcn or markers of oxidative stress and inflammation, indicates that other potential pathways might mediate the effects of telomere shortening on the predisposition to diabetes. Further investigations are necessary to uncover mechanisms and pathways underlying this association.

To the best of our knowledge, this is the first study to explore an independent effect of mtDNAcn on the risk of future diabetes. We did not detect an independent effect of this parameter on the risk of diabetes at the 3-year follow-up. Another prospective study observed that adding mtDNAcn as an adjunctive marker augmented the predictive power of HbA1c and OGTT for future diabetes, but its independent role was not investigated [41]. Considering the small sample size of our present study, larger population studies are needed to validate our conclusion.

Several previous cross-sectional and prospective studies reported associations between LTL and the fasting insulin sensitivity index HOMA-IR [33, 42, 43]. In our correlation analysis, we found that LTL positively correlated with the Matsuda index at baseline, but not with HOMA-IR either at baseline or at the 3-year follow-up, which conflicts with the previous reports. Although the small sample size of the present study might partially explain this discrepancy, it is worth noting that these two indices provide different information: the Matsuda index reflects whole-body insulin sensitivity, whereas HOMA-IR has been recognized as an index of peripheral or hepatic insulin sensitivities based on the assumption that they are equivalent [21]. In addition, the Matsuda index has been reported to highly correlate with the gold standard whole-body insulin sensitivity index derived from euglycemic insulin clamp, suggesting that the Matsuda index provides a reasonable approximation of the whole-body insulin sensitivity [21]. Unfortunately, during the OGTT at the 3-year follow-up, PG30min and PG60min were not obtained; therefore, the value of the Matsuda index could not be calculated. The prospective association of LTL and the Matsuda index needs to be explored in future studies.

Interestingly, we found that both mtDNAcn and TNF-α correlated with PG30min. The mediation model analysis revealed that TNF-α indirectly contributed to the elevation of PG30min via the decrease of mtDNAcn, whereas its direct effect was not significant. Our findings indicated that TNF-α likely contributes to hyperglycemia predominately via mitochondrial dysfunction rather than other pathways, but this remains to be confirmed. Notably, except for PG30min, parameters such as FPG, PG60min, and PG120min did not correlate with mtDNAcn and TNF-α. After glucose loading, the first-phase insulin secretion occurs within 30 min, suppressing glucose production in the liver [44]. Impaired first-phase insulin secretion leads to the elevation of PG30min. Therefore, our result suggests that mtDNAcn and TNF-α are possibly associated with impaired first-phase insulin secretion.

Overall, we made several important observations in the present study. First, to the best of our knowledge, our analysis for the first time demonstrated that LTL independently predicted 3-year diabetes risk after adjusting for the confounding effects of mtDNAcn, oxidative stress, and inflammation marker levels, as well as of the conventional diabetes risk factors. Second, we investigated the associations between the cellular aging markers on the one hand and PG30min and PG60min values on the other hand, which were rarely determined in previous studies. Intriguingly, this is the first study that reports correlations of mtDNAcn and TNF-α level with PG30min, although the exact mechanisms underlying these relationships need to be clarified in future studies.

However, several limitations also need to be mentioned. The relatively small sample size limited statistical power and thereby likely precluded confirmation of the associations between some parameters that could be revealed in a larger cohort. This circumstance partially explains the absence of correlations among the indicators. On the other hand, the follow-up time was relatively short; thus, the earlier predictive value of these cellular aging markers was not explored in the current study.

5. Conclusion

Collectively, our observations suggest that LTL is an independent predictor for the 3-year diabetes risk. Its predictive role, which was independent of mtDNAcn and markers of oxidative stress and inflammation, as well as of conventional diabetes risk factors, might provide novel insights into the potential pathways that mediate the effect of telomere shortening on diabetes development. These specific pathways will need to be identified in future studies. Importantly, identification of telomere length as an additional risk modifier could assist early prediction for future diabetes and its timely prevention. However, large population-based studies are needed to validate the clinical value of our present findings.

Data Availability

The SPSS Statistics data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the CAMS Innovation Fund for Medical Sciences (CIFMS) (grant number CIFMS2016-I2M-4-001) and the Nonprofit Central Research Institute Fund of the Chinese Academy of Medical Sciences (grant number No. 2019PT320007). We would like to thank Editage (https://www.editage.cn) for English language editing.

Supplementary Materials

Table S1: multivariate-adjusted OR (95% CI) for diabetes risk and oxidative stress and inflammation according to baseline LTL quintiles. (Supplementary Materials)