|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Epidemiology and Outcomes |





* Department of Global Epidemiology, Amgen, Inc., Thousand Oaks, California;
Department of Internal Medicine, Veterans Affairs Medical Center,
Dialysis Outcomes and Practice Patterns Study (DOPPS), University Renal Research and Education Association and
Department of Biostatistics, University of Michigan, Ann Arbor, Michigan
Address correspondence to: Dr. Brian D. Bradbury, 1 Amgen Center Drive, MS: 24-1-C, Thousand Oaks, CA 91320; Phone: 805-313-4343; Fax: 805-447-1984; E-mail: bradbury{at}amgen.com
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Studies that quantify mortality rates and identify mortality predictors among incident HD populations are useful for understanding the influence of predialysis care and length of time on dialysis. Data from the 2004 Annual Dialysis Report (1) indicate higher mortality rates for patients who have received HD for >5 yr as compared with <2 yr, suggesting that length of time on HD modifies mortality risk. Most patients begin HD with several comorbid conditions that may worsen or develop additional comorbidities with continued dialysis. We hypothesized that the influence on mortality of specific comorbid conditions may change over time, whereas the influence of others remains constant. This study identified a period of elevated mortality risk among incident HD patients and examined the magnitude of associations between various patient characteristics and mortality during this period of elevated risk and throughout the first year on HD using a representative sample of US incident HD patients from the Dialysis Outcomes and Practice Patterns Study (DOPPS).
| Materials and Methods |
|---|
|
|
|---|
Study Population
Information was collected at the time of enrollment on 15,500 and 13,000 dialysis patients in DOPPS I and II, respectively. For this analysis, only US patients who began dialysis <30 d before enrollment and had at least 1 d of follow-up were included (n = 4802). Limiting the study to incident patients reduced the amount of missing information on predialysis care and permitted analyses of mortality risk immediately after dialysis initiation. Demographic characteristics, medical history, predialysis medical care, and laboratory data that were collected at study enrollment were used for this analysis. Patients were followed from enrollment until the first of the following: Death, transplant, transfer out of the facility, switch to peritoneal or home dialysis, withdrawal from dialysis, or end of the study period. Patients who withdrew from HD were followed for an additional 60 d to capture patients who chose to stop dialysis and died within 2 mo of their last HD.
Outcomes
The primary outcomes in this analysis were all-cause and cause-specific mortality. Causes of death were grouped as follows: Cardiovascular; vascular; infection; a combined category of liver disease, gastrointestinal, and other known causes; and a final category of unknown cause. In addition, deaths that occurred subsequent to withdrawal, irrespective of the immediate cause of death, were examined. Deaths during the first year of follow-up were categorized as deaths that occurred within versus after the identified early period of elevated mortality risk. In analyses of deaths that occurred during this early period, patients who survived were censored at the end of the early period. Analyses of deaths that occurred beyond this early period of elevated risk included only patients who survived beyond this period; those who survived beyond 365 d were censored at 366 d.
Covariates
Patient information included four types: Demographics and dialysis care, laboratory values, medical history, and predialysis medical care. Demographic characteristics included age (18 to 44, 45 to 64, 65 to 74, 75+ yr), gender, race (white, nonwhite), body mass index ([kg/m2]; <20, 20 to <25, 25 to <30, 30+), and primary cause of ESRD (diabetes, hypertension, other [includes glomerulonephritis; secondary glomerulonephritis-vasculitis; interstitial nephritis; polycystic kidney disease, congenital; and renal tumor] and missing); dialysis care was measured by vascular access type in place at dialysis initiation (arteriovenous fistula, graft [bovine, synthetic], catheter [cuffed percutaneous, temporary], other, missing). Laboratory values were categorized according to the Kidney Disease Outcomes Quality Initiative (K/DOQI) guidelines (2729) and include the following: Albumin (<3.5,
3.5 g/dl, unknown), albumin-corrected calcium (<8.4, 8.4 to 9.5, >9.5 g/dl, unknown), hemoglobin (<11.0, 11.0 to <12.0, 12.0+ g/dl, unknown), phosphorus (<3.5, 3.5 to 5.5, >5.5 g/dl, unknown), and calcium-phosphorus product (<55, >55 mg2/ml2, unknown). The effect of parathyroid hormone (PTH) was not assessed because of large amount (>50%) of missing information. Composite variables that reflect patient medical history characteristics were described previously (25); each variable (dichotomized as yes/no) was captured at enrollment. Visit to a nephrologist >1 mo before enrollment in DOPPS (pre-ESRD nephrology care; yes/no) was used to assess pre-ESRD medical care.
Statistical Analyses
All analyses were conducted on the combined DOPPS I and DOPPS II samples unless stated otherwise. Descriptive statistics characterized the study population. Life-table analyses were used to examine the estimated hazard function, by month, for patients who were captured in the census (all HD patients in the facility) and separately for patients who were enrolled in the DOPPS I and II samples to identify a period of elevated risk. Discrete survival analysis using logistic regression with intervals of 1 mo was used to test for differences in mortality rates for early versus late periods (e.g., 0 to 90 versus 91 to 365, 0 to 120 versus 121 to 365, 0 to 150 versus 151 to 365). The t test and
2 analyses were used to compare characteristics of patients who were identified in the census but not included in the DOPPS II sample (nonparticipants) with those who were included in the DOPPS II sample (participants). This subanalysis was limited to DOPPS II patients because date of first dialysis was not collected in the DOPPS I census, and, therefore, incident patients, other than those in the DOPPS I sample, could not be identified. Cox proportional hazards regression was used to generate crude and adjusted hazard ratios (AHR) and 95% confidence intervals (CI) for the association between baseline patient characteristics and mortality in different periods during the first year of follow-up. The full models adjusted for all patient characteristics. Two-sided P < 0.05 was considered statistically significant. Time-dependent Cox regression also was used to examine the interaction between patient characteristics and time at risk for mortality, controlling for all other patient characteristics. In particular, tests for differences in the HR for a particular covariate before and after the period of elevated mortality were performed. All covariates initially were assessed for time-dependent effects using plots of ln(ln[S(t)]) versus time for each covariate in turn, with S(t) estimated from models that were adjusted for other baseline predictors but stratified on the covariate of interest. The final model included the significant time-dependent interaction terms that were identified in previous models (P < 0.05) along with all other main effects. Analyses were conducted using SAS 9.1 (SAS Institute, Cary, NC).
| Results |
|---|
|
|
|---|
Among the 4802 patients in the combined DOPPS I and DOPPS II sample, 4292 (89%) enrolled in DOPPS within 14 d of dialysis initiation; 1977 (41%) enrolled at their first dialysis session. The median time between first dialysis and DOPPS enrollment was 3 d. The distribution of patient characteristics for the population overall and according to mortality during the first year is presented in Table 1. Overall, 51% of the patients were 65 yr or older, 56% were male, and 65% were white.
|
|
|
HR estimates for the association between each predictor and death up to 120 d, 121 to 365 d, and during the first 365 d overall, adjusted for other patient characteristics, are provided in Table 2. Patient characteristics that were significantly associated with an increased risk for mortality during the first 120 d included age >75 yr (AHR 2.49), white race (AHR 1.40), catheter vascular access (AHR 1.71), albumin levels <3.5 g/dl (AHR 1.57), phosphorus levels <3.5 mg/dl (AHR 1.47), HIV/AIDS (AHR 2.85), history of cancer (AHR 1.41), history of lung disease (AHR 1.33), history of neurologic disease (AHR 1.50), and history of a psychiatric disorder (AHR 1.35). Factors that were associated with a significantly lower mortality risk during the first 120 d included the diagnosis of hypertension (AHR 0.55) and pre-ESRD nephrology care >1 mo before initiation of dialysis (AHR 0.65). Cause-specific death rates (per 100 person-years) were lower in the first 120 d among those who saw a nephrologist >1 mo before ESRD compared with those who did not (cardiac 9.24 versus 17.7; infection 2.2 versus 3.5; vascular 1.8 versus 3.2; liver/gastrointestinal/other 3.2 versus 6.3; unknown 5.3 versus 7.9 [data not shown]). Similarly, deaths subsequent to withdrawal were lower among those who saw a nephrologist >1 mo before ESRD compared with those who did not (4.5 versus 8.2 per 100 person-years, respectively).
|
Finally, we evaluated whether the association between each predictor and mortality differed between the <120-d and 121- to 365-d periods using predictor by time interaction terms while simultaneously adjusting for the other main effects. Of variables that were significantly associated with death in either the <120-d or 121- to 365-d time period, only serum calcium <8.4 mg/L (test for interaction with time, P = 0.009), HIV/AIDS status (P = 0.02), and pre-ESRD nephrologist care (P = 0.003) were observed to have a statistically significant different impact on mortality risk in the first 120-d than in the subsequent 121- to 365-d period.
| Discussion |
|---|
|
|
|---|
Few studies have assessed mortality rates and mortality predictors in incident HD patients (2,20,21). Evaluations that use large population-based databases such as Medicare are difficult because of the 90-d entitlement period (1), which precludes complete ascertainment of patients with ESRD during this period. Analyses of mortality during this initial period are prone to bias, particularly given the elevated mortality risk during the first 120 d that we observed, unless conducted in populations in which all patients who initiate dialysis can be identified and are guaranteed continued coverage. Studies that were conducted in nonpopulation-based samples of incident patients have reported 90-d mortality rates ranging from 12 (2022) to 16% (23). A Michigan study of nearly 2400 patients with ESRD reported deaths within 90 d of initiation of dialysis in 10% of patients (24). Most recently, Soucie and McClellan (2) reported 6% of HD patients in an ESRD registry in the southeastern United States died within 90 d after dialysis initiation. In a subanalysis (data not shown), we found that 5.6% of deaths occurred within 90 d of HD initiation; 8% of patients died within the first 120 d, representing 46% of all deaths during the first year of follow-up. Both our study and that of Soucie and McClellan (2) were population based, minimizing concerns about selection bias, and both include relatively recent data from patients who initiated HD. The latter point may account for the lower mortality rates as compared with findings from earlier studies, potentially reflecting improvements in pre-ESRD treatment, including earlier referral to a nephrologist, which we and others (30,31) suggest improves survival once on dialysis. It is interesting that patients with diabetes were more likely than patients without diabetes (78.6 versus 72.6%; P < 0.0001) to receive predialysis care from a nephrologist, but this greater access to care did not seem to explain the unexpected lower 120-d mortality rate that was observed among patients with diabetes.
Our findings with respect to factors that are associated with elevated mortality risk are consistent with previous studies among prevalent HD patient populations, which are composed of patients with a range of times since dialysis initiation. The strength of older age and white race as mortality risk factors throughout the first 365 d found in our study among incident patients and in other studies of prevalent populations (10,11) suggests a consistent effect of these factors on mortality risk. We also found that factors such as low serum albumin and phosphorus levels, catheter vascular access, and absence of referral to a nephrologist >1 mo before HD initiation were associated strongly with early mortality. These results support the importance of the ongoing recommendations for arteriovenous fistula placement in the predialysis period and avoidance of temporary or permanent catheter placement (27). Furthermore, these results suggest that greater medical attention before ESRD onset may lead to improved survival (20,31). Our findings are consistent with those from the RightStart Program (32), which demonstrated significantly lower mortality rates for patients who were new to dialysis when specific disease management protocols were used before ESRD. Moreover, our results remained largely unchanged even after adjustment for pre-ESRD erythropoietin (EPO) use, which reflects one component of pre-ESRD care. Previous studies have reported a reduction in mortality for patients who were treated with EPO before initiating HD (33), but in this analysis, when pre-ESRD nephrologist care and EPO use were included in the same model, only pre-ESRD nephrologist care remained strongly protective in the first 120 d. In the 121- to 365-d period, when both pre-ESRD nephrologist care and EPO use were included in the same model, the effect of pre-ESRD nephrologist care was completely attenuated.
We found that specific chronic comorbid conditions, including history of CHF, lung and neurologic disease (e.g., dementia, depression, peripheral neuropathy), and HIV/AIDS, were associated with increased mortality risk throughout the first year of follow-up, although the magnitude of their effects was not constant. With the exception of CHF, the strength of these mortality risk factors increased with longer follow-up. History of CHF was more strongly associated with death within the first 120 d, emphasizing the prognostic importance of cardiovascular disease in patients who initiate HD (34).
In our analyses, the crude mortality rates for both the first 120 d and the subsequent 121 to 365 d were substantially lower when patients who withdrew from dialysis were censored at the time of their last HD (were not included as deaths in the estimation of incidence rates), reflecting the greater numbers of deaths captured during the 60 d after withdrawal. Although the majority of withdrawals occurred within the first 120 d and a larger proportion of deaths subsequent to withdrawal occurred during this period, when patients who withdrew were censored at the time of withdrawal, the resulting effect estimates remained largely unchanged (data not shown), indicating that mortality risk factors were not simply predictors of withdrawal. However, these findings do provide further evidence that previous estimates of mortality risk among prevalent patient populations may be underestimates for incident patients, particularly given the incomplete ascertainment of HD patients within the first 90 d of dialysis initiation.
This study should be evaluated in light of the following limitations. First, the risk for mortality among the entire HD patient population (DOPPS II census) is slightly higher than among study participants, especially in the first 120 d after HD initiation, suggesting possible selection bias. In fact, the death rates in the first 120 d (versus 121 to 365 d) among patients who were randomly selected for inclusion in the sample were 28.9 versus 21.0 per 100 person-years, respectively, compared with the analogous death rates among those who actually were included in the sample (25.4 versus 19.8 deaths per 100 person-years, respectively). Although the 121- to 365-d death rates are similar, the early death rates are quite different, showing an underascertainment of deaths during the early period, consistent with the comparison of the DOPPS II census versus the sample in Figure 1. One explanation for this underascertainment is that in many cases, hospital records for patients who died were archived, which made accessing these records difficult at the time of study form completion. In addition, the Health Insurance Portability and Accountability Act, which was enacted in April 2003 during the DOPPS II period, required patient consent for any new DOPPS II patients before enrollment. Consequently, patients who were close to death at the time of HD initiation may not have been approached by a study coordinator or agreed to participate, which could have resulted in a slightly healthier population being enrolled. Therefore, the 120-d mortality rate that we report may be a slight underestimation of the true mortality rate in the incident patient population.
Second, the time-varying influence of mortality risk factors that was observed in this study may reflect depletion of patients with specific comorbidities over time. Because we assessed only the influence of patient characteristics that were collected at baseline and included only patients who were new to dialysis, our effect estimates may underestimate the association between mortality and characteristics that often change after beginning dialysis (e.g., levels of hemoglobin, calcium-phosphorus product). As such, a single measurement that is taken at baseline might not characterize adequately a patients relevant exposure. The choice to analyze only baseline characteristics was influenced partially by the frequency of DOPPS data collection (approximately every 4 mo). Patients who died within 120 d of initiating HD would have no available follow-up information, but those who survived beyond 120 d potentially would have additional time points at which data could be collected. More important, however, the intent of this analysis was to identify patient characteristics at dialysis initiation that predict mortality over some early period of follow-up. Patient information that is collected subsequent to starting dialysis, particularly laboratory values, is sensitive to concomitant treatments and may be modified by increased contact with a nephrologist, obscuring associations with specific characteristics at HD initiation.
The strengths of this study also should be acknowledged. This analysis represents the largest population of incident HD patients examined to date. Because data collection for patients who were enrolled in DOPPS was not connected to Medicare reimbursement, underascertainment of the population at risk as a result of the 90-d coordination-of-benefits period was not a concern. In addition, because most patients were enrolled before patient consent was required, consent issues were minimized. Substantial information on patient medical history as well as laboratory parameters was collected during the baseline assessment, facilitating a more comprehensive investigation of mortality risk factors.
| Conclusion |
|---|
|
|
|---|
| Disclosures |
|---|
|
|
|---|
| Acknowledgments |
|---|
| Footnotes |
|---|
See the related editorial, "Predialysis Nephrology Care Improves Dialysis Outcome: Now What? Or Chapter Two," on pages 143145.
Received September 29, 2005. Accepted September 26, 2006.
| References |
|---|
|
|
|---|
Related Article
This article has been cited by other articles:
![]() |
M. K. Sigrist, A. Levin, L. Er, and C. W. McIntyre Elevated osteoprotegerin is associated with all-cause mortality in CKD stage 4 and 5 patients in addition to vascular calcification Nephrol. Dial. Transplant., June 2, 2009; (2009) gfp253v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
I-W. Wu, S.-Y. Wang, K.-H. Hsu, C.-C. Lee, C.-Y. Sun, C.-J. Tsai, and M.-S. Wu Multidisciplinary predialysis education decreases the incidence of dialysis and reduces mortality--a controlled cohort study based on the NKF/DOQI guidelines Nephrol. Dial. Transplant., June 2, 2009; (2009) gfp259v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Couchoud, M. Labeeuw, O. Moranne, V. Allot, V. Esnault, L. Frimat, B. Stengel, and for the French Renal Epidemiology and Information A clinical score to predict 6-month prognosis in elderly patients starting dialysis for end-stage renal disease Nephrol. Dial. Transplant., May 1, 2009; 24(5): 1553 - 1561. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. M. McClellan, H. Wasse, A. C. McClellan, A. Kipp, L. A. Waller, and M. V. Rocco Treatment Center and Geographic Variability in Pre-ESRD Care Associate with Increased Mortality J. Am. Soc. Nephrol., May 1, 2009; 20(5): 1078 - 1085. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. D. Bradbury, M. D. Danese, M. Gleeson, and C. W. Critchlow Effect of Epoetin Alfa Dose Changes on Hemoglobin and Mortality in Hemodialysis Patients with Hemoglobin Levels Persistently below 11 g/dL Clin. J. Am. Soc. Nephrol., March 1, 2009; 4(3): 630 - 637. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. B. Patwardhan, K. Kawamoto, D. Lobach, U. D. Patel, and D. B. Matchar Recommendations for a Clinical Decision Support for the Management of Individuals with Chronic Kidney Disease Clin. J. Am. Soc. Nephrol., February 1, 2009; 4(2): 273 - 283. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Carrera and M. Burnier Use of darbepoetin alfa in the treatment of anaemia of chronic kidney disease: clinical and pharmacoeconomic considerations NDT Plus, January 1, 2009; 2(suppl_1): i9 - i17. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |