Hyperuricemia as a Predictor of Progression of Chronic Kidney Disease: A Matched Cohort Analysis ()
1. Introduction
Chronic kidney disease (CKD) is affected by multiple risk factors for disease progression [1] [2] , and it is extremely important to identify risk factors for the acceleration of the CKD progression. The impact of the serum uric acid (UA) on renal prognosis of CKD patients has recently been attracting attention [3] , and our group has reported finding that the serum UA value can serve as a predictor of progression of kidney disease in biopsy-proven nephrosclerosis patients [4] . However, because the serum UA elevation due to the decreased renal UA clearance in CKD patients makes it difficult to assess the role of UA in the progression of CKD, the predictive value of UA in regard to the progression of CKD remains a matter of controversy [5] - [11] .
A high serum UA level is associated with obesity [12] [13] and is a risk factor for carotid atherosclerosis [14] . The aim of the present study was to elucidate the predictive value of hyperuricemia (HU) in relation to kidney disease progression in CKD patients.
2. Methods
2.1. Patient Selection
Abdominal computed tomography (CT) and carotid ultrasonography were performed on 201 CKD outpatients who attended the Kidney Center of Tokyo Women’s Medical University Hospital between August 2006 and August 2007. CKD was diagnosed according to previously described criteria [15] . Among these 201 patients, the following patients were excluded from participation: 1) patients who had been treated with antihyperuricemic agents; 2) patients with no UA data. The remaining 122 patients were ultimately enrolled as subjects of the present study (Figure 1).
Figure 1. Flow chart of patient selection. The 79 patients who did not meet the entry criteria were excluded from the 201 patients screened, and the other 122 patients were enrolled in this study.
The subjects’ human rights and method of protecting personal information were well considered. All of the relevant and responsible staff adhered to the Helsinki Declaration (amended October 2013) and the Ethical Guidelines for Clinical Studies (revised July 31, 2008, referred to hereafter as the Clinical Studies Ethical Guidelines) in the execution of this study. The study was approved by the Medical Ethics Committee of Tokyo Women’s Medical University (#4188). All participants gave their informed consent to participate in the study.
2.2. Measurements of Covariates
At a regular ambulatory visit, each subject underwent anthropometric and physical examinations that included blood pressure (BP), height, body weight, visceral fat area (VFA), subcutaneous fat area (SFA), and maximum carotid intima- media thickness (IMT). BP was measured in triplicate with a mercury sphygmomanometer, and the average value was used in the analysis. VFA and SFA were measured by CT, and IMT was measured by carotid ultrasonography. Details of the measurement techniques are described below in this section. All biochemical analyses were performed on samples obtained from fasted subjects. Serum creatinine (Cr) levels were measured enzymatically. The estimated glomerular filtration rate (eGFR) was calculated by using the previously described formula for Japanese patients [16] . Urinary albumin was measured by a latex agglutination test [17] . Urinary albumin excretion was expressed as the urine albumin-to-creatinine ratio (UACR). We also assessed the impact of concomitant drug use and comorbidities at entry [18] [19] . The concomitant drugs being taken by the patients were: antihypertensive drugs, diuretics, and drugs for the treatment of hyperuricemia, dyslipidemia, and diabetes mellitus.
The baseline parameters assessed in this study were: age, gender, systolic blood pressure (SBP), diastolic blood pressure (DBP), mean blood pressure (MBP), BMI, VFA, SFA, IMT, eGFR, hemoglobin, serum albumin, UA, total cholesterol (TC), LDL cholesterol (LDL-C), HDL cholesterol (HDL-C), and triglyceride (TG) values, and UACR. Comorbidities are defined below. The subjects were followed up until July 2016.
2.3. Abdominal CT Examination
The multi-slice CT examination was performed by using a GE Light Speed. An index image was obtained before scanning, and the umbilicus to L4-5 level was identified [20] [21] . VFA at this level correlated extremely well with whole intra-abdominal fat content (r = 0.94). Horizontal images were obtained at 400 mA and 120 kVp with a scan time of 1.0 s. The range of CT values covered optimal CT numbers for fat tissue, i.e., from −150 to −40.14 Data were stored and analyzed with GE advantage workstation Ver.4.0.
2.4. Carotid Ultrasonography
An experienced examiner blinded to the patient data performed high-resolution duplex carotid ultrasonography with a 7.5-MHz duplex scanner (Aplio XG; Toshiba, Tokyo, Japan). The common and internal carotid arteries were scanned cross-sectionally and longitudinally to estimate the presence and distribution of atherosclerotic plaques. The entire length of both common carotid arteries and both internal carotid arteries up to approximately 20 mm distal to the tip of the carotid bifurcation were scanned. Maximum intima-media thickness (IMT) measurements were obtained to identify the thickest region of the arterial wall according to an international consensus report [22] .
2.5. Definition of Comorbidities
Hypertension: Being treated with an oral antihypertensive agent, SBP ≥ 140 mmHg, DBP ≥ 90 mmHg.
Hyperuricemia (HU): Being treated with an oral antihyperuricemic agent, serum UA level ≥ 7.1 mg/dL.
Hypercholesterolemia: Being treated with an oral antidyslipidemic agent, serum TC level ≥ 220 mg/dL, or serum LDL level ≥ 140 mg/dL.
Hypertriglyceridemia: Being treated with an oral antidyslipidemic agent or serum TG level ≥ 150 mg/dL.
Diabetes mellitus: Being treated with an antidiabetic agent or a history of diagnosis with diabetes mellitus.
2.6. Outcome Evaluation (Endpoint)
The outcome variable of interest was kidney disease progression, which was defined as a ≥50% decline in eGFR from baseline (≥50% eGFR decline) or ESRD requiring dialysis.
2.7. Statistical Analysis
Continuous variables are reported as the mean ± standard deviation (SD), and categorical variables are reported as percentages unless otherwise stated. We compared participant outcomes by performing an unpaired t-test or chi-square test, or Fisher’s exact test. The correlations between UA and the other variables were assessed by using Pearson’s correlation coefficients. The optimal cut-off serum UA value for discriminating ≥a 50% eGFR decline or ESRD during follow-up examinations was determined by performing a receiver operating characteristic (ROC) analysis. The patients were divided according to their UA level at baseline into an HU group, i.e., a group whose UA value was ≥7.1 mg/dL, and a non-HU group, i.e., a group whose UA value was <7.1 mg/dL. Prognostic variables for renal outcome were assessed by the univariate and multivariate Cox proportional hazards method. We included covariates for age, sex, VFA, IMT, eGFR, UACR, and comorbidities, including HU, at baseline in the Cox proportional hazards models. Variables with P-values less than 0.1 in the univariate model were included in the multivariate model. Renal outcome which was a ≥50% eGFR decline or ESRD, and interval estimates between the HU group and the non-HU group was calculated by the Kaplan?Meier method and evaluated by the log-rank test.
To further assess whether associations were consistent across clinically matched subgroups, we fit propensity score-matched models that included several potential modifying variables (sex, eGFR, MBP, VFA), and performed subgroup analyses of the groups. The caliper-matching method was used with a maximum tolerance level of 0.1. Standardized differences were calculated to assess the appropriateness of matching, and 95% confidence intervals (CIs) were calculated. P values <0.05 were considered statistically significant. All statistical analyses were performed by using the JMP Pro ver.12.1.0 software program (SAS Institute, Cary, NC, USA).
3. Results
3.1. Patients
The 122 subjects consisted of 48 males and 74 females, and their mean age at baseline was 58.7 ± 13.9 years (range 24 - 84 years). The mean SBP was 124.9 ± 8.1 mmHg, DBP 76.4 ± 6.2 mmHg, MBP 92.6 ± 6.6 mmHg, BMI 23.9 ± 4.1 kg/m2, VFA 117.1 ± 62.1 cm2, IMT 1.44 ± 0.81 mm, UACR 238.7 ± 548.1 mg/g Cr, and eGFR 63.7 ± 19.8 mL/min/1.73 m2. The concomitant drug data showed that 75 subjects were being treated with an antihypertensive agent, 43 with an antidyslipidemic agent, 17 with an antidiabetic agent, 20 with a corticosteroid, 7 with an immunosuppressing agent, and 19 with a diuretic. The comorbidity data showed that 76 patients had hypertension, 20 had HU (UA ≥ 7.1 mg/dL), 73 had hypercholesterolemia, 78 had hypertriglyceridemia, and 26 had diabetes mellitus. The overall follow-up period was 8.3 ± 2.4 years, and 16 patients had reached the endpoint (≥50% eGFR decline or ESRD) during the follow-up period.
3.2. Serum UA Cut-Off Value as a Predictor of Kidney Disease Progression
We performed ROC analyses to identify the optimal UA cut-off value for discriminating a ≥50% eGFR decline or ESRD during the follow-up examination, and the results showed that it was 7.1 mg/dL (AUC = 0.74, sensitivity = 62.5%, specificity = 90.6%, Figure 2).
3.3. Comparison between the Clinical Findings According to Serum UA Value Groups in the Total Cohort
We compared the clinical characteristics of the two groups established according to UA value at baseline (Table 1). The baseline serum UA (7.8 ± 0.6 vs. 5.0 ± 1.2 mg/dL, P ≤ 0.0001), SBP (130.0 ± 7.8 vs. 123.9 ± 7.8 mmHg, P = 0.002), DBP (80.5 ± 5.2 vs. 75.6 ± 6.1 mmHg, P = 0.001), MBP (97.0 ± 5.9 vs. 91.7 ± 6.4 mmHg, P = 0.0009), VFA (158.8 ± 75.4 vs. 108.9 ± 56.0 cm2, P = 0.0008), and UACR (649.7 ± 1035.8 vs. 158.1 ± 343.6 mg/g Cr, P = 0.0002) values of the HU group were significantly higher than in the non-HU group, and the eGFR (43.8 ± 19.0 vs. 67.6 ± 17.6 mL/min/1.73 m2, P ≤ 0.0001) and HDL-C (51.3 ± 15.5 vs. 60.7 ± 16.8 mg/dL, P = 0.02) values of the HU group were significantly lower
Table 1. Patient characteristics according to baseline serum uric acid levels (total cohort n = 122).
Continuous values are expressed as means ± standard deviation. Count data are expressed as percentages. Abbreviation: n, number; HU, Hyperuricemia; UA, Uric Acid; SBP, systolic blood pressure; DBP, diastolic blood pressure; MBP, mean blood pressure; BMI, body mass index; IMT, maximum carotid intima-media thickness; eGFR, estimated glomerular filtration rate; UACR, urine albumin-to-creatinine ratio.
Figure 2. Receiver operating characteristic analysis to identify the optimal serum uric acid cut-off value for predicting an eGFR decline by ≥50% from baseline or end-stage renal disease during the follow-up examination period.
than in the non-HU group. The percentage of males (70.0% vs. 33.3%, P = 0.002) was significantly higher in the HU group than in the non-HU group.
We performed a Kaplan-Meier analysis to assess kidney survival, with a ≥50% eGFR decline or ESRD used as the end-point. The kidney survival curves showed that the kidney survival rate of the CKD patients in the HU group was significantly lower than in the non-HU group (Figure 3(a)). At the 100-month follow- up examination, a decrease in eGFR value of at least a 50% was observed in 64.2% of the HU group (log rank, P ≤ 0.0001). The results of the sex-stratified analyses showed a significantly lower kidney survival rate in the HU group than in the non-HU group (Figure 3(b) and Figure 3(c)).
3.4. HU as a Prognostic Indicator in CKD Patients
To determine whether severe HU (UA ≥ 7.1 mg/dL) was an independent predictor of a decline in renal function, we performed univariate and multivariate regression analyses based on the Cox hazard model for associations between the clinical findings and a ≥50% eGFR decline or ESRD during the follow-up period (Table 2). The results showed significant associations between HU [hazard ratio (HR) = 8.07, 95% CI (2.25 - 33.3), P = 0.001], diabetes mellitus (HR = 8.61, P = 0.003), and UACR (HR = 1.02, P = 0.0005) and a ≥50% eGFR decline or ESRD (Table 2(a)), and sex-stratified analysis yielded similar results. Both in male and female, HU, diabetes mellitus, UACR and eGFR were significantly associated with a ≥50% eGFR decline or ESRD (Table 2(b) and Table 2(c)).
3.5. Correlations between Serum UA Levels and Other Parameters
Since the serum UA values may have been affected by confounders, the baseline serum UA values were tested for correlations with clinical and laboratory para-
(a) (b)(c)
Figure 3. Kidney survival rate of the group with serum uric acid levels ≥7.1 mg/dL and group with serum uric acid levels <7.1 mg/dL in the total cohort (a); the males (b) and the females (c).
(a) (b) (c)
Table 2. Univariate and multivariate analysis of risk factors associated with a ≥50% eGFR decline or ESRD. (a) Total cohort n = 122; (b) male n = 48; (c) female n = 74.
Variables with P-values of less than 0.1 in the univariate model were included in the multivariate model. Abbreviation: ESRD, end-stage renal disease; n, number; CI = confidence interval; IMT, maximum cartoid intima-media thickness; eGFR, estimated glomerular filtration rate; UACR, urine albumin-to- creatinine ratio.
meters at baseline. As shown in Table 3, the UA levels were significantly negatively correlated with both eGFR (r = −0.48, P < 0.0001, Figure 4(a)) and HDL-C (r = −0.29, P = 0.001), and significantly positively correlated with SBP (r = 0.22, P = 0.01), DBP (r = 0.25, P = 0.006), MBP (r = 0.25, P = 0.007), TG (r = 0.29, P = 0.002), UACR (r = 0.25, P = 0.005), VFA (r = 0.32, P = 0.0003, Figure 4(b)), and IMT (r = 0.32, P = 0.0003, Figure 4(c)).
3.6. Comparison between the Clinical and Laboratory Findings in the Groups Established According to Serum UA Levels in the Propensity Score-Matched Cohorts
Since UA was correlated with gender, blood pressure, VFA, and especially with eGFR, we fit propensity score-matched models that included potential modifying variables (sex, eGFR, MBP, VFA) and performed subgroup analyses of the groups. We created a propensity score-matched cohort of HU patients and non-HU patients, and the comparisons between the clinical and laboratory findings of the two groups at baseline are summarized in Table 4. There were no significant differences between the propensity score-matched groups in any of the parameters except those associated with UA levels. The serum UA levels of the HU group (7.6 ± 0.4 mg/dL) were significantly higher than in the non-HU group (5.6 ± 1.5 mg/dL, P = 0.0003). The results of the Kaplan-Meier analysis with a ≥50% eGFR decline or ESRD as the end-point showed that the kidney survival rate of the HU group of CKD patients was significantly lower than in the non-HU group (Figure 5(a)). At the 100-month follow-up examination,
(a)(b)(c)
Figure 4. Relationship between serum uric acid levels and eGFR (a); visceral fat area (b), and maximum carotid intima-media thickness (c).
(a)(b)(c)
Figure 5. Kidney survival rate of the group with serum uric acid levels ≥7.1 mg/dL and group with serum uric acid levels <7.1 mg/dL in the total cohort (a); the males (b) and the females (c) of the propensity score-matched cohort.
Table 3. Coefficients of the correlation between serum uric acid levels and continuous variables.
Abbreviation: n, number; SBP, systolic blood pressure; DBP, diastolic blood pressure; MBP, mean blood pressure; BMI, body mass index; IMT, maximum cartoid intima-media thickness; eGFR, estimated glomerular filtration rate; UACR, urine albumin-to-creatinine ratio.
Table 4. Patient characteristics according to baseline serum uric acid levels. (Propensity score matched cohort n = 24).
Continuous values are expressed as means ± standard deviation. Count data are expressed as percentages. Abbreviation: n, number; HU, Hyperuricemia; UA, Uric Acid; SBP, systolic blood pressure; DBP, diastolic blood pressure; MBP, mean blood pressure; BMI, body mass index; IMT, maximum cartoid intima-media thickness; eGFR, estimated glomerular filtration rate; UACR, urine albumin-to-creatinine ratio.
there was a decrease in eGFR value at least a 50% or ESRD in 48.4% of the HU patients (log rank, P = 0.02). Although there were no statistically significant differences due to the limitation imposed by the small number of cases in the matched cohort, the sex-stratified analyses showed that the kidney survival rate of the CKD patients in the HU group was lower than in the non-HU group (Figure 5(b) and Figure 5(c)).
4. Discussion
The predictive value of UA in regard to the progression of CKD is matter of controversy [5] - [11] . Although a causal role of HU in CKD has not been established [23] [24] , the results of several epidemiologic studies have indicated the existence of an association between HU and the development of CKD [25] [26] . HU has been found to independently predict the progression of kidney disease in nephrosclerosis [4] , diabetic nephropathy [27] [28] , IgA nephropathy [29] [30] [31] [32] , chronic allograft nephropathy [33] , and CKD [10] [26] , and HU has been reported to be associated with atherosclerosis [14] [34] and hypertension [35] . Thus, HU would seem to be associated with atherosclerosis and subsequent risk of hypertension and CKD, but the serum UA elevation in CKD patients makes it difficult to assess the pathophysiological role of UA in the progression of CKD. The present study focused on atherosclerosis/arteriosclerosis- related factors and examined the predictive value of HU in relation to CKD progression by using a propensity score-matching cohort.
The results of this study showed that even when eGFR was completely adjusted by propensity score matching (Table 4), the kidney survival rate of the HU group of CKD patients was significantly lower than in the non-HU group (Figure 5(a), log rank, P = 0.02). Multivariate analysis of the Cox proportional hazards model showed that HU (UA ≥ 7.1 mg/dL) (HR= 8.07, P = 0.001) was significantly associated with a ≥50% eGFR decline or ESRD in the CKD patients (Table 2(a)). Since the BP of our cohort was well controlled (mean SBP/DBP = 125/76 mmHg) with antihypertensive agents (61.5%), hypertension was not a significant prognostic risk factor in this cohort. Rather, not only diabetes mellitus, HU was demonstrated to be a significant prognostic risk factor of CKD progression in this BP-controlled cohort. Women generally have lower UA levels, but it has never been clearly elucidated whether gender difference is involved in kidney disease progression [10] [26] [36] [37] . The results of our ROC analyses showed that the optimal UA cut-off value for discriminating ≥50% eGFR decline or ESRD in the total cohort was 7.1 mg/dL and because of the limitation of the small numbers of cases it was impossible to detect a gender difference at the cut- off value. However, our sex-stratified Kaplan-Meier survival curves and multivariate regression analyses based on the Cox hazard model showed similar significant associations between HU (UA ≥ 7.1 mg/dL) and a ≥50% eGFR decline or ESRD in both males and females (Table 2, Figure 3).
The results of the present study showed that UA was significantly correlated with VFA and with IMT (Figure 4(b) and Figure 4(c)). Recent experimental studies have indicated that UA may play a role in the pathogenesis of atherosclerosis and metabolic syndrome via the inflammatory pathway in adipose tissue [38] [39] . Several clinical studies have reported finding that a significant association between HU and VFA [40] [41] and that HU is a risk factor for carotid atherosclerosis [14] . Nakagawa et al. [42] reported finding evidence suggesting that HU might have a causal role in the development of obesity, metabolic syndrome, and endothelial dysfunction in rats fed a high-fructose, atherogenic [43] , diet. The results of our study also indicate the existence of pathophysiological relationships between UA and VFA/IMT/CKD progression the same as in these other studies. We have to pay attention to the possibility that IMT and VFA are potential confounding factors or predictive markers for progression of kidney disease. The results of our multivariate analysis of the Cox proportional hazards model adjusted for eGFR, VFA and IMT indicated that HU was independently associated with a ≥50% eGFR decline or ESRD (Table 2(a)), and the propensity score matched cohort, in which both VFA and IMT were higher in non-HU group (Table 4), showed that the kidney survival rate of the HU group of CKD patients was significantly lower than in the non-HU group (Figure 5(a)). These results indicate the most important finding that UA is an independent kidney outcome predictor. Although the mechanism of the effects of HU on atherosclerosis has not been fully elucidated, recent data have shown a direct harmful effect of UA on endothelial cells [42] [44] and smooth muscle cells [45] [46] . The precise mechanism of the nephrotoxicity of HU also remains unknown. However, HU has been found to induce systemic hypertension and afferent arteriolar sclerosis in animal models [47] [48] , and high serum UA levels have been reported to be associated with renal arteriolosclerosis in humans [4] [49] [50] . These reports support the results of the present study.
The present study had several limitations. First, it was an observational study and had the typical limitations of a small number of subjects. Second, even if we used a propensity score-matched cohort as a control group, it is difficult to select a desirable control group.
5. Conclusion
The results obtained by using a propensity score-matched cohort in the present study showed that HU is an independent predictor of CKD progression. In addition to treatment of hypertension, treatment of HU in CKD patients will become more important.
Disclosure
All authors have no conflict of interest to disclose.