Skip to main content

Growth differentiation factor-15 and incident chronic kidney disease: a population-based cohort study

Abstract

Background

The relationship between growth differentiation factor 15 (GDF-15) and the development of chronic kidney disease (CKD) is still unclear. We sought to examine whether plasma GDF-15 was related to incident CKD and kidney function decline using a large prospective cohort study.

Methods

4318 participants of the Malmö Diet and Cancer Study-Cardiovascular Cohort were examined in 1991-1994. Incidence of CKD was followed prospectively by linkage with national patient registers. Estimated glomerular filtration rate (eGFR) was available for all participants at baseline, and was re-measured in a subgroup of 2744 subjects after 16.6 ± 1.49 years. Incidence of CKD was examined in relation to GDF-15 using Cox regression analysis. Logistic regression was used to examine the association of GDF-15 with eGFR change and eGFR-based CKD. Models were carefully corrected for potential confounders including baseline eGFR, N-terminal pro-B-type natriuretic peptide, and competing risk from death.

Results

165 patients developed CKD after 19.2 ± 4.04 years of follow-up. The adjusted hazard ratio (95% confidence interval, CI) for CKD in 4th versus 1st quartile of GDF-15 was 2.37 (1.33, 4.24) (p for trend < 0.01). Each per 1 standard deviation increase in GDF-15 was associated with a decline in eGFR of − 0.97 mL/min/1.73 m2 (95% CI, − 1.49 ~ − 0.45; p < 0.001). GDF-15 was also significantly associated eGFR-based CKD in 2713 subjects with baseline eGFR ≥60 mL/min/1.73 m2.

Conclusions

GDF-15 predicted incidence of CKD and eGFR decline in the general population, independent of a wide range of potential risk factors and competing risk of death.

Peer Review reports

Introduction

Around 9.1% of the global population are suffering from chronic kidney disease (CKD), accounting for 35.8 million disability-adjusted life-years and 1.2 million deaths annually [1]. Treatment of CKD can be costly, in particular for patients with end-stage renal disease (ESRD), which poses a considerable financial burden to families and health systems [2]. Fortunately, CKD is preventable with early detection and timely intervention [3, 4]. Efforts have thus been made to develop efficient screening strategies for CKD. In this regard, biomarkers have drawn increasing attention as they may not only help identify high-risk individuals but also provide insights into mechanism of kidney injury.

CKD is a major risk factor for CVD [1, 5], and vice versa, cardiac dysfunction could also lead to kidney injury [6]. Due to the interdependency of heart and kidney, several studies have explored kidney disease in relation to cardiovascular biomarkers, among which the predictive value of growth differentiation factor 15 (GDF-15) was recognized [7,8,9,10,11,12,13]. Most available evidence has focused on patients with existing kidney pathology [9,10,11,12]. For instance, GDF-15 was found to predict estimated glomerular filtration rate (eGFR) decline and mortality in type 1 diabetic patients with nephropathy [9], mortality in ESRD [10], eGFR decline and progression to ESRD in CKD [11], and progression to dialysis and mortality in light chain amyloidosis [12]. Cross-sectionally, GDF-15 was negatively associated with eGFR and was higher in the elderly with than without CKD [13]. However, community-based data regarding kidney function decline in relation to GDF-15 are scarce. To our knowledge, there were only two relevant studies. Participants in the study by Carlsson, et al. [8] were limited to elderly people. The authors demonstrated that GDF-15 did not predict decline of eGFR independently of baseline eGFR. In contrast, Ho, et al. [7] reported a positive association of GDF-15 with incident CKD, but some potential confounders (smoking, obesity, C-reactive protein, etc. [14]) were not adjusted for. In addition, incident CKD was identified by calculation of eGFR < 60 mL/min/1.73 m2, and no information about clinical diagnoses were included.

Therefore, we aimed to investigate the association of GDF-15 with incident CKD, as obtained from national registers, as well as from eGFR calculation. Analyses were performed in a prospective study with a large general population sample and a long-term follow-up, while taking into account potential baseline confounders and competing risk from death.

Methods

The Malmö diet and Cancer cardiovascular (MDC-CV) cohort study

The MDCS is a large prospective cohort study with participants recruited from Malmö, a city in southern Sweden [15]. During 1991-1994, a random sample of 6103 participants was taken from MDCS to investigate the epidemiology of carotid artery atherosclerosis (MDC-CV cohort study) [14]. Among them, 5540 donated fasting blood samples. We excluded participants with missing baseline data on eGFR (N = 302), waist circumference (N = 10), smoking (N = 281), low density lipoprotein (LDL, N = 2), C-reactive protein (CRP, N = 225), or GDF-15 (N = 395), or participants with previously diagnosed CKD (N = 6) or lost to follow-up (N = 1). Therefore, 4318 participants (Fig. 1, analysis 1, mean aged 57.5 ± 5.95 years, male 39.4%) remained for analyses of incident CKD, as detected by national registers of hospital inpatients and outpatients [16].

Fig. 1
figure 1

Population flow chart

During 2007-2012, MDC-CV participants who were still alive and living in the Malmö area were invited to a re-examination. A total of 3734 attended, which corresponds to 75.8% of the eligible population [17]. Among the 4318 individuals in this study, 2827 attended re-examination and 2744 had follow-up data available for eGFR. This sub-cohort study was then analyzed for decline in eGFR (Fig. 1, analysis 2). Incident CKD based on eGFR was further analyzed as the outcome in 2713 participants with baseline eGFR ≥60 mL/min/1.73 m2 (Fig. 1, analysis 3). Written informed consent was obtained from all included participants. The study conformed to the Declaration of Helsinki and was approved by the ethical committee at Lund University, Lund, Sweden (LU 51/90).

GDF-15 measurement

Fasting blood samples were collected from the cubital vein and stored at − 80 °C until assay. GDF-15 levels were measured by the SciLifeLab analysis service (Uppsala, Sweden) using Proseek® Multiplex CVD I96 × 96 reagent kit where a Proximity Extension Assay technology was applied [14, 18]. Briefly, the assay procedure consisted of three key steps: incubation, extension and detection. Raw Proseek data went through a pre-processing normalization procedure and were set relative to a fixed background level, after which Normalized Protein Expression (log2 scale) values were generated, measured in arbitrary units (AU). High AU values corresponded to a high protein concentration. GDF-15 levels in 987 subjects measured by Proseek assay closely correlated (r = 0.89 [L. Lind, unpublished data]) with the values by an electrochemiluminescence immunoassay (Roche Diagnostics, Mannheim, Germany) [14].

CKD based on the ICD codes from the national register

Information on CKD diagnosis was obtained from the Swedish patient register with nation-wide coverage. Moreover, the Swedish renal registry was searched for any additional cases of CKD [19]. CKD was defined as codes 585-586 according to ICD-9, and N18 and N19 according to ICD-10. All participants without any previous diagnosis of CKD were followed from baseline until the occurrence of a diagnosis of CKD (registry-based CKD), emigration from Sweden, death or December 31st, 2013, whichever came first. Participants were considered to have develop acute kidney injury (AKI) if they developed codes 580, 584 or 590 according to ICD-9, or N00, N10 or N17 according to ICD-10.

The CKD diagnosis in the Swedish patient register has been previously described and validated [20]. Briefly, for validation, CKD diagnoses were evaluated by two experienced specialists in nephrology. Patient records and laboratory data were reviewed and CKD cases were defined following the 2012 KDIGO criteria [21]. Validation showed that 94% of patients had a correct diagnosis of CKD [20].

CKD based on eGFR, and eGFR decline from baseline to follow-up

eGFR at baseline and follow-up was determined from plasma creatinine and cystatin C using the CKD-Epidemiology Collaboration 2012 eq. [22]. Single measurements of eGFR were assessed at each time point. A cut-off value of 60 mL/min/1.73 m2 was used to identify participants with eGFR-based CKD [20]. The difference between these two measurements was defined as eGFR change.

At baseline, creatinine and cystatin C were analyzed with the Jaffé method (Beckman Synchron LX20-4; Beckman-Coulter) and with a particle-enhanced immunonephelometry assay (N Latex Cystatin; Dade Behring, Deerfield, IL), respectively. Since the world calibrator was not introduced until 2010, cystatin C values were not standardized (reference value: 0.53 ~ 0.95 mg/L). During 2007-2012, creatinine was determined in follow-up samples using an enzymatic method (Cobas autoanalyzer; Roche Diagnostics) calibrated by isotope-dilution mass spectrometry-traceable (IDMS) creatinine [23], and cystatin C was analyzed using an automated particle-based immunoassay, adjusted to the international reference preparation ERM-DA 71/IFCC.38 T [24]. Therefore, values of creatinine and cystatin C could not be directly compared between baseline and follow-up.

Other variables and definitions

Baseline characteristics were obtained from self-administered questionnaire, physical examination, and blood measurements. Data on medication, history of kidney stone, smoking habits and alcohol consumption were collected by questionnaires. Participants were classified into current smokers, former smokers and never smokers. An average daily alcohol consumption > 40 g for males or > 30 g for females was considered as high alcohol consumption. Waist circumference was determined as being midpoint between the end of the 12th rib and the iliac crest. Blood pressure was measured with a mercury-column sphygmomanometer after 10 min of rest while the subject was in a supine position. Participant with a history of coronary event or stroke was considered to have CVD at baseline.

Glucose concentration was measured in fresh whole blood samples after an overnight fasting, following standard procedures at the Department of Clinical Chemistry, University Hospital Malmö. Diabetes was defined as self-reported physician diagnosis of diabetes, use of anti-diabetic drugs or fasting whole blood glucose ≥6.1 mmol/L (corresponding to plasma glucose ≥7.0 mmol/L). LDL concentration was estimated using the Friedewald’s formula. Measurements of biomarkers were conducted later using frozen (− 80 °C) plasma samples. CRP was measured with a Tina-quant® CRP latex assay (Roche Diagnostics, Basel, Switzerland). Methods to measure N-terminal pro-B-type natriuretic peptide (NT-proBNP) levels was the same way as that for GDF-15 [14, 18].

Statistical analyses

Baseline characteristics are presented for participants divided into quartiles (Q1-Q4) according to GDF-15 concentration, using sex-specific quartile limits. For skewed variables, log-transformation was performed to achieve a normal distribution. Differences across GDF-15 quartiles were examined using analysis of variance for continuous variables and logistic regression analysis for categorized variables.

Cox proportional hazard regression was used to analyze the association between baseline GDF-15 and incident CKD discovered by the national register. Hazard ratios (HRs) and 95% confidential intervals (CIs) were obtained. GDF-15 was treated both as a continuous variable (per standard deviation (SD) change) and as a categorized variable (in quartiles). In multivariate-adjusted models, potential covariates taken into consideration were age, sex, waist circumference, smoking, high alcohol consumption, systolic blood pressure, LDL, CRP, diabetes, CVD, anti-hypertensive drug medication, and baseline eGFR. Since GDF-15 has been frequently considered as a cardiovascular biomarker in recent years, NT-proBNP, a traditional cardiovascular marker was additionally adjusted for in a sensitivity analysis to explore whether the association of GDF-15 with CKD could be mediated by cardiac function. History of kidney stone was also included as a covariate in another sensitivity analysis (N = 4309). A restricted cubic spline function was incorporated into the Cox model to test for possible non-linearity, with knots placed at 20, 40, 60 and 80 percentages of GDF-15 concentration. Possible effect modifications were examined by introducing an interaction term between GDF-15 levels and risk factors into the multivariate model one by one. The competing risks of death was accounted for in a sensitivity analysis by the Fine and Gray proportional subdistribution hazards models method. In another sensitivity analysis, the association between GDF-15 and CKD was analyzed while participants with baseline eGFR < 60 mL/min/1.73 m2 were excluded. In addition, for participants with follow-up data available for eGFR, multiple linear regression was used to analyze the association between GDF-15 and eGFR change from baseline to the end of the follow-up. A multiple logistic regression analysis was conducted for the association between GDF-15 and eGFR-based CKD. Since CKD can be a consequence of AKI, the associations of GDF-15 with both registry- and eGFR- based CKD were re-estimated after excluding those with incident AKI before the development of CKD.

All analyses were performed using the Statistical Analysis System version 9.3 for Windows (SAS Institute Inc., Cary, NC, USA). A 2-tailed p < 0.05 was considered statistically significant.

Results

Baseline characteristics

The mean GDF-15 concentration in the cohort was 8.75 ± 0.56 AU and mean eGFR was 89.2 ± 0.56 mL/min/1.73 m2. The clinical and biochemical characteristics of the population across GDF-15 quartiles are shown in Table 1. As compared to participants with relatively low GDF-15 concentration, those with higher GDF-15 concentration tended to have decreased eGFR at baseline. An increasing trend was observed for most of the other covariates, except for sex and high alcohol consumption.

Table 1 Characteristics of individuals across quartiles (Q1-Q4) of growth differentiation factor-15 (GDF-15)

Incidence of register-based CKD in relation to GDF-15 (analysis 1)

During a mean of 19.2 ± 4.04 years of follow-up, a total of 165 subjects developed CKD. After multivariate adjustment (Table 2, Model 3), an increased GDF-15 level was associated with a higher risk of developing CKD. The HR (highest vs. lowest quartiles of GDF-15) for incident CKD was 2.74 (95% CI, 1.53 ~ 4.89; p for trend < 0.001). This value was slightly attenuated after additionally adjusting for baseline eGFR (HR, 2.37; 95% CI, 1.33 ~ 4.24; p for trend < 0.01). The adjusted HR for each 1 SD increase in GDF-15 was 1.39 (95% CI, 1.16 ~ 1.65; p < 0.001). In the sensitivity analysis when NT-proBNP was additionally added into the model, the association between GDF-15 and CKD hardly changed (HR: 1.35, 95% CI, 1.13 ~ 1.62; p < 0.01, per 1 SD increment of GDF-15). However, no statistical significance was observed for NT-proBNP (p = 0.10).

Table 2 Incident cases of register-based chronic kidney disease in relation to growth differentiation factor-15 (GDF-15)

No obvious evidence of non-linearity in the association between GDF-15 and CKD was detected by restricted cubic spline function (p for effect test < 0.0001, p for non-linearity test = 0.26). Meanwhile, no interaction between GDF-15 and covariates was found with respect to CKD. During the follow-up period, 952 individuals died from causes other than CKD. When competing risk of death was taken in consideration, the adjusted HR was 2.11 (95% CI, 1.19 ~ 3.76) for Q4 versus Q1 of GDF-15 (p for trend =0.01), and was 1.23 (95% CI, 1.03 ~ 1.48) for each 1 SD increase in GDF-15 (data not shown). After adjusting for history of kidney stone (395 cases out of 4309 subjects), the adjusted HR was 2.40 (95% CI, 1.34 ~ 4.29) for Q4 versus Q1 of GDF-15 (p for trend < 0.01), and was 1.39 (95% CI, 1.16 ~ 1.65) for each 1 SD increase in GDF-15 (data not shown). Among 4244 individuals with baseline eGFR ≥60 mL/min/1.73 m2, 145 developed CKD. The adjusted HR was 2.08 (95% CI, 1.14 ~ 3.78; p for trend =0.02), and 1.22 (95% CI, 1.01 ~ 1.48) for Q4 versus Q1 and per 1 SD increase of GDF-15, respectively (data not shown). A total of 182 subjects developed AKI during the follow-up, and 13 out of them (7.14%) were later re-classified as CKD. After excluding subjects with AKI prior to CKD, the adjusted HR was 2.32 (95% CI, 1.25 ~ 4.30; p for trend < 0.01), and 1.33 (95% CI, 1.10 ~ 1.60) for Q4 versus Q1 and per 1 SD increase of GDF-15, respectively (data not shown).

eGFR decline and incidence of eGFR-based CKD in relation to GDF-15

A total of 2744 had repeated eGFR values after a mean follow-up of 16.6 ± 1.49 years. The mean eGFR values at baseline and follow-up were 90.2 ± 12.8 and 66.2 ± 15.1 mL/min/1.73 m2, respectively (N = 2744). The association of eGFR change in relation to GDF-15, both in quartiles and per 1 SD increase, is presented in Table 3 (Analysis 2). As compared to Q1 of GDF-15, Q4 was associated with a greater eGFR decline during follow-up (− 2.42 mL/min/1.73 m2 (95% CI, − 3.91 ~ − 0.94); p for trend < 0.01) after multivariate adjustment including baseline eGFR. In addition, each 1 SD increase in GDF-15 was associated with a decline in eGFR of − 0.97 mL/min/1.73 m2 (95% CI, − 1.49 ~ − 0.45; p < 0.001) over the follow-up period. Results were consistent in participants with baseline eGFR ≥60 mL/min/1.73 m2 (N = 2713; B = -0.92; 95% CI, − 1.49 ~ − 0.36, p < 0.01) (data not shown).

Table 3 Change in estimated glomerular filtration rate in relation to baseline growth differentiation factor-15 (GDF-15)

Out of the 2713 individuals with baseline eGFR ≥60 mL/min/1.73 m2, 862 developed CKD, as defined by eGFR < 60 mL/min/1.73 m2. After adjusting for covariates including baseline eGFR, the odds ratio (OR) comparing Q4 vs Q1 of GDF-15 was 1.41 (95% CI, 1.06 ~ 1.89; p for trend =0.02) for developing eGFR-based CKD. The corresponding OR was 1.15 (95% CI, 1.04 ~ 1.28) per 1 SD increase in GDF-15 (Analysis 3: data not shown). Before re-examinations, 43 participants developed AKI based on ICD codes, of whom 26 (60.5%) were later identified to have eGFR-based CKD during re-examinations. After excluding them, the adjusted OR was 1.14 (95% CI, 1.03 ~ 1.27) for developing eGFR-based CKD per 1 SD increase in GDF-15.

Discussion

Our findings suggest that GDF-15 is associated with increased incidence of CKD and eGFR decline. The association was independent of baseline eGFR, smoking, waist circumference, CRP, etc., and remained after controlling for competing risk of death. GDF-15, therefore, may be a useful marker of increased risk of CKD.

GDF-15 is a distant member of the transforming growth factor-β (TGF-β) family [25]. It has consistently been associated with deterioration of kidney function or adverse outcomes in patients with existing kidney diseases [9,10,11,12]. GDF-15 has also been investigated in relation to the development of CKD [7, 8], yet the association has remained uncertain. Carlsson, et al. used a discovery cohort (the PIVUS Study; N = 687, mean age = 70 years) and a replication cohort (the ULSAM study; N = 360, mean age = 78 years) to identify predictors of eGFR decline from 80 CVD biomarkers [8]. Whereas GDF-15 was initially observed to be associated with eGFR decline during a 5-year follow-up period, the association disappeared after adjusting for eGFR at baseline. The validity of this study may be affected by the relatively small sample size, the elderly population and potential survival bias. Results from the Framingham Offspring cohort study (N = 2614, mean follow-up period = 9.5 years) [7] suggested a superior predictive value of GDF-15 for incident CKD (estimated by one eGFR measurement). The association of GDF-15 with CKD was independent of baseline eGFR.

Even though the exact production rate and the kidney clearance rate is currently not available in our understanding of GDF-15, as a 24.5-kDa active circulating dimeric protein [26], levels of GDF-15 may largely depend upon renal excretory function for elimination. In support of this view, GDF-15 concentration negatively correlated with eGFR in the current and previous studies [8, 13], and was decreased by kidney function improvement through kidney transplantation [27]. Nevertheless, a strong association of GDF-15 with incident CKD was observed in our study and the study by Ho, et al. [7], even after adjustment for baseline eGFR, which may abolish some of the doubts about the kidney elimination dependence. Noteworthily, some potential confounders were not taken into account in the study by Ho, et al. [7]. We confirmed previous findings [14], suggesting that smoking, obesity, hs-CRP, etc. are strongly correlated with GDF-15 levels, and these are factors that also have been associated with CKD [28,29,30]. Our finding provides added evidence that GDF-15 was significantly associated with incidence of CKD and eGFR decline even after adjustments for baseline eGFR as well as other potential risk factors.

Dysfunction of the heart or kidney could potentially induce dysfunction of the other organ [1, 5, 6], known as a pathologic condition termed the cardiorenal syndrome [31]. In our sensitivity analysis, GDF-15 and NT-proBNP were mutually adjusted for each other in the multivariate model. GDF-15 but not NT-proBNP remained significantly associated with CKD and eGFR decline, suggesting that GDF-15 might be more specific to kidney outcomes than other CVD biomarkers. Moreover, the competing risk of death is usually high in studies on geriatric populations or with long-term follow-ups, and may largely bias the findings [32, 33]. We demonstrated for the first time that the fully-adjusted association of GDF-15 with eGFR decline and CKD remained after controlling for death as a competing risk.

The kidney can be an important source of GDF-15. In kidney tissue from adult rats, GDF-15 mRNA expression was mainly detected in the S3 segment of the nephron and the collecting ducts by in situ hybridization [34]. In response to stimuli such as surgery, toxin, ischemia, and hyperoxia, GDF-15 can be immediately induced in kidney, possibly through TNF- and p53-dependent and -independent pathways, and acts as a regulator of inflammation, cell survival, proliferation, and apoptosis [35]. Duong et al., observed that when stimulated by metabolic acidosis, GDF-15 expression was strongly induced in mouse kidney outer medullary collecting duct [36]. The increased GDF-15 played a key role in collecting duct lengthening by triggering compensatory proliferation of acid-secreting intercalated cells [36]. In contrast, genetic deletion of GDF-15 aggravated tubular and interstitial injury, which resulted in glycosuria and polyuria in mice with diabetes [37]. Meanwhile, epidemiological evidence showed that circulating levels of GDF-15 were closely correlated with mRNA expression of GDF-15 in renal tubulointerstitium, and significantly predicted risk of disease progression in patients with CKD [11]. Similarly, urinary GDF-15 was tightly linked with proximal tubule damage and kidney function decline in diabetic patients [38]. Based on the above-mentioned evidence, it is speculated that the observed association between plasma GDF-15 and incident CKD might at least partly be due to enhanced GDF-15 expression in kidney as a protective response against early-stage renal damage.

Strengths of this study included a prospective study design and a long-term follow-up. In addition, endpoints were retrieved from hospital registers with national coverage and a high case validity of the CKD diagnosis, as well as based on eGFR determined from a combination of plasma creatinine and cystatin C. There were several limitations that need to be mentioned. Our findings were descriptive in nature and provided limited information on mechanisms. Urine albumin is a strong predictor of incident CKD but data on albuminuria was not available. This did not affect eGFR calculation, but unfortunately we could not examine whether GDF-15 was associated with CKD independently of albuminuria. Moreover, residual confounding could result from other known risk factors, e.g. serum uric acid, family history of CKD, viral hepatitis, autoimmune disease and use of nephrotoxic agents as well as from unknown risk factors, even though extensive adjustments have already been made. Thus, further studies are still needed to further confirm the relationship.

CKD cases identified among hospital registers may not include less severe cases which were treated in primary care or were asymptomatic, and thus tended to underestimate CKD incidence. On the contrary, eGFR-based CKD may overestimate CKD incidence by using one eGFR measurement instead of repeated testing 3 months apart [21]. As a result, huge difference existed in different measures of CKD incidence, both of which might not be exactly the same as the true incidence. However, since GDF-15 was consistently associated with eGFR decline, incidence of register-based CKD as well as eGFR-based criteria, the overall results strongly supported a relationship with clinical CKD outcomes. Another important limitation is that the calculation of eGFR change may be influenced since the methods of measuring cystatin C and creatinine at baseline were different from those at the follow-up. Nevertheless, we expect that any bias introduced by these different measurements should be non-differential in relation to GDF-15 levels, and a greater decline in eGFR from baseline to follow-up could still reflect deteriorated kidney function. Therefore, we consider that the association between GDF-15 and eGFR decline was valid.

Conclusions

In conclusion, this study indicated that elevated GDF-15 levels predicted incidence of CKD and eGFR decline in the general population, independent of a wide range of potential risk factors and competing risk of death.

Availability of data and materials

The data that support the findings of this study are available from Lund University, but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of Lund University.

Abbreviations

AU:

Arbitrary units

CI:

Confidence interval

CKD:

Chronic kidney disease

CRP:

C-reactive protein

eGFR:

Estimated glomerular filtration rate

ESRD:

End-stage renal disease

GDF-15:

Growth differentiation factor 15

HR:

Hazard ratio

IDMS:

Isotope-dilution mass spectrometry-traceable

LDL:

Low density lipoprotein

MDC-CV:

Malmö Diet and Cancer Cardiovascular

NT-proBNP:

N-terminal pro-B-type natriuretic peptide

OR:

Odds ratio

SD:

Standard deviation

TGF-β:

Transforming growth factor-β

References

  1. Collaboration GCKD. Global, regional, and national burden of chronic kidney disease, 1990-2017: a systematic analysis for the global burden of disease study 2017. Lancet. 2020;395:709–33.

    Article  Google Scholar 

  2. Himmelfarb J, Ikizler TA. Hemodialysis. N Engl J Med. 2010;363:1833–45.

    CAS  Article  Google Scholar 

  3. Stevens PE, Levin A, Members KDIGOCKDGDWG. Evaluation and management of chronic kidney disease: synopsis of the kidney disease: improving global outcomes 2012 clinical practice guideline. Ann Intern Med. 2013;158:825–30.

    Article  Google Scholar 

  4. Berns JS. Nutritional Management of Chronic Kidney Disease. N Engl J Med. 2018;378:584.

    PubMed  Google Scholar 

  5. Sarnak MJ, Levey AS, Schoolwerth AC, Coresh J, Culleton B, Hamm LL, et al. Kidney disease as a risk factor for development of cardiovascular disease: a statement from the American Heart Association councils on kidney in cardiovascular disease, high blood pressure research, clinical cardiology, and epidemiology and prevention. Circulation. 2003;108:2154–69.

    Article  Google Scholar 

  6. Damman K, Testani JM. The kidney in heart failure: an update. Eur Heart J. 2015;36:1437–44.

    Article  Google Scholar 

  7. Ho JE, Hwang SJ, Wollert KC, Larson MG, Cheng S, Kempf T, et al. Biomarkers of cardiovascular stress and incident chronic kidney disease. Clin Chem. 2013;59:1613–20.

    CAS  Article  Google Scholar 

  8. Carlsson AC, Ingelsson E, Sundstrom J, Carrero JJ, Gustafsson S, Feldreich T, et al. Use of proteomics to investigate kidney function decline over 5 years. Clin J Am Soc Nephrol. 2017;12:1226–35.

    Article  Google Scholar 

  9. Lajer M, Jorsal A, Tarnow L, Parving HH, Rossing P. Plasma growth differentiation factor-15 independently predicts all-cause and cardiovascular mortality as well as deterioration of kidney function in type 1 diabetic patients with nephropathy. Diabetes Care. 2010;33:1567–72.

    CAS  Article  Google Scholar 

  10. Breit SN, Carrero JJ, Tsai VW, Yagoutifam N, Luo W, Kuffner T, et al. Macrophage inhibitory cytokine-1 (MIC-1/GDF15) and mortality in end-stage renal disease. Nephrol Dial Transplant. 2012;27:70–5.

    CAS  Article  Google Scholar 

  11. Nair V, Robinson-Cohen C, Smith MR, Bellovich KA, Bhat ZY, Bobadilla M, et al. Growth differentiation Factor-15 and risk of CKD progression. J Am Soc Nephrol. 2017;28:2233–40.

    CAS  Article  Google Scholar 

  12. Kastritis E, Papassotiriou I, Merlini G, Milani P, Terpos E, Basset M, et al. Growth differentiation factor-15 is a new biomarker for survival and renal outcomes in light chain amyloidosis. Blood. 2018;131:1568–75.

    CAS  Article  Google Scholar 

  13. Kim JS, Kim S, Won CW, Jeong KH. Association between plasma levels of growth differentiation Factor-15 and renal function in the elderly: Korean frailty and aging cohort study. Kidney Blood Press Res. 2019;44:405–14.

    CAS  Article  Google Scholar 

  14. Bao X, Borne Y, Muhammad IF, Nilsson J, Lind L, Melander O, et al. Growth differentiation factor 15 is positively associated with incidence of diabetes mellitus: the Malmo diet and Cancer-cardiovascular cohort. Diabetologia. 2019;62:78–86.

    CAS  Article  Google Scholar 

  15. Berglund G, Elmstahl S, Janzon L, Larsson SA. The Malmo diet and Cancer study. Design and feasibility. J Intern Med. 1993;233:45–51.

    CAS  Article  Google Scholar 

  16. Ludvigsson JF, Andersson E, Ekbom A, Feychting M, Kim JL, Reuterwall C, et al. External review and validation of the Swedish national inpatient register. BMC Public Health. 2011;11:450.

    Article  Google Scholar 

  17. Rosvall M, Persson M, Ostling G, Nilsson PM, Melander O, Hedblad B, et al. Risk factors for the progression of carotid intima-media thickness over a 16-year follow-up period: the Malmo diet and Cancer study. Atherosclerosis. 2015;239:615–21.

    CAS  Article  Google Scholar 

  18. Lundberg M, Eriksson A, Tran B, Assarsson E, Fredriksson S. Homogeneous antibody-based proximity extension assays provide sensitive and specific detection of low-abundant proteins in human blood. Nucleic Acids Res. 2011;39:e102.

    CAS  Article  Google Scholar 

  19. Qureshi AR, Evans M, Stendahl M, Prutz KG, Elinder CG. The increase in renal replacement therapy (RRT) incidence has come to an end in Sweden-analysis of variations by region over the period 1991-2010. Clin Kidney J. 2013;6:352–7.

    Article  Google Scholar 

  20. Harari F, Sallsten G, Christensson A, Petkovic M, Hedblad B, Forsgard N, et al. Blood Lead levels and decreased kidney function in a population-based cohort. Am J Kidney Dis. 2018;72:381–9.

    CAS  Article  Google Scholar 

  21. Kidney Disease: Improving Global Outcomes (KDIGO) CKD Work Group. KDIGO 2012 clinical practice guideline for the evaluation and management of chronic kidney disease. Kidney Int Suppl. 2013;3:1–150.

    Article  Google Scholar 

  22. Inker LA, Schmid CH, Tighiouart H, Eckfeldt JH, Feldman HI, Greene T, et al. Estimating glomerular filtration rate from serum creatinine and cystatin C. N Engl J Med. 2012;367:20–9.

    CAS  Article  Google Scholar 

  23. Nyman U, Grubb A, Larsson A, Hansson LO, Flodin M, Nordin G, et al. The revised Lund-Malmo GFR estimating equation outperforms MDRD and CKD-EPI across GFR, age and BMI intervals in a large Swedish population. Clin Chem Lab Med. 2014;52:815–24.

    CAS  Article  Google Scholar 

  24. Grubb A, Horio M, Hansson LO, Bjork J, Nyman U, Flodin M, et al. Generation of a new cystatin C-based estimating equation for glomerular filtration rate by use of 7 assays standardized to the international calibrator. Clin Chem. 2014;60:974–86.

    CAS  Article  Google Scholar 

  25. Bootcov MR, Bauskin AR, Valenzuela SM, Moore AG, Bansal M, He XY, et al. MIC-1, a novel macrophage inhibitory cytokine, is a divergent member of the TGF-beta superfamily. Proc Natl Acad Sci U S A. 1997;94:11514–9.

    CAS  Article  Google Scholar 

  26. Desmedt S, Desmedt V, De Vos L, Delanghe JR, Speeckaert R, Speeckaert MM. Growth differentiation factor 15: a novel biomarker with high clinical potential. Crit Rev Clin Lab Sci. 2019;56:333–50.

    Article  Google Scholar 

  27. Connelly PW, Yan AT, Nash MM, Lok CE, Gunaratnam L, Prasad GVR. Growth differentiation factor 15 is decreased by kidney transplantation. Clin Biochem. 2019;73:57–61.

    CAS  Article  Google Scholar 

  28. Jones-Burton C, Seliger SL, Scherer RW, Mishra SI, Vessal G, Brown J, et al. Cigarette smoking and incident chronic kidney disease: a systematic review. Am J Nephrol. 2007;27:342–51.

    Article  Google Scholar 

  29. Garofalo C, Borrelli S, Minutolo R, Chiodini P, De Nicola L, Conte G. A systematic review and meta-analysis suggests obesity predicts onset of chronic kidney disease in the general population. Kidney Int. 2017;91:1224–35.

    Article  Google Scholar 

  30. Cachofeiro V, Goicochea M, de Vinuesa SG, Oubina P, Lahera V, Luno J. Oxidative stress and inflammation, a link between chronic kidney disease and cardiovascular disease. Kidney Int Suppl. 2008:S4–9.

  31. Wang J, Zhang W, Wu L, Mei Y, Cui S, Feng Z, et al. New insights into the pathophysiological mechanisms underlying cardiorenal syndrome. Aging (Albany NY). 2020;12:12422–31.

    CAS  Article  Google Scholar 

  32. Berry SD, Ngo L, Samelson EJ, Kiel DP. Competing risk of death: an important consideration in studies of older adults. J Am Geriatr Soc. 2010;58:783–7.

    Article  Google Scholar 

  33. Ensrud KE, Harrison SL, Cauley JA, Langsetmo L, Schousboe JT, Kado DM, et al. Impact of competing risk of mortality on Association of Weight Loss with Risk of central body fractures in older men: a prospective cohort study. J Bone Miner Res. 2017;32:624–32.

    Article  Google Scholar 

  34. Unsicker K, Spittau B, Krieglstein K. The multiple facets of the TGF-beta family cytokine growth/differentiation factor-15/macrophage inhibitory cytokine-1. Cytokine Growth Factor Rev. 2013;24:373–84.

    CAS  Article  Google Scholar 

  35. Zimmers TA, Jin X, Hsiao EC, McGrath SA, Esquela AF, Koniaris LG. Growth differentiation factor-15/macrophage inhibitory cytokine-1 induction after kidney and lung injury. Shock. 2005;23:543–8.

    CAS  PubMed  Google Scholar 

  36. Duong Van Huyen JP, Cheval L, Bloch-Faure M, Belair MF, Heudes D, Bruneval P, et al. GDF15 triggers homeostatic proliferation of acid-secreting collecting duct cells. J Am Soc Nephrol. 2008;19:1965–74.

    Article  Google Scholar 

  37. Mazagova M, Buikema H, van Buiten A, Duin M, Goris M, Sandovici M, et al. Genetic deletion of growth differentiation factor 15 augments renal damage in both type 1 and type 2 models of diabetes. Am J Physiol Ren Physiol. 2013;305:F1249–64.

    CAS  Article  Google Scholar 

  38. Simonson MS, Tiktin M, Debanne SM, Rahman M, Berger B, Hricik D, et al. The renal transcriptome of db/db mice identifies putative urinary biomarker proteins in patients with type 2 diabetes: a pilot study. Am J Physiol Ren Physiol. 2012;302:F820–9.

    Article  Google Scholar 

Download references

Acknowledgments

We thank all the participants in this study and the staff of the department of Clinical Chemistry at University Hospital, Malmö. We also would like to acknowledge the Swedish Patient Registry and the Swedish Renal Registry for providing information about renal outcomes and the Clinical biomarker facility at SciLifeLab Sweden for providing assistance in protein analyses.

Funding

This work was supported by the Swedish Heart–Lung Foundation (2016-0315, 2017-0626); Lund University Infrastructure grant “Malmö population-based cohorts” (STYR 2019/2046); the Natural Science Foundation of Jiangsu Province (BK20200128); and the Fundamental Research Funds for the Central Universities (0214-14380474). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Affiliations

Authors

Contributions

XB, BX and GE contributed to the conception or design of the work. XB and YB contributed to the acquisition, analysis, or interpretation of data for the work. XB drafted the manuscript. BX, YB, MOM, OM, JN, AC and GE critically revised the manuscript. All authors reviewed and approved the final manuscript.

Corresponding authors

Correspondence to Biao Xu or Gunnar Engström.

Ethics declarations

Ethics approval and consent to participate

All procedures performed in this study were approved by the ethics committee at Lund University Lund, Sweden (LU 51/90). Written informed consent was obtained from all individual participants included in the study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

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

Verify currency and authenticity via CrossMark

Cite this article

Bao, X., Xu, B., Borné, Y. et al. Growth differentiation factor-15 and incident chronic kidney disease: a population-based cohort study. BMC Nephrol 22, 351 (2021). https://doi.org/10.1186/s12882-021-02558-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12882-021-02558-w

Keywords

  • Growth differentiation factor 15
  • Chronic kidney disease
  • Estimated glomerular filtration rate
  • Competing risk
  • Cohort study