key: cord-0017935-3jx82c9l
authors: Matías-García, Pamela R.; Ward-Caviness, Cavin K.; Raffield, Laura M.; Gao, Xu; Zhang, Yan; Wilson, Rory; Gào, Xīn; Nano, Jana; Bostom, Andrew; Colicino, Elena; Correa, Adolfo; Coull, Brent; Eaton, Charles; Hou, Lifang; Just, Allan C.; Kunze, Sonja; Lange, Leslie; Lange, Ethan; Lin, Xihong; Liu, Simin; Nwanaji-Enwerem, Jamaji C.; Reiner, Alex; Shen, Jincheng; Schöttker, Ben; Vokonas, Pantel; Zheng, Yinan; Young, Bessie; Schwartz, Joel; Horvath, Steve; Lu, Ake; Whitsel, Eric A.; Koenig, Wolfgang; Adamski, Jerzy; Winkelmann, Juliane; Brenner, Hermann; Baccarelli, Andrea A.; Gieger, Christian; Peters, Annette; Franceschini, Nora; Waldenberger, Melanie
title: DNAm-based signatures of accelerated aging and mortality in blood are associated with low renal function
date: 2021-06-02
journal: Clin Epigenetics
DOI: 10.1186/s13148-021-01082-w
sha: 332f177b9a53f4938cf603d7b2dd73fd27768be0
doc_id: 17935
cord_uid: 3jx82c9l

BACKGROUND: The difference between an individual's chronological and DNA methylation predicted age (DNAmAge), termed DNAmAge acceleration (DNAmAA), can capture life-long environmental exposures and age-related physiological changes reflected in methylation status. Several studies have linked DNAmAA to morbidity and mortality, yet its relationship with kidney function has not been assessed. We evaluated the associations between seven DNAm aging and lifespan predictors (as well as GrimAge components) and five kidney traits (estimated glomerular filtration rate [eGFR], urine albumin-to-creatinine ratio [uACR], serum urate, microalbuminuria and chronic kidney disease [CKD]) in up to 9688 European, African American and Hispanic/Latino individuals from seven population-based studies. RESULTS: We identified 23 significant associations in our large trans-ethnic meta-analysis (p < 1.43E−03 and consistent direction of effect across studies). Age acceleration measured by the Extrinsic and PhenoAge estimators, as well as Zhang’s 10-CpG epigenetic mortality risk score (MRS), were associated with all parameters of poor kidney health (lower eGFR, prevalent CKD, higher uACR, microalbuminuria and higher serum urate). Six of these associations were independently observed in European and African American populations. MRS in particular was consistently associated with eGFR (β = − 0.12, 95% CI = [− 0.16, − 0.08] change in log-transformed eGFR per unit increase in MRS, p = 4.39E−08), prevalent CKD (odds ratio (OR) = 1.78 [1.47, 2.16], p = 2.71E-09) and higher serum urate levels (β = 0.12 [0.07, 0.16], p = 2.08E−06). The “first-generation” clocks (Hannum, Horvath) and GrimAge showed different patterns of association with the kidney traits. Three of the DNAm-estimated components of GrimAge, namely adrenomedullin, plasminogen-activation inhibition 1 and pack years, were positively associated with higher uACR, serum urate and microalbuminuria. CONCLUSION: DNAmAge acceleration and DNAm mortality predictors estimated in whole blood were associated with multiple kidney traits, including eGFR and CKD, in this multi-ethnic study. Epigenetic biomarkers which reflect the systemic effects of age-related mechanisms such as immunosenescence, inflammaging and oxidative stress may have important mechanistic or prognostic roles in kidney disease. Our study highlights new findings linking kidney disease to biological aging, and opportunities warranting future investigation into DNA methylation biomarkers for prognostic or risk stratification in kidney disease. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s13148-021-01082-w.

The kidneys are responsible for maintenance of homeostasis and blood filtration, and their function is most commonly clinically assessed by measuring serum creatinine levels to estimate glomerular filtration rate (eGFR) [1, 2] . Chronic kidney disease (CKD), defined by low eGFR (< 60 ml/min/m 2 ) and/or presence of protein in urine, is an increasingly prevalent non-communicable disease with a considerable burden worldwide [3] [4] [5] . Increased urinary albumin-to-creatinine ratio (uACR) is a marker of kidney injury, measured to identify early kidney damage which can precede eGFR decline (for example, diabetic nephropathy) [6] . Albuminuria is a predictor of CKD progression and mortality [6, 7] , whereas high serum levels of urate, a molecule of purine nucleotide metabolism excreted by the kidney, is a risk factor for incident cardiovascular and kidney disease, and also a biomarker of low eGFR [8] .

DNA methylation (DNAm), defined as the covalent addition of a methyl group to a DNA nucleotide (usually the cytosine of a cytosine-guanine dinucleotide [CpG] ), is the most extensively studied epigenetic mechanism, and its role in numerous conditions and diseases has been demonstrated [9] . Age-predicting algorithms based on the percentages of DNAm observed at sets of CpGs, such as the ones proposed by Hannum [10] and Horvath [11] , have been used to predict an individual's age (DNAmAge) and assess biological aging by calculating the difference between an individual's predicted and chronological age-a concept known as DNAmAge acceleration (DNAmAA) [12, 13] . Other measures have been derived to assess specific aspects of aging mechanisms, such as intrinsic epigenetic age acceleration (IEAA), which assesses aging independent of blood immune system changes [14] , or extrinsic epigenetic age acceleration (EEAA) [15, 16] , which specifically estimates aging as related to the immune system and reflected in changes in blood immune cell-type proportions. A "second generation" of DNAm-based aging signatures incorporated physiological markers to better capture changes in traditional biological aging biomarkers [12] . PhenoAge was developed as a marker meant to mirror physiological dysregulation as reflected in changes in age and 9 additional age-related features, such as C-reactive protein and serum glucose [17] . GrimAge is a mortality predictor based on mortality-related DNAm-estimated traits [18] . Another DNAm-based lifespan predictor, the 10-CpG epigenetic mortality risk score (MRS), stands out for its simplicity and its recent validation [19, 20] .

Multiple studies have shown, although with varying findings, a positive relationship between DNAmAge measured in blood and aging-related diseases and mortality [12, 13, 21] . DNAmAge is associated with all-cause mortality [22, 23] , frailty [24] , cognitive function and physical fitness [25] , body mass index (BMI) [26] and obesity [27] , lifetime stress [28] and a number of other age-related conditions [12] . The available evidence points to DNAmAge as a potential global biomarker of biological aging and health, though potential for publication bias must be considered [21] . Although some of the "second-generation" DNAm-based aging measures include proteins or markers known to be associated with kidney function [17, 18] , whether these and other DNAm-based predictors are correlated with different parameters of kidney aging and low function has not been investigated [13, 29, 30] .

We evaluated the association between five kidney traits (eGFR, prevalent CKD, uACR, microalbuminuria and serum urate) and seven DNAm-based age and/or lifespan predictors (HannumAA, HorvathAA, EEAA, IEAA, PhenoAA, GrimAA and MRS) in up to seven populationbased studies in a large trans-ethnic meta-analysis. We also evaluated kidney trait associations with secondary DNAm-based predictors: categorical epigenetic mortality risk score (MRS) variables, and eight DNAm-estimated traits underlying the GrimAge mortality predictor.

[0.07, 0.16], p = 2.08E−06). The "first-generation" clocks (Hannum, Horvath) and GrimAge showed different patterns of association with the kidney traits. Three of the DNAm-estimated components of GrimAge, namely adrenomedullin, plasminogen-activation inhibition 1 and pack years, were positively associated with higher uACR, serum urate and microalbuminuria.

Conclusion: DNAmAge acceleration and DNAm mortality predictors estimated in whole blood were associated with multiple kidney traits, including eGFR and CKD, in this multi-ethnic study. Epigenetic biomarkers which reflect the systemic effects of age-related mechanisms such as immunosenescence, inflammaging and oxidative stress may have important mechanistic or prognostic roles in kidney disease. Our study highlights new findings linking kidney disease to biological aging, and opportunities warranting future investigation into DNA methylation biomarkers for prognostic or risk stratification in kidney disease.

Keywords: Aging, Kidney function, Epigenetic age acceleration, DNAm age, Glomerular filtration rate, UACR , Serum urate

We additionally performed ethnic-specific meta-analyses to identify robust associations across cohorts of different ethnicities.

We conducted a trans-ethnic meta-analysis of seven DNAm-based age/lifespan predictors and five kidney traits using data from seven population-based cohorts (Fig. 1 , population characteristics in Table 1 ). Serum creatinine-based traits (eGFR and prevalent CKD) were available for all European ancestry (k = 5), African American (k = 2) and one additional Hispanic/Latino study (N = 9688). Sample sizes for the other traits were smaller: serum urate was available in four studies with European ancestry and one study with African American participants (N = 5903), while uACR and microalbuminuria were available in three of the European ancestry studies and one African American study (N = 4110). Additional details on the cohorts are provided in Additional file 1: Table S1 .

DNAm-based age and lifespan predictors were used as independent variables and kidney traits as dependent variables in covariate-adjusted regression models. Studylevel results showed associations were slightly attenuated after the inclusion of additional covariates in model 2 (namely BMI, log-transformed triglycerides, HDL, hypertension, smoking status and diabetes) in comparison with model 1 (basic model adjusting for chronological age and sex; Additional file 1: Table S2 ). Sensitivity analyses showed that "crude" bivariate correlations were largely attenuated after adjustment for chronological age; the introduction of additional variables did not significantly alter the observed effects, despite a slight increase in the coefficient after the introduction of smoking in models with MRS and GrimAge (Additional file 3: Note Fig. 1 Study design. Regression analyses were conducted in each study by modeling DNAm-based predictors as independent variables and kidney traits as dependent variables, adjusting for confounders. Results from the fully adjusted model (with chronological age, sex, BMI, log-transformed triglycerides, HDL, hypertension, smoking status, diabetes as covariates and baseline eGFR for serum urate analyses) were meta-analyzed using inverse-variance weighted fixed-effects and random-effects models. We based our main interpretations on the fixed-effects results; if heterogeneity was large (I 2 > 0.50 and Cochran's Q p het < 0.05), we based our interpretations on the random-effects results. The Venn diagram shows the set of 23 statistically significant associations between DNAm-based predictors and kidney traits identified in the trans-ethnic meta-analysis (p < 1.43E−03 and consistent direction of effect across studies) S1). Estimates from the fully adjusted model were metaanalyzed using inverse-variance weighted fixed-effects and random-effects models. We based our main interpretations on the fixed-effects results; if heterogeneity was large (I 2 > 0.50 and Cochran's Q p het < 0.05), we based our interpretations on the random-effects results (Methods).

We identified 23 significant DNAm-based predictorkidney trait associations (p < 1.43E−03 and concordant direction of effect across studies; Table 2 ). Three interesting groups worth further discussion are: (1) PhenoAA, MRS and EEAA were associated with all parameters of poor kidney health (lower eGFR, prevalent CKD, higher uACR or microalbuminuria and higher serum urate); (2) the "first-generation" epigenetic aging markers, where HannumAA was associated with all kidney traits but serum urate, and HorvathAA was only associated with lower eGFR; and (3) an analogous measure to age acceleration in GrimAge was associated with uACR, microalbuminuria and serum urate (Fig. 1) . Six associations between DNAm-based predictors and kidney traits were replicated across ethnic groups, and ethnic-specific replication was observed for 16 associations in total.

From these three DNAm-based predictors, associations with MRS had the smallest p-values: one MRS unit increase had a − 0.12 (95% CI = [− 0.16, − 0.08]) change in one standard deviation (SD) of log-transformed eGFR (p = 4.39E−08) and was associated with 78% [47-116%] increased odds of prevalent CKD (p = 2.71E−09) ( Table 2 ). Although high heterogeneity was identified in associations between MRS and CKD, uACR and Table 1 Population characteristics Population characteristics of all participating studies. The means and standards deviation (SD) are shown for continuous traits, and N (%) for categorical traits. **Skewed variables, for which median and (1st, 3rd quartile) are shown. The sample sizes presented here for each of the studies correspond to the number of observations with information on DNAm-predictors, serum-based creatinine kidney traits, chronological age and sex. Age was measured in years at time of participation in study; BMI in kg/m 2 ; serum creatinine in mg/dL; eGFR, serum-creatinine-based estimated glomerular filtration rate in mL/min/1.73 m 2 ; CKD: prevalent chronic kidney disease, defined as eGFR < 60 ml/min/1.73 m 2 ; uACR in mg/g; microalbuminuria was defined as uACR ≥ 30 mg/g; serum urate in mg/dl; prevalent diabetes was defined based on use of glucose lowering drugs or fasting plasma glucose ≥ 126 mg/dl; hypertension defined using the Joint National Committee (JNC) VII definition (blood pressure > 140/90 mm Hg or use of anti-hypertensive medications); HDL cholesterol and triglycerides in mg/dl; C-reactive protein in mg/L; NA denotes the trait was not available. Abbreviations used in ethnic background row: Eur., European ancestry; Afr.Am.., African American; His., Hispanic/Latino. An extended version of this table is shown in Additional file 1: Figure 2 shows the study-level regression coefficients of the association between the continuous kidney traits and DNAm-based predictors standardized to one SD deviation in both terms to allow for their comparison. While the strength of association with serum urate of these three DNAm-based predictors was similar, MRS and EEAA had the largest effects on eGFR (Fig. 2 , standardized effects in Additional file 1: Table S4 ). Similar observations were done for the binary kidney traits (Additional file 2: Fig. S1 ). Figure 3 shows results from ethnic-specific meta-analyses significant in both European ancestry and African American meta-analyses (p < 1.43E−03, Additional file 1: Table S3 ). MRS stands out for its replication across all subgroups, including the small Hispanic/Latino cohort (Fig. 3a , study-level results in Additional file 1: Table S2 ). The eGFR-EEAA and CKD-PhenoAA effects were also replicated at the Bonferroni-corrected significance level across ethnic-specific meta-analyses (Fig. 3b) , whereas Table 2 Trans-ethnic meta-analyses of associations between kidney traits and DNAm-based age and lifespan predictors Results from trans-ethnic meta-analyses of associations between kidney traits and DNAm-based age and lifespan predictors in up to seven population-based studies. Study-level associations were adjusted for chronological age, sex, BMI, blood lipids, hypertension, smoking and diabetes. Beta coefficients are given as changes in one standard deviation (SD) of the continuous kidney trait. Fully adjusted associations of serum-creatinine-based traits (eGFR, CKD) are based on N ≤ 9390 observations, whereas the sample size for urinary albumin-based traits (uACR, microalbuminuria) is N ≤ 4406 and for urate N ≤ 5769. I 2 is the heterogeneity statistic, and (Q) p het corresponds to Cochran's Q heterogeneity statistic Shown in bold are statistically significant associations (p < 1.43E-03 and consistent direction of effect across studies) with either no evidence of heterogeneity in the fixed-effects model or supporting findings from the random-effects model a p het of MRS with CKD, uACR and microalbuminuria < 0.05, therefore reported association based on significant random-effects models: MRS-uACR: their "complementary" associations (CKD-EEAA and eGFR-PhenoAA) were nominally significant (Additional file 2: Fig. S2A ). Further effects observed in both European ancestry meta-analysis and the African American cohort were the associations of higher serum urate with MRS and PhenoAA (Fig. 3c ). On the other hand, the associations between EEAA, PhenoAA and MRS with higher uACR (and prevalent microalbuminuria) identified in the trans-ethnic meta-analysis were mostly driven by the effects from the European ancestry cohorts (Additional file 2: Figs. S3 and S4), as well as that of serum urate and EEAA (Additional file 2: Fig. S5A ). Most notably, the association between uACR and MRS was replicated at the Bonferroni-corrected level in two of the European ancestry cohorts (KORA and ESTHER-II), and nominally significant in the third one (Additional file 1: Table S2 ).

HorvathAA, age acceleration measured by the "first-generation" DNAm-based predictor HorvathAge, was exclusively associated with low eGFR: a one year difference between DNAm-estimated age and chronological age was associated with a − 0.006 [− 0.01, − 0.003] change in log-transformed eGFR (p = 5.15E−04, Table 2 ). Han-numAA, another "first-generation" DNAm-based age predictor, was also associated with eGFR; the strength of the association with both "first-generation" DNAmbased predictors was similar (Fig. 2) . HannumAA was additionally associated with CKD, uACR and microalbuminuria (Additional file 2: Figs. S2-S4B). While the effects observed for HannumAA with eGFR and CKD were robustly replicated by cohorts with African American participants (Additional file 2: Fig. S2B ), the association between HorvathAA and eGFR was mostly driven by studies with European ancestry (Additional file 2: Fig.  S2C) .

A year of GrimAge acceleration was associated with an increase of 0.03 SD of log-transformed uACR and a 10.6% increase in the odds of having microalbuminuria, both early markers of renal damage (Additional file 2: Figs S3-S4C). Likewise, a one-year difference in GrimAge was associated with higher serum urate levels, a risk factor for cardiorenal disease (Additional file 2: Fig. S5B ). The GrimAA effects on uACR were the largest across DNAm-based predictors and (Fig. 2) , and a similar effect was observed for microalbuminuria (Additional file 2: Fig. S1 ).

In the categorical MRS, where risk groups were defined based on the number of CpGs with methylation levels beyond pre-defined high-risk thresholds, a 0.431 (95% CI = [0.29, 0.57]) SD increase in log-uACR was observed for individuals with > 5 "aberrantly" methylated CpGs (high risk MRS) compared to those with 0-1 aberrantly methylated loci (low risk MRS) (p = 4.01E-09, Additional file 2: Fig. S6A ). Similar associations were observed for microalbuminuria (OR = 2.075 [1.582, 2.722], p = 1.36E-07; Additional file 1: Table S3 ). The moderate MRS category (defined as those with 2-5 "aberrantly" methylated loci) was also associated with serum urate (Additional file 2: Fig. S6B ). These effects were mostly driven by European ancestry studies, where effects had a similar strength of association and replicated within this group Multivariate regression meta-analysis of association between kidney traits and MRS, EEAA and PhenoAA. Multivariate regression models were used to assess the relationship between kidney traits and age acceleration as measured by seven DNAm-based predictors of age and/or lifespan, including Zhang's 10-CpG mortality risk score (MRS), PhenoAge and extrinsic epigenetic age acceleration (EEAA). Results from the fully adjusted model included chronological age, sex, BMI, log-transformed triglycerides, HDL, hypertension, smoking status and diabetes as covariates, and baseline eGFR for serum urate analyses. Inverse-variance weighted fixed-effects models meta-analysis was conducted, where if heterogeneity was observed ([Q] p het < 0.05), a random-effects model was further interpreted. Individual panels show forest plots with the study-level and meta-analytic results for each association between kidney trait and DNAm-based predictor from the fully adjusted model. The rows correspond to the different studies and the sample size for each analysis (N). For eGFR and urate, the regression estimates represent the change in one standard deviation of the kidney trait per unit change in the MRS or per one year of age acceleration for PhenoAA and EEA. For CKD, the estimate column corresponds to the odds ratio (OR). The x-axis shows the estimates obtained from either the regression model (for single studies, data point shape is a black square) or the meta-analytic estimate (data point shape is a green diamond for serum creatinine-based traits and a red diamond for urate) with their 95% CI. Estimate: regression coefficient for the continuous traits and OR for CKD, 95% CI: 95% confidence interval of the estimate (Additional file 2: Fig. S7 ). None of these effects replicated in the African American cohort, perhaps due to the smaller sample size and reduced statistical power.

In regard to the secondary analyses with the eight DNAm-based components that constitute GrimAge, we found associations for three of them. DNAm-estimated adrenomedullin (DNAmADM), plasminogen-activation inhibition 1 (DNAmPAI1) and pack years (DNAm-PACKYRS) were positively associated with higher uACR, higher serum urate levels and microalbuminuria (Additional file 1: Table S5 ). Of note, the effects identified for serum urate replicated in both the European-ancestry analyses and in the African American cohort (Additional file 1: Table S2 ), whereas those in uACR and microalbuminuria were mostly replicated within the European studies.

We identified 23 associations between kidney traits and DNAm-based predictors of aging and/or mortality in a large trans-ethnic meta-analysis from up to seven multiethnic population-based cohorts. PhenoAA, MRS and EEAA were associated with all parameters of poor kidney health (lower eGFR, prevalent CKD, higher uACR or microalbuminuria and higher serum urate). Distinct patterns of association were observed with age acceleration in the "first-generation" clocks (Hannum and Horvath) and an analogous measure in GrimAge: HorvathAA was only associated with lower eGFR, while HannumAA was associated with all kidney traits but serum urate. Finally, GrimAA was associated with uACR, microalbuminuria and serum urate.

Sensitivity analyses showed that the inclusion of chronological age in a "crude" model lead to large attenuation of the correlation between DNAm-based predictors and kidney traits, as expected when adjusting for variables confounding the exposure-outcome association. The posterior introduction of additional variables did not further alter the observed effects. A slight increase in the coefficient after the introduction of smoking was observed for MRS and GrimAge, which may be explained in general by the strong correlation between smoking and DNAm/ DNAmAge acceleration in blood [31] [32] [33] [34] , and specifically in relation to the these two markers as they either directly incorporate cigarette smoking into its formulation [19] or capture the effects of cigarette smoking [20] .

From the three "universally" associated DNAmbased predictors, MRS was robustly associated with multiple kidney phenotypes, some of which replicated across ethnic-specific analyses (i.e., eGFR, prevalent CKD and urate). Together with EEAA and GrimAA, MRS showed the largest effects in comparison with the other DNAm-based predictors. Moreover, MRS moderate and high risk categories [20] were also associated with eGFR and markers of kidney injury in the secondary analyses. Recent evidence suggests the MRS is a DNAm-based biomarker which reflects mortality risks by capturing the effects of oxidative stress and systemic inflammation, as well as inflammation-driven changes in immune cell counts-all mechanisms shared by numerous chronic diseases [35] [36] [37] [38] . Particularly in kidney aging, mitochondrial dysfunction, uremia-induced epigenetic changes and the production of reactive oxygen species in the glomeruli that lead to barrier function impairment and albuminuria are all mechanisms promoting oxidative stress [8, [39] [40] [41] . These factors may thus explain the associations between kidney traits and MRS here identified. The grounding of the MRS in inflammation and oxidative stress mechanisms may also explain the predictive power of this predictor in regard to outcomes from cancer to cardiovascular disease mortality [20, 21] .

EEAA and PhenoAA, also associated with all kidney traits, are both extrinsic aging measures (i.e., tracking changes in blood cell composition) that also capture (intrinsic) aging-related physiological dysregulation [12, 18] . EEAA, identified as a measure of immune system aging, is considered better at predicting age-related decline of tissue performance [12, 14, 18] . Increased allostatic load, activation of stress and pro-aging pathways, impairment of protective pathways as well as exogenous lifestyle and environmental factors are all factors driving premature aging [17, 39] . Immunosenescence [42, 43] , systemic low-grade inflammation ('inflammaging') [44] and oxidative stress [37, 39] contribute to an increased allostatic load and are mechanisms present in kidney aging [17, 45] . In particular, the immune system plays an important role in the incidence, development and resolution of renal disease [43] . A signature of 447 genes involved in renal aging further confirmed the relevance of cellular pathways common to immune function and renal physiological decline [46] . Likewise, CKD patients show signs of premature immunological aging (such as poor naïve T-cell frequency and reduced thymic output), a status induced by ureamia (high concentrations of serum urate in blood) that is associated with poor clinical outcomes [47] . At the cellular level, telomere attrition [39, 48] and the cytokine secretory profile of senescent cells promote inflammation and lead to fibrotic damage [37, 39, 43, 45] , further linking inflammaging to renal dysfunction [9] . Our findings, in line with the positive association between PhenoAge and albumin excretion rate identified in 499 subjects with type 1 diabetes [49] , may be thus explained by the relationship between these DNAm-based predictors and the aging-related changes to the immune system, low-grade chronic inflammation and oxidative stress that impact renal disease [37, 47] . Moreover, EEAA, PhenoAge and the MRS have shown a stronger predictive association with time to death than HannumAA, HorvathAA and IEAA, suggesting they better reflect mortality risks associated with biological aging [16, 20] -and based on our findings, also better reflect immune system changes associated with kidney aging. The correlation between these DNAm-based predictors and poor kidney health may also explain their strong associations with mortality, as the ensuing contribution of renal disease to physiological dysregulation may increase mortality risk [6, 7, 39, 50] .

Distinct patterns of association were identified with two of the "first-generation" clocks, HorvathAA and HannumAA. HorvathAA, thought to track cell-intrinsic aging (e.g., epigenetic stability mechanisms, cell growth and survival as well as organismal development) [11, 12, 16, 18] , was exclusively associated with eGFR in our study. The statistical power derived from the larger sample size in our study likely explains this positive association, unlike prior studies reporting null findings with eGFR [49, 51] . Considering tissue from individuals with renal disease was included in the development of this pan-tissue marker [11] , this epigenetic marker reflects alterations in ubiquitous cell-intrinsic pathways that our findings suggest may be relevant for renal (glomerular) function. Genome-wide association studies of the Hannum and Horvath DNAm-based predictors have shown that, even though they capture different aspects of aging, both markers are influenced by genes associated with metabolic and immune system pathways [52, 53] . Consequently, these DNAm-based predictors may also somewhat be reflective of the immune molecular mechanisms previously described. An additional factor of potential relevance for the eGFR-HorvathAA association may be aberrant glucocorticoid signaling, given its potentially pathogenic role in renal function [54] and the enrichment of glucocorticoid response elements in this DNAm marker [29] .

HannumAA, associated with lower renal function and markers of early renal damage (uACR and microalbuminuria) in our study, is also strongly correlated with blood cell counts [16, 23] and is sensitive to environmental influences [15, 53] . Moreover, it has also been associated with higher levels of inflammatory biomarkers, creatinine and certain lipid classes in individuals of European ancestry [14, 15, 55] . Like other extrinsic measures, Hannum seems to be a better marker for later-life diseases and mortality than Horvath [22] and has even been proposed as a prognostic marker of pathological metabolic processes [56] . HannumAA in cancerous kidney tissue vs normal samples has also been reported [10] , thus offering further evidence on the relevance of our findings in blood to renal disease.

The association between one-year difference in the GrimAge predictor and higher uACR, serum urate and microalbuminuria-but not eGFR or prevalent CKDsuggests this DNAm-based predictor might be more sensitive to systemic inflammation and early renal damage. Albuminuria changes can occur before eGFR decline in early kidney disease [6] : early structural glomerular lesions in patients with normal eGFR are better correlated with changes in uACR than with GFR decline, where the latter might not be present yet [57] . Moreover, albuminuria is a predictor of CKD progression and mortality independently from eGFR changes, which suggests they represent two independent mechanisms underlying renal disease progression [6, 7] . Our findings are in line with prior reports of GrimAge association with albumin excretion in T1D patients and non-diabetic subjects [19, 49] , where lower power may explain the null findings reported in [51] . Asymptomatic hyperuricaemia, or high serum urate levels, is involved in pro-inflammatory mechanisms and is associated with a high risk of cardiovascular and renal disease [8] . Serum urate, both in its crystal and soluble forms, activates innate immunity and triggers DNAm epigenetic mechanisms (e.g., promoting cytokine secretion, pro-inflammatory pathways including oxidative stress) leading to persisting inflammation and an increased allostatic load [8, 39] . These effects may, in turn, explain our findings in relation to GrimAge. A high degree of heterogeneity in the eGFR-and CKD-GrimAA analyses was observed, similarly to reports in prior studies [58, 59] .

Overall, the observed associations did not show a clear pattern across kidney traits, thus supporting the proposed notion that the existing DNAm-based predictors might reflect different aspects of biological aging. This is in line with their differential association with risk factors, intermediate phenotypes and diseases [12, 58, 60, 61] , and their inclusion of non-overlapping CpGs sets [19, 20, 60, 62] . The CpG overlap between DNAm-based predictors was assessed in detail by Liu et.al, where the lack of CpG overlap may be explained by the redundancy of the methylome: CpGs selected in the construction of different DNAm predictors may represent different aging hallmarks or pathways, despite the potential biological similarities of their genomic regions [62] . Our findings are also consistent with the many associations observed with EEAA and other blood immune system correlated DNAm-based predictors [12] , and with the lack of associations with IEAA observed in prior studies [52, 61] .

DNAm-estimated adrenomedullin (DNAmADM), plasminogen activator inhibitor-1 (DNAmPA1) and smoking pack years (DNAmPackYears) were positively associated with decreased renal function. Consistent with our findings, patients with chronic cardiorenal diseases have higher blood levels of ADM [63, 64] and PAI-1 [65] [66] [67] , where the first is a potential biomarker of CKD progression [64, 68] and the latter a risk factor for cardiorenal disease [66, 69] . Several factors involved in kidney disease pathogenesis (e.g., oxidative stress, inflammation) induce PAI-1 expression (66, 70, 71) , which in turn has been linked to fibrosis, glomeruli damage and other pathogenic mechanisms in renal disease [65, 67] , as well as to thrombosis and an increased hypercoagulable state-a shared phenotype of inflammaging [37] and renal disease [72] . Moreover, the effects of smoking in DNAm [31, 32] and their association to DNAmAge acceleration in blood are well known [33, 34] . PhenoAge and MRS capture effects of cigarette smoking [18, 20, 33, 58, 60] , whereas GrimAge specifically incorporates cigarette smoking into its formulation [19] . Cigarette smoking has been associated with renal function decline and increased inflammation [73] , where several mechanisms explaining the negative effects of smoking on renal function (e.g., oxidative stress, endothelial dysfunction, immune function modulation) contribute to renal disease progression [74] . Overall, our findings are in line with the roles described in the literature for the studied DNAmestimated markers, and further support the notion that associations between epigenetic aging and health outcomes may be mediated by age-related pro-inflammatory mechanisms [75] . Moreover, they suggest DNAm-based estimates might prove to be valuable proxies in settings where such variables are not available (e.g., limitations in the clinical use of ADM [68] and self-reported smoking [34] ).

Strengths of this work are the large sample size and inclusion of multiple independent studies involving multi-ethnic populations. We comprehensively addressed biological aging and lifespan as predicted by DNAm and assessed multiple kidney traits reflecting different aspects of renal health. Our results in regard to PhenoAge and GrimAge represent confirmatory findings to some extent, given that renal function variables were included in the derivation of these algorithms [18, 19] . Of note, although the MRS was derived using data of two of the cohorts included in this study, it has been independently validated [20, 21] and the replication of its associations across multiple cohorts suggest our findings are not a product of data overfitting. All in all, our study meets the considerations proposed by a recent literature review and meta-analysis on the topic [22] .

Limitations of this study include the estimation of DNAm markers in blood samples rather than renal tissue, although there is currently no epigenetic age predictor derived in kidney tissue. Age-related methylation changes can be tissue-specific [76] and show inter-individual variation [46] , yet associations between eGFR and DNAm in blood have been demonstrated to be relevant to kidney traits [46, 77] . Moreover, this remains the only viable approach for research conducted in population-based studies, where taking renal biopsies from participants is not done due to practical and ethical considerations. Despite the bias inherent to the calculation of eGFR using equations that systematically produce higher values for individuals identified as black [2, 78, 79] , associations with trans-ethnic replication in our study featured lower eGFR (CKD). Nevertheless, future kidney research would benefit from the development and use of methods relying on filtration markers independent from muscle mass and/or moving beyond race as a variable [78] . The lack of trans-ethnic replication of all associations may be explained by multiple factors, most notably the smaller sample sizes from non-European studies, or ethnic biases in DNAm-based predictors (as those reported for PhenoAge in [80] ). Future studies using larger, homogeneous sample sizes from diverse ethnic groups are needed to address the generalizability of our findings, and to interrogate the contributions of environmental and social determinants of health disparities in epigenetic aging. Our models assumed a linear relationship in the age range studied here, and residual confounding after adjustment for the covariates included in our regression analyses is a possibility. Another potential limitation is that the cross-sectional nature of the study does not allow to draw conclusions on temporal relationships between DNAm and renal phenotypes, although DNAm patterns reflect lifetime environmental exposures and genetic factors. Our findings do not provide a mechanistic or causal explanation for renal aging and blood epigenetic aging markers, but should be considered hypothesis-generating research. Preliminary results from the largest genome-wide association (GWAS) study of DNAm-based aging and lifespan predictors found no evidence of causal effects on renal outcomes (uACR, eGFR, albuminuria and serum urate, among 150 studied traits) [81] . Nevertheless, DNAm-based aging/lifespan signatures could still be a valuable biomarker of kidney disease prognosis, risk stratification or kidney-related outcomes. Future research should aim to expand our understanding of epigenetic aging in chronic diseases and on the clinical utility of the DNAm-based predictors.

In this study of multi-ethnic population-based cohorts, kidney traits were robustly associated with DNAmbased aging and lifespan predictors measured in whole blood, as well as with some secondary DNAm-estimated markers. Our findings are consistent with a body of literature on the role immunosenescence, inflammaging and oxidative stress play in renal function and damage, as well as offer evidence on the relevance of cell-intrinsic aging mechanisms. DNAm age and lifespan predictors seem to capture the contribution of multiple CpGs to pathological changes common to systemic inflammation and renal disease, highlighting the systemic nature of age-related physiological functional decline. Future research in longitudinal studies is required to evaluate the translational value of our findings as either prognostic biomarkers for disease progression and mortality, or as means to enhance risk stratification; functional studies to explore the complex physiological interplay between epigenetic mechanisms and biological aging are also warranted.

The association between kidney traits as dependent variables and DNAm aging/lifespan predictors as independent variables was modeled using linear regression following a meta-analytic approach (Fig. 1) . Study-level results from up to seven studies were included in the meta-analyses: four studies with participants of European ancestry, one study of African American participants and three substudies from the WHI with European American, African American and Hispanic/Latino participants. The studies were KORA (Kooperative Gesundheitsforschung in der Region Augsburg), NAS (Normative Aging Study), ESTHER (Epidemiologische Studie zu Chancen der Verhütung, Früherkennung und optimierten THerapie chronischer ERkrankungen in der älteren Bevölkerung), WHI (Women's Health Initiative) and the Jackson Heart Study (JSH). The ESTHER, NAS and WHI studies contributed multiple sets of data that were analyzed separately: data sets from ESTHER corresponded to two surveys with non-overlapping sets of participants, whereas data from NAS were longitudinal and collected over consecutive examinations and analyzed taking into account these repeated measures. Further details on the data collection and methods used in each study are available in Additional file 3.

Serum creatinine values obtained with a Jaffé assay before 2009 were calibrated by multiplying by 0.95 [82] and used to calculate estimated glomerular filtration rate (eGFR) as per the CKD-EPI equation [2] in its implementation in the R package nephro [83] . Prevalent chronic kidney disease (CKD) was defined as eGFR < 60 ml/min/1.73 m 2 [84] . eGFR and urinary albumin-creatinine ratio (UACR) were log transformed prior to statistical analysis. Microalbuminuria was defined as uACR ≥ 30 mg/g. Serum urate was also studied.

Methylation was measured using the Illumina Infinium HumanMethylation450K or EPIC array in whole blood and used to estimate measures of DNAmAge and mortality (additional details on each predictor are described in Additional file 3: Note S2). Five DNAmAge and lifespan predictors were calculated using the online DNAm Age calculator (https:// dnama ge. genet ics. ucla. edu/) [11] : Hannum's estimate (HannumAge), ExtrinsicAge (EEAA) [10] , Horvath's estimate (HorvathAge) [11] , PhenoAge [18] and GrimAge [19] . IEAA, a marker capturing cellintrinsic aging properties that are independent of blood cell types, was derived by regressing HorvathAge on cell counts [10] . Quality control was conducted as in previous meta-analyses of epigenetic measures [16] , with exclusion of individuals with mismatching predicted and reported sex data. Age acceleration (AA) measures were calculated in each study as the difference between the predicted DNAmAge and chronological age, with chronological age included in all models as an adjustment for known chronological age effects across the lifespan. Defining age acceleration as the difference, rather than the residual of chronological age regressed on epigenetic age, has advantages as it is an individual measure as opposed to a population measure, and is not defined to have mean 0 in each population as is the case for the residual measure. A sixth measure, the 10-CpG-based epigenetic mortality risk score (MRS) in its continuous form was calculated as the sum of the methylation β values multiplied by the regression coefficients of each of the ten CpGs for all-cause mortality, as described in [20] .

Further measures of epigenetic aging were included in the secondary analysis: the risk-level MRS variable was built based on the cumulative number of "aberrantly" methylated CpG sites, defined by the cut-offs derived from the 4th quartile of the CpG positively correlated with mortality (cg08362785) and the 1st quartile of the other nine loci defined in [20] . Participants were then assigned to one of three risk levels based on the total number of "aberrantly" methylated CpGs: low risk, MRS = 0-1; moderate risk, MRS = 2-5; and high risk, MRS > 5. Finally, to better understand GrimAge, we also included in the analysis its eight underlying traits (smoking pack-years, adrenomedullin, beta-2 microglobulin, cystatin C, growth differentiation factor 15, leptin, plasminogen activation inhibitor 1, tissue inhibitor metalloproteinase 1) [19] if a renal phenotype was associated with GrimAge.

Linear and logistic regression models were run with kidney traits as outcomes and measures of DNAmbased age acceleration and lifespan as predictors, including covariates to adjust for potential confounding by biological and technical factors. Chronological age and sex were included in a basic model, whereas additional adjustment for BMI (kg/m 2 ), log transformed triglycerides, HDL, hypertension, smoking status (current/ever, never) and diabetes was done in the fully adjusted model. Linear regression models for serum urate additionally adjusted for baseline eGFR. Details on the study-specific definition or exclusion of variables, as well as additional information on all of the cohorts, are given in Additional file 3: Note S3. Leukocyte count (either measured or estimated by the Houseman approach [85] ) was additionally included in the regression models for Horvath's estimate as to obtain the intrinsic age acceleration measure (IEAA). In the secondary analyses, the aforementioned covariates from the basic and the fully adjusted models were used, with the exception of no smoking adjustment for DNAm-predicted pack years. All measures of association between epigenetic markers and continuous renal traits (eGFR, uACR, urate) were standardized to the standard deviation of the given renal trait as to obtain estimates comparable across renal traits.

All outcomes were available in at least one cohort of European ancestry and African American studies (Additional file 1: Table S1 ), although eGFR and CKD were the only outcomes reported by all participating studies. Each cohort provided regression estimates and standard errors, which were pooled using inverse-variance fixed-effects and random-effects models. Between-study heterogeneity was assessed using Cochran's Q and I 2 statistics. High heterogeneity was defined as I 2 > 0.50 and (Q) p het < 0.05. If high heterogeneity was detected in the fixed-effects model, the random-effects model was interpreted. All statistical analyses were conducted using R version 3.5.3 [86] , where meta-analyses were conducted using the metafor package v2.0 [87] . Multiple testing was addressed by correcting the significance level for the total number of statistical tests (i.e., Bonferroni correction, 0.05/7 epigenetic markers * 5 renal traits). Associations were considered significant if p < 1.43 E−03 in the transethnic meta-analysis and if they had consistent direction of effect across studies. Ethnic-specific replication was defined as associations with consistent direction of effect reaching nominal statistical significance (p < 0.05) in two or more studies. 

The KORA study was initiated and financed by the Helmholtz Zentrum München-German Research Center for Environmental Health, which is funded by the German Federal Ministry of Education and Research (BMBF) and by the State of Bavaria. Furthermore, KORA research has been supported within the Munich Center of Health Sciences (MC-Health), Ludwig-Maximilians-Universität, as part of LMUinnovativ. The ESTHER study was supported by the Baden-Württemberg State Ministry of Science, Research and Arts (Stuttgart, Germany), the Federal Ministry of Education and Research (Berlin, Germany), and the Federal Ministry of Family Affairs, Senior Citizens, Women and Youth (Berlin, Germany). The sponsors had no role in the study design, in the collection, analysis, and interpretation of data and preparation, review, or approval of the manuscript. The Normative Aging Study is supported by the National Institute of Environmental Health Sciences (grants P30ES009089, R01ES021733, R01ES025225, and R01ES027747). The VA Normative Aging Study is supported by the Cooperative Studies Program/Epidemiology Research and Information Center of the U.S. Department of Veterans Affairs and is a component of the Massachusetts Veterans Epidemiology Research and Information Center, Boston, Massachusetts. The WHI program is funded by the National Heart, Lung and Blood Institute, National Institutes of Health, U.S. Department of Health and Human Services through contracts HHSN268201600018C, HHSN268201600001C, HHSN268201600002C, HHSN268201600003C and HHSN268201600004C. The authors thank the WHI investigators and staff for their dedication, and the study participants for making the program possible. A full listing of WHI investigators can be found at: http:// www. whi. org/ resea rchers/ Docum ents% 20% 20Wri te% 20a% 20Pap er/ WHI% 20Inv estig ator% 20Long% 20List. pdf. This study was supported by the National Institutes of Health R01-MD012765, R01-DK117445 and R21-HL140385 to NF; by NIH/NIEHS R01-ES020836 to LH, AB and EAW; and by NIH/NHLBI contract 60442456 BAA23 to SH. The Jackson Heart Study (JHS) is supported and conducted in collaboration with Jackson State University (HHSN268201800013I), Tougaloo College (HHSN268201800014I), the Mississippi State Department of Health (HHSN268201800015I) and the University of Mississippi Medical Center (HHSN268201800010I, HHSN268201800011I and HHSN268201800012I) contracts from the National Heart, Lung and Blood Institute (NHLBI) and the National Institute on Minority Health and Health Disparities (NIMHD). The authors also wish to thank the staffs and participants of the JHS. The views expressed in this manuscript are those of the authors and do not necessarily represent the views of the National Heart, Lung and Blood Institute; the National Institutes of Health; or the US Department of Health and Human Services. The funders had no role in the design and conduct of the study, in the collection, analysis and interpretation of the data, and in the preparation, review or approval of the manuscript. The project described was supported by the National Center for Advancing Translational Sciences, National Institutes of Health, through Grant KL2TR002490 (LMR). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH. LMR was also funded by T32 HL129982 and R21-HL140385.

The dataset(s) supporting the conclusions of this article is(are) included within the article (and its additional file(s)). The informed consent given by the study participants does not cover posting of participant level phenotype data in public databases. Pre-existing data access policies for each of the five studies state research data requests can be submitted to each steering committee. Study-specific details regarding such requests are described in the next points: KORA data are available upon request from KORA Project Application Self-Service Tool (https:// epi. helmh oltz-muenc hen. de/); data requests can be submitted online and are subject to approval by the KORA Board. ESTHER data are not allowed to be publicly available due to restrictions of informed consent. However, the use of the data for collaboration projects has been and will remain the approach for data sharing. Data used in this analysis were produced and used in accordance with the policies of the Jackson Heart Study under contracts from the National, Heart, Lung and Blood Institute and are not the domain of the authors but that of the Jackson Heart Study. These data are available to other researchers for purposes of reproducing the results or replicating the procedures by submitting a manuscript proposal to the Jackson Heart Study at jhspub@umc.edu, as described at https:// www. jacks onhea rtstu dy. org/ Resea rch/ Publi catio ns# submi tmanu script. Data updates for the Jackson Heart Study are also deposited regularly in the National Institutes of Health data repositories, dbGaP (https:// www. ncbi. nlm. nih. gov/ gap/) and BioLincc (https:// bioli ncc. nhlbi. nih. gov/ home/). NAS data are available on request due to privacy/ethical restrictions. Data that support the findings of this study are available from AAB upon reasonable request. The DNA methylation datasets for WHI are publicly available under dbGAP access number phs000200.v10.p3 or upon request to www. whi. org. WHI datasets are also available through BioLINCC https:// bioli ncc. nhlbi. nih. gov/ studi es/ whi_ ctos/. Scripts used in data processing and statistical analyses have been made publicly accessible at https:// ascgi tlab. helmh oltz-muenc hen. de/ pamela. matias/ dnam-aging-in-kora-f4 and https:// ascgi tlab. helmh oltz-muenc hen. de/ pamela. matias/ kidney-dnama ge.

Ethics approval and consent to participate All participants provided written consent, and study-specific details are described next: ESTHER: The study was approved by the ethics committees of the University of Heidelberg and of the Medical Association of Saarland. All participants provided written informed consent. JHS: All participants included in this analysis provided written, informed consent for use of genetic data, and all study protocols conform to the 1975 Declaration of Helsinki guidelines. The study was approved by the Institutional Review Boards of the participating institutions (University of Mississippi Medical Center, Jackson State University and Tougaloo College). KORA: The KORA cohort ethical approval was granted by the ethics committee of the Bavarian Medical Association (REC reference numbers: F4: #06068) and all were carried out in accordance with the principles of the Declaration of Helsinki. All research participants have signed informed consent prior to taking part in any research activities. The KORA data protection procedures were approved by the responsible data protection officer of the Helmholtz Zentrum München. NAS: The NAS was approved by the Department of Veterans Affairs Boston Healthcare System and written informed consent was obtained from each subject before participation. WHI: All study participants have provided written consent to participate in genetic studies.

GFR estimation: from physiology to public health

A new equation to estimate glomerular filtration rate

Evolving importance of kidney disease: from subspecialty to global health burden

Global prevalence of chronic kidney disease: a systematic review and meta-analysis

Chronic kidney disease as a global public health problem: approaches and initiatives: a position statement from Kidney Disease Improving Global Outcomes

Association of estimated glomerular filtration rate and albuminuria with all-cause and cardiovascular mortality in general population cohorts: a collaborative meta-analysis

Lower estimated glomerular filtration rate and higher albuminuria are associated with mortality and end-stage renal disease. A collaborative meta-analysis of kidney disease population cohorts

Asymptomatic hyperuricaemia: a silent activator of the innate immune system

DNA methylation biomarkers in aging and age-related diseases

Genome-wide methylation profiles reveal quantitative views of human aging rates

DNA methylation age of human tissues and cell types

DNA methylation-based biomarkers and the epigenetic clock theory of ageing

DNA methylation age-environmental influences, health impacts, and its role in environmental epidemiology

An epigenetic clock analysis of race/ethnicity, sex, and coronary heart disease

Epigenetic clock analysis of diet, exercise, education, and lifestyle factors

DNA methylation-based measures of biological age: meta-analysis predicting time to death

An epigenetic biomarker of aging for lifespan and healthspan

DNA methylation GrimAge strongly predicts lifespan and healthspan

DNA methylation signatures in peripheral blood strongly predict all-cause mortality

Comparative validation of an epigenetic mortality risk score with three aging biomarkers for predicting mortality risks among older adult males

The epigenetic clock as a predictor of disease and mortality risk: a systematic review and meta-analysis

DNA methylation age of blood predicts all-cause mortality in later life

Epigenetic age acceleration predicts cancer, cardiovascular, and all-cause mortality in a German case cohort

Frailty is associated with the epigenetic clock but not with telomere length in a German cohort

The epigenetic clock is correlated with physical and cognitive fitness in the Lothian Birth Cohort

A systematic review and meta-analysis of environmental, lifestyle, and health factors associated with DNA methylation age

Obesity accelerates epigenetic aging of human liver

Lifetime stress accelerates epigenetic aging in an urban, African American cohort: relevance of glucocorticoid signaling

The role of epigenetics in renal ageing

Aging and DNA methylation

Epigenetic signatures of cigarette smoking

The dynamics of smoking-related disturbed methylation: a two time-point study of methylation change in smokers, non-smokers and former smokers

Smokingrelated DNA methylation is associated with DNA methylation phenotypic age acceleration: the veterans affairs normative aging study

Relationship of tobacco smoking and smoking-related DNA methylation with epigenetic age acceleration

Oxidative stress and epigenetic mortality risk score: associations with all-cause mortality among elderly people

Epigenetic predictors of all-cause mortality are associated with objective measures of neighborhood disadvantage in an urban population

Chronic inflammation (inflammaging) and its potential contribution to age-associated diseases

Biomarkers of aging and lung function in the normative aging study

Chronic kidney disease and premature ageing

KLOTHO methylation is linked to uremic toxins and chronic kidney disease

TNFalpha increases albumin permeability of isolated rat glomeruli through the generation of superoxide

Why does COVID-19 disproportionately affect older people? Aging

The role of the immune system in kidney disease

Inflamm-aging: an evolutionary perspective on immunosenescence

Ticktock chimes the kidney clock-from biology of renal ageing to clinical applications

A transcriptional profile of aging in the human kidney

Uraemia-induced immune senescence and clinical outcomes in chronic kidney disease patients

Telomere attrition is associated with inflammation, low fetuin-A levels and high mortality in prevalent haemodialysis patients

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

Age affects outcomes in chronic kidney disease

Intrinsic and extrinsic epigenetic age acceleration are associated with hypertensive target organ damage in older African Americans

GWAS of epigenetic aging rates in blood reveals a critical role for TERT

A meta-analysis of genome-wide association studies of epigenetic age acceleration

Glucocorticoids and renal Na(+) transport: implications for hypertension and salt sensitivity

Metabolic and inflammatory biomarkers are associated with epigenetic aging acceleration estimates in the GOLDN study

Investigation of bidirectional longitudinal associations between advanced epigenetic age and peripheral biomarkers of inflammation and metabolic syndrome

Changes in albuminuria but not GFR are associated with early changes in kidney structure in type 2 diabetes

Socioeconomic position, lifestyle habits and biomarkers of epigenetic aging: a multi-cohort analysis

Association of adult lung function with accelerated biological aging

Education and lifestyle factors are associated with dna methylation clocks in older African Americans

The role of epigenetic aging in education and racial/ethnic mortality disparities among older U

Underlying features of epigenetic aging clocks in vivo and in vitro

Adrenomedullin as a renal regulator peptide

Emerging risk factors and markers of chronic kidney disease progression

PAI-1 and kidney fibrosis

PAI-1 and atherothrombosis

Plasminogen activator inhibitor-1 in chronic kidney disease: evidence and mechanisms of action

Pro-A-type natriuretic peptide and pro-adrenomedullin predict progression of chronic kidney disease: the MMKD study

Aging and plasminogen activator inhibitor-1 (PAI-1) regulation: implication in the pathogenesis of thrombotic disorders in the elderly

convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year submit your research ?

Plasminogen activator inhibitor-1 and the kidney

Cardiovascular-renal complications and the possible role of plasminogen activator inhibitor: a review

Haemostasis in chronic kidney disease

Cigarette smoking and chronic kidney disease in African Americans in the Jackson heart study

Smoking: a risk factor for progression of chronic kidney disease and for cardiovascular morbidity and mortality in renal patientsabsence of evidence or evidence of absence?

DNA methylation age is associated with an altered hemostatic profile in a multi-ethnic meta-analysis

Age-related DNA methylation changes are tissue-specific with ELOVL2 promoter methylation as exception

Epigenome-wide association studies identify DNA methylation associated with kidney function

Reconsidering the consequences of using race to estimate kidney function

Hidden in plain sight-reconsidering the use of race correction in clinical algorithms

Array-based epigenetic aging indices may be racially biased

Genome-wide association studies identify 137 loci for DNA methylation biomarkers of ageing

Decline in estimated glomerular filtration rate and subsequent risk of end-stage renal disease and mortalityrisk of end-stage renal disease and mortalityrisk of end-stage renal disease and mortality

Estimating the glomerular filtration rate in the general population using different equations: effects on classification and association

Chronic kidney disease: global dimension and perspectives

DNA methylation arrays as surrogate measures of cell mixture distribution

R: a language and environment for statistical computing

Conducting meta-analyses in R with the metafor package

Publisher's Note

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

We are grateful to all study participants of ESTHER, JHS, KORA, NAS and WHI for their invaluable contributions to these studies, as well as all members of field staff conducting the studies.

The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s13148-021-01082-w.Additional file 1: Tables S1-S5 showing the results from all conducted analyses across renal parameters and DNAm-based predictors, providing both estimates from individual studies and meta-analyses.Additional file 2: Figures S1-S7 showing forest plots for all significant associations not presented in the main figures, as well as secondary traits (categorical MRS), comparison of MRS and DNAmAA effect sizes across renal traits.Additional file 3: Notes including additional information on sensitivity analyses, the derivation of the DNAm-based predictors, cohort-specific details (background and study design, as well as data collection and preprocessing), and all legends to supplemental tables and figures.