We have been committed to evaluating the most accurate and clinically practical equation for the Chinese population out of the internationally popular GFR estimation equations in recent years. The results of our study showed that the FASScr−cysC equation was most suitable for estimating GFR in Chinese adults because of its superior accuracy (P30 of 57.5%), excellent concordance with mGFR (ICC of 0.921) and good diagnostic accuracy (ROCAUC of 0.951, Kappa value was 0.364) for all subjects. Moreover, its advantage that it can be used for the full age spectrum increases its clinical practicality. However, we still found that the CKD-EPIcysC equation had advantages when applied to individuals with moderately and severely reduced GFR (GFR < 60 ml/min/1.73 m2) and adults over 70 years, as its low bias (-2.23), high accuracy (57.4% and 59.4%), and lowest mean bias (7.46) were shown in Bland-Altman plots. Otherwise, the XiangYa equation was confirmed to be appropriate for individuals with slightly reduced GFR (GFR ≥ 60 ml/min/1.73 m2) due to its low bias (-5.79), high precision (18.21), high accuracy (86.2%), and very good concordance with mGFR (ICC of 0.912).
The FAS equation was developed on the basis of the average GFR of the human body before the age of 40 being approximately 107.3 ml/min/1.73 m2. After 40 years of age, the GFR decreases at the level of 1 ml/min/1.73 m2 per year, and the rate of decline is approximately 0.988[26]. Chai L et al performed the first study to suggest the applicability of the FAS equation in Chinese CKD patients and showed that the FASScr equation had lower bias (1.28) and higher accuracy (P30 of 63.64, RMSE was 19.49 in adults under 70 years and 14.06 in elderly individuals over 70 years) than the CKD-EPI equation or the MDRD equation in a single-centre study of 396 subjects in eastern China [27]. Another single-centre study of 162 subjects conducted by Xie P et al in northern China found that the FASScr equation had excellent concordance with the result measured by the 99mTc-DTPA dual plasma sample clearance method (mGFR = -0.374 + 1.029eGFR (p < 0001)), even performing better than the results measured by the 99mTc-DTPA renal dynamic imaging method (bias-1.22 vs 8.92, precision 15.69 vs 18.36, P30 75.31% vs 59.26%) [28]. A Korean study by Jeong TD et al compared the FASScr equation with the CKD-EPIScr equation in a retrospective study of 1312 patients and found that the performance of the Korean FASScr equation was equivalent to that of the CKD-EPIScr equation with low bias (-0.2) and high accuracy (P30 of 75.8%, RMSE was 15.8) [20]. However, we did not find that the Korean calibrated FASScr equation performed better than the primitive FASScr equation in our study of a Chinese population (bias − 10.37 vs -10.16, IQR 21.38 vs 20.51, P30 50.4% vs 51.4%, RMSE 24.88 vs 23,82). Yong Z et al first conducted a multicentre validation study comparing the FAS and CKD-EPI equations in Chinese adults and verified that the FASScr−cysC equation performed favourably with the lowest bias (-2.87), highest precision (IQR of 19.01), highest accuracy (P30 of 74.16%, RMSE was 17.84) and an ROCAUC of 0.953. Furthermore, the FASScr−cysC equation exhibited an absolute advantage in the elderly population over 60 years with the highest accuracy (P30 of 70.37%, RMSE was 15.21)[29]. In our study, we found that the FASScr−cysC equation performed second only to the CKD-EPIcysC equation in individuals with mGFR < 60 ml/min/1.73 m2 or adults over 70 years. However, we also obtained similar results in that the FASScr−cysC equation performed best in all subjects and had good diagnostic accuracy. Several factors may explain the conclusion. First, the FAS equation was developed for a Caucasian population. The Q value in the FAS equation is mainly based on the characteristics of the white Belgian growth curve. The nutritional status and dietary structure of Chinese people are differentiated from those of European or American white or black individuals[30]. The ideal Q value that is applicable should match the different characteristics of the Chinese population[11]. Second, the FAS combined equation assigns Scr/QScr and cysC/QcysC according to weight. The weights of the two factors in this study each account for 50%. In fact, the optimal proportion of creatinine and cystatin C for the calculation results is not yet conclusive. However, regardless of the development of the CKD-EPI equation or FAS equation, the combination of creatinine and cystatin C could reduce the imprecision of either biomarker alone[19, 31]. Third, we found that the CKD-EPIcysC equation obviates the race/ethnicity coefficient. However, the calibration coefficients of the various equations in the Chinese race are all greater than 1[32], and there is a paradox between them. Therefore, it is necessary to carry out multicentre and multiethnic equation verification for Chinese people.
Creatinine is produced by muscle catabolism and secretion by kidney tubular cells. Its concentration is greatly affected by the body's nutritional status and muscle mass. Poor nutritional status, infection along with an inflammatory state, or disability will affect the production of creatinine, leading to an inability to accurately assess renal function, especially between sexes. Cystatin C is less affected by muscle protein content and is totally reabsorbed in the proximal tubules. The changes in the concentration of cystatin C are relatively stable. A retrospective cohort study conducted by Wang Y et al with 308 patients with diabetic nephropathy in southwestern China reclassified 39% of patients with poor prognosis in CKD stages 1–2 to CKD stages 3–5 by the CKD-EPIcysC equation compared with the results of the CKD-EPIScr equation (the median survival of the reclassified and not reclassified groups was 52 and 94 months, respectively)[33], which indicated that the CKD-EPIcysC equation could sensitively detect kidney injury, especially in individuals with moderately and severely reduced GFR; the same conclusion was reached in our study. Yang M et al performed a retrospective study of 632 CKD patients from two general hospitals in southeastern China and verified that the cystatin C-based equation had an advantage over the creatine-based equation, especially in patients with CKD stages 3–5; the CKD-EPIcysC equation had the lowest bias (-4,10), best precision (IQR of 17.35), and best accuracy (P30 of 50.5%, RMSE was 0.93) when GFR was < 60 ml/min/1.73 m2[34]. These results were similar to those in our study, proving the advantage of the CKD-EPIcysC equation in Chinese individuals with moderately and severely reduced GFR. Moreover, a cohort study of 70 patients who underwent kidney transplantation conducted in Korea revealed that the ROCAUC for cystatin C at the cut-off value of 45 ml/min/1.73 m2 was 0.800 (that for creatinine was 0.763), and the best predictive value for cystatin C was 1.27 mg/L, with a sensitivity of 77.8% and specificity of 78.8%[35]. The study also showed the diagnostic advantage for the cystatin C-based equation in predicting renal function, as GFR < 45 ml/min/1.73 m2 may indicate significant kidney injury among transplant recipients.
In this study, we also found that the CKD-EPIcysC equation performed well in elderly individuals over 70 years of age. According to the law of GFR decline with age, whether GFR < 60 ml/min/1.73 m2 is used as the diagnostic criterion for CKD in the elderly population is controversial, which was proposed in the BIS equation study[17]. Epidemiological studies have shown that elderly patients over 65 years old account for 53.07% of the CKD population in China[36]. A single-centre study conducted by Ye X et al among Chinese elderly CKD patients aged ≥ 60 years showed that the CKD-EPIcysC equation did not perform well (bias of -9.05, IQR of 19.61, P30 of 71.74%, RMSE was 18.50). However, the BIS-1Scr and BIS-2Scr − cysC equation performed better, with lower bias (-5.20 and − 8.65), better precision (IQR of 21.58 and 16.31), and higher accuracy (P30 of 79.1%, RMSE was 17.20 and 16.94) [37]. Another cross-sectional study of 218 elderly patients over 75 years of age conducted by Changjie G et al in southern China found that the BIS-2Scr − cysC equation performed best (bias of 0.63, IQR of 4.36, P30 of 94.50%, RMSE was 7.21) rather than the CKD-EPIcysC equation (bias of 3.86, IQR of 17.07, P30 of 67.89%, RMSE was 17.22)[38]. These results both showed the opposite conclusion compared to the conclusion of our study. We found that the FASScr−cysC equation performed equivalently to the BIS-2Scr − cysC equation but had slightly more bias in our study. However, considering the constitution of the equation, it was convenient to use the FAS equation in clinical practice. As a special population, the elderly population often combines various complications other than nephropathy with decreases in nutritional status and immune function. For cardiovascular disease, cerebrovascular diseases, and people who take vasoactive drugs, the applicability of eGFR equations remains to be further verified.
The MDRD equation performed best among the equations based on creatinine alone in all subjects, but the accuracy was significantly lower when GFR was < 60 ml/min/1.73 m2, which is consistent with the conclusion that the accuracy of the MDRD equation decreased when renal function was moderately reduced, as previously reported by Murata K et al. Even though the CKD-EPIScr equation enhanced the specificity of the detection of GFR < 60 ml/min/1.73 m2 compared with that of the MDRD Eq. (98% vs 94%), at the cost of reducing the sensitivity (50% vs 70%), the creatinine-based equation was still maligned in CKD patients [39]. The XiangYa equation is a new creatinine-based equation developed in recent years based on the central Chinese population. This equation has been externally verified in a Han population in central China and the Uighur population in western China. To our knowledge, we performed the first study to validate the XiangYa equation with four other main equations together in one centre. Surprisingly, we found that the XiangYa equation performed excellently in individuals with a GFR ≥ 60 ml/min/1.73 m2 and had a good ICC value with mGFR. However, we also made a similar conclusion that it performed terribly in individuals with GFR < 60 ml/min/1.73 m2 and the elderly population. This issue may be related to the original development subjects of the XiangYa equation causing significant bias, as the mean age was 52.34 with an SD of 13.23, and the mean mGFR was 71.320 ml/min/1.73 m2 with an SD of 23.96[21].
In this study, the accuracy P30 of each equation was found to be less than 70%. This is because the plasma biochemical indexes included in this study were from 48 hours before kidney nuclear radiograph, and it was not necessarily the same result as in the morning. This was in line with the current procedures of most hospitals in China. Due to the large population in China, the imbalance between doctors and patients, and limited medical resources, biochemical blood tests and imaging examinations are often unable to be synchronized. Using biochemical test results from within 48 hours to assess renal function is not unacceptable for CKD. In recent reports, the rate of change in biomarkers was shown to be smaller than that in kidney nuclear radiography. The eGFR computed by equations is at least as reliable as the mGFR measured by nuclear medicine for monitoring patients over time[10]. The comparative verification data adopted in this study are more applicable to China's national conditions, and its statistical results have more practical clinical significance.
We verified the effect of the current popular equation on the evaluation of renal function in adult patients in central China. We must take into account the limitations of this study. First, we only enrolled Chinese Han individuals in this single-centre experiment. China is a multi-ethnic and multi-regional country. Different living customs lead to diverse human physiques, and different medical conditions contribute to inconsistent diagnoses of CKD in various regions. We will consider involve more ethnic groups in our multi-center verification research. Second, the method of assaying serum creatinine and cystatin C in this study may be different from those used in the derivation studies for the above equations, so was the 99mTc-DTPA used for mGFR measurement. Although the difference may exist in assaying equipment and materials, the current assaying level can be traced back to international certification standards, as a consequence its results are acceptable. Systematic assessment can be performed to evaluate the influence of different assaying methods on the performance of eGFR equation forward. Third, this was a retrospective validation study, and the concentration of blood biochemical indexes came from one-time testing. We were not able to evaluate the accuracy of eGFR assessed by each equation in terms of a dynamic change in renal function. Considering that the renal function of CKD patients does not change significantly in the short term, it is advisable to evaluate the calculation results of patients in hospitalization and discharging together for a comprehensive analysis in the future study.