1. Introduction
Freezing of gait (FoG) is a debilitating movement disorder [
1,
2] present in roughly 50% of individuals with advanced Parkinson’s disease (PD) [
3]. It is characterized as a “brief, episodic absence or marked reduction of forward progression of the feet despite the intention to walk” [
2], leading to falls, increased dependence [
4], and reduced quality of life [
5]. FoG is associated with levodopa off-state [
6] and longer disease duration [
7]. Several factors can trigger FoG, including step initiation, changing direction and environment contexts, narrow doorways, anxiety, walking under temporal constraints, or stress [
2]. On the other hand, some situations relieve FoG symptoms: sensory cues and attentional focus on gait [
8]. Factors that trigger and relieve FoG depend highly on cortical mechanisms (emotion, attention, executive, and visuospatial functions), which can be explained by compensating cortical circuits on the affected subcortical regions [
9].
Recent studies found functionally deficient prefrontal areas in PD individuals with FoG, such as the anterior cingulate (error monitoring, inhibition, and volition); the orbitofrontal cortex (decision making, sensory integration, reward expectation); and the dorsolateral prefrontal cortex (motor planning, organization and changing of context) compared to individuals without PD (nFoG) [
10,
11]. Current evidence from radiological studies addresses FoG to structural disorders in gray and white matter. Voxel-based morphometry studies have shown atrophy in frontal [
10,
12] and parietal [
10,
13] cortices, with some studies also reporting gray matter loss in temporal [
12] and occipital [
12,
14] cortices. Tractography, utilizing diffusion tensor imaging (DTI), allows for the visualization and analysis of white matter tracts, providing insights into the integrity of neural pathways critical for motor control. Alterations in white matter tracts, such as the corpus callosum (CC), right inferior longitudinal fasciculus (RILF), and superior longitudinal fasciculus (SLF), have also been consistently observed in individuals with FoG [
15,
16,
17].
Given FoG’s complexity, it is imperative to explore the role of white matter integrity in step initiation and overall gait performance. The periventricular white matter, comprising longitudinal, thalamic, and callosal fibers, has been implicated in various studies as crucial for motor control and coordination. However, the specific contributions of these tracts to FoG, particularly during the anticipatory postural adjustments (APA) required for step initiation, still need to be explored. Step initiation requires APAs to shift the center of mass to the support leg, releasing the moving leg to take a step. Studies have shown that individuals with FoG have disorders of APA (lower amplitude and longer latency) compared to PD individuals without FoG [
18], which are related to disorders in cortical functioning [
19]. Recent studies demonstrated potential gray matter regions linked to impaired APA in step initiation in FoG [
19,
20,
21,
22]. The involvement of white matter integrity in step initiation behavior in individuals with FoG remains unclear. In the present study, we aim to investigate the integrity of the white matter and its relationship with the step initiation behavior in individuals with FoG.
2. Materials and Methods
2.1. Participants
Forty-three individuals with moderate idiopathic PD (Hoehn & Yahr = 3; 26 with FoG; age: 62.97 ± 10.59 years; disease duration: 10.03 ± 5.94 years; equivalent dose of levodopa: 611.54 ± 429.23 mg·day
−1; motor score Unified Parkinson’s Disease Rating Scale - UPDRS: 31.09 ± 10.39; Montreal Cognitive Assessment - MoCA: 18.25 ± 2.51) participated in this study. The individuals were recruited in the Movement Disorders Division of the Medical School of the University of São Paulo. All individuals were diagnosed by neurologists and were submitted to the Unified Parkinson’s Disease Rating Scale - UPDRS on the day of the exam. Inclusion criteria were: (1) Hoehn & Yahr score of 3; (2) availability to engage in task training, brain imaging, and neuropsychological assessments; (3) capacity to walk independently, without assistance devices; (4) ability to understand instructions to perform the experimental motor tasks; (5) absence of neurological diseases other than PD; (6) absence of musculoskeletal impairments possibly affecting the performance of the experimental tasks. The exclusion criteria were as follows: (1) difficulty performing the experimental motor tasks; (2) drop out of any assessment; (3) poor quality of the brain volumes acquired during the fMRI: head motion above 1 mm [
23]. Individuals were classified as having FoG if they answered affirmatively the first question of the New Freezing of Gait Questionnaire (nFoGq) following a video showing examples of individuals experiencing FoG. All the FoG individuals scored 3 or 4 in the NFoGq question: “How frequently do you experience freezing episodes when initiating the first step?” indicating moderate or severe FoG during step initiation. All evaluations were taken during the ON state of dopaminergic medications. All participants provided informed consent, and the University of São Paulo Institutional Review Board approved the experimental procedures.
2.2. Behavioral Task: Step Initiation
The individuals performed 30 trials of a step initiation task with the right leg. Once on the platform, the position of their feet was marked with tape on the floor to keep the same initial position in each trial. Firstly, the individual was to gaze at a cross on a TV monitor placed 3 m ahead and relax. The following command displayed was to prepare (a circle that could be open or closed, depending on the randomization of the trials), then an open or closed circle (if the first was closed, the second was open, and vice-versa) to take a step (
Figure 1).
The participants performed the task on a force platform (AMTI OR6-6). For kinematic analysis, we placed a reflective marker on the lateral malleolus of both subjects’ legs. Twelve optoelectronic cameras (Vicon, model T10) tracked the markers. Kinematic and kinetic data were acquired at a sampling frequency of 200 Hz using a motion capture system (Vicon Nexus). A trigger synchronizes the motion capture system and the fNIRS system.
The anticipatory postural adjustment (APA) began with the initial lateral displacement of the center of pressure (CoP) toward the leading limb. The end of the APA, and therefore the start of the step, occurred when the free vertical moment changed direction (i.e., became negative when initiating with the right limb) [
24]. The following variables were analyzed: (a) APA duration, the time between the beginning of the APA and the beginning of the step; (b) APA amplitude normalized, the maximum amplitude of the last peak of mediolateral displacement of the CoP before the start of the step normalized by the distance between the malleoli; (c) Step length normalized by each individual’s foot size.
2.3. Magnetic Resonance Imaging
To acquire diffusion images, the individuals underwent a Magnetic Resonance Imaging exam in a 3.0 T scanner (Intera Achieva, Philips Healthcare, Best, The Netherlands) in which every individual used a thirty-two-channel head coil (Philips Healthcare, Best, The Netherlands). T1-weighed axial 3D images were acquired with isotropic voxel size 1 mm, FOV 256 × 256, 180 sagittal slices, 60 mm thick each, covering the whole brain, phase encoding anterior-posterior, TR = 7 ms, TE = 3.2ms, TI = 1120 ms, flip angle 8°, acquisition time 4 min 12 s. FLAIR (Axial Fluid-Attenuated Inversion Recovery) images were acquired with isotropic voxel size 1 mm, FOV 256 × 256, 176 sagittal slices, phase encoding anterior-posterior, TR = 5,000 ms, TE = 397 ms, TI = 1,800 ms, flip angle 120°, acquisition time 4 min 42 s. Diffusion-weighted imaging (DWI) was obtained using a single-shot spin echo EPI sequence. Whole brain diffusion images received b = 700 s/mm2 with 32 directions. The repetition time (TR) and echo time (TE) were 13 and 101 ms for this image. Extra brain volumes received four no diffusion weighting (b = 0 s/mm2) with opposing phase-encode directions (Anterior-Posterior and opposite). The structural images consisted of single-shot T1-weighed axial 3D images containing 180 slices, 60 mm thick each, and covering the whole brain. They followed the parameters: flip angle = 8°, TR/TE = 7/3.2 s, FOV = 740 mm², NSA = 1, Fast Imaging Mode = TSE, needing a total of 58.7 s to be taken. DWI images consisted of 72 transversal slices covering the whole brain, with 2 mm³ voxels, 128 x 128 mm matrix, 256 x 256 mm field of view and no gap between slices, following the parameters: TR/TE = 8705/30 ms, NSA = 1, 844 volumes, 90-degree flip angle, b-factor = 2, max b-factor = 1000, diffusion gradient = 31.4/16.3 and EPI factor = 67. The total length of DWI acquisition was 6 minutes 11.5 seconds.
2.4. Image Data Processing
Subcortical segmentation and fiber tracking of the SLF were automated and obtained using FreeSurfer (
https://surfer.nmr.mgh.harvard.edu, software package version). Next, TRACULA, a tool available on FreeSurfer, was used to reconstruct the SLF tract based on each of the images obtained. The tool calculated the tract’s likelihood to be next to (or in) each of the cortical and subcortical segmentation labels from FreeSurfer, using manually labeled tracts on training subjects and an MNI 152 template as a reference.
The automated cortical and subcortical segmentation and parcellation were processed using the multimodal recon-all processing stream option, combining both the high-resolution T1 3D-weighted and T2-FLAIR images intended to improve dural surface segmentation [
25], which includes motion correction, skull-stripping, registration, segmentation, smoothing, and parcellation mapping in FreeSurfer 6.0 (
https://surfer.nmr.mgh.harvard.edu/fswiki). Segmentation results for each subject were manually inspected and edited for accuracy by the first author (MPN), per Freesurfer guidelines. For global probabilistic tractography analyses, diffusion tensor images (DTI) were processed using Freesurfer’s TRActs Constrained by UnderLying Anatomy (TRACULA) tool [
26], using a standard pipeline, including eddy-current correction via affine registration, rotation of the gradients to match eddy-current correction, extraction of the brain mask, registration from individual DTI to the individual T1 image (intra-subject registration), registration to MNI space (inter-subject registration), cortical and white matter mask creation, tensor fitting for the extraction of tensor-based measures, and computing anatomical priors for tracts form the TRACULA atlas.
This method reconstructs probabilistic distributions of 18 major white matter tracts for every individual. It computes the following DTI metrics: fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD), and axial diffusivity (AD).
The tract of interest consisted of commissural tracts (Forceps Major and Minor), anterior thalamic radiation, cingulum-angular (infracallosal) bundle, cingulum-cingulate (supracallosal) bundle, corticospinal tract, inferior longitudinal fasciculus, superior longitudinal fasciculus [temporal part, parietal part], and the uncinate fasciculus (UF) in each hemisphere. We checked all reconstructed pathways visually; no failed tracts had aberrant trajectories. At the end of the processing, the software’s output is the FA, AD, RD, and MD values of the most probable voxels in the bundle.
2.5. Statistical Analyses
Shapiro-Wilk and Levene’s tests were used to assess data normality and homogeneity of variance, respectively. For the two groups (FoG and nFoG) separately, we performed linear multiple regressions using the stepwise method to explain the variance of dependent variables (measurements of step initiation: APA duration, APA amplitude normalized, Step length normalized). First, the univariate analyses were used to test which tracts (Forceps Major and Minor, ATR, CAB, CCG, CST, SLFt, SLFp, UF in each hemisphere) would be associated with the dependent variables (APA duration, APA amplitude normalized, Step length normalized). Afterward, to explain the variance of the dependent variables, we included the factors in the linear multivariate analysis using the stepwise model if they presented a
P value ≤0.10 and a correlation of lower than 0.6 between them to avoid collinearity [
27]. The association between dependent and independent variables was assessed through the Spearman test. Statistical procedures were performed using the software SAS 9.4 (Institute Inc., Cary, NC, USA), and the significance level was set at
P≤0.05.
4. Discussion
Authors should discuss the results and how they can be interpreted from the perspective of previous studies and of the working hypotheses. The findings and their implications should be discussed in the broadest context possible. Future research directions may also be highlighted.
This study investigates whether white matter integrity can predict step initiation behavior in individuals with PD and FoG. To address this question, we included the tracts involved in the pathophysiology of PD and FoG in our regression model. This is the first investigation into the relationship between tract integrity and step initiation impairment in individuals with FoG. Our results show that the radial diffusivity in the RILF explained 31% of the variance in APA duration, while the RD in the CC Forceps Major explained 25% in step length in FoG.
Our study provides insights into the role of white matter tracts in modulating step initiation in FoG. These results align with previous research highlighting the importance of white matter tracts in motor control and cognitive functions in PD [
15,
16,
17]. Radial diffusivity (RD) refers to the water diffusion perpendicular to the main axis of the white matter fibers. Changes in RD provide insights into the microstructural integrity of white matter. Increased RD is generally interpreted as losing myelin integrity [
28] and is associated with motor and cognitive decline in PD [
29]. Specifically, increased RD in extensive white matter circuitry has been suggested to play a role in the pathophysiology of FoG [
30,
31,
32]. The results of Iseki et al. [
30] indicated a significant positive correlation between the nFoGq scores and RD values, suggesting that increases in RD reflect a degradation of myelin integrity. This increase in RD was interpreted as indicative of damage to both myelin and axons, which may contribute to the motor difficulties observed in PD, particularly those related to FoG. The analysis also demonstrated that changes in RD were more closely associated with FoG than with other cognitive or gait parameters, highlighting the specificity of these white matter abnormalities to the symptom of FoG in PD.
Step initiation impairments have also been linked with the pathophysiology of FoG. During step initiation, intact inhibitory control allows timely coordination between the interruption of the planning phase associated with the APA and the release of the step [
33]. This process is reflected in the APA duration. A longer time to transition from the automatic APA phase to the release of the voluntary step has been associated with FoG [
34,
35]. A noteworthy finding in our study is the association of the right hemisphere with increased APA duration in the FoG. Prior studies have also reported predominantly right hemisphere impairments in FoG patients, including white matter [
20,
36,
37,
38]. Functional brain imaging studies have shown that individuals with FoG often exhibit right-sided deficits [
39]. These regions are crucial for integrating sensory information necessary for gait control. Additionally, metabolic imaging studies have revealed lower glucose metabolism in the right parietal cortex [
39], suggesting impaired neural activity in this region among FoG. This right-hemisphere predominance highlights the importance of lateralized brain functions in the neural control of gait. It suggests that targeted interventions might need to consider these asymmetries to address FoG effectively.
The RILF’s association with APA duration in our study underscores its involvement in the preparatory phase of the movement, likely reflecting its role in integrating sensory and motor information crucial for step initiation. The RILF connects the occipital and temporal lobes and has been identified as one of the tracts involved in inhibitory circuitry [
40]. This finding aligns with existing literature that emphasizes the RILF’s connectivity with regions involved in visuospatial processing and inhibitory control, which are essential for planning and executing movements. Damage to the RILF can lead to difficulties integrating these sensory inputs, disrupting the initiation and smooth continuation of gait. This disruption is particularly evident in PD during movement initiation, where the inability to effectively process visual and spatial cues can result in FoG [
41]. Monaghan et al. [
42] investigated the association between white matter integrity and improvements in reactive stepping following a training program for individuals with PD. Improvements in step behavior were correlated with greater integrity of RILF and CC. Therefore, the RILF connection highlights how deficits in these areas can exacerbate FoG symptoms by impairing the coordination necessary for step initiation.
Regarding the CC, our results showed that those with smaller step lengths had higher RD in the FoG group. The CC’s role in step length regulation points to its importance in interhemispheric communication and coordination of bilateral motor activities. The step length depends on the integrity of the supporting and stepping leg coordination, which may be controlled by the CC [
43]. Dale et al. [
44] showed that lesions in the CC led to FoG, possibly due to disrupted callosal connections between the right and left pre-SMA. Accordingly, Wang et al. [
45] (2016) showed that FoG patients had more white matter abnormalities in the CC compared to controls. The correlation between reduced step length and RD in the FoG group underscores the importance of CC integrity in maintaining proper leg coordination.
Similar to our results, other studies found that RILF and CC were associated with gait disorders. Tan et al. [
46] showed that RILF and CC are more affected in individuals with PD and postural instability gait disorders (PIGD). Wen et al. [
47] found strong correlations of motor severity with higher RD of RILF and CC in PIGD compared to the tremor subtype of PD. Pietracupa [
37] showed that white matter changes in FoG compared with healthy subjects were prevalent in the RILF and the frontal radiations of the CC, and these abnormalities correlated with cognitive measures, including inhibitory control.
The current study’s results must be interpreted cautiously due to a relatively small sample size, limiting the findings’ generalizability. Also, the cross-sectional design limits the ability to infer causal relationships between white matter integrity and step initiation impairments. Although we controlled for some variables, we did not account for age, disease duration, L-dopa equivalent, and UPDRS, which differed between the groups. However, the FoG group participants were younger than those in the nFoG group, so the regression results cannot be attributed to brain aging. Despite the lack of covariate control, our results remain robust and trustworthy. The methodologies employed were rigorous, with precise measurements of white matter integrity and step initiation impairments, ensuring high reliability.
Our study also reinforces the potential utility of advanced neuroimaging techniques, such as DTI, in elucidating the microstructural changes associated with PD and FoG. DTI metrics provide a sensitive measure of white matter integrity that can reveal subtle but clinically significant changes in neural pathways. This approach not only enhances our understanding of the pathophysiological basis of FoG but also offers a valuable tool for monitoring disease progression and the effects of therapeutic interventions. By identifying specific biomarkers associated with motor impairments, we can improve the precision of clinical assessments and develop more targeted and effective treatment strategies.
Our findings are consistent with previous research, reinforcing their validity. This study underscores the critical role of white matter integrity in the pathophysiology of step initiation impairments in FoG. By demonstrating that increased RD in the RILF and CC significantly correlates with longer APA and reduced step length, respectively, the research highlights the importance of these white matter tracts in motor control. These findings provide valuable insights into the neural mechanisms underlying FoG, suggesting that targeting white matter integrity could be a potential therapeutic avenue. Furthermore, our results underscore the necessity for early detection of white matter deterioration in PD, which could allow for timely and more effective therapeutic interventions. By improving our understanding of the neural underpinnings of FoG, this study lays the groundwork for developing novel diagnostic tools and treatment modalities that specifically address the microstructural changes associated with this debilitating phenomenon. Future studies should focus on validating these findings across diverse populations and exploring the longitudinal impact of therapeutic interventions on white matter integrity and FoG symptoms.