Abstract
Background and objectives We hypothesized that longitudinal changes in uric acid (UA) may have independent associations with changes in nutritional parameters over time and consequently, long-term survival of patients on maintenance hemodialysis.
Design, setting, participants, & measurements We conducted a retrospective, longitudinal cohort study of a clinical database containing the medical records of patients on maintenance hemodialysis receiving dialysis between June of 1999 and December of 2012 in a single center; 200 patients (130 men and 70 women) with a median age of 69.0 (interquartile range, 59.3–77.0) years old were included in the study. Dietary intake, biochemical markers of nutrition, anthropometric measurements, and UA levels were recorded at 0, 6, 12, 18, 24, 30, and 36 months followed by 15 additional months of clinical observations. The patients were followed until January 31, 2015 (median follow-up was 38.0 [interquartile range, 30.0–46.8] months).
Results In a linear mixed effects model adjusted for baseline demographics and clinical parameters, each 1.0-mg/dl longitudinal increase in UA was associated with a 13.4% slower rate of decline in geriatric nutritional risk index (GNRI) levels over 3 years of observation (95% confidence interval [95% CI], 0.11 to 0.39; P<0.001 for UA × time interaction). UA remained associated with the rate of change in GNRI, even after controlling for C-reactive protein. During the follow-up, 87 (43.5%) all-cause and 38 (19.0%) cardiovascular deaths were reported. For each 1.0-mg/dl increase in serum UA over time, the multivariate adjusted all–cause mortality hazard ratio using Cox models with the effect of time-varying risk was 0.83 (95% CI, 0.74 to 0.95; P<0.01), which continued to be significant, even after including the baseline GNRI levels in this model: 0.89 (95% CI, 0.79 to 0.98; P=0.02).
Conclusions Longitudinal changes in serum UA seem to track with changes in nutritional status over time, and these changes are associated with survival of patients on maintenance hemodialysis. An increase in serum UA levels over time is accompanied by improvement of nutritional status and lower mortality rate.
- uric acid
- hemodialysis
- inflammation
- nutrition
- geriatric nutritional risk index
- mortality
- Cohort Studies
- Humans
- Longitudinal Studies
- Nutritional Status
Introduction
The importance of serum uric acid (SUA) levels for patients on maintenance hemodialysis (MHD) has not been well established. Uric acid (UA) is a purine metabolite that is formed in the liver, and 65%–75% of its daily produced amount is excreted by the kidneys (1). When significant amounts of UA accumulate in patients with ESRD because of lack of renal clearance, the mean UA removal is approximately 1 g per hemodialysis (HD) session, even with high-flux hemodialysers (2). Consequently, hyperuricemia is common in patients with ESRD undergoing treatment by either HD (3,4) or peritoneal dialysis (PD) (5–7). In the general population, hyperuricemia has been shown to be associated with an increased risk of hypertension (8,9), prehypertension (10), peripheral arterial disease (11), diabetes mellitus (12), CKD (13), and cardiovascular disease (CVD) (8,10,14). Conversely, in the HD population, hyperuricemia has been shown as a good nutritional marker (15) that can also predict lower rates of upcoming hospitalizations (15) as well as lower all–cause and cardiovascular death risks (15,16). Although the data about the mortality predictive role of higher SUA levels in the HD population are somewhat controversial (3,4,15–18), most of the studies conducted on this topic presented low SUA as a mortality risk factor in patients on MHD (3,4,15,16,18).
However, to our knowledge, there is no information in the literature on the clinical effect of longitudinal changes in SUA levels on the prognosis of patients on HD. We hypothesized that longitudinal changes in SUA may have independent associations with changes in nutritional parameters over time and consequently, long-term survival in patients on MHD. This study was, therefore, designed to determine the association of changes in nutritional and inflammatory parameters over time on changes in SUA and subsequent clinical outcomes in patients on MHD.
Materials and Methods
Patients
We conducted a retrospective, longitudinal cohort study of a clinical database containing the medical records of patients with ESRD receiving dialysis treatment from June 7, 1999 to July 31, 2015 in a single dialysis center (Assaf Harofeh Medical Center, Zerifin, Israel). This study was approved by our Institutional Ethics Committee.
The database includes patient demographics, medications, clinical laboratory values, dialysis treatment records, records of dietitian consultations, and all comorbidities available to the dialysis unit. Patients included in the analysis were those ≥18 years old who received thrice weekly in–center HD for ≥1 year from June 7, 1999 to December 5, 2012. The first (baseline) visit for each patient was the calendar date on which the patient’s dialysis duration was >60 days. Study measurements were selected from the database at 0, 6, 12, 18, 24, 30, and 36 months followed by 15 additional months of clinical observations. Follow-up time began on the date of entry into the cohort. The end date of the follow-up was January 31, 2015. A flowchart of the study is presented in Figure 1. Patients were censored at the time of death, renal transplantation, or transfer to PD or another dialysis unit (loss of follow-up). From 2432 patients receiving dialysis during the cohort period, 990 were ineligible (436 patients treated with PD and 554 patients on MHD for <1 year were excluded). From 1442 potential participants, 1242 patients were excluded because of incomplete data (if we were unable to record SUA, C-reactive protein [CRP], dietary intake calculation, and anthropometric measurements every 6 months during the longitudinal part of the study) (Figure 1). However, a comparison of included patients versus patients excluded because of incomplete data (Table 1) showed that the two groups were almost identical with regard to age, comorbidity index, BP, Kt/V, dietary intake, anthropometric measurements, and UA, creatinine, cholesterol, triglyceride, transferrin, CRP, and geriatric nutritional risk index (GNRI) levels. Although statistically significant, the differences in sex, dialysis vintage, diabetes status, vascular access type, serum albumin, hemoglobin levels, mortality, and use of diuretics or allopurinol were of small magnitude. In total, 200 patients (130 men and 70 women) with a median age of 69.0 (interquartile range, 59.3–77.0) years old who had available UA and pertinent covariate data were included in the study. These patients contributed a total of 634 patient-years of at-risk time, and the median follow-up was 38.0 (interquartile range, 30.0–46.8) months. During the follow-up, in total, 87 (43.5%) all-cause and 38 (19.0%) cardiovascular deaths were reported. The causes of renal failure included diabetic kidney disease (50.0%), hypertension (32.0%), GN (4.0%), autosomal dominant polycystic kidney disease (6.0%), obstructive uropathy (4.0%), and other diseases or unknown causes (4.0%). All patients underwent regular dialysis via their vascular access for 4–5 hours three times per week at a blood flow rate of 250–300 ml/min and a dialysis solution flow rate of 500 ml/min. Dialysis treatments were performed with low–flux biocompatible dialyzer membranes with surface areas of 1.4–1.8 m2.
Flow diagram of the study. PD, peritoneal dialysis.
Baseline characteristics of included patients (n=200) versus patients excluded because of missing data (n=1242)
Dietary Intake
Dietitian records on the basis of self–completed food diaries with continuous 3-day dietary histories (including a dialysis day, a weekend day, and a nondialysis day) were used to calculate the dietary intake. The dietary energy and protein intake were calculated and normalized for ideal body weight according to the European Best Practice Guidelines (19). Ideal weight in this study was calculated from the Lorentz equations differently for men and women. Dietary intake was calculated using computerized analysis (a Disk Operating System–based program MANA specially adapted for data entry and analysis of food intake records). This program was developed by the Israeli Food and Nutrition Administration on the basis of the Food Intake Analysis System of the US Department of Agriculture (20) and especially adapted to meet the needs of the Israeli population.
Dietary protein intake was also estimated by using the normalized protein catabolic rate (nPCR) calculation from the patient’s urea generation rate by urea kinetics modeling (21). Single–pool model urea kinetics were used to estimate the nPCR.
Anthropometric Measurements
Body mass index (BMI), triceps skinfold thickness (TSF), midarm circumference (MAC), and calculated midarm muscle circumference (MAMC) were measured as anthropometric variables. TSF was measured with a conventional skinfold caliper using standard techniques. MAC was measured with a plastic measuring tape. MAMC was estimated as follows: .
GNRI
GNRI was calculated from the patient’s serum albumin, body weight, and height by using the equation developed by Bouillanne et al. (22) and modifying the nutritional risk index for elderly patients.
Comorbidity Index and Clinical Outcomes
We determined the comorbidity index, which was recently developed by Liu et al. (23) and validated specifically for populations of patients on dialysis, as a measure of comorbid conditions.
Hospitalization data were obtained for all 200 patients on HD. Hospitalization was defined as any hospital admission that included at least one overnight stay in the hospital. The access-related hospitalizations were not included in the hospitalization data.
CVD was defined as myocardial infarction requiring coronary artery procedures, such as angioplasty or surgery, cerebrovascular accident, or peripheral vascular disease requiring angioplasty, bypass, or amputation. Cardiovascular mortality was defined as death resulting from coronary heart disease, sudden death, stroke, or complicated peripheral vascular disease.
Laboratory Evaluation
Blood samples were obtained from nonfasting patients on a midweek day predialysis, with the exception of postdialysis serum urea nitrogen to calculate urea kinetics. Albumin was measured using the bromocresol green method. All biochemical analyses, including UA, creatinine, urea, albumin, transferrin, complete blood count, triglycerides, and total cholesterol, were measured by an automatic analyzer. Additionally, serum high–sensitivity CRP was measured by a turbidimetric immunoassay.
Statistical Analyses
Data are expressed as means±SDs, medians and interquartile ranges (quartiles 1–3) for variables that did not follow a normal distribution, or frequencies as noted.
Longitudinal data were analyzed by using the Mixed model. Base models were adjusted for age, sex, diabetes status, dialysis vintage, comorbidity index, Kt/V, diuretic use, and vascular access type. F tests were used to assess the significance of the fixed effects, and P values of <0.05 were considered significant. To evaluate whether changes in time of UA associate the changes in various nutritional parameters, we included in each base model terms for individual UA by time interactions. We adjusted for baseline levels of UA when analyzing changes in nutritional parameters to minimize the likelihood of regression to the mean.
Normally distributed continuous variables were compared between the two groups using a two–sided t test, with chi-squared tests used for categorical variables and nonparametric Mann–Whitney U tests used for not normally distributed continuous variables. Associations between baseline SUA and nutritional parameters were assessed using Pearson correlation coefficients or Spearman rank order correlation coefficients in cases of skewed distribution of data. Survival analyses were examined in Cox models with UA as a time-varying predictor. Univariate and multivariate Cox regression analyses are presented as hazard ratios (HRs) and 95% confidence intervals (95% CIs).
All statistical analyses were performed using SPSS software, version 16.0 (IBM SPSS, Chicago, IL).
Results
The baseline characteristics of the cohort are shown in Table 1. The 200 patients on prevalent HD (median age, 69.0 years old; interquartile range, 59.3–77.0 years old) selected for this study included 35% women. One half of the participants (51.5%) had diabetes mellitus and a median dialysis vintage of 8.0 months. nPCR as a marker of daily protein intake, laboratory nutritional markers (albumin, creatinine, triglycerides, and transferrin), body composition parameters (BMI and TSF), and the overall nutritional status assessed by GNRI were positively correlated to SUA at baseline.
We examined the individual differences in the rate of change in GNRI using a linear mixed effects model controlling for age, sex, diabetes mellitus status, dialysis vintage, vascular access type, Kt/V, diuretic use, and comorbidity index. GNRI levels decreased by about 0.47 U per 3 years (95% CI, −0.65 to −0.29; P<0.001). We next added a term for UA and its interaction with time to the previous model (Table 2) to examine whether UA was associated with the rate of change in GNRI. In this model (model 2), both UA and GNRI levels decreased (UA by 0.38 mg/dl per 3 years; 95% CI, −0.77 to <0.01 and GNRI by 1.87 U per 3 years; 95% CI, −2.69 to −1.06), but each 1.0-mg/dl increase in UA at baseline was associated with a slower decrease in GNRI levels of 0.25 U per 3 years (P<0.001 for UA × time interaction). In practical terms, a person with UA 1 mg/dl above the mean at baseline would show about a 13.4% slower rate of decline in GNRI levels. Next, we examined whether controlling for CRP confounded the association of UA with the rate of change in GNRI over time. After adding CRP and its interaction with time to the previous model, UA remained associated with changes in GNRI over 3 years of observation. These results suggest that UA is independently associated with the rate of change in GNRI over time.
Associations of longitudinal uric acid changes with changes in geriatric nutritional risk index (slopes) over 36 months on the basis of a mixed effects model
We also examined the associations of UA with longitudinal changes in nutritional parameters (slopes) over 36 months including fixed parameters, such as age, sex, Kt/V, diuretic use, diabetes status, dialysis vintage, vascular access type, and comorbidity index (Table 3). Longitudinally, SUA levels were associated with most of the nutritional markers studied at any given time point, controlling for demographic and clinical parameters. In addition, there was the significant association of the UA × time interaction on the slope of albumin, creatinine, transferrin, cholesterol, and triglycerides over 36 months of observation. Specifically, in a person with SUA 1 mg/dl above the mean at baseline, the decline in serum albumin was reduced by almost 15.7%, and the decline in serum creatinine was reduced by 18.9%. Additionally, UA levels were associated with the rates of change in BMI and TSF (slopes) over time as shown in Table 3. We did not observe any association between SUA and CRP levels over time.
Associations of uric acid with nutritional parameters and their rates of changes (slopes) over 36 months on the basis of a mixed effects model
Next, we examined the trajectory of longitudinal changes in SUA and GNRI levels in groups of patients divided by tertiles of the absolute change in SUA (ΔSUA) from baseline to 36 months or from baseline to the final measurement of SUA for those dropped out during the longitudinal part of the study (Figure 2). GNRI in the first tertile of ΔSUA (the group of patients with a decrease in SUA to <−0.8 mg/dl, with a corresponding rate of change of −0.21; 95% CI, −0.25 to −0.16; P<0.001) exhibited a significant decline over time (−0.80; 95% CI, −1.11 to −0.48; P<0.001). However, in the second group (the group with SUA decline over time that was more gradual, with ΔSUA ranging from −0.8 to 0.2 mg/dl and a corresponding rate of change of −0.05; 95% CI, −0.08 to −0.02; P=0.002) and the third tertile of ΔSUA (the group with an increase in SUA to >0.2 mg/dl and a corresponding rate of change of 0.08; 95% CI, 0.03–0.13; P=0.002), GNRI did not decline significantly over the 36-month period.
Comparison of changes in geriatric nutritional risk index (GNRI) levels over a period of 36 months in the tertiles of absolute changes in serum uric acid (ΔSUA) over the same period. The first tertile is uric acid decrease to <−0.8 mg/dl (circles and solid line). In this group, SUA declined significantly over time: the estimate was −0.21 (95% confidence interval, −0.25 to −0.16; P<0.001). The second tertile is change of SUA between −0.8 and 0.2 mg/dl (circles and dashed line). SUA decline was statistically significant but more gradual (−0.05; 95% confidence interval, −0.08 to −0.02; P=0.002). The third tertile is uric acid increase to >0.2 mg/dl (triangles and dotted line); in this group, SUA increased significantly over time (0.08; 95% confidence interval, 0.03 to 0.13; P=0.002). In the whole cohort, GNRI declined during the study period (−0.33; 95% confidence interval, −0.49 to −0.17; P<0.001; crude analysis), whereas a statistically significant decline was observed only in the first tertile of the ΔSUA over time (−0.80; 95% confidence interval, −1.11 to −0.48; P<0.001), and GNRI did not change in the second (−0.15; 95% confidence interval, −0.43 to 0.12; P=0.27) and third (−0.06; 95% confidence interval, −0.30; 0.19; P=0.66) tertiles of ΔSUA. Data are presented as means with 95% confidence intervals.
Overall, for each 1.0-mg/dl increase in SUA, the crude all–cause mortality (87 of 200 participants died during the study period) HR using Cox models with the effect of time-varying risk was 0.91 (95% CI, 0.83 to 1.01; P=0.06), whereas after multivariate adjustments for all above–mentioned demographic and clinical covariates, HR was 0.83 (95% CI, 0.74 to 0.95; P<0.01), which continued to be significant even after including the baseline GNRI levels in this model (0.89; 95% CI, 0.79 to 0.98; P=0.02). Time-varying UA did not show an association with CVD-related hospitalization (75 cardiovascular events were reported during the follow-up) and mortality (38 cardiovascular deaths were reported) (Table 4).
Hazard ratios (95% confidence intervals) of all-cause mortality, cardiovascular disease–related mortality, and hospitalization associated with time–varying serum uric acid levels in unadjusted and multivariable–adjusted Cox models
Discussion
In this study, we wished to determine whether changes of SUA levels are associated with parallel changes in nutritional parameters over time in a cohort of patients on MHD. We provide novel evidence that changes in SUA levels independently reflect changes in nutritional status in our population. In addition, these changes can predict long-term survival in patients on MHD.
Regarding the effects of longitudinal associations, we found that SUA, most of the studied nutritional parameters, and nutritional status as a whole expressed by GNRI declined simultaneously over 36 months of observation. Only a few existing longitudinal studies have investigated changes in SUA levels in patients on MHD. In a 6-month longitudinal study with three 3-month repeated blood samplings, no statistically significant differences among the UA levels were found in 48 patients on HD (24). However, in agreement to our findings, in a 2-year study, on the basis of four 6-month blood samples, UA significantly decreased over the study period in 117 patients on HD (25). With regard to nutritional status, most of studies have shown that MHD associates with a significant decline in certain laboratory nutritional parameters, including albumin (26,27), as well as anthropometric measurements (27,28). GNRI, an objective and simple nutritional index, has also been reported to decrease significantly with time in patients on MHD in a recent prospective longitudinal study (29). The longitudinal behavior of SUA may be explained as follows: it is a proposed nutritional marker in the MHD population as reported in a recently published prospective observational study (15), and changes in SUA levels over time may simply follow changes in nutritional status. In addition, overproduction of UA is usually associated with visceral fat accumulation (30,31) as part of the metabolic syndrome, the prevalence of which is high in the ESRD population (32). Therefore, higher and stable over time SUA levels may express a proportion of metabolic syndrome in our cohort. Although not all markers of the metabolic syndrome were measured in this study, the significant associations observed between SUA changes over time and the rates of change (slopes) in BMI, TSF (as a surrogate of fat mass), and triglycerides (statistically significant SUA × time interactions) may support this assumption. Metabolic syndrome is associated with better nutritional status (expressed as better Subjective Global Assessment scores) in patients on MHD as reported by Xie et al. (33).
One must also consider whether changes in inflammatory markers observed in our study were relevant to the association with SUA changes and in parallel to GNRI changes over time. Although there is in vitro and experimental evidence of the proinflammatory effects of UA in human vascular smooth muscle and endothelial cells (34,35), the relationship of SUA with the inflammatory response in the HD population has not been sufficiently studied and is still controversial. SUA was positively associated with serum CRP levels in an incident HD population (4). Conversely, no significant association of SUA with CRP was found in the same population by another study (17). Of interest, monocytes from patients with ESRD purified from peripheral blood mononuclear cells and incubated with monosodium urate crystals showed a reduced secretion of proinflammatory cytokines, including IL-6, compared with healthy controls (36). We found no longitudinal study in the literature regarding the relationship between SUA and inflammation. We showed that SUA changes independently associate with the rate of change in GNRI levels over time, even after controlling for such covariates as CRP. We suggest, therefore, that changes in SUA over time associated with longitudinal changes in GNRI are independent from inflammation. However, because IL-6, the best inflammatory biomarker in the MHD population (37), was not measured in this study, we cannot completely rule out the confounding role of inflammation in the association between SUA and nutritional status.
In the context of associating SUA with prognosis, our study expands the cross-sectional observations (15,16) showing that changes of SUA over time may predict outcomes in patients on MHD. We did not find use of a longitudinal design in the few studies that we did find providing data about the association of SUA and clinical outcome in patients on MHD (3,4,15–18). Instead, a relatively large number of studies has used longitudinal data to analyze the association of changes in albumin (38), BMI (39), and TSF (28) over time with prognosis. This probably supports our data given the estimates of the regression coefficients in Table 3 expressing associations of SUA and the above–mentioned nutritional parameter changes over time. Thus, in our study, we have shown that the longitudinal increase in SUA over time is correlated with a more gradual decrease in GNRI and consequently, better all–cause mortality but not correlated with CVD mortality and first CVD event. Regarding the association between changes of SUA over time and CVD morbidity and mortality, prospective observational but not longitudinal studies present strong evidence on the prognostic power of SUA (8–12,10,14), even in patients on MHD (15,16). Therefore, this question should be examined in larger and prospective longitudinal studies.
Finally, multivariable adjustments in Cox models with the effect of time-varying risk, including baseline GNRI levels, did not abolish the prognostic value of longitudinal SUA changes on mortality. This suggests that, together with being associated with nutritional status, additional mechanisms contribute to the prognostic significance. In this context, the antioxidant properties of SUA (40) may play a role. Oxidative stress is highly linked to inflammation (41) and poor survival (42) in patients on MHD.
In interpreting the results of our study, we should be concerned with whether regression to the mean is relevant to the association of longitudinal SUA changes on changes in GNRI levels (slopes) over time. The linear mixed model used in our study to analyze longitudinal data models the covariance matrix, which gives a high statistical power and adjusts each subject’s follow-up measurement according to the baseline measurement. The use of analysis of covariance in such a way is accepted as one of the practical ways to deal with regression to the mean in statistical analysis (43).
A number of limitations inherent to observational research applies to our analysis. First, this study only used an observational approach without manipulating exposure factors, and therefore, no definitive cause and effect relationship can be derived for any of the risk factors analyzed. Second, selection bias is typical for retrospective studies. A difference in incidence of the outcome of interest between those who participated and those who did not would give biased results. Despite the similarity between the study sample and excluded patients (Table 1), selection bias can still affect the generalizability of our findings to the wider MHD population. Third, because of the retrospective design, we could not obtain a number of parameters of interest (e.g., residual renal function, waist circumference, hydration status, IL-6, and glutathione). Fourth, the overall nutritional status was assessed by GNRI, which may not be the ideal nutritional marker in the MHD population. Despite the observed significant associations of GNRI with dietary intake, laboratory markers of nutrition, body composition, and prognosis (29,44), the only validation study of GNRI as a nutritional marker in the MHD population was conducted in Japanese patients on HD (45). Fifth, the proportion of cardiovascular deaths may be underestimated, because the causes of death were extracted from patient records and were not confirmed by autopsies. Despite these limitations, the wide array of nutritional parameters, which included anthropometrics, biochemical markers, and dietary intake, longitudinal design, and long–term follow-up used in this study strengthens the results herein.
In summary, our study shows that longitudinal changes in SUA seem to track changes in the nutritional status over time and that these changes are associated with survival of patients on MHD. An increase in SUA levels over time is accompanied by improvement of nutritional status and lower mortality rate.
Disclosures
None.
Acknowledgments
The authors thank Dr. Mechael Kanovsky for his language editing, which has greatly improved the manuscript.
Footnotes
Published online ahead of print. Publication date available at www.cjasn.org.
- Received September 30, 2015.
- Accepted February 23, 2016.
- Copyright © 2016 by the American Society of Nephrology