Participants
This study has been approved by the Ethics Committee of Guangzhou First People's Hospital (the Second Affiliated Hospital of South China University of Technology). All subjects provided informed consent. The study was carried out in accordance with the principles of the Declaration of Helsinki and its later amendments. PD subjects were all patients admitted to the Department of Neurology of Guangzhou First People's Hospital from March 2020 to August 2024. The inclusion criteria were as follows: (1) PD was diagnosed for the first time; (2) No drug intervention has been carried out; (3) Completed imaging data collection and clinical evaluation. The exclusion criteria were as follows: (1) Poor image quality was not conducive to the analyst; (2) Patients combined with other nervous system disorders. Age - and sex-matched healthy individuals recruited from the community as healthy controls (HCs) (Fig. 1).
Following the inclusion and exclusion criteria, 39 HC and 122 PD were included. According to MMSE scores, PD patients were divided into two groups: MCI and NC.
Clinical assessments
Baseline clinical features included: Baseline age, sex, duration of disease from PD diagnosis to enrollment, Hoehn & Yahr (H-Y) staging, MDS-UPDRS sub-scores (UPDRS-I, UPDRS-II, UPDRS-III, and UPDRS-IV), and Mini-Mental State Examination (MMSE) score. H-Y grading method is used to record the severity of the disease for PD patients24. The first part of UPDRS is used to assess the degree of mental activity, behavior and emotional disorders of patients. The second part of UPDRS is used to assess patients' daily living ability. The third part of UPDRS is used to evaluate patients' motor function; and the fourth part of UPDRS is used to record the complications of patients after treatment25. Participants' cognitive function was scored using MMSE. The lower the MMSE score, the more serious the cognitive dysfunction26. According to the score, mild cognitive impairment was defined as an MMSE score of 21-26, and cognitive normal was defined as an MMSE score of 27 or above.
Imaging data acquisition and preprocessing
- MRI data acquisition
MRIs were acquired on a 3 Tesla Siemens system (Erlangen, Germany) equipped with a 16-channel phased-array head coil. 3D T1-weighted MPRAGE sequence was acquired using the following specifications: echo time (TE) = 2.19 ms, repetition time (TR) = 1,900 ms, field of view (FOV) = 250 mm x 250 mm, flip angle = 9°, slice thickness = 1.0 mm and slice gap = 0 mm. Other sequences included rs-fMRI acquisition with a single-shot echo-planar imaging sequence (TE = 21 ms, TR = 2,000 ms, FOV = 224 mm x 224 mm, flip angle = 78°, matrix = 64 x 64, slice thickness = 4.0 mm, voxel size = 3.5 mm x 3.5 mm x 4.0 mm, and comprising 220 time points) and DTI with a spin echo planar imaging sequence (TE = 102 ms, TR = 8700 ms, FOV = 230 mm x 230 mm², number of excitations (NEX) = 1, slice thickness = 2.5 mm, and matrix = 92 x 92).
- fMRI data preprocessing
All rs-fMRI data were meticulously processed utilizing the Data Processing & Analysis of Brain Imaging (DPABI) toolbox (version 8.2) 27 alongside Statistical Parametric Mapping (SPM, version 12)28. Pre-processing steps included excluding the first 10 time points to achieve magnetization equilibrium, slice-timing correction, motion correction, spatial smoothing (full width at maximum (FWHM) = 6mm), linear detrending and bandpass filtering (0.01–0.1 Hz). Given their relevance to the current study, no nuisance regression analysis was conducted on gBOLD signals, CSF signals, or motion parameters.
Evaluation of gBOLD–CSF coupling strength
After all preconditioning procedures, we extracted gBOLD, Yeo-7 network BOLD and CSF fMRI signals (Fig2a) from the gray matter and CSF regions to explore the coupling changes. Cortical gray matter regions of interest (ROI) were defined according to the Harvard-Oxford cortical structural atlas, and gBOLD were extracted from the ROI. We obtained a dice score map between Automated Anatomical Labelling (AAL) atlas and the Yeo-7 Networks atlas from the neuroparc OSF registered repository29, 30 (Fig3). According to the dice score, the ninety brain regions of AAL were corresponded to the seven brain networks proposed by Thomas Yeo. Then each of Yeo-7 networks was given its own mask generated based on the above map using the REST Image Calculator Toolkit in the Resting-State fMRI Data Analysis Toolkit (REST) 31, 32. Consistent with validation by Fultz et al, CSF signals surrounding the inferior cerebellum were chosen due to their sensitivity to cerebrospinal fluid influx in this region. The ROIs pertaining to CSF were delineated manually on rs-fMRI and subsequently validated with T1WI.
Following the extraction of rs-fMRI signals from both cortical gray matter and CSF, we proceeded to compute the cross-correlation function between the BOLD signals and the CSF signals. This analysis was conducted to evaluate the coupling strength across a range of time lags, specifically from -10 to 10 seconds, utilizing Pearson correlation methods for the assessment. Given the symmetrical amplitude of negative peaks at +6 seconds and positive peaks at -6 seconds, the negative peak's BOLD-CSF correlation at +6 seconds was used to quantify BOLD-CSF coupling strength (Fig2b).
Statistical analysis
Pearson test was used to check the normality of the data. Non-parametric tests were used to compare differences in non-normal distribution data (age, duration of illness, MMSE score and UPDRS scores). Chi-square test was used to examine differences in classification indicators (gender and H-Y staging). The differences of gBOLD-CSF among PD-MCI group, PD-NC group and HC group were compared by ANOVA. Spearman correlation analysis was used to analyze the relationship between gBOLD-CSF coupling and clinical indicators (age and MMSE score), and Bonferroni multiple comparison correction test was used. Statistical analysis and graph generation using GraphPad Prism 8.0. Non-normally distributed data are described by interquartile spacing. The difference was statistically significant with p value <0.05.