Skip to main content

DNA methylation age calculators reveal association with diabetic neuropathy in type 1 diabetes

Abstract

Background

Many CpGs become hyper or hypo-methylated with age. Multiple methods have been developed by Horvath et al. to estimate DNA methylation (DNAm) age including Pan-tissue, Skin & Blood, PhenoAge, and GrimAge. Pan-tissue and Skin & Blood try to estimate chronological age in the normal population whereas PhenoAge and GrimAge use surrogate markers associated with mortality to estimate biological age and its departure from chronological age. Here, we applied Horvath’s four methods to calculate and compare DNAm age in 499 subjects with type 1 diabetes (T1D) from the Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications (DCCT/EDIC) study using DNAm data measured by Illumina EPIC array in the whole blood. Association of the four DNAm ages with development of diabetic complications including cardiovascular diseases (CVD), nephropathy, retinopathy, and neuropathy, and their risk factors were investigated.

Results

Pan-tissue and GrimAge were higher whereas Skin & Blood and PhenoAge were lower than chronological age (p < 0.0001). DNAm age was not associated with the risk of CVD or retinopathy over 18–20 years after DNAm measurement. However, higher PhenoAge (β = 0.023, p = 0.007) and GrimAge (β = 0.029, p = 0.002) were associated with higher albumin excretion rate (AER), an indicator of diabetic renal disease, measured over time. GrimAge was also associated with development of both diabetic peripheral neuropathy (OR = 1.07, p = 9.24E−3) and cardiovascular autonomic neuropathy (OR = 1.06, p = 0.011). Both HbA1c (β = 0.38, p = 0.026) and T1D duration (β = 0.01, p = 0.043) were associated with higher PhenoAge. Employment (β = − 1.99, p = 0.045) and leisure time (β = − 0.81, p = 0.022) physical activity were associated with lower Pan-tissue and Skin & Blood, respectively. BMI (β = 0.09, p = 0.048) and current smoking (β = 7.13, p = 9.03E−50) were positively associated with Skin & Blood and GrimAge, respectively. Blood pressure, lipid levels, pulse rate, and alcohol consumption were not associated with DNAm age regardless of the method used.

Conclusions

Various methods of measuring DNAm age are sub-optimal in detecting people at higher risk of developing diabetic complications although some work better than the others.

Background

CpGs are regions of DNA where a cytosine is followed by a guanine nucleotide. Cytosines within CpGs can be methylated, and CpG methylation levels affect gene expression. Many CpGs become hyper or hypo-methylated with age [1,2,3,4]. In 2013, Horvath used publicly available DNA methylation (DNAm) data to define and evaluate a DNAm age predictor, Pan-tissue, which is accurate across most tissues and cell types. Chronological age was regressed on CpG methylation levels using a penalized regression model (elastic net) which selected 353 CpGs [1] (Supplementary Table 1). Pan-tissue has been widely used and has shown that faster epigenetic aging is associated with multiple age-related diseases and conditions (e.g., Alzheimer’s, cancer, cardiovascular diseases (CVD)) indicating that epigenetic age is an indicator of health status. Some risk factors for type 2 diabetes (T2D) including BMI, waist circumference, and fasting glucose have been associated with higher Pan-tissue epigenetic age acceleration (EAA = epigenetic age − chronological age) [5, 6]. On average liver Pan-tissue EAA increased significantly by 0.33 years per BMI unit [7].

Pan-tissue performed sub-optimally estimating fibroblast age in in vitro studies. Therefore, Horvath et al. described Skin & Blood DNAm age using a similar method to Pan-tissue which performed remarkably well across a wide spectrum of cells that are most frequently used in in vitro studies, including the blood (Supplementary Table 1). Skin & Blood EAA was highly predictive of time to all-cause mortality. It was positively correlated with waist/hip ratio (WHR), blood insulin, glucose, triglyceride, systolic blood pressure, and BMI and negatively correlated with HDL and physical exercise. However, the respective correlation coefficients were weak (|r| < 0.11) [2].

Pan-tissue and Skin & Blood were both developed using chronological age as a surrogate for biological age. Therefore, they may not capture CpGs that signal departure of biological age from chronological age. In a newer method called PhenoAge, chronological age was replaced with a surrogate measure of phenotypic age developed using clinical data. A Cox penalized regression model was applied where the hazard of mortality was regressed on 42 clinical markers and chronological age to select variables for inclusion in the phenotypic age score. Nine clinical markers and chronological age were selected and used to estimate the 10-year mortality risk score which was then converted into units of years. Finally, the resulting phenotypic age estimate was regressed on CpG methylation levels using an elastic net regression model (Supplementary Table 1). A 1-year increase in PhenoAge was associated with a 4.5% increase in the risk of all-cause mortality in independent populations without diabetes. PhenoAge predicted mortality significantly better than Pan-tissue. It was also associated with increased risk of CVD and differed significantly between never, current, and former smokers. It was also positively correlated with blood insulin, glucose, triglyceride, and WHR and negatively correlated with HDL and physical exercise [8].

Most recently, Horvath et al. defined another DNAm age calculator called GrimAge. To develop this method, 88 plasma proteins and smoking pack-years were individually regressed on chronological age, sex, and the CpGs methylation levels using an elastic net regression model. Twelve plasma proteins and pack-years had high correlations between their DNAm estimation and the corresponding measured levels. These twelve DNAm estimated plasma proteins and smoking pack-years as well as chronological age and sex were regressed on the hazard of aging-related mortality using a Cox penalized regression model. This selected DNAm pack-years, age, sex, and the predicted DNAm for 7 of the 12 plasma proteins. The combined estimate of these factors was then transformed into GrimAge which has the same mean and variance as chronological age (Supplementary Table 1). GrimAge was highly predictive of lifespan and time-to-CVD even after adjustment for known risk factors and outperformed Pan-tissue and PhenoAge. It was also associated with hypertension and T2D. In addition, GrimAge was correlated with BMI, WHR, and physical exercise; blood insulin, glucose, HbA1c, triglyceride, and HDL; and albuminuria (all correlations were in the expected directions) [9].

Only a small proportion of CpGs are common among three of the four epigenetic ages (Supplementary Figure 1), and there is only weak/moderate correlation among them [2, 8, 9].

To our knowledge, there has been no study of DNAm age in type 1 diabetes (T1D). However, telomere length, another indicator of aging, has been investigated in T1D, and it was found to be significantly shorter in T1D compared to non-diabetic subjects [10]. Shorter telomere length has also been associated with T1D duration, pulse pressure [10], BMI [11], systolic blood pressure, all-cause mortality [12], and diabetic nephropathy [13] in subjects with T1D.

Here, we investigated DNAm age calculated by all four methods in 499 subjects with T1D from the Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications (DCCT/EDIC) study and its association with diabetic complications (CVD, nephropathy, retinopathy, and neuropathy) and their risk factors [14,15,16]. We also examined DNAm age in a smaller subset at two time points, 1617 years apart, to investigate changes in DNAm age over time.

Results

Illumina whole blood EPIC data

Comparison of the four DNAm ages with chronological age and with each other

Characteristics of the subjects with EPIC DNAm data are summarized in Table 1. All four epigenetic ages were highly correlated with chronological age and with each other. However, there were significant difference among them: GrimAge was higher than Pan-tissue, and both were higher than chronological age whereas Skin & Blood and PhenoAge were both lower than chronological age, and PhenoAge was lower than Skin & Blood (GrimAge > Pan-tissue > chronological age > Skin & Blood > PhenoAge) (all p < 0.0001) (Table 2, Fig. 1, Supplementary Figure 2).

Table 1 Characteristics of the subjects with Illumina EPIC array data
Table 2 Correlation and mean difference between epigenetic ages and chronological age in EPIC data
Fig. 1
figure 1

DNAm age vs. chronological age in EPIC data. Corr, Spearman correlation coefficient. The line is X = Y

Changes in EAA by chronological age

The differences between chronological age and each of Pan-tissue, Skin & Blood, and GrimAge decreased by 0.13, 0.07, and 0.14 years per 1-year increase in chronological age (p = 1.18E−4, 6.49E−3, and 1.03E−4, respectively). The difference between chronological age and PhenoAge did not differ significantly by chronological age (Fig. 2, Supplementary Table 2).

Fig. 2
figure 2

Epigenetic age acceleration vs. chronological age in EPIC data. The dash line is epigenetic age acceleration = 0. The solid line is the line fitting linear regression model with chronological age as predictor and epigenetic age acceleration as outcome and is present when chronological age is significantly associated with epigenetic age acceleration (p < 0.05)

Association of the four DNAm ages with development of diabetic complications

Epigenetic ages were not significantly associated with the development of CVD or retinopathy. Although DNAm ages were not associated with estimated glomerular filtration (eGFR), both PhenoAge (β = 0.023, p = 0.007) and GrimAge (β = 0.029, p = 0.002) were positively associated with repeated measures of albumin excretion rate (AER, natural log transformed) which remained significant after adjustment of HbA1c levels. The effect of GrimAge on AER increased over time (GrimAge × EDIC year interaction β (SE) = 0.0013 (0.0005), p = 5.10E−3) whereas effect of PhenoAge was not significantly different over time (PhenoAge × EDIC year interaction p = 0.85) (Table 3, Supplementary Figure 3-4). Pan-tissue, Skin & Blood, and PhenoAge were not associated with neuropathy, but GrimAge was positively associated with both diabetic peripheral neuropathy (DPN: OR = 1.07, p = 9.24E−3) and cardiovascular autonomic neuropathy (CAN: OR = 1.06, p = 0.011). These associations also remained significant after further adjustment for time-weighted HbA1c (Table 3).

Table 3 Association of four DNAm ages with development of T1D complication

Association of risk factors of diabetic complication with the four DNAm ages

The univariable associations of different factors with EAAs are shown in Supplementary Table 2. The results regarding multivariable associations of only sex, age, T1D duration, and time-weighted HbA1c with epigenetic ages (minimal model) are shown in Supplementary Table 3. In the multivariable analysis when all factors were included, males had on average 1.5 years higher Pan-tissue (p = 8.00E−4) and GrimAge (p = 9.99E−5) compared to females whereas females had on average 1.5 years higher PhenoAge compared to males (p = 0.005) (Table 4). PhenoAge increased 0.4 years per 1% increase in time-weighted HbA1c (p = 0.026) and 0.01 years per 1-month increase in T1D duration (p = 0.043) (Table 4). The effect of time-weighted HbA1c on PhenoAge was not significantly different in either conventional (β (SD) = 0.45 (0.29), p = 0.12) or intensive (β (SD) = 0.15 (0.48), p = 0.76) therapy group (interaction p = 0.56). Skin & Blood increased 0.09 years per one-unit increase in BMI (p = 0.048) (Table 4). Those with strenuous physical activity at work on average had 1.99 years lower Pan-tissue compared to those with sedentary jobs (p = 0.045), and those who achieved one to two times the recommended level of physical activity during leisure time on average had 0.8 years lower Skin & Blood compared to those who did not achieve the recommended level (p = 0.022) (Table 4). Current smokers had on average 7.1 years higher GrimAge compared to non-smokers (p = 9.03E−50). The other factors were not significantly associated with epigenetic ages (Table 4).

Table 4 Multivariable association of different factors with DNAm ages in the EPIC dataset

Illumina 450K whole blood data

Characteristics of the subjects with 450K data are summarized in Supplementary Table 4. Chronological age and all four epigenetic ages were highly correlated. However, there were significant differences among them: Pan-tissue > (greater than) GrimAge > Skin & Blood > chronological age > PhenoAge (all p < 0.0001) (Supplementary Table 5, Supplementary Figure 5-6). EAAs were not significantly different in the two treatment groups (Supplementary Table 6).

The difference between chronological age and GrimAge decreased by 0.13 years per 1-year increase in chronological age (p = 0.01) (Supplementary Table 6 and Figure 9).

Sex, cohort, treatment group, T1D duration, stimulated C-peptide, and time-weighted HbA1c were not significantly associated with epigenetic ages (Supplementary Table 6, Table 5).

Table 5 Multivariable association of different factors with DNAm ages in 450K dataset

Illumina 450K monocyte data

Chronological age and all four epigenetic ages were highly correlated (Supplementary Table 7, Supplementary Figure 8-9). However, except for Pan-tissue and GrimAge, there were significant differences among them: PhenoAge > Pan-tissue ≈ GrimAge > Skin & Blood > chronological age (p < 0.001).

Although Pan-tissue EAA was on average 2.9 years lower in the former DCCT intensive versus conventional treatment group, the difference was not significant in the multivariable analysis. The other EAAs were also not significantly different between the two treatment groups (Table 5, Supplementary Table 8).

The differences between chronological age and each of Pan-tissue, Skin & Blood, and GrimAge decreased by 0.3, 0.2, and 0.2 years per 1-year increase in chronological age (p = 0.005, 0.003, and 0.004, respectively) (Supplementary Table 8 and Figure 10).

All four EAAs were highly correlated between whole blood and monocyte (Fig. 3).

Fig. 3
figure 3

Monocyte vs. whole blood epigenetic age acceleration in 450K data. EAA, epigenetic age acceleration; Corr, Spearman correlation coefficient. The solid line is X = Y. The dash lines are EAA = 0

In the multivariable analysis, Pan-tissue on average decreased 13.05 years per 1 pmol/ml increase in stimulated C-peptide at DCCT eligibility (p = 0.016). Association of C-peptide with Pan-tissue was not significantly different in the two treatment groups (p = 0.61). Sex, cohort, treatment group, T1D duration, and time-weighted HbA1c were not significantly associated with epigenetic ages (Table 5).

Discussion

We used DNAm data in whole blood measured by Illumina EPIC array and four different methods to estimate epigenetic ages in 499 subjects with T1D. Subsequently, we compared estimated epigenetic ages with chronological age and investigated if the epigenetic ages were associated with development of T1D complications (CVD, nephropathy/decreased renal function, retinopathy, and neuropathy) and their risk factors.

All four epigenetic ages were correlated with chronological age and with each other, but there were significant differences between them. Pan-tissue and Skin & Blood were developed to predict chronological age in healthy individuals whereas PhenoAge and GrimAge used biological biomarkers associated with time-to-death to predict differences in life expectancy of individuals. In addition, they were developed and tested in datasets measured in different tissues by different arrays (Illumina 27K, 450K, and EPIC) and used different statistical methods. As a result, there have been only low to moderate correlations between them [2, 8, 9]. Pan-tissue and GrimAge were significantly higher than chronological age, but Skin & Blood and PhenoAge were significantly lower than chronological age. T1D usually has a negative impact on general health [19, 20]. Therefore, we expected DNAm age to be higher than chronological age in subjects with T1D [21]. However, DCCT subjects were a relatively healthy group of subjects with T1D at baseline due to the extensive inclusion/exclusion criteria applied. Out of ~ 7000 individuals who made initial contact, only 1441 subjects aged 13–39 years with 1–15 years of T1D and no serious long-term complications of diabetes were included in DCCT. Subjects were excluded if they were at risk for adverse effects (e.g., history of frequent ketoacidosis, hypoglycemic coma, or seizure) had known risk factors for vascular complications, were unlikely to comply with the demands of treatment protocols, did not demonstrate an adequate understanding of the DCCT’s purpose, or had drug addiction, chronic alcoholism, or major mental illness [22]. In addition, individuals with EPIC DNAm data are not an entirely random sample of DCCT/EDIC subjects (Supplementary Tables 9-13).

With the exception of PhenoAge, the difference between the epigenetic ages and chronological age decreased in older subjects. This is consistent with previous findings where longitudinal changes in Pan-tissue EAA were tracked using linear mixed models (LMMs) within five different cohorts and showed that epigenetic age increases at a slower rate than chronological age across the life span especially in older populations [23].

We did not find significant association between epigenetic ages and development of CVD. Studies have investigated the association of Pan-tissue and risk of developing CVD in non-diabetic subjects. In 2543 African Americans, the hazard ratio of fatal coronary heart disease increased by 1.03 per year increase in Pan-tissue over a ≈ 17-year follow-up period (144 events) [21]. Similarly, in a cohort of 1863 subjects aged 50–75 years from Germany, the risk of CVD mortality increased by 1.04 per year increase in Pan-tissue over 13 years of follow-up (194 events) [24]. However, our study had only 22% power to detect an effect with this size, due to the relatively small number (N = 58) of CVD events. Another study did not find any significant association between Pan-tissue and incidence of CVD in women from African, Caucasian, and Hispanic ancestry [25]. This discrepancy could be due to race, gender, age, etc. differences between the studied populations as well as statistical power. Increases in PhenoAge and GrimAge have also been reported in association with increased risk of CVD in non-diabetic subjects (β = 1.10 and HR = 1.07, respectively) [8, 9]. However, our study was again under-powered to detect these effects (power = 0.55 and 0.72, respectively). It is noteworthy that in all these studies including the current study, Pan-tissue was calculated in whole blood which is not the target tissue for diabetic complication. Although it has been shown that DNAm profiles are quite similar in different tissues, there are tissue-specific differentially methylated regions which can affect DNAm calculation [26]. However, accessing target tissues are not feasible especially in large scale studies of living subjects.

To our knowledge, association of DNAm age with risk of developing nephropathy, retinopathy, or neuropathy has not been investigated before. We did not observe any significant association between DNAm age and development of nephropathy and retinopathy. However, we found higher GrimAge and PhenoAge to be associated with higher levels of repeated measures of AER, an indicator of decreasing renal function. These associations remained significant even after adjusting for time-weighed HbA1c indicating that at least part of the effects is independent of HbA1c levels. This result is consistent with previous findings where GrimAge has been associated with albuminuria in non-diabetic subjects [9]. Of the four epigenetic ages, higher GrimAge was also associated with development of both DPN and CAN. Larger number of events and higher statistical power along with relatively larger effect sizes could be among the reasons that we detected association of GrimeAge with neuropathy but not the rest of the complications, although CSME has similar number of cases.

Intensive therapy and keeping HbA1c levels close to the normal range has been associated with decreased risk of developing diabetic complications [27]. In our study, both time-weighted HbA1c (but not treatment group) and T1D duration were associated with higher PhenoAge but not with the other epigenetic ages.

Males had higher Pan-tissue and GrimAge compared to females consistent with previous findings [9, 21, 25] whereas PhenoAge was higher in women and sex was not associated with Skin & Blood. All Pan-tissue and PhenoAge CpGs are on autosomal chromosome. Only one of Skin & Blood CpG (cg01892695) is on the X chromosome. Therefore, sex differences in Pan-tissue and PhenoAge are not due to CpGs being on the X chromosome. The CpGs included in GrimAge are not publicly available.

The majority of known risk factors for diabetic complications were not associated with the four epigenetic ages. Two studies have investigated the association of CVD risk factors with Pan-tissue in non-diabetic subjects, but they reported conflicting results [21, 25]. In our study, individuals with physically strenuous jobs had significantly lower Pan-tissue compared to individuals with sedentary jobs, and those who achieved one to two times the recommended level of physical activity during leisure time had lower Skin & Blood compared to those who did not achieve the recommended level. Consistent with these findings, a twin study found that although only a small amount of variance in Pan-tissue is explained by non-shared environmental factors in younger individuals, leisure time physical activity can affect Pan-tissue during adult years [28]. Higher BMI was associated with higher Skin & Blood consistent with what has been observed in non-diabetic subjects before [2]. Current smoking had a large effect on GrimAge; GrimAge increased on average 7 years in current smokers compared to non-smokers. This finding was expected since pack-years were one of the surrogate biomarkers used to generate GrimAge [9]. Triglyceride, HDL, BMI, and physical exercise have been correlated with PhenoAge and/or GrimAge [8, 9]. In our study, although some of these factors were associated with PhenoAge and/or GrimAge in the univariable analysis, these associations did not remain significant in the multivariable analysis [22].

We also investigated epigenetic age in a smaller subset (N = 63) using Illumina 450K array in whole blood and 16–17 years later in monocytes which gave us the opportunity to investigate the difference between the two arrays and the change in epigenetic age over time. In 450K whole blood data, on average, Pan-tissue had the highest estimated value whereas in EPIC data GrimAge was estimated higher than the other epigenetic ages. In addition, Skin & Blood, which was on average lower than chronological age in EPIC data, was on average higher than chronological age in 450K data. Therefore, it appears that the type of array can affect the estimated epigenetic age. However, this could also be due to subjects being highly selected (two extremes of HbA1c and complications risk) in 450K data. Epigenetic ages in whole blood and monocyte measured 16–17 years apart were always correlated indicating that DNAm age is consistent over time and in multiple cell types. This is consistent with previous findings showing that Pan-tissue is consistent across life span [23], and a substantial amount (over 70%) of its changes are shared between different tissue/cell types [29] and also PhenoAge being correlated in different tissues and cell types [8]. In monocyte 450K data, we found a new association which we did not observe in EPIC data: higher stimulated serum C-peptide at DCCT eligibility was associated with lower monocyte Pan-tissue measured decades later at EDIC year 16–17. The association was in the expected direction since preserved beta cell function as measured by stimulated C-peptide is associated with better clinical outcomes (i.e., better glycemic control [30] and lower risk for diabetic complications [31,32,33,34]). The fact that this association was not observed in EPIC data with a much larger sample size could partly be due to the fact that the two sample populations are different: the 450K sample was selected from two extremes of HbA1c and complications risk, whereas the EPIC sample was selected randomly from each cohort/treatment group. In addition, probes may perform slightly differently between the two arrays (450K and EPIC), and as a result measured levels of methylation can be different [35]. In addition, whole blood is a mixture of different cells with different half-lives which could dilute the association(s) of individual cell types especially if they make up only a small proportion (e.g., monocytes, 2–8%).

We investigated the physical distance between all CpGs that are included in Pan-tissue, Skin & Blood, and PhenoAge epigenetic age calculations and T1D GWAS loci [36, 37]: all CpGs are > 25 Mb away from them. Therefore, it is unlikely that estimated DNAm ages were confounded by methylation levels of the CpGs associated with T1D. However, the CpGs included in GrimAge are not publicly available. There have been two epigenome-wide association studies of T1D [38, 39]; however, associated CpGs (N = 132) are available for one of them [38]. Of these, only two CpGs are common with Pan-tissue (cg02047577 (Chr20, 62,587,702 (HG19)) and cg16494477 (Chr5, 170,847,251)): the later CpG and another site are in common with PhenoAge (cg11903057 (Chr4, 40,198,776) and cg16494477 (Chr5, 170,847,251)). Therefore, it appears unlikely that they can have major effect on DNAm age.

Conclusions

We found that although all four epigenetic ages are correlated with each other and with chronological age, however, there are significant differences between them in subjects with T1D. We also found that DNAm age is consistent over time, but type of array (450K vs. EPIC) and cell type can affect the estimated epigenetic age. None of the epigenetic ages were associated with CVD or retinopathy, but PhenoAge and GrimAge were both associated with decreasing renal function as measured by AER. GrimAge was also associated with both DPN and CAN. Only PhenoAge was positively associated with HbA1c levels and T1D duration, two major risk factors for diabetic complications. Some of the other risk factors of diabetic complications were associated with individual epigenetic ages. Therefore, it seems that the investigated epigenetic ages all work sub-optimally in detecting subjects with T1D who are at higher risk to develop complications. However, PhenoAge and specifically GrimAge performed better that Pan-tissue and Skin & Blood suggesting that including biomarkers associated with aging-related mortality improves the accuracy of DNAm measurement. Nevertheless, only some of the risk factors of diabetic complication which are also among the main factors associated with aging-related mortality in the general population were considered in their development (serum creatinine and glucose in PhenoAge and smoking pack-years in GrimAge), and major factors such as hypertension, lipid levels, BMI, and HbA1c which is a better indicator of glycemic levels in long-term compared to serum glucose were not considered [8, 9]. Including these factors could potentially improve epigenetic age estimation in the general population and specifically in subjects with T1D.

Methods

Subjects

The study subjects were from the DCCT/EDIC study. Subjects with T1D aged 13–39 years were recruited into DCCT in 1983–1989 in two cohorts. The primary prevention cohort included participants with 1–5 years of diabetes and no pre-existing retinopathy or nephropathy. The secondary cohort included participants with 1–15 years of diabetes and pre-existing mild retinopathy. Subjects were randomly assigned to receive intensive or conventional diabetes therapy [40]. The DCCT ended in 1993, and subjects subsequently have been followed annually through the EDIC study (Supplementary Figure 11).

Genome-wide DNAm measurement, QC, and normalization

Genome-wide DNAm was measured in whole blood by Illumina Infinium Human Methylation EPIC BeadChip array in a subset of DCCT participants (N = 499) representing about 125 randomly selected adult subjects from each cohort-treatment group who provided informed consent and had sufficient DNA through a 2-year time period prior to DCCT closeout (Zhuo Chen, Feng Miao, Barbara H Braffett, John M Lachin, Lingxiao Zhang, XiweiWu, Delnaz Roshandel, Melanie Carless, Xuejun Arthur Li, Joshua D Tompkins, JohnKaddis, Arthur D Riggs, Andrew D Paterson, DCCT/EDIC Study Group, RamaNatarajan: DNA methylation: a mediator of HbA1c-associated complicationdevelopment in Type 1 diabetes, Submitted) . DNAm was also measured by Infinium Human Methylation 450K BeadChip array in a smaller subset from the DCCT/EDIC (N = 63, 22 overlap with EPIC data) in the whole blood during the same 2-year time period prior to DCCT closeout and in monocytes at EDIC follow-up year 16–17. These included 32 subjects from the conventional treatment group with mean DCCT HbA1c > 9.1% (76 mmol/mol) and significant progression of retinopathy and/or nephropathy from the DCCT closeout to EDIC year 10, and 31 subjects from the intensive treatment group with mean DCCT HbA1c < 7.3% (56 mmol/mol) and no development of retinopathy and nephropathy until EDIC year 10 [41]. Three subjects had missing methylation data for monocytes including two subjects from the conventional and one subject from the intensive treatment group.

The R package meffil (https://github.com/perishky/meffil; accessed on December 2018) [42] was used to perform QC and normalization. Samples were removed if their predicted sex based on DNAm did not match their recorded sex or had > 10% probes with detection p value > 0.01 or > 10% probes with bead number < 3. Samples were also excluded if their SNP genotypes did not match with those from GWAS-array (concordance threshold = 0.8). Illumina HumanCoreExome BeadArray (Illumina, San Diego, CA, USA) data imputed to 1000 Genomes (phase 3, v5) was used for this comparison [43]. One sample from EPIC and one sample from 450K monocyte data did not pass QC. Functional normalization (“noob” for dye-bias and background correction followed by “quantile” normalization implemented in meffil) was then performed to account for technical variation in the data [42]. Blood cell proportions were estimated using the method [44] implemented in meffil with “blood gse35069 complete” as reference [42].

DNAm age calculation

Pan-tissue [1], Skin & Blood [2], PhenoAge [8], and GrimAge [9] were calculated by uploading the data to https://dnamage.genetics.ucla.edu/ (accessed December 2018) with the normalization option selected.

T1D complications

CVD was described as any CVD from DNA collection date to EDIC follow-up year 20 [45].

Nephropathy was described as first occurrence of sustained (2 consecutive) eGFR < 60 ml/min/1.73 m2 [46, 47] or AER ≥ 300 mg/24 h [48] from DCCT closeout to EDIC year 18. The Chronic Kidney Disease Epidemiology Collaboration equation [49] was used to calculate eGFR.

Retinopathy was defined as severe non-proliferative diabetic retinopathy (SNPDR), proliferative diabetic retinopathy (PDR), or clinically significant macular edema (CSME) from DCCT closeout to EDIC follow-up year 18 [50].

Neuropathy was defined as DPN at EDIC year 13/14 and CAN at EDIC year 13/14 and/or EDIC year 16/17 [51].

Statistical analysis

Spearman correlation was used to test for correlation between chronological age and DNAm age and between the four epigenetic ages. Paired sample T tests were used to determine if epigenetic ages were significantly different from chronological age and if the four epigenetic ages differed significantly.

EAA was calculated by subtracting chronological age from epigenetic age (EAA = epigenetic age − chronological age). EPIC array was performed in seven batches. Therefore, batch was included in all multivariable analyses of EPIC data. Since whole blood is a combination of 7 different cell types (neutrophil, B cell, CD4T, CD8T, natural killer cell, eosinophil, and monocyte) and their proportions affect DNAm and epigenetic ages, six predicted cell proportions were included as covariates in all multivariable analyses of whole blood. Neutrophils were excluded as the seven cell proportions sum to one.

Cox proportional hazard models were used to test association of DNAm age with development of CVD, nephropathy, and retinopathy using EPIC data. Logistic regression was used to test association of DNAm age with development of neuropathy during EDIC. Subjects who developed complications during DCCT (before DNAm measurement) were excluded from their corresponding analysis. We also investigated association of DNAm age with repeated measures of renal function, eGFR (annual), and AER (biannual, natural log transformed) from DCCT closeout to EDIC follow-up year 18 using LMMs. Two models were fit for both Cox and LMMs. Model 1 included sex, age, and T1D duration at the time of DNAm measurement along with batch and cell proportions. Model 2 included all covariates in model 1 plus repeated cross-sectional measures of HbA1c (Supplementary Table 12).

In univariable analysis, linear regression was used to test the association of different factors with EAA one at a time. In multivariable analysis, all factors plus chronological age were included in the model, and their associations were tested with DNAm age using linear regression. These factors included sex, cohort, treatment group, and stimulated C-peptide at DCCT eligibility as well as T1D duration, time-weighted HbA1c, BMI, systolic and diastolic blood pressure, HDL, LDL, triglyceride, total cholesterol, pulse rate, current smoking, drinking status, physical activity during work, and leisure time at the time of DNAm measurement (Supplementary Table 3). Since HbA1c and treatment group are highly associated, only HbA1c was included in the multivariable analysis of EPIC data. Due to small sample size and being highly selected on multiple traits, associations of only a subset of factors (sex, age, T1D duration, stimulated C-peptide, HbA1c, cohort, and treatment group) were tested in 450K data.

All the statistical analyses were performed using SAS 9.4 (Cary, NC).

Availability of data and materials

DCCT data are available to authorized users at https://repository.niddk.nih.gov and www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000086.v3.p1.

The 450K whole blood and monocyte DNAm datasets are available in the Gene Expression Omnibus (GEO) database, www.ncbi.nlm.nih.gov/geo (accession no. GSE76171).

The EPIC DNAm dataset will be uploaded into dbGaP https://www.ncbi.nlm.nih.gov/gap/.

Abbreviations

AER:

Albumin excretion rate

CAN:

Cardiovascular autonomic neuropathy

CSME:

Clinically significant macular edema

CVD:

Cardiovascular diseases

DCCT/EDIC:

Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications

DNAm:

DNA methylation

DPN:

Diabetic peripheral neuropathy

EAA:

Epigenetic age acceleration

eGFR:

estimated glomerular filtration

LMMs:

Linear mixed models

PDR:

Proliferative diabetic retinopathy

SNPDR:

Severe non-proliferative diabetic retinopathy

T1D:

Type 1 diabetes

T2D:

Type 2 diabetes

WHR:

Waist/hip ratio

References

  1. Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14(10):R115.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Horvath S, Oshima J, Martin GM, Lu AT, Quach A, Cohen H, et al. Epigenetic clock for skin and blood cells applied to Hutchinson Gilford Progeria Syndrome and ex vivo studies. Aging (Albany NY). 2018;10(7):1758–75.

    Article  CAS  Google Scholar 

  3. Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49(2):359–67.

    Article  CAS  PubMed  Google Scholar 

  4. Bocklandt S, Lin W, Sehl ME, Sanchez FJ, Sinsheimer JS, Horvath S, et al. Epigenetic predictor of age. PLoS One. 2011;6(6):e14821.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Grant CD, Jafari N, Hou L, Li Y, Stewart JD, Zhang G, et al. A longitudinal study of DNA methylation as a potential mediator of age-related diabetes risk. Geroscience. 2017;39(5-6):475–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Nevalainen T, Kananen L, Marttila S, Jylhava J, Mononen N, Kahonen M, et al. Obesity accelerates epigenetic aging in middle-aged but not in elderly individuals. Clin Epigenetics. 2017;9:20.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Horvath S, Erhart W, Brosch M, Ammerpohl O, von Schonfels W, Ahrens M, et al. Obesity accelerates epigenetic aging of human liver. Proc Natl Acad Sci U S A. 2014;111(43):15538–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Levine ME, Lu AT, Quach A, Chen BH, Assimes TL, Bandinelli S, et al. An epigenetic biomarker of aging for lifespan and healthspan. Aging (Albany NY). 2018;10(4):573–91.

    Article  Google Scholar 

  9. Lu AT, Quach A, Wilson JG, Reiner AP, Aviv A, Raj K, et al. DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging (Albany NY). 2019;11(2):303–27.

    Article  CAS  Google Scholar 

  10. Januszewski AS, Sutanto SS, McLennan S, O'Neal DN, Keech AC, Twigg SM, et al. Shorter telomeres in adults with type 1 diabetes correlate with diabetes duration, but only weakly with vascular function and risk factors. Diabetes Res Clin Pract. 2016;117:4–11.

    Article  CAS  PubMed  Google Scholar 

  11. Tesovnik T, Kovac J, Hovnik T, Kotnik P, Battelino T, Trebusak PK. Association of average telomere length with body-mass index and vitamin D status in juvenile population with type 1 diabetes. Zdr Varst. 2015;54(2):74–8.

    PubMed  PubMed Central  Google Scholar 

  12. Astrup AS, Tarnow L, Jorsal A, Lajer M, Nzietchueng R, Benetos A, et al. Telomere length predicts all-cause mortality in patients with type 1 diabetes. Diabetologia. 2010;53(1):45–8.

    Article  CAS  PubMed  Google Scholar 

  13. Fyhrquist F, Tiitu A, Saijonmaa O, Forsblom C, Groop PH. Telomere length and progression of diabetic nephropathy in patients with type 1 diabetes. J Intern Med. 2010;267(3):278–86.

    Article  CAS  PubMed  Google Scholar 

  14. Risk factors for cardiovascular disease in type 1 diabetes. Diabetes. 2016;65(5):1370–9.

  15. Hainsworth DP, Bebu I, Aiello LP, Sivitz W, Gubitosi-Klug R, Malone J, et al. Risk Factors for retinopathy in type 1 diabetes: The DCCT/EDIC Study. Diabetes Care. 2019;42(5):875–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Perkins BA, Bebu I, de Boer IH, Molitch M, Tamborlane W, Lorenzi G, et al. Risk Factors for kidney disease in type 1 diabetes. Diabetes Care. 2019;42(5):883–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Makura CB, Nirantharakumar K, Girling AJ, Saravanan P, Narendran P. Effects of physical activity on the development and progression of microvascular complications in type 1 diabetes: retrospective analysis of the DCCT study. BMC Endocr Disord. 2013;13:37.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Haskell WL, Lee IM, Pate RR, Powell KE, Blair SN, Franklin BA, et al. Physical activity and public health: updated recommendation for adults from the American College of Sports Medicine and the American Heart Association. Med Sci Sports Exerc. 2007;39(8):1423–34.

    Article  PubMed  Google Scholar 

  19. Soedamah-Muthu SS, Fuller JH, Mulnier HE, Raleigh VS, Lawrenson RA, Colhoun HM. All-cause mortality rates in patients with type 1 diabetes mellitus compared with a non-diabetic population from the UK general practice research database, 1992-1999. Diabetologia. 2006;49(4):660–6.

    Article  CAS  PubMed  Google Scholar 

  20. Harjutsalo V, Forsblom C, Groop PH. Time trends in mortality in patients with type 1 diabetes: nationwide population based cohort study. Bmj. 2011;343:d5364.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Roetker NS, Pankow JS, Bressler J, Morrison AC, Boerwinkle E. Prospective study of epigenetic age acceleration and incidence of cardiovascular disease outcomes in the ARIC study (Atherosclerosis Risk in Communities). Circ Genom Precis Med. 2018;11(3):e001937.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Adverse events and their association with treatment regimens in the diabetes control and complications trial. Diabetes Care. 1995;18(11):1415–27.

  23. Marioni RE, Suderman M, Chen BH, Horvath S, Bandinelli S, Morris T, et al. Tracking the epigenetic clock across the human life course: a meta-analysis of longitudinal cohort data. J Gerontol A Biol Sci Med Sci. 2019;74(1):57–61.

    Article  PubMed  Google Scholar 

  24. Perna L, Zhang Y, Mons U, Holleczek B, Saum KU, Brenner H. Epigenetic age acceleration predicts cancer, cardiovascular, and all-cause mortality in a German case cohort. Clin Epigenetics. 2016;8:64.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Horvath S, Gurven M, Levine ME, Trumble BC, Kaplan H, Allayee H, et al. An epigenetic clock analysis of race/ethnicity, sex, and coronary heart disease. Genome Biol. 2016;17(1):171.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Lokk K, Modhukur V, Rajashekar B, Martens K, Magi R, Kolde R, et al. DNA methylome profiling of human tissues identifies global and tissue-specific methylation patterns. Genome Biol. 2014;15(4):r54.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. The Diabetes Control and Complications Trial Research Group. The effect of intensive treatment of diabetes on the development and progression of long-term complications in insulin-dependent diabetes mellitus. N Engl J Med. 1993;329(14):977–86.

    Article  Google Scholar 

  28. Sillanpaa E, Ollikainen M, Kaprio J, Wang X, Leskinen T, Kujala UM, et al. Leisure-time physical activity and DNA methylation age-a twin study. Clin Epigenetics. 2019;11(1):12.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Zhu T, Zheng SC, Paul DS, Horvath S, Teschendorff AE. Cell and tissue type independent age-associated DNA methylation changes are not rare but common. Aging (Albany NY). 2018;10(11):3541–57.

    Article  Google Scholar 

  30. Lachin JM, McGee P, Palmer JP. Impact of C-peptide preservation on metabolic and clinical outcomes in the Diabetes Control and Complications Trial. Diabetes. 2014;63(2):739–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Steffes MW, Sibley S, Jackson M, Thomas W. Beta-cell function and the development of diabetes-related complications in the diabetes control and complications trial. Diabetes Care. 2003;26(3):832–6.

    Article  PubMed  Google Scholar 

  32. Sorensen JS, Johannesen J, Pociot F, Kristensen K, Thomsen J, Hertel NT, et al. Residual beta-Cell function 3-6 years after onset of type 1 diabetes reduces risk of severe hypoglycemia in children and adolescents. Diabetes Care. 2013;36(11):3454–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Kuhtreiber WM, Washer SL, Hsu E, Zhao M, Reinhold P, 3rd, Burger D, et al. Low levels of C-peptide have clinical significance for established type 1 diabetes. Diabet Med. 2015.

  34. Hoeldtke RD, Hampe CS, Bekris LM, Hobbs G, Bryner KD, Lernmark A. Antibodies to GAD65 and peripheral nerve function in the DCCT. J Neuroimmunol. 2007;185(1-2):182–9.

    Article  CAS  PubMed  Google Scholar 

  35. Logue MW, Smith AK, Wolf EJ, Maniates H, Stone A, Schichman SA, et al. The correlation of methylation levels measured using Illumina 450K and EPIC BeadChips in blood samples. Epigenomics. 2017;9(11):1363–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Oram RA, Patel K, Hill A, Shields B, McDonald TJ, Jones A, et al. A type 1 diabetes genetic risk score can aid discrimination between type 1 and type 2 diabetes in young adults. Diabetes Care. 2016;39(3):337–44.

    Article  CAS  PubMed  Google Scholar 

  37. Onengut-Gumuscu S, Chen WM, Burren O, Cooper NJ, Quinlan AR, Mychaleckyj JC, et al. Fine mapping of type 1 diabetes susceptibility loci and evidence for colocalization of causal variants with lymphoid gene enhancers. Nat Genet. 2015;47(4):381–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Rakyan VK, Beyan H, Down TA, Hawa MI, Maslau S, Aden D, et al. Identification of type 1 diabetes-associated DNA methylation variable positions that precede disease diagnosis. PLoS Genet. 2011;7(9):e1002300.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Stefan M, Zhang W, Concepcion E, Yi Z, Tomer Y. DNA methylation profiles in type 1 diabetes twins point to strong epigenetic effects on etiology. J Autoimmun. 2014;50:33–7.

    Article  CAS  PubMed  Google Scholar 

  40. The Diabetes Control and Complications Trial (DCCT). Design and methodologic considerations for the feasibility phase. The DCCT Research Group. Diabetes. 1986;35(5):530–45.

    Article  Google Scholar 

  41. Chen Z, Miao F, Paterson AD, Lachin JM, Zhang L, Schones DE, et al. Epigenomic profiling reveals an association between persistence of DNA methylation and metabolic memory in the DCCT/EDIC type 1 diabetes cohort. Proc Natl Acad Sci U S A. 2016;113(21):E3002–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Min JL, Hemani G, Davey Smith G, Relton C, Suderman M. Meffil: efficient normalization and analysis of very large DNA methylation datasets. Bioinformatics. 2018.

  43. Roshandel D, Gubitosi-Klug R, Bull SB, Canty AJ, Pezzolesi MG, King GL, et al. Meta-genome-wide association studies identify a locus on chromosome 1 and multiple variants in the MHC region for serum C-peptide in type 1 diabetes. Diabetologia. 2018;2018 May;61(5):1098-1111.

  44. Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13:86.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Braffett BH, Dagogo-Jack S, Bebu I, Sivitz WI, Larkin M, Kolterman O, et al. Association of insulin dose, cardiometabolic risk factors, and cardiovascular disease in type 1 diabetes during 30 years of follow-up in the DCCT/EDIC study. Diabetes Care. 2019.

  46. de Boer IH, Sun W, Cleary PA, Lachin JM, Molitch ME, Steffes MW, et al. Intensive diabetes therapy and glomerular filtration rate in type 1 diabetes. N Engl J Med. 2011;365(25):2366–76.

    Article  PubMed  CAS  Google Scholar 

  47. Molitch ME, Gao X, Bebu I, de Boer IH, Lachin J, Paterson A, et al. Early glomerular hyperfiltration and long-term kidney outcomes in type 1 diabetes: the DCCT/EDIC experience. Clin J Am Soc Nephrol. 2019;14(6):854–61.

    Article  PubMed  PubMed Central  Google Scholar 

  48. DCCT/EDIC research group. Effect of intensive diabetes treatment on albuminuria in type 1 diabetes: long-term follow-up of the Diabetes Control and Complications Trial and Epidemiology of Diabetes Interventions and Complications study. Lancet Diabetes Endocrinol. 2014;2(10):793–800.

    Article  Google Scholar 

  49. Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF 3rd, Feldman HI, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–12.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Lachin JM, White NH, Hainsworth DP, Sun W, Cleary PA, Nathan DM. Effect of intensive diabetes therapy on the progression of diabetic retinopathy in patients with type 1 diabetes: 18 years of follow-up in the DCCT/EDIC. Diabetes. 2015;64(2):631–42.

    Article  CAS  PubMed  Google Scholar 

  51. Braffett BH, Gubitosi-Klug RA, Albers JW, Feldman EL, Martin CL, White NH, et al. Risk factors for diabetic peripheral neuropathy and cardiovascular autonomic neuropathy in the Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications (DCCT/EDIC) study. Diabetes. 2020.

Download references

Acknowledgements

The authors acknowledge the participants and researchers of the DCCT/EDIC. Some of these results were presented as abstracts at the American Society of Human Genetics annual meeting in 2018.

DCCT/EDIC Research Group:

S. Scherer1, F. Miao2, L. Zhang2, J. Brown-Friday6, J. Crandall6, H. Engel6, S. Engel6, H. Martinez6 , M. Phillips6 , M. Reid6 , H. Shamoon6 , J. Sheindlin6 , R. Gubitosi-Klug7, J. Wood7, L. Mayer7, D. Miller7, A. Nayate7, M. Novak7, S. Pendegast7, L. Singerman7, D. Weiss7, H. Zegarra7, E. Brown7, P. Crawford7, S. Genuth7, M. Palmert7, P. Pugsley7, J. Quin7, S. Smith-Brewer7, W. Dahms7, J. McConnell7, C. Beck7, K. Farrell7, P. Gaston7, R. Trail7, V. Monnier7, D. Sell7, C. Strauch7, N.S. Gregory8, R. Hanna8, R. Chan8, S. Kiss8, A. Orlin8, M. Rubin8, S. Barron8, B. Bosco8, D. Brillon8, S. Chang8, A. Dwoskin8, M. Heinemann8, L. Jovanovic8, M.E. Lackaye8, T. Lee8, B. Levy8, V. Reppucci8, M. Richardson8, R. Campbell8, A. Bhan9, J.K. Jones9, D. Kruger9, P.A. Edwards9, H. Remtema9, E. Angus9, A. Galprin9, M. McLellan9, A. Thomas9, J.D. Carey9, F. Whitehouse9, R. Bergenstal10, S. Dunnigan10, K. Gunyou10, M. Johnson10, A. Carlson10, L. Thomas10, R. Birk10, P. Callahan10, G. Castle10, R. Cuddihy10, M. Franz10, D. Freking10, L. Gill10, J. Gott10, P. Hollander10, D. Kendall10, J. Laechelt10, S. List10, W. Mestrezat10, J. Nelson10, B. Olson10, N. Rude10, M. Spencer10, D. Etzwiler10, K. Morgan10, L.P. Aiello11, E. Golden11, P. Arrigg11, R. Beaser11, L. Bestourous11, J. Cavallerano11, R. Cavicchi11, O. Ganda11, O. Hamdy11, T. Murtha11, D. Schlossman11, S. Shah11, G. Sharuk11, P. Silva11, P. Silver11, M. Stockman11, J. Sun11, E. Weimann11, L.M. Aiello11, V. Asuquo11, A. Jacobson11, R. Kirby11, L. Rand11, J. Rosenzwieg11, H. Wolpert11, D.M. Nathan12, M.E. Larkin12, M. Cayford12, A. deManbey12, L . Gurry12, J. Heier12, A. Joseph12, F. Leandre12, K. Martin12, C. Shah12, C. Stevens12, N. Thangthaeng12, E. Anderson12, H. Bode12, S. Brink12, M. Christofi12, C. Cornish12, D. Cros12, S. Crowell12, L. Delahanty12, K. Folino12, S. Fritz12, C. Gauthier-Kelly12, J. Godine12, C. Haggan12, K. Hansen12, P. Lou12, J. Lynch12, C. McKitrick12, D. Moore12, D. Norman12, M. Ong12, E. Ryan12, C. Taylor12, D. Zimbler12, A. Vella13, G. Ziegler13, A. Barkmeier13, B. French13, M. Haymond13, J. Mortenson13, J. Pach13, R. Rizza13, L. Schmidt13, W.F. Schwenk13, F.J. Service13, R. Woodwick13, R. Colligan13, A. Lucas13, B. Zimmerman13, K. Nickander13, P. Low13, C. Sommer13, M. Lopes-Virella14, L. Luttrell14, L. Spillers14, J. Fernandes14, K. Hermayer14, S. Kwon14, K. Lee14, T. Lyons14, M. Nutaitis14, A. Blevins14, M. Bracey14, S. Caulder14, J. Colwell14, S. Elsing14, A. Farr14, D. Lee14, P. Lindsey14, R. Mayfield14, J. Parker14, N. Patel14, C. Pittman14, J. Selby14, J. Soule14, M. Szpiech14, T. Thompson14, D. Wood14, S. Yacoub-Wasef 14, A. Jenkins14, R. Klein14, G. Virella14, A. Jaffa14, R. Carter14, J. Stoner14, W.T. Garvey14, D. Lackland14, M. Brabham14, D. McGee14, D. Zheng14, R. K. Mayfield14, M. Molitch15, A. Wallia15, D. Adelman15, M. Hartmuller15, S. Ajroud-Driss15, P. Astelford15, A. Degillio15, M. Gill15, L. Jampol15, C. Johnson15, L. Kaminski15, N. Leloudes15, A. Lyon15, R. Mirza15, D. Ryan15, E. Simjanoski15, Z. Strugula15, S. Colson15, B. Schaefer15, O. Kolterman16, S. Mudaliar16, G. Lorenzi16, M. Goldbaum16, T. Clark16, M. Giotta16, I. Grant16, K. Jones16, R. Lyon16, M. Prince16 , R. Reed16, M. Swenson16, G. Friedenberg16, W.I. Sivitz17, B. Vittetoe17, J. Kramer17, M. Bayless17, C. Fountain17, R. Hoffman17, J. MacIndoe17, N. Olson17, H. Schrott17, L. Snetselaar17, T. Weingeist17, R. Zeitler17, R. Miller18, S. Johnsonbaugh18, M. Carney18, D. Counts18, T. Donner18, J. Gordon18, M. Hebdon18, R. Hemady18, B. Jones18, A. Kowarski18, R. Liss18, S. Mendley18, D. Ostrowski18, M. Patronas18, P. Salemi18, S. Steidl18, W.H. Herman19, C.L. Martin19, P. Lee19, J. W. Albers19, R. Pop-Busui19, E.L. Feldman19. D.A. Greene19, M.J. Stevens19, N. Burkhart19, T. Sandford19, J. Floyd19, A. Sarma19, R. Dunn19, T. Speigelberg19, J. Bantle20, N. Flaherty20, D. Koozekanani20, S. Montezuma20, J. Terry20, F. Goetz20, C. Kwong20, L. McKenzie20, M. Mech20, J. Olson20, B. Rogness20, T. Strand20, R. Warhol20, N. Wimmergren20, M. Steffes20, A. Karger20, J. Seegmiller20, V. Arends20, J. Bucksa20, B. Chavers20, A. Killeen20, M. Nowicki20, A. Saenger20, R. Crow20, B. Gloeb20, C. O’Donnell20, S. Thomas20, D. Goldstein21, D. Hainsworth21, S. Hitt21, J. Giangiacomo21, D.S. Schade22, J.L. Canady22, R.B. Avery22, M.R. Burge22, J.E. Chapin22, A. Das22, L.H. Ketai22, D. Hornbeck22, C. Johannes22, J. Rich22, M.L Schluter22, M. Schutta23, P.A. Bourne23, A. Brucker23, S. Braunstein23, B.J. Maschak-Carey23, S. Schwartz23, L. Baker23, T. Orchard24, L. Cimino24, D. Rubinstein24, D. Becker24, B. Doft24, D. Finegold24, K. Kelly24, L. Lobes24, N. Silvers24, T. Songer24, D. Steinberg24, L. Steranchak24, J.Wesche24, A. Drash24, C. Ryan24, B. Burzuk24, E. Cupelli24, M. Geckle24, D. Sandstrom24, F. Thoma24, T. Williams24, T. Woodfill24, J.I. Malone25, A. Morrison25, M.L. Bernal25, P.R. Pavan25, L. Babbione25, T.J. DeClue25, N. Grove25, D. McMillan25, H. Solc25, E.A. Tanaka25, J. Vaccaro-Kish25, S. Dagogo-Jack26, C. Wigley26, S. Huddleston26, A. Patel26, M. Bryer-Ash26, E. Chaum26, A. Iannacone26, H. Lambeth26, D. Meyer26, S. Moser26, M.B. Murphy26, H. Ricks26, S. Schussler26, S. Yoser26, A. Kitabchi26, P. Raskin27, S. Strowig27, YG. He27, E. Mendelson27, RL. Ufret-Vincenty27, M. Basco27, S. Cercone27, L. Sun28, A. Boright28, B. A. Perkins28, B. Zinman28, A. Barnie28, N. Bakshi28, M. Brent28, R. Devenyi28, K. Koushan28, M. Mandelcorn28, F. Perdikaris28, L. Tuason28, D. Daneman28, R. Ehrlich28, S. Ferguson28, A. Gordon28, K. Perlman28, S. Rogers28, J. Palmer29, R. Fahlstrom29, IH. de Boer29, L. Olmos de Koo29, L. Van Ottingham29, S. Catton29, J. Ginsberg29, J. Kinyoun29, J. Brunzell29, J. Purnell29, H. Wessells29, S. Holt29, J. Hotaling29, C. Kim29, Q. Clemens29, J. Brown29, C. McDonald30, M. Driscoll30, J. Bylsma30, T. Sheidow30, W. Brown30, C. Canny30, P. Colby30, S. Debrabandere30, J. Dupre30, J. Harth30, I. Hramiak30, M. Jenner30, J. Mahon30, D. Nicolle30, N.W. Rodger30, T. Smith30, O. Crofford31, M. May31, J. Lipps Hagan31, R. Ramker31, T. Adkins31, A. Agarwal31, C. Lovell31, S. Feman31, R. Lorenz31, L. Survant31, N.H. White32, L. Levandoski32, I. Boniuk32, J. Santiago32, W. Tamborlane33, P. Gatcomb33, K. Stoessel33, J. Ahern33, K. Fong33, P. Ossorio33, P. Ramos33, J. Lachin34, I. Bebu34, B. Braffett34, J. Backlund34, L. Diminick34, X. Gao34, D. Kenny34, K. Klumpp34, M. Lin34, V. Trapani34, K. Anderson34, K. Chan34, P. Cleary34, A. Determan34, L. Dews34, W. Hsu34, P. McGee34, H. Pan34, B. Petty34, D. Rosenberg34, B. Rutledge34, W. Sun34, S. Villavicencio34, N. Younes34, C. Williams34, C. Cowie35, C. Siebert35, A. Harrington36, D. O’Leary36, J. Polak36, L. Funk36, E.Z. Soliman37, M. Barr37, C. Campbell37, S. Hensley37, J. Hu37, L. Keasler37, Y. Li37, T. Taylor37, Z.M. Zhang37, R. Prineas37, B. Blodi38, R. Danis38, D. Lawrence38, H. Wabers38, M. Burger38, M. Davis38, J. Dingledine38, V. Gama38, S. Gangaputra38, L. Hubbard38, S. Neill38, R. Sussman38, K. Cruickshanks38, D. Dalton38, K. Bainbridge38, M. Budoff, S39. Darabian39, P. Rezaeian39, R. Detrano39, M. Fox39, L. Kim39, R. Oudiz39, N. Wong39, J. Lima40, D. Bluemke40, E. Turkbey40, H. Chahal40, A. Jain40, R. Jarboe40, C. Liu40, A. Malayeri40, C. Miao40, R. J. van der Geest40, S. Hazen41, A. Pratt41, W. Tang41, J. Maynard42, J. Farquhar43, M. Moran43, M. Pfiefer43, M. Schumer43, K. McVary44

6 Albert Einstein College of Medicine

7 Case Western Reserve University

8 Weill Cornell Medical College

9 Henry Ford Health System

10 International Diabetes Center

11 Joslin Diabetes Center

12 Massachusetts General Hospital

13 Mayo Clinic

14 Medical University of South Carolina

15 Northwestern University

16 University of California, San Diego

17 University of Iowa

18 University of Maryland, School of Medicine

19 University of Michigan

20 University of Minnesota

21 University of Missouri

22 University of New Mexico

23 University of Pennsylvania

24 University of Pittsburgh

25 University of South Florida

26 University of Tennessee

27 University of Texas Southwestern Medical Center at Dallas

28 University of Toronto

29 University of Washington

30 University of Western Ontario

31 Vanderbilt University

32 Washington University, St. Louis

33 Yale University School of Medicine

34 George Washington University, The Biostatistics Center

35 National Institute of Diabetes and Digestive and Kidney Disease

36 New England Medical Center

37 Wake Forest School of Medicine

38 University of Wisconsin

39 Harbor UCLA Research and Education Institute

40 Johns Hopkins University, National Heart Lung and Blood Institute

41 Cleveland Clinic

42 Veralight

43 Southern Illinois University

44 Northwest University

Funding

Industry contributors have had no role in the DCCT/EDIC study but have provided free or discounted supplies or equipment to support participants’ adherence to the study: Abbott Diabetes Care (Alameda, CA, USA), Animas (Westchester, PA, USA), Bayer Diabetes Care (Tarrytown, NY, USA), Becton Dickinson (Franklin Lakes, NJ, USA), Eli Lilly (Indianapolis, IN, USA), Extend Nutrition (St Louis, MO, USA), Insulet Corporation (Bedford, MA, USA), LifeScan (Milpitas, CA, USA), Medtronic Diabetes (Minneapolis, MN, USA), Nipro Home Diagnostics (Ft Lauderdale, FL, USA), Nova Diabetes Care (Billerica, MA, USA), Omron (Shelton, CT, USA), Perrigo Diabetes Care (Allegan, MI, USA), Roche Diabetes Care (Indianapolis, IN, USA), and Sanofi (Bridgewater, NJ, USA). The DCCT/EDIC has been supported by cooperative agreement grants (1982–1993 and 2012–2017) and contracts (1982–2012) with the Division of Diabetes Endocrinology and Metabolic Diseases of the NIDKK (current grant nos. U01-DK-094176 and U01-DK-094157) and through support by the National Eye Institute, the National Institute of Neurologic Disorders and Stroke, the General Clinical Research Centers Program (1993–2007), and the Clinical Translational Science Center Program (2006–present), Bethesda, MD, USA. Funding for genotyping by Illumina HumanCoreExome was provided by JDRF grant 17-2013-9. The study sponsor was not involved in the design of the study; the collection, analysis, and interpretation of data; writing the report; or the decision to submit the report for publication.

Methylation study is supported by grants from the National Institutes of Health (NIH): DP3 DK106917-01 and R01 DK065073 (to R.N) and the Wanek family project at the City of Hope (to R.N). Research reported in this publication included work performed in the Integrative Genomics Core at the City of Hope (supported by the National Cancer Institute of the NIH under award number P30CA33572) and the Genomics Core at the University of S. California (Dr. Dan Weisenberger, Director).

Author information

Authors and Affiliations

Authors

Consortia

Contributions

All named authors affirm that authorship is merited based on the International Committee of Medical Journal Editors (ICMJE) authorship criteria. DR designed the study, performed the QC/normalization of the methylation data, analyzed the data, and wrote the manuscript. ZC and RN provided the raw methylation data, contributed to interpretation of the results, and reviewed the manuscript. AJC and SBB provided statistical consult through data analysis and interpretation of the results and reviewed the manuscript. ADP designed the study, researched the data, and reviewed/edited the manuscript and is the guarantor of this work. All authors approve the final version to be published.

Corresponding author

Correspondence to Andrew D. Paterson.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was obtained at each center separately and in George Washington University, The Biostatistics Center. All participants consent to participate in the study.

Consent for publication

Not applicable, our manuscript does not contain any individual person’s data in any form.

Competing interests

The authors declare that there is no competing interest associated with this manuscript.

Additional information

Publisher’s Note

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

Supplementary information

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

Roshandel, D., Chen, Z., Canty, A.J. et al. DNA methylation age calculators reveal association with diabetic neuropathy in type 1 diabetes. Clin Epigenet 12, 52 (2020). https://doi.org/10.1186/s13148-020-00840-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-020-00840-6

Keywords