1,852
Views
0
CrossRef citations to date
0
Altmetric
Clinical Study

Validation of the EKFC equation for glomerular filtration rate estimation and comparison with the Asian-modified CKD-EPI equation in Chinese chronic kidney disease patients in an external study

, , , , , , & show all
Article: 2150217 | Received 29 Jun 2022, Accepted 02 Nov 2022, Published online: 12 Jan 2023

Abstract

Objectives

The aim of this study is to determine whether new European Kidney Function Consortium (EKFC) equation is more applicable than Asian-modified CKD-EPI equation in clinical practice, having a higher accuracy in estimating GFR in our external CKD population.

Methods

We calculated estimated GFREKFC and GFRCKD-EPI independently using the EKFC and Asian-modified CKD-EPI formulas, respectively. The clinical diagnostic performance of the two equations was assessed and compared by median bias, precision, accuracy (P30) and so on, using 99mTc-DTPA dual plasma sample clearance method as a reference method for GFR measurement (mGFR). The equation that met the following targets was superior: (1) median bias within ± 3 mL/min/1.73 m2; (2) P30 > 75%; and (3) better precision and 95% limits of agreement in Bland–Altman analysis.

Results

Totally, 160 CKD patients were recruited in our external cohort. GFREKFC was highly related to mGFR, with a regression equation of GFREKFC=mGFR × 0.87 + 5.27. Compared with the Asian-modified CKD-EPI equation, EKFC equation demonstrated a wider median bias (–1.64 vs. 0.84 mL/min/1.73 m2, p < 0.01) that was within 3 mL/min/1.73 m2 and not clinically meaningful. Furthermore, the precision (12.69 vs. 12.72 mL/min/1.73 m2, p = 0.42), 95% limits of agreement in Bland–Altman analysis (42.4 vs. 44.4 mL/min/1.73 m2) and incorrect reclassification index of the two target equations were almost identical. Although, EKFC equation had a slightly better P30 (80.0% vs. 74.4%, p = 0.01).

Conclusions

The overall performance of EKFC equation is acceptable. There is no clinically meaningful difference in the performance of the Asian-modified CKD-EPI and EKFC equations within the limits imposed by the small sample size.

Introduction

It is especially crucial for clinical laboratories to report an accurate glomerular filtration rate (GFR), which is a positive alternative for renal function assessment [Citation1,Citation2]. Accurate renal function assessment plays an important role in the early diagnosis, treatment adjustment, and prognostic management of chronic kidney disease (CKD) patients. GFR has been estimated by various simple creatinine-based equations in recent decades [Citation3–5] as an alternative to inulin as the gold standard. We focus on how accurate the developed equations are in our external CKD population.

Our previous study also assessed the clinical utility of various creatinine-based equations using the 99mTc-diethylenetriaminepentaacetic acid (99mTc-DTPA) dual plasma sample clearance method as a reference method, which was accepted as the standard GFR assessment method at the 21st International Annual Meeting of the Society of Nuclear Medicine Europe, and demonstrated that the Asian-modified CKD-EPI equation and full age spectrum (FAS) formula were clinically more acceptable for GFR estimation with lower median bias, more superior precision and better P30 [Citation6–9].

In 2021, Pottel et al. developed and validated a new creatinine-based equation, the new European Kidney Function Consortium (EKFC) equation [Citation5], which can be applied to the full spectrum of age and renal function by combining the properties of the FAS [Citation3] and CKD-EPI equations [Citation10]. The internal development and external validation datasets suggested improved accuracy and precision compared with commonly used equations for estimating GFR from Scr levels. However, whether the EKFC equation is more accurate than the Asian-modified EPI-CKD equation for GFR estimation remains unknown, which was verified to be more accurate and convenient to use in our external CKD patients. The current study validated the EKFC equation and compared it with the Asian-modified CKD-EPI formula to recommend a more suitable equation for clinical GFR assessment in our external CKD population.

Methods

Ethics statement

All procedures performed in study involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. The study protocol was approved by the Hebei Medical University ethical committee (No. 2017-028-1), and written informed consent was obtained from each participant.

Study subjects

This study was a prospective cohort study. Chinese patients who were diagnosed with chronic kidney disease (CKD) according to the National Kidney Foundation-Kidney Disease Outcomes Quality Initiative (K/DOQI) clinical practice guidelines [Citation11] and who underwent 99mTc-DTPA dual plasma clearance measurement were enrolled in the study cohort. We collected serum samples for the serum creatinine examination. Patients with acute kidney function deterioration, edema, cardiac insufficiency, pleural or abdominal effusion, or disabled limbs or who underwent treatment with cimetidine, trimethoprim or replacement therapy were excluded [Citation12].

Staging criteria [Citation13]: According to the mGFR determined by the 99mTc-DTPA dual plasma sample clearance method, all patients were divided into 5 groups: normal renal function group (GFR > 90 mL/min/1.73 m2); mildly decreased renal function group (GFR 60–90 mL/min/1.73 m2); moderately decreased renal function group (GFR 30–59 mL/min/1.73 m2); severely decreased renal function group (GFR 15–29 mL/min/1.73 m2); and renal failure group (GFR <15 mL/min/1.73 m2).

Grouping criteria: The patients were categorized into 2 subgroups based on the mGFR calculated by the reference method: the lower-GFR subgroup (mGFR ≤ 60 mL/min/1.73 m2) and the higher-GFR subgroup (mGFR > 60 mL/min/1.73 m2).

Laboratory measurements

mGFR measurement by the 99mTc-DTPA dual plasma sample clearance method

At 2 and 4 h after intravenous injection of 99mTc-DTPA, 3 mL heparin anti-coagulated blood samples were collected from the elbow vein, and the plasma radioactivity was measured by a multifunction well counter. Then, the clearance of 99mTc-DTPA (Cl′) was calculated from a single exponential: Cl′ = [D×ln(P1/P2)]/(t2t1) ×exp[(t1×ln(P2)–t2×ln(P1))/(t2t1)]. Cl′ was corrected to GFR by Brochner–Mortensen’s formula [Citation14]: GFR= 0.990778Cl′ − 0.001218Cl′Citation2. The GFR was also standardized for a BSA of 1.73 m2, mGFR = GFR× (1.73/BSA), according to the Haycock formula [Citation15] of BSA (m2) = 0.024265 × Wt0.5378 × Ht0.3964, using the patients’ height (cm) and weight (kg) parameters. Detailed procedures are described in a previous work [Citation12].

GFR measurement by the EKFC equation

The serum creatinine (Scr) level was automatically measured by the enzymatic IDMS-traceable method using a biochemical analyzer (AU5821, Beckman, USA).

The EKFC equation is shown in [Citation5]. Normalized serum creatinine (Scr/Q) was mathematically obtained, where Q is the median Scr from healthy populations accounting for age and sex.

Table 1. EKFC equation.

GFR measurement by the Asian-modified CKD-EPI equation (GFRCKD-EPI)

The Asian-modified CKD-EPI equation is shown in [Citation4].

Table 2. Asian-modified CKD-EPI equation.

Statistical analysis

Normally distributed continuous variables were described as the mean ± standard deviation (SD); otherwise, they were described as the median and interquartile range (P25P75). Categorical variables were described as frequencies and percentages (%).

The relationship between GFREKFC/GFRCKD-EPI and mGFR was assessed with Spearman correlation analysis and the linear regression method. The Bland–Altman method was used to evaluate the degree of agreement between GFREKFC/GFRCKD-EPI and mGFR. The comparative performance indicators of GFR estimation for the EKFC equation and the Asian-modified CKD-EPI equation included median bias, precision, accuracy (P30), and incorrect reclassification index. The equation that met the following targets was considered superior: (1) median bias between ± 3 mL/min/1.73 m2; (2) P30 > 75%, which is considered sufficient for good clinical decision-making according to the 2002 K/DOQI benchmark [Citation13]; and (3) better precision and 95% limits of agreement in Bland–Altman analysis. Bias and precision were defined as the difference and the interquartile range (IQR) of GFREKFC/GFRCKD-EPI minus mGFR, respectively. The percentage of GFR results within 30% deviation of mGFR (P30) was used as accuracy. The respective 95% confidence intervals (95% CIs) were calculated by means of bootstrap methods (2000 bootstraps) [Citation16].

According to the patients’ basic characteristics and serum creatine values, we calculated the bias of the GFR of the EKFC and Asian-modified CKD-EPI equations using the 99mTc-diethylenetriaminepentaacetic acid (99mTc-DTPA) dual plasma sample clearance method as a reference method. We defined values more than 2 times the IQR of all the bias data of the two equations as outliers. We calculated the incorrect reclassification index, which was the total percentage of patients incorrectly reclassified into a different CKD stage by the EKFC equation and the Asian-modified CKD-EPI equation. The Wilcoxon signed rank test was performed to compare the median bias between the two models, the bootstrap method was used for precision comparison, and the McNemar test was used for comparison of P30 and the incorrect reclassification index. All statistical analyses were performed using IBM SPSS 21.0 (IBM Corp., Armonk, NY, USA), MATLAB software (version 2020 b, MathWorks) and the MedCalc application (version 4.3, MedCalc software, Mariekerke, Belgium). The p value was two-sided, and p < 0.05 was considered to be statistically significant.

Results

Characteristics of the study populations

We collected a total of 192 CKD patients with GFR estimated by the 99mTc-DTPA dual plasma sample clearance method and recollected serum samples for serum creatinine examination. However, 32 patients were excluded, 7 patients lacking Scr data due to reluctance to have their blood collected again, 8 patients undergoing dialysis, 3 patients taking drugs affecting the serum creatine value, 4 patients with edema and cardiac insufficiency, and 10 patients determined to be outliers after outlier analysis. In total, 160 patients were enrolled in our study cohort; 52 patients had chronic glomerulonephritis, 36 patients had diabetic nephropathy, 30 patients had chronic pyelonephritis, 21 patients had hypertensive nephropathy, and the remaining 21 patients had CKD due to other causes or unknown causes. The basic characteristics of the patients are shown in and .

Figure 1. Strobe diagram of the included patients.

Figure 1. Strobe diagram of the included patients.

Table 3. Basic characteristics of the study populations.

Table 4. Performance comparison of the EKFC and Asian-modified CKD-EPI equations.

Clinical assessment of the EKFC equation

GFREKFC was highly related to mGFR, with a correlation coefficient of 0.95 (95% CI, 0.93–0.96) and a regression equation of GFREKFC=mGFR × 0.87 + 5.27 (). The Bland–Altman plot showed that the 95% limits of agreement for the EKFC equation were −23.2 to 19.2 mL/min/1.73 m2 ().

Figure 2. Scatter plots and regression equation of GFREKFC and mGFR (mL/min/1.73 m2). mGFR is located on the X axis, and GFREKFC is located on the Y axis. The solid blue line represents the regression line between GFREKFC and mGFR, and dashed red lines represent 95% confidence intervals for the regression line. The solid red line represents the identity line of y = x.

Figure 2. Scatter plots and regression equation of GFREKFC and mGFR (mL/min/1.73 m2). mGFR is located on the X axis, and GFREKFC is located on the Y axis. The solid blue line represents the regression line between GFREKFC and mGFR, and dashed red lines represent 95% confidence intervals for the regression line. The solid red line represents the identity line of y = x.

Figure 3. Bland–Altman plot of GFREKFC and mGFR (mL/min/1.73 m2). The mean of mGFR plus GFREKFC is located on the X axis, and the value of mGFR minus GFREKFC is located on the Y axis. The solid blue line represents the mean difference between methods, the dashed dark red lines represent 95% limits of agreement of the mean difference between methods, the dotted red line represents the regression line of the difference between methods against mGFR, and the solid orange lines represent 95% confidence intervals for the regression line.

Figure 3. Bland–Altman plot of GFREKFC and mGFR (mL/min/1.73 m2). The mean of mGFR plus GFREKFC is located on the X axis, and the value of mGFR minus GFREKFC is located on the Y axis. The solid blue line represents the mean difference between methods, the dashed dark red lines represent 95% limits of agreement of the mean difference between methods, the dotted red line represents the regression line of the difference between methods against mGFR, and the solid orange lines represent 95% confidence intervals for the regression line.

Performance comparison of the EKFC equation and the Asian-modified CKD-EPI equation for GFR estimation

The linear relationship between GFRCKD-EPI and mGFR was GFRCKD-EPI=0.98 × mGFR +3.75. The Bland–Altman plot demonstrated that the 95% limits of agreement for the Asian-modified CKD-EPI equation were −18.5 to 23.9 mL/min/1.73 m2 (). The 95% limits of agreement in Bland–Altman analysis of the EKFC and Asian-modified CKD-EPI equations were almost the same (42.4 vs. 44.4 mL/min/1.73 m2). However, compared with the EKFC equation, the slope of the regression equation of the Asian-modified CKD-EPI equation was closer to 1, and the intercept approached zero.

Figure 4. Bland–Altman plot of GFRCKD-EPI and mGFR (mL/min/1.73 m2). The mean mGFR plus GFRCKD-EPI is located on the X axis, and the value of mGFR minus GFRCKD-EPI is located on the Y axis. The solid blue line represents the mean difference between methods, the dashed dark red lines represent 95% limits of agreement of the mean difference between methods, the dotted red line represents the regression line of the difference between methods against mGFR, and the solid orange lines represent 95% confidence intervals for the regression line.

Figure 4. Bland–Altman plot of GFRCKD-EPI and mGFR (mL/min/1.73 m2). The mean mGFR plus GFRCKD-EPI is located on the X axis, and the value of mGFR minus GFRCKD-EPI is located on the Y axis. The solid blue line represents the mean difference between methods, the dashed dark red lines represent 95% limits of agreement of the mean difference between methods, the dotted red line represents the regression line of the difference between methods against mGFR, and the solid orange lines represent 95% confidence intervals for the regression line.

Compared with the Asian-modified CKD-EPI equation, the EKFC equation demonstrated a wider median bias (–1.64 vs. 0.84 mL/min/1.73 m2, p < 0.01), which was still within 3 mL/min/1.73 m2 and was clinically not meaningful, although it was statistically significant. Furthermore, precision (12.69 vs. 12.72 mL/min/1.73 m2, p = 0.42) and the incorrect reclassification index in the different stages of CKD between the two models were not significantly different. However, the EKFC equation had a slightly better P30 (80.0% vs. 74.4%, p = 0.01).

In the mGFR ≤ 60 mL/min/1.73 m2 subgroup analysis, the precision (9.78 vs. 10.25 mL/min/1.73 m2, p < 0.01) and P30 (68.1% vs. 59.3%, p = 0.02) of the EKFC equation were better than those of the Asian-modified CKD-EPI equation, but the median bias was not significantly different (1.48 vs. −0.22 mL/min/1.73 m2, p = 0.18). However, in the mGFR > 60 mL/min/1.73 m2 subgroup, the EKFC equation did not perform better than the Asian-modified CKD-EPI equation, with a wider median bias (–8.60 vs. −0.42 mL/min/1.73 m2, p < 0.01), an inferior precision (16.39 vs. 15.86 mL/min/1.73 m2, p < 0.01), and a statistically nonsignificant P30 (95.7% vs. 94.2%, p = 1.00).

Furthermore, we compared the performance of the two equations in males and females and found that the results of the two equations were not different, with a nonmeaningful median bias and insignificantly different P30. However, the precision of the Asian-modified CKD-EPI equation was better in females than males, and the EKFC equation was superior in males. The performance comparison between the EKFC equation and the Asian-modified CKD-EPI equation is summarized in and .

Table 5. Incorrect reclassification results between the two models.

Discussion

The Kidney Disease Outcomes Quality Initiative (K/DOQI) guidelines from 2002 clearly indicated that the patient’s GFR was the primary determinant of the diagnostic rationale and staging criteria for CKD patients. Several formulas in current clinical practice were constructed to estimate GFR. In 2021, Pottel et al. developed and validated the EKFC equation, which was a modified full age spectrum, Scr-based GFR estimation equation. How accurate the EKFC equation is in external populations was the main purpose of our study.

In the current study, we assessed the clinical utility of the EKFC equation and demonstrated that it was acceptable in clinical practice. GFREKFC was highly related to mGFR; the regression coefficient was near 1, and the intercept was close to zero. Furthermore, the Bland–Altman plot showed that the 95% limits of agreement for the EKFC equation were satisfactory in our external population. As a full age spectrum equation, the new EKFC equation provides continuity across the entire age range and facilitates GFR reporting in the clinical laboratory by simple calculation. Compared with the Asian-modified CKD-EPI equation, the EKFC equation demonstrated a wider median bias within the range of ±3 mL/min/1.73 m2, and the difference had no clinical value. Despite a slightly better P30 of 80%, the EKFC equation demonstrated nearly identical precision and was not significantly different in the incorrect reclassification index compared with the Asian-modified CKD-EPI equation. In total, there is no clinically meaningful difference in the performance of the two equations in our external Chinese chronic kidney disease population. In the mGFR ≤ 60 mL/min/1.73 m2 subgroup, the EKFC equation produced excellent results. However, its performance was relatively disappointing in the mGFR > 60 mL/min/1.73 m2 subgroup, with a wider median bias and inferior precision compared with the Asian-modified CKD-EPI formula in our external CKD patients.

Furthermore, we compared the performance of the two equations in males and females and found that there were no differences in sex. The reason for this could be that the EKFC equation used Q values that were independently different for males and females, and the Asian-modified equation was developed for females and males separately.

Overall, the EKFC equation was clinically acceptable in our external Chinese CKD patient cohort. The overall bias of −1.64 mL/min/1.73 m2 was not different from the bias of −0.9 mL/min/1.73 m2 that was reported in the external validation dataset in the original study. However, the P30 of 80% in our study was lower than the overall 88% in the original study. The reason why the P30 of the EKFC equation was worse in our external population may be because the EKFC equation was developed based on a cross-sectional analysis of pooled data, from which many algorithms were developed by different reference methods, for example, 99mTc-DTPA renal dynamic imaging. In fact, a large number of studies have shown that 99mTc-DTPA renal dynamic imaging for GFR is not ideal and should not be used as a reference standard [Citation17] because it may introduce errors in the process of developing and assessing the EKFC equation. In our validation cohort, we used the 99mTc-DTPA dual plasma clearance method as the reference standard, which was accepted as the standard GFR method at the 21st International Annual Meeting of the Society of Nuclear Medicine Europe. The utilization of a different gold standard may have caused the different validation results in our external cohort. Additionally, in the mGFR > 60 mL/min/1.73 m2 subgroup, the EKFC equation had a wider median bias than those of the Asian-modified CKD-EPI equation and the original external validation dataset. We found that the median bias of the EKFC equation was higher in the mGFR > 60 mL/min/1.73 m2 subgroup than in the mGFR ≤ 60 mL/min/1.73 m2 subgroup. It is possible that the higher the mGFR is, the more obvious the underestimation of the EKFC equation, as shown in . Furthermore, the EKFC equation was developed and validated in white populations. Our study population included participants of Asian ethnicity; the median Scr concentration is different among the various ethnicities, which leads to a different Q value, where Q represents the median creatine value of healthy persons, and thus distinct GFR according to the EKFC equation. In future work, we will calculate the Q value according to the Asian median Scr concentration and further validate the EKFC equation with a larger sample size.

Our study has some limitations. First, our population’s sample size was not sufficient. Second, our study included Asian adult CKD patients; therefore, we could not determine the accuracy of the EKFC equation in children’s populations. More external validating studies are needed to verify the clinical application of the EKFC equation.

Conclusions

Overall, for adult CKD patients, there is no clinically meaningful difference in the performance of the two equations within the limits imposed by the small sample size. Clinicians should assess the GFR of adult CKD patients based on the patient’s actual situation in combination with the equation in clinical practice. For adolescents and children with CKD, more external validation studies are needed to assess the accuracy of the EKFC equation.

Author Contributions

L.Z. and H.-l.L. analyzed the data and wrote the main manuscript. H.-j.L. and J.-M. collected basic characteristics, the serum specimen of study subjects and performed the experimental operation. W.L. conducted the document retrieval and data entry. J.-m.H. and L.-g.W. performed the 99mTc-diethylenetriaminepentaacetic acid (99mTc-DTPA) dual plasma sample experiment. P.X. designed the experiments and revised the manuscript.

Acknowledgement

This manuscript has been preprinted at the Research Square, Version 1, posted 15 March 2022. https://doi.org/10.21203/rs.3.rs-1437549/v1

Disclosure statement

No potential conflict of interest was reported by the author(s).

Data availability statement

The datasets used and/or analyzed during the current study are available from the correspondence author on reasonable request.

Additional information

Funding

This work was partially supported by the Natural Science Foundation of Hebei Province under Grant number H2020206422.

References

  • Liu X, Li N, Lv L, et al. Improving precision of glomerular filtration rate estimating model by ensemble learning. J Transl Med. 2017;15(1):231.
  • Levey A, Inker S, Lesley A. GFR as the "gold standard": estimated, measured, and true. Am J Kidney Dis. 2016;67(1):9–12.
  • Pottel H, Hoste L, Dubourg L, et al. An estimated glomerular filtration rate equation for the full age spectrum. Nephrol Dial Transplant. 2016;31(5):798–806.
  • Stevens LA, Claybon MA, Schmid CH, et al. Evaluation of the chronic kidney disease epidemiology collaboration equation for estimating the glomerular filtration rate in multiple ethnicities. Kidney Int. 2011;79(5):555–562.
  • Pottel H, Björk J, Courbebaisse M, et al. Development and validation of a modified full age spectrum creatinine-based equation to estimate glomerular filtration rate: a cross-sectional analysis of pooled data. Ann Intern Med. 2021;174(2):183–191.
  • Xie P, Li HL, Huang JM, et al. Validation of the full-age spectrum equation in the approximation of glomerular filtration rate in Chinese patients with chronic kidney disease. Ren Fail. 2019;41(1):467–472.
  • Xie P, Huang JM, Lin HY, et al. CDK-EPI equation may be the most proper formula based on creatinine in determining glomerular filtration rate in Chinese patients with chronic kidney disease. Int Urol Nephrol. 2013;45(4):1057–1064.
  • Xie P, Huang JM, Li Y, et al. The modified CKD-EPI equation may be not more accurate than CKD-EPI equation in determining glomerular filtration rate in Chinese patients with chronic kidney disease. J Nephrol. 2017;30(3):397–402.
  • Wang J, Xie P, Huang JM, et al. The new Asian modified CKD-EPI equation leads to more accurate GFR estimation in Chinese patients with CKD. Int Urol Nephrol. 2016;48(12):2077–2081.
  • Andrew SL, Lesley AS, Christopher HS, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–612.
  • Andrassy KM. Comments on KDIGO 2012 clinical practice guideline for the evaluation and management of chronic kidney disease. Kidney Int. 2013;84(3):622–623.
  • Zhao L, Zhang J, J, Tian X, et al. The ensemble learning model is not better than the Asian modified CKD-EPI equation for glomerular filtration rate estimation in Chinese CKD patients in the external validation study. BMC Nephrol. 2021;22(1):372.
  • Kidney FN. K/doqi clinical practice guidelines for chronic kidney disease: evaluation, classification, and stratification. Am J Kidney Dis. 2002;39: S1–S266.
  • Bröchner-Mortensen J. A simple method for the determination of glomerular filtration rate. Scand J Clin Lab Invest. 1972;30(3):271–274.
  • Haycock GB, Schwartz GJ, Wisotsky DH. Geometric method for measuring body surface area: a height-weight formula validated in infants, children, and adults. J Pediatr. 1978;93(1):62–66.
  • Carpenter J, Bithell J. Bootstrap confidence intervals: when, which, what? A practical guide for medical statisticians. Statist Med. 2000;19(9):1141–1164.
  • Xie P, Huang JM, Liu XM, et al. 99mTc-DTPA renal dynamic imaging method may be unsuitable to be used as the reference method in investigating the validity of CDK-EPI equation for determining glomerular filtration rate. PLoS One. 2013;8(5):e62328.