Associations of the Disrupted Functional Brain Network and Cognitive Function in End-Stage Renal Disease Patients on Maintenance Hemodialysis: A Graph Theory-Based Study of Resting-State Functional Magnetic Resonance Imaging
Die Zhang
Yingying Chen1,3 
Jianlin Wu- 1Department of Radiology, Affiliated Zhongshan Hospital of Dalian University, Dalian, China
- 2Department of Radiology, Shenzhen Third People’s Hospital, Shenzhen, China
- 3Department of Radiology, National Cancer Center, National Clinical Research Center for Cancer, Cancer Hospital, Shenzhen Hospital, Chinese Academy of Medical Sciences, Peking Union Medical College, Shenzhen, China
- 4Department of Nephrology, Affiliated Zhongshan Hospital of Dalian University, Dalian, China
Objective: Cognitive impairment (CI) is a common neurological complication in patients with end-stage renal disease undergoing maintenance hemodialysis (MHD). Brain network analysis based on graph theory is a promising tool for studying CI. Therefore, the purpose of this study was to analyze the changes of functional brain networks in patients on MHD with and without CI by using graph theory and further explore the underlying neuropathological mechanism of CI in these patients.
Methods: A total of 39 patients on MHD (19 cases with CI and 20 without) and 25 healthy controls (HCs) matched for age, sex, and years of education were enrolled in the study. Resting-state functional magnetic resonance imaging (rs-fMRI) and T1-weighted high-resolution anatomical data were obtained, and functional brain networks for each subject were constructed. The brain network parameters at the global and regional levels were calculated, and a one-way analysis of covariance was used to compare the differences across the three groups. The associations between the changed graph-theory parameters and cognitive function scores in patients on MHD were evaluated using Spearman correlation analysis.
Results: Compared with HCs, the global parameters [sigma, gamma, and local efficiency (Eloc)] in both patient groups decreased significantly (p < 0.05, Bonferroni corrected). The clustering coefficient (Cp) in patients with CI was significantly lower than that in the other two groups (p < 0.05, Bonferroni corrected). The regional parameters were significantly lower in the right superior frontal gyrus, dorsolateral (SFGdor) and gyrus rectus (REC) of patients with CI than those of patients without CI; however the nodal local efficiency in the left amygdala was significantly increased (all p < 0.05, Bonferroni corrected). The global Cp and regional parameters in the three brain regions (right SFGdor, REC, and left amygdala) were significantly correlated with the cognitive function scores (all FDR q < 0.05).
Conclusion: This study confirmed that the topology of the functional brain network was disrupted in patients on MHD with and without CI and the disruption of brain network was more severe in patients with CI. The abnormal brain network parameters are closely related to cognitive function in patients on MHD.
Introduction
End-stage renal disease (ESRD) is defined as chronic kidney disease with a glomerular filtration rate below 15 ml/min//1.73 m2 or requiring permanent renal replacement therapy (Brouns and De Deyn, 2004). Maintenance hemodialysis (MHD) is one of the most important treatment methods.
Patients with ESRD are often affected by various neurological complications, such as cognitive impairment (CI). CI affects quality of life and mortality in patients with ESRD (Murray et al., 2006; Sarnak et al., 2013; Van Zwieten et al., 2019). In particular, patients with ESRD undergoing hemodialysis suffer from more severe CI (Murray et al., 2006). Moreover, a large number of unrecognized cognitive impairments in patients on hemodialysis have a significant impact on nutritional status, compliance with medical treatment, and hospitalization (Sehgal et al., 1997; Murray, 2008; Kwan et al., 2020). Therefore, exploration of the mechanism of hemodialysis-related cognitive decline in patients with ESRD is of great significance.
Neuroimaging is an important tool for studying cognitive decline in patients with ESRD (Kong et al., 2014; Chen et al., 2015; Zhang et al., 2015; Luo et al., 2016; Lu et al., 2019). These studies usually focus on changes in brain function and structure. Abnormal brain activity in regions related to cognitive function was found using resting-state functional magnetic resonance imaging (rs-fMRI) (Chen et al., 2015; Luo et al., 2016). Another study showed that altered functional connectivity of the frontal lobe is related to cognitive decline (Lu et al., 2019). Further, structural damage to the white matter is also associated with the performance of cognitive function (Kong et al., 2014; Zhang et al., 2015).
The human brain is a complex network system that can efficiently integrate and transfer information. Recently, a graph theory-based method revealed that the functional brain network in patients with ESRD was damaged, which may affect the cognitive function of patients (Park et al., 2020; Wu et al., 2020). However, the difference in brain network architecture between patients on MHD with and without CI remains unclear.
Therefore, in this study, we hypothesized that although the topology of functional brain networks in patients on MHD changes, these changes are not consistent between those with and without CI. To test our hypothesis, we used graph theory-based methods of brain network analysis to detect potential differences in the topology of functional brain networks between healthy controls (HCs) and patients on MHD with and without CI, and to explore the relationship between these possible alterations and cognitive function.
Materials and Methods
The ethics committee of our hospital approved this study, and informed consent was obtained from all the subjects.
Subjects
Between December 2016 and October 2018, 39 patients with ESRD on regular hemodialysis treatment (dialysis duration ≥ 3 months and 3–4 times a week) in our hospital’s hemodialysis center were enrolled in this study. The etiological factors of ESRD included: hypertension (12), chronic glomerulonephritis (10), diabetes (9), nephrotic syndrome (1), purpura nephritis (1), and unknown causes (6). The patients were divided into two groups: patients with CI [Montreal Cognitive Assessment Scale (MoCA) scores < 26] and those without CI (MoCA scores ≥ 26). Clinical data including age, sex, years of education, body mass index (BMI), and dialysis duration were collected from the electronic medical record system. During this cohort period, 25 healthy controls (HCs) who had no kidney disease, diabetes, hypertension, or other neuropsychiatric diseases were recruited from the local community.
All subjects were right-handed and the exclusion criteria were as follows: (1) age < 18 or >70 years; (2) obvious lesions (>1.0 cm) in the brain, such as cerebral infarction and brain tumor, confirmed by clinical history and conventional MRI; (3) severe heart failure and liver disease; (4) a history of drug or alcohol abuse; (5) contraindications to MRI scanning and inability to complete the neuropsychological test; and (6) translational motion greater than 2.5 mm, rotation greater than 2.5°, and mean framewise displacement values greater than 0.5 in data processing (Power et al., 2012).
Blood Biochemistry and Neuropsychological Test
Blood samples for biochemistry tests were collected up to 3 days before MRI scanning and results for blood urea, serum creatinine, uric acid, potassium (K), calcium (Ca), phosphorus (P), red blood cell (RBC), hemoglobin (HGB), and hematocrit (HCT) were checked. Blood biochemical parameters were not tested in the HCs.
To effectively perform the neuropsychological assessments for each subject, a neuropsychological test was conducted in a relatively quiet and comfortable environment and attempts were made to ensure that the results are not affected by the subjects’ uncoordinated behavior or emotional conflict. The test was commenced within one hour of performing the MRI scan. Three scales, including the MoCA (Drew et al., 2020), Trail Making Test A (TMT-A) (Llinàs-Reglà et al., 2017), and Digit Span Task (DST) (Anil et al., 2003), were administered to assess global cognitive function, executive function, and short-term memory function, respectively. Anxiety and depression were assessed using the Zung’s Self-rating Depression Scale (SDS) and Zung’s Self-rating Anxiety Scale (SAS), respectively.
Magnetic Resonance Imaging Data Acquisition
A 3.0-T MRI scanner equipped with a 12-channel head coil (Siemens, Verio, Germany) was used to obtain the imaging data. The participants were asked to lie quietly on the examination bed and the head was fixed with a foam sponge to reduce the impact of head movement. The scan sequences and parameters were as follows: (1) T2 fluid-attenuated inversion-recovery sequence (T2 FLAIR): repetition time (TR) = 4000 ms; echo time (TE) = 77 ms; flip angle (FA) = 150°; field of view (FOV) = 250 mm × 226 mm; slices = 20; and slice thickness = 5.0 mm. (2) T1-weighted three-dimensional (T1-3D) high-resolution: TR = 2530 ms; TE, 2.22 ms; FA, 7°; FOV, 224 × 224 mm; slices = 192; slice thickness = 1 mm. T2-FLAIR and T1-3D high-resolution images can be used to detect severe structural brain lesions (e.g., severe white matter lesions and tumors). (3) fMRI: TR = 2000 ms; TE = 30 ms; FG = 90°; FOV = 224 × 224 mm; slices = 31; slice thickness = 3.5 mm; 240 volumes were obtained for each participant. Before the scan of the functional series, all subjects were asked to remain still with their eyes closed and keep their mind clear without thinking about anything.
Data Preprocessing
Magnetic resonance imaging image processing tools (Data Processing & Analysis for Brain Imaging toolbox1 (Yan et al., 2016) based on MATLAB 2013b (The MathWorks, Inc., Natick, MA, United States)2 software were used to preprocess the MRI data. The preprocessing steps were as follows: (1) removal of the first 10 volumes; (2) slice timing and realignment of the remaining 230 volumes; (3) regression of nuisance covariates including white matter signal, cerebral spinal fluid signal, Friston-24 parameters of head motions, and linear detrend; (4) spatial normalization using the DARTEL algorithm (Chen et al., 2020); and (5) rmoval of temporal filtering (0.01–0.08 Hz).
Network Construction and Threshold Selection
The functional connection matrix for each subject was constructed using the GRETNA 2.0 package3 (Wang et al., 2015) in MATLAB 2013b. First, the whole brain was parcellated into 90 regions (nodes) based on automated anatomical labeling (AAL) atlas, which contains cortical and subcortical regions (Tzourio-Mazoyer et al., 2002). Next, the mean time series for each node was obtained, and the partial correlation coefficients of the mean time series between each node pair were calculated. Finally, a 90 × 90 functional connectivity matrix was acquired for each subject, and this matrix was converted into a binary matrix (i.e., adjacency matrices) based on a predefined threshold (described below), where the entry aij was 1, if the absolute value of the partial correlation coefficients between nodes i and j was greater than the predefined threshold. Otherwise, aij was 0.
Currently, there is no optimal or single strategy for threshold selection. Thus, a wide range of sparsity (S) was employed instead of the choice of a single threshold to meet a small-worldness as well as to perform sparse properties for each functional network (Rubinov and Sporns, 2010; Zhang et al., 2011). According to previous studies (Wu et al., 2020), the wide sparsity for brain network construction in this study ranges from 0.1 to 0.4, with an interval of 0.01.
Network Parameters
The brain network parameters at the global and regional levels were calculated at each sparsity (Rubinov and Sporns, 2010; Smith et al., 2011). In this study, global parameters consisted of the clustering coefficient (Cp), characteristic path length (Lp), normalized Cp (gamma), normalized Lp (lambda), small-worldness scalar (sigma), global efficiency (Eg), and local efficiency (Eloc). Regional parameters included nodal efficiency (Ne), nodal local efficiency (NLe), and nodal clustering coefficient (NCp).
In a network, clusters (including Cp, gamma, and NCp) and Eloc measured functional segregation (Rubinov and Sporns, 2010). Paths (including Lp and lambda) and Eg measured functional integration (Achard and Bullmore, 2007; Rubinov and Sporns, 2010). Sigma, which was defined as the ratio of gamma divided by lambda, was a scalar to quantify small-worldness. The three nodal parameters (Ne, NLe, and NCp) measured the local information processing and transmission at the regional level (details in the Supplementary Material).
Repeatability Analysis
The brain network parameters calculated according to different brain templates were not completely consistent (Wang et al., 2009; Fornito et al., 2010). The Harvard Oxford atlas (HOA) (Smith et al., 2004) with 112 brain regions and Brainnetome atlas with 246 brain regions (Fan et al., 2016) were selected to further evaluate the stability and repeatability of the results of this study.
Statistical Analysis
SPSS software (SPSS 24.0, IBM Corporation) was used to analyze the demographic and clinical data. For the patient group, blood biochemistry and dialysis duration were analyzed using a two-sample t-test or Mann-Whitney U test. For other clinical data, a one-way analysis of variance (ANOVA) or the Kruskal-Wallis test was used to detect differences among the three groups. Categorical variables were compared using the Pearson’s χ2 test.
The area under the curve (AUC), as a summarized scalar, is a highly sensitive scalar independent of a single threshold selection (Lei et al., 2015; Suo et al., 2015). Therefore, a one-way analysis of covariance (ANCOVA) based on the AUC was used to determine the inter-group differences in the brain network measures. The possible confounding factors included age, sex, years of education, and SDS scores. The problem of multiple comparisons at the nodal level was addressed using the Benjamini-Hochberg false discovery rate correction method (FDR q < 0.05) (Genovese et al., 2002; Chen et al., 2018).
The Spearman correlation analysis was used to explore the relationship between the changed network parameters and cognitive performance in patients on MHD, with age, sex, years of education, and SDS scores as covariates.
The significance level of this study was set as 0.05, two-tailed, and post hoc tests (p < 0.05, Bonferroni corrected) were further used for pairwise comparison (Genovese et al., 2002; Chen et al., 2018).
Results
Demographic and Clinical Characteristics
Demographic and clinical characteristics are summarized in Table 1. MoCA and DST scores were significantly lower in patients with CI than in patients without CI and HCs (p < 0.05, Bonferroni corrected). The time of TMT-A in patients with CI was significantly longer than that in the other two groups (p < 0.05, Bonferroni corrected). SDS scores significantly increased in patients with CI, compared with HCs (p < 0.05, Bonferroni corrected).
Table 1. Participant demographic and clinical characteristics.
Alterations in Global Parameters
All participants demonstrated a small-worldness, with gamma higher than 1 and lambda approximately equal to 1 (Figure 1). According to the post hoc test (Table 2), sigma, gamma, and Eloc significantly decreased in both patient groups compared with those in HCs (p < 0.05, Bonferroni corrected). Cp significantly decreased in patients with CI compared with that in the other two groups (p < 0.05, Bonferroni corrected) (Figure 2).
Figure 1. The three groups have small-worldness in a range of sparsity (0.1–0.4), i.e., gamma > 1 (A), lambda ≈ 1 (B). HCs, healthy controls; CI, cognitive impairment.
Table 2. Area under the curve (AUC) under all thresholds at the global level.
Figure 2. Global parameters with differences among the three groups. (A–D) Clustering coefficient (Cp), normalized Cp (Gamma), small-worldness scalar (Sigma), and local efficiency (Eloc), respectively. * indicates p < 0.05, ** indicates p < 0.01, *** indicates p < 0.001 (Bonferroni corrected). HCs, healthy controls; CI, cognitive impairment.
Alterations in Regional Parameters
Analysis of covariances revealed that there was no significant difference in Ne among the three groups (FDR q > 0.05); however, several brain regions with significant alterations in NLe and NCp were observed in the frontal and temporal lobes (FDR q < 0.05) (Table 3 and Figures 3, 4).
Table 3. Brain regions with significant group differences in regional characteristics across the three groups.
Figure 3. This figure shows significant differences in nodal local efficiency (NLe) among the three groups. Dark blue indicates the nodes with differences among the three groups using ANCOVA. Red represents a significant increase and green represents a significant decrease (p < 0.05, Bonferroni corrected). SFGdor, superior frontal gyrus dorsolateral; REC, gyrus rectus; DCG, median cingulate and paracingulate gyri; AMYG, amygdala; HES, Heschl’s gyrus.
Figure 4. This figure shows significant differences in the nodal clustering coefficient (NCp) among the three groups. Dark blue indicates the nodes with differences among the three groups using ANCOVA. Red represents a significant increase and green represents a significant decrease (p < 0.05, Bonferroni corrected). SFGdor, superior frontal gyrus dorsolateral; REC, gyrus rectus; DCG, median cingulate and paracingulate gyri; AMYG, amygdala; HES, Heschl’s gyrus.
Compared with HCs, regional parameters of the right superior frontal gyrus, dorsolateral (SFGdor), and bilateral median cingulate and paracingulate gyri (DCG) increased significantly in patients without CI, and significantly decreased in the left amygdala (AMYG) and right Heschel (HES) (p < 0.05, Bonferroni corrected).
Compared with HCs, the regional parameters of the right DCG were significantly increased in patients with CI, while the right gyrus rectus (REC) and HES were significantly decreased (p < 0.05, Bonferroni corrected).
Relative to patients without CI, significantly increased and decreased nodal parameters were located in the left AMYG and right frontal lobes (right SFGdor and REC), respectively (p < 0.05, Bonferroni corrected).
Correlation Analysis
The detailed results of the correlation analysis are shown in Table 4. Cp was associated with the MoCA score (r = 0.505; FDR q = 0.004). For the right SFGdor, the MoCA score was positively correlated with the NLe (r = 0.403; FDR q = 0.03) and NCp (r = 0.425; FDR q = 0.021). For the right REC, the MoCA, DST-B, and DST-T scores were positively associated with the NLe (r = 0.4, FDR q = 0.03; r = 0.494, FDR q = 0.005; and r = 0.49, FDR q = 0.01, respectively) and NCp (r = 0.448, FDR q = 0.021; r = 0.485, FDR q = 0.012; and r = 0.446, FDR q = 0.024, respectively). For the left AMYG, the MoCA score was negatively correlated with the NCp (r = –0.378; FDR q = 0.036), and the TMT-A time was positively correlated with both the NLe (r = 0.417; FDR q = 0.04) and NCp (r = 0.477; FDR q = 0.012).
Table 4. Correlation analysis of changed brain network parameters and cognitive function in patients on MHD.
Results of Repeatability
The global parameters using the HOA and Brainnetome atlas were similar to the results from AAL templates (Supplementary Table 1). According to the HOA and Brainnetome atlas, significant changes in nodal parameters were not only found in the frontal and temporal lobes but were also found in other brain regions, such as the insula and thalamus (Supplementary Table 2).
Discussion
This study, based on rs-fMRI, revealed abnormal changes in functional brain networks in patients on MHD by using graph theory analysis, and the main findings are as follows: (1) in comparison with HCs, global properties, including sigma, gamma, and Eloc, were significantly decreased in patients with and without CI. Cp was significantly lower in patients with CI than in the other two groups; (2) in comparison with HCs, the altered nodal parameters (NCp and NLe) of patients on MHD were mainly located in the frontal and temporal lobes; (3) in comparison with patients without CI, the nodal parameters (NCp and NLe) of the prefrontal cortex (SFGdor and REC) were significantly decreased and NLe of the amygdala was significantly increased in patients with CI; and (4) several altered global and regional parameters of patients on MHD were significantly correlated with cognition scores. These results contribute to the understanding of the underlying neuropathological mechanism of cognitive decline in patients on MHD.
Altered Global Parameters
The human brain network is a network model with small-worldness (Zhang et al., 2011; Suo et al., 2015), which combines a higher Cp and shorter Lp to satisfy the opposing demands of global and local processing and ensure highly efficient transformation and processing of information at the global and regional levels (Kaiser and Hilgetag, 2006; Bullmore and Sporns, 2009). The more the brain network structure deviates from the optimal state, the more prone it is to abnormal brain function (Rubinov and Sporns, 2010; Xu et al., 2019). The results of this study demonstrated that the three groups had small-worldness (Figure 1); however, the small-worldness scalar (sigma) of both patient group was significantly lower than that of HCs. This reflects the fact that the overall structure of the small-world network does not reach the optimal state, which further affects the transmission and processing of whole brain information. A study (Wu et al., 2020) revealed that compared with HCs, the sigma of patients on hemodialysis was also significantly reduced, but the reduced Lp and lambda in patients was not consistent with our results. The cognitive function of patients with ESRD is affected by dialysis duration (Brouhard et al., 2000), and the brain network architecture is related to cognitive level (Liu et al., 2012). Given that the duration of hemodialysis among patients with ESRD in this study was longer than that in previous studies (Wu et al., 2020; Yue et al., 2021), it is postulated that these inconsistent results may be due to hemodialysis duration. According to the definition of Eloc and gamma (Sporns and Zwi, 2004; Rubinov and Sporns, 2010), their decreased values in the patient group reflect the decreased information transmission ability of patients on MHD at the regional level. In addition, these results further confirm that hemodialysis can adversely affect the brain networks of patients (Chou et al., 2019; Wu et al., 2020).
At the global level, an important finding was that the Cp of patients with CI decreased significantly compared with patients without CI (p < 0.05, Bonferroni corrected), which implies that the disruption of the brain network was more serious in patients with CI. Cp measures the extent of local cliquishness of the network, i.e., functional segregation (Rubinov and Sporns, 2010). A decrease in the Cp reflects a decrease in the information transmission efficiency at the regional level. In this study, patients with CI had significantly reduced Cp compared with patients without CI, which implies local information processing was significantly reduced in those patients. Previous studies (Park et al., 2020; Xiong et al., 2020) have confirmed a close relationship between Cp and cognitive function. Moreover, the correlation analysis in this study revealed Cp was significantly and positively correlated with the global cognitive function (MoCA) in patients on MHD. Therefore, the reduction of local information processing may be an important neuropathological mechanism for the decline of global cognitive function in patients on MHD.
Altered Regional Parameters
Consistent with the changes of the Cp and Eloc, NCp, and NLe were the significantly altered nodal parameters. They mainly affected the frontal and temporal lobes. The frontal and temporal lobes play an important role in several cognitive functions, such as executive function, attention, and memory (Chayer and Freedman, 2001; Oh and Jagust, 2013; Fang et al., 2016). Consequently, the abnormal node attributes in the frontal-temporal lobe may be associated with potential cognitive decline in patients on MHD.
Further analysis found that, except for the amygdala, the nodal parameters of the frontal and temporal cortices in patients with CI were significantly or slightly lower than those in patients without CI. This indicated that the more severe brain network damage in patients on MHD with CI mainly occurred in the frontal-temporal lobe, particularly in the prefrontal lobe [NCp and NLe of the prefrontal lobe in patients with CI were significantly lower than those in patients without CI (p < 0.05, Bonferroni corrected)]. A study confirmed that the impairment of the brain network in the prefrontal lobe in patients with ESRD with CI was more serious than in patients without CI (Yue et al., 2021). However, unlike in the previous study, patients on MHD with CI in our study did not show a significant reduction in node attributes in the precuneus and posterior cingulate gyrus. This may have been due to the use of different nodal parameters for evaluation of the changes in node attributes; a more important reason being the different hemodialysis duration, which was similar to the difference at the global level. In addition, compared with HCs, the nodal parameters of the frontal lobe increased and those of the temporal lobe decreased in patients on MHD without CI. This result is similar to the conclusion of a study by Wu et al. (2020) that suggested that the nodal degree of the frontal cortex in patients on hemodialysis is higher than that in HCs, and the nodal degree of the temporal lobe is lower than that in HCs. The increase in degree centrality in multiple brain regions, including the middle frontal gyrus, may be an important compensatory mechanism for patients with Alzheimer’s disease to delay dementia (Behfar et al., 2020). Similarly, the enhancement of frontal cortex activity is also a compensatory mechanism for healthy elderly people to resist neurodegeneration and maintain normal cognitive levels (Scheller et al., 2018). Therefore, the increased nodal parameters of the frontal lobes in patients on MHD without CI in this study may reflect an adaptive compensatory mechanism of delaying cognitive decline in patients on MHD. Taken together, these results suggest that the occurrence of CI in patients on MHD may be related to the impairment of the frontotemporal brain network and the weakening of frontal nerve compensation.
The nodal parameters of the amygdala (including NCp and NLe) in patients on MHD without CI were significantly lower than those in HCs, which is consistent with the results of the previous studies (Wu et al., 2020; Yue et al., 2021). The amygdala is an important structure for the regulation and processing of emotions (Dalton et al., 2005). Decreased functional connectivity between the amygdala and multiple brain regions may be related to depression in patients with ESRD (Chen et al., 2017). However, further analysis revealed that the NLe of the amygdala in patients with CI was significantly higher than that in patients without CI, even if SDS scores were taken as a covariate. This indicated that abnormal brain connectivity in the amygdala of patients with CI may not be related to abnormal emotion regulation only. We speculate that the abnormal brain connectivity of the amygdala in patients with CI may be the result of interaction between cognition and emotion regulation since there is an interaction effect between emotion and cognitive function (Brosch et al., 2013).
Relationship Between Altered Network Parameter and Cognitive Function in Patients on Maintenance Hemodialysis
Clustering coefficient was significantly and positively correlated with the MoCA score. Although more notably, the nodal parameters of the brain regions (right SFGdor and REC, and left AMYG), with significant differences among the patient groups, were significantly correlated with the cognitive performance of patients on MHD. Specifically, both the NCp and NLe of the prefrontal lobe (right SFGdor and REC) were significantly and positively correlated with the MoCA score. In patients with ESRD, the prefrontal lobe is an important brain region related to cognition, and functional and structural abnormalities of the prefrontal lobe have been observed in previous studies (Ma et al., 2015; Jin et al., 2020). Other studies have confirmed that for patients with ESRD, prefrontal connectivity is significantly correlated with global cognitive function (i.e., the MoCA score) (Wu et al., 2020; Huang et al., 2021). These results suggest that weakening of local information transmission, especially in the prefrontal lobe, may play an important role in the global cognitive decline of patients on MHD and may become an imaging marker for early evaluation of cognitive decline in such patients. Additionally, there was a significant positive correlation between the nodal parameters (NCp and NLe) of the right REC and short-term memory function scores (DST-B and DST-T scores). The REC plays an important role in memory function; however, both memory and language function would be negatively affected should the REC be resected (Joo et al., 2016). The REC is also an important part of the orbitofrontal cortex, which is directly involved in memory processing (Petrides, 2007). Accordingly, the functional disconnect of the REC may be involved, at least in part, in the disruption of short-term memory function in patients on MHD. The nodal parameters of the left AMYG were negatively correlated with the MoCA score and positively correlated with TMT-A time (a test for executive function). The AMYG can participate in and coordinate executive attention through the amygdala–prefrontal circuit (Salzman and Fusi, 2010). Consequently, the abnormal node attribute of the AMYG in this study may be related to cognitive deficits in patients on MHD, especially executive function.
Repeatability Analysis
The global and regional attributes of the three different anatomical atlases have high similarity, which indicates that the results of this study have good repeatability and stability. Particularly at the global level, ANCOVA results showed that the three atlases obtained more consistent results.
Limitations
This study had some limitations. First, the sample size of this study was small. The cross-sectional study design of this study may make it difficult to reveal the longitudinal changes in functional brain networks and cognitive function in patients on MHD. This needs to be further verified by increasing the sample size and conducting longitudinal studies in the future. Second, the evaluation of cognitive function in this study lacks a systematic approach and objectivity, which may have resulted in some bias. Third, the lack of detailed grouping in this study to distinguish the primary and secondary diseases of ESRD may also affect the results.
Conclusion
This study confirmed that the topological architecture of the functional brain network in patients on MHD was disrupted. Furthermore, the impairment of the brain network in patients with CI was more severe than that in patients without CI, particularly in the prefrontal cortex. Changes in brain network parameters at global and regional levels are related to the cognitive function of patients on MHD. Functional brain network analysis based on graph theory can provide an important imaging basis for the neuropathological mechanism of cognitive dysfunction in patients on MHD.
Data Availability Statement
The data analyzed in this study is subject to the following licenses/restrictions: Raw data will not be made public to protect the subjects’ privacy. Requests to access these datasets should be directed to JW, cjr.wujianlin@vip.163.com.
Ethics Statement
The studies involving human participants were reviewed and approved by Ethics Committee of the Affiliated Zhongshan Hospital of Dalian University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JW and DZ contributed to the conception and design of the study. DZ, YC, and LL performed the statistical analysis and wrote the manuscript. HW, QX, CC, and LJ performed the experiments. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
The authors declare that this study was received funding from the National Natural Science Foundation of China (grant number 81371526).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank all patients and volunteers for their participation.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnhum.2021.716719/full#supplementary-material
Footnotes
References
Achard, S., and Bullmore, E. (2007). Efficiency and cost of economical brain functional networks. PLoS Comput. Biol. 3:e17. doi: 10.1371/journal.pcbi.0030017
Anil, A. E., Kivircik, B. B., Batur, S., Kabakçi, E., Kitiş, A., Güven, E., et al. (2003). The Turkish version of the auditory consonant trigram test as a measure of working memory: a normative study. Clin. Neuropsychol. 17, 159–169. doi: 10.1076/clin.17.2.159.16510
Behfar, Q., Behfar, S. K., von Reutern, B., Richter, N., Dronse, J., Fassbender, R., et al. (2020). Graph theory analysis reveals resting-state compensatory mechanisms in healthy aging and prodromal Alzheimer’s disease. Front. Aging Neurosci. 12:576627. doi: 10.3389/fnagi.2020.576627
Brosch, T., Scherer, K. R., Grandjean, D., and Sander, D. (2013). The impact of emotion on perception, attention, memory, and decision-making. Swiss Med. Wkly. 143:w13786. doi: 10.4414/smw.2013.13786
Brouhard, B. H., Donaldson, L. A., Lawry, K. W., McGowan, K. R., Drotar, D., Davis, I., et al. (2000). Cognitive functioning in children on dialysis and post-transplantation. Pediatr. Transplant. 4, 261–267. doi: 10.1034/j.1399-3046.2000.00121.x
Brouns, R., and De Deyn, P. P. (2004). Neurological complications in renal failure: a review. Clin. Neurol. Neurosurg. 107, 1–16. doi: 10.1016/j.clineuro.2004.07.012
Bullmore, E., and Sporns, O. (2009). Complex brain networks: graph theoretical analysis of structural and functional systems. Nat. Rev. Neurosci. 10, 186–198. doi: 10.1038/nrn2575
Chayer, C., and Freedman, M. (2001). Frontal lobe functions. Curr. Neurol. Neurosci. Rep. 1, 547–552. doi: 10.1007/s11910-001-0060-4
Chen, H. J., Qi, R., Kong, X., Wen, J., Liang, X., Zhang, Z., et al. (2015). The impact of hemodialysis on cognitive dysfunction in patients with end-stage renal disease: a resting-state functional MRI study. Metab. Brain Dis. 30, 1247–1256. doi: 10.1007/s11011-015-9702-0
Chen, H. J., Qiu, J., Fu, Q., and Chen, F. (2020). Alterations of spontaneous brain activity in hemodialysis patients. Front. Hum. Neurosci. 14:278. doi: 10.3389/fnhum.2020.00278
Chen, H. J., Shi, H. B., Jiang, L. F., Li, L., and Chen, R. (2018). Disrupted topological organization of brain structural network associated with prior overt hepatic encephalopathy in cirrhotic patients. Eur. Radiol. 28, 85–95. doi: 10.1007/s00330-017-4887-8
Chen, H. J., Wang, Y. F., Qi, R., Schoepf, U. J., Varga-Szemes, A., Ball, B. D., et al. (2017). Altered amygdala resting-state functional connectivity in maintenance hemodialysis end-stage renal disease patients with depressive mood. Mol. Neurobiol. 54, 2223–2233. doi: 10.1007/s12035-016-9811-8
Chou, M. C., Ko, C. H., Chang, J. M., and Hsieh, T. J. (2019). Disruptions of brain structural network in end-stage renal disease patients with long-term hemodialysis and normal-appearing brain tissues. J. Neuroradiol. 46, 256–262. doi: 10.1016/j.neurad.2018.04.004
Dalton, K. M., Nacewicz, B. M., Johnstone, T., Schaefer, H. S., Gernsbacher, M. A., Goldsmith, H. H., et al. (2005). Gaze fixation and the neural circuitry of face processing in autism. Nat. Neurosci. 8, 519–526. doi: 10.1038/nn1421
Drew, D. A., Tighiouart, H., Rollins, J., Duncan, S., Babroudi, S., Scott, T., et al. (2020). Evaluation of screening tests for cognitive impairment in patients receiving maintenance hemodialysis. J. Am. Soc. Nephrol. 31, 855–864. doi: 10.1681/ASN.2019100988
Fan, L., Li, H., Zhuo, J., Zhang, Y., Wang, J., Chen, L., et al. (2016). The human Brainnetome atlas: a new brain atlas based on connectional architecture. Cereb. Cortex 26, 3508–3526. doi: 10.1093/cercor/bhw157
Fang, S., Wang, Y., and Jiang, T. (2016). The influence of frontal lobe tumors and surgical treatment on advanced cognitive functions. World Neurosurg. 91, 340–346. doi: 10.1016/j.wneu.2016.04.006
Fornito, A., Zalesky, A., and Bullmore, E. T. (2010). Network scaling effects in graph analytic studies of human resting-state FMRI data. Front. Syst. Neurosci. 4:22. doi: 10.3389/fnsys.2010.00022
Genovese, C. R., Lazar, N. A., and Nichols, T. (2002). Thresholding of statistical maps in functional neuroimaging using the false discovery rate. Neuroimage 15, 870–878. doi: 10.1006/nimg.2001.1037
Huang, J., Xie, L., Guo, R., Wang, J., Lin, J., Sun, Z., et al. (2021). Abnormal brain activity patterns during spatial working memory task in patients with end-stage renal disease on maintenance hemodialysis: a fMRI study. Brain Imaging Behav. 15, 1898–1911. doi: 10.1007/s11682-020-00383-7
Jin, M., Wang, L., Wang, H., Han, X., Diao, Z., Guo, W., et al. (2020). Structural and functional alterations in hemodialysis patients: a voxel-based morphometry and functional connectivity study. Front. Hum. Neurosci. 14:80. doi: 10.3389/fnhum.2020.00080
Joo, M. S., Park, D. S., Moon, C. T., Chun, Y. I., Song, S. W., and Roh, H. G. (2016). Relationship between gyrus rectus resection and cognitive impairment after surgery for ruptured anterior communicating artery aneurysms. J. Cerebrovasc. Endovasc. Neurosurg. 18, 223–228. doi: 10.7461/jcen.2016.18.3.223
Kaiser, M., and Hilgetag, C. C. (2006). Nonoptimal component placement, but short processing paths, due to long-distance projections in neural systems. PLoS Comput. Biol. 2:e95. doi: 10.1371/journal.pcbi.0020095
Kong, X., Wen, J. Q., Qi, R. F., Luo, S., Zhong, J. H., Chen, H. J., et al. (2014). Diffuse interstitial brain edema in patients with end-stage renal disease undergoing hemodialysis: a tract-based spatial statistics study. Medicine 93:e313. doi: 10.1097/MD.0000000000000313
Kwan, E., Draper, B., Endre, Z. H., Harvey, S. B., and Brown, M. A. (2020). Prevalence, types and recognition of cognitive impairment in dialysis patients in south eastern Sydney. Intern. Med. J. doi: 10.1111/imj.14976 [Epub ahead of print].
Lei, D., Li, K., Li, L., Chen, F., Huang, X., Lui, S., et al. (2015). Disrupted functional brain connectome in patients with posttraumatic stress disorder. Radiology 276, 818–827. doi: 10.1148/radiol.15141700
Liu, Z., Zhang, Y., Yan, H., Bai, L., Dai, R., Wei, W., et al. (2012). Altered topological patterns of brain networks in mild cognitive impairment and Alzheimer’s disease: a resting-state fMRI study. Psychiatry Res. 202, 118–125. doi: 10.1016/j.pscychresns.2012.03.002
Llinàs-Reglà, J., Vilalta-Franch, J., López-Pousa, S., Calvó-Perxas, L., Torrents Rodas, D. T., and Garre-Olmo, J. (2017). The trail making test. Assessment 24, 183–196. doi: 10.1177/1073191115602552
Lu, H., Gu, Z., Xing, W., Han, S., Wu, J., Zhou, H., et al. (2019). Alterations of default mode functional connectivity in individuals with end-stage renal disease and mild cognitive impairment. BMC Nephrol. 20:246. doi: 10.1186/s12882-019-1435-6
Luo, S., Qi, R. F., Wen, J. Q., Zhong, J. H., Kong, X., Liang, X., et al. (2016). Abnormal intrinsic brain activity patterns in patients with end-stage renal disease undergoing peritoneal dialysis: a resting-state functional MR imaging study. Radiology 278, 181–189. doi: 10.1148/radiol.2015141913
Ma, X., Jiang, G., Li, S., Wang, J., Zhan, W., Zeng, S., et al. (2015). Aberrant functional connectome in neurologically asymptomatic patients with end-stage renal disease. PLoS One 10:e0121085. doi: 10.1371/journal.pone.0121085
Murray, A. M. (2008). Cognitive impairment in the aging dialysis and chronic kidney disease populations: an occult burden. Adv. Chronic Kidney Dis. 15, 123–132. doi: 10.1053/j.ackd.2008.01.010
Murray, A. M., Tupper, D. E., Knopman, D. S., Gilbertson, D. T., Pederson, S. L., Li, S., et al. (2006). Cognitive impairment in hemodialysis patients is common. Neurology 67, 216–223. doi: 10.1212/01.wnl.0000225182.15532.40
Oh, H., and Jagust, W. J. (2013). Frontotemporal network connectivity during memory encoding is increased with aging and disrupted by beta-amyloid. J. Neurosci. 33, 18425–18437. doi: 10.1523/JNEUROSCI.2775-13.2013
Park, B. S., Kim, S. E., Lee, H. J., Kim, Y. W., Kim, I. H., Park, J. H., et al. (2020). Alterations in structural and functional connectivities in patients with end-stage renal disease. J. Clin. Neurol. 16, 390–400. doi: 10.3988/jcn.2020.16.3.390
Petrides, M. (2007). The orbitofrontal cortex: novelty, deviation from expectation, and memory. Ann. N. Y. Acad. Sci. 1121, 33–53. doi: 10.1196/annals.1401.035
Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L., and Peterson, S. E. (2012). Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59, 2142–2154. doi: 10.1016/j.neuroimage.2011.10.018
Rubinov, M., and Sporns, O. (2010). Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52, 1059–1069. doi: 10.1016/j.neuroimage.2009.10.003
Salzman, C. D., and Fusi, S. (2010). Emotion, cognition, and mental state representation in amygdala and prefrontal cortex. Annu. Rev. Neurosci. 33, 173–202. doi: 10.1146/annurev.neuro.051508.135256
Sarnak, M. J., Tighiouart, H., Scott, T. M., Lou, K. V., Sorensen, E. P., Giang, L. M., et al. (2013). Frequency of and risk factors for poor cognitive performance in hemodialysis patients. Neurology 80, 471–480. doi: 10.1212/WNL.0b013e31827f0f7f
Scheller, E., Schumacher, L. V., Peter, J., Lahr, J., Wehrle, J., Kaller, C. P., et al. (2018). Brain aging and APOE ε4 interact to reveal potential neuronal compensation in healthy older adults. Front. Aging Neurosci. 10:74. doi: 10.3389/fnagi.2018.00074
Sehgal, A. R., Grey, S. F., Deoreo, P. B., and Whitehouse, P. J. (1997). Prevalence, recognition, and implications of mental impairment among hemodialysis patients. Am. J. Kidney Dis. 30, 41–49. doi: 10.1016/S0272-6386(97)90563-1
Smith, S. M., Jenkinson, M., Woolrich, M. W., Beckmann, C. F., Behrens, T. E., Johansen-Berg, H., et al. (2004). Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23, S208–S219. doi: 10.1016/j.neuroimage.2004.07.051
Smith, S. M., Miller, K. L., Salimi-Khorshidi, G., Webster, M., Beckmann, C. F., Nichols, T. E., et al. (2011). Network modelling methods for FMRI. Neuroimage 54, 875–891. doi: 10.1016/j.neuroimage.2010.08.063
Sporns, O., and Zwi, J. D. (2004). The small world of the cerebral cortex. Neuroinformatics 2, 145–162.
Suo, X., Lei, D., Li, K., Chen, F., Li, F., Li, L., et al. (2015). Disrupted brain network topology in pediatric posttraumatic stress disorder: a resting-state fMRI study. Hum. Brain Mapp. 36, 3677–3686. doi: 10.1002/hbm.22871
Tzourio-Mazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., et al. (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 15, 273–289. doi: 10.1006/nimg.2001.0978
Van Zwieten, A., Wong, G., Ruospo, M., Palmer, S. C., Teixeira-Pinto, A., Barulli, M. R., et al. (2019). Associations of cognitive function and education level with all-cause mortality in adults on hemodialysis: findings from the COGNITIVE-HD study. Am. J. Kidney Dis. 74, 452–462. doi: 10.1053/j.ajkd.2019.03.424
Wang, J., Wang, L., Zang, Y., Yang, H., Tang, H., Gong, Q., et al. (2009). Parcellation-dependent small-world brain functional networks: a resting-state fMRI study. Hum. Brain Mapp. 30, 1511–1523. doi: 10.1002/hbm.20623
Wang, J., Wang, X., Xia, M., Liao, X., Evans, A., and He, Y. (2015). Gretna: a graph theoretical network analysis toolbox for imaging connectomics. Front. Hum. Neurosci. 9:386. doi: 10.3389/fnhum.2015.00386
Wu, B., Li, X., Zhang, M., Zhang, F., Long, X., Gong, Q., et al. (2020). Disrupted brain functional networks in patients with end-stage renal disease undergoing hemodialysis. J. Neurosci. Res. 98, 2566–2578. doi: 10.1002/jnr.24725
Xiong, Y., Chen, X., Zhao, X., Fan, Y., Zhang, Q., and Zhu, W. (2020). Altered regional homogeneity and functional brain networks in Type 2 diabetes with and without mild cognitive impairment. Sci. Rep. 10:21254. doi: 10.1038/s41598-020-76495-3
Xu, J., Chen, F., Liu, T., Wang, T., Zhang, J., Yuan, H., et al. (2019). Brain functional networks in type 2 diabetes mellitus patients: a resting-state functional MRI study. Front. Neurosci. 13:239. doi: 10.3389/fnins.2019.00239
Yan, C. G., Wang, X. D., Zuo, X. N., and Zang, Y. F. (2016). DPABI: data processing & analysis for (resting-state) brain imaging. Neuroinformatics 14, 339–351. doi: 10.1007/s12021-016-9299-4
Yue, Z., Wang, P., Li, X., Ren, J., and Wu, B. (2021). Abnormal brain functional networks in end-stage renal disease patients with cognitive impairment. Brain Behav. 11:e02076. doi: 10.1002/brb3.2076
Zhang, J., Wang, J., Wu, Q., Kuang, W., Huang, X., He, Y., et al. (2011). Disrupted brain connectivity networks in drug-naive, first-episode major depressive disorder. Biol. Psychiatry 70, 334–342. doi: 10.1016/j.biopsych.2011.05.018
Keywords: end-stage renal disease, maintenance hemodialysis, graph theory, resting-state functional magnetic resonance imaging (rs-fMRI), cognitive impairment
Citation: Zhang D, Chen Y, Wu H, Lin L, Xie Q, Chen C, Jing L and Wu J (2021) Associations of the Disrupted Functional Brain Network and Cognitive Function in End-Stage Renal Disease Patients on Maintenance Hemodialysis: A Graph Theory-Based Study of Resting-State Functional Magnetic Resonance Imaging. Front. Hum. Neurosci. 15:716719. doi: 10.3389/fnhum.2021.716719
Received: 01 June 2021; Accepted: 22 November 2021;
Published: 13 December 2021.
Edited by:
Raffaella Franciotti, G. d’Annunzio University of Chieti-Pescara, ItalyReviewed by:
Jie Lu, Capital Medical University, ChinaShao-Wei Xue, Affiliated Hospital of Hangzhou Normal University, China
Copyright © 2021 Zhang, Chen, Wu, Lin, Xie, Chen, Jing and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianlin Wu, cjr.wujianlin@vip.163.com