The use of absolute values improves performance of estimation formulae: a retrospective cross sectional study

Abstract

Background

Estimation of Glomerular Filtration Rate (GFR) by equations such as Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) or Modification of Diet in Renal Disease (MDRD) is usually expressed as a Body Surface Area (BSA) indexed value (ml/min per 1.73 m²). This can have severe clinical consequences in patients with extreme body sizes, resulting in an underestimation in the case of obesity or an overestimation of GFR in the case of underweight patients. The aim of this study was to compare the performance of both estimation formula expressed in ml/min, instead of ml/min per 1.73 m², with a reference method.

Methods

Retrospective single centre cross sectional study of 185 patients. GFR was measured with ⁵¹Cr-EDTA and estimated with CKD-EPI and MDRD. Bias, precision and accuracy of absolute estimated GFR was calculated.

Results

Bias of CKD-EPI and MDRD formulae expressed as an absolute value was 0.49 and 0.27 ml/min respectively, which is lower than previously reported. Precision was 12.95 and 16.33 and accuracy expressed as P30 was over 92.43% for CKD-EPI. There were no significant differences in GFR between the reference method and the estimation formulae.

Conclusions

The performance of CKD-EPI and MDRD formulae can be significantly improved in the individual patient if the absolute values are used by removing the BSA normalization factor. Absolute estimated GFR by CKD-EPI is comparable to measured GFR, improving the performance of this formula in the assessment of individual kidney function, thus providing clinicians with an alternative to reference methods.

Background

Glomerular filtration rate (GFR) is accepted as the best indicator of kidney function in clinical practice, and its measurement using urinary or plasma clearance of exogenous filtration markers is considered the gold standard [1]. In clinical practice, measuring GFR by these methods is impractical, expensive and in some cases unnecessary. This is one of the reasons why estimating equations have emerged.

Estimated GFR (e-GFR) is nearly always expressed as a Body Surface Area (BSA) indexed value (ml/min per 1.73 m²).

Measured GFR (m-GFR) can be expressed as an absolute value (ml/min) or a BSA indexed value (ml/min per 1.73 m²). These results are comparable as long as the patient has a BSA around 1.73 m², which is the reference value in use over the last century. Whether this value of 1.73 m² reflects a normal BSA nowadays, is disputable [2].

It is important for clinicians to be aware of the implications of using BSA indexed GFR values in populations with extreme BSA values, as suggested by an increasing amount of scientific papers investigating relationship between GFR and BSA over the last ten years [3–5].

Several studies have demonstrated that normalization of GFR with BSA in those cases with extreme body sizes might be misleading [6] and other studies have documented that obese patients have higher absolute GFR values than lean patients, independent of diabetes and hypertension [7–10].

There are different reasons why we keep using BSA as a normalization index of GFR, as already elegantly discussed by Delanaye [11], but the most important of them is to allow comparisons of different populations.

The impact of indexing GFR by BSA can be rather important if the results are to be used for drug dosing [12] or evaluation of kidney living donor in those patients outside the normal range of BSA [13]. In both cases, knowing the real individual kidney function is rather important. This is already recommended by the National Kidney Education Program (NKDEP) that suggests using the absolute GFR values in individual drug dosing of patients with extreme body sizes, as well as to consider measuring renal function with exogenous filtration markers when prescribing drugs with narrow therapeutic indices [14].

In practice very few are aware of these recommendations and most professionals use the BSA indexed GFR values.

For these reasons we decided to investigate whether the use of estimated absolute GFR values obtained with Modification of diet in renal disease (MDRD) and chronic kidney disease epidemiology collaboration (CKD-EPI) formula is comparable to absolute GFR values obtained with a gold standard method, thus avoiding professionals the need of investigating patients with expensive isotope methods, without compromising the quality of the results.

Methods

Retrospective single centre cross-sectional study of patients who had their kidney function measured with a reference method, 51 chromium ethylene diamine tetraacetic acid clearance (⁵¹Cr-EDTA).

An analysis of differences between absolute GFR values obtained with ⁵¹Cr-EDTA and two creatinine-based estimation formulae: CKD-EPI and MDRD was performed.

The study was approved by the Danish Data Protection Agency (j.nr 2007-41-1006) and conducted in agreement with the principles of the Declaration of Helsinki.

Participants

A total of 362 patients were referred to the local Department of Nuclear and Clinical Physiology in order to determine GFR by ⁵¹Cr-EDTA over a 2 year period. All patients were examined at the same laboratory.

The inclusion criteria were: Age ≥ 18; a minimum of two serum creatinine values with less than 15% variation in range within a period of three months and a serum creatinine value measured a maximum of three weeks before the ⁵¹Cr-EDTA measurement. We only used the first value of ⁵¹Cr-EDTA in those cases where this test was performed more than once.

Exclusion criteria: Only one creatinine value determination; acute renal insufficiency; diabetic ketoacidosis; increased values of protein (plasma (P) protein > 90 g/l), glucose (P glucose > 17 mmol/l), bilirubin (P bilirubin > 65 micromole/l) or uric acid (P uric acid > 0.45 mmol/l); renal replacement therapy; pregnancy; amputation; treatment with cephalosporin, cimetidine, methyldopa or trimethoprim.

The study population had not been diagnosed with end-stage renal disease (ESRD), as renal replacement therapy was an exclusion criterion. None of the patients were undergoing treatment with erythropoietin (EPO) or active D vitamin at the time of inclusion.

A total of 189 consecutive patients fulfilled the inclusion criteria.

Four were excluded due to a high serum level of uric acid [15], leaving 185 patients included for further analysis.

This patient population is now included in a larger meta-analysis [16].

Estimated glomerular filtration rate

Analysis of serum creatinine

Serum creatinine was measured with the Jaffé method (Abbott Architect C systems TMc8000, reagent 7D64) which was calibrated with Isotope Dilution Mass Spectrometry (IDMS).

Calculation of estimated glomerular filtration rate

The MDRD and CKD-EPI formulas have previously been validated in patients with CKD stage 3 or lower, defined as GFR ≤ 59 ml/min per 1.73 m². Validation was set at this level because higher values are associated with increased imprecision.

CKD-EPI equation: e ‒ GFR = 141 × min (Scr/k, 1)ª × max (Scr/k, 1)^‒ 1.209 × 0.993^Age × 1.018 [if female] × 1.159 [if black].

Where Scr is serum creatinine, k is 0.7 for females and 0.9 for males, a is -0.329 for females and -0.411 for males, min indicates the minimum of Scr/k or 1, and max indicates the maximum of Scr/k or 1.

Modification of diet in renal disease formula (14) based on serum creatinine (S_c), age, gender and race

$$\begin{array}{ll}\mathrm{e}‒\mathrm{GFR}\left(\mathrm{ml},/,\min ,,\mathrm{per},,1.73,{\mathrm{m}}^2\right.& =175*{\left({\mathrm{S}}_{\mathrm{cr}}/88.4\right)}^{‒1.154}*{\mathrm{age}}^{‒0.203}\\ \left.*,\left(0.742\mathrm{if}\mathrm{female}\right),*,\left(1.21\mathrm{if}\mathrm{black}\right)\right).\end{array}$$

Where S_cr is in micromoles/l; 88.4 is the molecular weight of S_cr and age is stated in years.

Body surface area (BSA) was calculated from the DuBois & DuBois formula with height in centimeters rather than meters, as BSA(m²) = 0.007184 × height (cm)^0.725 × weight (kg)^0.425.

We looked at absolute values and indexed values; therefore we used these formulae to convert one to another.

Conversion of BSA indexed GFR values to absolute GFR values

$$\begin{array}{ll}\mathrm{Absolute}\mathrm{GFR}=& \left(\mathrm{BSA}\mathrm{indexed}\mathrm{estimated}\mathrm{GFR}\times \mathrm{Patients}\mathrm{BSA}\right)\\ /1.73{\mathrm{m}}^2=\mathrm{ml}/\min .\end{array}$$

Conversion of absolute GFR values to BSA indexed GFR values

$$\begin{array}{ll}\mathrm{BSA}\mathrm{Indexed}\mathrm{GFR}& =\left(\mathrm{absolute}\mathrm{GFR}\times 1.73{\mathrm{m}}^2\right)/\mathrm{Patients}\mathrm{BSA}\\ =\mathrm{ml}/\min \mathrm{per}1.73{\mathrm{m}}^2.\end{array}$$

Measured Glomerular Filtration Rate (m-GFR)

GFR was determined from the total (renal and extra-renal) ⁵¹Cr-EDTA plasma clearance by a simplified single-injection technique with a single-plasma sample [17].

Statistics

Data was analyzed with Excel and Statistica 99 version, by StatSoft inc. All GFR measurements and estimates follow a normal distribution.

Bias, precision and accuracy of e-GFR

Bias in this study is an expression of systemic error of estimated absolute GFR and was defined as the mean difference between estimated and measured absolute GFR.

Precision is an expression of the random variation or spread of estimated GFR around the measured GFR. It was expressed as the standard deviation (SD) of the bias. Large width equals low precision. Both precision and bias were expressed as ml/min.

Accuracy is affected by both bias and imprecision and was defined as the percentage of estimated GFR values within 10% and 30% limits of the measured GFR. Values over 90% are considered as indicators of high accuracy.

Comparison of absolute m-GFR versus absolute e-GFR

We investigated whether there were differences between measured and estimated GFR whenever results were expressed independently from BSA, in other words, GFR was expressed as an absolute value (ml/min) rather than as a BSA indexed value (ml/min per 1.73 m²).

Our hypothesis was that there were no differences in between the absolute m-GFR by ⁵¹Cr-EDTA or absolute e-GFR by CKD-EPI or MDRD.A T test for unpaired data was used, values of P <0.05 were considered significant.

We also performed an analysis where the BSA indexed values where compared against the absolute values for all the three methods to ascertain whether the bias found with ⁵¹Cr-EDTA, was also confirmed with CKD-EPI and MDRD formulae. We also analyzed data with a Bland Altman plot as this is the recommended analysis to compare methods.

Results

In this study we looked at the performance of these two formulae whenever the results are expressed in absolute values, to avoid a possible influence of BSA in the results. In a previous similar study based on the same set of data we looked at similar results but only using the BSA indexed values [17].

Patient characteristics

As shown in Table 1, the mean patient age exceeded 60 years, with a range of 26 to 83 years. There were more females than males. The mean absolute m-GFR was 88 ± 25 ml/min, (Range 33–156 ml/min) which means that there were no patients with CKD IV and V included in this study.

Table 1 Patient characteristics

32 (17%) patients had m-GFR > 110 ml/min, with a median of 120 ml/min (111–156).

There were 23 (12%) patients with a GFR less than 60 ml/min (Range 33–59), median 49 ml/min.

The average BSA in this study was 1.78 ± 0.20 m², higher than the standard BSA of 1.73 m², with a range of 1.23 to 2.30 m². In comparison BMI mean was 24.2 kg/m² and range 15.6 to 36.4 kg/m².

Bias, precision and accuracy results of e-GFR with CKD-EPI and MDRD formulae

With regards to CKD-EPI, Table 2 shows that bias was lower than in the previous study based on BSA indexed data (0.49 vs. 1.16), precision was also slightly better (12.95 versus 13.37) but the best results arose from the higher P30 which crossed the level of 90%, considered as the reference value for a good accuracy (92.43% versus 89.73%) [18].

Table 2 Bias, precision and accuracy results of e-GFR with CKD-EPI and MDRD

The bias for the whole group was even lower for MDRD, but the P30 was not above 90% and precision was lower than for CKD-EPI.

When looking at the same analysis stratified by different levels of BSA, we found that in the group of patients with BSA more than 1.80 m² with CKD-EPI formula, there is an improvement in the results as compared to MDRD, with a bias of -0.68 ml/min, a high precision (11.85 vs. 12.95) and a higher P10 (56.79 vs. 52.43%) as well as P30 (98.77 vs. 92.43%).

Comparison of measured absolute values versus estimated absolute values

Table 3 shows the bias between absolute values obtained with the gold standard method and both estimation formulae. In other words, the normalization factor was removed.

Table 3 Comparison of CKD-EPI and MDRD versus 51Cr-EDTA as absolute GFR values (ml/min)

This analysis showed that there were no differences in between GFR values obtained with ⁵¹Cr-EDTA and CKD-EPI or MDRD estimation formulae, as long as the absolute values are used.

An analysis of data stratified by BSA showed that there was a clinical significant difference in the group of males with BSA < 1.60 m² for both formula, but there were only 2 participants in the group. This difference was disregarded.

By looking at the data in a graphic manner with Bland-Altman plot it is quite clear to see that the comparison between ⁵¹Cr-EDTA and CKD-EPI (Figure 1), demonstrates that most values are around 0 with a lesser spreading of data and no outliners. Data for MDRD is more spread, showing that this method is less comparable to ⁵¹Cr-EDTA (Figure 2).

Figure 1

Bland altman plot 51cr-EDTA versus CKD-EPI formula. 51cr-EDTA: 51 chromium ethylene diamine tretraacetic acid clearance; CKD-EPI: Chronic kidney disease epidemiology collaboration.

Figure 2

Bland altman plot 51cr-EDTA versus MDRD formula. 51cr-EDTA: 51 chromium ethylene diamine tretraacetic acid clearance; MDRD: Modification of diet in renal disease.

Comparison of BSA indexed and absolute GFR values for all three methods

When differences in between indexed and absolute values for all three methods were analyzed, as shown in Table 4, a very similar bias was found across all three methods, suggesting that the bias is related to the BSA normalization factor.

Table 4 Differences in between BSA indexed GFR (ml/min per 1.73 m ² ) and absolute GFR values (ml/min) for 51Cr-EDTA, CKD-EPI and MDRD

There was a significant bias in all three groups as a whole as well as whenever BSA was > 1.80 m².

These results showed that in patients with the lowest BSA under 1.60 m² there was a significant overestimation of GFR patent with all three methods, demonstrating that the tendency can be found with the gold standard method as well as the estimation formula.

There were no clinical or statistical differences between mean bias for patients with at BSA <1.60 m² measured with ⁵¹Cr-EDTA (10.23 ± 4.76 ml/min) and mean bias for the same group calculated with CKD-EPI (11.08 ± 5.82 ml/min) as p value was 0.51.

The same results applied to the comparison between ⁵¹Cr-EDTA and MDRD where p value was 0.53.

The opposite effect was observed in the group of patients with BSA over 1.80 m², as there was an underestimation of the measured GFR with all three methods if the BSA indexed GFR value was used.

Again there were no differences in the mean bias between all three methods, as underestimation of GFR by ⁵¹Cr-EDTA of -12.53 ± 7.84 ml/min compared to underestimation by CKD-EPI of -12.31 ± 7.37 ml/min (n = 81) showed a p value of 0.85. Same analysis with MDRD values showed a p value of 0.78.

Table 5 shows the analysis of bias between ⁵¹Cr-EDTA and the two estimation methods stratified by age. The chosen range is 25 to >70 years, as the youngest patient was 25.64 and the oldest was 83.42 years old. As demonstrated in Table 3 by BSA stratification, no differences are found in the 3 chosen groups. The best results can be observed with CKD-EPI formula in the group of patients with an age range between 51 and 70 years, as bias is as low as 0.07 and no significant clinical or statistical differences between the reference method and the estimation formula are found.

Table 5 Comparison of CKD-EPI and MDRD versus 51Cr-EDTA as absolute GFR values (ml/min) stratified by age

Discussion

To our knowledge this is the first clinical study of a group of patients in a creatinine steady state, comparing the performance of MDRD and CKD-EPI GFR estimation equation expressed in absolute values and not BSA indexed values, with a reference method for measuring GFR.

The absolute GFR values obtained with both estimation formulae are very similar to the results obtained with ⁵¹Cr-EDTA, suggesting that the performance of both formulae can be improved if the BSA normalization factor is removed, but in general CKD-EPI performance was slightly better measured with different parameters.

This study also shows what other studies have suggested, with increasing BSA there is an underestimation of GFR as long as the BSA indexed values are used, what is new in this study is that this underestimation is eliminated if the absolute values are used, improving the performance of CKD-EPI and MDRD formulae in all groups regardless BSA.

The best results regarding bias, precision and accuracy were obtained for CKD-EPI formula in the group of patients with BSA >1.80, showing a P30 of 98.77%, and an underestimation of only 0.68 ml/min with regards to ⁵¹Cr-EDTA.

By looking at the different age groups, the best results were also found with the CKD-EPI formula in the range of 51 to 70 years with a bias as low as 0.07 ml/min.

The Bland Altman plots show that the best results are seen with CKD-EPI, with a good concentration of values around the bias line.

A review about BSA indexing renal function parameters by Delanaye et al. [19] has argued against using BSA indexed measured GFR values, especially if the results are to be used for drug dosing in an obese population or GFR longitudinal follow-up purposes, but they do not recommend to “de-index” MDRD in patients with a BSA over 2.4 m², because the MDRD formula was based on patients with a BSA range between 1.5 and 2.4 m².

In this study we demonstrate that de-indexation of estimated GFR by MDRD or CKD-EPI follows the same trend as absolute measured GFR by ⁵¹Cr-EDTA, but our population has a BSA range of 1.23 to 2.30 m².

By expressing the raw estimated GFR values, factors such as cachexia or obesity are removed, and that might avoid under and overestimation largely observed in many studies, which could contribute to a more accurate dose calculation of GFR based drugs, such as carboplatin and a reduction of the inaccuracies in drug dosing reported by other studies [5, 20].

There are many studies looking at the problems related to BSA based GFR estimation formula, suggesting other parameters such as lean body mass [3] but if the goal of estimating GFR is drug dosing or an individual assessment of kidney function, the absolute values without further correction might be more accurate and easy to calculate.

In those cases where GFR is used to compare populations, the use of extracellular volume (ECV) as normalization factor seems to be more adequate as BSA, body mass index (BMI) and lean body mass are all related to weight and height, but measuring ECV is not feasible in daily practice [4].

Conclusions

In view of these facts we recommend clinicians to convert the estimated GFR values obtained with CKD-EPI to the absolute GFR values whenever the individual kidney function needs to be investigated, as bias with regards to the reference method is rather low. These considerations might also apply to other groups of patients when calculating dose of nephrotoxic drugs, as recommended by NKDEP or in assessing GFR for other purposes such as suitability for living kidney donation [13] or indication for renal replacement therapy.

Relevant organizations might consider expressing GFR estimation equations as absolute values as well as BSA indexed values, leaving clinicians the opportunity to use the most relevant value for a given patient in different clinical situations [21].

A prospective study looking at the impact in chemotherapy drug dosing intention by using GFR estimation formula expressed in absolute values is desirable, as such an study might show whether there are clinical consequences or not.