Ability of a Set of Trunk Inertial Indexes of Gait to Identify Gait Instability and Recurrent Fallers in Parkinson's Disease

Differences in Trunk Acceleration-Derived Gait Indexes in Stroke Subjects with and without Stroke-Induced Immunosuppression

Background: Stroke-induced immunosuppression (SII) represents a negative rehabilitative prognostic factor associated with poor motor performance at discharge from a neurorehabilitation unit (NRB). This study aims to evaluate the association between SII and gait impairment at NRB admission. Methods: Forty-six stroke patients (65.4 ± 15.8 years, 28 males) and 42 healthy subjects (HS), matched for age, sex, and gait speed, underwent gait analysis using an inertial measurement unit at the lumbar level. Stroke patients were divided into two groups: (i) the SII group was defined using a neutrophil-to-lymphocyte ratio ≥ 5, and (ii) the immunocompetent (IC) group. Harmonic ratio (HR) and short-term largest Lyapunov's exponent (sLLE) were calculated as measures of gait symmetry and stability, respectively. Results: Out of 46 patients, 14 (30.4%) had SII. HR was higher in HS when compared to SII and IC groups (p < 0.01). HR values were lower in SII when compared to IC subjects (p < 0.01). sLLE was lower in HS when compared to SII and IC groups in the vertical and medio-lateral planes (p ≤ 0.01 for all comparisons). sLLE in the medio-lateral plane was higher in SII when compared to IC subjects (p = 0.04). Conclusions: SII individuals are characterized by a pronounced asymmetric gait and a more impaired dynamic gait stability. Our findings underline the importance of devising tailored rehabilitation programs in patients with SII. Further studies are needed to assess the long-term outcomes and the role of other clinical features on gait pattern.

Local Dynamic Stability of Trunk During Gait is Responsive to Rehabilitation in Subjects with Primary Degenerative Cerebellar Ataxia

This study aimed to assess the responsiveness to the rehabilitation of three trunk acceleration-derived gait indexes, namely the harmonic ratio (HR), the short-term longest Lyapunov's exponent (sLLE), and the step-to-step coefficient of variation (CV), in a sample of subjects with primary degenerative cerebellar ataxia (swCA), and investigate the correlations between their improvements (∆), clinical characteristics, and spatio-temporal and kinematic gait features. The trunk acceleration patterns in the antero-posterior (AP), medio-lateral (ML), and vertical (V) directions during gait of 21 swCA were recorded using a magneto-inertial measurement unit placed at the lower back before (T0) and after (T1) a period of inpatient rehabilitation. For comparison, a sample of 21 age- and gait speed-matched healthy subjects (HSmatched) was also included. At T1, sLLE in the AP (sLLEAP) and ML (sLLEML) directions significantly improved with moderate to large effect sizes, as well as SARA scores, stride length, and pelvic rotation. sLLEML and pelvic rotation also approached the HSmatched values at T1, suggesting a normalization of the parameter. HRs and CV did not significantly modify after rehabilitation. ∆sLLEML correlated with ∆ of the gait subscore of the SARA scale (SARAGAIT) and ∆stride length and ∆sLLEAP correlated with ∆pelvic rotation and ∆SARAGAIT. The minimal clinically important differences for sLLEML and sLLEAP were ≥ 36.16% and ≥ 28.19%, respectively, as the minimal score reflects a clinical improvement in SARA scores. When using inertial measurement units, sLLEAP and sLLEML can be considered responsive outcome measures for assessing the effectiveness of rehabilitation on trunk stability during walking in swCA.

Optimizing Rare Disease Gait Classification through Data Balancing and Generative AI: Insights from Hereditary Cerebellar Ataxia

The interpretability of gait analysis studies in people with rare diseases, such as those with primary hereditary cerebellar ataxia (pwCA), is frequently limited by the small sample sizes and unbalanced datasets. The purpose of this study was to assess the effectiveness of data balancing and generative artificial intelligence (AI) algorithms in generating synthetic data reflecting the actual gait abnormalities of pwCA. Gait data of 30 pwCA (age: 51.6 ± 12.2 years; 13 females, 17 males) and 100 healthy subjects (age: 57.1 ± 10.4; 60 females, 40 males) were collected at the lumbar level with an inertial measurement unit. Subsampling, oversampling, synthetic minority oversampling, generative adversarial networks, and conditional tabular generative adversarial networks (ctGAN) were applied to generate datasets to be input to a random forest classifier. Consistency and explainability metrics were also calculated to assess the coherence of the generated dataset with known gait abnormalities of pwCA. ctGAN significantly improved the classification performance compared with the original dataset and traditional data augmentation methods. ctGAN are effective methods for balancing tabular datasets from populations with rare diseases, owing to their ability to improve diagnostic models with consistent explainability.

Dynamic Stability, Symmetry, and Smoothness of Gait in People with Neurological Health Conditions

Neurological disorders such as stroke, Parkinson's disease (PD), and severe traumatic brain injury (sTBI) are leading global causes of disability and mortality. This study aimed to assess the ability to walk of patients with sTBI, stroke, and PD, identifying the differences in dynamic postural stability, symmetry, and smoothness during various dynamic motor tasks. Sixty people with neurological disorders and 20 healthy participants were recruited. Inertial measurement unit (IMU) sensors were employed to measure spatiotemporal parameters and gait quality indices during different motor tasks. The Mini-BESTest, Berg Balance Scale, and Dynamic Gait Index Scoring were also used to evaluate balance and gait. People with stroke exhibited the most compromised biomechanical patterns, with lower walking speed, increased stride duration, and decreased stride frequency. They also showed higher upper body instability and greater variability in gait stability indices, as well as less gait symmetry and smoothness. PD and sTBI patients displayed significantly different temporal parameters and differences in stability parameters only at the pelvis level and in the smoothness index during both linear and curved paths. This study provides a biomechanical characterization of dynamic stability, symmetry, and smoothness in people with stroke, sTBI, and PD using an IMU-based ecological assessment.

Neurophysiology of cerebellar ataxias and gait disorders

There are numerous forms of cerebellar disorders from sporadic to genetic diseases. The aim of this chapter is to provide an overview of the advances and emerging techniques during these last 2 decades in the neurophysiological tests useful in cerebellar patients for clinical and research purposes. Clinically, patients exhibit various combinations of a vestibulocerebellar syndrome, a cerebellar cognitive affective syndrome and a cerebellar motor syndrome which will be discussed throughout this chapter. Cerebellar patients show abnormal Bereitschaftpotentials (BPs) and mismatch negativity. Cerebellar EEG is now being applied in cerebellar disorders to unravel impaired electrophysiological patterns associated within disorders of the cerebellar cortex. Eyeblink conditioning is significantly impaired in cerebellar disorders: the ability to acquire conditioned eyeblink responses is reduced in hereditary ataxias, in cerebellar stroke and after tumor surgery of the cerebellum. Furthermore, impaired eyeblink conditioning is an early marker of cerebellar degenerative disease. General rules of motor control suggest that optimal strategies are needed to execute voluntary movements in the complex environment of daily life. A high degree of adaptability is required for learning procedures underlying motor control as sensorimotor adaptation is essential to perform accurate goal-directed movements. Cerebellar patients show impairments during online visuomotor adaptation tasks. Cerebellum-motor cortex inhibition (CBI) is a neurophysiological biomarker showing an inverse association between cerebellothalamocortical tract integrity and ataxia severity. Ataxic gait is characterized by increased step width, reduced ankle joint range of motion, increased gait variability, lack of intra-limb inter-joint and inter-segmental coordination, impaired foot ground placement and loss of trunk control. Taken together, these techniques provide a neurophysiological framework for a better appraisal of cerebellar disorders.

Predicting the fatigue in Parkinson's disease using inertial sensor gait data and clinical characteristics

Objectives

The study aimed to analyze the clinical features and gait characteristics of patients with Parkinson's disease (PD) who also suffer from fatigue and to develop a model that can help identify fatigue states in the early stages of PD.

Methodology

A total of 81 PD patients have been enrolled for the Parkinson's Fatigue Scale (PFS-16) assessment and divided into two groups: patients with or without fatigue. Neuropsychological assessments of the two groups, including motor and non-motor symptoms, were collected. The patient's gait characteristics were collected using a wearable inertial sensor device.

Results

PD patients who experienced fatigue had a more significant impairment of motor symptoms than those who did not, and the experience of fatigue became more pronounced as the disease progressed. Patients with fatigue had more significant mood disorders and sleep disturbances, which can lead to a poorer quality of life. PD patients with fatigue had shorter step lengths, lower velocity, and stride length and increased stride length variability. As for kinematic parameters, PD patients with fatigue had lower shank-forward swing max, trunk-max sagittal angular velocity, and lumbar-max coronal angular velocity than PD patients without fatigue. The binary logistic analysis found that Movement Disorder Society-Unified Parkinson's Disease Rating Scale-I (MDS-UPDRS-I) scores, Hamilton Depression Scale (HAMD) scores, and stride length variability independently predicted fatigue in PD patients. The area under the curve (AUC) of these selected factors in the receiver operating characteristic (ROC) analysis was 0.900. Moreover, HAMD might completely mediate the association between Hamilton Anxiety Scale (HAMA) scores and fatigue (indirect effect: β = 0.032, 95% confidence interval: 0.001-0.062), with a percentage of mediation of 55.46%.

Conclusion

Combining clinical characteristics and gait cycle parameters, including MDS-UPDRS-I scores, HAMD scores, and stride length variability, can identify PD patients with a high fatigue risk.

Machine learning based estimation of dynamic balance and gait adaptability in persons with neurological diseases using inertial sensors

Poor dynamic balance and impaired gait adaptation to different contexts are hallmarks of people with neurological disorders (PwND), leading to difficulties in daily life and increased fall risk. Frequent assessment of dynamic balance and gait adaptability is therefore essential for monitoring the evolution of these impairments and/or the long-term effects of rehabilitation. The modified dynamic gait index (mDGI) is a validated clinical test specifically devoted to evaluating gait facets in clinical settings under a physiotherapist's supervision. The need of a clinical environment, consequently, limits the number of assessments. Wearable sensors are increasingly used to measure balance and locomotion in real-world contexts and may permit an increase in monitoring frequency. This study aims to provide a preliminary test of this opportunity by using nested cross-validated machine learning regressors to predict the mDGI scores of 95 PwND via inertial signals collected from short steady-state walking bouts derived from the 6-minute walk test. Four different models were compared, one for each pathology (multiple sclerosis, Parkinson's disease, and stroke) and one for the pooled multipathological cohort. Model explanations were computed on the best-performing solution; the model trained on the multipathological cohort yielded a median (interquartile range) absolute test error of 3.58 (5.38) points. In total, 76% of the predictions were within the mDGI's minimal detectable change of 5 points. These results confirm that steady-state walking measurements provide information about dynamic balance and gait adaptability and can help clinicians identify important features to improve upon during rehabilitation. Future developments will include training of the method using short steady-state walking bouts in real-world settings, analysing the feasibility of this solution to intensify performance monitoring, providing prompt detection of worsening/improvements, and complementing clinical assessments.

Multiscale Entropy Algorithms to Analyze Complexity and Variability of Trunk Accelerations Time Series in Subjects with Parkinson's Disease

The aim of this study was to assess the ability of multiscale sample entropy (MSE), refined composite multiscale entropy (RCMSE), and complexity index (CI) to characterize gait complexity through trunk acceleration patterns in subjects with Parkinson's disease (swPD) and healthy subjects, regardless of age or gait speed. The trunk acceleration patterns of 51 swPD and 50 healthy subjects (HS) were acquired using a lumbar-mounted magneto-inertial measurement unit during their walking. MSE, RCMSE, and CI were calculated on 2000 data points, using scale factors (τ) 1-6. Differences between swPD and HS were calculated at each τ, and the area under the receiver operating characteristics, optimal cutoff points, post-test probabilities, and diagnostic odds ratios were calculated. MSE, RCMSE, and CIs showed to differentiate swPD from HS. MSE in the anteroposterior direction at τ4 and τ5, and MSE in the ML direction at τ4 showed to characterize the gait disorders of swPD with the best trade-off between positive and negative posttest probabilities and correlated with the motor disability, pelvic kinematics, and stance phase. Using a time series of 2000 data points, a scale factor of 4 or 5 in the MSE procedure can yield the best trade-off in terms of post-test probabilities when compared to other scale factors for detecting gait variability and complexity in swPD.

Identification of a Gait Pattern for Detecting Mild Cognitive Impairment in Parkinson's Disease

The aim of this study was to determine a gait pattern, i.e., a subset of spatial and temporal parameters, through a supervised machine learning (ML) approach, which could be used to reliably distinguish Parkinson's Disease (PD) patients with and without mild cognitive impairment (MCI). Thus, 80 PD patients underwent gait analysis and spatial-temporal parameters were acquired in three different conditions (normal gait, motor dual task and cognitive dual task). Statistical analysis was performed to investigate the data and, then, five ML algorithms and the wrapper method were implemented: Decision Tree (DT), Random Forest (RF), Naïve Bayes (NB), Support Vector Machine (SVM) and K-Nearest Neighbour (KNN). First, the algorithms for classifying PD patients with MCI were trained and validated on an internal dataset (sixty patients) and, then, the performance was tested by using an external dataset (twenty patients). Specificity, sensitivity, precision, accuracy and area under the receiver operating characteristic curve were calculated. SVM and RF showed the best performance and detected MCI with an accuracy of over 80.0%. The key features emerging from this study are stance phase, mean velocity, step length and cycle length; moreover, the major number of features selected by the wrapper belonged to the cognitive dual task, thus, supporting the close relationship between gait dysfunction and MCI in PD.

Identification of Gait Unbalance and Fallers Among Subjects with Cerebellar Ataxia by a Set of Trunk Acceleration-Derived Indices of Gait

This study aimed to assess the ability of 25 gait indices to characterize gait instability and recurrent fallers among persons with primary degenerative cerebellar ataxia (pwCA), regardless of gait speed, and investigate their correlation with clinical and kinematic variables. Trunk acceleration patterns were acquired during the gait of 34 pwCA, and 34 age- and speed-matched healthy subjects (HS_matched) using an inertial measurement unit. We calculated harmonic ratios (HR), percent recurrence, percent determinism, step length coefficient of variation, short-time largest Lyapunov exponent (sLLE), normalized jerk score, log-dimensionless jerk (LDLJ-A), root mean square (RMS), and root mean square ratio of accelerations (RMSR) in each spatial direction for each participant. Unpaired t-tests or Mann-Whitney tests were performed to identify significant differences between the pwCA and HS_matched groups. Receiver operating characteristics were plotted to assess the ability to characterize gait alterations in pwCA and fallers. Optimal cutoff points were identified, and post-test probabilities were calculated. The HRs showed to characterize gait instability and pwCA fallers with high probabilities. They were correlated with disease severity and stance, swing, and double support duration, regardless of gait speed. sLLEs, RMSs, RMSRs, and LDLJ-A were slightly able to characterize the gait of pwCA but failed to characterize fallers.

The kinectome: A comprehensive kinematic map of human motion in health and disease

Human voluntary movement stems from the coordinated activations in space and time of many musculoskeletal segments. However, the current methodological approaches to study human movement are still limited to the evaluation of the synergies among a few body elements. Network science can be a useful approach to describe movement as a whole and to extract features that are relevant to understanding both its complex physiology and the pathophysiology of movement disorders. Here, we propose to represent human movement as a network (that we named the kinectome), where nodes represent body points, and edges are defined as the correlations of the accelerations between each pair of them. We applied this framework to healthy individuals and patients with Parkinson's disease, observing that the patients' kinectomes display less symmetrical patterns as compared to healthy controls. Furthermore, we used the kinectomes to successfully identify both healthy and diseased subjects using short gait recordings. Finally, we highlighted topological features that predict the individual clinical impairment in patients. Our results define a novel approach to study human movement. While deceptively simple, this approach is well-grounded, and represents a powerful tool that may be applied to a wide spectrum of frameworks.

Harmonic ratio is the most responsive trunk-acceleration derived gait index to rehabilitation in people with Parkinson's disease at moderate disease stages

BACKGROUND

Harmonic ratios (HRs), recurrence quantification analysis in the antero-posterior direction (RQA_detAP), and stride length coefficient of variation (CV) have recently been shown to characterize gait abnormalities and fall risk in people with Parkinson's disease (pwPD) at moderate disease stages.

RESEARCH QUESTION

This study aimed to i) assess the internal and external responsiveness to rehabilitation of HR, RQA_detAP, and CV, ii) identify the baseline predictors of normalization of the gait stability indexes, and iii) investigate the correlations between the gait indexes modifications (∆) and clinical and kinematic ∆s in pwPD at Hoehn and Yahr disease staging classification 3.

METHODS

The trunk acceleration patterns of 21 pwPD and 21 age- and speed-matched healthy subjects (HSmatched) were acquired during gait using an inertial measurement unit at baseline (T0). pwPD were also assessed after a 4-week rehabilitation period (T1). Each participant's HR in the antero-posterior (HR_AP), medio-lateral (HR_ML), and vertical directions, RQA_detAP, CV, spatio-temporal, and kinematic variables were calculated.

RESULTS

At T1, HR_AP and HR_ML improved to normative values and showed high internal and external responsiveness. Lower HRs and higher pelvic rotation values at baseline were predictors of ∆HRs. A minimal clinically important difference (MCID) ≥ 21.5 % is required to normalize HR_AP with 95 % probability. MCID ≥ 36.9 % is required to normalize HR_ML with 92 % probability. ∆HR_AP correlated with ∆HR_ML and both correlated with ∆stride length and ∆pelvic rotation, regardless of ∆gait speed. RQAdetAP and step length CV were not responsive to rehabilitation.

SIGNIFICANCE

When using inertial measurement units, HR_AP and HR_ML can be considered as responsive outcome measures for assessing the effectiveness of rehabilitation on trunk smoothness during walking in pwPD at moderate disease stages.

Machine Learning Approach to Support the Detection of Parkinson's Disease in IMU-Based Gait Analysis

The aim of this study was to determine which supervised machine learning (ML) algorithm can most accurately classify people with Parkinson’s disease (pwPD) from speed-matched healthy subjects (HS) based on a selected minimum set of IMU-derived gait features. Twenty-two gait features were extrapolated from the trunk acceleration patterns of 81 pwPD and 80 HS, including spatiotemporal, pelvic kinematics, and acceleration-derived gait stability indexes. After a three-level feature selection procedure, seven gait features were considered for implementing five ML algorithms: support vector machine (SVM), artificial neural network, decision trees (DT), random forest (RF), and K-nearest neighbors. Accuracy, precision, recall, and F1 score were calculated. SVM, DT, and RF showed the best classification performances, with prediction accuracy higher than 80% on the test set. The conceptual model of approaching ML that we proposed could reduce the risk of overrepresenting multicollinear gait features in the model, reducing the risk of overfitting in the test performances while fostering the explainability of the results.

Sensitivity to gait improvement after levodopa intake in Parkinson's disease: A comparison study among synthetic kinematic indices

The synthetic indices are widely used to describe balance and stability during gait. Some of these are employed to describe the gait features in Parkinson's disease (PD). However, the results are sometimes inconsistent, and the same indices are rarely used to compare the individuals affected by PD before and after levodopa intake (OFF and ON condition, respectively). Our aim was to investigate which synthetic measure among Harmonic Ratio, Jerk Ratio, Golden Ratio and Trunk Displacement Index is representative of gait stability and harmony, and which of these are more sensitive to the variations between OFF and ON condition. We found that all indices, except the Jerk Ratio, significantly improve after levodopa. Only the improvement of the Trunk Displacement Index showed a direct correlation with the motor improvement measured through the clinical scale UPDRS-III (Unified Parkinson's Disease Rating Scale-part III). In conclusion, we suggest that the synthetic indices can be useful to detect motor changes induced by, but not all of them clearly correlate with the clinical changes achieved with the levodopa administration. In our analysis, only the Trunk Displacement Index was able to show a clear relationship with the PD clinical motor improvement.

Quantitative Characterization of Motor Control during Gait in Dravet Syndrome Using Wearable Sensors: A Preliminary Study

Dravet syndrome (DS) is a rare and severe form of genetic epilepsy characterized by cognitive and behavioural impairments and progressive gait deterioration. The characterization of gait parameters in DS needs efficient, non-invasive quantification. The aim of the present study is to apply nonlinear indexes calculated from inertial measurements to describe the dynamics of DS gait. Twenty participants (7 M, age 9-33 years) diagnosed with DS were enrolled. Three wearable inertial measurement units (OPAL, Apdm, Portland, OR, USA; Miniwave, Cometa s.r.l., Italy) were attached to the lower back and ankles and 3D acceleration and angular velocity were acquired while participants walked back and forth along a straight path. Segmental kinematics were acquired by means of stereophotogrammetry (SMART, BTS). Community functioning data were collected using the functional independence measure (FIM). Mean velocity and step width were calculated from stereophotogrammetric data; fundamental frequency, harmonic ratio, recurrence quantification analysis, and multiscale entropy (τ = 1...6) indexes along anteroposterior (AP), mediolateral (ML), and vertical (V) axes were calculated from trunk acceleration. Results were compared to a reference age-matched control group (112 subjects, 6-25 years old). All nonlinear indexes show a disruption of the cyclic pattern of the centre of mass in the sagittal plane, quantitatively supporting the clinical observation of ataxic gait. Indexes in the ML direction were less altered, suggesting the efficacy of the compensatory strategy (widening the base of support). Nonlinear indexes correlated significantly with functional scores (i.e., FIM and speed), confirming their effectiveness in capturing clinically meaningful biomarkers of gait.

Predicting Axial Impairment in Parkinson's Disease through a Single Inertial Sensor

Background: Current telemedicine approaches lack standardised procedures for the remote assessment of axial impairment in Parkinson’s disease (PD). Unobtrusive wearable sensors may be a feasible tool to provide clinicians with practical medical indices reflecting axial dysfunction in PD. This study aims to predict the postural instability/gait difficulty (PIGD) score in PD patients by monitoring gait through a single inertial measurement unit (IMU) and machine-learning algorithms. Methods: Thirty-one PD patients underwent a 7-m timed-up-and-go test while monitored through an IMU placed on the thigh, both under (ON) and not under (OFF) dopaminergic therapy. After pre-processing procedures and feature selection, a support vector regression model was implemented to predict PIGD scores and to investigate the impact of L-Dopa and freezing of gait (FOG) on regression models. Results: Specific time- and frequency-domain features correlated with PIGD scores. After optimizing the dimensionality reduction methods and the model parameters, regression algorithms demonstrated different performance in the PIGD prediction in patients OFF and ON therapy (r = 0.79 and 0.75 and RMSE = 0.19 and 0.20, respectively). Similarly, regression models showed different performances in the PIGD prediction, in patients with FOG, ON and OFF therapy (r = 0.71 and RMSE = 0.27; r = 0.83 and RMSE = 0.22, respectively) and in those without FOG, ON and OFF therapy (r = 0.85 and RMSE = 0.19; r = 0.79 and RMSE = 0.21, respectively). Conclusions: Optimized support vector regression models have high feasibility in predicting PIGD scores in PD. L-Dopa and FOG affect regression model performances. Overall, a single inertial sensor may help to remotely assess axial motor impairment in PD patients.

The Effects of Intensive Neurorehabilitation on Sequence Effect in Parkinson's Disease Patients With and Without Freezing of Gait

Background: The sequence effect (SE), defined as a reduction in amplitude of repetitive movements, is a common clinical feature of Parkinson's disease (PD) and is supposed to be a major contributor to freezing of gait (FOG). During walking, SE manifests as a step-by-step reduction in step length when approaching a turning point or gait destination, resulting in the so-called destination sequence effect (dSE). Previous studies explored the therapeutic effects of several strategies on SE, but none of them evaluated the role of an intensive rehabilitative program. Objectives: Here we aim to study the effects of a 4-week rehabilitative program on dSE in patients with PD with and without FOG. Methods: Forty-three patients (30 males, 70.6 ± 7.5 years old) with idiopathic PD were enrolled. The subjects were divided into two groups: patients with (PD + FOG, n = 23) and without FOG (PD - FOG, n = 20). All patients underwent a standardized 4-week intensive rehabilitation in-hospital program. At hospital admission (T0) and discharge (T1), all subjects were evaluated with an inertial gait analysis for dSE recording. Results: At T0, the dSE was more negative in the PD + FOG group (-0.80 ± 0.6) when compared to the PD - FOG group (-0.39 ± 0.3) (p = 0.007), even when controlling for several clinical and demographic features. At T1, the dSE was reduced in the overall study population (p = 0.001), with a more pronounced improvement in the PD + FOG group (T0: -0.80 ± 0.6; T1: -0.23 ± 0.4) when compared to the PD - FOG group (T0: -0.39 ± 0.3; T1: -0.22 ± 0.5) (p = 0.012). At T1, we described in the overall study population an improvement in speed, cadence, stride duration, and stride length (p = 0.001 for all variables). Conclusions: dSE is a core feature of PD gait dysfunction, specifically in patients with FOG. A 4-week intensive rehabilitative program improved dSE in PD patients, exerting a more notable beneficial effect in the PD + FOG group.