Hemodialysis and its effects on glymphatic system function in end-stage kidney disease patients
Introduction
Chronic kidney disease (CKD) is a growing global health concern, affecting millions of individuals worldwide. With an incidence rate of approximately 24.2 per 100,000 people and an annual growth rate of around 8%, CKD represents a significant public health challenge (1). End-stage kidney disease (ESKD), the final stage of CKD, is typically diagnosed when the glomerular filtration rate (GFR) falls below 15 mL/min/1.73 m2. One of the major complications of ESKD is cognitive impairment, with prevalence rates ranging from approximately 16% to 38%, nearly three times higher than those observed in the age-matched general population (2). Although the mechanisms underlying cognitive impairment in ESKD are not fully understood, growing evidence suggests that disruptions in brain waste clearance pathways, such as the glymphatic system, may contribute to this phenomenon (3).
The glymphatic system, a type of cerebrospinal fluid (CSF) transport system, is employed to remove waste from the brain (4,5). The glymphatic system facilitates the movement of CSF from the arterial perivascular spaces into the brain parenchyma through aquaporin-4 (AQP4) channels expressed in astrocytic endfeet. From there, interstitial fluid and waste products are transported along venous perivascular routes and cleared through meningeal and cervical lymphatic vessels (6). Currently, methods for evaluating the glymphatic system are relatively limited. Magnetic resonance imaging (MRI) with intrathecal gadolinium-based contrast agents is the most representative technique for assessing glymphatic function (7). However, its clinical application is limited by its invasiveness and procedural complexity, particularly due to the contraindication of contrast agents in patients with renal impairment. Recently, a diffusion tensor image analysis along the perivascular space (DTI-ALPS) index has been developed, which utilizes non-invasive diffusion imaging to assess glymphatic function (8). The DTI-ALPS method measures water diffusion along perivascular spaces, reflecting the efficiency of CSF and interstitial fluid exchange, which is a key component of glymphatic clearance. A higher ALPS index may reflect more efficient glymphatic system function, whereas a lower index may suggest reduced clearance efficiency. Multiple studies have demonstrated that a reduced ALPS index is negatively correlated with diseases such as Alzheimer’s disease, Parkinson’s disease, epilepsy, stroke, and obstructive sleep apnea–hypopnea syndrome, where glymphatic dysfunction contributes to waste accumulation and neurodegeneration (9-13).
Studies have demonstrated that both early and end stages of CKD are linked with abnormalities in the glymphatic system (14,15). This glymphatic dysfunction is associated with the cognitive impairment observed in CKD patients (16). However, there is limited research on whether short-term improvements in glymphatic function occur following hemodialysis in ESKD patients.
We hypothesized that acute hemodynamic shifts and rapid clearance of small water-soluble uremic toxins during hemodialysis might transiently modulate glymphatic function. Therefore, the aim of this study was to assess the glymphatic system in pre- and post-hemodialysis ESKD patients as well as compare them to healthy controls (HCs) using the DTI-ALPS method. This non-invasive and sensitive technique allows for a detailed assessment of glymphatic function and provides insight into the potential short-term effects of hemodialysis on this system. The findings of this study may provide valuable insights into the mechanisms underlying cognitive impairment in ESKD patients. We present this article in accordance with the STROBE reporting checklist (available at https://qims.amegroups.com/article/view/10.21037/qims-2025-74/rc).
Methods
Participants
This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendment and was approved by the Ethics Committee of Tongren People’s Hospital (No. 202304). Informed consent was provided by all individual participants. We enrolled 25 patients with ESKD from Tongren People’s Hospital and recruited 27 age- and gender-matched HCs from November 2023 to June 2024. The study flow chart is illustrated in Figure 1.
The inclusion criteria for this study were as follows: adult ESKD patients aged 18–75 years, undergoing regular hemodialysis with a stable clinical condition, no acute illness or infection within the past month, and the ability to undergo MRI scans both before and after hemodialysis. Age- and sex-matched HCs with normal renal function, no history of kidney disease or chronic conditions affecting brain health, and no recent infections or acute illnesses within the past month were also recruited. Exclusion criteria for all participants included neurological or psychiatric disorders (e.g., stroke, dementia, major depressive disorder) that could impact brain function, a history of traumatic brain injury, MRI contraindications (e.g., metal implants, claustrophobia, pacemakers), pregnancy or breastfeeding, and current alcohol or substance abuse. Additional exclusion criteria for ESKD patients included unstable medical conditions, such as poorly controlled hypertension or diabetes, and recent changes in hemodialysis regimen within the past month.
Laboratory data collection and analysis
Blood samples were collected from ESKD patients at two time points to evaluate a panel of laboratory variables: (I) immediately before hemodialysis (pre-hemodialysis samples) and (II) immediately after completion of the hemodialysis session (post-hemodialysis samples). Pre-hemodialysis samples were obtained from the dialysis tubing immediately following the initiation of hemodialysis sessions to ensure consistency in sampling conditions.
All laboratory analyses were performed using standardized, automated analytical platforms following manufacturer protocols. Biochemical assays, including tests for renal function [blood urea nitrogen (BUN), creatinine], liver function [aspartate aminotransferase (AST), alanine aminotransferase (ALT)], electrolytes (potassium, sodium, chloride, calcium), iron metabolism (iron, ferritin), lipid profile (total cholesterol), and nutritional markers (albumin, total protein, prealbumin), were conducted using the Beckman Coulter AU5800 automated biochemistry analyzer (Beckman Coulter, Brea, CA, USA). Hematological parameters (hemoglobin, hematocrit) were measured using the Mindray BC-7500[NR]CRP automated hematology analyzer (Mindray, Shenzhen, China). Parathyroid hormone (PTH) levels were measured using an electrochemiluminescence immunoassay (ECLIA) on a Roche Cobas e601 analyzer (Roche Diagnostics, Basel, Switzerland). The estimated glomerular filtration rate (eGFR) was calculated based on serum creatinine levels using the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation (17).
To ensure analytical accuracy, all blood samples were processed within 2 hours of collection. Samples were transported under controlled conditions to the laboratory, where they were centrifuged, and plasma or serum was separated as required for specific assays. Quality control measures were implemented following the manufacturers’ guidelines, and all instruments were calibrated prior to analysis.
MRI acquisition
MRI data were acquired using a Philips Ingenia CX 3.0T scanner (Philips, Amsterdam, Netherlands). During MRI scanning, participants were instructed to remain still in a supine position with their eyes closed, maintain steady breathing, stay awake, and refrain from engaging in deliberate thought. To reduce noise, earplugs or headphones were provided, and foam cushions were placed on either side of the head to minimize head movement.
HCs underwent a single MRI scan, whereas ESKD patients underwent two scans: one before and one after hemodialysis. The pre-hemodialysis MRI scan was performed 2–24 hours prior to hemodialysis, and the post-hemodialysis MRI scan occurred 24–48 hours afterward. DTI data were acquired using a single-shot echo planar imaging (EPI) sequence with the following parameters: repetition time/echo time (TR/TE) =9,700/74.2 ms, b-values =0 and 1,000 s/mm2, 32 diffusion directions, field of view (FOV) =256 mm × 256 mm, matrix size =128×128, 48 slices, 3 mm slice thickness without an intersection gap, and in-plane resolution 2×2 mm2.
Participants with excessive motion or artifacts were excluded. In addition, conventional MRI scans, including axial and sagittal T2, as well as axial T2 fluid-attenuated inversion recovery (FLAIR), were performed to rule out organic brain abnormalities.
DTI-ALPS processing
DTI-ALPS data processing was conducted following the methodology outlined in our previous study (9). Digital Imaging and Communications in Medicine (DICOM) images were initially converted to Neuroimaging Informatics Technology Initiative (NIfTI) format using the MRIcron tool (https://www.nitrc.org/projects/mricron). During this conversion, the corresponding b-vector and b-value files, essential for DTI-ALPS processing, were generated. Preprocessing of the DTI images was conducted using the FSL pipeline in version 6.0 of the FMRIB Software Library (FSL; https://fsl.fmrib.ox.ac.uk/fsl) (18,19).
The preprocessing workflow included several key steps. Artifact correction was performed using the Marchenko-Pastur Principal Component (MP-PCA) denoising algorithm and Gibbs-unringing, implemented with MRtrix3 commands “dwidenoise” and “mrdegibbs”. Eddy current and motion corrections were then applied using the “eddy” command in FSL. Head motion was quantified using the mean absolute intervolume displacement (MAID) (20,21), calculated as the average Euclidean distance of each volume’s translational displacement relative to the first volume, derived from the transformation matrices in eddy correction outputs. Skull stripping was conducted to remove non-brain tissue from the images, followed by the generation of fractional anisotropy (FA) and diffusivity maps along x-, y-, and z-axes using the FSL function “dtifit”.
Next, each participant’s FA map was co-registered to the Johns Hopkins University-International Consortium for Brain Mapping (JHU-ICBM-FA) template using FSL’s “flirt” command, and the transformation matrix was applied to the corresponding diffusivity maps. The projection and association fibers at the level of the lateral ventricle body, specifically the superior corona radiata (SCR) and the superior longitudinal fasciculus (SLF), were identified based on the JHU-ICBM-DTI-81-white-matter labeled atlas. Regions of interest (ROIs) were automatically defined as 5 mm-diameter spheres within the bilateral SLF and SCR on all participants’ diffusivity maps, as illustrated in Figure 2. The central coordinates of these ROIs, referenced to the JHU-ICBM-FA template, were as follows: left SCR (116, 110, 99), left SLF (128, 110, 99), right SCR (64, 110, 99), and right SLF (51, 110, 99) (19).
Diffusivity values (Dxx, Dyy, and Dzz) for the bilateral SLF and SCR were automatically extracted for the ALPS index calculation. Automatic extraction was achieved by registering each participant’s FA map to the JHU-ICBM-FA template and programmatically placing spherical ROIs at predefined template-based coordinates for the SLF and SCR, ensuring consistent and reproducible ROI positioning across participants. Manual verification was performed to ensure the accuracy of the registration and ROI placement for each participant.
The ALPS index was calculated as the ratio of the mean x-axis diffusivity in the projection fibers regions (Dxxproj) and the association fibers regions (Dxxassoc) to the mean of the y-axis diffusivity in the projection fibers regions (Dyyproj) and the z-axis diffusivity in the association fibers regions (Dzzassoc), as shown in Eq. [1] (8).
The final ALPS index was obtained by averaging the bilateral ALPS index values.
Sensitivity analysis
To address potential confounding effects of comorbidities on the ALPS index, we conducted a sensitivity analysis excluding ESKD patients with documented cardiovascular or systemic comorbidities. This subgroup analysis aimed to isolate the impact of kidney dysfunction itself on the ALPS index. Of the 25 ESKD patients, 13 had ≥1 comorbidity and were excluded, leaving 12 “comorbidity-free” ESKD patients for comparison with the original HC group (HCs, n=27). Statistical methods remained consistent with the primary analysis.
Statistical analysis
Statistical comparisons for group differences were conducted by applying the chi-square test for categorical variables. For continuous variables, the Shapiro-Wilk test was first used to assess the normality of distributions in each group. When both groups met the assumption of normality, Levene’s test was then performed to evaluate the homogeneity of variances. Based on the results, an appropriate statistical test was selected. Specifically, if both normality and homogeneity of variances were confirmed, an independent two-sample t-test was used. If normality was met, but variances were unequal, Welch’s t-test was applied to account for the variance difference. If either group did not satisfy the normality assumption, the Mann-Whitney U test was employed as a non-parametric alternative. In the analysis of ALPS index differences between groups, analysis of variance (ANOVA) was used with sex and age included as covariates to adjust for their potential effects.
Paired comparisons for measurements before and after hemodialysis were analyzed by first assessing the normality of paired differences with the Shapiro-Wilk test. If the paired differences followed a normal distribution, a paired t-test was used to compare means between samples; if normality was not met, the Wilcoxon signed-rank test was employed as a non-parametric alternative, allowing for deviations from normality.
The correlation between the ALPS index and laboratory metrics and clinical data was evaluated by using either Pearson or Spearman correlation coefficients, depending on distribution characteristics. The normality of both the ALPS index and laboratory metrics was first examined with the Shapiro-Wilk test. For pairs of variables with normal distributions, the Pearson correlation coefficient was used to measure the strength and direction of their linear relationship. For cases where one or both variables did not meet the normality assumption, the Spearman rank correlation coefficient was employed to assess the strength and direction of their monotonic relationship.
All statistical analyses were conducted using R software (version 4.3.2, R Foundation for Statistical Computing, Vienna, Austria) and GraphPad Prism (version 10.1.2, GraphPad Software, San Diego, CA, USA). A two-tailed P value of less than 0.05 was considered statistically significant. The false discovery rate (FDR) correction method was applied for multiple comparisons.
Results
Demographics and clinical characteristics
A total of 52 participants were included in this study, comprising 25 ESKD patients (mean age 49.8±13.01 years; 14 males) and 27 age- and sex-matched HCs (mean age 48.67±10.16 years; 12 males) (Table 1). Shapiro-Wilk tests were conducted to evaluate the normality of the relevant continuous variables, with the results presented in Table S1. There were no significant differences between ESKD patients and HCs in terms of age (U=335.5, P=0.978) and sex (χ2=0.308, P=0.579). In the laboratory data, calcium and eGFR levels increased in ESKD patients after hemodialysis compared to pre-hemodialysis levels, whereas BUN, creatinine, and potassium levels decreased following hemodialysis (Figure 3).
Table 1
| Clinical data | Pre-hemodialysis ESKD (n=25) | Post-hemodialysis ESKD (n=25) | HCs (n=27) | Test statistics | P value |
|---|---|---|---|---|---|
| Demographic data | |||||
| Age, years, mean (SD) | 49.80 (13.01) | 48.67 (10.16) | 335.5a | 0.978 | |
| Male, N (%) | 14 (56.0) | 12 (44.44) | 0.308b | 0.579 | |
| Hemodialysis duration, months, mean (SD) | 57.87 (40.38) | – | – | – | |
| Comorbidities | |||||
| Diabetes mellitus, N (%) | 4 (16.0) | 0 (0.0) | – | – | |
| Hypertension, N (%) | 11 (44.0) | 0 (0.0) | – | – | |
| Coronary artery disease, N (%) | 1 (4.0) | 0 (0.0) | – | – | |
| Atrial fibrillation, N (%) | 1 (4.0) | 0 (0.0) | – | – | |
| IgA nephropathy, N (%) | 1 (4.0) | 0 (0.0) | – | – | |
| Total number with ≥1 comorbidity, N (%) | 13 (52.0) | 0 (0.0) | – | – | |
| Head motion [MAID, mm, mean (SD)] | |||||
| Three-group comparison | 3.64 (1.07) | 3.60 (1.03) | 3.30 (0.83) | 1.07c | 0.351 |
| Pre- vs. post-hemodialysis | 3.64 (1.07) | 3.60 (1.03) | – | 0.31d | 0.761 |
| Laboratory data | |||||
| BUN, mmol/L, mean (SD) | 23.07 (9.06) | 7.33 (3.35) | – | 153e | <0.001* |
| Creatinine, μmol/L, mean (SD) | 1,036.65 (241.42) | 414.75 (109.83) | – | 16.75d | <0.001* |
| Hemoglobin, g/L, mean (SD) | 116.09 (22.41) | – | – | – | – |
| Hematocrit, %, mean (SD) | 35.15 (6.84) | – | – | – | – |
| Albumin, g/L, mean (SD) | 39.16 (4.03) | – | – | – | – |
| Total protein, g/L, mean (SD) | 65.18 (5.41) | – | – | – | – |
| Prealbumin, mg/L, mean (SD) | 323.32 (67.7) | – | – | – | – |
| Potassium, mmol/L, mean (SD) | 4.83 (0.89) | 3.63 (0.56) | – | 5.476d | 0.001* |
| Sodium, mmol/L, mean (SD) | 137.59 (3.54) | 136.43 (2.33) | – | 30e | 0.426 |
| Chloride, mmol/L, mean (SD) | 102.66 (4.62) | 100.91 (2.7) | – | 1.845d | 0.102 |
| Calcium, mmol/L, mean (SD) | 2.09 (0.28) | 2.20 (0.2) | – | –4.971d | 0.003* |
| Iron, μmol/L, mean (SD) | 6.87 (3.43) | – | – | – | – |
| Ferritin, ng/mL, mean (SD) | 234.01 (245.93) | – | – | – | – |
| eGFR, mL/min/1.73 m2, mean (SD) | 5.83 (2.2) | 12.93 (3.94) | – | 0e | <0.001* |
| AST, U/L, mean (SD) | 17.20 (5.79) | – | – | – | – |
| ALT, U/L, mean (SD) | 12.61 (6.87) | – | – | – | – |
| PTH, pg/mL, mean (SD) | 499.83 (355.99) | – | – | – | – |
| Total cholesterol, mmol/L, mean (SD) | 5.61 (4.57) | – | – | – | – |
a, Mann-Whitney U test; b, chi-squared test; c, analysis of variance; d, paired t-test; e, Wilcoxon signed-rank test. *, statistical significance (P<0.05). ALT, alanine aminotransferase; AST, aspartate aminotransferase; BUN, blood urea nitrogen; eGFR, estimated glomerular filtration rate; ESKD, end-stage kidney disease; HCs, healthy controls; MAID, mean absolute intervolume displacement; PTH, parathyroid hormone; SD, standard deviation.
Head motion analysis revealed no significant differences in MAID between pre-hemodialysis ESKD patients (MAID =3.64±1.07 mm), post-hemodialysis ESKD patients (MAID =3.60±1.03 mm), and HCs (MAID =3.30±0.83 mm) (ANOVA, P=0.351) (Table 1). Similarly, no significant difference was observed between pre- and post-hemodialysis scans in ESKD patients (paired t-test, P=0.761). Given the lack of significant differences, we did not include motion as a nuisance regressor in the statistical analysis.
Comparison of ALPS index
The DTI-ALPS values for the participant groups are summarized in Table 2 and Figure 4. The ALPS index was significantly lower in pre-hemodialysis ESKD patients compared to HCs (1.379 vs. 1.469, P=0.023). There was no significant difference in the ALPS index between pre- and post-hemodialysis in ESKD patients within a single session (1.379 vs. 1.387, P=0.493). Additionally, diffusivity along the x-axis and z-axis in the association fibers, as well as along the y-axis in the projection fibers, which are not considered to reflect glymphatic function based on the DTI-ALPS method (8), also showed significant differences between pre-hemodialysis ESKD patients and HCs (P<0.05).
Table 2
| Parameters | Pre-hemodialysis ESKD (n=25) | Post-hemodialysis ESKD (n=25) | HCs (n=27) | P value† | P value‡ | |||||
|---|---|---|---|---|---|---|---|---|---|---|
| Mean | SD | Mean | SD | Mean | SD | |||||
| Association fibers | ||||||||||
| Dxx | 0.00068 | 0.00008 | 0.00068 | 0.00009 | 0.00063 | 0.00005 | 0.553a | 0.043c | ||
| Dzz | 0.00043 | 0.00007 | 0.00044 | 0.00008 | 0.00038 | 0.00004 | 0.52a | 0.005c | ||
| Projection fibers | ||||||||||
| Dxx | 0.00061 | 0.00006 | 0.00061 | 0.00006 | 0.00060 | 0.00003 | 0.931b | 0.232c | ||
| Dyy | 0.00052 | 0.00009 | 0.00051 | 0.00008 | 0.00046 | 0.00004 | 0.543a | 0.011c | ||
| ALPS index | 1.37900 | 0.15370 | 1.38700 | 0.14330 | 1.46900 | 0.09020 | 0.493b | 0.023c | ||
a, Wilcoxon signed-rank test; b, Paired t-test; c, analysis of variance. †, comparisons of ALPS index between ESKD patients in the pre- and post-hemodialysis; ‡, comparisons of ALPS index between ESKD patients in the pre-hemodialysis and HCs, with sex and age included as covariates. ALPS index, diffusion tensor image analysis along the perivascular space index; Dxx, diffusivity along the x-axis; Dyy, diffusivity along the y-axis; Dzz, diffusivity along the z-axis; ESKD, end-stage kidney disease; HCs, healthy controls; SD, standard deviation.
After excluding 13 ESKD patients with comorbidities, the remaining subgroup (n=12) showed no significant difference in ALPS index compared to HCs (1.410 vs. 1.469, P=0.1297). Similarly, within this subgroup, no significant change was observed between pre- and post-hemodialysis measurements (1.410 vs. 1.417, P=0.5731) (Figure S1).
Correlation between ALPS index and clinical characteristics
In pre-hemodialysis ESKD patients, a significant positive correlation was observed between the ALPS index and prealbumin levels after FDR correction (r=0.663, uncorrected P=0.002, corrected P=0.036, FDR-corrected), as illustrated in Figure 5. No significant correlations were observed between the ALPS index and other laboratory parameters.
Exploratory analysis showed no significant correlations between pre-hemodialysis ALPS index and dialysis duration (months) (r=−0.180, P=0.390), age (r=−0.370, P=0.071), or cumulative number of hemodialysis sessions (r=−0.147, P=0.484). Similarly, the change in ALPS index (post- minus pre-hemodialysis) showed no significant associations with dialysis duration (months) (r=0.022, P=0.919), age (r=−0.185, P=0.387), or cumulative number of hemodialysis sessions (r=0.026, P=0.904).
Discussion
This study explored glymphatic system alterations in patients with ESKD and examined the potential impact of hemodialysis using the DTI-ALPS method. Our findings showed that ESKD patients exhibited lower ALPS indices compared to HCs, suggesting altered perivascular diffusivity. However, no significant change in the ALPS index was observed following a single hemodialysis session. Although these results may reflect impaired interstitial fluid dynamics in ESKD, the ALPS index remains an indirect imaging marker and does not directly quantify glymphatic clearance. Therefore, interpretations regarding glymphatic dysfunction should be made cautiously and considered within the broader context of cerebrovascular and systemic influences.
The post-hemodialysis changes observed in the laboratory results of this study, including increases in calcium and eGFR, as well as reductions in BUN, creatinine, and potassium, align with the expected biochemical responses to hemodialysis. These changes highlight the procedure’s effectiveness in waste clearance and electrolyte balance, key functions for maintaining systemic homeostasis in ESKD patients. However, although these biochemical markers improve with dialysis, they primarily reflect systemic rather than central nervous system (CNS) clearance. This distinction may explain the lack of improvement in ALPS index observed in our study. Hemodialysis effectively clears small, water-soluble molecules from the blood but is less efficient at removing larger, protein-bound uremic toxins that can cross the blood-brain barrier and accumulate in the brain (22). These retained toxins may continue to impair the glymphatic system by promoting neuroinflammation and oxidative stress, which can hinder interstitial fluid flow and reduce perivascular clearance (23). This distinction underscores the need for CNS-specific markers to evaluate dialysis efficacy fully and suggests that modified or adjunctive approaches could benefit CNS clearance in ESKD patients.
The reduced ALPS index in ESKD patients may indicate glymphatic system dysfunction, suggesting that CKD could profoundly impact brain health. In ESKD, the accumulation of uremic toxins and chronic inflammation may impair glymphatic clearance by crossing the blood-brain barrier, disrupting cerebral circulation, and limiting interstitial fluid movement through perivascular pathways (24). This dysfunction in brain waste clearance could lead to the buildup of neurotoxic proteins such as amyloid-beta, which are normally cleared by the glymphatic system and may partly explain the higher prevalence of cognitive decline in ESKD patients (25). This study confirmed that a single hemodialysis session does not enhance the ALPS index, potentially because hemodialysis cannot effectively remove certain toxins, such as middle molecules and protein-bound solutes, that contribute to chronic inflammation and oxidative stress, further damaging glymphatic function (26). Given that impaired glymphatic function has been associated with neurodegenerative conditions such as Alzheimer’s disease (27), these findings emphasize the potential role of glymphatic impairment as a neurological risk pathway in kidney disease. The lack of improvement in the ALPS index highlights the need for additional therapeutic strategies that target neuroinflammation and support glymphatic activity. Adjustments to hemodialysis protocols or complementary treatments may be critical for mitigating cognitive risks in this population. Further research is essential to uncover the mechanisms linking renal dysfunction to brain health and to develop integrated therapeutic strategies for ESKD patients.
The correlations between the ALPS index and clinical indicators enhance our understanding of the relationship between kidney function and glymphatic activity. Prealbumin, an important biomarker of nutritional status, is critical for waste clearance and neuroprotection against toxicity (28). In CKD and ESKD patients, prealbumin levels are commonly affected by malnutrition and chronic low-grade inflammation. The positive correlation between prealbumin and the ALPS index suggests that nutritional status may influence brain waste clearance mechanisms, indicating that better-nourished patients may exhibit more effective clearance capabilities. Low prealbumin levels often indicate malnutrition, which is prevalent in these populations and crucial for disease management and prognosis. Additionally, inflammation can suppress prealbumin production, meaning that low levels may reflect both nutritional deficiencies and an inflammatory state. This interplay between chronic inflammation and malnutrition, referred to as “cachexia syndrome”, can increase the risk of disease progression in CKD patients (29). Exploratory analyses also revealed no significant relationships between ALPS index or its change and dialysis duration, age, or cumulative number of hemodialysis sessions, suggesting that these factors may not be primary predictors of glymphatic function in this context; however, these null findings may reflect the small sample size and warrant further investigation with larger cohorts.
To further assess the influence of comorbidities, our sensitivity analysis excluding ESKD patients with comorbidities demonstrated that the difference in ALPS indices between ESKD patients and HCs was no longer significant. This suggests that cardiovascular and systemic comorbidities may partly contribute to impaired cerebral waste clearance in ESKD. Additionally, the reduced sample size after exclusion may have limited statistical power. Therefore, the ALPS index observed in ESKD patients should be interpreted as potentially reflecting the combined effects of kidney failure and coexisting comorbid conditions. Future studies with larger, comorbidity-stratified cohorts are warranted to better delineate these relationships.
This study has several limitations. First, the relatively small sample size and the use of a single time point for HCs may limit the generalizability of our findings and prevent the detection of inter-individual variability over time. Second, the timing of MRI acquisition relative to hemodialysis sessions (2–24 hours pre-hemodialysis and 24–48 hours post-hemodialysis) was constrained by clinical and logistical factors. This variability may have obscured transient dialysis-related changes in glymphatic activity. Third, we did not assess cognitive performance, which limits the interpretation of clinical relevance regarding neurological outcomes. Additionally, the ALPS index is an indirect measure of glymphatic activity based on perivascular water diffusivity and does not capture the full complexity of glymphatic clearance (30-32). Moreover, the use of spherical ROIs may have introduced partial volume effects, affecting measurement reliability. Ulloa et al. recommended using squared or cubic ROIs to achieve more consistent and reliable results (33). Future studies should adopt larger, stratified samples, tighter imaging schedules, and multimodal neuroimaging approaches (including direct functional assessments) to better elucidate the brain-kidney axis and its implications for cognitive health in ESKD.
Conclusions
This study found that patients with ESKD exhibited lower ALPS indices compared to HCs, suggesting potential alterations in perivascular diffusivity. However, this difference was no longer statistically significant after excluding patients with cardiovascular and systemic comorbidities, indicating that such conditions may also contribute to the observed changes. No significant change in the ALPS index was observed after a single hemodialysis session under the current scanning schedule; further studies with optimized scan timing are needed to clarify acute dialysis effects. Importantly, the ALPS index is an indirect imaging marker and does not directly reflect glymphatic clearance function. Therefore, our findings should be interpreted with caution. Future studies using larger, comorbidity-stratified cohorts and multimodal imaging approaches are needed to better understand the relationship between kidney dysfunction, systemic comorbidities, and brain fluid dynamics.
Acknowledgments
None.
Footnote
Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at https://qims.amegroups.com/article/view/10.21037/qims-2025-74/rc
Data Sharing Statement: Available at https://qims.amegroups.com/article/view/10.21037/qims-2025-74/dss
Funding: This study was supported by
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://qims.amegroups.com/article/view/10.21037/qims-2025-74/coif). X.Z. reports that this study was supported by the Guizhou Provincial High-level Innovative Talent Training Program (QianKeHe Platform Talent-GCC [2023] 083) and the National Natural Science Foundation of China (Grants Nos. 82460344 and 82060314). The other authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendment and was approved by the Ethics Committee of Tongren People’s Hospital (No. 202304). Informed consent was provided by all individual participants.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Kim GH. Pharmacologic Treatment of Chronic Hyperkalemia in Patients with Chronic Kidney Disease. Electrolyte Blood Press 2019;17:1-6. [Crossref] [PubMed]
- Kurella Tamura M, Yaffe K. Dementia and cognitive impairment in ESRD: diagnostic and therapeutic strategies. Kidney Int 2011;79:14-22. [Crossref] [PubMed]
- Xu S, Wang J, Sun K, Meng L, Qin C, Feng R, Tian Y, Zhai Y, Liang D, Zhang R, Tian H, Liu H, Chen Y, Fu Y, Chen P, Zhu Q, Teng J, Wang X. Cognitive Impairment in Chronic Kidney Disease Is Associated with Glymphatic System Dysfunction. Kidney Dis (Basel) 2023;9:384-97. [Crossref] [PubMed]
- Iliff JJ, Wang M, Liao Y, Plogg BA, Peng W, Gundersen GA, Benveniste H, Vates GE, Deane R, Goldman SA, Nagelhus EA, Nedergaard M. A paravascular pathway facilitates CSF flow through the brain parenchyma and the clearance of interstitial solutes, including amyloid β. Sci Transl Med 2012;4:147ra111. [Crossref] [PubMed]
- Jessen NA, Munk AS, Lundgaard I, Nedergaard M. The Glymphatic System: A Beginner's Guide. Neurochem Res 2015;40:2583-99. [Crossref] [PubMed]
- Rasmussen MK, Mestre H, Nedergaard M. The glymphatic pathway in neurological disorders. Lancet Neurol 2018;17:1016-24. [Crossref] [PubMed]
- Taoka T, Naganawa S. Glymphatic imaging using MRI. J Magn Reson Imaging 2020;51:11-24. [Crossref] [PubMed]
- Taoka T, Masutani Y, Kawai H, Nakane T, Matsuoka K, Yasuno F, Kishimoto T, Naganawa S. Evaluation of glymphatic system activity with the diffusion MR technique: diffusion tensor image analysis along the perivascular space (DTI-ALPS) in Alzheimer's disease cases. Jpn J Radiol 2017;35:172-8. [Crossref] [PubMed]
- Xiong Z, Bai M, Wang Z, Wang R, Tian C, Wang L, Nie L, Zeng X. Resting-state fMRI network efficiency as a mediator in the relationship between the glymphatic system and cognitive function in obstructive sleep apnea hypopnea syndrome: Insights from a DTI-ALPS investigation. Sleep Med 2024;119:250-7. [Crossref] [PubMed]
- Shen T, Yue Y, Ba F, He T, Tang X, Hu X, Pu J, Huang C, Lv W, Zhang B, Lai HY. Diffusion along perivascular spaces as marker for impairment of glymphatic system in Parkinson's disease. NPJ Parkinsons Dis 2022;8:174. [Crossref] [PubMed]
- Lee DA, Park BS, Ko J, Park SH, Lee YJ, Kim IH, Park JH, Park KM. Glymphatic system dysfunction in temporal lobe epilepsy patients with hippocampal sclerosis. Epilepsia Open 2022;7:306-14. [Crossref] [PubMed]
- Toh CH, Siow TY. Glymphatic Dysfunction in Patients With Ischemic Stroke. Front Aging Neurosci 2021;13:756249. [Crossref] [PubMed]
- Chang HI, Huang CW, Hsu SW, Huang SH, Lin KJ, Ho TY, Ma MC, Hsiao WC, Chang CC. Gray matter reserve determines glymphatic system function in young-onset Alzheimer's disease: Evidenced by DTI-ALPS and compared with age-matched controls. Psychiatry Clin Neurosci 2023;77:401-9. [Crossref] [PubMed]
- Heo CM, Lee DA, Park KM, Lee YJ, Park S, Kim YW, Ko J, Yoo BC, Park BS. Glymphatic system dysfunction in patients with early chronic kidney disease. Front Neurol 2022;13:976089. [Crossref] [PubMed]
- Heo CM, Lee WH, Park BS, Lee YJ, Park S, Kim YW, Lee DA, Yoo BC, Park KM. Glymphatic Dysfunction in Patients With End-Stage Renal Disease. Front Neurol 2021;12:809438. [Crossref] [PubMed]
- Viggiano D, Wagner CA, Martino G, Nedergaard M, Zoccali C, Unwin R, Capasso G. Mechanisms of cognitive dysfunction in CKD. Nat Rev Nephrol 2020;16:452-69. [Crossref] [PubMed]
- Inker LA, Eneanya ND, Coresh J, Tighiouart H, Wang D, Sang Y, et al. New Creatinine- and Cystatin C-Based Equations to Estimate GFR without Race. N Engl J Med 2021;385:1737-49. [Crossref] [PubMed]
- Jenkinson M, Beckmann CF, Behrens TE, Woolrich MW, Smith SM. FSL. Neuroimage 2012;62:782-90. [Crossref] [PubMed]
- Liu X, Barisano G, Shao X, Jann K, Ringman JM, Lu H, et al. Cross-Vendor Test-Retest Validation of Diffusion Tensor Image Analysis along the Perivascular Space (DTI-ALPS) for Evaluating Glymphatic System Function. Aging Dis 2024;15:1885-98. [Crossref] [PubMed]
- Power JD, Barnes KA, Snyder AZ, Schlaggar BL, Petersen SE. Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 2012;59:2142-54. [Crossref] [PubMed]
- Aoki Y, Yoncheva YN, Chen B, Nath T, Sharp D, Lazar M, Velasco P, Milham MP, Di Martino A. Association of White Matter Structure With Autism Spectrum Disorder and Attention-Deficit/Hyperactivity Disorder. JAMA Psychiatry 2017;74:1120-8. [Crossref] [PubMed]
- Cai H, Chen L. Advances in Research Related to Cognitive Impairment in Maintenance Hemodialysis Patients. Academic Journal of Medicine & Health Sciences 2023. doi:
10.25236/AJMHS.2023.040404 . - Cai Y, Zhang Y, Leng S, Ma Y, Jiang Q, Wen Q, Ju S, Hu J. The relationship between inflammation, impaired glymphatic system, and neurodegenerative disorders: A vicious cycle. Neurobiol Dis 2024;192:106426. [Crossref] [PubMed]
- Ko J, Park BS, Heo CM, Yi J, Lee DA, Park KM. Effect of glymphatic system function on cognitive function in patients with chronic kidney disease. Front Neurol 2024;15:1480536. [Crossref] [PubMed]
- Pépin M, Klimkowicz-Mrowiec A, Godefroy O, Delgado P, Carriazo S, Ferreira AC, et al. Cognitive disorders in patients with chronic kidney disease: Approaches to prevention and treatment. Eur J Neurol 2023;30:2899-911. [Crossref] [PubMed]
- Zhao H, Sun M, Zhang Y, Kong W, Fan L, Wang K, et al. Connecting the Dots: The Cerebral Lymphatic System as a Bridge Between the Central Nervous System and Peripheral System in Health and Disease. Aging Dis 2024;15:115-52. [Crossref] [PubMed]
- Zhang X, Wang Y, Jiao B, Wang Z, Shi J, Zhang Y, Bai X, Li Z, Li S, Bai R, Sui B. Glymphatic system impairment in Alzheimer's disease: associations with perivascular space volume and cognitive function. Eur Radiol 2024;34:1314-23. [Crossref] [PubMed]
- Uchida K. Waste Clearance in the Brain and Neuroinflammation: A Novel Perspective on Biomarker and Drug Target Discovery in Alzheimer's Disease. Cells 2022;11:919. [Crossref] [PubMed]
- Graterol Torres F, Molina M, Soler-Majoral J, Romero-González G, Rodríguez Chitiva N, Troya-Saborido M, Socias Rullan G, Burgos E, Paúl Martínez J, Urrutia Jou M, Cañameras C, Riera Sadurní J, Vila A, Bover J. Evolving Concepts on Inflammatory Biomarkers and Malnutrition in Chronic Kidney Disease. Nutrients 2022;14:4297. [Crossref] [PubMed]
- Taoka T, Ito R, Nakamichi R, Nakane T, Kawai H, Naganawa S. Diffusion Tensor Image Analysis ALong the Perivascular Space (DTI-ALPS): Revisiting the Meaning and Significance of the Method. Magn Reson Med Sci 2024;23:268-90. [Crossref] [PubMed]
- Wright AM, Wu YC, Chen NK, Wen Q. Exploring Radial Asymmetry in MR Diffusion Tensor Imaging and Its Impact on the Interpretation of Glymphatic Mechanisms. J Magn Reson Imaging 2024;60:1432-41. [Crossref] [PubMed]
- Drenthen GS, van der Thiel MM. Editorial for "Exploring Radial Asymmetry in MR Diffusion Tensor Imaging and Its Impact on the Interpretation of Glymphatic Mechanisms". J Magn Reson Imaging 2024;60:1442-3. [Crossref] [PubMed]
- Ulloa P, Rudolf JC, Kremer J, Schmidt A, Schramm P. Influence of orientation, size and shape of the region of interest in diffusion MRI along perivascular spaces index. MAGMA 2025;38:761-9. [Crossref] [PubMed]
(English Language Editor: J. Jones)


