Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Estimated Glomerular Filtration Rate Is a Poor Predictor of the Concentration of Middle Molecular Weight Uremic Solutes in Chronic Kidney Disease

  • Nathalie Neirynck ,

    Contributed equally to this work with: Nathalie Neirynck, Sunny Eloot

    n.neirynck@ugent.be

    Affiliation Nephrology Section, Department of Internal Medicine, Ghent University Hospital, Gent, Belgium

  • Sunny Eloot ,

    Contributed equally to this work with: Nathalie Neirynck, Sunny Eloot

    Affiliation Nephrology Section, Department of Internal Medicine, Ghent University Hospital, Gent, Belgium

  • Griet Glorieux,

    Affiliation Nephrology Section, Department of Internal Medicine, Ghent University Hospital, Gent, Belgium

  • Daniela V. Barreto,

    Affiliations INSERM U-1088, Amiens, France, Clinical Research Centre, Division of Clinical Pharmacology, Amiens University Hospital, Amiens, France, and the Jules Verne University of Picardy, Amiens, France

  • Fellype C. Barreto,

    Affiliations INSERM U-1088, Amiens, France, Clinical Research Centre, Division of Clinical Pharmacology, Amiens University Hospital, Amiens, France, and the Jules Verne University of Picardy, Amiens, France

  • Sophie Liabeuf,

    Affiliations INSERM U-1088, Amiens, France, Clinical Research Centre, Division of Clinical Pharmacology, Amiens University Hospital, Amiens, France, and the Jules Verne University of Picardy, Amiens, France

  • Aurélie Lenglet,

    Affiliations INSERM U-1088, Amiens, France, Division of Nephrology, Amiens University Hospital, Amiens, France

  • Horst D. Lemke,

    Affiliation EXcorLab GmbH, Obernburg, Germany

  • Ziad A. Massy,

    Affiliations INSERM U-1088, Amiens, France, Clinical Research Centre, Division of Clinical Pharmacology, Amiens University Hospital, Amiens, France, and the Jules Verne University of Picardy, Amiens, France, Division of Nephrology, Amiens University Hospital, Amiens, France

  • Raymond Vanholder

    Affiliation Nephrology Section, Department of Internal Medicine, Ghent University Hospital, Gent, Belgium

Abstract

Background

Uremic solute concentration increases as Glomerular Filtration Rate (GFR) declines. Weak associations were demonstrated between estimated GFR (eGFR) and the concentrations of several small water-soluble and protein-bound uremic solutes (MW<500Da). Since also middle molecular weight proteins have been associated with mortality and cardiovascular damage in Chronic Kidney Disease (CKD), we investigated the association between several eGFR formulae and the concentration of Low Molecular Weight Proteins (LMWP) (MW>500Da).

Materials and Methods

In 95 CKD-patients (CKD-stage 2–5 not on dialysis), associations between different eGFR-formulae (creatinine, CystatinC-based or both) and the natural logarithm of the concentration of several LMWP’s were analyzed: i.e. parathyroid hormone (PTH), Cystatin C (CystC), interleukin-6 (IL-6), tumor necrosis factor-alpha (TNF-α), leptin, retinol binding protein (RbP), immunoglobin light chains kappa and lambda (Ig-κ and Ig-λ), beta-2-microglobulin (β2M), myoglobin and fibroblast growth factor-23 (FGF-23)).

Results

The regression coefficients (R2) between eGFR, based on the CKD-EPI-Crea-CystC-formula as reference, and the examined LMWP’s could be divided into three groups. Most of the LMWP’s associated weakly (R2 <0.2) (FGF-23, leptin, IL-6, TNF-α, Ig-κ, Ig-λ) or intermediately (R2 0.2–0.7) (RbP, myoglobin, PTH). Only β2M and CystC showed a strong association (R2 >0.7). Almost identical R2-values were found per LMWP for all eGFR-formulae, with exception of CystC and β2M which showed weaker associations with creatinine-based than with CystC-based eGFR.

Conclusion

The association between eGFR and the concentration of several LMWP’s is inconsistent, with in general low R2-values. Thus, the use of eGFR to evaluate kidney function does not reflect the concentration of several LMWP’s with proven toxic impact in CKD.

Introduction

Chronic Kidney Disease (CKD) is an independent risk factor for mortality and cardiovascular disease (CVD) [1]. As Framingham risk calculation cannot correctly predict this risk [2], [3], other than traditional risk factors are at play. When kidney function declines, retention of uremic solutes with potential to cause vessel damage and other toxic effects, conceivably plays a role in this [4], [5].

Glomerular Filtration Rate (GFR) is used to express kidney function and this can accurately be measured by time-consuming and labor-intensive methods [6]. In clinical practice, serum creatinine (Crea) based formulae are used to calculate estimated GFR (eGFR), which offer an acceptable estimate of measured GFR (mGFR) [6][9]. However, if possible, mGFR is to be preferred as it may differ from eGRF especially in the lower GFR range in a CKD population or in patients with a body constitution that deviates from the average [6], [10]. On the other hand, measuring GFR by one of these techniques is more costly and labor-intensive than to determine eGFR. Also, current guidelines classify CKD based on the Modification of Diet in Renal Disease study (MDRD) formula [11], [12]. More recently, the CKD-EPI-Crea formula [13] has been proposed as a valid alternative, especially if eGFR is >60 ml/min/1.73 m2 [9], so that it possibly will be incorporated into the upcoming KDIGO guideline [14].

Since concentrations of uremic solutes rise when GFR deteriorates, it has been thought that GFR reflects the retention state of the patient and that the elevation of individual solute concentration of uremic toxins is closely related to the gradual deterioration of GFR. However, Eloot et al. [15] found very low regression coefficients between eGFR and several low molecular weight retention solutes in a CKD population.

The low molecular weight proteins (LMWP) are among the main representatives of the middle molecules, the third family of uremic retention solutes [16], and are interesting to study for their relationship with eGFR as with normal kidney function they are freely filtered through the glomerular basement membrane (GBM) and then mainly degraded into amino acids by the proximal tubules [17]. Furthermore the concentrations of several of the investigated LMWP’s, such as inflammatory parameters and FGF-23, are already elevated in patients with a moderate reduction in GFR [18][21] or in more advanced CKD [22]. As a consequence, associations between these solutes and eGFR are often assumed. Assessing the predictive value of eGFR for their concentration is furthermore also relevant, because several LMWP’s, such as interleukin-6 (IL-6) [23][26], tumor necrosis factor-alpha (TNF-α) [27], [28], beta-2-microglobulin (β2M) [29], [30], and fibroblast growth factor-23 (FGF-23) [31][33], have been linked to mortality or surrogate outcomes like vascular damage or progression of kidney failure. In addition, active removal of middle molecules by dialysis has been associated with better outcome [34].

Therefore, we investigated in a CKD population whether the concentration of several LMWP’s would associate with eGFR, calculated by several eGFR formulae.

Materials and Methods

Ethics Statement

The study was approved by the local ethical committee (Comité Consultatif de Protection des Personnes dans la Recherche Biomédicale (CCPPRB) de Picardie, CHU Amiens, Amiens, France) and performed in accordance to the Declaration of Helsinki. Written informed consent was obtained from all patients.

Study Population

This evaluation is a planned sub-analysis of a study undertaken over an 18-month period (January ’06- June ’07), which screened 150 Caucasian patients with prevalent CKD stage 2–5D from the Nephrology Department at Amiens University Hospital, in which uremic retention solutes in relation to clinical outcomes were analyzed [23], [35][37].

All patients were over 40 years old and had a confirmed diagnosis of CKD (two previous eGFR measures of <90 ml/min, calculated by the Cockcroft-Gault formula with an interval of 6–9 months) [38]. Exclusion criteria were chronic inflammatory disease, atrial fibrillation, complete heart block, abdominal aorta aneurysm, an aortic and/or femoral artery prosthesis, primary hyperparathyroidism, kidney transplantation, and any acute cardiovascular event in a 3 month period prior to screening for inclusion.

From the 140 patients who met the inclusion criteria, 45 were excluded from the current study because of hemodialysis treatment, which has an impact on solute concentration and on eGFR. The 95 patients, included in this study, were classified in CKD stages according to the CKD-EPI-Crea-CystC formula for further analysis [8].

Sampling and Laboratory Methods

Blood samples of all patients were collected in the morning from 9 a.m. on, centrifuged, aliquoted, frozen and stored at −80°C. Cystatin C (CystC) (MW: 13.3 kDa) concentration was determined by immune-nephelometry (N latex Cystatin C®, Siemens Healthcare, Dade Behring, Marburg, Germany) and that of intact parathyroid hormone (PTH) (MW: 9.5kDa) with a chemiluminometric immunoassay (Liaison N-tact PTH CLIA®, Diasorin, Stillwater, MN, USA). The determination of retinol binding protein (RbP) (MW: 21kDa), beta-2-microglobulin (β2M) (MW: 11.8kDa), myoglobin (MW: 17kDa) and total immunoglobulin light chains kappa (Ig-κ) and lambda (Ig-λ) (MW: 23kDa) was performed by laser nephelometry (BNProSpec®, Siemens Healthcare, Dade Behring, Marburg, Germany). ELISA’s were used to determine the levels of interleukin-6 (IL-6) (MW: 23kDa), tumor necrosis factor-alpha (TNF-α) (MW: 17kDa) (R&D Systems, Wiesbaden, Germany), and leptin (MW: 16kDa) (DRG diagnostics, Marburg, Germany). Intact fibroblast growth factor-23 (FGF-23) (MW: 32kDa), was measured by a two-site (N-terminal and C-terminal) ELISA (Immunotopics, San Clemente, CA, USA). Serum creatinine (Crea) (MW: 113Da) was measured colorimetrically by standard laboratory methods.

eGFR- Calculation

Six different formulae were used to estimate GFR: the CKD-EPI formula, based on Crea and CystC (CKD-EPI-Crea-CystC) eGFR = 177.6·Crea−0.65·CystC− 0.57·age−0.20·0.82 (if female) [8]; two formulae based on Crea alone: the MDRD eGFR = 175·Crea−1.154·age−0.203·(0.742 if female)·(1.21 if black) [7] and the CKD-EPI creatinine (CKD-EPI-Crea) eGFR = 141.min(Crea/κ,1)α·max(Crea/κ,1) −1.209·0.993Age·1.018 (if female)·1.159 (if black) (κ: 0.7 if female, 0.9 if male; α: −0.329 if female, −0.411 if male) [13]; and three formulae based on CystC alone: Stevens eGFR = 127.7·CystC−1.17·age−0.13·0.91 (if female)·1.06 (if black) [8], Le Bricon eGFR = [78·(1/CystC)]+4 [39] and Rule eGFR = 66.8·(CystC)−1.3 [40].

Statistical Analysis

The data are expressed as mean ± standard deviation and analysed by ANOVA if they were normally distributed. For data that were not normally distributed, median with interquartile range and Kruskall-Wallis test were used. Linear regressions and Pearson correlations were calculated on semi-logarithmic (LN) concentrations as a function of eGFR. Multifactorial analysis was performed to correct for well-known influencing factors for the concentration of several solutes. The regression model of CystC, β2M, IL-6, TNF-α, Ig-κ and Ig-λ was adjusted for C-reactive protein (CRP), the one of FGF-23 and PTH for calcium, phosphorus and vitamin D-supplementation, the one of leptin for body mass index (BMI) and gender, and the one of RbP for BMI, 1/CRP and diabetes mellitus. A P<0.05 was considered as statistically significant. All statistical analyses were performed using SPSS Statistics 19 (SPSS Inc, Chicago, IL) for Windows (Microsoft Corp, Redmond, WA).

Results

Ninety-five patients at different stages of CKD were included: 11.5% CKD stage 2, 39.0% CKD stage 3, 39.0% CKD stage 4, and 10.5% CKD stage 5 not on dialysis. Table 1 summarizes the demographic and clinical characteristics of the study population.

thumbnail
Table 1. Main demographic and clinical characteristics of the study population (n = 95).

https://doi.org/10.1371/journal.pone.0044201.t001

The concentrations of the studied LMWP’s, except for the immunoglobulin light chains, increased progressively with declining kidney function (Table 2).

thumbnail
Table 2. Concentrations of uremic solutes ± standard deviation according to CKD-stage (CKD-EPI-Crea-CystC).

https://doi.org/10.1371/journal.pone.0044201.t002

Our analysis primarily focused on the linear regression analysis with the natural logarithm (LN) of the concentration of each studied uremic retention solute concentration as dependent variable and the CKD-EPI-Crea-CystC eGFR as independent variable. This formula was chosen as reference because it is considered as one of the most accurate ones at this time while it incorporates both Crea and CystC, in contrast to all other studied formulae which are based on either Crea or CystC [8]. Associations between eGFR and LMWP’s were expressed as regression coefficients (R2) and are summarized in table 3 and figure 1.

thumbnail
Figure 1. Regression coefficients between LN of studies LMWP’s and eGFR.

The coefficients of the linear regression analysis between the natural logarithm of the studied low molecular weight protein concentrations and estimated Glomerular Filtration Rate, according to CKD-EPI-Crea-CystC, can be divided into 3 groups: strong (R2 >0.7), moderate (R2 0.2–0.7) and weak (R2 <0.2). The dashed lines indicate R2 = 0.2 and 0.7. All correlations were significant except for Ig-κ and Ig-λ. LN: natural logarithm, LMWP: low molecular weight protein, eGFR estimated glomerular filtration rate, R2: regression coefficient, Cyst C: Cystatin C, β2M: beta-2-microglobulin, RbP: retinol binding protein, PTH: parathyroid hormone, IL-6: interleukin-6, FGF-23: fibroblast growth factor-23, TNF-α: tumor necrosis factor-alpha, Ig-κ: immunoglobulin light chain kappa, Ig-λ: immunoglobulin light chain lambda.

https://doi.org/10.1371/journal.pone.0044201.g001

thumbnail
Table 3. Regression coefficients of LMWP’s and different eGFR formulae.

https://doi.org/10.1371/journal.pone.0044201.t003

The R2-values per individual solute were divergent; according to these, associations could be arbitrarily divided into three groups: strong (R2 >0.7), moderate (R2 0.2–0.7) and weak (R2 <0.2) (Figure 1). As expected, CystC (R2 = 0.828) was strongly associated as it is one of the used parameters in the formula. Only β2M showed a similar association (R2 = 0.770). Retinol binding protein (RbP), myoglobin and parathyroid hormone (PTH) were moderately associated to eGFR with R2-values of 0.423, 0.303 and 0.231, respectively. The association with eGFR was only weak for IL-6 (R2 = 0.117), leptin (R2 = 0.084), FGF-23 (R2 = 0.058) and TNF-α (R2 = 0.056). There was even no association for immunoglobulin light chain kappa (Ig-κ) (R2 = 0.021) and immunoglobulin light chain lambda (Ig-λ) (R2 = 0) (P = N.S.).

Figure 2 shows the dot plots of solute concentrations as a function of eGFR. Whereas the relation is strong with little scatter around the linear regression line for β2M (Panel A), the scatter is much larger for solutes with moderate to weak R2-values, as illustrated for myoglobin (Panel B), IL-6 (Panel C) and especially Ig-λ (Panel D) for which there is no association at all. The large standard deviations or wide interquartile ranges of the individual solute concentrations per CKD-stage also illustrate the large inter-individual variability of LMWP concentration within the same eGFR-range (Table 2).

thumbnail
Figure 2. Dot plots with best fit linear regression lines for LN of LMWP’s in function of eGFR.

Dot plots with best fit linear regression lines for natural logarithms of β2M, myoglobin, IL-6 and Ig-λ, as examples of strongly, moderately and weakly correlating low molecular weight proteins, in function of estimated Glomerular Filtration Rate, calculated by CKD-EPI-Crea-CystC. The dots represent the individual concentrations and the lines the best fit linear regression line with the 95% confidence interval. LN: natural logarithm, LMWP: low molecular weight protein, β2M : beta-2- microglobulin, IL-6: Interleukin-6, Ig-λ: immunoglobulin light chain lambda, EPI-Crea-CystC: CKD-EPI formula based on serum creatinine and Cystatin C, R2: regression coefficient, LN: natural logarithm.

https://doi.org/10.1371/journal.pone.0044201.g002

In addition we analyzed the correlation coefficients between the concentrations of the different LMWP and eGFR (CKD-EPI-Crea-CystC) in the group CKD stage 2–3 versus CKD stage 4–5. The correlation between RbP, PTH, myoglobin, FGF-23 and eGFR was significant in CKD-stage 2–3 while not in CKD stage 4–5. For the other investigated solutes, the R2-values in CKD stage 2–3 and CKD stage 4–5 were more confirm to each other. The respective R2-values are summarized in table 4.

thumbnail
Table 4. Regression coefficients (R2) of the concentration of LMWP’s and eGFR (CKD-EPI-Crea-CystC) comparing CKD stage 2–3 versus CKD stage 4–5.

https://doi.org/10.1371/journal.pone.0044201.t004

In a second step, the same analysis was performed with the other formulae under evaluation and compared to the results with the CKD-EPI-Crea-CystC-formula. The R2-values between all eGFR formula and individual solutes were strikingly similar with only one exception (Table 3). β2M was only moderately associated to Crea-based eGFR, with R2-values of approximately 0.55, as compared to CystC-based eGFR (R2 >0.8). In this way β2M followed the same pattern as CystC, for which this discrepancy could be attributed to whether CystC was included as a factor in the formula or not. Considering the other studied LMWP’s, only PTH showed a moderately similar trend, with R2 approximately 0.27 compared to R2 approximately 0.13, with CystC- or Crea- based eGFR-formulae, respectively.

Finally, we performed multifactorial regression analysis for the different LMWP’s with adjustment for several relevant variables. However, only two models induced a marked increase in R2-value: for leptin the association rose from weak to moderate when BMI was added to the regression model (R2 from 0.084 to 0.346), with BMI as an independent predictor for the leptin concentration. Likewise, after adjustment for CRP, IL-6 became moderately associated with eGFR (R2 from 0.117 to 0.305 after adjustment). For all other solutes there was no change in R2. (Data not shown).

Discussion

We analyzed the linear regression coefficients between the concentrations of several LMWP’s retained in CKD and different eGFR-formulae in a CKD population, stage 2–5 not on dialysis. As a main finding, the R2-values diverged considerably, ranging from high, R2 >0.7, to low, R2 <0.2. The majority of the evaluated LMWP’s associated weakly (R2 <0.2 for IL-6, TNF-α, FGF-23 and leptin) or moderately (R2: 0.2–0.7 for RbP, myoglobin and PTH). There was no correlation at all for the immunoglobulin light chains. Only CystC and β2M showed a strong association with eGFR (R2 >0.7) (Figure 1, Table 3). Although in some studies a correlation was sought for individual LMWP’s and eGFR of mGFR, this present study seeked out the association of the concentration of several LMWP’s and eGFR formulae together allowing their comparison.

The R2-values for the weakly and moderately associating LMWP’s did not differ substantially whether eGFR was calculated with the CKD-EPI-Crea-CystC-formula [8], the Crea-based formulae (MDRD [7] and CKD-EPI-Crea [13]), or the three different CystC-based formulae, (Stevens [8], Rule [40] and Le Bricon [39]) (Table 3). These low regression coefficients can partially be attributed to the known limitations of eGFR, as an index of mGFR [6], [41]. However in at least four other studies, almost identical low regression coefficients were found between mGFR, assessed with different techniques, and the concentration of RbP (R2 0.16) [42], myoglobin (R2 0.38) [43], leptin (R2 0.0004) [44] and FGF-23 (R2 0.09) [45] as in our study, be it that transformation of the concentrations varied from study to study. In addition, the imperfect reflection of true GFR by eGFR can explain that regression coefficients are substantially lower than 1, but not that the range in between individual molecules is so discordant, whereas per molecule they are almost identical (Table 3, Figure 1). There was also an unpredictable and large variability in concentrations of different solutes within each eGFR stratum (Table 2). These data suggest another reason for the sometimes deceiving associations than a discrepancy among mGFR and eGFR, namely that uremic solute concentration depends on other factors than GFR as well. In this way, our study corroborates findings in an earlier study with small water-soluble and protein-bound compounds [15], [46].

These results are somewhat unexpected from a physiological point of view as the renal clearance of these LMWP’s depends to a large extent on GFR alone. All these LMWP’s are freely filtered through the GBM, followed, at normal physiological concentrations, by an almost entire uptake by the proximal tubules via a receptor-mediated process to be degraded subsequently into amino acids in the tubular lysosomes [17], [47]. In this way, the proximal tubulus plays an important role in LMWP metabolism but without a direct contribution to their renal clearance; regarding the latter, GFR is the rate limiting step. This probably explains why we did not find any association between eGFR and the total (free plus bound) immunoglobulin light chains, in contrast to Hutchison et al who evaluated only free light chains for their association to eGFR (free Ig-κ: R2: 0.52; free Ig-λ: R2: 0.44) [48], as only the free fraction passes the GBM. This is also in contrast to the small water-soluble and protein-bound uremic toxins, for which tubular secretion and/or reabsorption play an important role in renal clearance [15], [46].

However, the concentration of small water-soluble and protein-bound solutes may be further influenced by many other factors as well, such as enzymatic metabolism, intestinal secretion/absorption, generation by intestinal flora, diet and changes in distribution volume [15], [46]. It is conceivable that also the concentration of the weakly and moderately correlating LMWP’s depends on other mechanisms than GFR, which even seem to have more important weight than GFR. Some known influencing factors such as changes in generation, homeostatic mechanisms and extra-renal clearance are summarized in table 5. Multifactorial regression analysis for the respective LMWP’s including some of these parameters, increased the R2-value as expected. E.g. for leptin, R2 rose from 0.084 to 0.346 when corrected for BMI,which was an independent covariate for leptin concentration in a model with eGFR and BMI. The R2-value between IL-6 and eGFR became 0.305 instead of 0.117, when adjusted for CRP, which, in contrast to BMI for leptin concentration, did however not independently predict IL-6 concentration in a model with eGFR and CRP. The majority of potentially influencing factors did however not importantly affect the R2-values. This suggests that other than well known mediators may influence these LMWP concentrations as well. In more advanced CKD, the influence of confounders, for example bone metabolism, is probably more important, which could explain partially that no significant associations were observed between PTH or FGF-23 and eGFR in CKD stage 4–5, while they were present in CKD stage 2–3 (table 4). Another contributing factor to this discrepancy in associations between CKD stage 2–3 versus CKD stage 4–5 might be purely mathematical, as the GFR-range in CKD 2–3 (30–90 ml/min/1.73 m2) is much larger than CKD 4–5 (±10–30 ml/min/1.73 m2).

thumbnail
Table 5. Main factors influencing the concentrations of the studied LMWP’s, other than GFR.

https://doi.org/10.1371/journal.pone.0044201.t005

This study demonstrates that eGFR is not a reliable predictor of the concentration of most of the evaluated LMWP’s, although several of them such as IL-6 [23][26], TNF-α [27], [28] and FGF-23 [31][33], [49] have been associated with mortality or with intermediate endpoints, such as vascular dysfunction or progression to end stage renal disease (ESRD) in CKD- or hemodialysis patients. Presumably, some of these solutes, especially if they would be representative for a cluster of other solutes, might by themselves become useful predictors of morbidity or mortality in CKD independently from eGFR. Based on the data collected in the present study, we investigated the mutual correlations between the concentrations of the different LMWP’s; however, we could not identify such a marker, correlating strongly to other LMWP’s without correlating to eGFR, among the investigated solutes (data not shown). This question, however, would be worthwhile to be investigated in larger populations.

In contrast to these weakly and moderately correlating LMWP’s, there is a remarkable similarity in regression coefficients between CystC and β2M. First, these molecules are the sole LMWP’s studied that result in acceptably high associations with eGFR (Table 3). Second, they associate better with CystC-based eGFR formulae [8], [39], [40] than with Crea-based ones [7], [13], the CKD-EPI-Crea-CystC [8] which contains both factors being intermediate (Table 3). Whereas this is no surprise for CystC which is included in some formulae and not in others, the pattern for β2M seems to be identical. This suggests that the kinetics of both molecules during progression of CKD depend on similar factors or at least factors with a similar impact on solute concentration. Of note, some of the non-renal elements with impact on both concentrations [50], [51], like chronic inflammatory disease or malignancy were among the exclusion criteria of this study. CystC was a superior marker of the association of GFR with outcome in a study by Peralta et al., who showed that the predictive value for mortality or CVD of eGFR <60 ml/min/1.73 m2, based on a CystC-based eGFR, was better than eGFR based on the CKD-EPI-Crea-formula [45]. Recently, in a general population, CystC and β2M were stronger predictors of mortality, CVD and evolution to ESRD than eGFR based on the CKD-EPI-Crea [52].

The present study has some shortcomings. First, the study population was rather small, with even smaller subgroups per CKD-stage. Second, we used eGFR which gives only an approximate value for glomerular filtration in comparison to more exact methods such as EDTA-clearance. We preferred to use methods which are applied on a day to day basis. As the differences in correlations are so striking, it is very likely that these findings can be extrapolated to GFR in general. The strengths of this study lie in the fact that several LMWP’s are evaluated together in the same population for different eGFR formulae based on Crea, CystC or both.

Our present data, together with the previous ones [15], showing extremely variable associations between uremic retention solutes and a surrogate of GFR, suggests that eGFR per se is an inadequate indicator of the uremic status. This is also suggested by other studies. In a CKD population, Lilitkarntakul et al. [53] demonstrated that renal function did not independently predict arterial stiffness or endothelial dysfunction while the uremic retention solutes asymmetric dimethylarginine (ADMA), isoprostanes or endothelin-A did. In the Initiating Dialysis Early and Late (IDEAL) trial [54], approximately 75% of the patients randomized to start dialysis at low eGFR (5–7 ml/min/1.73 m2), initiated dialysis earlier, mainly because of uremic symptoms.

In this study, the regression coefficients of different LMWP’s in relation to eGFR are diverse and in general low. This shows that other factors than GFR are important for the development of the ‘uremic status’. Further research is needed to evaluate whether these uremic toxins can be used as biomarkers for the risk stratification associated to uremic toxicity within the different CKD-stages and beyond eGFR.

Supporting Information

Table S1.

Main factors influencing the concentrations of the studied LMWP’s, other than GFR.

https://doi.org/10.1371/journal.pone.0044201.s001

(PDF)

Author Contributions

Conceived and designed the experiments: ZAM RV. Performed the experiments: DVB FCB SL AL ZAM. Analyzed the data: SE NN GG. Contributed reagents/materials/analysis tools: HDL ZAM SL. Wrote the paper: NN RV.

References

  1. 1. Weiner DE, Tighiouart H, Amin MG, Stark PC, MacLeod B, et al. (2004) Chronic Kidney Disease as a Risk Factor for Cardiovascular Disease and All-Cause Mortality: A Pooled Analysis of Community-Based Studies. J Am Soc Nephrol 15: 1307–1315.
  2. 2. Weiner DE, Tighiouart H, Elsayed EF, Griffith JL, Salem DN, et al. (2007) The Framingham predictive instrument in chronic kidney disease. J Am Coll Cardiol 50: 217–224.
  3. 3. Sciarretta S, Valenti V, Tocci G, Pontremoli R, Rosei EA, et al. (2010) Association of renal damage with cardiovascular diseases is independent of individual cardiovascular risk profile in hypertension: data from the Italy - Developing Education and awareness on MicroAlbuminuria in patients with hypertensive Disease study. J Hypertens 28: 251–258.
  4. 4. Vanholder R, De Smet R (1999) Pathophysiologic effects of uremic retention solutes. J Am Soc Nephrol 10: 1815–1823.
  5. 5. Vanholder R, Baurmeister U, Brunet P, Cohen G, Glorieux G, et al. (2008) A bench to bedside view of uremic toxins. J Am Soc Nephrol 19: 863–870.
  6. 6. Stevens LA, Levey AS (2009) Measured GFR as a Confirmatory Test for Estimated GFR. J Am Soc Nephrol 20: 2305–2313.
  7. 7. Levey AS, Coresh J, Greene T, Stevens LA, Zhang YL, et al. (2006) Using standardized serum creatinine values in the modification of diet in renal disease study equation for estimating glomerular filtration rate. Ann Intern Med 145: 247–254.
  8. 8. Stevens LA, Coresh J, Schmid CH, Feldman HI (2008) Estimating GFR using Cystatin C alone or in combination with creatinine: a pooled analysis of 3418 individuals with CKD. Am J Kidney Dis 51: 395–406.
  9. 9. Stevens LA, Schmid CH, Greene T, Zhang YP, Beck GJ, et al. (2010) Comparative Performance of the CKD Epidemiology Collaboration (CKD-EPI) and the Modification of Diet in Renal Disease (MDRD) Study Equations for Estimating GFR Levels Above 60 mL/min/1.73 m(2). Am J Kidney Dis 56: 486–495.
  10. 10. Kuan Y, Hossain M, Surman J, El Nahas AM, Haylor J (2005) GFR prediction using the MDRD and Cockcroft and Gault equations in patients with end-stage renal disease. Nephrol Dial Transplant 20: 2394–2401.
  11. 11. National Kidney Foundation (2002) K/DOQI clinical practice guidelines for chronic kidney disease: evaluation, classification, and stratification. Am J Kidney Dis 39: S1–266.
  12. 12. Levey AS, Eckardt KU, Tsukamoto Y, Levin A, Coresh J, et al. (2005) Definition and classification of chronic kidney disease: a position statement from Kidney Disease: Improving Global Outcomes (KDIGO). Kidney Int 67: 2089–2100.
  13. 13. Levey AS, Stevens LA, Schmid CH, Zhang YP, Castro AF, et al. (2009) A New Equation to Estimate Glomerular Filtration Rate. Ann Intern Med 150: 604–607.
  14. 14. Levey AS, de Jong PE, Coresh J, El Nahas M, Astor BC, et al. (2011) The definition, classification, and prognosis of chronic kidney disease: a KDIGO Controversies Conference report. Kidney Int 80: 17–28.
  15. 15. Eloot S, Schepers E, Barreto DV, Barreto FC, Liabeuf S, et al. (2011) Estimated glomerular filtration rate is a poor predictor of concentration for a broad range of uremic toxins. Clin J Am Soc Nephrol 6: 1266–1273.
  16. 16. Vanholder R, De Smet R, Glorieux G, Argiles A, Baurmeister U, et al. (2003) Review on uremic toxins: Classification, concentration, and interindividual variability. Kidney Int 63: 1934–1943.
  17. 17. Maack T, Johnson V, Kau ST, Fugueiredo J, Sigulem D (1979) Renal filtration, transport, and metabolism of low-molecular- weight proteins: A review. Kidney Int 16: 251–270.
  18. 18. Ix JH, Shlipak MG, Wassel CL, Whooley MA (2010) Fibroblast growth factor-23 and early decrements in kidney function: the Heart and Soul Study. Nephrol Dial Transplant 25: 993–997.
  19. 19. Isakova T, Wahl P, Vargas GS, Gutierrez OM, Scialla J, et al. (2011) Fibroblast growth factor 23 is elevated before parathyroid hormone and phosphate in chronic kidney disease. Kidney Int 79: 1370–1378.
  20. 20. Keller C, Katz R, Cushman M, Fried LF, Shlipak M (2008) Association of kidney function with inflammatory and procoagulant markers in a diverse cohort: a cross-sectional analysis from the Multi-Ethnic Study of Atherosclerosis (MESA). BMC Nephrol 9: 9 doi:10.1186/1471-2369-9-9.
  21. 21. Upadhyay A, Larson MG, Guo CY, Vasan RS, Lipinska I, et al. (2011) Inflammation, kidney function and albuminuria in the Framingham Offspring cohort. Nephrol Dial Transplant 26: 920–926.
  22. 22. Pecoits R, Heimburger O, Barany P, Suliman M, Fehrman-Ekholm I, et al. (2003) Associations between circulating inflammatory markers and residual renal function in CRF patients. Am J Kidney Dis 41: 1212–1218.
  23. 23. Barreto DV, Barreto FC, Liabeuf S, Temmar M, Lemke HD, et al. (2010) Plasma interleukin-6 is independently associated with mortality in both hemodialysis and pre-dialysis patients with chronic kidney disease. Kidney Int 77: 550–556.
  24. 24. Tripepi G, Mallamaci F, Zoccali C (2005) Inflammation markers, adhesion molecules, and all-cause and cardiovascular mortality in patients with ESRD: searching for the best risk marker by multivariate modeling. J Am Soc Nephrol 16 Suppl 1: S83–S88.
  25. 25. Panichi V, Maggiore U, Taccola D, Migliori M, Rizza GM, et al. (2004) Interleukin-6 is a stronger predictor of total and cardiovascular mortality than C-reactive protein in haemodialysis patients. Nephrol Dial Transplant 19: 1154–1160.
  26. 26. Pecoits-Filho R, Barany P, Lindholm B, Heimburger O, Stenvinkel P (2002) Interleukin-6 is an independent predictor of mortality in patients starting dialysis treatment. Nephrol Dial Transplant 17: 1684–1688.
  27. 27. Kimmel PL, Phillips TM, Simmens SJ, Peterson RA, Weihs KL, et al. (1998) Immunologic function and survival in hemodialysis patients. Kidney Int 54: 236–244.
  28. 28. Futh R, Herder C, Forster S, Muller-Scholze S, Kruse N, et al. (2004) Evaluation of diagnostic relevance of mRNA levels in peripheral blood: predictive value for mortality in hemodialysis patients. Cytokine 27: 166–172.
  29. 29. Okuno S, Ishimura E, Kohno K, Fujino-Katoh Y, Maeno Y, et al. (2009) Serum beta(2)-microglobulin level is a significant predictor of mortality in maintenance haemodialysis patients. Nephrol Dial Transplant 24: 571–577.
  30. 30. Cheung AK, Rocco MV, Yan GF, Leypoldt JK, Levin NW, et al. (2006) Serum beta-2 microglobulin levels predict mortality in dialysis patients: Results of the HEMO study. J Am Soc Nephrol 17: 546–555.
  31. 31. Isakova T, Xie HL, Yang W, Xie DW, Anderson AH, et al. (2011) Fibroblast Growth Factor 23 and Risks of Mortality and End-Stage Renal Disease in Patients With Chronic Kidney Disease. J Am Med Assoc 305: 2432–2439.
  32. 32. Jean G, Terrat JC, Vanel T, Hurot JM, Lorriaux C, et al. (2009) High levels of serum fibroblast growth factor (FGF)-23 are associated with increased mortality in long haemodialysis patients. Nephrol Dial Transplant 24: 2792–2796.
  33. 33. Yilmaz MI, Sonmez A, Saglam M, Yaman H, Kilic S, et al. (2010) FGF-23 and vascular dysfunction in patients with stage 3 and 4 chronic kidney disease. Kidney Int 78: 679–685.
  34. 34. Locatelli F, Martin-Malo A, Hannedouche T, Loureiro A, Papadimitriou M, et al., for the Membrane Permeability Outcome (MPO) Study Group (2009) Effect of Membrane Permeability on Survival of Hemodialysis Patients. J Am Soc Nephrol 20: 645–654.
  35. 35. Barreto FC, Barreto DV, Liabeuf S, Meert N, Glorieux G, et al. (2009) Serum indoxyl sulfate is associated with vascular disease and mortality in chronic kidney disease patients. Clin J Am Soc Nephrol 4: 1551–1558.
  36. 36. Liabeuf S, Barreto DV, Barreto FC, Meert N, Glorieux G, et al. (2010) Free p-cresylsulphate is a predictor of mortality in patients at different stages of chronic kidney disease. Nephrol Dial Transplant 25: 1183–1191.
  37. 37. Schepers E, Barreto DV, Liabeuf S, Glorieux G, Eloot S, et al., On behalf of the European Uremic Toxin Work Group (EUTox) (2011) Symmetric Dimethylarginine as a Proinflammatory Agent in Chronic Kidney Disease. Clin J Am Soc Nephrol 6: 2374–2383.
  38. 38. Cockcroft DW, Gault MH (1976) Prediction of creatinine clearance from serum creatinine. Nephron 16: 31–41.
  39. 39. Le Bricon T, Thervet E, Froissart M, Benlakehal M, Bousquet B, et al. (2000) Plasma cystatin C is superior to 24-h creatinine clearance and plasma creatinine for estimation of glomerular filtration rate 3 months after kidney transplantation. Clin Chem 46: 1206–1207.
  40. 40. Rule AD, Bergstralh EJ, Slezak JM, Bergert J, Larson TS (2006) Glomerular filtration rate estimated by cystatin C among different clinical presentations. Kidney Int 69: 399–405.
  41. 41. Stevens LA, Coresh J, Greene T, Levey AS (2006) Assessing kidney function–measured and estimated glomerular filtration rate. N Engl J Med 354: 2473–2483.
  42. 42. Donadio C, Lucchesi A, Ardini M, Giordani R (2001) Cystatin C, beta 2-microglobulin, and retinol-binding protein as indicators of glomerular filtration rate: comparison with plasma creatinine. J Pharm Biomed Anal 24: 835–842.
  43. 43. Hallgren R, Karlsson FA, Roxin LE, Venge P (1978) Myoglobin Turnover - Influence of Renal and Extra-Renal Factors. J Lab Clin Med 91: 246–254.
  44. 44. Menon V, Wang X, Greene I, Beck GJ, Kusek JW, et al. (2004) Factors associated with serum leptin in patients with chronic kidney disease. Clin Nephrol 61: 163–169.
  45. 45. Bacchetta J, Dubourg L, Harambat J, Ranchin B, bou-Jaoude P, et al. (2010) The Influence of Glomerular Filtration Rate and Age on Fibroblast Growth Factor 23 Serum Levels in Pediatric Chronic Kidney Disease. J Clin Endocrinol Metabol 95: 1741–1748.
  46. 46. Vanholder R, Eloot S, Schepers E, Neirynck N, Glorieux G, et al. (2012) An Obituary for GFR as the Main Marker for Kidney Function? Semin Dial 25: 9–14.
  47. 47. Christensen EI, Birn H (2001) Megalin and cubilin: synergistic endocytic receptors in renal proximal tubule. Am J Physiol - Renal Physiol 280: F562–F573.
  48. 48. Hutchison CA, Harding S, Hewins P, Mead GP, Townsend J, et al. (2008) Quantitative Assessment of Serum and Urinary Polyclonal Free Light Chains in Patients with Chronic Kidney Disease. Clin J Am Soc Nephrol 3: 1684–1690.
  49. 49. Gutierrez OM, Mannstadt M, Isakova T, Rauh-Hain JA, Tamez H, et al. (2008) Fibroblast growth factor 23 and mortality among patients undergoing hemodialysis. N Engl J Med 359: 584–592.
  50. 50. Drueke TB, Massy ZA (2009) Beta2-microglobulin. Semin Dial 22: 378–380.
  51. 51. Seronie-Vivien S, Delanaye P, Pieroni L, Mariat C, Froissart M, et al. (2008) Cystatin C: current position and future prospects. Clin Chem Lab Med 46: 1664–1686.
  52. 52. Astor BC, Shafi T, Hoogeveen RC, Matsushita K, Ballantyne CM, et al. (2012) Novel Markers of Kidney Function as Predictors of ESRD, Cardiovascular Disease, and Mortality in the General Population. Am J Kidney Dis 653–662.
  53. 53. Lilitkarntakul P, Dhaun N, Melville V, Blackwell S, Talwar DK, et al. (2011) Blood pressure and not uraemia is the major determinant of arterial stiffness and endothelial dysfunction in patients with chronic kidney disease and minimal co-morbidity. Atherosclerosis 216: 217–225.
  54. 54. Cooper BA, Branley P, Bulfone L, Collins JF, Craig JC, et al. (2010) A randomized, controlled trial of early versus late initiation of dialysis. N Engl J Med 363: 609–619.