Abstract
Background and objectives Estimates of the effect of estimated GFR (eGFR) decline on mortality have focused on populations with normal kidney function, or have included limited information on factors previously shown to influence the risk of death among patients with CKD.
Design, setting, participants, & measurements We retrospectively assessed the effect of rate of eGFR decline on survival of patients with CKD receiving primary care through a large integrated health care system in central Pennsylvania between January 1, 2004, and December 31, 2009.
Results A total of 15,465 patients were followed for a median of 3.4 years. Median rates of eGFR change by those in the lower, middle, and upper tertiles of eGFR slope were −4.8, −0.6, and 3.5 ml/min per 1.73 m2/yr, respectively. In Cox proportional hazard modeling for time to death, adjusted for baseline proteinuria, changes in nutritional parameters, and episodes of acute kidney injury during follow-up (among other covariates), the hazard ratio for those in the lower (declining) and upper (increasing) eGFR tertiles (relative to the middle, or stable, tertile) was 1.84 and 1.42, respectively. Longitudinal changes in nutritional status as well as episodes of acute kidney injury attenuated the risk only modestly. These findings were consistent across subgroups.
Conclusions eGFR change over time adds prognostic information to traditional mortality risk predictors among patients with CKD. The utility of incorporating eGFR trends into patient-risk assessment should be further investigated.
Introduction
Non–dialysis-dependent chronic kidney disease (NDD-CKD), affecting up to 15% of the adult US population, increases the risk of death compared with those without CKD, predominantly because of an increased burden of cardiovascular disease (1–3). Traditionally, fixed estimates of GFR (eGFR) at a single point in time have been used to define this risk (4). Several factors—including proteinuria, malnutrition, inflammation, and episodes of acute kidney injury (AKI)—have been proposed as mediators of the cardiovascular morbidity and mortality observed in this population (5–7).
Recently, several retrospective analyses have identified the rate of eGFR decline as an independent risk factor for all-cause and cardiovascular mortality and morbidity (8–11). These studies focused largely on older populations, and/or predominantly included patients with normal, or near normal, kidney function. None of the reported risk models included information about proteinuria or AKI during follow-up, and only one included longitudinal markers of nutritional status—all important mortality predictors in the NDD-CKD population (9).
Using a longitudinal, single health system electronic database with detailed demographic, pharmaceutical, clinical, laboratory, claims, and outcomes information, we investigated the effect of change in eGFR on mortality among a large, primary care population with NDD-CKD in central Pennsylvania between 2004 and 2009. We hypothesized that among those with a higher rate of eGFR decline, when compared with those with more stable longitudinal renal function, mortality risk is increased.
Materials and Methods
This retrospective cohort study was reviewed and approved under “exempt” status by the Geisinger Medical Center IRB in February 2010. The data source was EpicCare, Geisinger Medical Center's electronic health record (12).
Study Population
Patients eligible for the study included the following: those between the ages of 18 and 88 years at the time of first qualifying eGFR value; enrollment for primary care at any Geisinger facility at any time during the period February 1, 2004, through March 31, 2009; baseline NDD-CKD (defined by two outpatient GFR values estimated using the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation between 15 and 59 ml/min per 1.73 m2 on at least 2 dates separated by at least 90 days but no more than 365 days—the CKD-defining window period); and at least one additional eGFR value a minimum of 6 months after the first qualifying value (13). The study index date was the date of the first qualifying eGFR value, once all other eGFR criteria were met. Follow-up for study outcomes occurred through December 31, 2009.
Patients were excluded if they suffered hospital-associated AKI during the CKD-defining window period (defined by an increase of 50% or more in serum creatinine during a hospitalization for any cause, with the baseline serum creatinine defined as the lowest recorded value between 90 days before the index admission and the hospital discharge date) (7). Community-associated AKI (defined as a 50% increase in serum creatinine occurring in the outpatient setting and within 30 days of another outpatient serum creatinine result) during the CKD-defining window period was likewise an excluding event. Additional exclusions included any prior history of solid-organ transplantation (ICD-9-CM diagnostic code V42, V42.0, V42.6, V42.7, V42.8, V42.83, V42.84, V42.89, V42.9, V43, 996.8, 996.80 through 996.87, or 996.89); prior hemo- or peritoneal dialysis (ICD-9-CM diagnostic code 585.5, 585.6, V45.11, V45.12, or V56.31); a diagnosis of metastatic cancer (ICD-9-CM diagnostic code 191.0 or 191.1); or an active prescription for cytotoxic or immunosuppressive therapy (cyclophosphomide, chlorambucil, prednisone, solumedrol, rituximab, tacrolimus, cyclosporine, or mycophenolate mofetil). Patients were censored for ESRD, defined as an eGFR <10 ml/min per 1.73 m2, renal transplantation, or initiation of maintenance hemo- or peritoneal dialysis (ICD-9-CM diagnostic code 585.5, 585.6, V42.0, 996.81, V45.11, V45.12, or V56.31), or for any 18-month period of time without a serum creatinine result. Thus, qualifying patients were followed from the date of first qualifying eGFR result until ESRD, death, the date of the last recorded serum creatinine value if no additional lab testing was performed during the subsequent 18 months, or end of study (December 31, 2009).
Study Covariates
Demographic and clinical information, along with medication prescription history and laboratory and vital sign data, were extracted from the database. The Charlson Comorbidity Index score was calculated for each patient (14). A diagnosis was coded if present by ICD-9-CM code on the medical history on at least two distinct outpatient encounter dates during the 12 months before the index date. Medication prescriptions active at the time of the index date were coded. The laboratory values and vital sign data closest in time and before (but within 12 months of) the index date were used for all analyses. For proteinuria and hemoglobin A1C, because more than 10% of the cohort did not have a result recorded for these parameters, we created additional dummy variables to indicate whether or not the test was ordered, and included these dummy variables in the univariate and adjusted models, as appropriate (15). We included changes in serum albumin and body mass index (BMI) during study follow-up as covariates in survival modeling (16). In both cases, a linear regression model was fit for each patient to estimate annual change in the nutritional parameter, using a smoothed estimate based on the change from baseline for each additional measurement of the parameter. This estimate was included in the Cox survival model along with the baseline value. Hospital- and community-associated AKI (defined as above) during follow-up were also included as study covariates in a time-dependent fashion. Detailed definitions of study covariates are available as Supplemental Information.
Serum creatinine was measured using the isotope dilution/mass spectroscopy (IDMS)–traceable Roche enzymatic method throughout the entirety of the study period (17). Instrument calibration at Geisinger laboratories is performed according to manufacturer's specifications. Proteinuria was coded as present if any semiquantitative dipstick analysis result was 30 mg/dl or greater, if assessment by urinary protein-to-creatinine ratio was greater than 0.2 mg protein/g creatinine, or if a 24-hour urine collection included at least 300 mg of protein.
Exposure
All eGFR values were obtained using outpatient serum creatinine values and the CKD-EPI estimating formula (13). The slope of the eGFR versus time curve was calculated using all available outpatient values. This value was then used to stratify patients by tertile of rate of eGFR change in univariate and Cox proportional hazards models. This stratification scheme resulted in three groups roughly corresponding to declining (lower tertile), stable (middle tertile), and increasing (upper tertile) eGFR trends over time; the terms “declining,” “stable,” and “increasing” are used throughout the remainder of this report to represent the three tertiles. eGFR change was also analyzed as a continuous variable in a cubic spline model examining the relationship between GFR variability and death.
Outcome
The primary study outcome was death. Information on vital status for Geisinger primary care recipients is updated monthly by query of the Social Security Administration's data set through the National Technical Information Service (18).
Statistical Analyses
Data are presented as mean and SD for continuous variables and as frequency and percentage for categorical variables. Baseline comparisons across the three eGFR groups were made using the Kruskal–Wallis nonparametric and Pearson chi-squared tests, as appropriate. A Poisson regression model was used to estimate the mortality rate across groups and expressed as incident rate ratios.
Time until death was analyzed using the Cox proportional hazard regression model. Variables significant at P < 0.10 in the unadjusted analysis, or previously associated with death among patients with NDD-CKD, were included in the model. Estimated changes in BMI and serum albumin by linear regression were used in the Cox model to account for changes in nutritional status during follow-up. Both hospital- and community-associated AKI were treated as time-dependent variables by sequentially accounting for such events as they occurred before death in survival models. Subgroup analyses were performed to examine the stability of the observed association between eGFR change and death across CKD subpopulations. Interactions also were tested to examine for differential effects of eGFR change across subpopulations in the full Cox proportional hazards regression model.
eGFR change was also analyzed as a continuous variable in unadjusted and fully adjusted models to estimate the relationship between GFR variability and death. A restricted cubic spline was created to parameterize this relationship as it allows for detecting nonlinear effects of exposure-outcome relationship.
Results
A total of 15,465 patients with NDD-CKD were enrolled in the study (Figure 1). The median (interquartile range [IQR]) change in eGFR in the decreasing, stable, and increasing groups was −4.8 (−8.2 to −3.2), −0.6 (−1.4 to 0.0), and 3.5 (1.9 to 6.7) ml/min per 1.73 m2/yr, respectively. “Cut” points for tertiles were −2.2 and 0.8 ml/min per 1.73 m2. In this study population, 41.3% had an eGFR that rose over time. The median (IQR) follow-up time was 2.9 (1.8 to 4.7), 4.5 (2.9 to 5.5), and 3.0 (1.9 to 4.5) years in those with decreasing, stable, and increasing eGFR groups over time, respectively. The median (IQR) number of serum creatinine values used to calculate eGFR values was 10 (6 to 17), 11 (7 to 17), and 8 (5 to 14) in the decreasing, stable, and increasing groups, respectively (P < 0.001 for the comparison of decreasing and stable eGFR, and increasing and stable eGFR).
Cohort development. CA-AKI, community-associated acute kidney injury; eGFR, estimated GFR; HA-AKI, hospital-associated acute kidney injury; NDD-CKD, non–dialysis-dependent chronic kidney disease.
The baseline characteristics of the study population are described in Table 1. Relative to the stable group, those in the decreasing eGFR group were more likely to be men, have a smoking history, and have diabetes mellitus, heart failure, and other comorbid disease conditions (as reflected in higher Charlson Comorbidity Index scores). Interestingly, they were less likely to have hypertension or hyperlipidemia. Proteinuria was more prevalent in decliners than in those with stable kidney function.
Baseline characteristics of adult patients with NDD-CKD, by tertiles of eGFR change
During follow-up, AKI was common in all three groups; rates of hospital-associated AKI were 76.2, 34.8, and 60.6 per 1000 patient-years in the declining, stable, and increasing eGFR groups, respectively, whereas rates of community-associated AKI were 38.6, 13.7, and 23.6 per 1000 patient-years, respectively. Nutritional status, as reflected by longitudinal changes in BMI and serum albumin levels, varied over time. Across all three groups, both BMI and serum albumin levels fell during follow-up (Table 2).
Change in nutritional parameters (BMI, serum albumin) over time among patients with NDD-CKD, by tertile of eGFR changea
There were 560, 379, and 388 deaths (with corresponding mortality rates of 34.7, 17.5, and 24.0 deaths per 1000 patient-years, P < 0.001 for all intergroup comparisons) in the declining, stable, and increasing eGFR groups, respectively.
In univariate analyses, baseline covariates associated with death included the following: older age, male gender, and smoking history; a history of heart failure, peripheral vascular disease, dementia, chronic liver disease, or higher Charlson Comorbidity Index score; proteinuria; and a prescription for a beta blocker, loop diuretic, aldosterone antagonist, coumadin, aspirin, calcium acetate, and insulin. Both hospital- and community-associated AKI were likewise associated with a higher likelihood of death. White race, hypertension, hyperlipidemia, increasing BMI over time, and the use of nonsteroidal anti-inflammatory agents at baseline were associated with a lower risk of death. Both declining and increasing eGFR over time were associated with a higher risk of death, relative to those in the stable eGFR group.
In unadjusted and multivariate adjusted Cox proportional hazards models examining time to death, both decreasing and increasing eGFR were independently associated with death, relative to those in the stable group (Table 3). Other factors independently associated with death include older age, male gender, chronic liver disease, a higher Charlson Comorbidity Index score, use of beta blockers, aspirin, or coumadin, and baseline proteinuria. Higher baseline BMI, diastolic BP, eGFR, and HDL cholesterol levels all were associated with a lower risk of death. With the sequential addition of markers of changing nutritional status over time, as well as the inclusion of AKI events (both hospital- and community-associated) in a time-dependent fashion, the association between eGFR change and death was attenuated but qualitatively unchanged. Interactions between beta blockers and heart failure, as well as angiotensin-converting enzyme inhibitors and/or angiotensin receptor blocker use and proteinuria, were not significant. Vital status ascertainment was complete on patients censored for 18-month serum creatinine gaps; the risk of death associated with both decreasing (1.88, 95% confidence interval [CI] 1.68 to 2.10) and increasing (1.64, 95% CI 1.46 to 1.84) eGFR change over time relative to those in the stable eGFR group was not qualitatively different when these censored patients were included and analyzed in the fully adjusted Cox proportional hazard model.
Cox proportional hazards for time to death among patients with non–dialysis-dependent CKD and changing kidney function over time
In fully adjusted subgroup analyses (Figure 2), the point estimate for the risk of death associated with decreasing eGFR over time was relatively consistent across groups. For those with increasing eGFR, the adjusted risk was likewise consistent across subgroups, with the exception of those younger than 60 years of age. Interactions were tested for each subpopulation of interest. Among those with declining (and relative to those with stable) eGFR, the mortality risk for those without diabetes was significantly larger than the risk for those with diabetes (P = 0.02). Among those with increasing (and relative to those with stable) eGFR, the mortality risk for those with proteinuria was significantly larger than the risk for those without proteinuria (P = 0.003). Interactions for all other subgroups with either declining or increasing eGFR were nonsignificant.
Risk of death among subgroups of patients with baseline nondialysis dependent CKD, by tertile of eGFR change over time. For the BMI/ALB subgroup analysis, patients with both BMI and serum albumin results (or with either one alone if the other was not assessed in follow-up) at the time of last follow-up equal to or greater than the baseline value were classified as nutritionally stable; those with one or both values less than baseline value were classified as “malnourished.” 439 (2.8%) patients had no follow-up BMI or serum albumin; these patients were not included in the analysis of nutritional change. ALB, serum albumin; BMI, body mass index; decl, declining GFR group; incr, increasing GFR group; neg, negative; pos, positive.
Figure 3 details the association between eGFR change as a continuous variable and the risk of death.
Functional point estimate (and 95% pointwise confidence band) between estimated GFR change and the hazard ratio for death among patients with non–dialysis-dependent chronic kidney disease, estimated by a restricted cubic spline function. Adjusted model includes age, gender, race, smoking history, hypertension, heart failure, peripheral vascular disease, hypertension, Charlson Comorbidity Index score, and dementia present at the time of cohort entry; a prescription for beta blocker, loop diuretic, aldosterone antagonist, calcium acetate, insulin, coumadin, or aspirin at time of cohort entry; systolic and diastolic BP, proteinuria, HDL and LDL cholesterol levels, estimated GFR at time of cohort entry, hospital- and community-associated acute kidney injury during follow-up, change in body mass index during follow-up, and baseline serum albumin. Bold lines represent point estimates; faint lines represent 95% confidence bands.
Discussion
In this study, declining eGFR over time was independently associated with an increased risk of death compared with those with more stable renal function. This risk persisted despite adjusting for factors previously known to be associated with death among patients with CKD, and was consistent across population subgroups.
Rifkin et al. evaluated eGFR decline among older community-dwelling adults using a retrospective analysis of 4380 participants with predominantly normal kidney function in the Cardiovascular Health Study (9). Rapid decline was associated with a 50% to 80% increased risk of death, depending on the method used to estimate GFR (19,20). This study included changes in weight over time as a model covariate, but did not include other measures of nutritional status, assessment of proteinuria, or AKI inj in risk modeling.
Unexpectedly, we found that increasing eGFR over time was independently associated with death. Al-Aly et al., reporting on the association between kidney function decline and mortality among a group (n = 4171) of military veterans with early stage 3 CKD and rheumatoid arthritis, found that approximately one third of their cohort had rising eGFR values over time, an observation which is inconsistent with consensus statements about progression of CKD (11,21). In addition to the observation that those with >4 ml/min per year of eGFR decline had approximately a 50% increased risk of death at a median follow-up of nearly 6 years, relative to those with stable kidney function, the investigators found that for those with a rising eGFR over time, there was a nonsignificant trend toward an increased risk of death relative to those with normal, age-related decline in kidney function. The investigators speculate that this observation is explained by malnutrition among those with increasing eGFR over time.
The reasons for our finding that increasing eGFR over time associates with an increased risk of death are not immediately apparent. Loss of lean body mass would be potentially explanatory, as poor nutritional status has been previously identified as a predictor of death among patients with CKD (6). We accounted for nutritional status at baseline (weight, BMI, and serum albumin levels) and during follow-up (changes in BMI and serum albumin); inclusion of these nutritional markers in our risk models only slightly attenuated the mortality risk associated with decline or gain in eGFR. Furthermore, in the subgroup of patients with declining nutritional parameters during follow-up, the adjusted mortality risk associated with an increasing eGFR over time remained robust, at 50% relative to those with stable function. Resolving AKI, while a recognized mortality risk factor in this population, might also account for a rising eGFR over time; however, inclusion of AKI events during follow-up did not qualitatively affect our study findings (and patients with both community- and hospital-associated AKI at entry were excluded). Improving intrinsic renal function is theoretically explanatory, yet is physiologically unlikely to occur in more than 40% of a general geriatric population with a substantial burden of cardiovascular disease and established CKD, and should not exacerbate mortality risk.
Changes in intravascular volume status, due to conditions like heart failure, liver disease, or other hemodynamic insults, might also explain an increase in eGFR over time, due to a dilutional effect on serum creatinine values. Subclinical physiologic changes in intravascular volume status might manifest as a rising eGFR over time, and in turn influence mortality risk among patients with CKD. Fluctuations in serum creatinine (and thus eGFR) due to subtle intravascular volume shifts is particularly common among those with chronic tubulointerstitial disease from small-vessel renal vascular disease. Finally, the statistical phenomenon of mean regression might theoretically explain some degree of eGFR rise over time among those with reduced baseline values. It is important to recognize that this study was not designed to assess the mortality risk associated with improving eGFR over time; thus, the observations reported in this regard remain hypothesis-generating and require further investigation.
Study limitations include the possible persistence of confounders that may have influenced the observed outcomes. These would include, among others, duration of CKD before study entry. The censoring of patients without an eGFR measurement after 18 months may have biased our sample toward a sicker, more resource-intense population with inherently higher risk of death; however, mortality risk did not qualitatively change after inclusion of these censored patients in the sensitivity analysis.
We did not use a gold standard measure of GFR, but measured GFR is not practical for typical clinical use or large epidemiologic studies. Our model assumes a constant rate of eGFR change over time and does not account for the influence of time-dependent rate changes. Given that eGFR slope calculations included data points that qualify as community-associated AKI, the influence of such events relative to the non–AKI associated decline or increase in eGFR over time is difficult to ascertain definitively in this study, despite the inclusion of all types of AKI in survival modeling. Finally, the largely Caucasian population receiving care in this hospital system precludes extrapolation of study findings to non-Caucasian populations.
The results of this study would suggest that identification of eGFR decline over time may enhance risk prognostication for patients with established CKD. Further research testing interventions to address this risk are warranted, as are investigations targeting a better understanding of the relationship between rising eGFR and outcomes among patients with CKD.
Disclosures
Dr. Perkins reports having received honoraria and research support from American Regent and research support from Amgen.
Acknowledgments
We thank Ms. Amanda Bengier, Mr. Raymond Menapace, and Mr. Joseph Leader for assistance with data extraction and programming, Ms. Haiyan Sun for assistance with analysis, and Dr. S. Mark Poler for early critical feedback on methods and design. In addition, we thank Dr. Linda Fried for her helpful critique of a draft of the manuscript. The study was supported by the Administrative Committee for Research at Geisinger Medical Center.
This study was presented in abstract form at the annual meeting of the American Society of Nephrology; November 16 through 21, 2010; Denver, CO.
Footnotes
Published online ahead of print. Publication date available at www.cjasn.org.
See related editorial, “Change in Kidney Function over Time and Risk for Adverse Outcomes: Is an Increasing Estimated GFR Harmful?,” on pages 1805–1806.
- Received January 18, 2011.
- Accepted April 25, 2011.
- Copyright © 2011 by the American Society of Nephrology