Objective and study design
A cross-sectional and unblinded study assessing the effect of levodopa and stimulation in FOG emerging after STN-DBS, in five different experimental conditions.
Primary Outcome: the difference in the number of FOG episodes between Med OFF/Stim OFF and Med ON/Stim ON 130 Hz conditions.
Study Population
We reviewed the electronical medical records off all PD patients operated at our center since the beginning of the DBS program (2006), to identify those with Med On/Stim ON FOG. The presence of FOG was defined by a score ≥ 2 on item 3.11 (freezing) of the MDS-UPDRS part III on the Med ON/Stim ON state. FOG was considered only after 6 months into the post-surgical period, as we considered 6 months the time required for treatment stabilization. We called into clinic for assessment those still alive and available for evaluation. Patients unable to walk without ambulatory aids in the Med ON/Stim ON state, and those with severe osteo-articular or other non-neurological issues affecting gait, were excluded. All patients had been selected for surgery according to the CAPSIT-PD protocol, and none had levodopa-resistant axial signs before surgery during best ON state. DBS surgery followed standard stereotactic techniques, and postoperative neuroimaging confirmed lead position. DBS programming parameters were optimized regularly by a movement disorders specialist.
Clinical Evaluation
Clinical evaluations were conducted under five different conditions following a pre-determined order: 1) Medication OFF/Stimulation ON (Med OFF/Stim ON),, 2) Medication OFF/Stimulation OFF (Med OFF/Stim OFF) 3) Medication ON/Stimulation OFF (Med ON/ Stim OFF), 4) Medication ON/Stimulation ON at 130 Hz High-frequency stimulation (Med ON/Stim ON HFS), and 5) Medication ON/Stimulation ON at 60 Hz Low-frequency stimulation (Med ON/Stim ON LFS).
The full assessment lasted approximately 4 hours in a single morning session. Patients were evaluated on the "practical OFF drug" condition after 12 hours of medication withdrawal. A levodopa challenge test (LCT) was performed using the same dose of LD as used on the pre-surgery LCT. The Stim ON condition used the stimulation settings that had shown the best clinical results over the previous 6 months for each patient, while we adjusted the total energy delivered when testing LFS according to the total electrical energy delivered (TEED) formula: TEED (1 second) = voltage2 X frequency x pulse width/impedance).15,73. No patient was on LFS at the time of evaluation. A 30-minute interval was maintained between frequency changes and after stimulation arrest.
Motor evaluations using the MDS-UPDRS part III74, gait/axial sub-score (MDS-UPDRS part III items 3.9–3.12)74, non-gait/axial sub-score (MDS-UPDRS III – gait/axial sub-score), Hoehn and Yahr scale75 and the Stand-Walk and Sit test (SWS-test)15 were performed on each of the five conditions. The SWS test is a standardized, timed test where subjects walk a 14-meter distance between sitting and standing. The overall duration of the test (time to walk 14 meters, SWS time) and the number of FOG episodes (#FOG episodes) were recorded. FOG was defined as a transient incapacity to move forward, despite the intention to walk, including both akinetic and "trembling in place" forms.1,76 Patients performed the SWS test three times, and the data was averaged.
We defined responsiveness to LD and stimulation as any decrease on the number of FOG episodes in the SWS (i.e, a clinical improvement). Accordingly, patients who worsened or had no change with LD or stimulation were classified as LD-resistant or stimulation-resistant, respectively.
Imaging procedures - Electrode reconstruction and localization
To ensure electrode placement across subjects, we reconstructed the DBS electrodes of all patients whose neuroimaging data was possible to retrospectively retrieve (n = 11). This was performed through the advanced processing pipeline in Lead-DBS77. For each patient, a post-operative CT scan was linearly coregistered to a pre-operative structural MRI (T1w), and then transformed into the ICBM 2009b NLIN asymmetric MNI space78 both using Advanced Normalization Tools79. A brain shift correction was employed and electrode trajectories were reconstructed automatically using the PaCER algorithm80, following manual readjustments. Group visualization was performed using Lead Group with the DISTAL Atlas segmentation81.
Development of an inertial sensor-based classifier for FOG
We used data collected from previous studies 62 to develop a model for automatic detection of FOG. In summary, 20 PD patients (freezers and non-freezers) wearing seven wearable sensors (inertial measurement units, IMUs) fixed to different body parts performed a self-paced gait task. 180° degrees turns were removed from the analysis leaving only straight-line walking. Clinical annotation of FOG episodes (presence and duration) was made by a PD expert clinician (RB) based on video recordings.
Each IMU consisted of a tri-axial accelerometer, a gyroscope and a magnetometer (Xsens Technologies, Enschede, The Netherlands), that was fixed to the patient’s body using a Velcro elastic band. The inertial sensors were positioned in pelvis, right and left thighs, legs and feet. We used the raw signal collected by the IMUs and applied a supervised learning methodology on labelled data. This strategy trained a model to distinguish between different labels: FOG and non-FOG.
Continuous IMU data was processed into a one-second dataset, in which all 9 dimensions were split into non-overlapping sequential 1-second slices, each corresponding to 100 time-points, given the collection rate of 100Hz. Additionally, each resulting slice was then associated with the clinician label of belonging to a FOG episode or not. Considering the continuous nature of sensor data, an algorithm for identifying FOG events was envisaged as a time series (TS) classification problem. Deep learning (DL) has been increasingly used in the TS area, especially for multivariate time series problems.82 Convolutional Neural Network (CNN) based architecture was used to build a FOG detection model, consisting of three sequential convolution modules 83The model was built in 13 patients and 7 were used to test the model on unseen data. FOG % was the main output of the model, representing the percentage of FOG presenting during straight gait.
A detailed description of the development and validation of this model is provided in Supplementary Methods.
Wearable sensor-based gait analysis:
During the SWS test, patients wore the same seven sensors as previously described in the same locations (pelvis, right and left thighs, legs and feet), during each of the five stimulation/medication conditions. The data collected from the IMUs were processed using the KINETIKOS (Coimbra, Portugal) cloud-based platform to reconstruct each subject's body motion using a 3D kinematic biomechanical model. Each trial out of the 3 SWS trials were individually computed and results were averaged. A final dataset consisted of 30 variables organized into 4 domains (spatio-temporal, asymmetry, variability, and non-linear metrics) selected based on their relevance in the literature (Supplementary Table S12). Spatiotemporal features were measured on both sides of the body, and the “worst-side” score was selected based on clinical assessment. For spatiotemporal variables, coefficients of variation (e.g. standard deviation of variable X score / mean of variable X score) and asymmetries [(variable X score on the right side – variable X score on the left side) / (variable X score on the right side + variable X score on the left side)] were calculated.
Statistical Analysis
The primary outcome of the study was the difference in the number of FOG episodes between Med OFF/Stim OFF and Med ON/Stim ON 130 Hz conditions. Secondary outcomes were: a) the change in the number of FOG episodes between the Med OFF/Stim OFF and the Med ON/Stim OFF and the Med OFF/Stim ON condition; b) the differences in the number of FOG episodes between the 130Hz and the 60 HZ conditions c) the association between the automatic detection of FOG and FOG assessed by clinical gait metrics. Exploratory analysis was performed to dissect kinematic dimensions related to FOG.
The statistical analysis was conducted in alignment with the prespecified hypothesis defined in our primary and secondary outcomes. This was done to mitigate the risk of Type II errors. Summary statistics were presented as mean (± SD) and median (IQR). To compare groups, paired non-parametric statistics were used (Wilcoxon signed rank test for comparisons of 2 groups and Friedman test for more than 2 groups with corrections to multiple comparisons performed using Dunn’s test). Spearman correlation was used to assess the association between kinematic and clinical variables. A Principal Component Analysis (PCA) used a dimensionality reduction method on the individual gait metrics to facilitate interpretation of data. Significance level was set at 0.05. Data processing was conducted using Matlab R2021a, Python 3.8 and R 4.2.2. Statistical analysis used Graphpad Prism 10.1.1.