# Longitudinal uric acid has nonlinear association with kidney failure and mortality in chronic kidney disease

Mathilde Prezelin-Reydit, Christian Combe, Denis Fouque, Luc Frimat, Christian Jacquelinet, Maurice Laville, Ziad Massy, Céline Lange, Carole Ayav, Roberto Pecoits-Filho, et al.

► To cite this version:

Mathilde Prezelin-Reydit, Christian Combe, Denis Fouque, Luc Frimat, Christian Jacquelinet, et al. Longitudinal uric acid has nonlinear association with kidney failure and mortality in chronic kidney disease. Scientific Reports, 2023, 13 (1), pp.3952. 10.1038/s41598-023-30902-7. hal-04032437

HAL Id: hal-04032437

https://u-picardie.hal.science/hal-04032437

Submitted on 6 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Longitudinal uric acid has nonlinear association with kidney failure and mortality in chronic kidney disease

Mathilde Prezelin-Reydit<sup>1,2,3✉</sup>, Christian Combe<sup>4,5</sup>, Denis Fouque<sup>6</sup>, Luc Frimat<sup>7,8</sup>, Christian Jacquelinet<sup>9,10</sup>, Maurice Laville<sup>11,12</sup>, Ziad A. Massy<sup>9,13</sup>, Céline Lange<sup>9</sup>, Carole Ayav<sup>14</sup>, Roberto Pecoits-Filho<sup>15,16</sup>, Sophie Liabeuf<sup>9,17,18</sup>, Bénédicte Stengel<sup>9</sup>, Jérôme Harambat<sup>1,3,19</sup>, Karen Leffondré<sup>1,3</sup> & CKD-REIN study group<sup>9*</sup>

We investigated the shape of the relationship between longitudinal uric acid (UA) and the hazard of kidney failure and death in chronic kidney disease (CKD) patients, and attempted to identify thresholds associated with increased hazards. We included CKD stage 3–5 patients from the CKD-REIN cohort with one serum UA measurement at cohort entry. We used cause-specific multivariate Cox models including a spline function of current values of UA (cUA), estimated from a separate linear mixed model. We followed 2781 patients (66% men, median age, 69 years) for a median of 3.2 years with a median of five longitudinal UA measures per patient. The hazard of kidney failure increased with increasing cUA, with a plateau between 6 and 10 mg/dl and a sharp increase above 11 mg/dl. The hazard of death had a U-shape relationship with cUA, with a hazard twice higher for 3 or 11 mg/dl, compared to 5 mg/dl. In CKD patients, our results indicate that UA above 10 mg/dl is a strong risk marker for kidney failure and death and that low UA levels below 5 mg/dl are associated with death before kidney failure.

Chronic kidney disease (CKD) is recognized as a major public health problem and identification of modifiable determinants of CKD progression, other than hypertension or proteinuria<fnref n="1" />, is essential to develop effective strategies to slow disease progression. In the past two decades, uric acid (UA) has drawn attention in the nephrology community.

Previous cohort studies found conflicting results on the relationship between UA and CKD progression<sup>2–11</sup>, and did not well characterize the shape of this longitudinal relationship<sup>12</sup>. The discrepancy between results may be due to differences in CKD stages at baseline, duration of follow-up, definitions of CKD progression,

<footnote n="1">¹INSERM, Bordeaux Population Health Research Center, UMR1219, Univ Bordeaux, Bordeaux, France. ²Maison du REIN AURAD Aquitaine, 2 allée des demoiselles, 33170 Gradignan, France. ³INSERM, CIC1401-EC, Univ Bordeaux, Bordeaux, France. ⁴Department of Nephrology Dialysis Transplantation, Centre Hospitalier Universitaire de Bordeaux, Bordeaux, France. ⁵INSERM U1026, Univ Bordeaux, Bordeaux, France. ⁶Nephrology Department, Centre Hospitalier Lyon Sud, Université de Lyon, Carmen, Pierre-Bénite, France. ⁷Nephrology Department, Centre Hospitalier Régional Universitaire de Nancy, Vandoeuvre-lès-Nancy, France. ⁸APEMAC, Lorraine University, Vandoeuvre-lès-Nancy, France. ⁹Center for Research in Epidemiology and Population Health (CESP), Clinical Epidemiology Team, Paris-Saclay University, Paris-Sud University, Versailles Saint-Quentin-en-Yvelines University, Villejuif, France. ¹⁰Renal Epidemiology and Information Network Registry, Biomedicine Agency, Saint Denis, France. ¹¹Carmen INSERM U1060, Université Claude Bernard Lyon 1, Pierre-Bénite, France. ¹²AURAL, Lyon, France. ¹³Division of Nephrology, Ambroise Paré University Hospital, Assistance publique - Hôpitaux de Paris, Boulogne-Billancourt/Paris, France. ¹⁴CHRU-Nancy, INSERM, Université de Lorraine, CIC, Epidémiologie Clinique, 54000 Nancy, France. ¹⁵Arbor Research Collaborative for Health, Ann Arbor, MI, USA. ¹⁶School of Medicine, Pontifícia Universidade Catolica do Parana, Curitiba, Brazil. ¹⁷Department of Clinical Pharmacology, Centre Hospitalier Universitaire, Amiens, France. ¹⁸Laboratoire MP3CV, EA7517, Université de Picardie Jules Verne, 80000 Amiens, France. ¹⁹Pediatric Nephrology Unit, Pellegrin-Enfants Hospital, Centre Hospitalier Universitaire de Bordeaux, Centre de Référence Maladies rénales rares Sorare, Bordeaux, France. ²A list of authors and their affiliations appears at the end of the paper.✉email: mathilde.reydit@maisondurein-aurad.fr</footnote>
the choice of adjustment factors¹³, or to the use of a single measure of UA assessed at baseline, i.e. entry into the cohort²⁻⁵,⁸⁻¹⁰. Indeed, entry into a cohort does not usually correspond to any relevant time point in patient's course of CKD and varies across studies and from patient to patient. To estimate the association between UA and CKD progression, most previous studies thus compared the baseline uric acid values between two patients who potentially did not enter the cohort at the same moment of their CKD history, and fully ignored their respective subsequent evolution of uric acid between baseline and the current time when the risk was assessed. Yet, UA may change over time, and from a clinical point of view, it might be more relevant to compare the level of UA reached at the time when the risk is assessed. To our knowledge, only one single-center cohort study in Taiwan used longitudinal measures of UA¹⁴. They found that elevated UA trajectories had an increased hazard of dialysis initiation and death. However, the analysis based on classes of trajectories did not allow the examination of the shape of the relationship between the level of UA reached at a given time and the hazard of kidney failure or death at the same time, and thus the identification of potential critical thresholds of UA, which yet might be useful in monitoring patients with CKD.

The objective of the present study was to investigate, in the French Chronic Kidney Disease-Renal Epidemiology and Information Network (CKD-REIN) cohort¹⁵, the shape of the relationship between longitudinal UA and the hazard of kidney failure and death in CKD stage 3–5 patients, and identify potential thresholds of longitudinal UA associated with increased hazards.

## Methods

**Study population.** CKD-REIN is an ongoing prospective cohort of CKD Stage 3–5 patients receiving nephrologist-led care, without prior chronic dialysis or kidney transplantation¹⁵. The study included 3033 patients over 18 years of age. They were first selected and recruited in an enrolment phase in 40 nephrology centers located over all metropolitan France and representative of all centers with respect to legal status (public, private non-for-profit, and private for-profit), and then actively included into the cohort at a first visit (baseline) between 2013 and 2016, in the same centers, after obtaining inform consent¹⁵. The study protocol was conducted with adherence to the Declaration of Helsinki and approved by the institutional review board at the French National Institute of Health and Medical Research (INSERM; reference: IRB00003888). The study was registered at ClinicalTrials.gov (NCT03381950). In the present study, we included patients who had at least one serum UA and one creatinine measurement within 6 months of their inclusion into the cohort.

**Exposure and outcomes.** The exposure was serum UA concentration, assessed in each center along with routine laboratory investigations, at baseline, annually per protocol and more frequently if considered necessary by the nephrologist.

The outcomes of interest were (1) kidney failure assessed by initiation of chronic dialysis or pre-emptive transplantation, and (2) death before kidney failure. To ensure complete collection of kidney failure, a record linkage was performed with the national REIN registry. Administrative censoring was performed on July 30, 2018.

**Covariates.** Baseline characteristics were recorded by clinical research associates from medical files or by interview. Data included age, sex, body mass index (BMI), hypertension (patients having an office blood pressure greater than or equal to 140/90 mmHg or an antihypertensive treatment), cardiovascular history (coronary artery disease, arrhythmic disorders, congestive heart failure, stroke, peripheral vascular disease and/or valvulopathy), diabetes (diabetes history or antidiabetic treatment or glycated hemoglobin ≥ 6.5% or fasting glycemia ≥ 7 mmol/l or non-fasting glycemia ≥ 11), gout history, dyslipidemia, primary kidney disease, time since CKD diagnosis (time elapsed from the date of CKD diagnosis found in the medical record and the cohort entry), number of consultation in the previous year with nephrologist and dietician, treatment (urate-lowering therapy (ULT), diuretics, antiplatelet agents, renin-angiotensin system inhibitors (RASi)), laboratory data (serum creatinine, eGFR estimated by the CKD-EPi equation, serum UA, albuminemia, C-reactive protein and, albuminuria—or equivalent—classified according to the KDIGO 2012 guidelines¹⁶), salt intake (estimated by 24-h natriuresis) and protein intake (estimated by 24-h urinary urea)¹⁷, medication adherence according to the Girerd score in categories (good, minimal and poor)¹⁸, health literacy according to their need for help reading medical documents (never need vs always or partly need)¹⁹ and type of center (university, non-university hospital, private non-profit and private for-profit clinic).

**Statistical analyses.** We described characteristics of all the patients in the CKD-REIN cohort, those included in the present study, as well as those who contributed to the estimation of the different statistical models. We also estimated the crude association between patients' characteristics and UA at baseline.

To estimate the shape of the relationship between the level of UA reached at a given time, i.e. current level of UA (cUA), and the hazard of kidney failure or death at the same time, we used a two-stage statistical approach which have been evaluated previously²⁰⁻²². At Stage 1, we estimated for each patient its full longitudinal trajectory of UA from a linear mixed model estimated using information on all patients. This allowed us to get cUA at any time point of follow-up for each patient. At Stage 2, we estimated the association between cUA and the hazards of kidney failure or death using cause-specific Cox models including cUA as a continuous time-dependent variable. More specifically, the linear mixed model (at Stage 1) included a 3-knot natural cubic spline function of time with random effects on each coefficient, and some selected baseline factors (age, sex, hypertension, eGFR, BMI, use of ULT and diuretics) to make the estimation of all individual UA trajectories more accurate. The cause-specific Cox models (at Stage 2) included a 2-knot natural cubic spline function of cUA (derived from Stage 1) to allow the estimation of nonlinear association and thus the detection of potential thresholds. The number of knots were selected using the Akaike Information Criterion (AIC)²³. The Cox models were adjusted for a set of
confounders which were selected from a directed acyclic graph (DAG), i.e. a diagram of causal pathways sum-
marizing a priori hypothetical causal relationship between variables (Figs. S1 and S2). This approach allows the
selection of optimal set of adjustment factors (i.e. factors associated both with the exposure and the outcome),
avoiding adjustment for unnecessary factors, mediators, and colliders²⁴. For kidney failure, this included age (in
years), sex, primary kidney disease (diabetic, glomerular, hypertensive, vascular, tubulo-interstitial, polycystic or
unknown nephropathy), hypertension (yes/no), diabetes (yes/no), cardiovascular disease (yes/no), dyslipidemia
(yes/no), BMI (< 25, 25–30, ≥ 30 kg/m²), albuminuria (A1, A2, A3), CKD stages (5, 4, 3B, 3A or less), medica-
tion adherence (good, minimal, poor), RASi (yes/no), and ULT (yes/no) (Fig. S1, Model 1 for kidney failure).
For death, we added spironolactone (yes/no) and anti-platelet agents (yes/no) (Fig. S2, Model 1 for death). All
confounders (including renal function represented by CKD stages) were taken at baseline only to respect the
temporal sequence between confounders and the exposure, i.e. the subsequent cUA which was the only time-
dependent variable. The analyses thus accounted for the fact that UA may be consequence of a decreased renal
function by adjusting the effect of cUA for CKD stage at baseline.

To investigate if the association between cUA and the hazard of kidney failure or death differed according to
sex, we included interaction terms with the spline function of cUA, and test the interaction using the likelihood
ratio test.

In a first sensitivity analysis, we further adjusted Model 1 for salt intake (< 95, 95–127, 128–170 and ≥ 170 mmol/day) and protein intake (in mmol/day) at baseline (called “eating habits” in the DAG) (Model 2) because patients with high salt or high protein intake may be at higher risk of hyperuricemia²⁵,²⁶ and CKD progression²⁷,²⁸. We perform this further adjustment in the subsample of patients having 24-h natriuresis and urinary urea at baseline. In a second sensitivity analysis, for comparison with previous studies, we estimated the association between baseline value of UA and outcomes using the same set of adjustment factors as Model 1 (Model 3).

In all Cox models, we accounted for correlation between patients of the same type of center using robust standard errors based on the sandwich estimator²⁹. Proportional hazards assumption was checked using Schoenfeld residuals. Linearity of the effect of all adjusting quantitative variables was checked using 4-df penalized spline functions³⁰, which was kept in the model if the effect was nonlinear³¹. All analyses were performed using R version 3.6.0³².

Results

**Patients' selection and characteristics, UA distribution, and number of events.** Among the 3033 patients enrolled in the CKD-REIN study, 2781 patients had a UA measurement within six months of inclusion into the cohort and were thus included in the present study (Fig. 1). At baseline, they had a median age of 69.0 years (interquartile range (IQR): 60.0–76.0), 65.5% were men, 96.1% had hypertension, 42.5% diabetes, 53.4% cardiovascular disease, 73.2% dyslipidemia, and 22% a gout history (Table 1, Included Population). About 40% of patients had either hypertensive nephropathy or diabetic nephropathy. Median eGFR was 32 ml/min/1.73 m² (IQR: 23–41), with 94% of CKD stage 3 or 4 (6% of patients who had progressed to CKD stage 2 or

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Included population, CKD-REIN, France, 2013–2018.</figcaption></figure>
<table><thead><tr><th rowspan="2"></th><th colspan="2">CKD-REIN population (n=3033)</th><th colspan="2">Included population (n=2781)</th><th colspan="2">Population for Models 1 and 3 (n=2344)</th><th colspan="2">Population for Model 2 (n=1212)</th></tr><tr><th>N</th><th>n (%) or median (IQR)</th><th>N</th><th>n (%) or median (IQR)</th><th>N</th><th>n (%) or median (IQR)</th><th>N</th><th>n (%) or median (IQR)</th></tr></thead><tbody><tr><td>Age (years)</td><td>3033</td><td>69.0 (60.0-76.0)</td><td>2781</td><td>69.0 (60.0, 76.0)</td><td>2344</td><td>68.0 (60.0-76.0)</td><td>1212</td><td>68.5 (61.0-76.0)</td></tr><tr><td>Male gender</td><td>3033</td><td>1982 (65.3)</td><td>2781</td><td>1821 (65.5)</td><td>2344</td><td>1549 (66.1)</td><td>1212</td><td>817 (67.4)</td></tr><tr><td>Body mass index (kg/m²)</td><td>2968</td><td>27.8 (24.6-31.8)</td><td>2726</td><td>27.8 (24.6, 31.6)</td><td>2344</td><td>27.9 (24.7-31.8)</td><td>1212</td><td>28.2 (25.0-32.1)</td></tr><tr><td>Hypertension</td><td>3026</td><td>2915 (96.3)</td><td>2774</td><td>2666 (96.1)</td><td>2344</td><td>2260 (96.4)</td><td>1212</td><td>1168 (96.4)</td></tr><tr><td>Cardiovascular history</td><td>2991</td><td>1594 (53.3)</td><td>2740</td><td>1464 (53.4)</td><td>2344</td><td>1230 (52.5)</td><td>1212</td><td>640 (52.8)</td></tr><tr><td>Diabetes</td><td>3026</td><td>1301 (43.0)</td><td>2778</td><td>1180 (42.5)</td><td>2344</td><td>1025 (43.7)</td><td>1212</td><td>520 (42.9)</td></tr><tr><td>Dyslipidemia</td><td>3019</td><td>2223 (73.6)</td><td>2719</td><td>2025 (73.2)</td><td>2344</td><td>1759 (75.0)</td><td>1212</td><td>920 (75.9)</td></tr><tr><td>Gout history</td><td>2968</td><td>618 (20.8)</td><td>2768</td><td>597 (22.0)</td><td>2325</td><td>525 (22.6)</td><td>1205</td><td>283 (23.5)</td></tr><tr><td>Primary kidney disease</td><td>3033</td><td></td><td>2781</td><td></td><td>2344</td><td></td><td>1212</td><td></td></tr><tr><td>Diabetic nephropathy</td><td></td><td>611 (20.1)</td><td></td><td>545 (19.6)</td><td></td><td>475 (20.3)</td><td></td><td>219 (18.1)</td></tr><tr><td>Glomerulopathy</td><td></td><td>532 (17.5)</td><td></td><td>505 (18.2)</td><td></td><td>455 (19.4)</td><td></td><td>257 (21.2)</td></tr><tr><td>Hypertensive nephropathy</td><td></td><td>633 (20.9)</td><td></td><td>570 (20.5)</td><td></td><td>476 (20.3)</td><td></td><td>239 (19.7)</td></tr><tr><td>Vascular nephropathy</td><td></td><td>216 (7.1)</td><td></td><td>203 (7.3)</td><td></td><td>156 (6.7)</td><td></td><td>93 (7.7)</td></tr><tr><td>Tubulo-interstitial nephropathy</td><td></td><td>377 (12.4)</td><td></td><td>349 (12.5)</td><td></td><td>297 (12.7)</td><td></td><td>159 (13.1)</td></tr><tr><td>Polykystic renal disease</td><td></td><td>166 (5.5)</td><td></td><td>157 (5.6)</td><td></td><td>131 (5.6)</td><td></td><td>74 (6.1)</td></tr><tr><td>Other or unknown</td><td></td><td>498 (16.4)</td><td></td><td>452 (16.2)</td><td></td><td>354 (15.1)</td><td></td><td>171 (14.1)</td></tr><tr><td>Glomerular filtration rate (ml/min/1.73 m²)</td><td>3027</td><td>32.0 (23.2,41.4)</td><td>2781</td><td>31.8 (23.2, 41.4)</td><td>2344</td><td>31.9 (23.3-41.5)</td><td>1212</td><td>31.5 (23.2-40.8)</td></tr><tr><td>Chronic Kidney Disease stage</td><td>3027</td><td></td><td>2781</td><td></td><td>2344</td><td></td><td></td><td></td></tr><tr><td>1</td><td></td><td>0 (0.0)</td><td></td><td>0 (0.0)</td><td></td><td>0 (0.0)</td><td></td><td>0 (0.0)</td></tr><tr><td>2</td><td></td><td>65 (2.1)</td><td></td><td>58 (2.1)</td><td></td><td>46 (2.0)</td><td></td><td>25 (2.1)</td></tr><tr><td>3</td><td></td><td>1612 (53.3)</td><td></td><td>1464 (52.6)</td><td></td><td>1234 (52.6)</td><td></td><td>619 (51.1)</td></tr><tr><td>4</td><td></td><td>1233 (40.7)</td><td></td><td>1146 (41.2)</td><td></td><td>977 (41.7)</td><td></td><td>522 (43.1)</td></tr><tr><td>5</td><td></td><td>117 (3.9)</td><td></td><td>113 (4.0)</td><td></td><td>87 (3.7)</td><td></td><td>46 (3.8)</td></tr><tr><td>Time since Chronic Kidney Disease diagnosis (years)</td><td>2891</td><td>5.1 (2.4, 8.0)</td><td>2655</td><td>5.2 (2.5, 10.3)</td><td>2266</td><td>5.2 (2.5-10.4)</td><td>1168</td><td>5.4 (2.6-11.5)</td></tr><tr><td>Uric acid (mg/dl)</td><td>2725</td><td>7.1 (5.8, 8.4)</td><td>2722</td><td>7.1 (5.8, 8.4)</td><td>2310</td><td>7.2 (5.9-8.4)</td><td>1204</td><td>7.2 (5.9-8.5)</td></tr><tr><td>Uric acid (in categories)</td><td>2725</td><td></td><td>2722</td><td></td><td>2310</td><td></td><td>1204</td><td></td></tr><tr><td>&lt;4 mg/dl</td><td></td><td>119 (4.4)</td><td></td><td>119 (4.4)</td><td></td><td>100 (4.3)</td><td></td><td>43 (3.6)</td></tr><tr><td>4-5 mg/dl</td><td></td><td>224 (8.2)</td><td></td><td>223 (8.2)</td><td></td><td>185 (8.0)</td><td></td><td>105 (8.7)</td></tr><tr><td>5-6 mg/dl</td><td></td><td>416 (15.3)</td><td></td><td>416 (15.3)</td><td></td><td>354 (15.3)</td><td></td><td>173 (14.4)</td></tr><tr><td>6-8 mg/dl</td><td></td><td>1048 (38.5)</td><td></td><td>1048 (38.5)</td><td></td><td>891 (38.6)</td><td></td><td>455 (37.8)</td></tr><tr><td>8-10 mg/dl</td><td></td><td>669 (24.6)</td><td></td><td>668 (24.5)</td><td></td><td>569 (24.6)</td><td></td><td>319 (26.5)</td></tr><tr><td>&gt; 10 mg/dl</td><td></td><td>249 (9.1)</td><td></td><td>248 (9.1)</td><td></td><td>211 (9.1)</td><td></td><td>109 (9.1)</td></tr><tr><td>Protein-to-creatinine ratio (mg/mmol)</td><td>1849</td><td>35.7 (13.7, 115.7)</td><td>1724</td><td>35.7 (13.5, 115.7)</td><td>1611</td><td>35.8 (13.1-114.4)</td><td>743</td><td>36.2 (13.5-126.7)</td></tr><tr><td>Albumin-to-creatinine ratio</td><td>2693</td><td></td><td>2507</td><td></td><td>2344</td><td></td><td></td><td></td></tr><tr><td>&lt;3 mg/mmol</td><td></td><td>742 (27.6)</td><td></td><td>684 (27.3)</td><td></td><td>642 (27.4)</td><td></td><td>318 (26.2)</td></tr><tr><td>3-30 mg/mmol</td><td></td><td>847 (31.5)</td><td></td><td>788 (31.4)</td><td></td><td>729 (31.1)</td><td></td><td>391 (32.3)</td></tr><tr><td>&gt; 30 mg/mmol</td><td></td><td>1104 (41.0)</td><td></td><td>1035 (41.3)</td><td></td><td>973 (41.5)</td><td></td><td>503 (41.5)</td></tr><tr><td>Natriuresis (mmol/day)</td><td>1663</td><td>128.0 (95.0-170.5)</td><td>1602</td><td>128.0 (95.0-170.0)</td><td>1458</td><td>128.0 (96.0-169.0)</td><td>1212</td><td>128.0 (96.0-167.2)</td></tr><tr><td>Urinary urea (mmol/day)</td><td>1403</td><td>304.7 (233.2-389.5)</td><td>1361</td><td>305.0 (234.8-389.6)</td><td>1233</td><td>306.6 (236.4-389.6)</td><td>1212</td><td>306.2 (236.4-388.1)</td></tr><tr><td>Albumin (µmol/l)</td><td>2459</td><td>585.5 (550.7-623.2)</td><td>2354</td><td>585.5 (550.7-623.2)</td><td>2058</td><td>584.1 (550.7-623.2)</td><td>1119</td><td>588.4 (550.7-623.2)</td></tr><tr><td>C-reactive Protein (mg/l)</td><td>1197</td><td>3.8 (1.7-7.3)</td><td>1149</td><td>3.9 (1.7-7.3)</td><td>1007</td><td>3.9 (1.8-7.4)</td><td>520</td><td>3.9 (1.7-7.3)</td></tr><tr><td>Urate lowering therapy (febuxostat or allopurinol)</td><td>3024</td><td>999 (33.0)</td><td>2774</td><td>936 (33.7)</td><td>2344</td><td>831 (35.5)</td><td></td><td>424 (35.0)</td></tr><tr><td>Diuretics (all types)</td><td>3024</td><td>1605 (53.1)</td><td>2774</td><td>1455 (52.5)</td><td>2344</td><td>1240 (52.9)</td><td>1212</td><td>618 (51.0)</td></tr><tr><td>Spironolactone</td><td>3024</td><td>109 (3.6)</td><td>2774</td><td>98 (3.5)</td><td>2344</td><td>87 (3.7)</td><td>1212</td><td>43 (3.5)</td></tr><tr><td>Antiplatelet agents</td><td>3024</td><td>1238 (40.7)</td><td>2774</td><td>1133 (40.8)</td><td>2344</td><td>978 (41.7)</td><td>1212</td><td>493 (40.5)</td></tr><tr><td>Renin-angiotensin inhibitors</td><td>3024</td><td>2294 (75.9)</td><td>2774</td><td>2091 (75.4)</td><td>2344</td><td>1808 (77.1)</td><td>1212</td><td>956 (78.9)</td></tr><tr><td>Medication adherence according to the Gیرerd score</td><td>3002</td><td></td><td>2753</td><td></td><td>2344</td><td></td><td>1212</td><td></td></tr><tr><td>Good (score equal to 0)</td><td></td><td>1129 (37.6)</td><td></td><td>1028 (37.3)</td><td></td><td>875 (37.3)</td><td></td><td>418 (34.5)</td></tr><tr><td>Minimal (score equal to 1 or 2)</td><td></td><td>1651 (55.0)</td><td></td><td>1528 (55.5)</td><td></td><td>1313 (56.0)</td><td></td><td>723 (59.7)</td></tr><tr><td>Poor (score ≥ 3)</td><td></td><td>222 (7.4)</td><td></td><td>197 (7.2)</td><td></td><td>156 (6.7)</td><td></td><td>71 (5.8)</td></tr><tr><td>Health literacy according to their need for help reading medical documents</td><td>3033</td><td></td><td>2781</td><td></td><td>2344</td><td></td><td>1212</td><td></td></tr></tbody><tfoot><tr><td colspan="9">Continued</td></tr></tfoot></table>
<table><caption>**Table 1.** Characteristics of included population compared to the characteristics of CKD-REIN population. N: available data. IQR: InterQuartile Range 25–75. Uric acid in mg/dl to µmol/l, × 59.48. Cardiovascular history defined as patients having coronary artery disease, arrhythmic disorders, congestive heart failure, stroke, peripheral vascular disease and/or valvulopathy. Diabetes defined as patients having diabetes history or antidiabetic treatment or glycated hemoglobin ≥ 6.5% or fasting glycemia ≥ 7 mmol/l or non-fasting glycemia ≥ 11 mmol/l. Hypertension defined as patients having an office blood pressure greater than or equal to 140/90 mmHg or an antihypertensive treatment.</caption><thead><tr><th rowspan="2"></th><th colspan="2">CKD-REIN population (n=3033)</th><th colspan="2">Included population (n=2781)</th><th colspan="2">Population for Models 1 and 3 (n=2344)</th><th colspan="2">Population for Model 2 (n=1212)</th></tr><tr><th>N</th><th>n (%) or median (IQR)</th><th>N</th><th>n (%) or median (IQR)</th><th>N</th><th>n (%) or median (IQR)</th><th>N</th><th>n (%) or median (IQR)</th></tr></thead><tbody><tr><td>Never</td><td></td><td>575 (19.0)</td><td></td><td>519 (18.7)</td><td></td><td>420 (17.9)</td><td></td><td>219 (18.1)</td></tr><tr><td>Rarely, sometimes, often, or always</td><td></td><td>2458 (81.0)</td><td></td><td>2262 (81.3)</td><td></td><td>1924 (82.1)</td><td></td><td>993 (81.9)</td></tr><tr><td>Number of nephrological consultations in the year before inclusion</td><td>2612</td><td></td><td>2414</td><td></td><td>2056</td><td></td><td>1080</td><td></td></tr><tr><td>0</td><td></td><td>52 (2.0)</td><td></td><td>48 (2.0)</td><td></td><td>41 (2.0)</td><td></td><td>20 (1.9)</td></tr><tr><td>1 or 2</td><td></td><td>1660 (63.6)</td><td></td><td>1516 (62.8)</td><td></td><td>1277 (62.1)</td><td></td><td>676 (62.6)</td></tr><tr><td>3</td><td></td><td>624 (23.9)</td><td></td><td>587 (24.3)</td><td></td><td>519 (25.2)</td><td></td><td>269 (24.9)</td></tr><tr><td>4 or more</td><td></td><td>276 (10.6)</td><td></td><td>263 (10.9)</td><td></td><td>219 (10.7)</td><td></td><td>115 (10.6)</td></tr><tr><td>Number of dietary consultations in the year before inclusion</td><td>2476</td><td></td><td>2295</td><td></td><td>1956</td><td></td><td>1017</td><td></td></tr><tr><td>0</td><td></td><td>1862 (75.2)</td><td></td><td>1723 (75.1)</td><td></td><td>1449 (74.1)</td><td></td><td>730 (71.8)</td></tr><tr><td>1</td><td></td><td>424 (17.1)</td><td></td><td>394 (17.1)</td><td></td><td>348 (17.8)</td><td></td><td>200 (19.7)</td></tr><tr><td>2</td><td></td><td>112 (4.5)</td><td></td><td>105 (4.6)</td><td></td><td>91 (4.7)</td><td></td><td>50 (4.9)</td></tr><tr><td>3 or more</td><td></td><td>78 (3.2)</td><td></td><td>73 (3.2)</td><td></td><td>68 (3.5)</td><td></td><td>37 (3.6)</td></tr><tr><td>Type of center</td><td>2892</td><td></td><td>2707</td><td></td><td>2344</td><td></td><td>1212</td><td></td></tr><tr><td>University center</td><td></td><td>1734 (60.0)</td><td></td><td>1609 (59.4)</td><td></td><td>1395 (59.5)</td><td></td><td>764 (63.0)</td></tr><tr><td>Hospital center</td><td></td><td>577 (20.0)</td><td></td><td>554 (20.5)</td><td></td><td>482 (20.6)</td><td></td><td>260 (21.5)</td></tr><tr><td>Non profit institution</td><td></td><td>119 (4.1)</td><td></td><td>115 (4.2)</td><td></td><td>113 (4.7)</td><td></td><td>36 (3.0)</td></tr><tr><td>For-profit institution</td><td></td><td>462 (15.9)</td><td></td><td>429 (15.8)</td><td></td><td>354 (15.1)</td><td></td><td>152 (12.5)</td></tr></tbody></table>

5 between their selection and their actual inclusion into the cohort). Median salt intake was 6 g/day (natriuresis of 128 mmol/day) and estimated median protein intake was 61.1 g/day (urinary urea of 305 mmol/day). A total of 938 patients (33.7%) were prescribed ULT (Table 1, Included Population). Median UA at baseline was 7.1 mg/dl (IQR: 5.8–8.5), and was statistically significantly higher in younger men and in patients with higher BMI, advanced CKD stages, diabetes, cardiovascular history, no gout history, and in those receiving diuretics and no ULT (Table S1).

Over the follow-up (Median 3.2, IQR: 2.6–3.8 years), the 2781 included patients had a median of 5 UA measures, with a median of 120 days between two consecutive measures (Table 2, Fig. S3A). This led to a total of 16 947 measures of UA which were mostly observed during the first three years of follow-up (Figure S3B) and were normally distributed in both men and women (Figure S4). Individual observed values of UA, as well as the true UA trajectory estimated from the linear mixed model, are shown in Figure S5 for some selected patients with extreme UA values. The overall goodness of fit of the linear mixed model is described in Figure S6, which shows how on average the individual predicted values of UA were close to the observed values all along the follow-up. Of the 2781 patients, 439 (15.8%) initiated dialysis (n = 375) or received a pre-emptive transplant (n = 64) and 264 (9.5%) died before kidney failure during the follow-up (Table 2).

Compared to the 3033 patients participating in the CKD-REIN cohort, the 2781 included patients (with UA measurement within six months of inclusion into the cohort) had similar baseline characteristics (Table 1). The 2344 patients used for Models 1 and 3 (with no missing data on adjustment factors) had a higher proportion of ULT use at baseline (Table 1). The 1212 patients used for Model 2 (with further 24-h urine collection at baseline) were more often followed-up in a University Hospital and tended to have more dietary consultations before inclusion (Table 1).

**UA and risk of kidney failure.** The hazard of kidney failure increased with increasing cUA, with a plateau for cUA between 6 and 10 mg/dl (Fig. 2A). At any time after inclusion into the cohort, patients with a cUA of precisely 3 mg/dl, had a 59% decreased hazard of kidney failure at that time compared to patients with a cUA of precisely 5 mg/dl at the same time (HR 0.41, 95% confidence interval (CI): 0.31, 0.54, Model 1 in Table 3). The hazard of kidney failure was increased by 70% for patients with a cUA of precisely 11 mg/dl compared to patients with a cUA of precisely 5 mg/dl at the same time (HR 1.70, 95% CI: 1.18, 2.47, Model 1 in Table 3). The association between cUA and the hazard of kidney failure tended to be similar in men and women (p-value for interaction of 0.07). However, an increase in cUA from 3 to 7 mg/dl was associated with a moderately higher increase in the hazard of kidney failure in females than in males. Above 10 mg/dl, the increase in the hazard of kidney failure was detectable in men only (Figure S7), the data being too sparse in women over that range (Figure S4). Further adjustment for salt and protein intake at baseline weakly affected the association in the subsample of the
<table><caption>**Table 2.** Description of repeated uric acid measures and outcomes over the follow-up in the population included in our analysis (N = 2781). UA uric acid, IQR interquartile range.</caption><thead><tr><th></th><th>Median (IQR 25-75)</th><th>n (%)</th></tr></thead><tbody><tr><td>Total number of UA measures over follow-up</td><td></td><td>16,947</td></tr><tr><td>Number of UA measures by patient</td><td>5 (3-8)</td><td></td></tr><tr><td>Time interval between two consecutive UA measures (days)</td><td>120 (63-197)</td><td></td></tr><tr><td>Number of patients with</td><td></td><td></td></tr><tr><td>Only one UA measure</td><td></td><td>231 (8.3)</td></tr><tr><td>Two UA measures</td><td></td><td>272 (9.8)</td></tr><tr><td>Three UA measures</td><td></td><td>323 (11.6)</td></tr><tr><td>More than three UA measures</td><td></td><td>1955 (70.5)</td></tr><tr><td>Number of UA measures over follow-up in each of the following UA value classes</td><td></td><td></td></tr><tr><td>&lt;4 mg/dl</td><td></td><td>808 (4.8)</td></tr><tr><td>4-6 mg/dl</td><td></td><td>4203 (24.8)</td></tr><tr><td>6-8 mg/dl</td><td></td><td>6509 (38.4)</td></tr><tr><td>8-10 mg/dl</td><td></td><td>4022 (23.7)</td></tr><tr><td>&gt;10 mg/dl</td><td></td><td>1405 (8.3)</td></tr><tr><td>Number of patients with</td><td></td><td></td></tr><tr><td>Kidney failure</td><td></td><td>439 (15.8%)</td></tr><tr><td>Death before kidney failure</td><td></td><td>264 (9.5%)</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** (A) Estimated effect of *current* uric acid value on the hazard of kidney failure in all patients (n = 2344 including 382 KRT, Model 1 in Table 3). (B) Estimated effect of *baseline* uric acid value on the hazard of kidney failure in all patients (n = 2344, Model 3 in Table 3). Results from cause-specific Cox models using a spline function for uric acid, adjusted for age, sex, primary kidney disease, hypertension, diabetes, cardiovascular disease, dyslipidemia, body mass index, albuminuria, medication adherence, use of renin-angiotensin system inhibitors and urate lowering therapy, all at baseline. The reference value of uric acid for the HR indicated in the y-axis was arbitrarily chosen at 5 mg/dl, which corresponds to the midpoint of the normal range of uric acid (uric acid in mg/dl to µmol/l: ×59.48). CKD-REIN cohort, France, 2013–2018.</figcaption></figure>

1212 patients with available information (Model 2 vs. Model 1 in Table 3) but adjusting or not for eating habits in this subsample produced very similar results (Table S2). Finally, the association was much weaker with baseline UA than with cUA (Model 3 in Table 3, Fig. 2B).

**UA and risk of death before kidney failure.** The hazard of death before kidney failure had a U-shape relationship with cUA (Fig. 3A), with the lowest mortality for a cUA of 6 mg/dl. At any time after inclusion into
<table><thead><tr><th rowspan="2"></th><th rowspan="2">Current value of UA*</th><th colspan="2">Model 1<br>(N=2344)</th><th colspan="2">Model 2<br>(N=1212)</th><th colspan="2">Model 3<br>(N=2344)</th></tr><tr><th>HR</th><th>95% CI</th><th>HR</th><th>95% CI</th><th>HR</th><th>95% CI</th></tr></thead><tbody><tr><td rowspan="5">Kidney failure</td><td>3 mg/dl</td><td>0.41</td><td>0.31-0.54</td><td>0.36</td><td>0.18-0.69</td><td>0.98</td><td>0.77-1.25</td></tr><tr><td>5 mg/dl</td><td>1</td><td></td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>7 mg/dl</td><td>1.28</td><td>1.15-1.43</td><td>1.16</td><td>0.78-1.72</td><td>1.03</td><td>0.97-1.09</td></tr><tr><td>9 mg/dl</td><td>1.23</td><td>1.03-1.46</td><td>0.98</td><td>0.64-1.51</td><td>1.06</td><td>0.88-1.28</td></tr><tr><td>11 mg/dl</td><td>1.70</td><td>1.18-2.47</td><td>1.64</td><td>1.20-2.24</td><td>1.09</td><td>0.83-1.42</td></tr><tr><td rowspan="5">Death</td><td>3 mg/dl</td><td>1.80</td><td>1.47-2.19</td><td>1.43</td><td>0.93-2.20</td><td>0.92</td><td>0.62-1.38</td></tr><tr><td>5 mg/dl</td><td>1</td><td></td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>7 mg/dl</td><td>0.88</td><td>0.75-1.04</td><td>0.92</td><td>0.89-0.95</td><td>1.11</td><td>1.00-1.23</td></tr><tr><td>9 mg/dl</td><td>1.25</td><td>0.90-1.76</td><td>1.23</td><td>0.95-1.60</td><td>1.22</td><td>0.94-1.57</td></tr><tr><td>11 mg/dl</td><td>2.30</td><td>1.84-2.87</td><td>2.23</td><td>1.73-2.89</td><td>1.23</td><td>0.97-1.56</td></tr></tbody></table>

**Table 3.** Association between current or baseline value of uric acid and the hazard of kidney failure or death before kidney failure. Results from time-dependent cause-specific Cox models accounting for nonlinear effect uric acid. CKD-REIN cohort, France, 2013–2018. UA, uric acid; HR, hazard ratio; CI, confidence intervals. Uric acid in mg/dl to µmol/l, × 59.48. Model 1: Cox model with UA as a continuous time-dependent variable and adjusted for age, sex, CKD stage, primary kidney disease, hypertension, diabetes, cardiovascular disease, dyslipidemia, body mass index, albuminuria, medication adherence, use of renin-angiotensin system inhibitors and urate lowering therapy, all at baseline. HR of death were further adjusted for spironolactone and antiplatelet agents at baseline. Model 2: Model 1 further adjusted for salt and protein intake at baseline. Model 3: Cox model with UA as a continuous variable measured only at baseline and adjusted for the same factors as Model 1. *The listed values of UA are precise current values since uric acid was taken as a continuous time-dependent covariate in the Cox model. HR of 1.70 for example means that a patient with a current value of uric acid of precisely 11 mg/dl had a 70% increased hazard of kidney failure at that time of follow-up compared to a patient with a value of uric acid of precisely 5 mg/dl at the same time.

<figure><img src="image_5.png" /></figure>

**Figure 3.** Estimated effect of current uric acid value on the hazard of death before kidney failure, Panel (A) Estimated effect of current uric acid value on the hazard of death before kidney failure adjusted for age, sex, primary kidney disease, hypertension, diabetes, cardiovascular disease, dyslipidemia, body mass index, albuminuria, CKD stage, medication adherence, use of renin-angiotensin system inhibitors, urate lowering therapy, spironolactone and anti-platelet agents, all at baseline (Model 1 in Table 3); Panel (B) Estimated effect of baseline uric acid value on the hazard of death before kidney failure, adjusted for the same factors as in (A) (Model 3 in Table 3). The reference value of uric acid for the HR indicated in the y-axis was arbitrarily chosen at 5 mg/gl, which corresponds to the midpoint of the normal range of uric acid (uric acid in mg/dl to µmol/l: × 59.48). Results from cause-specific Cox models using a spline function for uric acid. CKD-REIN cohort, France, 2013–2018 (n = 2344, including 218 death before kidney failure).
the cohort, patients with a cUA of precisely 3 mg/dl, had an 80% increased mortality at that time compared to patients with a cUA of precisely 5 mg/dl at the same time (HR 1.80, 95%CI 1.47, 2.19, Model 1 in Table 3). Mortality was twice higher for patients with a cUA of precisely 11 mg/dl compared to patients with a cUA of precisely 5 mg/dl at the same time (HR 2.30, 95%CI 1.84, 2.87, Fig. 3A, Model 1 in Table 3). The U-shaped association was similar in men and women (p-value for interaction of 0.68), but the confidence intervals were large in women (Figure S8) because of much less deaths before kidney failure in them (80 vs. 184 in men). As for kidney failure, further adjustment for salt and protein intake at baseline weakly affected the association (Model 2 vs. Model 1 in Table 3) and the results were very similar after adjusting or not for eating habits in this subsample (Table S2). Finally, the association was also much weaker with baseline UA than with cUA (Model 3 in Table 3, Fig. 3B).

Discussion

Using longitudinal data analysis, our results highlight the strong nonlinear association between longitudinal UA and both kidney failure and mortality in CKD patients. After adjustment for major risk factors for CKD progression, the hazard of kidney failure increased with increasing cUA, with a plateau between 6 and 10 mg/dl. By contrast, mortality before kidney failure had a U-shape relationship with cUA, with a minimum for cUA of 6 mg/dl, and a mortality twice higher for cUA of 3 or 11 mg/dl, compared to 5 mg/dl. The association with UA at inclusion was much weaker for both kidney failure and mortality.

Previous experimental studies showed that hyperuricemia may cause and accelerate CKD<sup>25,33</sup>, by mitochondrial dysfunction<sup>34</sup>, activation of the renin-angiotensin-aldosterone system<sup>35</sup>, induction of afferent arteriolar sclerosis<sup>36,37</sup>, pro-inflammation, or urate crystals deposition in the tubules<sup>38,39</sup>. Several epidemiologic studies also found that hyperuricemia was associated with CKD progression<sup>6,7,9-11,14,40,41</sup>, but most of them used baseline UA only, and only two examined the shape of the relationship between baseline UA and kidney failure and death<sup>10,14</sup>. As us, they found no or weak association between baseline UA and the hazard of kidney failure after adjustment for baseline eGFR. These results for baseline UA contrasted with our results for cUA suggesting a strong increased hazard of kidney failure for cUA above 11 mg/dl. The contrast of results between baseline and current UA value may explain why others studies did not find any association between UA and kidney failure<sup>2-5</sup>. The stronger association with cUA than with baseline UA may also suggest a stronger short term association than long term association, but this should be further explored using specific lags in the statistical analysis. Furthermore, the moderately stronger increase in the hazard of kidney failure in women than in men, associated with any increase of cUA till 7 mg/dl is consistent with previous studies which found a significant association between UA and CKD progression in women only<sup>41</sup>.

The U-shape of association between UA and death has already been found in a population of dialysis patients<sup>42</sup> or in a Korean population without CKD at baseline<sup>43</sup>, but this was investigated using baseline UA and not cUA. Similarly, another study found a significant increased all-cause mortality in non-diabetic patients with severe CKD and UA below 5 mg/dl<sup>44</sup>. Two other studies rather found a J-shape association in non-dialysis CKD patients, with an increased mortality for any baseline UA above 9 mg/dl<sup>10</sup> or above 11 mg/dl<sup>14</sup>, after adjustment for baseline eGFR. A potential explanation of our results suggesting an increased mortality at low cUA values, is that UA is involved in reducing oxidative stress and that a moderate increase in UA is needed to counteract oxidative damage, particularly in the context of arteriosclerosis. Indeed, in patients with atherosclerotic risks, hypouricemia has been associated with a higher hazard of all-cause and cardiovascular mortality<sup>45</sup>. Similarly, dialysis patients with hypouricemia are also at greater risk of all-cause and cardiovascular death than patients with normal uric acid levels<sup>46,47</sup>. In contrast, in non-diabetic CKD patients, Lee et al. did not find any association between hypouricemia and cardiovascular mortality, using an unique measure of UA at baseline<sup>44</sup>. It may therefore be of interest to replicate our analyzes in CKD patients by focusing on cardiovascular mortality. Another explanation for an increased mortality associated with low UA could be that patients with low UA had a poor nutritional status, as in hemodialysis patients<sup>48</sup> or in elderly patients<sup>49</sup>. However, further adjustment for salt and protein intake did not change the magnitude of association between cUA and mortality.

Our study has several strengths. First, the CKD-REIN cohort is a large multicentric prospective cohort, based on a nationally representative sample of nephrology clinics which is likely to enable adequate statistical power and generalizability of our findings to all French patients with CKD under nephrology care<sup>15</sup>. Second, thanks to frequent measures of UA and the use of advanced statistical methods, we were able to account for changes in UA over time and to estimate the nonlinear association between the reached level of UA and the hazard of death and kidney failure, which was much stronger than with baseline UA only.

However, our study has also limitations. The major limitation is that it is based on observational data, so we cannot exclude residual confounding, despite the use of DAGs to identify appropriate sets of adjusting factors. In particular, we adjusted for diet (salt and protein intake in particular) in a sensitivity analysis only because it was available for only half of patients. The estimated effect of cUA was weaker after adjustment, but this was more likely due to a selection bias in the subsample of half patients, than to a strong confounding bias. Indeed, adjusting or not for salt and protein intake in the subsample produced very similar results. However, it would be important to further investigate the role of diet in general in the association. It would also be of interest to investigate if the association is similar in patients with and without mutations or common variants in *UMOD*. Another limitation of our study is the use of a two-stage statistical approach instead of a joint analysis of UA trajectory and hazards of kidney failure and death<sup>21,50</sup>. Such a joint analysis would have accounted for potential informative dropout from the study in the estimation of the current value of UA and for its uncertainty. However, at the time of our analysis, the R packages that allowed the estimation of joint models assumed linear effects of biomarkers and thus did not allow nonlinear effects to be investigated. Because of these limitations and the complex and bidirectional relationship between changes in UA and kidney function during the course of CKD, we acknowledge that our study does not fully elucidate the causal role of UA in CKD progression. This remains
all the more unclear that a recent meta-analysis of published placebo-controlled clinical trials⁵¹, including the two most recent⁵²,⁵³, found no evidence of benefits of ULT on the risk of kidney failure. These findings, as well as our results confirming the increased mortality for low UA⁴⁴,⁴⁶,⁴⁷,⁵⁴, question the use of ULT to slow progression of CKD in patients with UA below 10 mg/dl.

To conclude, our study investigating longitudinal UA rather than baseline UA only, indicates a strong non-linear monotonic association between the reached level of UA and the hazard of kidney failure, and confirms the U-shape relationship with mortality. Although the use of ULT to slow progression of CKD is not recommended, we believe that a UA above 10 mg/dl may be considered as a strong risk marker for kidney failure and death, and thus should encourage nephrologists to be stricter in controlling cardiovascular and nephroprotective factors.

## Data availability

The datasets generated and/or analysed during the current study are not publicly available due to privacy/ethical restrictions but are available from the study principal investigator (Benedicte Stengel: benedicte.stengel@inserm.fr) on reasonable request.

Received: 1 February 2022; Accepted: 3 March 2023

Published online: 09 March 2023

## References

1. Modification of Diet in Renal Disease Study Group, Hunsicker LG, Adler S, Caggiula A, England BK, Greene T, et al. Predictors of the progression of renal disease in the Modification of Diet in Renal Disease Study. *Kidney Int.* **51**(6), 1908–1919 (1997).

2. Sturm G, Kollerits B, Neyer U, Ritz E, Kronenberg F, MMKD Study Group. Uric acid as a risk factor for progression of non-diabetic chronic kidney disease? The mild to Moderate Kidney Disease (MMKD) Study. *Exp. Gerontol.* **43**(4), 347–352 (2008).

3. Madero, M. et al. Uric acid and long-term outcomes in CKD. *Am. J. Kidney Dis.* **53**(5), 796–803 (2009).

4. Liu, W. C. et al. Association of hyperuricemia with renal outcomes, cardiovascular disease, and mortality. *Clin. J. Am. Soc. Nephrol.* **7**(4), 541–548 (2012).

5. Nacak, H. et al. Uric acid is not associated with decline in renal function or time to renal replacement therapy initiation in a referred cohort of patients with Stage III, IV and V chronic kidney disease. *Nephrol. Dial. Transpl.* **30**(12), 2039–2045 (2015).

6. Chang, W. X. et al. Predictors and the subsequent risk of end-stage renal disease—usefulness of 30% decline in estimated GFR over 2 years. *PLoS ONE* **10**(7), e0132927 (2015).

7. Uchida, S. et al. Targeting uric acid and the inhibition of progression to end-stage renal disease—a propensity score analysis. *PLoS ONE* **10**(12), e0145506 (2015).

8. Altemtam, N., Russell, J. & El Nahas, M. A study of the natural history of diabetic kidney disease (DKD). *Nephrol Dial Transpl.* **27**(5), 1847–1854 (2012).

9. Nacak, H., van Diepen, M., de Goeij, M. C., Rotmans, J. I. & Dekker, F. W. Uric acid: association with rate of renal function decline and time until start of dialysis in incident pre-dialysis patients. *BMC Nephrol.* **15**(1), 91 (2014).

10. Srivastava, A., Kaze, A. D., McMullan, C. J., Isakova, T. & Waikar, S. S. Uric acid and the risks of kidney failure and death in individuals with CKD. *Am. J. Kidney Dis.* **71**(3), 362–370 (2018).

11. Shi, Y. et al. Clinical outcome of hyperuricemia in IgA nephropathy: A Retrospective Cohort Study and Randomized Controlled Trial. *Kidney Blood Press Res.* **35**(3), 153–160 (2012).

12. Zoccali, C. & Mallamaci, F. Uric acid in chronic kidney disease: The quest for causality continues. *Nephrol Dial Transpl.* **33**(2), 193–195 (2018).

13. Bonino, B., Leoncini, G., Russo, E., Pontremoli, R. & Viazzi, F. Uric acid in CKD: Has the jury come to the verdict?. *J. Nephrol.* **33**(4), 715–724 (2020).

14. Tsai, C. W. et al. Uric acid predicts adverse outcomes in chronic kidney disease: a novel insight from trajectory analyses. *Nephrol. Dial Transpl.* **33**(2), 231–241 (2018).

15. Stengel, B. et al. The French chronic kidney disease-renal epidemiology and information network (CKD-REIN) cohort study. *Nephrol Dial Transpl.* **29**(8), 1500–1507 (2014).

16. KDIGO. KDIGO. Clinical practice guideline for the evaluation and management of chronic kidney disease. *Kidney Int. Suppl.* **2013**(3), 1–150 (2012).

17. Ikizler, T. A. et al. KDOQI clinical practice guideline for nutrition in CKD: 2020 update. *Am. J. Kidney Dis.* **76**(3), S1–S107 (2020).

18. Girerd, X., et al. Evaluation de l'observance par l'interrogatoire au cours du suivi des hypertendus dans des consultations spécialisées. *Arch. Mal. Coeur. Vaiss.* **94**(8), 839–842 (2001).

19. Chew, L. D. et al. Validation of screening questions for limited health literacy in a large VA outpatient population. *J. Gen. Intern. Med.* **23**(5), 561–566 (2008).

20. Dafni, U. & Tsiatis, A. Evaluating surrogate markers of clinical outcome when measured with error. *Biometrics* **54**(4), 1445–1462 (1998).

21. Sweeting, M. J. & Thompson, S. G. Joint modelling of longitudinal and time-to-event data with application to predicting abdominal aortic aneurysm growth and rupture. *Biom J* **53**(5), 750–763 (2011).

22. Ye, W., Lin, X. & Taylor, J. M. G. Semiparametric modeling of longitudinal measurements and time-to-event data—a two-stage regression calibration approach. *Biometrics* **64**(4), 1238–1246 (2008).

23. Harrell, F. E. General aspects of fitting regression models. In *Regression Modeling Strategies: With Applications to Linear Models, Logistic and Ordinal Regression, and Survival Analysis* (eds Harrell, J. & Frank, E.) 13–44 (Springer, New York, 2015).

24. Suttorp, M. M., Siegerink, B., Jager, K. J., Zoccali, C. & Dekker, F. W. Graphical presentation of confounding in directed acyclic graphs. *Nephrol. Dial. Transpl.* **30**(9), 1418–1423 (2015).

25. Sanchez-Lozada, L. G. et al. Uric acid and hypertension: An update with recommendations. *Am. J. Hypertens.* **33**(7), 583–594 (2020).

26. Hediger, M. A., Johnson, R. J., Miyazaki, H. & Endou, H. Molecular physiology of urate transport. *Physiology* **20**(2), 125–133 (2005).

27. Kalantar-Zadeh, K. & Fouque, D. Nutritional management of chronic kidney disease. *N. Engl. J. Med.* **377**(18), 1765–1776 (2017).

28. Ko, G. J., Rhee, C. M., Kalantar-Zadeh, K. & Joshi, S. The effects of high-protein diets on kidney health and longevity. *J. Am. Soc. Nephrol.* **31**(8), 1667–1679 (2020).

29. Paul, R. & Zaihra, T. Interval estimation of risk difference for data sampled from clusters. *Stat. Med.* **27**(21), 4207–4220 (2008).

30. Eilers, P. H. C. & Marx, B. D. Flexible smoothing with B-splines and penalties. *Stat. Sci.* **11**, 89–121 (1996).

31. Leffondre, K., Jager, K. J., Boucquemont, J., Stel, V. S. & Heinze, G. Representation of exposures in regression analysis and interpretation of regression coefficients: Basic concepts and pitfalls. *Nephrol. Dial. Transpl.* **29**(10), 1806–1814 (2014).
32. R Core Team (2019). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available from: https://www.R-project.org/

33. Johnson, R. J., Kivlighn, S. D., Kim, Y. G., Suga, S. & Fogo, A. B. Reappraisal of the pathogenesis and consequences of hyperuricemia in hypertension, cardiovascular disease, and renal disease. *Am. J. Kidney Dis.* **33**(2), 225–234 (1999).

34. Sánchez-Lozada, L. G. et al. Uric acid-induced endothelial dysfunction is associated with mitochondrial alterations and decreased intracellular ATP concentrations. *Nephron. Exp. Nephrol.* **121**(3–4), e71–78 (2012).

35. Mazzali, M. et al. Elevated uric acid increases blood pressure in the rat by a novel crystal-independent mechanism. *Hypertension* **38**(5), 1101–1106 (2001).

36. Mazzali, M. et al. Hyperuricemia induces a primary renal arteriolopathy in rats by a blood pressure-independent mechanism. *Am. J. Physiol. Renal. Physiol.* **282**(6), F991–997 (2002).

37. Sánchez-Lozada, L. G. et al. Mild hyperuricemia induces vasoconstriction and maintains glomerular hypertension in normal and remnant kidney rats. *Kidney Int.* **67**(1), 237–247 (2005).

38. Bjornstad, P. et al. Role of bicarbonate supplementation on urine uric acid crystals and diabetic tubulopathy in adults with type 1 diabetes. *Diabetes Obes. Metab.* **20**(7), 1776–1780 (2018).

39. Bjornstad, P. et al. Hyperfiltration and uricosuria in adolescents with type 1 diabetes. *Pediatr. Nephrol. Berl. Ger.* **31**(5), 787–793 (2016).

40. Iseki, K. et al. Significance of hyperuricemia as a risk factor for developing ESRD in a screened cohort. *Am. J. Kidney Dis.* **44**(4), 642–650 (2004).

41. Tsai, C. W., Lin, S. Y., Kuo, C. C. & Huang, C. C. Serum uric acid and progression of kidney disease: A longitudinal analysis and mini-review. *PLoS ONE* **12**(1), e0170393 (2017).

42. Zawada, A. M. et al. Serum uric acid and mortality risk among hemodialysis patients. *Kidney Int. Rep.* **5**(8), 1196–1206 (2020).

43. Cho, S. K., Chang, Y., Kim, I. & Ryu, S. U-shaped association between serum uric acid level and risk of mortality. *Arthritis Rheumatol.* **70**(7), 1122–1132 (2018).

44. Lee, C. L. & Tsai, S. F. Association between mortality and serum uric acid levels in non-diabetes-related chronic kidney disease: An analysis of the National Health and Nutrition Examination Survey, USA, 1999–2010. *Sci. Rep.* **10**(1), 17585 (2020).

45. Cang, Y. et al. Serum uric acid revealed a U-shaped relationship with all-cause mortality and cardiovascular mortality in high atherosclerosis risk patients: The ASSURE study. *Front. Cardiovasc. Med.* **8**, 641513 (2021).

46. Li, Q., Wu, C., Kuang, W., Zhan, X. & Zhou, J. Correlation analysis of low-level serum uric acid and cardiovascular events in patients on peritoneal dialysis. *Int. Urol. Nephrol.* **53**(11), 2399–2408 (2021).

47. Li, M. et al. Low serum uric acid levels increase the risk of all-cause death and cardiovascular death in hemodialysis patients. *Ren. Fail.* **42**(1), 315–322 (2020).

48. Beberashvili, I. et al. Longitudinal study of serum uric acid, nutritional status, and mortality in maintenance hemodialysis patients. *Clin. J. Am. Soc. Nephrol.* **11**(6), 1015–1023 (2016).

49. Tseng, W. C. et al. Taiwan Geriatric Kidney Disease (TGKD) Research Group. U-shaped association between serum uric acid levels with cardiovascular and all-cause mortality in the elderly: The role of malnourishment. *J. Am. Heart Assoc.* **7**(4), e007523 (2018).

50. Rizopoulos, D. *Joint Models for Longitudinal and Time-to-Event Data: With Applications in R* 1st edn. (CRC Press, New-York, 2012).

51. Chen, Q. et al. Effect of urate-lowering therapy on cardiovascular and kidney outcomes: A systematic review and meta-analysis. *Clin. J. Am. Soc. Nephrol.* **15**(11), 1576–1586 (2020).

52. Badve, S. V. et al. Effects of allopurinol on the progression of chronic kidney disease. *N. Engl. J. Med.* **382**(26), 2504–2513 (2020).

53. Doria, A. et al. Serum urate lowering with allopurinol and kidney function in type 1 diabetes. *N. Engl. J. Med.* **382**(26), 2493–2503 (2020).

54. Perez-Gomez, M. V. et al. Potential dangers of serum urate-lowering therapy. *Am. J. Med.* **132**(4), 457–467 (2019).

## Acknowledgements

We acknowledge the CKD-REIN study coordination staff for their efforts in setting up the CKD-REIN cohort, including Julie de Chefdebien, Reine Ketchemin and all the clinical research associates. We also thank all the patients and CKD-REIN clinical sites and investigators, by region: Alsace: Prs T. Hannedouche et B. Moulin (CHU, Strasbourg) and Dr A. Klein (CH Colmar); Aquitaine: Pr C. Combe (CHU, Bordeaux), Dr J.P. Bourdenx (Clinique St Augustin, Bordeaux), Dr A. Keller, Dr C. Delclaux (CH, Libourne), Dr B. Vendrely (Clinique St Martin, Pessac), Dr B. Deroure (Clinique Delay, Bayonne) and Dr A. Lacraz (CH, Bayonne); Basse Normandie: Dr T. Lobbedez (CHU, Caen) and Dr I. Landru (CH, Lisieux); Ile de France: Pr Z. Massy (CHU, Boulogne-Bilancourt), Pr P. Lang (CHU, Créteil), Dr X. Belenfant (CH, Montreuil), Pr E. Thervet (CHU, Paris), Dr P. Urena (Clinique du Landy, St Ouen) and Dr M. Delahousse (Hôpital Foch, Suresnes); Languedoc-Roussillon: Dr C. Vela (CH, Perpignan); Limousin: Dr M. Essig (CHU, Limoges); Lorraine: Dr H. Sekhri, Dr M.Smati (CH, Epinal); Dr M. Jamali and Dr B. Hacq (Clinique Louis Pasteur, Essey-les-Nancy), Dr V. Panescu, Dr M. Bellou (Polyclinique de Gentilly, Nancy) and Pr Luc Frimat (CHU, Vandœuvre-les-Nancy); Midi-Pyrénées: Pr L. Rostaing and Pr D. Chauveau (CHU, Toulouse); Nord-Pas-de-Calais: Prs C. Noël et F. Glowacki (CHU, Lille), Dr L. Vrigneaud (CH, Valenciennes), Dr R. Azar (CH, Dunkerque) and Dr M. Hoffmann (Hôpital privé La Louvière, Lille); Pays-de-la Loire: Pr M. Hourmant (CHU, Nantes), Dr A. Testa (Centre de Dialyse, Rezé) and Dr D. Besnier (CH, St Nazaire); Picardie: Pr G. Choukroun (CHU, Amiens) and Dr G. Lambrey (CH, Beauvais); Provence-Alpes-Côte d'Azur: Pr S. Burtey (CHU, Marseille), Dr G. Lebrun (CH, Aix-en-Provence) and Dr E. Magnant (Polyclinique du Parc Rambot, Aix-en-Provence); Rhône-Alpes: Pr M. Laville, Pr D. Fouque (CHU, Lyon-Sud) et L. Juillard (CHU Edouard Herriot, Lyon), Dr C. Chazot (Centre de rein artificiel Tassin Charcot, Ste Foy-les-Lyon), Pr P. Zaoui (CHU, Grenoble) and Dr F. Kuentz (Centre de Santé Rénale, Grenoble).

## Author contributions

M.P-R., J.H. and K.L. drafted the article. C.C., D.F., L.F., C.J., M.L., Z.A.M., C.L., C.A., R.P.-C., S.L. and B.S. and all the members of the CKD-REIN study group contributed to the study design and/or the conception of questionnaires and/or the biosample collection design. M.P-R., J.H. and K.L. performed the statistical analyses. All authors approved the final version of the manuscript.

## Competing interests

CKD-REIN is supported by a public-private partnership with funding from seven pharmaceutical companies as listed above. The partners did not participate in this work nor in the writing of the manuscript. CKD-REIN
is funded by the 'Agence Nationale de la Recherche' through the 2010 'Cohortes-Investissements d'Avenir' programme and by the 2010 national 'Programme Hospitalier de Recherche Clinique'. CKD-REIN is also supported through a public-private partnership with Amgen, Baxter, Fresenius Medical Care, GlaxoSmithKline (GSK), Merck Sharp & Dohme-Chibret (MSD France) since 2012, Lilly France since 2013, Otsuka Pharmaceutical since 2015 and Sanofi-Genzyme from 2012 to 2015. Inserm Transfert set up and has managed this partnership since 2011. The remaining authors do not have any conflict of interest to declare.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-023-30902-7.

**Correspondence** and requests for materials should be addressed to M.P.-R.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023

**CKD-REIN study group**

Natalia Alencar de Pinho⁹, Carole Ayav¹⁴, Christian Combe⁴,⁵, Denis Fouque⁶, Luc Frimat⁷,⁸, Yves-Edouard Herpe<fnref n="20" />, Christian Jacquelinet⁹,¹⁰, Maurice Laville¹¹,¹², Ziad A. Massy⁹,¹³, Christophe Pascal²¹, Roberto Pecoits-Filho¹⁵,¹⁶, Joost Schanstra²²,²³, Bénédicte Stengel⁹, Céline Lange⁹, Oriane Lambert⁹, Sophie Liabeuf⁹,¹⁷,¹⁸, Marie Metzger⁹ & Elodie Speyer⁹

<footnote n="20">²⁰Biobanque de Picardie, Département de la Recherche, CHU Amiens-Picardie, 80000 Amiens, France. ²¹Institut de formation et de recherche sur les organisations sanitaires et sociales (IFROSS), Université Jean Moulin Lyon 3, 69007 Lyon, France. ²²Institut National de la Santé et de la Recherche Médicale (INSERM), U1297, Institute of Cardiovascular and Metabolic Disease, Toulouse, France. ²³Université Toulouse III Paul-Sabatier, Toulouse, France.</footnote>