Participants
The present study was approved by the Ethics Committee of the Second Hospital of Anhui Medical University. Written informed consent was provided from all participants. Fifty peritoneal dialysis ESRD patients (25 patients with SHPT, 14 males and 11 females, mean age 42.20±7.53 years (SHPT group); 25 patients without SHPT, 15 males and 10 females, mean age 41.96±6.17 years (Non-SHPT group)) were enrolled at the Nephrology Department in our hospital from January 2017 to July 2018. Inclusion criteria for patients as followings: (a) ESRD patients were diagnosed by (GRE (estimated glomerular filtration rate)≤15 mL/min/1.73 m2), (b) receiving regular peritoneal dialysis three times a week at our hospital, (c) SHPT for iPTH level > 600 pg/ml (9 times of normal upper level). 25 healthy controls (HC group, 12 males and 13 females, mean age 41.80 ± 7.15 years) were recruited. The healthy controls had no kidney disease and did not undergo laboratory tests.
All subjects were right-handed, and more than 18 years old. All participants received the neuropsychological test and MRI examination on the same day, receiving regular hemodialysis three times a week in our hospital. Exclusion criteria were as follows: (1) obvious encephalopathy (cerebrovascular disease) revealed by clinical or imaging, (2) any history of drug/alcohol abuse, (3) leukoencephalopathy on a previous MRI, and (4) MRI contraindications.
Neuropsychological test
Before the MRI scan, we used the Montreal Cognitive Assessment (MoCA) to assess the participants’ global cognitive abilities. A score less than 26 was defined as a diagnosis of cognitive impairment [20].
Laboratory test
Serum creatinine and urea levels within 24 hours prior to the MRI scan were conducted to evaluate renal function of all patients.
MRI Acquisition
MR imaging data were acquired in a Siemens Verio 3.0T MR scanner using an 8-channel coil. Use foam padding to reduce head motion. MRI sequence include 3D T1 weighted structural sagittal images (TE=2.98 ms, TR=1900 ms, thickness=1 mm, voxel size=1×1×1 mm3, FOV=256×256, 176 slices) and diffusion weighted images (30 directions, TE=84 ms, TR=8400 ms, slice thickness=3 mm, FOV=256×256, b=0, 1000 s/mm2, no gap) in all subjects. All subjects were told to keep their eyes closed and remain awake within the MRI scanning process.
Data Processing
Data preprocessing and network construction are accomplished by using PANDA toolbox (www.nitrc.org/projects/panda) [21]. Briefly, the image preprocessing procedure included format conversion from DICOM to NIFTI format, BET (skull removal), extract brain tissue and structure, eddy current and motion artifact correction, calculation of fractional anisotropy (FA) and diffusion tensor tractography. Fiber tracts were terminated if the FA value was lower than 0.2, or exceed an angular threshold of 45°.
WM network construction
Network node definition. The automated anatomical labeling (AAL) atlas [22] were used to parcel 90 cortical and subcortical regions (45 for each brain hemisphere) to define the network nodes [19] (Table 1, Figure 1).
Network edge definition. To define the edges of the structural network, we selected a threshold value for the fiber bundles with two end points located in both two regions. To remove spurious connections, we used a minimum of three fibers as a threshold [23]. We used fiber number (FN) of the connected fibers between two regions for the weights of the network edges. Finally, for each participant, the FN weighted 90×90 matrix structural networks were constructed [24] (Figure 1).
Network analysis
We performed by GRETNA (www.nitrc.org/projects/gretna) toolbox to analyze the WM network topological properties [25], included clustering coefficient (Cp), shortest path length (Lp), small-worldness (σ) (normalized clustering coefficient (γ) and normalized shortest path length (λ), local efficiency (Eloc), and global efficiency (Eg). The Cp of a network is the average of the clustering coefficient of each node in the network, which indicates the local efficiency for the transformation of the information. The Lp of the entire network refers to the average shortest travel distance across all nodes and indicates the most efficiency information transfer between the two nodes.
In order to test γ and λ in this study, 100 matched random networks are generated by Markov-chain Monte Carlo method. These networks have the same number of nodes, edges and degree distributions, but retain the weighted distribution of the real network. We computed the γ (γ=Cp/Cprandom) and λ (λ=Lp/Lprandom), where Cprandom and Lprandom are the average clustering coefficient and shortest path length over the random network, respectively. A small world network should meet the γ>1 and λ≈1.
Eg is the average of the inverse of the shortest path length of all node pairs in the network and usually reflects the ability of the network in parallel information processing. Eloc is the average of the global efficiency of the community neighboring all nodes in the network and represents the fault tolerance level of the network.
Identification of Hubs
Hubs of networks are key nodes that are estimated by different methods [26, 27]. Consistent with previous study [28], we identify the hub regions by regional efficiency. The hubs defined with regional efficiency at least one SD greater than the mean nodal efficiency of the network.
Statistical analysis
Demographic factors and clinical data were analyzed by SPSS 23.0 (SPSS, Chicago, IL, USA). We used mean ± SD for continuous variables, and frequencies for categorical variables in descriptive analysis. The differences among three groups including age, education and neuropsychological tests were analyzed with one-way analysis of variance (ANOVA). Post hoc pairwise comparisons were then performed using t tests. The sex data were analyzed by the chi-square test. In addition, two-sample t test was used to compare the clinical data between two groups for normality and homoscedasticity of the distributions, otherwise the Mann-Whitney U test was used.
For group effects in global and regional network metrics, comparisons were conducted among the three groups by one-way ANOVA and Nonparametric test (Kruskal-Wallis Test) for normal distribution and abnormal distribution, respectively. Levels of significance were set at p<0.05 and a Bonferroni correction for multiple testing was applied for group differences (α=0.05/3=0.017). We also detected the relationship between network metrics, neuropsychological test scores and the clinical index by partial correlation analysis. To investigate the correlation between the neuropsychological test scores and specific brain regions, the nodes with a significant group differences were taken.