Benefits of Incremental Hemodialysis Seen in a Historical Cohort Study.

PURPOSE: Previous research on incremental hemodialysis transition has mainly focused on one or two benefits or prognoses. We aimed to conduct a comprehensive analysis by investigating whether incremental hemodialysis was simultaneously associated with adequate dialysis therapy, stable complication indicators, long-lasting arteriovenous vascular access, and long-lasting preservation of residual kidney function (RKF) without increasing mortality or hospitalization. PATIENTS AND METHODS: Incident hemodialysis patients from Huashan Hospital in Shanghai, China, over the period of 2012 to 2019, were enrolled and followed every three months until death or the time of censoring. Changes in complication indicators from baseline to all post-baseline visits were analyzed by mixed-effects models. The outcomes of RKF loss, arteriovenous vascular access complications, and the composite of all-cause mortality and cardiovascular events were compared between incremental and conventional hemodialysis by Cox proportional hazards model.
RESULTS: Of the 113 patients enrolled in the study, 45 underwent incremental and 68 conventional hemodialysis. There were no significant differences in the changes from baseline to post-baseline visits in complication indicators between the two groups. Incremental hemodialysis reduced the risks of RKF loss (HR, 0.33; 95% CI, 0.14-0.82), de novo arteriovenous access complication (HR, 0.26; 95% CI, 0.08-0.82), and recurrent arteriovenous access complications under the Andersen-Gill (AG) model (HR, 0.27; 95% CI, 0.10-0.74) and the Prentice, Williams and Peterson Total Time (PWP-TT) model (HR, 0.31; 95% CI, 0.12-0.80). There were no significant differences in all-cause hospitalization or the composite outcome between groups.
CONCLUSION: Incremental hemodialysis is an effective dialysis transition strategy that preserves RKF and arteriovenous access without affecting dialysis adequacy, patient stability, hospitalization risk and mortality risk. Randomized controlled trials are warranted.

Introduction

Thrice-weekly hemodialysis has been considered a standard renal replacement therapy (RRT) since the 1970s.1,2 Recently, an increasing number of patients have progressed to end-stage renal disease (ESRD) and initiated hemodialysis therapy with residual kidney function (RKF).3,4 Thus, controversy has arisen concerning the abrupt transition from non-dialysis-dependent chronic kidney disease to a thrice-weekly fixed hemodialysis schedule.2,4–8 Incremental hemodialysis, a strategy in which the hemodialysis dose is gradually intensified according to the degree of RKF among incident hemodialysis patients, was proposed.2,5,8 In most cases, incremental hemodialysis refers to incident hemodialysis patients being prescribed hemodialysis from once weekly to twice weekly and then thrice weekly according to RKF.5,8

It has been reported that incremental hemodialysis was associated with preservation of RKF,6,9,10 lower serum β2-microglobulin,11,12 and lower erythropoiesis-stimulating agent consumption and erythropoietin resistance index.10,11 Recently, the clinical outcomes of incremental hemodialysis were studied, and mortality was the main focus.6,10,13 To date, a comprehensive analysis of whether incremental hemodialysis could simultaneously maintain sufficient adequacy of dialysis and the stability of complication indicators while delaying the loss of residual kidney function without increasing mortality or cardiovascular disease has not been reported. Moreover, protection of arteriovenous vascular access has been proposed but not examined in a cohort of incremental hemodialysis patients.14,15

In China, a total of 26% of ESRD patients have been dialyzing twice weekly, compared with fewer than 5% in other Dialysis Outcomes and Practice Pattern Study (DOPPS) regions.16 Thus, evidence on the safety and benefits of incremental hemodialysis is even more urgent in China than in any other country. Our center has implemented incremental hemodialysis for nearly 10 years with careful and regular monitoring of the symptoms and signs of patients. Therefore, we conducted a historical cohort study to comprehensively investigate the association of incremental hemodialysis with RKF loss, complications of arteriovenous access, all-cause hospitalization and the composite outcome of all-cause mortality and cardiovascular events among incident hemodialysis patients. We also explored the changes in complication-related clinical and laboratory values over the first year among patients on incremental hemodialysis and then compared all the values between the incremental and the conventional hemodialysis.

Patients and Methods

Study Population

We performed a historical cohort study using the data of incident hemodialysis patients aged over 18 years in the Hemodialysis Centre, Huashan Hospital, Shanghai, China, between 2012 and 2019. Of the 137 incident hemodialysis patients who were treated with in-center hemodialysis, 9 were excluded for being treated temporarily, and 15 were excluded for having missing data on urine output or having urine output less than 200 mL/day at baseline (). Our final study population consisted of 113 patients. Patients who initiated once- or twice-weekly hemodialysis for at least six months and who subsequently transitioned to a thrice-weekly regimen as RKF declined were categorized as the incremental group. Those who underwent thrice-weekly hemodialysis or transitioned from infrequent hemodialysis to thrice-weekly hemodialysis within six months were categorized as the conventional group. Six months was decided all from clinical practice for it covered the transition from non-dialysis CKD to dialysis and the evaluation for the hemodialysis adequacy and all the complications. It is an enough period to evaluate whether hemodialysis regimens fitted patients. Those who had undergone once-weekly or twice-weekly hemodialysis for palliative purposes or during recovery from acute kidney injury were not included in our cohort. Patients were followed up from the first day of hemodialysis to the end of the second year or the day of censoring for death, transferring to another center, transplantation, or the end of the study (January 17, 2020).

Demographic, Clinical, and Laboratory Values

Information for age, sex, body mass index (BMI), comorbid conditions, and the primary causes of kidney disease at baseline was obtained from the electronic database of the Hemodialysis Centre. The baseline estimated glomerular filtration rate (eGFR) before hemodialysis initiation was calculated using the Modification of Diet in Renal Disease (MDRD) equation, and the urine output was measured at the same time. All patients visited the doctor with laboratory tests for the regular assessment of complications in January, April, July, and October every year. The results of the first regular visit, which was within three months after hemodialysis initiation, were considered as the baseline.

Blood samples were drawn using uniform techniques in our center at the second hemodialysis session or the first hemodialysis session when patients receiving more-than once-weekly hemodialysis or once-weekly hemodialysis, respectively, and transported to the Department of Clinical Laboratory, typically within two hours. Laboratory values were measured by automated and standardized methods. Furthermore, patients were encouraged to collect 24-hour urine samples on the last day of the interdialytic interval before the hemodialysis session when laboratory tests were to be performed at the upcoming regular visit, from which residual kidney urea clearance (KRU) was calculated. In addition, the values of ultrafiltration volume, pre-dialysis weight, systolic blood pressure, diastolic blood pressure, and dialysis hours for each patient during each regular-visit month were averaged for the following analysis of this study. The single-pool Kt/V delivered by hemodialysis was calculated by the second-generation Daugirdas’ equation.17 The urea reduction ratio (URR) was the fractional reduction of urea during hemodialysis.

Clinical Outcomes

Outcomes of interest were as follows: 1) RKF loss; 2) time to first complication and time to recurrent complications of arteriovenous vascular access; 3) all-cause hospitalization; and 4) the composite outcome of mortality and cardiovascular events. Cardiovascular events were defined as attacks due to myocardial infarction, congestive heart failure, hospitalization for myocardial ischemia, arrhythmia, or stroke. RKF loss was defined as a urine output less than 200 mL per 24 hours. Types of arteriovenous access included arteriovenous fistula and graft. Vascular access complications referred to failure, stenosis, or thrombosis of the arteriovenous vascular access that necessitated surgical interventions, including percutaneous transluminal angioplasty. All medical events that happened in either our hospital or other hospitals were documented in our center’s database in routine clinical practice.

Statistical Analysis

Continuous variables were expressed as mean ± standard deviation (SD) or median with interquartile range (IQR), and categorical variables were expressed as percentages or ratios. The differences between groups were examined using the t-test, Wilcoxon’s rank sum test, the χ2 test or Fisher’s exact test, as appropriate. The trends of laboratory values across the first four regular visits were examined by Cuzick’s test. The incidence rate was expressed as the number of cases per 100 person-years. The decline rate of urine output was defined as the slope between the urine output at baseline and the last follow-up urine volume over time.

Changes from baseline to all post-baseline visits in blood pressure, log-transformed NT-proBNP, hemoglobin, potassium, carbon dioxide combining power (CO2CP), calcium, phosphate and log-transformed intact parathyroid hormone (iPTH) were analyzed by a maximum-likelihood, mixed-effects repeated-measures models (MMRM). The model included the above parameter as a response variables and treatment, visit, and the treatment × visit interaction, age, sex, BMI, diabetes, baseline urine output and eGFR as covariates.

The risk (hazard ratio, HR) of RKF loss, the risk of the composite outcome of all-cause mortality and cardiovascular events, and the time to the first complication of arteriovenous access were compared between the incremental and the conventional hemodialysis by the Kaplan–Meier curve, the Log rank test, and the Cox proportional hazards model. Arteriovenous access complications could occur more than once for an individual patient. The time to recurrent complications of arteriovenous access was examined by the Andersen–Gill model (AG model) and the Prentice, Williams and Peterson Total Time model (PWP-TT model).18 The incidence rate differences in all-cause hospitalization or arteriovenous access complications were examined by the negative binomial model. To adjust the cofounders between groups, two levels of adjustment were performed in all models: 1) Model 1, which included age, sex, BMI, diabetes, urine output, eGFR and systolic blood pressure at baseline and 2) Model 2, which included covariates in Model 1 plus log-transformed NT-proBNP, albumin and C-reactive protein at baseline.

All statistical analyses were performed with Stata, version 16.0 (StataCorp LLC). A P-value less than 0.05 was considered statistically significant.

Results

Baseline Characteristics

The cohort included 45 patients undergoing incremental hemodialysis and 68 undergoing conventional hemodialysis. Baseline characteristics between the two groups are shown in Table 1. Compared with the conventional group, patients in the incremental group were younger and had a lower prevalence of diabetes. The baseline ultrafiltration volume, systolic blood pressure and NT-proBNP were lower, while serum albumin and urine output at baseline were higher in the incremental group.

Table 1

Baseline Characteristics of Incident Hemodialysis Patients by Hemodialysis Regimen

Characteristics	Incremental Group (N = 45)	Conventional Group (N = 68)	P
Age (years)	56.3 ± 14.3	61.9 ± 13.8	0.04
Male, n (%)	22 (48.9)	41 (60.3)	0.23
Body mass index (kg·m−2)	22.1 ± 2.8	22.5 ± 3.9	0.64
Primary kidney disease			0.01
Diabetes, n (%)	4 (8.9)	25 (36.8)
Hypertension, n (%)	5 (11.1)	7 (10.3)
Glomerulonephritis, n (%)	25 (55.6)	21 (30.9)
Polycystic kidney disease, n (%)	6 (13.3)	6 (8.8)
Others/Unknown, n (%)	5 (11.1)	9 (13.2)
Comorbidity
Diabetes mellitus, n (%)	8 (17.8)	26 (38.2)	0.02
Cardiovascular disease, n (%)	6 (13.3)	15 (22.1)	0.24
Dialysis-related index
Urine output (mL/24h)	1566.2 ± 533.9	1129.4 ± 521.9	<0.001
eGFR (mL/min/1.73m2)	6.71 (5.11, 8.97)	6.39 (4.72, 9.45)	0.90
Ultrafiltration (mL)	1268.3 ± 730.3	1975.6 ± 700.5	<0.001
Predialysis weight (kg)	60.7 ± 10.3	61.6 ± 13.6	0.71
SBP (mmHg)	131.3 ± 14.9	140.8 ± 16.0	0.002
DBP (mmHg)	76.4 ± 10.6	78.1 ± 10.6	0.42
Dialytic treatment time (hours/session)	3.98 ± 0.32	3.94 ± 0.12	0.31
spKt/V	1.35 ± 0.26	1.27 ± 0.33	0.13
URR	0.69 ± 0.07	0.65 ± 0.10	0.03
Laboratory measures
Albumin (g/dL)	4.02 ± 3.37	3.72 ± 4.35	<0.001
Hemoglobin (g/dL)	9.69 ± 1.11	9.54 ± 1.33	0.55
Calcium (mg/dL)	8.77 ± 0.91	8.79 ± 1.02	0.90
Phosphate (mg/dL)	4.85 ± 1.09	5.29 ± 1.68	0.13
iPTH (pg/mL)	204 (107, 336)	220 (119, 336)	0.85
CO2CP (mmol/L)	21.6 ± 4.2	21.5 ± 3.9	0.85
Potassium (mg/dL)	4.81 ± 0.75	4.67 ± 0.74	0.34
NT-proBNP (pg/mL)	1161 (733, 4021)	2867 (1124, 6196)	<0.001
Creatinine (mg/dL)	8.46 ± 2.12	8.20 ± 2.44	0.55
SUN (mg/dL)	77.2 ± 18.1	71.0 ± 18.8	0.09
C-reactive protein			0.20
<10 mg/dL (%)	43 (95.6)	59 (86.8)
≥10 mg/dL (%)	2 (4.4)	9 (13.2)

Notes: Continuous data presented as mean ± standard deviation or median (interquartile range); categorical data presented as number (percent).

Abbreviations: eGFR, estimated glomerular filtration rate; iPTH, intact parathyroid hormone; CO2CP, carbon dioxide combining power; NT-proBNP, N-terminal pro–B-type natriuretic peptide; SUN, serum urea nitrogen; URR, urea reduction ratio.

Comparisons of Clinical and Laboratory Values Between Groups

In the incremental group, the percentages of those who transitioned to thrice-weekly hemodialysis were 17.1% (95% CI, 8.5%–32.6%) and 31.5% (95% CI, 19.1%–49.1%) at the end of the first and second years, respectively, estimated by Kaplan–Meier curve (). In the first year of incremental hemodialysis, most values including blood pressure, serum potassium, iPTH, CO2CP, and NT-proBNP remained stable across the four regular visits (P for trend > 0.05; Table 2). The exceptions were that ultrafiltration volume and creatinine had upward trends (P for trend < 0.05; Table 2). Moreover, there were no significant differences in the overall changes from baseline to all post-baseline visits in NT-proBNP, blood pressure, calcium, phosphate, iPTH, hemoglobin, potassium and CO2CP (P > 0.05; Table 3). In the subgroup analysis, no significant difference was found in stdKt/V over the first four visits (Pfor group = 0.14; Pfor group×time interaction = 0.59; ).

Table 2

The Changes in Clinical and Laboratory Values Over the First Four Regularly Visits After Initiating Hemodialysis Therapy Among Patients with Incremental Hemodialysis

Characteristics	Q1 (n = 45)	Q2 (n = 45)	Q3 (n = 43)	Q4 (n = 41)	P for trend
Ultrafiltration (mL)	1278 (670, 1686)	1339 (752, 2314)	1771 (748, 2814)	1721 (826, 2375)	0.02
Predialysis weight (kg)	60.7 ± 10.3	61.5 ± 10.5	62.7 ± 10.4	62.3 ± 10.3	0.42
SBP (mmHg)	131.3 ± 14.9	131.5 ± 15.2	130.9 ± 14.6	132.9 ± 14.4	0.78
DBP (mmHg)	76.4 ± 10.6	76.3 ± 12.0	78.0 ± 12.1	77.2 ± 9.7	0.58
Dialytic treatment time (hours)	3.98 ± 0.32	4.04 ± 0.32	4.03 ± 0.37	4.00 ± 0.33	0.78
spKt/V	1.35 ± 0.26	1.35 ± 0.25	1.35 ± 0.25	1.38 ± 0.27	0.76
URR	0.69 ± 0.07	0.68 ± 0.07	0.68 ± 0.06	0.69 ± 0.07	0.97
Laboratory variables
Albumin (g/dL)	4.02 ± 0.34	4.04 ± 0.34	4.03 ± 0.34	4.06 ± 0.28	0.71
Hemoglobin (g/dL)	9.69 ± 1.11	10.64 ± 1.26	10.58 ± 1.39	10.32 ± 1.16	0.05
Calcium (mg/dL)	8.77 ± 0.91	9.05 ± 1.05	8.86 ± 1.01	9.06 ± 1.29	0.42
Phosphate (mg/dL)	4.85 ± 1.09	5.11 ± 1.45	5.15 ± 1.44	4.67 ± 1.53	0.63
iPTH (pg/mL)	204 (107, 336)	174 (64, 391)	181 (72, 302)	154 (93, 222)	0.07
CO2CP (mmol/L)	21.6 ± 4.2	20.8 ± 4.1	20.6 ± 3.2	21.7 ± 3.3	0.77
Potassium (mg/dL)	4.81 ± 0.75	4.69 ± 0.60	4.89 ± 0.72	4.74 ± 0.75	0.80
NT-proBNP (pg/mL)	1161 (733, 2595)	1100 (607, 1931)	1252 (738, 2183)	1454 (918, 3378)	0.17
Creatinine (mg/dL)	8.46 ± 2.12	9.29 ± 2.51	10.25 ± 3.08	10.15 ± 3.00	0.001
SUN (mg/dL)	77.2 ± 18.1	79.1 ± 16.5	84.9 ± 20.2	84.9 ± 20.2	0.36
C-reactive protein					0.24
<10 mg/dL (%)	43 (95.6)	44 (97.8)	41 (95.3)	37 (90.2)
≥10 mg/dL (%)	2 (4.4)	1 (2.2)	2 (4.7)	4 (9.8)

Notes: Continuous data presented as mean ± standard deviation or median (interquartile range); categorical data presented as number (percent). Q1–4, the first to fourth quarterly visits after hemodialysis initiation.

Abbreviations: eGFR, estimated glomerular filtration rate; iPTH, intact parathyroid hormone; CO2CP, carbon dioxide combining power; NT-proBNP, N-terminal pro–B-type natriuretic peptide; SUN, serum urea nitrogen; URR, urea reduction ratio.

Table 3

Repeated-Measures Analysis of Changes in Laboratory Values and Blood Pressure Between Groups from Baseline to Post-Baseline Visits

	Second Visit		Third Visit		Fourth Visit
	Diff (95% CI)a	Pb	Diff (95% CI)a	Pb	Diff (95% CI)a	Pb	Pc
Log-transformed NT-proBNP
Conventional group (n = 68)	0.10 (−0.09, 0.30)	0.44	−0.01 (−0.21, 0.19)	0.75	0.10 (−0.11, 0.30)	0.46	0.50
Incremental group (n = 45)	−0.02 (−0.26, 0.22)		0.04 (0.21, 0.29)		0.22 (−0.03, 0.47)
Systolic blood pressure
Conventional group (n = 68)	−2.82 (−5.84, 0.20)	0.22	−2.09 (−5.12, 0.95)	0.59	−3.09 (−6.23, 0.05)	0.07	0.30
Incremental group (n = 45)	0.20 (−3.51, 3.91)		−0.77 (−4.53, 2.99)		1.43 (−2.39, 5.25)
Diastolic blood pressure
Conventional group (n = 68)	−1.74 (−3.64, 0.16)	0.28	−1.31 (−3.22, 0.60)	0.08	−1.69 (−3.66, 0.28)	0.09	0.25
Incremental group (n = 45)	−0.08 (−2.41, 2.25)		1.45 (−0.92, 3.82)		1.03 (−1.37, 3.43)
Potassium
Conventional group (n = 68)	0.14 (−0.02, 0.31)	0.05	0.24 (0.07, 0.42)	0.23	0.18 (0.01, 0.35)	0.04	0.15
Incremental group (n = 45)	−0.12 (−0.32, 0.09)		0.08 (−0.13, 0.29)		−0.11 (−0.32, 0.10)
CO2CP
Conventional group (n = 68)	−0.49 (−0.52, 1.49)	0.10	−0.03 (−1.04, 0.98)	0.23	0.29 (−0.75, 1.33)	0.76	0.33
Incremental group (n = 45)	−0.85 (−2.08, 0.38)		−1.01 (−2.26, 0.25)		0.04 (−1.23, 1.30)
Hemoglobin
Conventional group (n = 68)	0.92 (0.51, 1.33)	0.92	1.10 (0.70, 1.51)	0.56	1.03 (0.61, 1.45)	0.23	0.55
Incremental group (n = 45)	0.95 (0.45, 1.45)		0.91 (0.40, 1.41)		0.63 (0.12, 1.14)
Calcium
Conventional group (n = 68)	0.33 (0.04, 0.62)	0.84	0.41 (0.12, 0.70)	0.20	0.26 (−0.04, 0.56)	0.89	0.48
Incremental group (n = 45)	0.28 (−0.07, 0.64)		0.11 (−0.25, 0.47)		0.29 (−0.07, 0.65)
Phosphate
Conventional group (n = 68)	−0.08 (−0.48, 0.33)	0.31	−0.08 (−0.49, 0.33)	0.31	0.10 (−0.33, 0.52)	0.32	0.15
Incremental group (n = 45)	0.26 (−0.24, 0.76)		0.26 (−0.25, 0.77)		−0.24 (−0.76, 0.28)
Log-transformed iPTH
Conventional group (n = 68)	−0.48 (−0.81, −0.15)	0.44	−0.46 (−0.79, −0.13)	0.32	−0.22 (−0.57, 0.12)	0.57	0.39
Incremental group (n = 45)	−0.27 (−0.68, 0.14)		−0.19 (−0.60, 0.22)		−0.38 (−0.80, 0.03)

Notes: aValues are adjusted least-squares means and 95% CI estimated from the mixed-effects model. bTest of significance of treatment group differences by visit (changes from baseline visit to the second visit, third visit or fourth visit between group) from the mixed-effects model. cTest of significance between treatment groups for the overall effect (group×visit) from the mixed-effects model.

Abbreviations: iPTH, intact parathyroid hormone; CO2CP, carbon dioxide combining power; NT-proBNP, N-terminal pro–B-type natriuretic peptide.

Association of Incremental Hemodialysis with Residual Kidney Function

The decline rate of urine output in the incremental group was significantly lower than that in the conventional group [−33.4 mL/month (IQR, −76.7 to −19.4 mL/month) vs −70.4 mL/month (IQR, −155.5 to −33.0 mL/month); P < 0.001]. The mean follow-up duration was 13.4 ± 7.9 months; during this period, 46 patients occurred outcomes of RKF loss, 37 occurred in the conventional group. A decrease in unadjusted risk for RKF loss was observed for incremental hemodialysis (P log-rank < 0.001, Figure 1A; HR, 0.23, 95% CI, 0.11–0.47, P < 0.001, Table 4). After controlling for covariates, this association persisted in Model 1 (HR, 0.37; 95% CI, 0.16–0.82; P = 0.01) and Model 2 (HR, 0.33; 95% CI, 0.14–0.82; P = 0.02) (Table 4). The association between hemodialysis regimens and RKF was diminished after ultrafiltration volume at baseline was added to the adjustment model.

Figure 1

Kaplan–Meier survival curve for (A) RKF loss, (B) the first arteriovenous access complication, (C) death; (D) the composite outcome of mortality and cardiovascular events. Numbers without parentheses in the risk table show the absolute numbers of patients at risk over time and numbers within parentheses show the absolute numbers of patients with events during a period.

Table 4

Unadjusted and Adjusted Hazard Ratio of RKF Loss by Hemodialysis Regimen

Model	Hazard Ratio	95% Confidential Interval	P
Univariate	0.23	0.11, 0.47	<0.001
Model 1	0.37	0.16, 0.82	0.01
Model 2	0.33	0.14, 0.82	0.02
Model 1 + UFV	0.45	0.20, 1.04	0.06
Model 2 + UFV	0.40	0.16, 1.05	0.06

Notes: Model 1, adjusted for age, sex, BMI, diabetes, urine output, eGFR, and systolic blood pressure at baseline; Model 2, adjusted for Model 1 plus NT-proBNP, albumin, and C-reactive protein at baseline.

Abbreviation: UFV, ultrafiltration volume at baseline.

Association of Incremental Hemodialysis with Complications of Arteriovenous Access

Three patients in the incremental group and six patients in the conventional group were excluded from the main cohort for failure of arteriovenous access maturation or censoring before arteriovenous access maturation. The mean follow-up duration was 17.0 ± 7.9 months; during this period, 36 arteriovenous outcome complications occurred, 30 occurred in the conventional group. The incidence rates of arteriovenous access complications in the two groups were 8.5 and 29.8 per 100 person-years, respectively. Incremental hemodialysis was associated with a lower incidence rate of arteriovenous access complications under the unadjusted and adjusted negative binomial models (). A lower risk of having an arteriovenous access complication was also observed for incremental hemodialysis by the Log rank test and the Cox model, with an adjusted HR of 0.26 (95% CI, 0.08–0.82) in Model 2 (Figure 1B; Table 5). Moreover, the adjusted risk of recurrent complications of arteriovenous access was lower for incremental hemodialysis under the AG model (HR, 0.27; 95% CI, 0.10–0.74; P = 0.01) and the PWP-TT model (HR, 0.31; 95% CI, 0.12–0.80; P = 0.02) (Table 5).

Table 5

Unadjusted and Adjusted Hazard Ratio of Arteriovenous Access Complications by Hemodialysis Regimen

Model	Hazard Ratio	95% Confidential Interval	P
Cox Model
Univariate	0.36	0.13, 0.98	0.04
Model 1	0.27	0.09, 0.81	0.02
Model 2	0.26	0.08, 0.82	0.02
AG Model
Univariate	0.29	0.11, 0.79	0.02
Model 1	0.25	0.10, 0.63	0.003
Model 2	0.27	0.10, 0.74	0.01
PWP-TT Model
Univariate	0.37	0.16, 0.86	0.02
Model 1	0.30	0.14, 0.66	0.003
Model 2	0.31	0.12, 0.80	0.02

Notes: Model 1, adjusted for age, sex, BMI, diabetes, urine output, eGFR, and systolic blood pressure at baseline; Model 2, adjusted for Model 1 plus NT-proBNP, albumin, and C-reactive protein at baseline.

Abbreviations: AG Model, Andersen and Gill Cox model; PWP-TT Model, Prentice, Williams and Peterson-Total Time Model.

Association of Incremental Hemodialysis with Hospitalization and the Composite Outcome

The mean follow-up duration was 18.9 ± 6.9 months; during this period, 18 patients occurred outcomes of the composite outcome, 14 occurred in the conventional group. The incidence rate of all-cause hospitalization was lower in the incremental group than in the conventional group under the unadjusted negative binomial model (IRR, 0.57; 95% CI, 0.34–0.93; P = 0.03), but this association was not significant under the adjusted negative binomial model (). Moreover, the incremental group was not associated with risk of death or the composite outcome by the unadjusted or adjusted Cox model (Figure 1C and D; Table 6).

Table 6

Unadjusted and Adjusted Hazard Ratio of the Death and the Composite Outcome of All-Cause Mortality and Cardiovascular Events by Hemodialysis Regimen

Model	Hazard Ratio	95% Confidential Interval	P
Composite Outcome
Univariate	0.41	0.13, 1.24	0.11
Model 1	0.80	0.23, 2.80	0.72
Model 2	0.88	0.23, 3.43	0.85
Death
Univariate	0.30	0.04, 2.59	0.28
Model 1	0.96	0.08, 12.02	0.99
Model 2	1.83	0.06, 52.25	0.72

Notes: Model 1, adjusted for age, sex, BMI, diabetes, urine output, eGFR, and systolic blood pressure at baseline; Model 2, adjusted for Model 1 plus NT-proBNP, albumin, and C-reactive protein at baseline.

Discussion

In the present study, we provide evidence that patients receiving incremental hemodialysis can maintain most of their complication-related clinical and laboratory values stable in the first year. We also found that incremental hemodialysis was significantly associated with preservation of RKF and fewer complications of arteriovenous access. Furthermore, patients receiving incremental hemodialysis were found to have a similar risk of hospitalization, death and the composite outcome as patients receiving conventional hemodialysis.

Most clinical and laboratory values in the first year remained stable among patients with incremental hemodialysis. In a subgroup analysis, we found that the weekly clearance of urea (stdKt/V) between the two groups was comparable, which means RKF compensated for the reduction in urea clearance when dialysis frequency decreased (). RKF has been reported to have a strong ability to remove fluid and salt, excrete phosphorus, clear middle molecules, and produce endogenous vitamin D and erythropoietin.4,19–21 The increasing trends in ultrafiltration volume and creatinine might still suggest an underlying decline in RKF during the first year. Unlike a gradual or mild decline in RKF, a significant decline in RKF usually results in uncontrolled hypertension or hypervolemia, and frequent hyperpotassemia and acidosis, which are the indexes in our center to increase hemodialysis frequency among patients receiving incremental hemodialysis. Taken together, RKF plays a significant role in implementing incremental hemodialysis.

So far, there has been no evidence that incremental hemodialysis was associated with prolonged longevity of arteriovenous access. In the present study, we first report the association of incremental hemodialysis with a decreased risk of complications of arteriovenous access. This association may arise due to lower intradialytic ultrafiltration and a lower incidence of intradialytic hypotension (IDH). Patients with incremental hemodialysis had better preserved RKF, which may slow the accumulation of fluids and avoid larger-volume intradialytic ultrafiltration. Large-volume ultrafiltration has been directly associated with the frequency of episodes of IDH, which may cause extreme hemodynamic derangements and thus increased the risk of arteriovenous access complications.22–24 In addition, more frequent venipunctures in conventional hemodialysis result in repeated trauma at venipuncture sites and subsequent turbulent blood flow. Both of these consequences may promote profibrotic cytokine production, local inflammation, and neointimal proliferation, which together increase arteriovenous access complications.25,26

Several studies have compared mortality between the incremental and the conventional hemodialysis. Some studies reported that the incremental hemodialysis did not increase mortality.6,13 In the present study, the outcome of all-cause hospitalization and the composite outcome of all-cause mortality and cardiovascular events were investigated. We found that the composite outcome, death and all-cause hospitalization were not significantly different between groups. This may be related to a stable control of complications, adequate hemodialysis and better RKF preservation as demonstrated in our study. Taken together, the above findings show that patients with RKF may safely be prescribed incremental hemodialysis in their transition to renal replacement therapy.

It is plausible that the observed safety of transition was moderated by longer preservation of RKF.27–29 In a previous study, we showed that thrice-weekly hemodialysis in the first year of hemodialysis was associated with a 7.2-fold increased likelihood of RKF loss compared to twice-weekly hemodialysis.9 Here, we demonstrated that the decline rate of RKF was slower in the incremental hemodialysis than in the conventional hemodialysis. Notably, this benefit was not traded off against a poorer prognosis as showed in this study. One of the mechanisms whereby RKF among patients with incremental hemodialysis declines more slowly may be that renal perfusion acutely decreases during hemodialysis therapy, resulting in renal ischemia, and that an incremental hemodialysis schedule decreases the frequency of ischemic insults and alleviates renal tissue damage at the initial stage.30,31 Membrane bio-incompatibility is one of the causes of inflammation in hemodialysis patients and is an independent risk factor for the progression of kidney function.32–34 Incremental hemodialysis may reduce exposure to those factors, which might partially account for RKF preservation.

Overall, it was demonstrated in our study that incremental hemodialysis reduces complications of arteriovenous access and slows the decline of RKF. As we known, hemodialysis patients have referred to their vascular access as their lifeline, which means the vascular access complications may have an adverse impact on hemodialysis experience and health-related quality of life (HRQOL), which may be one of the reasons why HRQOL impairs among patients on hemodialysis. In addition, high RKF was significantly associated with fewer symptom burden.35 Overall, incremental hemodialysis may improve the health-related quality of life and symptom burden among hemodialysis patients, whereby the protection of arteriovenous access and residual kidney function. Subsequent studies may need more consideration in HRQOL and symptom burden.

There were some limitations to our study. This was a historical cohort study, in which two hemodialysis regimens were not assigned randomly; thus, the results may suffer from some confounders such as age, comorbidities and so on. However, we adjusted potential confounders and tried to investigate the independent associations of incremental hemodialysis with clinical outcomes. We also obtained consistent results using the inverse probability weighting method, which is another method to adjust for the baseline imbalance between groups ( and ). However, residual confounding factors, such as medications, may still exist. Small sample study was another limitation and low statistical power of small sample size may result in poor ability to detect clinical outcomes of significance. Multicenter study with larger sample size may be warranted in the near future. Moreover, not all patients had KRU records at every regular visit, but estimated urine output was monitored. The standard weekly clearance of urea (stdKt/V) thus cannot be calculated for all patients. In addition, trajectories of renal function before initiation of long-term hemodialysis have not been considered, which might have an association with the clinical outcome after initiating hemodialysis. Nevertheless, we are the first to report a protective effect of incremental hemodialysis on arteriovenous access. Our comparisons of clinical outcomes between groups are more comprehensive than other studies. We also have summarized a set of clinical criteria in our center to implement incremental hemodialysis, which is based on the experience of 10-year clinical practices and need further verification from cohort studies ().

Conclusion

Incremental hemodialysis is a strategy to keeping hemodialysis adequate and complications stable and to protect RKF and arteriovenous access without increasing all-cause hospitalization and the composite outcome of mortality and cardiovascular events. This study provides better evidence supporting the implementation of incremental hemodialysis through comprehensive analyses. Our findings lead to an urgent call for randomized controlled trials to support incremental hemodialysis schedule.