Associations between the erythropoiesis-stimulating agent resistance index and the geriatric nutritional risk index of maintenance hemodialysis patients and increased mortality

Abstract

Background

There may be an association between the responsiveness to erythropoiesis-stimulating agent (ESA) and the nutritional status on mortality of maintenance hemodialysis (MHD) patients. Therefore, we assessed the associations of geriatric nutritional risk index (GNRI) and ESA resistance index (ERI) in mortality in MHD patients.

Methods

A total of 298 MHD patients were enrolled. Unadjusted and adjusted Cox regression analyses for all-cause mortality were performed. Odds ratios were estimated using a logistic regression model for the association between the GNRI and ERI and mortality, adjusting for potential confounders.

Results

During the follow-up period of 34.6 ± 6.1 months, 51 died with 31 cardiovascular deaths. Receiver operating curve analyses indicated the cutoff values of GNRI and ERI for mortality were 94.9 and 13.5, respectively. When subjects were stratified by ERI and GNRI values into four groups, those who had low GNRI and high ERI were associated with the highest risk of mortality among the four groups (log-rank p < 0.001).

Conclusions

High ERI and low GNRI were associated with an increased risk of all-cause mortality.

Background

Protein-energy wasting (PEW) is a major problem in end-stage renal disease (ESRD) patients, and its reported prevalence has ranged from about 20 to 78 % [1]. Since nutritional risk is one of the strongest predictors of increased mortality in maintenance hemodialysis (MHD) patients, and evaluation of nutritional status and nutritional management are essential clinical procedures for managing MHD patients [2, 3], it is necessary to regularly screen patients for PEW by making a variety of measurement with a variety of measures [4]. The geriatric nutritional risk index (GNRI) is a very simple and objective method of assessing nutritional status in a number of pathological conditions that is based on body weight, height, and serum albumin levels. Some studies have investigated the reliability of the GNRI as a tool for assessing malnutrition in HD patients [5] and predicting their overall mortality [6, 7].

The development of erythropoiesis-stimulating agents (ESAs) has resulted in substantial health benefits for ESRD patients, including improved quality of life, reduced need for blood transfusion, and increased exercise capacity. Unfortunately, however, a considerable proportion of ESRD patients exhibit a suboptimal hematologic response to ESAs as evidenced by the persistence of anemia despite adequate dosing or by the need for high ESA doses to achieve the recommended hemoglobin (Hb) target [8]. ESA doses and the Hb levels achieved in HD patients have been found to be associated with mortality [9], and hyporesponsiveness to ESA therapy is known to predict a poorer outcome [10]. Both high-dose ESA therapy and low Hb levels may contribute to poorer outcomes of HD patients who are hyporesponse to ESAs [11, 12].

The aim of the present study was to clarify the relationship between both the GNRI and responsiveness to ESA of MHD patients and their long-term mortality.

Methods

Study design

This was a prospective, observational cohort study conducted at a single center in Japan. The subjects were recruited from among patients who had been routinely dialyzed via an arteriovenous fistula in the dialysis unit of the Jyoban Hospital for at least 6 months. The Institutional Review Board of the Jyoban Hospital approved all study protocols (no. 11), and the protocols were carried out in accordance with the Declaration of Helsinki guidelines regarding ethical principles for medical research involving human subjects. Written informed consent was obtained from every subject.

HD patients with a malignancy, active inflammation, liver cirrhosis, gastrointestinal bleeding, cardiac valvular disease, or severe illness were excluded from participation and were transferred to another dialysis unit for intensive care. The patients who were enrolled as subjects (n = 298) underwent stable regular HD with a bicarbonate dialysate. Their underlying diseases of ESRD were as follows: diabetic nephropathy (n = 144), chronic glomerulonephritis (n = 76), hypertensive nephrosclerosis (n = 70), chronic pyelonephritis (n = 3), and unknown (n = 5). Autosomal-dominant polycystic kidney disease patients often have higher Hb levels and lower ESA requirements than ESRD patients with other underlying diseases, and they were excluded.

All patients were on thrice-weekly HD sessions. Blood pressure (BP) was measured with a mercury sphygmomanometer with the patient in the supine position after resting for 10 to 15 min, and mean values for the 1-month period preceding enrollment were used in the statistical analysis. Dry weight was targeted to achieve a normotensive edema-free state. History of previous cardiovascular disease was collected from the medical records. Diabetes was recorded as a history or presence of diabetes and/or a fasting plasma glucose concentration >126 mg/dl or HbA1c concentration >6.5 % or prescription of glucose-lowering agents.

The study period was divided into two phases to evaluate the time-dependent risk of mortality based on subjects’ longitudinal dialysis status. The first phase (12 months) was used to collect laboratory data to evaluate longitudinal dialysis status. The second phase (the subsequent 2 years) was used as a means of evaluating the time-dependent risk of mortality. In order to use longitudinal dialysis condition as a predictive value, laboratory data were not assessed during the second phase. The ESA resistance index (ERI) was calculated by dividing the weekly weight-adjusted ESA dose (IU/kg/week) by the Hb concentration (g/dL), as previously described [12]. Patients were divided into four groups by quartiles (Q) of ERI as follows: Q1 ERI <4.96, Q2 ERI =4.97–8.30, Q3 ERI =8.31–13.1, and Q4 >13.2 IU/kg/week/g/L.

Laboratory and nutrition parameters

Blood sampling after an overnight fast was performed before the mid-week dialysis session day. Serum urea nitrogen, creatinine, calcium, phosphorous, albumin, total cholesterol, ferritin, transferrin saturation, and C-reactive protein (CRP) levels and the Hb concentration were measured by standard laboratory methods. Total calcium values were corrected by the patient’s albumin level. Intact parathyroid hormone (iPTH) was measured by an immunoradiometric assay. Body mass index (BMI) was expressed in kilograms per square meter. Body weight was calculated as dry weight defined as the post-dialysis weight of normotensive patients with no signs of overhydration. Urea kinetics were assessed by measuring a blood-based dialysis parameter, Kt/V [13], and the mean value of the three measurements during each of the 3 months before the start of the study was used in the analysis. The normalized protein catabolism rate (nPCR) was used as an indirect indicator of protein intake and was calculated by using the formula previously described [14].

The GNRI was calculated by modifying the nutritional risk index for elderly patients reported by Yamada et al. [5] as follows: GNRI =[14.89 × albumin (g/dl)] + [41.7 × (body weight/ideal body weight)]. We used the body weight value at the end of the dialysis session to calculate the BMI. When a body weight exceeded the ideal body weight, body weight/ideal body weight was set equal to 1 [6]. The ideal body weight in the present study was calculated using height and a BMI of 22, which has been reported to be associated with the lowest morbidity rate in the Asian population [15].

Anemia treatment strategy

Subjects’ anemia was treated according to the guidelines published in 2004 by the Japanese Society of Dialysis Therapy (JSDT) [16]. The target Hb level was 10–11 g/dL, and epoetin beta was administered to achieve the target Hb level. Intravenous iron was administered whenever a subject’s transferrin saturation was found to be <20 % or the serum ferritin level was found to be <100 ng/mL.

Study outcome

Data for endpoints were obtained from hospital charts and by telephone interviews of patients that were conducted by trained interviewers who were blinded to the study protocol. The primary endpoint of the study was all-cause mortality during the follow-up period from 1 July 2011 to 31 July 2014. The secondary endpoint was cardiovascular death, including death due to heart failure, myocardial infarction, arrhythmia, and stroke. The vital status of the subjects was assessed by searching their electronic dialysis records. Patients were censored if they were alive on 31 July 2014.

Statistical analysis

Normally distributed, unpaired continuous values were expressed as means ± SD and compared by performing an analysis of variance. Nonparametric values were expressed as median values and compared by performing the Kruskal-Wallis test. Categorical values were expressed as percentages and compared by performing the Fisher’s exact test. Logistic regression analysis was used in a model to determine how independent variables predicted the ERI value. We used Spearman’s rank-order correlation coefficients for selected analyses where indicated. Multivariate regression analyses including linear and logistic regression were performed to analyze the data for associations between clinical parameters and ERI. Our fully adjusted model was adjusted for age, gender, dialysis vintage, presence of diabetes, history of cardiovascular disease, Kt/V (single pool), serum total cholesterol and phosphorus levels, transferrin saturation (TSAT), and GNRI.

The survival analysis was based on the Kaplan-Meier curve with subjects censored for death. A log-rank test was used to compare the survival rates of the two groups. Receiver operating characteristic (ROC) analysis was used to select cutoff values of GNRI and ERI for all-cause mortality. We considered some variables that possess P value <0.10 in univariate logistic regression analysis as independent variables for multivariate analyses. A multivariate logistic regression model with adjustment for multiple factors was used to evaluate mortality risk. Results were expressed as odds ratio (OR) and 95 % confidence interval (CI). A p value of <0.05 was considered statistically significant. All statistical analyses were performed by using the SAS version 9.2 software program (SAS Institute Inc., Cary, NC, USA) for Windows personal computers.

Results

Table 1 shows the baseline characteristics of the 298 HD patients according to their ERI quartiles. The mean age was 66.0 ± 13.1 years, 62.8 % were males, 48.3 % were diabetic, and mean dialysis vintage was 5.4 ± 6.2 years. Patients with higher ERI were older and had lower GNRI, BMI, Hb, phosphorus, and TSAT. There was a significant relationship between ERI quartiles and both BMI (p = 0.0002) and GNRI (p = 0.0003).

Table 1 Baseline characteristics of the study population

Table 2 shows the results of the multivariate linear regression analyses between ESA responsiveness and clinical parameters after adjustment for significant confounders by univariate analyses (p < 0.05). Both older age (β = 0.149) and female gender (β = 0.233) were significant and independent predictors of ERI using multivariate regression analyses. In addition, serum total cholesterol (β = −0.194), TSAT (β = −0.193), and GNRI (β = −0.205) were found to be independent negative predictors of ERI. As shown in Fig. 1, the ERI was inversely correlated with the GNRI (r = 0.287, p < 0.0001).

Table 2 Predictors of ERI using univariate and multivariate regression analyses

Fig. 1

Scatter diagram showing the association between the GNRI and ERI

Table 3 shows the likelihood of ESA at the worst hyporesponsiveness using multivariate logistic regression analyses to compare the highest (worst) versus three lower quartiles (as reference) of responsiveness to ESA. The likelihood of belonging to the highest ERI quartile in comparison with the three lower quartiles decreased for every 1 SD higher GNRI (OR = 0.681; 95 % CI = 0.494 –0.930), every 1 SD higher total cholesterol (OR = 0.643; 95 % CI = 0.456–0.892), and every 1 SD higher TSAT (OR = 0.582; 95 % CI = 0.415–0.797).

Table 3 Likelihood of the worst ESA hyporesponsiveness using multivariate logistic regression analyses to compare the highest versus lowest quartiles of ERI

The mean follow-up period was 34.6 ± 6.1 months. Of the 298 enrolled subjects, 51 died during the follow-up period, and the causes of death were cardiovascular disease in 31 subjects (60.7 %) and non-cardiovascular diseases including infection and malignancy in 20 subjects (39.3 %).

ROC analyses were performed to select the optimal cutoff values of GNRI and ERI for increased risk of mortality. Figure 2 shows that the cutoff value of GNRI was 94.9 (area under curve = 0.664, p = 0.0001) and that of ERI was 13.5 (area under curve = 0.574, p = 0.1022). When the subjects were stratified according to their ERI and GNRI values into four groups, those who had low GNRI and high ERI were associated with the highest risk of mortality among the four groups (Fig. 3; log-rank p < 0.001).

Fig. 2

Receiver operating characteristic analysis showing cutoff values of GNRI and ERI for increased all-cause mortality

Fig. 3

Kaplan-Meier curve of survival probability according to the GNRI and ERI. High GNRI and ERI high (GHEH) (n = 32); high GNRI and low ERI (GHEL) (n = 142); low GNRI and high ERI (GLEH) (n = 36); and low GNRI and low ERI (GLEL) (n = 88)

In an attempt to clarify whether the combination of ERI and poor nutritional status could be a better predictor for all-cause mortality of MHD patients, we examined ORs of the combination factor (four groups stratified by ERI and GNRI) by multivariate analyses and compared the OR of GNRI only. As shown in Table 4, GNRI was an independent predictor of all-cause mortality (OR = 0.929; 95 % CI = 0.894–0.964; p < 0.001). In addition, low GNRI and high ERI (GLEH) was the worst predictor of all-cause mortality (OR = 5.325; 95 % CI = 1.888–15.512; p = 0.0016).

Table 4 Multivariate logistic regression analysis of all-cause mortality according to GNRI and combination factor (four groups stratified by ERI and GVRI)

Finally, we performed multivariate analyses of all-cause mortality using confounders including quartiles of ERI and GNRI. As shown in Table 5, Q1 and Q2 of GNRI were significant predictors of all-cause mortality, whereas Q3 of ERI only was an independent factor for mortality.

Table 5 Multivariate logistic regression analysis of all-cause mortality according to quartiles of ERI or GNRI

Discussion

The results of this study showed that ESA responsiveness and nutritional status were associated with increased risk of death after adjustment for potential confounding factors, and, in particular, they showed that patients with high ERI and low GNRI had the highest risk of mortality. These findings suggest that ESA responsiveness and nutritional status should be considered simultaneously when evaluating the relationship between renal anemia and malnutrition and mortality risk in MHD patients.

ESA responsiveness is relative, not absolute, and the two determinants of responsiveness are ESA dose and Hb level. Previous studies have demonstrated associations between ESA resistance and high all-cause mortality and high cardiovascular mortality in ESRD patients [17–19]. ERI is considered a useful means of ESA resistance, but the follow-up periods in the previous studies were relatively short (i.e., 1 or 2 years). The results of this study showed that a combination of high ERI and low GNRI was associated with increased mortality risk during a mean follow-up period of 34.6 ± 6.1 months. Because the ERI and GNRI were calculated from the baseline data, the ERI and GNRI may have changed during the long-term follow-up period.

ESA responsiveness varies widely according to a number of clinical factors. The results of the multivariate logistic regression analysis showed that the likelihood of belonging to the highest ERI quartile in comparison with the three lower quartiles decreased for every 1 SD higher GNRI (OR = 0.681; 95 % CI = 0.494–0.930), every 1 SD higher total cholesterol (OR = 0.643; 95 % CI = 0.456–0.892), and every 1 SD higher TSAT (OR = 0.582; 95 % CI = 0.415–0.797). These results suggest that nutritional markers are independent predictors of the ESA responsiveness of MHD patients.

Assessment and monitoring of protein and energy nutritional status are essential to prevent, diagnose, and treat uremic malnutrition in dialysis patients [20]. Many nutritional screening tools have been developed for use in the elderly, children, hospitalized patients, patients in the community, and patients with cancer or infections, and some of them can be safely and easily applied to MHD patients as well [21]. The GNRI was one of the parameters that was a predictor of increased overall mortality in the present study. The GNRI is calculated from only a few objective components, including serum albumin and BMI, and is a simple nutritional screening tool [5]. Kobayashi et al. [6] reported finding that HD patients whose GNRI was <90 had a poorer nutritional status in terms of lower serum levels of albumin and creatinine than HD patients whose GNRI was >90. However, it is very difficult to determine the ideal cutoff value. Park et al. [7] found that GNRI may be a significant predictor of increased mortality in Korean HD patients. More recently, Panichi et al. [22] obtained results closer to those obtained by Yamada et al. and indicated that values below 92 may be a strong indicator of an unfavorable outcome.

Various mechanisms, including derangements in muscle, gastrointestinal, hematopoietic and immune systems, and abnormal activation of the inflammatory process may explain the link between PEW and mortality in renal failure [23]. Low serum albumin and cholesterol levels are indicators of PEW but may not be causally responsible for negative outcome. The reduction in muscle mass, namely sarcopenia, observed in PEW may be due to uremic toxins or procatabolic conditions (metabolic, hormonal, or neuropathic abnormalities). Muscle wasting may in turn lead to reduced skeletal, respiratory, and cardiac muscle function and cause functional insufficiency [24].

The present study had several limitations. First, baseline data were used to define exposure categories in this cohort, and we were unable to examine the effect of changes from the baseline category during the follow-up period. Second, since this was a single center cohort study, the data and results may not be representative of MHD patients in other countries, particularly since the ERI and GNRI are not widely used in Japan. Third, the cohort size was small and the mortality rate was low. However, the follow-up period of our cohort was longer than in previous studies.

Conclusions

A combination of a high ERI and low GNRI was an independent predictor of increased all-cause mortality in MHD patients. In future studies, greater attention should be focused on background factors that increase responsiveness to ESA in addition to efforts to improve nutritional status.