Subjects
In this study, 41 patients with advanced Parkinson's disease (Hoehn-Yahr stage Ⅲ ~ IV, age 58.6 ±7.99 years, duration of disease 8.95 ±3.78 years) were enrolled. The selection criteria included more than 5 years of diagnosis of PD, eligibility for DBS surgery, and no taboos in MRI. All patients were treated with STN-DBS surgery. The details of the subjects' demographic and clinical characteristics are shown in Table 1. All patients were given informed written consent before participating the experiment and the study was approved by the local ethics board.
Table 1
Relevant values are reported as mean ± standard deviation. OFF-medication was defined as the practically defined OFF state with overnight medication withdrawal. Post-operative UPDRS-III values are ON-stimulation with bilateral clinical stimulation settings over first year after DBS.
Patient demographics and clinical characteristics
|
|
|
pre-surgery
|
1 year follow up
|
P value
|
Age(years)
|
59.02±7.80
|
-
|
-
|
Disease duration(years)
|
9.10±3.72
|
-
|
-
|
UPDRS III baseline
|
57.04±12.28
|
52.31±15.98
|
0.137
|
UPDRS III "on" dopamine
|
29.73±11.12
|
36.12±19.97
|
0.056
|
UPDRS III "on" DBS
|
-
|
29.02±12.91
|
-
|
H-Y stage
|
3.57±0.78
|
3.20±0.53
|
<0.05
|
LEDD(mg)
|
642.68±384.40
|
467.57±240.85
|
<0.05
|
PDQ39
|
69.08±16.92
|
48.29±24.47
|
<0.0001
|
fMRI experiment
Functional and structural MRI data were acquired with a 3T GE (Achieva TX) MRI scanner on OFF medication state before DBS surgery, following 8-h period of medication withdrawal. Structural images were acquired using a sagittal magnetization prepared rapid gradient echo three-dimensional T1-weighted sequence (repetition time [TR] = 8.5 ms, echo time [TE] = 3.2 ms, inversion time [TI] = 450 ms, flip angle [FA] = 12◦). Functional MRI images were obtained using the following SE-EPI sequence: repetition time = 2,000 ms, echo time = 30 ms, slice thickness/gap = 3.6/0 mm, axial slices = 38 layers, flip angle = 90◦, FOV = 256 × 256 mm, matrix size = 64 × 64, and scanning time = 242s.
DBS surgery
All patients underwent standard frame-based DBS implantation. Two quadrupole DBS electrodes (L301C electrode of Beijing Pins Medical Co., Ltd., China) and a LFP sensing nerve stimulator (G102R, Beijing Pins Medical Co., Ltd.) were received. Accurate electrode placement was confirmed by intraoperative microelectrode recording, intraoperative motor response to stimulation, and postoperative CT. The stimulation was started one month after the operation, and regular return visits were conducted to adjust the stimulation parameters.
Clinical evaluation
Each patient underwent two clinical evaluations: UPDRS-III was measured during the pre-surgical inpatient on and during postoperative DBS programming clinical visits. The preoperative score was collected within 1 week prior to DBS, and the postoperative ON stimulation score was collected 12 months after DBS placement. All UPDRS-III scores were evaluated on OFF medication state, after at least 8 hours of medication withdrawal. At the end of each visit, the ON-stimulation rating was performed after the completion of DBS settings. We calculated the preoperative off-medication and postoperative off-medicine and on-stimulation score changes to determine the effect of DBS. The average score during the entire postoperative period was used to mitigate random fluctuations and to represent the actual performance of each subject after receiving DBS.
Date preprocessing
MRI data were analyzed using Data Processing & Analysis for Brain Imaging open-source package(DPABI V6.0_210501, based on SPM8 and MATLAB)[25]. The standard pre-processing included:(a) Remove the first five volumes to stabilize magnetization;(b) Slice timing correction and realignment to correct the time differences existed in each slice on the same volume;(c) Reorient functional and structural images by manually alignement;(d) Co-register of structural images and functional images, and segment the structural images into gray matter, white matter, and cerebrospinal fluid; (e) Regress nuisance covariates including Friston 24 head motion parameters, white matter, and cerebrospinal fluid signals;(f) Normalize the functional images to Montreal Neurological Institute (MNI) space based on the diffeomorphic anatomical registration by exponentiated lie algebra, and resample the images to 3.0x3.0x3.0 mm; (g) Spatial smoothing using a Gaussian kernel with 6mm FWHM; (h) Temporal filtering to 0.01-0.08Hz.
No participant had a head motion of more than 1.5 mm maximum displacement in any of the x, y, or z directions or 1.5 of any angular motion throughout the course of the scan, and data greater than this threshold were excluded from further analysis.
Brain connectivity estimation
PD25 template with 16 Region of Interests (ROIs) which was specifically defined based on a Parkinson’s disease cohort was used to construct the FC matrix. The average time course of each ROI was calculated, and the functional connection was generated by computing the Pearson’s correlation coefficients between each pair of ROIs. Fisher z-transformation was used to convert r value to z value to increase the normality of correlation distribution. Connectivity between the putamen and other 15 ROIs was chosen for further analysis.
Statistical Analysis
Connectivity between the putamen and other ROIs was correlated with the overall changes in UPDRS-III in order to explore the relationship between preoperative resting-state functional connectivity and postoperative motor improvement. The level of significance was set at p<0.05. The False Discovery Rate (FDR) correction was used for multiple comparison correction. FDR correction was performed using Matlab. Demographic information, such as age at surgery, gender, and disease duration, were not assessed using a t-test, since our data were conducted for the same subject pre-surgery and post-surgery.
To investigate the predictive power of the corrections between FC and overall change in UPDRS-III, Leave-One-Out Cross-Validator (LOOCV) was adopted. Linear regression models were performed with Scikit-learning software using FC and other covariates, which included medication improvement rate and preoperative UPDRS-III. Among them, the improvement value means the difference between preoperative OFF-medication and postoperative OFF-medication and ON-stimulation scores. We also carried out a series of ablation studies on various factors. For example, generate linear models with and without medication improvement rate or baseline UPDRS-III. Mean absolute error was used to evaluate the validity of the linear regression model.