Visual Overview
Abstract
Iron deficiency (ID) is independently associated with an increased risk of death in renal transplant recipients (RTRs). ID promotes production and cleavage of intact fibroblast growth factor 23 (iFGF23) into C-terminal fibroblast growth factor 23 (cFGF23), elevated levels of which are also prospectively associated with adverse outcomes. We hypothesized that in RTRs, the relationship between ID and mortality is mediated by FGF23. We measured plasma iFGF23 and cFGF23 levels in 700 stable RTRs at a median of 5.4 years after transplant. RTRs with ID had median (interquartile range) cFGF23 concentrations higher than those of RTRs without ID (223 [131–361] versus 124 [88–180] RU/ml; P<0.001), whereas iFGF23 concentrations were similar between groups. In multivariable-adjusted Cox regression analyses, ID associated with increased mortality (81 events; hazard ratio, 1.95; 95% confidence interval, 1.22 to 3.10; P<0.01). However, this association lost significance after additional adjustment for cFGF23 levels (hazard ratio, 1.45; 95% confidence interval, 0.87 to 2.51; P=0.15). In further mediation analysis, cFGF23 explained 46% of the association between ID and mortality, whereas iFGF23 did not mediate this association. In conclusion, we found that cFGF23 levels are increased in iron-deficient RTRs and that the underlying biologic process driving production and cleavage of iFGF23, or alternatively the increased level of cFGF23 fragments, probably is an important mediator of the association between ID and mortality. Our results underline the strong relationship between iron and FGF23 physiology, and provide a potential mechanism explaining the relationship between ID and adverse outcome in RTRs.
Iron deficiency (ID) is highly prevalent among renal transplant recipients (RTRs) and an important contributor to post-transplant anemia.1,2 In addition to its role in hemoglobin synthesis, iron also plays a pivotal role in oxygen sensing, synthesis of DNA, electron transport, and cellular proliferation and differentiation.3 Independent of anemia, ID is a known risk factor for mortality in RTRs, although the underlying mechanisms are unclear.4
Recent studies suggest that ID is crucially involved in fibroblast growth factor 23 (FGF23) production and metabolism. FGF23 is an osteocyte-derived hormone, and an essential regulator of phosphate metabolism, among others, by influencing vitamin D homeostasis. Higher plasma FGF23 levels are associated with an increased risk of mortality in RTRs,5,6 which may be explained, at least in part, by off-target effects of high FGF23 levels on the heart and other organs.7,8 Experimental animal models revealed that ID stimulates FGF23 transcription accompanied by increased intracellular cleavage of FGF23, which results in elevated circulating C-terminal FGF23 (cFGF23) concentrations, but relatively normal intact FGF23 (iFGF23) concentrations.9 Observational studies in humans have also demonstrated an inverse relationship between markers of iron status and cFGF23 levels.10–13 Moreover, administration of intravenous iron in iron-deficient anemic women resulted in markedly decreased circulating cFGF23 concentrations,14 consistent with the notion that ID is an important physiologic regulator of increased cFGF23 levels.
To date, it is unknown whether an association exists between FGF23 and ID in RTRs and whether FGF23 modulates the increased mortality risk observed in RTRs with ID. Therefore, the aim of this study was to elucidate whether ID influences FGF23 levels, and whether the association of ID with all-cause mortality is mediated by FGF23.
Results
Baseline Characteristics
We included 700 RTRs at a median of 5.4 (interquartile range [IQR], 1.9–12.0) years after transplantation. Mean age was 53±13 years; 57% of participants were men; mean body mass index (BMI) was 26.7±4.8 kg/m2. Additional baseline characteristics are shown in Table 1.
Baseline characteristics of RTRs according to ID
Median plasma iFGF23 and cFGF23 concentrations were 62 (IQR, 43–99) pg/ml and 140 (95–233) RU/ml, respectively. Mean hemoglobin concentration was 13.2±1.7 g/dl; median ferritin concentration was 118 (54–222) µg/L; and mean TSAT was 25.4%±10.8%. ID, defined as transferrin saturation (TSAT) <20% and ferritin<300 µg/L, was present in 208 (30%) patients. Significant differences in baseline characteristics between RTRs with versus without ID were noted with respect to sex, BMI, smoking status, time since transplantation, diabetes mellitus, hemoglobin concentration, mean corpuscular volume, ferritin, TSAT, proteinuria, high-sensitivity C-reactive protein (hs-CRP), and use of angiotensin converting enzyme (ACE) inhibitors and diuretics (Table 1). Increased cFGF23 levels were noted in the iron-deficient compared with non–iron-deficient RTRs (223 [131–361] versus 124 [88–180] RU/ml; P<0.001), whereas iFGF23 levels were similar between groups (Table 1).
ID, FGF23, and All-Cause Mortality
During a median follow-up of 3.1 (2.7–3.9) years, 81 (12%) RTRs died, of whom 38 (47%) died from cardiovascular causes. Other causes of death were infection (24%), malignancy (16%), and miscellaneous causes (14%).
In unadjusted Cox regression analysis, ID (hazard ratio [HR], 2.04 for present versus absent; 95% confidence interval [95% CI], 1.31 to 3.16; P<0.001), cFGF23 (HR, 1.61 per SD; 95% CI, 1.35 to 1.91; P<0.001), and iFGF23 (HR, 1.33 per SD; 95% CI, 1.11 to 1.60; P=0.002) were associated with all-cause mortality. The simulation study showed that this difference in HRs could not be explained by the differences in intra-assay coefficients of variation (CVs) (HR cFGF23, 1.60 per SD; 95% CI, 1.35 to 1.90). We observed no effect modification of the association between ID and mortality by age, sex, eGFR, proteinuria, history of cardiovascular disease, time since transplantation, smoking status, BMI, presence of diabetes, hs-CRP, serum calcium, phosphate, parathyroid hormone (PTH), use of ACE inhibitors, use of diuretics, cFGF23, and iFGF23 concentrations (Pinteraction>0.10 for all). In multivariable Cox regression analysis, the association between ID and mortality remained significant after adjustment for age, sex, eGFR, proteinuria, time since transplantation, primary renal disease, history of cardiovascular disease, and smoking status (HR, 1.95; 95% CI, 1.22 to 3.10; P<0.01). Further adjustment for iFGF23 did not affect the association between ID and mortality (HR, 1.94; 95% CI, 1.22 to 3.10; P<0.01). In contrast, adjustment for cFGF23 abolished the association between ID and mortality such that it was no longer significant (HR, 1.45; 95% CI, 0.87 to 2.51; P=0.15) (Table 2). The correlation coefficient of ID with cFGF23 levels was r=0.35, P<0.001, and the variance inflation factor (VIF) was 1.1, indicating absence of relevant colinearity.
Univariate and multivariate-adjusted associations between ID and all-cause mortality
Mediation Analysis
In mediation analyses, iFGF23 was not found to be a significant mediator (P value for indirect effect >0.05) of the effect between ID and mortality. In contrast, cFGF23 was a significant mediator (P value for indirect effect <0.05; 46% of the association between ID and mortality was explained by cFGF23) (Table 3). Differences in intra-assay CVs could not explain this difference in mediation effects either (proportion of mediation from simulated data, 45.7%).
Mediation analysis of ID with all-cause mortality through iFGF23, cFGF23, and hs-CRP
Sensitivity Analyses
In sensitivity analyses, we first assessed the prospective association between ID and mortality using an alternative definition of ID: TSAT<20% and ferritin<200 µg/L. The association between ID and mortality remained significant independent of adjustment for age, sex, eGFR, proteinuria, time since transplantation, primary renal disease, history of cardiovascular disease, and smoking status (HR, 1.87; 95% CI, 1.16 to 3.01; P=0.01). Similar to the primarily used definition of ID, the association persisted after adjustment for iFGF23 (HR, 1.87; 95% CI, 1.16 to 3.02; P=0.01), whereas adjustment for cFGF23 abolished the association (HR, 1.40; 95% CI, 0.81 to 2.40; P=0.23).
Second, we assessed whether the association between ID and mortality remained independent of adjustment for markers related to mineral metabolism other than FGF23, i.e., serum calcium and serum PTH. In these analyses, the association between ID and mortality remained materially unchanged independent of further adjustment for serum calcium in addition to adjustment for age, sex, eGFR, proteinuria, time since transplantation, primary renal disease, and smoking status (HR, 1.90; 95% CI 1.19 to 3.03; P<0.01). The same was true with further adjustment for serum PTH (HR, 1.96; 95% CI, 1.23 to 3.13; P<0.01 [see Table 2, model 2 for comparison]).
Furthermore, we assessed whether the association between ID and mortality remained independent of adjustment for ACE inhibitors and angiotensin II (AII) antagonists, hemoglobin levels, acute rejection, presence of diabetes, use of iron supplements, or hs-CRP. In these analyses, the association between ID and mortality remained materially unchanged independent of further adjustment for use of ACE inhibitors and AII antagonists (HR, 1.97; 95% CI, 1.23 to 3.16; P<0.01 [see Table 2, model 2 for comparison]), in addition to adjustment for age, sex, eGFR, proteinuria, time since transplantation, primary renal disease, and smoking status. The same was true with further adjustment for hemoglobin (HR, 1.76; 95% CI, 1.09 to 2.84; P=0.02), acute rejection (HR, 1.93; 95% CI, 1.21 to 3.09; P<0.01), presence of diabetes (HR, 1.83; 95% CI, 1.13 to 2.94; P=0.01), use of iron supplements (HR, 2.06; 95% CI, 1.29 to 3.29; P=0.003), and hs-CRP (HR, 1.86; 95% CI, 1.15 to 2.99; P=0.01).
Finally, we assessed the associations between the individual iron status components and mortality. The relationship between serum ferritin and mortality demonstrated a nonlinear relationship, as shown by cubic restricted splines (Figure 1A). When divided into quintiles and adjusted according to model 2, the lowest and the highest quintiles of serum ferritin were associated with an increased risk of mortality (HR, 2.82; 95% CI, 1.20 to 6.63; P=0.02; and HR, 2.49; 95% CI, 1.10 to 5.65; P=0.03, respectively), compared with the fourth quintile. Adjustment for iFGF23 did not materially alter the association of the lowest quintile (HR, 2.79; 95% CI, 1.19 to 6.59; P=0.02) and the highest quintile (HR, 2.54; 95% CI, 1.12 to 5.77; P=0.03) of ferritin with mortality. In contrast, after adjustment for cFGF23, the association of the lowest quintile with mortality was markedly weakened and no longer significant (HR, 1.67; 95% CI, 0.66 to 4.22; P=0.27), but the highest quintile of ferritin was still significantly associated with mortality (HR, 2.78; 95% CI, 1.22 to 6.37; P=0.02; Figure 1B).
Associations of low serum ferritin and low TSAT with all-cause mortality are abrogated by adjustment for cFGF23. Data were fit by a Cox proportional hazard regression model on the basis of restricted cubic splines. For both variables, the median was utilized as reference (4.8 µg/L for naturally logarithmic ferritin, 24% for TSAT). The black line represents the HR, the gray area the 95% CI. Panels (A and C) show univariate analyses between serum ferritin and TSAT and all-cause mortality. Panels (B and D) show adjustment for cFGF23.
TSAT was inversely associated with mortality (HR, per 5% increase, 0.86; 95% CI, 0.77 to 0.96; P<0.01; Figure 1C). In multivariable analysis, TSAT remained significantly associated with mortality (HR, 0.87; 95% CI, 0.78 to 0.98; P=0.02). Adjustment for iFGF23 did not alter the association of TSAT with mortality (HR, 0.88; 95% CI, 0.78 to 0.98; P=0.02). In contrast, further adjustment for cFGF23 abolished the association (HR, 0.96; 95% CI, 0.84 to 1.08; P=0.45; Figure 1D).
Discussion
In this study, we show that, in RTRs, ID is accompanied by markedly higher cFGF23 levels than in RTRs without ID, whereas iFGF23 levels were similar in both groups. Importantly, variation in cFGF23 explained a considerable part of the association between ID and mortality. Similar findings were obtained in sensitivity analyses using an alternative definition of ID, and using TSAT and ferritin as individual components of ID. Thus, this study confirms earlier findings that iron plays an essential role in FGF23 production and metabolism, extends these findings to a novel patient setting (i.e., RTRs), and supports the notion that the biologic process by which ID simultaneously upregulates FGF23 production and its cleavage, or alternatively the increased level of cFGF23 fragments, is an important mediator of the association between ID and mortality.
Post-transplantation anemia is highly prevalent in RTRs, affecting approximately one third of the population.2,4 Furthermore, it has been widely documented that post-transplant anemia is associated with poor outcomes.1,15,16 Recently, we showed that ID, independent of anemia, is associated with an increased risk of mortality, shifting the focus from anemia to ID.4 The etiology of ID-related risk of mortality is unknown. In this study, we have shown that the association of ID with mortality is largely explained by cFGF23. The extent of this effect is illustrated by the restricted cubic splines depicting the individual components of ID (Figure 1). The association between low serum ferritin and low TSAT and mortality is abrogated by adjustment for cFGF23, whereas the association between high serum ferritin levels and mortality was not altered. In mediation analyses, cFGF23 was a prominent statistical mediator of the association between ID and mortality. Therefore, it seems that an ID-induced rise in cFGF23 levels plays an important role in the outcome of RTRs. Indeed, cFGF23 has been shown to be an independent risk factor for death in various patient groups, including postoperative AKI, nondialysis CKD, ESRD, and RTRs.5,6,17–20 It seems likely that a so far unknown, underlying process driving FGF23 production and cleavage plays a role in adverse outcomes of ID, but not iron overload or inflammation. The fact that adjustment for cFGF23 did not materially alter the upper part of the nonlinear association that was present in sensitivity analyses on the association of serum ferritin with mortality is supportive of our hypothesis that inflammation is not the driving force behind the association of ID and mortality. In line with this observation, the association of ID with mortality persisted upon adjustment for hs-CRP. Similarly, renal function in itself also does not seem to be the sole underlying factor, because adjustment for eGFR did not materially change the results either in this study or in previous studies addressing the association between cFGF23 and mortality.5,6 As an alternative to an underlying factor or process mediating both FGF23 cleavage and adverse outcome, it may be that C-terminal fragments in themselves drive adverse outcome. Previously, C-terminal fragments have been found to compete with iFGF23 for binding to its receptor complex and function as a competitive inhibitor, which may impair phosphaturia and aggravate soft tissue calcification.21
The present data confirm the concept that iron plays an important role in regulating FGF23 production and metabolism. FGF23 is regulated by a complex, partly unrevealed, interplay between local bone factors that modulate bone turnover and mineralization, and systemic factors that regulate mineral metabolism.22 PTH, 1,25-dihydroxyvitamin D, klotho, glucocorticoids, calcium, and phosphate are all known to regulate FGF23 production.23,24 It has been established that systemic ID stabilizes hypoxia-inducible factor 1-α (HIF-1α), which increases FGF23 transcription9 and simultaneously upregulates furin, which cleaves FGF23.25,26 Normal osteocytes couple increased FGF23 production with commensurately increased FGF23 cleavage, which ensures normal phosphate homeostasis in the event of ID26 because the intact, biologically active iFGF23 levels remain relatively unchanged. In keeping with this hypothesis, in our iron-deficient RTRs’ cFGF23 levels were upregulated, whereas the levels of iFGF23 were similar in patients with ID and without ID.
To our knowledge, this is the first study to investigate the interplay between ID and FGF23 in RTRs. In animal models of CKD, it has been shown that ID stimulates FGF23 production but also upregulates cleavage, which leads to increased levels of circulating cFGF23, but normal levels of iFGF23.9,27,28 In humans, Wolf and colleagues showed that ID anemia is associated with normal iFGF23 but markedly elevated cFGF23 levels to an extent that is only seen in advanced CKD or in hereditary diseases of FGF23 overproduction.14 Moreover, rapid correction of ID with different intravenous iron preparations reduced cFGF23 levels by approximately 80% within 24 hours,14 consistent with ID as an important stimulus for elevated cFGF23 levels.
Our study has several strengths as well as limitations. The main strength is the large prospectively followed cohort of stable RTRs, in which several markers of ID as well as both iFGF23 and cFGF23 levels were measured. Moreover, end point evaluation was complete in all participants despite a considerable follow-up period. We also acknowledge several limitations of this study. First, because of the observational status of our single center study, we cannot exclude the possibility of residual confounding. Second, a limitation is that no gold standard for the definition of ID exists.29 In this study, the definition of ID was on the basis of a combination of two commonly used and clinically relevant markers, namely ferritin (iron load) and TSAT (iron transport availability). To increase robustness of our findings, we performed sensitivity analyses where we used an alternative definition of ID and also assessed the association of the individual iron status parameters with mortality, and demonstrated similar results.
In conclusion, we identified that iron-deficient RTRs have elevated levels of cFGF23, but not iFGF23, compared with non–iron-deficient RTRs. Importantly, ID was independently associated with mortality, and this association was to a large extent explained by variation in cFGF23 levels. Future studies are needed to unravel the complex interplay between ID, FGF23, and adverse outcomes in RTRs and other populations.
Concise Methods
Patient Population
We approached all RTRs (aged≥18 years) who were at least 1 year post-transplantation for participation in this study. RTRs were approached during outpatient clinic visits between 2008 and 2011, as described previously.30 All kidney transplantations took place in the University Medical Center Groningen (Groningen, The Netherlands). Among 817 RTRs who were approached, 707 (87%) chose to participate. We excluded patients with missing data on ID (n=7), resulting in 700 RTRs eligible for analyses. All patients provided written informed consent. All study protocols were approved by the Medical Ethical Committee (METc 2008/186) and adhered to the principles of the Declaration of Helsinki. The primary end point of this study was all-cause mortality. The continuous surveillance system of the outpatient program ensures up-to-date information on patient status. General practitioners or referring nephrologists were contacted in the case that the status of a patient was unknown. End points were recorded until the end of May of 2013 with no loss to follow-up.
Data Collection
The measurement of clinical parameters has been described in detail previously.31 In brief, information on medical history and medication use was obtained from patient records. Participants’ height and weight were measured with participants wearing indoor clothing without shoes. BP was measured according to a strict protocol as previously described.31 Information on alcohol consumption and smoking behavior was gathered using a questionnaire. Smoking behavior was classified as never, former, or current smoker.
Laboratory Procedures
Blood was drawn in the morning after an 8–12-hour overnight fast. Plasma cFGF23 levels were measured with human FGF23 (C-terminal) ELISA (Immutopics, Inc., San Clemente, CA) in stored plasma samples with intra-assay CVs between 2.2% and 4.4%, and interassay CVs between 9% and 16%.32 iFGF23 levels were measured in stored plasma samples by ELISA (Kainos Laboratories, Inc., Tokyo, Japan) with intra-assay CVs between 5.3% and 9.7%, and interassay CVs between 5.7% and 14%.32 The cFGF23 immunometric assay uses two antibodies directed against different epitopes within the C-terminal part of FGF23, which therefore detects both the intact hormone and C-terminal cleavage products. In contrast, the iFGF23 assay detects only the intact molecule.26 Transferrin was measured using an immunoturbidimetric assay (Cobas c analyzer, Modular P system; Roche Diagnostics, Mannheim, Germany). Serum ferritin concentrations were determined using the electrochemiluminescence immunoassay (Modular analytics E170; Roche Diagnostics). Serum iron was measured using photometry (Modular P800 system; Roche Diagnostics). Serum creatinine was measured using an enzymatic, isotope dilution mass spectrometry–traceable assay on a Roche P-Modular automated analyzer (Roche Diagnostics). TSAT (%) was calculated as 100× serum iron (µmol/L)/25× transferrin (g/L).33 ID was defined as TSAT<20% and ferritin<300 µg/L.34 Renal function was determined by estimating GFR by applying the CKD Epidemiology Collaboration equation.35 Proteinuria was defined as urinary protein excretion ≥0.5 g/24 h.
Statistical Analyses
Data were analyzed using IBM SPSS software, version 22.0 (SPSS Inc., Chicago, IL), R version 3.2.3 (Vienna, Austria), and STATA 14.1 (STATA Corp., College Station, TX). Data are expressed as mean±SD for normally distributed variables and as median (25th–75th IQR) for variables with a skewed distribution. Categoric data are expressed as number (percentage). We evaluated between-group differences comparing RTRs with versus without ID using t test, Mann–Whitney U test, or chi-squared test, as appropriate. To study the association between ID and all-cause mortality, we used Cox proportional hazards regression analysis. We performed analyses in which we first adjusted for age and sex (model 1); and additionally for eGFR, proteinuria, time since transplantation, primary renal disease, history of cardiovascular disease, and smoking status (model 2). In further models, we adjusted for iFGF23 (model 3) and cFGF23 (model 4). Because of skewed distribution, iFGF23, cFGF23, and time since transplantation were natural log–transformed. We tested for colinearity between ID and cFGF23 by calculating a correlation coefficient and a VIF score. A correlation coefficient of <0.7 and a VIF<5 indicate no evidence for colinearity.36 Potential effect modifications by age, sex, eGFR, proteinuria, history of cardiovascular disease, time since transplantation, smoking status, BMI, presence of diabetes, hs-CRP, serum calcium, phosphate, PTH, use of ACE-inhibitors, use of diuretics, cFGF23, and iFGF23 concentrations were tested by fitting models containing both main effects and their crossproduct terms. In sensitivity analyses, we repeated the Cox regression analysis by using an alternative, frequently used definition of ID: TSAT<20% and ferritin<200 µg/L.37 We also evaluated ID using the iron status components (i.e., TSAT and ferritin) individually. Splines of individual iron status components with all-cause mortality were fit using a Cox proportional hazards regression model on the basis of restricted cubic splines in univariate analyses and after adjustment for cFGF23. As sensitivity analyses, we performed adjustments of the association of ID with mortality as in Table 2, model 2, for serum calcium, serum PTH, use of ACE inhibitors and AII antagonists, hemoglobin levels, acute rejection events, presence of diabetes, use of iron supplements, and hs-CRP levels, each time in addition to existing adjustment for age, sex, eGFR, proteinuria, time since transplantation, primary renal disease, and smoking status. Finally, we performed mediation analyses with the method described by Preacher and Hayes, which is on the basis of logistic regression. These analyses allow for testing significance and magnitude of mediation.38,39 In all analyses, a two-sided P value <0.05 was considered significant.
Simulation Study
In order to investigate the effect of differences in intra-assay CVs between cFG23 and iFGF23 we simulated values of cFGF23 100 times. In each simulation, normally distributed noise with a mean of 0 and an SD equal to 4.3% of the cFGF23 value was added to the original cFGF23 value in order to simulate an intra-assay CV between 9.5% and 10%, which is similar to that of iFGF23. The Cox regression and the mediation analyses were repeated for each of these simulated cFGF23 variables and the means of the HRs, their 95% CIs, and the proportion of mediation were calculated.
Disclosures
C.A.J.M.G. received speaking fees and research funding from Vifor Pharma. The other authors have declared that no conflict of interest exists.
Acknowledgments
We thank Wendy Dam for measuring cFGF23 levels in the plasma samples of the cohort.
Generation of the TransplantLines Food and Nutrition Biobank and Cohort Study (TxL-FN) was funded by Top Institute Food and Nutrition. M.H.d.B. is supported by a Veni grant (The Netherlands Organization for Scientific Research, grant no 016.146.014).
Footnotes
M.H.d.B. and C.A.J.M.G. are the cosenior authors.
Published online ahead of print. Publication date available at www.jasn.org.
- Copyright © 2017 by the American Society of Nephrology