- Research article
- Open Access
- Open Peer Review
Variability in estimated glomerular filtration rate values is a risk factor in chronic kidney disease progression among patients with diabetes
BMC Nephrologyvolume 16, Article number: 34 (2015)
It is unknown whether variability of estimated Glomerular Filtration Rate (eGFR) is a risk factor for dialysis or death in patients with chronic kidney disease (CKD). This study aimed to evaluate variability of estimated Glomerular Filtration Rate (eGFR) as a risk factor for dialysis or death to facilitate optimum care among high risk patients.
A longitudinal retrospective cohort study of 70,598 Veterans Health Administration veteran patients with diabetes and CKD (stage 3–4) in 2000 with up to 5 years of follow-up. VHA and Medicare files were linked to derive study variables. We used Cox proportional hazards models to evaluate association between time to initial dialysis/death and key independent variables: time-varying eGFR variability (measured by standard deviation (SD)) and eGFR means and slopes while adjusting for prior hospitalizations, and comorbidities.
There were 76.7% older than 65 years, 97.5% men, and 81.9% Whites. Patients were largely in early stage 3 (61.2%), followed by late stage 3 (28.9%), and stage 4 (9.9%); 29.1%, 46.8%, and 73.3%, respectively, died or had dialysis during the follow-up. eGFR SDs (median: 5.8, 5.1, and 4.0 ml/min/1.73 m2 ) and means (median: 54.1, 41.0, 27.2 ml/min/1.73 m2) from all two-year moving intervals decreased as CKD advanced; eGFR variability (relative to the mean) increased when CKD progressed (median coefficient of variation: 10.9, 12.8, and 15.4). Cox regressions revealed that one unit increase in a patient’s standard deviation of eGFRs from prior two years was significantly associated with about 7% increase in risk of dialysis/death in the current year, similarly in all three CKD stages. This was after adjusting for concurrent means and slopes of eGFRs, demographics, prior hospitalization, and comorbidities. For example, the hazard of dialysis/death increased by 7.2% (hazard ratio:1.072; 95% CI = 1.067, 1.080) in early stage 3.
eGFR variability was independently associated with elevated risk of dialysis/death even after controlling for eGFR means and slopes.
Chronic kidney disease (CKD) affects 15% of the U.S. adult population and is associated with a high burden of morbidities and mortality [1,2]. Patients with CKD are a heterogeneous group characterized by varying rates of progression to a devastating outcome such as end-stage renal disease (ESRD), likely due to presence/absence of various risk factors [3-6]. Identifying risk factors of CKD progression or death will directly benefit patients by helping clinicians optimize health care and provide proper referral to multi-disciplinary clinics based on patients’ risk profiles, hence may help reduce patients’ risk of morbidities and mortality. Existing literature has provided support for many risk factors of CKD progression or death. A greater degree of estimated glomerular filtration rate (eGFR) downward trend was found to be associated with higher mortality . Other identified predictors of CKD progression included demographics (age, sex, race), laboratory test results (proteinuria, phosphorus, albumin, hemoglobin), and comorbidities (diabetes, vascular disease, and hypertension) [3,7,8]. Recently, an algorithm was proposed to combine those factors in predicting a patient’s progression . However, CKD progression is usually not constant over time [10-12]. For example, acute kidney injury—sudden deterioration of kidney function, is a critical adverse event in CKD progression [13-15].
eGFR is often used in tracking CKD progression. Although it is well recognized by clinicians that a single eGFR may not provide sufficient and complete information to one’s kidney function. For instance, two patients with the same eGFR value may have different resilience to kidney offense. Namely, other risk factors being equal, an adverse event or an acute disease (i.e., a stimulus) such as an infection may deteriorate kidney function (less renal resilience) in one patient but not in the other (greater resilience). In this study, we have proposed to use variability of eGFR values to measure kidney resilience, with larger eGFR variability indicating reduced kidney resilience to stimuli, and considered the state of kidney resilience as a potential risk factor of CKD progression. There are studies that have attempted to use repeatedly measured eGFRs to model CKD progression using different statistical approaches, including modeling of non-linear trajectories [12,16-18]. A recent study measured variability of kidney function using repeated measured eGFRs  and found it to be associated with increased risk of death in stage 3–5 CKD patients.
In this study, we also aimed to evaluate eGFR variability as a risk factor; however, we have used different study designs and methods, and included not only death but also dialysis as the adverse CKD outcomes of interest. Our hypothesis is that greater eGFR variability is associated with a higher risk of incurring dialysis or dialysis-free death, independently of the CKD progression trend. In essence, we believe although renal function decline may be inevitable for patients with CKD, those with similar decline rate of eGFR may be further differentiated by degree of renal resilience, as reflected by eGFR variability.
We chose to focus on CKD patients with coexisting diabetes in the study. Patients with diabetes are at higher risk of developing CKD. The risk of morbidity and mortality increase in CKD patients with diabetes , highlighting the importance of assessing and managing risk factors associated with progression of CKD. Furthermore, both conditions have considerable costs and negative impacts to individuals, health care systems, as well as the society as a whole, and they are increasingly prevalent. In the Veteran Health Administration (VHA) health care system, about one in five veterans have diabetes  and 30% of veterans with diabetes have CKD [22,23].
Study design, data sources and study population
This was a longitudinal retrospective cohort study of VHA veteran patients with diabetes in fiscal year (FY) 2000 (9/30/1999-10/1/2000) from a national research database containing linked longitudinal patient records of inpatient and outpatient services . We determine the baseline CKD status using two outpatient eGFR values less than 60 ml/min/1.73 m2: the first (index) eGFR value was identified in FY2000, followed by a subsequent qualifying eGFR 90 to 365 days apart from the index eGFR.  eGFR was calculated with the four-variable version of the Modification of Diet in Renal Disease equation [24,25]. The baseline period was the 12 months prior to (and including) the date of the qualifying eGFR. Individuals were followed through FY2004, with follow-up time ranging between less than a year and close to five years.
Of the 653,064 VHA clinic users with diabetes in FY2000, we excluded individuals according to different criteria detailed in Figure 1. Patients enrolled in Medicare Health Maintenance Organization (HMO) were excluded because detailed Medicare HMO data is not routinely collected by Center for Medicare and Medicaid Services. We also excluded individuals with evidence of prior dialysis, using International Classification of Diseases, Ninth Edition, Clinical Modification (ICD-9-CM) diagnosis and procedure codes, Current Procedural Terminology codes, revenue codes (Medicare Part A), and VHA clinic stop codes (; Additional file 1). We determined CKD stage using the qualifying eGFR based on the National Kidney Foundation Kidney Disease Outcome Quality Initiative (K/DOQI) criteria , and divided stage 3 into early (eGFR 45–59 ml/min/1.73 m2) and late stage (eGFR 30–44 ml/min/1.73 m2) [1,23]. We excluded individuals in stage 5 CKD, who are close to developing ESRD.
The VHA New Jersey Healthcare System Institutional Review Board approved the study. This study received a waiver of consent from the facility Internal Review Board.
The outcome was time from the date of the qualifying eGFR to the event of dialysis or dialysis-free death, whichever came first. The event time was censored if the individual was enrolled in Medicare HMO or alive at the end of FY 2004. Death information came from the VHA Vital Status File . We ascertained dialysis using the method described above, but selected codes explicitly referring to dialysis treatment (see Additional file 1).
The primary exposure variable was variability in the outpatient eGFR values. For each patient, we calculated standard deviations (SDs), concurrent means, and slopes of eGFR progression (derived from a linear regression) of repeatedly measured eGFRs in a series of (five) two-year moving intervals. In our analyses, the number of actual eGFR measurements per patient was at least two in each of the five intervals. These variables were later used as time-varying variables in Cox regression models. Other independent variables included demographic variables (age, sex, and race/ethnicity), VHA priority status, and clinical factors: prior hospitalization and medical comorbidities, all derived from the baseline period. The VHA categorizes veteran enrollees into 1 of 8 mutually exclusive priority groups  based on poverty level and service connected disability. We regrouped patients into four levels (severely disabled, moderately disabled, poverty, and co-payment required). We applied the Selim Comorbidity Index  to determine health status using ICD-9-CM codes, including indicators of physical (cardiovascular and non-cardiovascular conditions separately), and mental comorbidities. Cardiovascular conditions included angina pectoris, myocardial infarction, congestive heart failure (CHF), peripheral vascular diseases, and stroke. Mental comorbidities included schizophrenia, depression, bipolar disorder, anxiety, post-traumatic stress disorder, and alcohol abuse.
We fitted the Cox proportional hazards models to evaluate the association between eGFR variability and risk of dialysis or dialysis-free death by each CKD stage. Prior to statistical modeling, we evaluated possible correlation/association among independent variables and found little indication of co-linearity among variables within each model.
Primary analysis: We divided the follow-up time into yearly intervals and included the eGFR SDs, means and slopes (derived from the two-year interval prior to each yearly interval in the follow-up) as lagged time-dependent covariates. Because higher values (hence greater means) are usually associated with greater SDs, and lower values (smaller means) are associated with smaller SDs, we included both SDs and means (derived from the same time intervals) as time-dependent covariates in the Cox models. We also included interactions of eGFR SDs and slopes to assess if the influence of eGFR variability differs with changes in the slopes.
Sensitivity analysis: We conducted four sensitivity analyses to evaluate robustness of the findings from the primary analysis. First, regarding the length of the time interval used to derive the time-varying eGFR SDs, means and slopes, we re-fitted Cox regression models using re-calculated time-varying eGFR SDs, means and slopes from a one-year interval (prior to each yearly interval in the follow-up). Second, we re-divided the follow-up time into quarterly intervals and included eGFR SDs, means and slopes from prior eight quarters (two years) as time-dependent variables in Cox regression models. Third, we examined the effect of extreme slopes on our main findings by removing those with absolute values greater than 36.5 or 25.0 ml/min/1.73 m2 per year; this also simultaneously removed extreme SDs. Fourth, as an alternative, for each patient in each defined time interval we obtained the standard deviation of eGFR residuals from the linear regression model used to obtained a eGFR slope , and refitted Cox regression models. Results of these sensitivity analyses were very similar and we chose to present only findings of the main analysis.
Of the 70,598 study patients, 61.2%, 28.9%, and 9.9% were in early stage 3, late stage 3, and stage 4 CKD, respectively. Table 1 also shows that they were largely men (97.5%), and the mean age was 70(±9) years. There were 81.9% white, 14.2% African American, 2.3% Hispanics, and 1.6% were of other racial groups. At the baseline, 36.6% had prior hospitalizations, 69.5% had cardiovascular conditions, 98.3% had other non-cardiovascular conditions, and 19.0% had mental comorbidities. Individuals in stage 4 were more likely to be non-White, and had higher rates of prior hospitalizations and comorbidities than those in stage 3. During follow-up, 19,939 (28.3%) had dialysis-free death, and 7,193 (10.3%) had dialysis. The rate of dialysis or dialysis-free death was higher for patients with more severe CKD: 29.1% for early stage 3, 46.8% for late stage 3, and 73.3% for stage 4 CKD. Overall, the median follow-up time was 3.8 years (range: <0.01-4.8).
Generally, eGFR SDs and means decreased as CKD advanced (Table 2). To account for the effect of mean in interpreting SDs as a measure of variability, we have supplemented the statistic of coefficient of variation (defined as standard deviation divided by mean). The summary statistics show that eGFR variability (relative to the mean) increased when CKD progressed. The median coefficient of variation was 10.9 for early stage 3, 12.8 for late stage 3, and 15.4 for stage 4 CKD. The median eGFR SDs from all of the two-year intervals were 5.8, 5.1, and 4.0 ml/min/1.73 m2 for early stage 3, late stage 3, and stage 4, respectively. A downward trend of eGFRs was observed for all CKD stages, with the magnitude increasing as CKD progressed. The median slopes went from −2.2 in early stage 3 to −3.5 ml/min/1.73 m2/year in stage 4.
Figure 2 delineates the relationship between variability in eGFRs (SDs) and presence/absence of the events of dialysis or dialysis-free death. Across all three CKD stages and for all follow-up years, individuals with events always had significantly (p < 0.001 from t-tests) greater average SDs than those without events, with the exception of the last four follow-up years in stage 4 (p > 0.05).
In Table 3, results from Cox regression models show that a greater eGFR SD was associated with elevated risk of dialysis or dialysis-free death, with or without adjusting for other independent variables. For example, in early stage 3, 1 ml/min/1.73 m2 increase in eGFR SD in prior two years would increase the risk of dialysis or dialysis-free death in the current year by 7.6% after adjusting for the eGFR means and slopes (Model 1, adjusted hazard ratio (AHR) = 1.076; 95% CI = (1.071, 1.080)), or by 7.2% after additionally adjusting for other demographic and clinical variables (Model 2, AHR = 1.072; 95% CI = (1.067, 1.080)). Similar associations were found in the other CKD stages. Interactions of eGFR SDs and slopes were not statistically significant.
Other independent variables were also significantly related to the risk of dialysis or dialysis-free death. If the eGFR declined by 10 ml/min/1.73 m2 per year in previous two years, the risk would increase by 1.0% (Model 2, AHR = 1.010; 95% CI = (1.003, 1.017)) in early stage 3 CKD, and 6.7% (Model 2, AHR = 1.067; 95% CI = (1.038, 1.096) in stage 4 CKD. Compared to patients younger than 45 years old, patients 65 or older in early stage 3, 75 or older in late stage 3, and 85 or older in CKD 4 were at greater risk. Women tended to have lower risk (AHR range: 0.63 to 0.66) than men. African Americans, but not other racial/ethnic minorities were at higher risk (1.07 to 1.14) than Whites. Individuals with severe service connected disability also had higher risk (1.17 to 1.26) compared to others. As regards clinical factors, individuals with hospitalization history (1.28 to 1.56), presence of cardiovascular conditions (1.21 to 1.52) and mental comorbidities (early stage 3: 1.12; late stage 3: 1.10) at the baseline had elevated risks. The magnitude of the AHR for these clinical factors somehow decreased when CKD advanced.
In this large cohort study of longitudinal outpatient eGFR values, we found that greater variability in eGFRs, measured by SDs in a dynamic/time-varying fashion over the follow-up, was associated with elevated risk of dialysis or dialysis-free death, even after controlling for downward trends (slopes) in CKD progression and other demographic and clinical factors. In all three CKD stages, one ml/min/1.73 m2 increase in a patient’s SD of eGFRs from prior two years was significantly associated with about 7% increase in risk of dialysis/death in current year (e.g., adjusted hazard ratio in early stage 3:1.072; 95% CI = (1.067, 1.080)). A recent study  of stage 3-5 CKD patients also found the variability in repeatedly measured eGFRs to be associated with death, although eGFR variability was measured differently and only death was evaluated as the outcome in the study. Both studies used VHA patients and have detailed individual-level electronic health records in the VHA health care system. However, our study had Medicare information for totality of data rather than VHA data alone (most VHA patients are also Medicare eligible). Additionally, we simultaneously evaluated eGFR slopes and eGFR variability as potential risk factors by including both as time-varying independent variables in the statistical models. We found that a faster downward eGFR progression was associated with increased risk of CKD outcomes. This is consistent with the literature  and the common knowledge among clinicians caring for CKD patients.
Since our findings support the hypothesis of eGFR variability (measured by standard deviation) being associated with increased risk of dialysis or dialysis free death, the key message for clinicians or health policy makers is to consider using a simple statistic—standard deviation, to assist in clinical care of CKD patients. We have proposed to measure eGFR variability using the statistic of standard deviation because it reflects the essence of the concept and is familiar among clinicians and policy makers. The derivation of this statistic can be easily built into any existing systems of electronic health records and the result (value) can be instantly made available. Because this statistic is a summary measure of the variability (non-stability) of renal function, it can be used as a tool by clinicians in assessing CKD management for a patient and assist in clinical decision making for CKD care. For example, a clinician can compare a patient’s current statistic with earlier results or the norm (to be established as future work) and an elevated value (indicating increased variability of renal function) may alert the clinician to an acute kidney injury or other reasons for such increased variability. This might direct clinicians to look into the patient’s history of eGFR values and health conditions, and to engage in active dialogues with the patient regarding the patient’s recent changes of health conditions and/or medications that may affect kidney function. Standard deviation decreases as mean decreases (also as seen in Table 2). Since the renal function of patients with CKD is expected to decline over time (hence decreased eGFR values), one would expect their standard deviations of eGFR values decrease as well. Therefore, it is important to note that in comparing a patient’s past and current eGFR SDs, an increased SD would not be expected and may warrant more detailed assessment.
Although variability in eGFR in CKD progression may have not been well studied, the concept of variability has been applied to different physiological markers in medicine and is supported by various biological mechanisms. For example, within-visit variability of blood pressure was associated with stroke and all-cause mortality . Heart rate variability predicted worst outcomes in perioperative and intensive care settings . Furthermore, home blood pressure and heart rate variability were associated with cardiovascular events , and hemoglobin (Hb) variability was associated with higher mortality among patients with ESRD  or type 2 diabetes . In patients with type 1 diabetes, HbA1C variability was predictive of development of microalbuminuria, progression of kidney disease and cardiovascular events .
In this study, we propose to measure eGFR variability as a proxy for kidney resilience to stimuli. Greater variability of eGFR values over time means decreased renal function and more stable eGFR values over time means less damaging renal function. We caution not to confuse this with the renal reserve (RR) test. The renal reserve measures ability of a kidney to increase basal GFR after a protein overload ; a larger value means a better kidney ability/function. In CKD patients, renal reserve is usually preserved but with a lower magnitude, suggesting that remaining nephrons are able to maintain some response to protein overload despite already being in a state of hyperfiltration [35,36].
Understanding predictors and pathways to ESRD/dialysis in patients with CKD is important. A recent study reported that patients were highly concerned with the uncertainty of CKD progression and that physicians were frustrated with their inability to predict the course of CKD . Uncertainty in predicting CKD progression may limit clinicians’ ability in effective CKD management and impair preparation for dialysis such as vascular access placement for high-risk patients . Understanding CKD progression and identifying patients at high risk of adverse outcomes is also valuable for health care resource and program planning to improve CKD care. The uncertainty rises mainly due to acute changes, eGFR variability, and/or other characteristics of the course of CKD progression, already stated above. Pathways other than linear slopes have been described, but require a rather difficult application of statistical modeling . Prospective clinical studies may help to address these issues. We suggest that regional health maintenance organizations with electronic health records attempt to replicate our findings. We also recognize that the general unavailability of electronic medical records makes implementation of eGFR variability problematic. However, the evaluation of individual level longitudinal data using electronic records are consistent with Federal policies to establish meaningful use criteria; it provides the opportunity to simultaneously improve quality, safety, and efficiency and the results can provide decision support for national high-priority conditions.
We would like to note that there are existing studies addressing the question of changes in repeatedly measured eGFR values. It was reported that patients in the quartile with the greatest eGFR change were more likely to incur death [39-41] and ESRD in one study  although not in another . In these studies, eGFR change was defined as annual decline (measured as change in eGFR values from two visits relative to the time elapsed between these visits)  or percentage of change in eGFR values (measured as change in eGFR values from two visits relative to the first eGFR value)  or both , or mean absolute eGFR residual values [41,43]; patients then were assigned into quartiles based on the value of eGFR change. The findings of these studies were largely consistent with ours, although different approaches were used.
There are a few unique features of this study: 1) the advanced modeling of variability, slopes, and means of repeated eGFR values over time as time-dependent variables (calculated in each two-year interval with sensitivity analyses) allows us to effectively test the hypothesis and address the research question; 2) The explicit adjustment of eGFR slopes (as a time-dependent covariate) allows us to evaluate whether there is an interaction between eGFR variability and CKD progression or only existence of their main effects; 3) inclusion of both dialysis and death as adverse CKD outcomes for evaluation provides us overall understanding of risks of two critical adverse outcomes of CKD progression in relation to eGFR variability; and 4) evaluation of the associations by CKD stage permits us to compare whether findings depend on baseline renal status.
Our study has limitations. First, the findings of the study may not be generalized beyond populations of largely male, elderly and with diabetes; and CKD patients identified based on two eGFR values less than 60 ml/min/1.73 m2. Second, we were not able to identify transient or temporary dialysis, resulting in censoring some patients whose kidney function was resilient and their creatinine values later improved during the follow-up. Lastly, the primary goal of this study was to evaluate eGFR variability as a new risk factor in CKD progression in order to enhance clinical decisions by clinicians, and resources planning by policy makers/administrators; nonetheless, assessing factors affecting eGFR variability as future work is warranted. We suspect variability in eGFRs may be related to medication changes, acute diseases (infections, cardiac events, etc.) and hospitalizations (although we included only outpatient creatinine values in this study). Indeed, in an effort to assess eGFR variability as a risk factor without the influence of some severe cardiac conditions, we conducted two separate additional analyses where we excluded 1) patients with a prior history of CHF at baseline or 2) patients with CHF anytime during the study period (but prior to CKD outcomes if applicable). The findings remained similar.
In conclusion, this study provides support for evaluating eGFR variability as an independent risk factor of dialysis/death, regardless of CKD stage. It has important clinical implication because patients may benefit from monitoring of eGFR variability in addition to their rates of eGFR decline.
Chronic kidney disease
Estimated glomerular filtration rate
Veterans health administration
Health maintenance organization
End-stage renal disease
Congestive heart failure
Adjusted hazard ratio
Go AS, Chertow GM, Fan D, McCulloch CE, Hsu C. Chronic kidney disease and the risks of death, cardiovascular events, and hospitalization. N Engl J Med. 2004;351(13):1296–305.
U.S. Renal Data System. USRDS 2011 Annual Data Report: Atlas of Chronic Kidney Disease and End-Stage Renal Disease in the United States. Bethesda, MD: National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases; 2011.
Levin A, Djurdjev O, Beaulieu M, Er L. Variability and risk factors for kidney disease progression and death following attainment of stage 4 CKD in a referred cohort. Am J Kidney Dis. 2008;52(4):661–71.
Onuigbo MA. The natural history of chronic kidney disease revisited–a 72-month Mayo Health System Hypertension Clinic practice-based research network prospective report on end-stage renal disease and death rates in 100 high-risk chronic kidney disease patients: a call for circumspection. Adv Perit Dial. 2009;25:85–8.
Zhang AH, Tam P, LeBlanc D, Zhong H, Chan CT, Bargman JM, et al. Natural history of CKD stage 4 and 5 patients following referral to renal management clinic. Int Urol Nephrol. 2009;41(4):977–82.
Barrett BJ, Garg AX, Goeree R, Levin A, Molzahn A, Rigatto C, et al. A nurse-coordinated model of care versus usual care for stage 3/4 chronic kidney disease in the community: a randomized controlled trial. Clin J Am Soc Nephrol. 2011;6(6):1241–7.
Al-Aly Z, Zeringue A, Fu J, Rauchman MI, McDonald JR, El-Achkar TM, et al. Rate of kidney function decline associates with mortality. J Am Soc Nephrol. 2010;21(11):1961–9.
Hoefield RA, Kalra PA, Baker P, Lane B, New JP, O’Donoghue DJ, et al. Factors associated with kidney disease progression and mortality in a referred CKD population. Am J Kidney Dis. 2010;56(6):1072–81.
Tangri N, Stevens LA, Griffith J, Tighiouart H, Djurdjev O, Naimark D, et al. A predictive model for progression of chronic kidney disease to kidney failure. JAMA. 2011;305(15):1553–9.
Ambrogi V, Thilly N, Boini S, Virion JM, Kessler M, Briancon S, et al. Patterns and predictors of kidney function decline in the last year prior to dialysis. Nephron Clin Pract. 2009;111(12):c95–101.
Heaf JG, Mortensen LS. Uraemia progression in chronic kidney disease stages 3–5 is not constant. Nephron Clin Pract. 2011;118(4):c367–74.
Li L, Astor BC, Lewis J, Appel LJ, Lipkowitz MS, Toto RD, et al. Longitudinal Progression Trajectory of GFR Among Patients With CKD. Am J Kidney Dis. 2012;59(4):504–12.
Ali T, Khan I, Simpson W, Prescott G, Townend J, Smith W, et al. Incidence and outcomes in acute kidney injury: a comprehensive population-based study. J Am Soc Nephrol. 2007;18(4):1292–8.
Ishani A, Xue JL, Himmelfarb J, Prescott G, Townend J, Smith W, et al. Acute kidney injury increases risk of ESRD among elderly. J Am Soc Nephrol. 2009;20(1):223–8.
Lafrance JP, Djurdjev O, Levin A. Incidence and outcomes of acute kidney injury in a referred chronic kidney disease cohort. Nephrol Dial Transplant. 2010;25(7):2203–9.
Klahr S, Levey AS, Beck GJ, Caggiula AW, Hunsicker L, Kusek JW, et al. The effects of dietary protein restriction and blood-pressure control on the progression of chronic renal disease. Modification of Diet in Renal Disease Study Group. N Engl J Med. 1994;330(13):877–84.
Wright J, Barkis G, Greene T, Agodoa LY, Appel LJ, Charleston J, et al. African American Study of Kidney Disease and Hypertension Study Group. Effect of blood pressure lowering and antihypertensive drug class on progression of hypertensive kidney disease: results from the AASK trial. JAMA. 2002;288(19):2421–31.
Gassman JJ, Greene T, Wright Jr JT, Agodoa L, Bakris G, Beck GJ, et al. Design and statistical aspects of the African American Study of Kidney Disease and Hypertension (AASK). J Am Soc Nephrol. 2003;14 Suppl 2:S154–65.
Al-Aly Z, Balasubramanian S, McDonald JR, Scherrer JF, O’Hare AM. Greater variability in kidney function is associated with an increased risk of death. Kidney Int. 2012;82(11):1208–2014.
Moen MF, Zhan M, Hsu VD, Walker LD, Einhorn LM, Seliger SL, et al. Frequency of hypoglycemia and its significance in chronic kidney disease. Clin J Am Soc Nephrol. 2009;4(6):1121–7.
Miller DR, Safford MM, Pogach LM. Who has diabetes? Best estimates of diabetes prevalence in the Department of Veterans Affairs based on computerized patient data. Diabetes Care. 2004;27 Suppl 2:B10–21.
Kern EFO, Maney M, Miller DR, Tseng CL, Tiwari A, Rajan M, et al. Failure of ICD-9-CM codes to identify patients with comorbid chronic kidney disease in diabetes. Health Serv Res. 2006;41(2):564–80.
Tseng CL, Kern EFO, Miller D, Tiwari A, Maney M, Rajan M, et al. Survival Benefit of Nephrologic Care in Patients with Diabetes and Chronic Kidney Disease. Arch Intern Med. 2008;168(1):55–62.
National Kidney Foundation. K/DOQI clinical practice guidelines for chronic kidney disease: evaluation, classification, and stratification. Am J Kidney Dis. 2002;39(2 Suppl 1):S1–266.
Levey AS, Coresh J, Greene T, Stevens LA, Zhang YL, Hendriksen S, et al. Chronic Kidney Disease Epidemiology Collaboration. Using standardized serum creatinine values in the modification of diet in renal disease study equation for estimating glomerular filtration rate. Ann Intern Med. 2006;145(4):247–54.
U.S. Department of Veterans Affairs, VA Information Resource Center (VIReC). VHA Vital Status File. (http://www.hsrd.research.va.gov/for_researchers/cyber_seminars/archives/694-notes.pdf). (Accessed January 20, 2015).
U.S. Department of Veterans Affairs. Health Benefits: priority group table. (http://www.va.gov/healthbenefits/resources/priority_groups.asp). (Accessed January 20, 2015).
Selim AJ, Fincke G, Ren XS, Lee A, Rogers WH, Miller DR, et al. Comorbidity assessments based on patient report: results from the Veterans Health Study. J Ambul Care Manage. 2004;27(3):281–95.
Yang W, Israni RK, Brunelli SM. Hemoglobin variability and mortality in ESRD. J Am Soc Nephrol. 2007;18(12):3164–70.
Muntner P, Shimbo D, Tonelli M, Reynolds K, Arnett DK, Oparil S. The relationship between visit-to-visit variability in systolic blood pressure and all-cause mortality in the general population: findings from NHANES III, 1988 to 1994. Hypertension. 2011;57(2):160–6.
Mazzeo AT, La Monaca E, Di Leo R, Vita G, Santamaria LB. Heart rate variability: a diagnostic and prognostic tool in anesthesia and intensive care. Acta Anaesthesiol Scand. 2011;55(7):797–811.
Johansson JK, Niiranen TJ, Puukka PJ, Jula AM. Prognostic value of the variability in home-measured blood pressure and heart rate: the Finn-Home Study. Hypertension. 2012;59(2):212–8.
Takao T, Matsuyama Y, Yanagisawa H, Kikuchi M, Kawazu S. Association between HbA1c variability and mortality in patients with type 2 diabetes. J Diabetes Complications. 2014;28(4):494–9.
Wadén J, Forsblom C, Thorn LM, Gordin D, Saraheimo M, Groop PH. Finnish Diabetic Nephropathy Study. A1C variability predicts incident cardiovascular events, microalbuminuria, and overt diabetic nephropathy in patients with type 1 diabetes. Diabetes. 2009;58(11):2649–55.
Musso CG. Renal reserve test: its methodology and significance. Saudi J Kidney Dis Transpl. 2011;22(5):990–3.
Barai S, Gambhir S, Prasad N, Sharma RK, Ora M. Functional renal reserve capacity in different stages of chronic kidney disease. Nephrology (Carlton). 2010;15(3):350–3.
Schell JO, Patel UD, Steinhauser KE, Ammarell N, Tulsky JA. Discussions of the kidney disease trajectory by elderly patients and nephrologists: a qualitative study. Am J Kidney Dis. 2012;59(4):495–503.
Levin A, Djurdjev O. On being better kidney doctors: understanding trajectories, probabilities, predictability, and people. Am J Kidney Dis. 2012;59(4):475–7.
Matsushita K, Selvin E, Bash LD, Franceschini N, Astor BC, Coresh J. Change in estimated GFR associates with coronary heart disease and mortality. J Am Soc Nephrol. 2009;20:2617–24.
Turin TC, Coresh J, Tonelli M, Stevens PE, de Jong PE, Farmer CK, et al. Change in the estimated glomerular filtration rate over time and risk of all-cause mortality. Kidney Int. 2013;83(4):684–91.
Perkins RM, Tang X, Bengier AC, Kirchner HL, Bucaloiu ID. Variability in estimated glomerular filtration rate is an independent risk factor for death among patients with stage 3 chronic kidney disease. Kidney Int. 2012;82(12):1332–8.
Turin TC, Coresh J, Tonelli M, Stevens PE, de Jong PE, Farmer CK, et al. Short-term change in kidney function and risk of end-stage renal disease. Nephrol Dial Transplant. 2012;27(10):3835–43.
Perkins RM, Kirchner HL, Hartle JE, Bucaloiu ID. Estimated glomerular filtration rate variability and risk of end-stage renal disease among patients with Stage 3 chronic kidney disease. Clin Nephrol. 2013;80(4):256–62.
This research was supported by the United States of America Department of Veteran Affairs (VA) Medical Service Epidemiology Grant (Drs. Pogach & Miller), the VA Diabetes Quality enhancement Research Initiative Coordinating Center (Dr. Pogach), the VA Health Services Research and Development (HSR&D) grant IIR 04–205 (Dr. Pogach), and the VA HSR&D grant IIR 11–077 (Dr. Tseng).
The authors declare that they have no competing interests.
CT had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis; CT contributed in the design, conduct, and completion of the study, as well as interpretation of the results and writing of the manuscript. JL, SL, LP, and DM contributed in the design of the study, interpretation of results, and writing of the manuscript. OS contributed in data derivation and analysis, interpretation of results, and writing of the manuscript. MM contributed in data acquisition, interpretation of results, and writing of the manuscript. All authors read and approved the manuscript.
The views expressed in this article are those of the authors and do not necessarily reflect the position or policy of the Department of Veterans Affairs or the United States government.
Codes used in the study for determination of dialysis and access (fistula procedure). Description: This file lists codes used in the study for determination of dialysis in the baseline period (to remove study subjects) and in the outcome years (for study of dialysis as an outcome).