- Research article
- Open Access
- Open Peer Review
Risk scores to predict decreased glomerular filtration rate at 10 years in an Asian general population
BMC Nephrologyvolume 18, Article number: 240 (2017)
Asians have among the highest prevalence of chronic kidney disease (CKD) or end-stage renal disease in the world. A risk score capable of identifying high risk individuals at the primary care level could allow targeted therapy to prevent future development of CKD. Risk scores for new CKD have been developed in US general populations, but the impact of various risks factors for development of CKD may differ in Asian subjects. In this study, we aimed to develop risk models and simplified risk scores to predict the development of decreased glomerular filtration rate (GFR) at 10 years in an Asian general population using readily obtainable clinical and laboratory parameters.
Employees of EGAT (The Electric Generating Authority of Thailand) were studied prospectively. Multivariable logistic regression models were used to assess risk factors and used to derive risk models and risk scores for developing decreased GFR at 10 years: Model 1 (Clinical only), Model 2 (Clinical + Limited laboratory tests), and Model 3 (Clinical + Full laboratory tests). The performance of the risk models or risk scores to predict incident cases with decreased GFR were evaluated by tests of calibration and discrimination.
Of 3186 subjects with preserved GFR (eGFR ≥60) at baseline, 271 (8.5%) developed decreased GFR (eGFR < 60) at 10 years. Model 1 (Age, sex, systolic blood pressure, history of diabetes, and waist circumference) had good performance (χ2 = 9.02; AUC = 0.72). Model 2 (Age, Sex, systolic blood pressure, diabetes, glomerular filtration rate) had better discrimination (χ2 = 10.87, AUC = 0.79) than Model 1. Model 3 (Model 2+ Uric acid, Hemoglobin) did not provide significant improvement over Model 2. Based on these findings, simplified categorical risk scores were developed for Models 1 and 2.
Clinical or combined clinical and laboratory risk models or risk scores using tests readily available in a resource-limited setting had good accuracy and discrimination power to estimate the 10-year probability of developing decreased GFR in a Thai general population. The benefits of the risk scores in identifying high risk individuals in the Thai or other Asian communities for special intervention requires further studies.
Chronic kidney disease (CKD) and decreased glomerular filtration rate (GFR) are both associated with elevated risks for end-stage renal disease (ESRD), cardiovascular disease, and death [1, 2]. The incidence of CKD has increased worldwide with important public health and economic implications especially in developing countries where resources are limited. Despite the high prevalence rate, CKD awareness rates are often very low in the community as early CKD is usually asymptomatic [3, 4]. As a consequence, CKD is not frequently detected until it has already advanced, and opportunities for intervention are lost. A risk score that identifies those at higher risk for future CKD has been proposed as a stratification and prediction device . Cardiovascular risk scores, such as the Framingham score,  have influenced public health policy in the primary prevention of cardiovascular diseases and have been tested in many populations . A simple and accurate renal risk score would lead to targeted medical management at the primary care level to individuals at the highest risk for future CKD.
Asians represent 60% of the global population and have among the highest prevalence of CKD in the world . Asian countries include developed and low to middle-income countries with different risk factors for CKD development [9, 10]. In low to middle-income countries, the burden of disease is changing from infections towards chronic lifestyle-related diseases as a result of demographic transition and urbanization. In 2011, Thailand was reclassified by the World Bank from a lower-middle income to a higher-middle-income country. Over the last decade, the numbers of patients with ESRD have increased by an order of magnitude. Data from the Thailand Renal Replacement Therapy registry reported that the numbers of patients on renal replacement increased from 30 per million people in 1997 to 1199 per million people in 2014 (file:///E:/TRT%202007–2017/1.TRT-report−2014-_3–11-59_-final_pdf). This staggering increase is both due to public health policy changes as well as due to higher rates of CKD.
Risk scores for incident CKD have been developed in the US general populations [11, 12], but there are limited information on risk models that allow long term predictions of incident CKD in Asians. Since risks factors for CKD and decreased GFR in populations from Asian countries including Thailand may differ from Caucasian populations, risk scores developed in an Asian community may be more appropriate to evaluate risks in Asian populations. In this study, we aimed to evaluate risk predictors and develop risk models and risk scores for developing decreased GFR at 10 years follow-up in a Thai general population cohort. We hypothesized that new cases of decreased GFR may be predicted by a risk score containing a subset of clinical variables. Because subjects from low to middle-income Asian countries often have limited resource and less access to routine medical checkups than their Western counterparts, we have developed a risk score based on clinical and simple laboratory parameters easily assessed in the primary care setting and compared them to a score based on more expanded, but still standard laboratory work up.
The subjects were employees of EGAT (The Electric Generating Authority of Thailand) who volunteered to participate in a health survey. All participants completed a medical evaluation and had routine laboratory investigations including urinalysis. Blood was drawn after a 12 h fast. The details of the EGAT study cohorts and the study protocols have previously been described . In summary, there were 3 EGAT cohorts in total. In 1985, 3499 workers of EGAT (half of the total employees) were randomly enrolled as EGAT 1 cohort. In 1998, 2999 employees were randomly enrolled as EGAT 2 cohort. The age range of both EGAT 1 and EGAT 2 was selected as 35–54 years old. In 2009, 2584 participants were recruited to the EGAT 3 cohort, but the age range was expanded to 25–54 years old. Participants in different EGAT cohorts did not overlap. The response rates in the first and second cohorts were >95%, but in EGAT 3 it dropped to 76%. The major reasons for this reduced response rate are thought to be that some selected participants were due to leave for periodic work in the sites outside Bangkok, and the fear of contracting swine flu in examination centers in 2009. The mean ages were not different between responders and non-responders (each was 41 years; P = 0.78), but responders were less likely to be male (74 vs 78%; P = 0.004). The EGAT 1 cohort was resurveyed in 1997, 2002, 2007, and 2012. The EGAT 2 cohort was resurveyed in 2003, 2008, and 2013. EGAT 3 was resurveyed in 2014. Each time, the same individuals were contacted by telephone and invitation letter to attend the follow-up examination, or else information about the cause of death was sought for those known to have died during the interim period. At each follow-up visit, the subjects underwent similar medical evaluations and had routine laboratory investigations as the baseline visit. The EGAT studies were conducted in accordance with the Helsinki Declarations and approved by the Ethics Committee, Ramathibodi Hospital, Mahidol University, Thailand. Written informed consent was obtained.
The Derivation cohort included subjects from EGAT 1 and EGAT 2 cohorts. The 2002–2003 period (EGAT 1 3rd examination and EGAT 2 2nd examination) was selected to be the baseline of the study to ensure a broad starting age range of approximately 40–70 years old. The outcome was evaluated 10 years later in 2012–13 (EGAT 1 5th examination and EGAT 2 4th examination). Of EGAT 1 (n = 2360) and EGAT 2 (n = 2651) participants who attended the baseline 2002–2003 examination, subjects (n = 37) with missing baseline creatinine data were excluded leaving a total of 5010 subjects with serum creatinine at baseline in 2002–2003. The Derivation cohort dataset was derived from participants with preserved GFR (estimate glomerular filtration rate (eGFR) ≥ 60 mL/min/1.73m2) at baseline who attended both the examinations in 2002–2003 and 2012–13 (EGAT 1 5rd examination and EGAT 2 4nd examination) (Fig. 1).
The Validation cohort dataset was derived from EGAT 3 cohort participants with preserved GFR (eGFR ≥ 60) at baseline in 2009 (EGAT 3 1st examination) who were followed up 5 years later in 2014 (EGAT 3 2nd examination). Participants younger than 40 years old at baseline were excluded to maintain similar cut-off age as the Derivation cohort (Fig. 1).
Incident cases with decreased eGFR (Decreased GFR) refers to subjects with preserved GFR (eGFR ≥60) at baseline who subsequently developed decreased GFR (eGFR < 60 mL/min/1.73m2) at the 10 year follow-up. The outcome is a modification from the KDIGO definition of CKD stage 3–5. The difference being that the kidney function was measured once and the presence of decreased eGFR was not confirmed to be present for greater than 3 months .
An average of two blood pressure measurements in seated position were used. Hypertension was defined as systolic blood pressure (SBP) ≥ 140 mmHg or diastolic blood pressure (DBP) ≥ 90 mmHg or use of oral antihypertensive medication . A positive history of diabetes was ascertained from subject testimony or the use of medications. Diabetes mellitus (DM) was defined as a fasting glucose of ≥126 mg/dl or a positive history of diabetes  Smokers were defined according to current status within the last 12 months. History of cardiovascular disease included subjects with previous peripheral vascular disease, heart muscle, heart attack and stroke. Waist circumference was measured midway between the lowest ribs and the iliac crest. Body mass index was defined as weight in kilograms divided by the square of height in meters. There were less than 5% missing data for the covariates and these were regarded as missing.
Laboratory measurements and GFR estimation
Cholesterol, triglycerides, and HDL cholesterol were determined using enzymatic methods as published  LDL cholesterol was calculated using the Friedewald equation . Serum creatinine (sCr) was measured by the enzymatic assay on the Vitros 350 analyzer (Ortho-Clinical Diagnostics, USA) using IDMS-Standard Reference Material (SRM) 967 as the standard . The coefficient of variations for low, high sCr were 1.64%, 0.41%. Urine protein was detected by urinalysis reagent strip (Bayer, Indiana, USA). Estimate glomerular filtration rate (eGFR) was calculated according to two-level race variable Chronic Kidney Disease–Epidemiology Collaboration (CKD-EPI) equation  using the non-black coefficient as this is now the preferred equation in the general population .
Continuous data is reported as mean (± SD). Categorical data is presented as percentages. Continuous variables were compared using Student’s t test. Categorical variable were compared using Chi-square test. The baseline data of the Derivation cohort were used to predict risk for Decreased GFR defined by CKD-EPI after excluding patients with eGFR <60 at baseline.
Multivariable logistic regression models were used to assess risk factors for Decreased GFR, selected from earlier reports [11, 21] Variables were sequentially added in a pre-specified order and incorporated using a P < .05 threshold for entry and retention in the final model. Three analysis models were used:
Model 1 contained only readily available clinical variables previously linked to CKD without the use of laboratory parameters. In order of entry, the candidates were: age, sex, history of diabetes, systolic blood pressure, waist circumference, current smoking (yes/no), and history of cardiovascular disease
Model 2 comprised of Model 1 plus baseline eGFR and blood sugar as linear parameters.
Model 3 comprised Model 2 and other biochemical parameters including uric acid, hemoglobin and HDL.
Because urine screening is not routinely performed in Thailand, we evaluated alternative models to assess the benefit of urinalysis in which dipstick proteinuria (defined as 1+ or higher) was also calculated in addition to existing models. Additional analyses were also performed in which hypertension was substituted for SBP and BMI was substituted for waist circumference.
Variables were dichotomized to enhance clinical utility. A risk scoring system for CKD was developed using previously established methods. In brief, each variable was assigned points proportional to the product of its regression coefficient from the multiple logistic regression model for CKD (described above) and the measured value of that variable [6, 22].
The performance of the risk models or risk scores to predict Decreased GFR were evaluated by tests of calibration and discrimination. Calibration measures how closely predicted outcomes agree with actual outcomes. The risk models were used to divide subjects into deciles of predicted risk for Decreased GFR. Differences of predicted and observed rate were compared using Hosmer-Lemeshow chi-square test . Chi-square values <20 generally support adequate fit. Discrimination is the ability of a prediction model to separate those who develop Decreased GFR from those who do not. We quantified this by calculating the c statistic, analogous to the area under a receiver operating characteristic curve (AUC) . This value represents an estimate of the probability that a model assigns a higher risk to those who develop Decreased GFR than to those who do not. AUC varies from 0.5 (no discrimination) to 1.0 (perfect discrimination) and is reported as AUC (95% confidence interval). AUC of different models were compared using Delong test .
As there is a day to day variation in eGFR, persons with borderline eGFR may at reevaluation have eGFR slightly lower than 60 just due to random variation. We tested the performance of the risk models after excluding Decreased GFR cases with an eGFR decline of less than 5 ml/min over 10 years as a sensitivity analysis.
We also evaluated the performance of the risk models using the four-level race variable CKD-EPI equation with Asians coefficients (×1.052) to calculate eGFR as this equation has been proposed by some investigators as a valid equation in Asian subjects .
Internal and External validation
Internal validation of the c-statistics was performed using bootstrapping with 1000 replications of individuals sampled with replacement. The Somer’D correlation was used to estimate the correlation between the observed and predicted values for Decreased GFR, called Dboot . Calibration bias of the model was assessed by subtracting the original Somer’D correlation from the mean Dboot.
External validation of the performance of the Risk Score developed in the Derivation Cohort from EGAT 1 and EGAT 2 participants for predicting Decreased GFR was evaluated in the Validation Cohort composed of EGAT 3 participants.
Analyses were performed with SPSS version 11.5 (SPSS Inc.; Chicago, IL). Comparisons of AUCs were performed with Stata 14 (Stata Corp LP, College Station, TX).
Baseline characteristics of the Derivation cohort and the incidence of Decreased GFR
Of EGAT 1 and EGAT 2 participants with serum creatinine (n = 5010) at baseline in 2002–2003, 428 subjects died, 222 retired and moved, 915 did not want to participate and four subjects had missing serum creatinine leaving 3441 participants with complete data for both 2002–2003 and 2012–2013 visits (Fig. 1). These subjects with complete follow-up data were 1 year younger, had slightly lower proportions of males, diabetes, hypertension, proteinuria compared to the total participants with serum creatinine in 2002–2003. The mean eGFR differed by 0.9 ml/min/1.73 m2. Although these differences were statistically significant, the actual differences were generally not large (Additional file 1: Table S1).
Of the 3441 subjects with complete data for both visits, subjects with eGFR < 60 at baseline (n = 255) were excluded leaving 3186 patients for analysis as the Derivation cohort. Baseline characteristics of the Derivation cohort are shown in Table 1. All subjects are Thais or Thai-Chinese. The age range was 39–71 years with distribution as follows: <45 years, 24.4%; 45–54 years; 40.9%; ≥ 55 years, 34.7%. Proterinuria (defined as 1+ or higher.) was present in 15.9% of subjects at baseline. Of the Derivation cohort, 271 (8.5%) developed Decreased GFR at follow up.
Performance of the risk prediction models
By univariate analysis, smoking and history of cardiovascular disease were not significant. The factors significant in the multivariate model for the simple clinical model were: Age, sex, systolic blood pressure, history of diabetes, waist circumference (Table 2) and the performance of Model 1 (Fig. 2) was good (χ2 = 9.02, p = 0.34; AUC = 0.72 (0.69–0.75), p < 0.001).
Substitution of waist circumference with overweight as a category (BMI ≥ 25) produced similar results (Additional file 1: Table S2) with the same factors remaining predictive (χ2 = 8.87, p = 0.35; AUC = 0.72 (0.69–0.75), p < 0.001). Substitution of hypertension for SBP also showed similar results (χ2 = 10.71, p = 0.22; AUC = 0.71 (0.68–0.74), p < 0.001).
Model 2: (Clinical + limited laboratory tests)
In model 2, the following factors were significant: Age, Sex, SBP, DM, GFR category. Waist circumference was not significant (Table 2). Substitution of waist with BMI ≥ 25 was also not significant. Model 2 (Fig. 2) had good performance (χ2 = 10.87, p = 0.21, AUC = 0.79 (0.76–0.82), p < 0.001). Substitution of SBP with hypertension also produced similar results (χ2 = 5.62, p = 0.69, AUC = 0.78 (0.76–0.81), p < 0.001).
The AUC of Model 2 was significantly different than AUC of Model 1 (p = <0.001).
Model 3: (Clinical +full laboratory tests)
Using the full laboratory model, the following parameters were significant predictors by multivariate analysis: Age, Sex, SBP, DM, GFR, Uric acid, Hemoglobin. Waist circumference and low HDL were not significant (Table 2). Model 3 (Fig. 2) performed well (χ2 = 8.28, p = 0.41, AUC = 0.80 (0.77–0.82), p < 0.001). Substitution of hypertension for SBP produced similar results (χ2 = 5.53, p = 0.70, AUC = 0.79 (0.760–0.82), p < 0.001).
The AUC of models 2 vs model 3 were not significantly different (p = 0.39).
Addition of proteinuria to risk models
When proteinuria was added to the risk factors used in Model 1, proteinuria was a significant risk factor, but waist circumference was no longer significant if proteinuria was added (Additional file 1: Table S3). The performance of Model 1+ proteinuria was good (χ2 = 6.07, p = 0.64, AUC =0.75 (0.70–0.80), p < 0.001). However, although AUC of Model 1+ proteinuria tended to be better than Model 1 without proteinuria, this was not statistically significant (p = 0.08). Proteinuria was not a significant factor when added to model 2 or 3. (Data not shown).
Exclusion of Decreased GFR subjects with eGFR decline of less than 5 ml/min over 10 years
Among 271 subjects who developed Decreased GFR, 26 cases (9.6%) had eGFR decline of less than 5 ml/min/1.73m2 over 10 years. Exclusion of these subjects did not alter the odds ratios of the risk factors in the multivariate analysis (data not shown) or alter the predictive values of the Models as follows: Model 1: AUC = 0.72 (0.69–0.75), χ2 = 4.11, p = 0.847; Model 2: AUC = 0.77 (0.74–0.80), χ2 = 11.71, p = 0.165; Model 3: AUC = 0.78 (0.75–0.81), χ2 = 6.89, p = 0.549.
Use of four-level race variable CKD-EPI equation with Asian coefficient
When the CKD-EPI Asian formula was used to calculate eGFR, the cohort consisted of 3271 subjects with preserved GFR (eGFR ≥60) at baseline. Subsequently, 204 (6.2%) developed Decreased GFR at 10 years. The performance of the models using the CKD-EPI Asian equation was comparable to the two level race-variable CKD-EPI equation using coefficient for non-blacks, although there was a slight reduction in calibration for Model 2 as follows: Model 1: AUC, 0.75 (0.72–0.79), χ2 = 3.25, p = 0.918. Model 2: AUC, 0.80 (0.76–0.83), χ2 = 23.15, p = 0.003.Model 3: AUC, 0.81 (0.78–0.84); χ2 = 6.60; p = 0.580.
Development of Risk scores
Risk algorithms for Model 1 (waist), Model 1 (BMI) and Model 2 were converted to risk scores using dichotomized values (Tables 3 and 4). Because the performance of Model 2 and Model 3 were very similar, we did not develop a risk score for model 3 given the additional complexity.
The performance of the Clinical Risk Score (derived from Model 1) were: χ2 = 9.19, p = 0.33; AUC = 0.71 (0.68–0.74), p < .001, and for the Clinical + Limited laboratory test Score (derived from Model 2) were: χ2 = 11.05, p = 0.19, AUC =0.79 (0.76–0.82), p < 0.001.
The difference in discrimination (AUC) between the Clinical and the Clinical + Limited laboratory test Score was significant (p < 0.001).
For Model 1, the calibration bias was 0.00004 (p = 0.60). The C statistics of the bootstrap validation model 1 (AUC: 0.71 (0.68–0.74)) was similar to the developed Model 1.
For Model 2, the calibration bias was 0.06 (p = 0.12). The C statistics of the bootstrap validated model 2 (AUC: 0.75 (0.72–0.78)) was similar to the developed Model 2.
From a total of 2584 participants from EGAT 3 cohort recruited in 2009, we excluded participants missing serum creatinine at baseline (n = 9) or at follow-up in 2014 (n = 14), age at baseline less than 40 years (n = 904), participants who did not attend follow-up (n = 234) or died (n = 14) or those with eGFR < 60 at baseline (n = 14). The Validation cohort comprised of 1395 subjects (Fig. 1). The baseline characteristics of the participants included in the Validation cohort are shown in Additional file 1: Table S4. At 5 years, 1.9% of the Validation cohort developed Decreased GFR. The performance of the Risk Scores for predicting Decreased GFR in the Validation cohort were as follows: Model 1 (waist): AUC, 0.66 (0.55–0.78); χ2,=4.31, p = 0.229; Model 2: AUC,0.88 (0.80–0.95); χ2,=2.29, p = 0.514.
We developed risk prediction models for developing decreased eGFR at 10 years in a middle-age to older Thai general population using standard clinical parameters and routine laboratory tests. The predictors for the clinical model were: age, sex, systolic blood pressure, waist circumference or body mass index, and history of diabetes. The predictors for clinical and limited laboratory tests comprised of age, sex, systolic blood pressure, diabetes mellitus and baseline eGFR. The risk models demonstrated good discrimination and calibration with good internal validation. The addition of more laboratory tests of hemoglogin concentration, uric acid, HDL did not increase the performance of the clinical and limited laboratory tests significantly. Based on these results, we developed 2 simplified risk scores: a clinical risk score and a combined clinical and limited laboratory risk score. External validation using a separate cohort confirmed good performances of these scores. The parameters used in the scores are readily available for self–testing or evaluation by medical personnel in the primary care settings appropriate for a resource-limited setting such as Thailand or other parts of Asia.
Improved clinical prediction is an essential component of personalized medicine. Clinical prediction tools such as the Framingham cardiovascular risk score  have helped shape public health policy in the primary prevention of cardiovascular disease in many countries. However, despite the identification of several key renal risk factors,  similarly useful risk scores for predicting long term risk of new CKD has not been developed in an Asian general population. We are aware of 3 prior published risk prediction scores for incident CKD from community-based cohorts. The first study was derived and validated using data from middle-aged and older adults in the US community . From this study, the final model included 8 variables: age, sex, anemia, hypertension, diabetes mellitus, cardiovascular disease, history of heart failure, and peripheral vascular disease. This risk score had moderate discriminatory power (c-statistic 0.70) and did not contain data on baseline GFR or proteinuria. The second study evaluated Taiwanese subjects, but was compromised by poor discriminatory power (c-statistic 0.67) and short follow-up (median 2.2 years) . Because the follow-up of this study was very short, only those subjects with very rapid decline in GFR would be detected and the cumulative effects of risk factors on CKD development would be underestimated. The most recent risk score was derived from Caucasian subjects from the Framingham cohort . This study shared several elements to our study including a similar follow-up period of 10 years with similar, but not identical risk predictors. In the Framingham study, age, diabetes and hypertension were significant predictors for the clinical score, and age diabetes hypertension, GFR, proteinuria were predictors in the combined clinical and laboratory model. Both the clinical and the combined clinical and laboratory tests had a high degree of accuracy and discriminatory power in US Caucasian or Black subjects (AUC 0.78–0.83). However, the Framingham risk score had low discriminatory power and accuracy when tested in our cohort (χ2 = 256.5, p < 0.001 and AUC 0.63 for model 2). Compared to our scoring system, age was the most significant contributor to the combined risk score in the Framingham study with diabetes, and hypertension only contributing in a minor role. In our score, both diabetes, hypertension and baseline eGFR were more important contributors to the score. In addition, being overweight was an important clinical predictor of CKD whereas this was not included in the Framingham score. The importance of obesity and diabetes highlighted in our score is especially striking given the rise in obesity and diabetes across low to middle-income populations in Asia with increasing globalization [8, 10].
Previous prediction scores for incident CKD were derived using the modification of diet in renal disease study (MDRD) equation . The MDRD formula was first developed in US patients with established CKD (6). CKD-EPI equation, which was derived from both CKD and normal subject cohorts has been shown to be more accurate than MDRD especially in subjects with preserved GFR  Furthermore, CKD-EPI has been shown to be superior to MDRD at predicting adverse outcomes and improved the accuracy in outcome prediction in Caucasian and Black US subjects  There is considerable controversy on the optimum eGFR equation in Asian populations [30, 31]. There are as much as 20–30 ml/min/1.73m2 differences in GFR estimates between various Asian formulae. These discrepancies results in as much as 10 fold variations in CKD prevalence rate, and alter the prognostic significance attributable to the presence of CKD . Differences in the reference GFR methods, and the proportion of non-CKD subjects in the development cohort likely account for these discrepancies as much as any biological differences between Asian subjects of various ethnicities [30, 31]. A Thai eGFR equation has been developed with Thai patients with established CKD  using a short plasma clearance of 99Tc DTPA as the GFR measurement method. Given the lack of inclusion of normal subjects in the development cohort and methodological issues used to develop the Thai eGFR equation , we elected to use CKD-EPI for the sake of generalizability of the score to other Asian cohorts and for comparisons with other global populations [34, 35]. The rationale of our choice is supported by the fact that the CKD-EPI equation-based CKD staging has been shown to result in similar risk predictions for adverse outcomes in Asians, Whites and Blacks in a large meta-analysis . In addition, we also tested the performance of our risk scores using the Asian coefficient of the four-level race variable CKD-EPI equation to calculate eGFR. Although this equation was developed in Asian populations, its role remain uncertain as the accuracy can vary in different Asian populations  Changing to the CKD-EPI Asian equation resulted in lower incident cases with Deceased GFR, but the performance of risk models were largely similar to using the original CKD-EPI equation. This suggests that the risk scores may be used to predict Decreased GFR when the Asian coefficient for CKD-EPI is used to calculate eGFR but the performance might be slightly reduced.
There are several potential implications of this work. First, by allowing physicians to determine an individual’s estimated risk for Decreased GFR, the score may inform clinical decision-making, for example to modify treatment, frequency of follow-up or institute renal primary prevention measures in high risk patients. Secondly, the use of the score may raise the profile of kidney disease among the general population, a key goal as the current CKD awareness rates is only 1.9% in Thailand . Thirdly, it is noteworthy that the discrimination of the clinical risk score is already fairly high and this score could be used for focused renal screening, identifying individuals in whom creatinine measurement would be most cost-effective. Of note, although proteinuria tended to improve the discrimination of the clinical model, the improvement was slight and not statistically significant. Given the increased cost, our study suggests that routine population-based dipstick testing proteinuria may not be worthwhile. Of course, the risk score should not be used as a substitute for established urinalysis-screening intervals in people with diabetes or have other high risks. Finally, our score may be useful in estimating the individual risk and future prevalence of CKD in middle age to older subjects from other Asian general population [8, 10]. The risk factors used in our score such as older age, diabetes and hypertension are universal risk factors for CKD . Many low to middle-income Asian countries are exposed to similar health impact of globalization as Thailand, and share similar prevalence of many CKD risk factors that may be considerably different from the West . Combined with a closer genetic background, risk scores developed in one Asian population may be more accurate at predicting CKD in another Asian population. Although our scores have been validated externally using a separate cohort, the cohort used for validation consisted of younger Thai subjects from a similar employment background as the Derivation cohort. Participants also had shorter follow-up period of five years and a lower incidence of Decreased GFR. Further studies in diverse Asian cohorts with longer follow-up duration are necessary to confirm the usefulness of our score in predicting the long term risk of CKD in other Asian populations.
There are several strengths to this study. To our knowledge, this is the first prospective risk score to predict incident cases of Decreased GFR with follow-up of up to 10 years in an Asian population. We employed a community-based cohort with detailed assessment of risk factors and standardized calibrated enzymatic creatinine measurements. Several limitations also should be acknowledged. Baseline and follow-up creatinine were measured on a single occasion. According to KDIGO guidelines, the diagnosis of CKD requires two estimates of GFR separated by 3 months . As such, the outcome we evaluated in this study does not fulfill the criteria of incident CKD, but rather, the outcome represented incident cases with decreased eGFR (eGFR < 60). This is a limitation our study shares with most published studies involved in developing incident CKD score including the Framingham heart study . Multiple measurements in cohort studies are costly to perform especially in a resource-limited setting such as ours. By measuring the follow-up eGFR only once, we cannot exclude the fact that some subjects may have reversible acute kidney injury rather than persistent CKD. In addition, some subjects with a rather low borderline eGFR may have a follow-up eGFR slightly lower than 60 just due to random variation of serum creatinine. In practice, the EGAT subjects who attended the follow-up examination were not acutely ill and significant acute kidney injury was probably not frequent. Exclusion of subjects with less than 5 ml/min/1.73 m2 change in eGFR in our sensitivity analysis did not alter the performance of the risk score. Thus we expect that the risk factors identified in our study and our risk score should be valuable in identifying subjects at risk of developing incident CKD in the general population. Nonetheless, a single measurement of eGFR may lead to an overestimation of incident CKD. Future studies with repeated creatinine measurements that can confirm the presence of decreased eGFR after 3 months should provide a more accurate prediction of risks of developing incident CKD, although such a study would be more expensive to perform.
We used dipstick proteinuria rather than urine albumin creatinine ratio. It is possible that quantitative proteinuria measurement might have provided better prediction for Decreased GFR and improved our prediction models. The aim of this study was to devise scores to screen subjects at risk of Decreased GFR in a resource-limited setting and urine albumin creatinine ratio was not performed at baseline because of the higher expense. The KDIGO 2012 guidelines  suggested that urine dipstick might be substituted for albumin creatinine ratio when the latter is not available. In other scores e.g. Framingham, substitution of urine dipstick for urine albumin creatinine ratio did not alter the results .
A number of the participants who attended the visit in 2002–2003 did not attend the follow-up visit in 2012–13. It is not surprising that the health risk profile of those who attended both visits were statistically better than the total study population at baseline since death or retirement were common reasons for non-attendance. Although the differences in these risk factors were statistically significant, the actual differences were clinically quite small for most variables. Nonetheless, it is possible that the subjects who did not attend the follow-up examination were sicker and the true incidence of Decreased GFR might have been underestimated. Although the EGAT cohort is a community-based cohort, there may be some differences in the participant profiles from the Thai population as a whole. All participants were Thais and Thai-Chinese who represent the vast majority (over 95%) of the Thai Census population. Our study included only middle-age to older subjects and had a higher percentage of males than females compared to the total Thai population. EGAT employees come from all regions of Thailand and cover a wide-range of sociodemographic backgrounds . Nonetheless, the socio-economic status of EGAT employees is probably better than some of the most severely economically disadvantaged Thais, and the study did not include the severely ill or disabled subjects excluded from employment. The prevalence of subjects with decreased eGFR and the CKD risk factors in our study are comparable to other a representative cross-sectional population surveys from Thailand . Therefore, our risk score should applicable in assessing the risks of developing decreased GFR in community-based Thai subjects, although caution may be necessary in extrapolating findings to groups not represented in our study (for example younger or very old subjects or those who are institutionalized).
We showed that clinical or combined clinical and laboratory risk models based on simple parameters available in the primary care had good accuracy and discrimination power to estimate the 10-year probability of developing Decreased GFR in a middle age to elderly non-institutionalized Thai general population. In addition, we validated these prediction models in a separate group of Thai subjects from a similar employment background. Additional studies are necessary to determine the validity of these scores in other Asian populations. The benefits of the two risk scores derived from these models in increasing self-awareness for CKD risks and for targeting individuals in the Thai or other Asian communities at high risk of CKD for specific interventions requires further studies.
Area under the curve
Chronic Kidney Disease–Epidemiology Collaboration
Diastolic blood pressure
The Electric Generating Authority of Thailand
End-Stage renal disease
Glomerular filtration rate
Modification of diet in renal disease study.
Systolic blood pressure
Jha V, Garcia-Garcia G, Iseki K, Li Z, Naicker S, Plattner B, Saran R, Wang AY, Yang CW. Chronic kidney disease: global dimension and perspectives. Lancet. 2013;382(9888):260–72.
Levey AS, Coresh J. Chronic kidney disease. Lancet. 2012;379(9811):165–80.
Hsu CC, Hwang SJ, Wen CP, Chang HY, Chen T, Shiu RS, et al. High prevalence and low awareness of CKD in Taiwan: a study on the relationship between serum creatinine and awareness from a nationally representative survey. Am J Kidney Dis. 2006;48:727–38.
Ingsathit A, Thakkinstian A, Chaiprasert A, Sangthawan P, Gojaseni P, Kiattisunthorn K, Ongaiyooth L, Vanavanan S, Sirivongs D, Thirakhupt P, et al. Prevalence and risk factors of chronic kidney disease in the Thai adult population: Thai SEEK study. Nephrol Dial Transplant. 2010;25(5):1567–75.
Taal MW, Brenner BM. Predicting initiation and progression of chronic kidney disease: Developing renal risk scores. Kidney Int. 2006;70(10):1694–705.
Wilson PW, D'Agostino RB, Levy D, Belanger AM, Silbershatz H, Kannel WB. Prediction of coronary heart disease using risk factor categories. Circulation. 1998;97(18):1837–47.
D'Agostino RB Sr, Grundy S, Sullivan LM, Wilson P, Group CHDRP. Validation of the Framingham coronary heart disease prediction scores: results of a multiple ethnic groups investigation. JAMA. 2001;286(2):180–7.
Stanifer JW, Muiru A, Jafar TH, Patel UD. Chronic kidney disease in low- and middle-income countries. Nephrol Dial Transplant. 2016;31(6):868–74.
Chang JM, Hwang SJ, Tsukamoto Y, Chen HC. Chronic kidney disease prevention--a challenge for Asian countries: report of the Third Asian Forum of Chronic Kidney Disease Initiatives. Clin Exp Nephrol. 2012;16(2):187–94.
Jha V, Prasad N. CKD and Infectious Diseases in Asia Pacific: Challenges and Opportunities. Am J Kidney Dis. 2016;68(1):148–60.
O'Seaghdha CM, Lyass A, Massaro JM, Meigs JB, Coresh J, D'Agostino RB Sr, Astor BC, Fox CS. A Risk Score for Chronic Kidney Disease in the General Population. Am J Med. 2012;125(3):270–7.
Kshirsagar AV, Bang H, Bomback AS, Vupputuri S, Shoham DA, Kern LM, Klemmer PJ, Mazumdar M, August PA. A simple algorithm to predict incident kidney disease. Arch Intern Med. 2008;168(22):2466–73.
Vathesatogkit P, Woodward M, Tanomsup S, Ratanachaiwong W, Vanavanan S, Yamwong S, Sritara P. Cohort Profile: The electricity generating authority of Thailand study. Int J Epidemiol. 2011;41(2):359–65.
Kidney Disease: Improving Global Outcomes (KDIGO) CKD Work Group. KDIGO clinical practice guideline for the evaluation and management of chronic kidney disease. Kidney Int Suppl. 2013;3:1–150.
National High Blood Pressure Education P. The Seventh Report of the Joint National Committee on Prevention, Detection, Evaluation, and Treatment of High Blood Pressure. Bethesda (MD): National Heart, Lung, and Blood Institute (US); 2004.
Organization WH. Definition and diagnosis of diabetes mellitus and intermediate hyperglycemia: report of a WHO/IDF consultation; 2006. p. 29.
Sritara P, Cheepudomwit S, Chapman N, Woodward M, Kositchaiwat C, Tunlayadechanont S, Sura T, Hengprasith B, Tanphaichitr V, Lochaya S, et al. Twelve-year changes in vascular risk factors and their associations with mortality in a cohort of 3499 Thais: the Electricity Generating Authority of Thailand Study. Int J Epidemiol. 2003;32(3):461–8.
Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem. 1972;18(6):499–502.
Levey AS, Coresh J, Greene T. Expressing the Modification of Diet in Renal Disease Study equation for estimating glomerular filtration rate with standardized serum creatinine values. Clin Chem. 2007;53:766–72.
Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF 3rd, Feldman HI, Kusek JW, Eggers P, Van Lente F, Greene T, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–12.
Tsai WC, Wu HY, Peng YS, Ko MJ, Wu MS, Hung KY, Wu KD, Chu TS, Chien KL. Risk Factors for Development and Progression of Chronic Kidney Disease: A Systematic Review and Exploratory Meta-Analysis. Medicine (Baltimore). 2016;95(11):e3013.
Sullivan LM, Massaro JM, D'Agostino RB Sr. Presentation of multivariate data for clinical use: The Framingham Study risk score functions. Stat Med. 2004;23(10):1631–60.
Hosmer DWHT, Le Cessie S, Lemeshow S. A comparison of goodness-of-fit tests for the logistic regression model. Stat Med. 1997;16(9):965–80.
Zou KH, O’Malley AJ, Mauri L. Receiver-Operating Characteristic Analysis for Evaluating Diagnostic Tests and Predictive Models. Circulation. 2007;115(5):654–7.
DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. 1988;44(3):837–45.
Stevens LA, Claybon MA, Schmid CH, Chen J, Horio M, Imai E, Nelson RG, Van Deventer M, Wang HY, Zuo L, 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–62.
Harrell FE Jr, Lee KL, Mark DB. Multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat Med. 1996;15(4):361–87.
Chien KL, Lin HJ, Lee BC, Hsu HC, Lee YT, Chen MF. A prediction model for the risk of incident chronic kidney disease. Am J Med. 2010;123(9):836–46. e832
Matsushita K, Selvin E, Bash LD, Astor BC, Coresh J. Risk implications of the new CKD Epidemiology Collaboration (CKD-EPI) equation compared with the MDRD Study equation for estimated GFR: the Atherosclerosis Risk in Communities (ARIC) Study. Am J Kidney Dis. 2010;55(4):648–59.
Delanaye P, Cavalier E, Mariat C, Krzesinski JM, Rule AD. Estimating glomerular filtration rate in Asian subjects: where do we stand? Kidney Int. 2011;80(5):439–40.
Rule AD, Teo BW. GFR estimation in Japan and China: what accounts for the difference? Am J Kidney Dis. 2009;53(6):932–5.
Yamwong S, Kitiyakara C, Vathesatogkit P, Saranburut K, Chittamma A, Cheepudomwit S, Vanavanan S, Akrawichien T, Sritara P. Effects of glomerular filtration rate estimating equations derived from different reference methods on staging and long term mortality risks of chronic kidney disease in a Southeast Asian cohort. Nephrology (Carlton). 2016;21(8):678–86.
Praditpornsilpa K, Lekhyananda S, Premasathian N, Kingwatanakul P, Lumpaopong A, Gojaseni P, et al. Prevalence Trend of Renal Replacement Therapy in Thailand: Impact of Health Economics Policy. Med Assoc Thai. 2011;94(Suppl. 4):S1–6.
Levin A, Tonelli M, Bonventre J, Coresh J, Donner JA, Fogo AB, Fox CS, Gansevoort RT, Heerspink HJL, Jardine M, et al. Global kidney health 2017 and beyond: a roadmap for closing gaps in care, research, and policy. Lancet. 2017. doi:10.1016/S0140-6736(17)30788-2. [Epub ahead of print]
Ene-Iordache B, Perico N, Bikbov B, Carminati S, Remuzzi A, Perna A, Islam N, Bravo RF, Aleckovic-Halilovic M, Zou H, et al. Chronic kidney disease and cardiovascular risk in six regions of the world (ISN-KDDC): a cross-sectional study. Lancet Glob Health. 2016;4(5):e307–19.
Wen CP, Matsushita K, Coresh J, Iseki K, Islam M, Katz R, McClellan W, Peralta CA, Wang H, de Zeeuw D, et al. Relative risks of chronic kidney disease for mortality and end-stage renal disease across races are similar. Kidney Int. 2014;86(4):819–27.
Collaboration NCDRF. Trends in adult body-mass index in 200 countries from 1975 to 2014: a pooled analysis of 1698 population-based measurement studies with 19.2 million participants. Lancet. 2016;387(10026):1377–96.
We would like to thank the staff at EGAT and Ramathibodi for clinical and administrative assistance.
The EGAT study is sponsored by the Faculty of Medicine, Ramathibodi Hospital, Mahidol University, the National Research Council, the Thailand Research Fund, the Thai Heart Association, the Thai Health Foundation, and the Higher Education Research Promotion and National Research University Development, Office of the Higher Education Commission.
Availability of data and materials
The datasets used and analysed during the current study available from the corresponding author on reasonable request.
Ethics approval and consent to participate
This study was conducted in accordance with the Helsinki Declaration and approved by the Ethics Committee, Ramathibodi Hospital, Mahidol University, Thailand. Written informed consent was obtained.
Consent for publication
All authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Baseline characteristics of the participants with serum creatinine at both baseline (2002-2003) and follow-up (2012-2013) compared to all subjects with serum creatinine at the baseline visit. Table S2. Alternative clinical model (Model 1) with body mass index. Table S3. Clinical Model (Model 1) with proteinuria. Table S4. Baseline characteristics of EGAT 3 participants in the Validation dataset. (DOCX 20 kb)