- 1Institute of Biomedical and Health Engineering, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
- 2Department of Neurology, Guangzhou First People’s Hospital, School of Medicine, South China University of Technology, Guangzhou, China
- 3Department of Radiology, Guangzhou First People’s Hospital, School of Medicine, South China University of Technology, Guangzhou, China
- 4Department of Radiation Oncology, Southwestern Medical Center, University of Texas, Dallas, TX, United States
- 5GYENNO Technologies Co., Ltd., Shenzhen, China
Objective: The purposes of this study are to investigate the regional homogeneity (ReHo) of spontaneous brain activities in Parkinson’s disease (PD) patients with freeze of gait (FOG) and to investigate the neural correlation of movement function through resting-state functional magnetic resonance imaging (RS-fMRI).
Methods: A total of 35 normal controls (NC), 33 PD patients with FOG (FOG+), and 35 PD patients without FOG (FOG−) were enrolled. ReHo was applied to evaluate the regional synchronization of spontaneous brain activities. Analysis of covariance (ANCOVA) was performed on ReHo maps of the three groups, followed by post hoc two-sample t-tests between every two groups. Moreover, the ReHo signals of FOG+ and FOG− were extracted across the whole brain and correlated with movement scores (FOGQ, FOG questionnaire; GFQ, gait and falls questionnaire).
Results: Significant ReHo differences were observed in the left cerebrum. Compared to NC subjects, the ReHo of PD subjects was increased in the left angular gyrus (AG) and decreased in the left rolandic operculum/postcentral gyrus (Rol/PostC), left inferior opercular-frontal cortex, left middle occipital gyrus, and supramarginal gyrus (SMG). Compared to that of FOG−, the ReHo of FOG+ was increased in the left caudate and decreased in the left Rol/PostC. Within the significant regions, the ReHo of FOG+ was negatively correlated with FOGQ in the left SMG/PostC (r = −0.39, p < 0.05). Negative correlations were also observed between ReHo and GFQ/FOGQ (r = −0.36/−0.38, p < 0.05) in the left superior temporal gyrus (STG) of the whole brain analysis based on AAL templates.
Conclusion: The ReHo analysis suggested that the regional signal synchronization of brain activities in FOG+ subjects was most active in the left caudate and most hypoactive in the left Rol/PostC. It also indicated that ReHo in the left caudate and left Rol/PostC was critical for discriminating the three groups. The correlation between ReHo and movement scores (GFQ/FOGQ) in the STG has the potential to differentiate FOG+ from FOG−. This study provided new insight into the understanding of PD with and without FOG.
Introduction
Parkinson’s disease (PD) is a kind of neurodegenerative disease characterized by motor deficits (Villarreal et al., 2018). Freezing of gait (FOG) is a disabling symptom characterized by brief episodes of an inability to take a step or taking extremely short steps that typically occurs on initiating gait or on turning while walking (Nutt et al., 2011). Although FOG is more commonly observed in PD patients with advanced disease stages and old age, it may also occur in the early stage of idiopathic PD. Nearly 50% of PD patients suffer from FOG (Macht et al., 2007). Though FOG is transient and lasts for only a few seconds, it greatly impacts the quality of life of affected patients.
In the most recent decade, an increasing number of neuroimaging studies have focused on exploring the pathophysiology of PD patients with FOG by using different imaging modalities (Bartels and Leenders, 2008; Fasano et al., 2015) such as positron emission tomography (PET) (Park et al., 2009; Bohnen et al., 2014), functional near-infrared spectroscopy (fNIRS) (Maidan et al., 2015), diffusion tensor imaging (DTI) (Schweder et al., 2010; Herman et al., 2013), and functional magnetic resonance imaging (fMRI) (Tessitore et al., 2012b; Wang et al., 2016; Li et al., 2018). FOG in PD patients is suggested to be associated with abnormalities in motor, executive, cognitive, and affective functions (Amboni et al., 2008; Shine et al., 2013b). It has been synthesized that freezing occurs through a neural pathway in which the transient increase in inhibitory basal ganglia output leads to hypoactivity within the gait-coordinated brainstem, which may be caused by dopaminergic depletion in the striatum and hyperactivity in the subthalamic nucleus (Lewis and Shine, 2016). However, the pathophysiological mechanisms of FOG are not yet fully understood.
Resting-state fMRI (RS-fMRI) reflects alterations in spontaneous brain activities by measuring blood-oxygen-level-dependent (BOLD) signals (Fox and Raichle, 2007), and regional homogeneity (ReHo) evaluates signal synchronization by calculating the concordance of temporal change in BOLD signals within local brain regions (Zang et al., 2004). ReHo has been used to evaluate the differences between PD patients and normal controls (NC) (Li et al., 2016; Pan et al., 2017), and significant differences have been observed in the motor- and executive-related brain regions of PD patients, including the prefrontal cortex (Choe et al., 2013; Borroni et al., 2015), inferior parietal lobule (Wu et al., 2009), basal ganglia (Wang et al., 2018), precentral gyrus (Li et al., 2016; Wang et al., 2016), and cerebellum (Jiang et al., 2016). PD patients with gait disturbance are observed to have impaired coordination of movement and locomotion (Plotnik et al., 2008; Peterson et al., 2012). Therefore, it is hypothesized that PD with FOG (FOG+) may exhibit altered local signal synchronizations of neural activities in comparison to PD without FOG (FOG−) and NC. In this study, ReHo was employed to investigate the regional synchronizations of spontaneous brain activities in FOG+. Moreover, due to the altered movement performance in PD, the neural correlations of movement function in clinical assessment were also explored within the whole brain for both FOG+ and FOG−.
Materials and Methods
Subjects and Clinical Assessments
In this study, 37 NC and 72 subjects with a diagnosis of PD were investigated. The NC were healthy subjects with no history of neurological disease, no symptoms of PD, and no disorder of cognitive function. PD patients were diagnosed according to the clinical criteria of the Movement Disorder Society (Postuma et al., 2015). The exclusion criteria for PD patients: severe comorbidity disease (cardiovascular disease, respiratory disease, and malignant tumor, etc.), a history of surgical operations (thalamotomy and posteroventral pallidotomy, deep brain stimulation (DBS), and organ transplantation, etc.), or a pacemaker/metal implanted in their body, which is forbidden in MRI scanning. Among the 72 PD patients, 35 patients were included as FOG+ according to two criteria: (1) rating scores >0 for the third item in the freezing of gait questionnaire (FOGQ), which was described by Giladi et al. (2000) as “Do you feel that your feet get glued to the floor while walking, making a turn or when trying to initiate walking (freezing)?”; (2) based on the former criteria, OFF-FOG patients whose symptoms of FOG were improved after drug therapy were included. The other 37 patients were grouped as PD without freezing of gait (FOG−). In addition, FOGQ was employed to evaluate the severity of FOG performance (Giladi et al., 2000). Other clinical assessments were also made across all subjects with PD. Ratings on the Hoehn and Yahr system (H&Y) (Hoehn and Yahr, 1998) were collected to evaluate the severity of PD symptoms. The gait and falls questionnaire (GFQ) was applied to evaluate the gait and falls risk (Giladi et al., 2000). The motor part of the Unified Parkinson’s Disease Rating Scale (UPDRS-III) was also applied. For the motor assessments of PD patients, both FOGQ and GFQ assessments were made for the most severe OFF medication state based on their experience over the last week, and the UPDRS-III rating was assessed for the ON medication state. Non-motor symptoms of cognitive function were evaluated by montreal cognitive assessment (MOCA) (Nasreddine et al., 2005) and mini-mental state examination (MMSE) (Folstein et al., 1975). The levodopa equivalent daily dose (LEDD) of all of the PD patients was also collected. All subjects were recruited by the Guangzhou First People’s Hospital from May 2017 to September 2018.
Data Acquisition
All subjects (37 NC, 37 FOG−, and 35 FOG+) were enrolled in 3.0T SIEMENS MRI scanning and were required to lie quietly in the scanner, staying awake with eyes closed. All of the PD patients were in the ON medication state when the MRI scanning was performed. Both functional and structural MRI images were obtained. The RS-fMRI was obtained by echo-planar imaging (EPI) with the following parameters: repetition time (TR) = 2000 ms; echo time (TE) = 21 ms; slice thickness/gap = 4 mm/0.6 mm; acquisition matrix = 64 × 64; flip angle = 78°; in-plane resolution = 3.5 mm × 3.5 mm; FOV = 224 × 224 mm2. Sagittal T1-weighted images were obtained with the following parameters: TR/TE = 1900 ms/2.22 ms; acquisition matrix = 256 × 215; flip angle = 9°; in-plane resolution = 0.488 mm × 0.488 mm; slice thickness/gap = 1 mm/0.5 mm.
Data Preprocessing
The functional images were preprocessed using the toolkits of DPABI (Yan et al., 2016), the RS-fMRI Data Analysis Toolkit (REST)1, and Statistical Parametric Mapping (SPM12)2, implemented on a MATLAB platform. Data preprocessing included removal of the first 10 of the 220 time points in case of unstable signal quality, slice-timing adjustment (33 slices), head-motion correction, segmentation using a new segment (Ashburner and Friston, 2005) and diffeomorphic anatomical registration through Exponentiated Lie Algebra (DARTEL) (Ashburner, 2007), regression of nuisance covariates (including white matter, cerebrospinal fluid, and Friston’s 24 parameters of head motion) (Friston et al., 1996; Satterthwaite et al., 2013; Yan et al., 2013), spatial normalization to Montreal Neurological Institute (MNI) space by resampling to 3mm × 3mm × 3mm by DARTEL (Ashburner, 2007), a temporal filter with a bandpass of 0.01–0.1 Hz, and removal of linear detrending. Six parameters of head motion (three directions each of rotation and translation) were recorded during the scanning. Subjects with maximal translations exceeding 2.5 mm or rotations over 2.5 degrees were excluded. According to this exclusion criterion, a total of six subjects were excluded from three groups, leaving 35/35/33 subjects for NC/FOG−/FOG+, respectively. Additionally, the mean frame-wise displacement (FD) (Jenkinson et al., 2002) was calculated, which represents the head motion. The mean FD was added as a covariate in the statistical analysis.
Regional Homogeneity
ReHo evaluates local signal synchronizations by assessing the similarity between the time series of a chosen voxel and those of its neighboring voxels, and Kendall’s coefficient concordance (KCC) is applied to ReHo calculation between a voxel and its 26 neighboring voxels (Zang et al., 2004). KCC-ReHo is a value between 0 and 1. Higher values indicate better local synchronization. Voxel-wise ReHo maps were calculated, and the ReHo maps were then spatially smoothed with a full width at half maximum (FWHM) of 4 mm. Additionally, Z-transformation was applied to the ReHo maps for standardization by subtracting the global mean value and then dividing by the global standard deviation. The standardized ReHo Z-maps were applied to the subsequent statistical and correlative analysis.
Statistical and Correlative Analysis
Analysis of covariance (ANCOVA) was applied to explore the ReHo differences among NC, FOG−, and FOG+, with age, sex, and mean FD Jenkinson as covariates. The resultant F-map was corrected by multiple comparisons of the Gaussian Random Field (GRF) with voxel p < 0.05 and cluster p < 0.05 within a gray matter mask, two-tailed (F > 3.83 and cluster size >6750 mm3). The surviving voxels were then extracted as a mask in the post hoc analysis of ReHo differences between every two groups by two-sample t-test with the covariates of age, sex, and mean FD. For the ReHo differences between FOG+ and FOG−, the clinical assessments that demonstrated significant group difference were also controlled as covariates. The resultant T-maps were further corrected by GRF with voxel p < 0.001 and cluster p < 0.05, which is beneficial for avoiding false positives (Woo et al., 2014).
In addition to the whole-brain gray matter, ReHo differences were also investigated within certain brain regions that have frequently been reported in previous motor- and gait-related PD studies, including the basal ganglia (caudate, putamen, and pallidum), sensorimotor cortices, cerebellum, hippocampus, para-hippocampus, and fusiform gyrus (Camicioli et al., 2003; Herman et al., 2014; Pan et al., 2017; Wang et al., 2018; Li et al., 2019). These brain regions were, respectively, extracted from the Automated Anatomical Labeling (AAL) template, which contains 116 brain regions, including 90 cerebrum regions and 26 cerebellum regions (Tzourio-Mazoyer et al., 2002). The corresponding AAL atlas regions were: the bilateral precentral gyrus (PreC, AAL-1,2), postcentral gyrus (PostC, AAL-57,58), supplementary motor area (SMA, AAL-19,20), cerebellum (AAL-91 to 116), caudate (AAL-71,72), putamen (AAL-73,74), pallidum (AAL-75,76), hippocampus (AAL-37,38), para-hippocampus (AAL-39,40), and fusiform gyrus (AAL-55,56). Both ANCOVA and post hoc analysis were performed within these regions, respectively, with covariates as in the analysis within whole-brain gray matter. Both ANCOVA and post hoc two-sample t-test were corrected by GRF with voxel p < 0.05 and cluster p < 0.05, two-tailed.
Brain regions showing significant differences were extracted as regions of interest (ROIs) for exploring the correlative relationship between signal synchronization and movement function (FOGQ, GFQ). The ReHo value was extracted from ROIs by averaging the values of all voxels within ROI. The Pearson correlation coefficient (statistical significance level p < 0.05) was used to quantify the correlation between ReHo and FOGQ/GFQ. Moreover, correlations were also analyzed within AAL templates to examine the neural interactions between regional signal synchronization and movement function across the whole brain.
Results
Demographic Characteristics and Clinical Assessments
After exclusion of subjects with excessive head motion, 35 NC, 35 FOG−, and 33 FOG+ remained. There was no significant difference in head motion (mean FD) among the three groups (p = 0.2096). Table 1 summarizes the demographic characteristics and clinical assessments of the subjects. A significant difference in age was observed among the three groups (p = 0.0001). Specifically, FOG+ showed higher ages than FOG− and NC, while no age difference was found between FOG− and NC (p > 0.05). Compared to FOG−, FOG+ had experienced longer disease durations and had higher severity of PD symptoms (H&Y scores) and higher ratings for GFQ and FOGQ. However, FOG+ and FOG− demonstrated no significant differences (p > 0.05) on UPDRS-III, MMSE, and MOCA. Note that some PD subjects refused to answer the GFQ and FOGQ and thereby, 35/32 GFQ results remained for FOG−/FOG+ and 34/31 FOGQ results for FOG−/FOG+, respectively.
Group Differences of Regional Homogeneity
The ACONVA analysis demonstrated significant ReHo differences among NC, FOG−, and FOG+. The identification of significant brain regions was based on XjView3. It was observed that all significant regions were located in the left cerebrum, including the sensorimotor area of PostC and rolandic operculum (Rol), the posterior middle temporal gyrus (MTG), the parietal gyrus with the angular gyrus (AG) and supramarginal gyrus (SMG), the middle occipital gyrus (MOG), and the opercular part of the inferior frontal gyrus (IFGoper) (Figure 1). Within these significant regions, post hoc analysis of two-sample t-test was performed between every two groups, corrected by GRF with voxel p < 0.001 and cluster p < 0.05, two-tailed (T > 3.45, cluster size >135 mm3). The T-map between FOG+ and FOG− was thresholded with an uncorrected p < 0.05 (Figure 2C), which failed to survive under GRF correction. The results of between-group ReHo differences are shown in Table 2 and Figure 2. Both FOG+ and FOG− showed higher ReHo than the NC in the left AG, and lower ReHo in the left IFGoper, left MOG, left Rol/PostC, and left SMG/PostC (Table 2 and Figures 2A,B). Moreover, FOG+ achieved higher ReHo than the FOG− in the left MTG and left AG (Table 2 and Figure 2A) and lower ReHo than the FOG− in the left Rol/PostC (Table 2 and Figure 2C).
Figure 1. ACONVA analysis of regional homogeneity (ReHo) among NC, FOG–, and FOG+ within gray matter. The F-map was corrected by GRF with voxel p < 0.05 and cluster p < 0.05 (F > 3.83, cluster size >6750 mm3) within a gray matter mask. L/R, left/right hemisphere; NC, normal controls; FOG+/FOG–, Parkinson’s disease with/without freezing of gait.
Figure 2. Post hoc analysis of ReHo differences among NC, FOG–, and FOG+. (A) ReHo differences between FOG+ and NC. (B) ReHo differences between FOG– and NC. (C) ReHo differences between FOG+ and FOG– (voxel p < 0.05, uncorrected). T-maps (A,B) were corrected by GRF with voxel p < 0.001 and cluster p < 0.05 (T > 3.45, cluster size >135 mm3) within the mask of significant brain regions identified by ANCOVA. ReHo signals were extracted from the left MTG/AG (D–F) and left Rol/PostC (G–I), which were identified by post hoc analysis. L/R, left/right hemisphere; NC, normal controls; FOG+/FOG–, Parkinson’s disease with/without freezing of gait; MTG, middle temporal gyrus; AG, angular gyrus; Rol, rolandic operculum; PostC, postcentral gyrus; ∗∗p < 0.005; ∗∗∗p < 0.0005.
Table 2. Brain regions showing significant regional homogeneity (ReHo) differences among NC, FOG−, and FOG+ within gray matter.
Furthermore, ReHo differences were also examined in the bilateral PreC, PostC, SMA, cerebellum, caudate, putamen, pallidum, hippocampus, para-hippocampus, and fusiform gyrus, respectively. However, significant results were only observed in the caudate (Table 3 and Figure 3). The F-map acquired by ANCOVA among the three groups was corrected by GRF with voxel p < 0.05 and cluster p < 0.05 (threshold of F > 3.84 and cluster size >1566 mm3) (Figure 3B). ReHo differences were demonstrated between the three groups in the left caudate (Figures 3B,C). Post hoc two-sample t-testing on the pair groups were performed within the significant left caudate, and the resultant T-maps were corrected by GRF with voxel p < 0.05 and cluster p < 0.05 (threshold of T > 2.00 and cluster size >486 mm3). FOG+ was observed to have higher ReHo than FOG− (Figures 3F,I) and NC (Figures 3D,G) in the left caudate. While no cluster survived the GRF correction for the T-map between FOG− and NC, thresholding with an uncorrected voxel p < 0.05 and cluster size of 351 mm3 demonstrated higher ReHo in the left caudate in FOG− than in NC (Figures 3E,H).
Figure 3. Regional homogeneity differences among NC, FOG–, and FOG+ within the bilateral caudate. (A) Bilateral caudate. (B) F-map of ANCOVA analysis of the three groups, corrected by GRF with voxel p < 0.05, threshold of F > 3.84 and cluster size >1566 mm3. (C) FOG+ subjects showed higher ReHo within the left caudate than did FOG– and NC. Post hoc two-sample t-tests on FOG+ and NC (D,G), FOG– and NC (E,H), and FOG+ and FOG– (F,I). The T-map in (E) was with a threshold of an uncorrected voxel p < 0.05, (D,F) were corrected by GRF with voxel p < 0.05, threshold of T > 2.00 and cluster size >486 mm3. L/R, eft/right hemisphere; NC, normal controls; FOG+/FOG–, Parkinson’s disease with/without freezing of gait.
Table 3. Brain regions showing significant ReHo differences among NC, FOG−, and FOG+ within the bilateral caudate.
Correlative Analysis
The correlative analysis between local signal synchronization (ReHo) and movement function (GFQ/FOGQ) was performed within significant brain regions identified by post hoc two-sample t-tests. Among the regions, significant results were observed in the left SMG/PostC (Figure 4A), which was visualized by BrainNet Viewer (Xia et al., 2013). The ReHo values of FOG+ and FOG− were significantly lower than those of NC (p < 0.0001), while no difference was found between those of FOG+ and FOG− (p = 0.2408) (Figure 4B). The ReHo of FOG+ was negatively correlated with FOGQ in the left SMG/PostC (r = −0.39, p < 0.05) (Figure 4D), while no significant correlation was found between the ReHo of FOG− and FOGQ (Figure 4C). In addition, neither the ReHo of FOG+ nor the ReHo of FOG− was correlated with GFQ.
Figure 4. Correlation between ReHo and movement function in the left supramarginal gyrus/postcentral gyrus (SMG/PostC). (A) Visualization of left SMG/PostC. (B) ReHo of left SMG/PostC. ∗∗∗p < 0.0001 between FOG+/FOG– and NC, p = 0.2408 between FOG+ and FOG–. (C) Correlation between ReHo of FOG– and movement function. (D) Correlation between ReHo of FOG+ and movement function.
As well as from the ROIs identified by the post hoc analysis (Tables 2, 3), ReHo signals were also extracted from AAL templates and correlated with GFQ and FOGQ, respectively. Significant results were observed in the left superior temporal gyrus (STG, AAL-81), left orbital part of the medial prefrontal cortex (MPFCorb, AAL-25), left gyrus rectus (AAL-27), right inferior temporal gyrus (ITG, AAL-90), and left cerebellum anterior lobe (CAL, AAL-97) (Figure 5). The ReHo values of FOG+ in the left STG were observed to have negative correlations (p < 0.05) with GFQ (r = −0.36) and FOGQ (r = −0.38) (Figure 5A). The ReHo values of FOG−, meanwhile, were observed to have negative correlations with GFQ/FOGQ in the left MPFCorb (Figure 5B), left gyrus rectus (Figure 5C), and right ITG (Figure 5D) and to be positively correlated with GFQ/FOGQ in the left CAL (Figure 5E). Though significant correlations were observed in the above five brain regions, the ReHo signals extracted from the five sub-templates showed no significant difference (p > 0.05) between FOG− and FOG+ according to two-sample t-test (Figure 5).
Figure 5. Correlative analysis between ReHo and movement function within automated anatomical labeling (AAL) templates. Significant correlation was found in the (A) left superior temporal gyrus, (B) left orbital part of the medial prefrontal cortex, (C) left gyrus rectus, (D) right inferior temporal gyrus, and (E) left cerebellum anterior lobe. FOGQ, freezing of gait questionnaire; GFQ, gait and falls questionnaire. Statistical significance notation: ∗p < 0.05, ∗∗∗p < 0.005.
Discussion
In this study, we assessed ReHo to explore differences in local synchronization among NC, FOG+, and FOG− subjects, and significant results were observed in the left cerebrum. In comparison to NC, FOG+, and FOG− showed increased ReHo in the left AG but decreased ReHo in the left IFGoper, left Rol/PostC, left MOG, and left SMG (Table 2 and Figures 2A,B). Compared to FOG−, the ReHo in FOG+ was decreased in the left Rol/PostC (Figure 2C) and increased in the left caudate (Figure 3). Neural correlation analysis showed that the ReHo of FOG+ was negatively correlated with FOGQ (r = −0.39) in the left SMG/PostC (Figure 4) and negatively correlated with GFQ/FOGQ (r = −0.36/−0.38) in the STG (Figure 5A).
Previous studies have reported structural changes in FOG+ vs. FOG−, with gray matter volume reductions in the posterior cingulate cortex/precuneus (PCC/PCu) (Tessitore et al., 2012a), left IFG (Pan et al., 2012), inferior parietal lobule (IPL) (Kostić et al., 2012), and AG (Herman et al., 2014). The PCC/PCu, medial prefrontal cortex (MPFC), AG, and posterior MTG are major parts of the default mode network (DMN), where brain activities are active in the resting state but passive in the task-induced state (Raichle et al., 2001). Gray-matter atrophy of the DMN in FOG+ probably leads to dysfunctional brain activities in the resting state. Meta-analysis of the ReHo of PD has been reported to give consistent results in the DMN and motor networks, with increased ReHo in the bilateral IPL, bilateral MPFC, and left cerebellum and decreased ReHo in the right putamen and right PreC (Pan et al., 2017; Wang et al., 2018). Our findings are consistent with the results showing increased ReHo in the DMN, observing increased ReHo in the left MTG/AG (Figures 2A,B,D,E). However, no gait-specificity of ReHo in the DMN was found in PD patients. Under a threshold of an uncorrected voxel p < 0.05, FOG+ were observed to have increased ReHo in the left MTG/AG (Figure 2C), while no difference was found between FOG+ and FOG− when the averaged ReHo was extracted from these regions (Figure 2F). A previous study on FOG+ also reported decreased ReHo in the frontal cortex and motor area (Zhou et al., 2018), which suggests that FOG+ subjects have decreased cognitive function and motor function. Our findings also demonstrated decreased ReHo in the IFG and MOG. A resting-state study with connectivity analysis of FOG+ observed decreased connectivity in the frontal and occipital lobes, which correspond to the executive and visual networks (Tessitore et al., 2012b). These findings suggest the disruption of executive and visual functions in FOG+.
One of the most critical findings of this study is that for ReHo in the left Rol/PostC, FOG+ <FOG−<NC, which suggests that this region is most hypoactive in FOG+ (Figures 2G–I). However, when using a threshold of an uncorrected voxel p < 0.05 when comparing the ReHo of FOG+ and FOG−, significant difference (p < 0.005) was observed in the left Rol/PostC by extracting the averaged ReHo from the surviving voxels (Figure 2I). This suggests that ReHo values in the left Rol/PostC, the sensorimotor areas, are critical features for discriminating the three groups. One interesting finding comes from the fact that all of the significant ROIs among FOG+, FOG−, and NC were located in the left hemisphere of the cerebrum (Figures 1,2 and Table 2). An fMRI study of emotional picture stimuli reported activations of left-brain activities in response to positive pictures and hence suggested that the left hemisphere of the brain is associated with positive emotions (Canli et al., 1998). Altered ReHo in the left hemisphere may reveal decreased positive emotion in FOG+. Previous studies also demonstrate that emotional state affects the motor control of gait (Naugle et al., 2011). Depression and anxiety are also the major emotional symptoms in PD, which exacerbate poor motor performance (Macht et al., 2005; Avanzino et al., 2018). Therefore, cognitive behavioral therapy is recommended for the treatment of patients with FOG (Berardelli et al., 2015).
The basal ganglia, sensorimotor area, cerebellum, hippocampus, parahippocampus, and fusiform gyrus are critical regions involved in the pathophysiology of PD. ReHo differences were also examined within these brain regions. However, significant results were only observed in the left caudate: FOG+ >FOG−>NC (Table 3 and Figure 3). Gray matter atrophy in the left caudate (Jia et al., 2019) and consistently decreased ReHo in the putamen (Pan et al., 2017; Wang et al., 2018) have been observed in PD patients. In our results, meanwhile, PD showed increased ReHo in the left caudate, with FOG+ being the most active. Decreased ReHo may indicate functional deficits caused by diseases, while increased ReHo may be related to a compensatory mechanism for maintaining normal function (Pan et al., 2017). It is suggested that FOG in PD is associated with functional decoupling between the cognitive control network and the basal ganglia (Shine et al., 2013a). Increased ReHo in the left caudate may reflect a compensation of the cognitive control function in PD patients.
Within the significant brain regions, the ReHo of FOG+ was negatively correlated (r = −0.39) with FOGQ in the left SMG/PostC (Figure 4). Negative correlation was also observed between ReHo and GFQ/FOGQ (r = −0.36/−0.38) in the STG based on the whole-brain analysis (Figure 5A). A DBS study on the effect of sub-thalamic nucleus stimulation in PD reported a positive correlation between motor scores and metabolic activity in parietal-temporal sensory-related brain areas (Chul et al., 2007). Activation of brain activity was found in the temporal lobe during a memory paradigm fMRI study after rehabilitation (Díez-Cirarda et al., 2017). These findings indicate that brain activities in the motor area and temporal lobe are associated with an improvement in motor performance.
There are two limitations to this study. The demographic characteristics of age and sex are not properly matched among FOG+, FOG−, and NC. The age of FOG+ is significantly higher than that of FOG− and NC. PD patients with advanced age and a higher stage of disease progression are more likely to experience FOG (Zhang et al., 2016). In this study, age and sex were included as covariates in the statistical analysis to regress out the unmatched confounds. The second limitation comes from the insufficient sample size. In our future work, more subjects will be enrolled for a better understanding of the neuroimaging features of PD patients with freezing of gait.
Conclusion
ReHo was used to explore the regional signal synchronization of brain activities in PD patients. The results suggest that the brain activities of PD patients with FOG were the most active in the left caudate and the most hypoactive in the left Rol/PostC. The correlation analysis between ReHo and movement scores (GFQ/FOGQ) in the left STG provides the potential to stratify PD patients with and without FOG. This study provides new insight for understanding PD patients.
Data Availability Statement
The datasets generated for this study are available on request to the corresponding author.
Ethics Statement
This study was approved by the Institutional Review Board (IRB) of Guangzhou First People’s Hospital. Written informed consent was obtained from all subjects.
Author Contributions
YL wrote the manuscript. YL, ML, and XW conceived of the idea and performed the literature review. YL, ML, XW, HC, GH, and SY performed the data analysis. XW, XR, ZLuo, and JZ contributed to the data collection. All authors interpreted the results, reviewed the manuscript and joined the discussion of the manuscript.
Funding
This work was supported in part by grants from the Shenzhen Matching Project (GJHS20170314155751703), the National Key Research and Development Program of China (2016YFC0105102), the National Science Foundation of China (61871374 and 81871846), the Science Foundation of Guangdong (2017B020229002, 2015B020233011, and 2014A030312006), the Leading Talent of Special Support Project in Guangdong (2016TX03R139), the Science and Technology Planning Project of Guangzhou (201804010032), the Guangzhou Municipal Health Bureau Project (20171A010247), the Guangzhou Key Project of R&D Innovation (2016201604030018), the Shenzhen Key Technical Research Project (JSGG20160229203812944), the Shenzhen Basic Technology Research Project (JCYJ20170413162458312), and the Shenzhen Engineering Laboratory for Key Technologies on Intervention Diagnosis and Treatment Integration.
Conflict of Interest
ZLuo was employed by the company of GYENNO Technologies Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Footnotes
- ^ http://www.restfmri.net/forum/REST
- ^ http://www.fil.ion.ucl.ac.uk/spm/
- ^ http://www.alivelearn.net/xjview/
References
Amboni, M., Cozzolino, A., Longo, K., Picillo, M., and Barone, P. (2008). Freezing of gait and executive functions in patients with Parkinson’s disease. Mov. Disord. 23, 395–400. doi: 10.1002/mds.21850
Ashburner, J. (2007). A fast diffeomorphic image registration algorithm. Neuroimage 38, 95–113. doi: 10.1016/j.neuroimage.2007.07.007
Ashburner, J., and Friston, K. J. (2005). Unified segmentation. Neuroimage 26, 839–851. doi: 10.1016/j.neuroimage.2005.02.018
Avanzino, L., Lagravinese, G., Abbruzzese, G., and Pelosin, E. (2018). Relationships between gait and emotion in Parkinson’s disease: a narrative review. Gait Posture 65, 57–64. doi: 10.1016/j.gaitpost.2018.06.171
Bartels, A. L., and Leenders, K. L. (2008). Brain imaging in patients with freezing of gait. Mov. Disord. 23, S461–S467. doi: 10.1002/mds.21912
Berardelli, I., Fabbrini, G., Pasquini, M., Biondi, M., Berardelli, A., and Roselli, V. (2015). Cognitive behavioral therapy in movement disorders: a review. Mov. Disord. Clin. Pract. 2, 107–115. doi: 10.1002/mdc3.12160
Bohnen, N. I., Frey, K. A., Studenski, S., Kotagal, V., Koeppe, R. A., Constantine, G. M., et al. (2014). Extra-nigral pathological conditions are common in Parkinson’s disease with freezing of gait: an in vivo positron emission tomography study. Mov. Disord. 29, 1118–1124. doi: 10.1002/mds.25929
Borroni, B., Premi, E., Formenti, A., Turrone, R., Alberici, A., Cottini, E., et al. (2015). Structural and functional imaging study in dementia with lewy bodies and Parkinson’s disease dementia. Park. Relat. Disord. 21, 1049–1055. doi: 10.1016/j.parkreldis.2015.06.013
Camicioli, R., Moore, M. M., Kinney, A., Corbridge, E., Glassberg, K., and Kaye, J. A. (2003). Parkinson’s disease is associated with hippocampal atrophy. Mov. Disord. 18, 784–790. doi: 10.1002/mds.10444
Canli, T., Desmond, J. E., Zhao, Z., Glover, G., and Gabrieli, J. D. E. (1998). Hemispheric asymmetry for emotional stimuli detected with fMRI. Neuroreport 9, 3233–3239. doi: 10.1097/00001756-199810050-199810019
Choe, I. H., Yeo, S., Chung, K. C., Kim, S. H., and Lim, S. (2013). Decreased and increased cerebral regional homogeneity in early Parkinson’s disease. Brain Res. 1527, 230–237. doi: 10.1016/j.brainres.2013.06.027
Chul, H. L., Aalto, S., Rinne, J. O., Ki, O. L., Seung, H. O., Jin, W. C., et al. (2007). Different cerebral cortical areas influence the effect of subthalamic nucleus stimulation on Parkinsonian motor deficits and freezing of gait. Mov. Disord. 22, 2176–2182. doi: 10.1002/mds.21609
Díez-Cirarda, M., Ojeda, N., Peña, J., Cabrera-Zubizarreta, A., Lucas-Jiménez, O., Gómez-Esteban, J. C., et al. (2017). Increased brain connectivity and activation after cognitive rehabilitation in Parkinson’s disease: a randomized controlled trial. Brain Imaging Behav. 11, 1640–1651. doi: 10.1007/s11682-016-9639-x
Fasano, A., Herman, T., Tessitore, A., Strafella, A. P., and Bohnen, N. I. (2015). Neuroimaging of freezing of gait. J. Parkinsons Dis. 5, 241–254. doi: 10.3233/JPD-150536
Folstein, M. F., Folstein, S. E., and McHugh, P. R. (1975). Mini-mental state”. a practical method for grading the cognitive state of patients for the clinician. J. Psychiatr. Res. 12, 189–198. doi: 10.1016/0022-3956(75)90026-90026
Fox, M. D., and Raichle, M. E. (2007). Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat. Rev. Neurosci. 8, 700–711. doi: 10.1038/nrn2201
Friston, K. J., Williams, S., Howard, R., Frackowiak, R. S. J., and Turner, R. (1996). Movement-related effects in fMRI time-series. Magn. Reson. Med. 35, 346–355. doi: 10.1002/mrm.1910350312
Giladi, N., Shabtai, H., Simon, E. S., Biran, S., Tal, J., and Korczyn, A. D. (2000). Construction of freezing of gait questionnaire for patients with Parkinsonism. Park. Relat. Disord. 6, 165–170. doi: 10.1016/S1353-8020(99)00062-60
Herman, T., Rosenberg-Katz, K., Jacob, Y., Auriel, E., Gurevich, T., Giladi, N., et al. (2013). White matter hyperintensities in Parkinson’s disease: do they explain the disparity between the postural instability gait difficulty and tremor dominant subtypes? PLoS One 8:e55193. doi: 10.1371/journal.pone.0055193
Herman, T., Rosenberg-Katz, K., Jacob, Y., Giladi, N., and Hausdorff, J. M. (2014). Gray matter atrophy and freezing of gait in Parkinson’s disease: is the evidence black-on-white? Mov. Disord. 29, 134–139. doi: 10.1002/mds.25697
Hoehn, M. M., and Yahr, M. D. (1998). Parkinsonism: onset, progression, and mortality. Neurology 50, 318–318. doi: 10.1212/WNL.50.2.318
Jenkinson, M., Bannister, P., Brady, M., and Smith, S. (2002). Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage 17, 825–841. doi: 10.1016/S1053-8119(02)91132-91138
Jia, X., Li, K., Fu, X., Liang, P., and Li, Y. (2019). Precuneus dysfunction in Parkinson’s disease with mild cognitive impairment. Front. Aging Neurosci. 10:427. doi: 10.3389/fnagi.2018.00427
Jiang, S., Wang, M., Zhang, L., Yuan, Y., Tong, Q., Ding, J., et al. (2016). Regional homogeneity alterations differentiate between tremor dominant and postural instability gait difficulty subtypes of Parkinson’s disease. J. Neural. Transm. 123, 219–229. doi: 10.1007/s00702-015-1490-1495
Kostić, V. S., Agosta, F., Pievani, M., Stefanova, E., Ječmenica-Lukić, M., Scarale, A., et al. (2012). Pattern of brain tissue loss associated with freezing of gait in Parkinson disease. Neurology 78, 409–416. doi: 10.1212/WNL.0b013e318245d23c
Lewis, S. J. G., and Shine, J. M. (2016). The next step: a common neural mechanism for freezing of gait. Neuroscientist 22, 72–82. doi: 10.1177/1073858414559101
Li, J., Yuan, Y., Wang, M., Zhang, J., Zhang, L., Jiang, S., et al. (2018). Decreased interhemispheric homotopic connectivity in Parkinson’s disease patients with freezing of gait: a resting state fMRI study. Park. Relat. Disord. 52, 30–36. doi: 10.1016/j.parkreldis.2018.03.015
Li, M., Liu, Y., Chen, H., Hu, G., Yu, S., Ruan, X., et al. (2019). Altered global synchronizations in patients with Parkinson’s disease: a resting-state fMRI study. Front. Aging Neurosci. 11:139. doi: 10.3389/fnagi.2019.00139
Li, Y., Liang, P., Jia, X., and Li, K. (2016). Abnormal regional homogeneity in Parkinson’s disease: a resting state fMRI study. Clin. Radiol. 71, 28–34. doi: 10.1016/j.crad.2015.10.006
Macht, M., Kaussner, Y., Möller, J. C., Stiasny-Kolster, K., Eggert, K. M., Krüger, H. P., et al. (2007). Predictors of freezing in Parkinson’s disease: a survey of 6,620 patients. Mov. Disord. 22, 953–956. doi: 10.1002/mds.21458
Macht, M., Schwarz, R., and Ellgring, H. (2005). Patterns of psychological problems in Parkinson’s disease. Acta Neurol. Scand. 111, 95–101. doi: 10.1111/j.1600-0404.2005.00375.x
Maidan, I., Bernad-Elazari, H., Gazit, E., Giladi, N., Hausdorff, J. M., and Mirelman, A. (2015). Changes in oxygenated hemoglobin link freezing of gait to frontal activation in patients with Parkinson disease: an fNIRS study of transient motor-cognitive failures. J. Neurol. 262, 899–908. doi: 10.1007/s00415-015-7650-7656
Nasreddine, Z. S., Phillips, N. A., Bédirian, V., Charbonneau, S., Whitehead, V., Collin, I., et al. (2005). The montreal cognitive assessment, MoCA: a brief screening tool for mild cognitive impairment. J. Am. Geriatr. Soc. 53, 695–699. doi: 10.1111/j.1532-5415.2005.53221.x
Naugle, K. M., Hass, C. J., Joyner, J., Coombes, S. A., and Janelle, C. M. (2011). Emotional state affects the initiation of forward gait. Emotion 11, 267–277. doi: 10.1037/a0022577
Nutt, J. G., Bloem, B. R., Giladi, N., Hallett, M., Horak, F. B., and Nieuwboer, A. (2011). Freezing of gait: moving forward on a mysterious clinical phenomenon. Lancet Neurol. 10, 734–744. doi: 10.1016/S1474-4422(11)70143-70140
Pan, P. L., Song, W., and Shang, H. F. (2012). Voxel-wise meta-analysis of gray matter abnormalities in idiopathic Parkinson’s disease. Eur. J. Neurol. 19, 199–206. doi: 10.1111/j.1468-1331.2011.03474.x
Pan, P. L., Zhan, H., Xia, M. X., Zhang, Y., Guan, D. N., and Xu, Y. (2017). Aberrant regional homogeneity in Parkinson’s disease: a voxel-wise meta-analysis of resting-state functional magnetic resonance imaging studies. Neurosci. Biobehav. Rev. 72, 223–231. doi: 10.1016/j.neubiorev.2016.11.018
Park, H. K., Kim, J. S., Im, K. C., Oh, S. J., Kim, M. J., Lee, J. H., et al. (2009). Functional brain imaging in pure Akinesia with gait freezing: [18F] FDG PET and [18F] FP-CIT PET analyses. Mov. Disord. 24, 237–245. doi: 10.1002/mds.22347
Peterson, D. S., Plotnik, M., Hausdorff, J. M., and Earhart, G. M. (2012). Evidence for a relationship between bilateral coordination during complex gait tasks and freezing of gait in Parkinson’s disease. Park. Relat. Disord. 18, 1022–1026. doi: 10.1016/j.parkreldis.2012.05.019
Plotnik, M., Giladi, N., and Hausdorff, J. M. (2008). Bilateral coordination of walking and freezing of gait in Parkinson’s disease. Eur. J. Neurosci. 27, 1999–2006. doi: 10.1111/j.1460-9568.2008.06167.x
Postuma, R. B., Berg, D., Stern, M., Poewe, W., Olanow, C. W., Oertel, W., et al. (2015). MDS clinical diagnostic criteria for Parkinson’s disease. Mov. Disord. 30, 1591–1601. doi: 10.1002/mds.26424
Raichle, M. E., MacLeod, A. M., Snyder, A. Z., Powers, W. J., Gusnard, D. A., and Shulman, G. L. (2001). A default mode of brain function. Proc. Natl. Acad. Sci. U.S.A. 98, 676–682. doi: 10.1073/pnas.98.2.676
Satterthwaite, T. D., Elliott, M. A., Gerraty, R. T., Ruparel, K., Loughead, J., Calkins, M. E., et al. (2013). An improved framework for confound regression and filtering for control of motion artifact in the preprocessing of resting-state functional connectivity data. Neuroimage 64, 240–256. doi: 10.1016/j.neuroimage.2012.08.052
Schweder, P. M., Hansen, P. C., Green, A. L., Quaghebeur, G., Stein, J., and Aziz, T. Z. (2010). Connectivity of the pedunculopontine nucleus in parkinsonian freezing of gait. Neuroreport 21, 914–916. doi: 10.1097/WNR.0b013e32833ce5f1
Shine, J. M., Matar, E., Ward, P. B., Frank, M. J., Moustafa, A. A., Pearson, M., et al. (2013a). Freezing of gait in Parkinson’s disease is associated with functional decoupling between the cognitive control network and the basal ganglia. Brain 136, 3671–3681. doi: 10.1093/brain/awt272
Shine, J. M., Naismith, S. L., and Lewis, S. J. G. (2013b). The differential yet concurrent contributions of motor, cognitive and affective disturbance to freezing of gait in Parkinson’s disease. Clin. Neurol. Neurosurg. 115, 542–545. doi: 10.1016/j.clineuro.2012.06.027
Tessitore, A., Amboni, M., Cirillo, G., Corbo, D., Picillo, M., Russo, A., et al. (2012a). Regional gray matter atrophy in patients with Parkinson disease and freezing of gait. Am. J. Neuroradiol. 33, 1804–1809. doi: 10.3174/ajnr.A3066
Tessitore, A., Amboni, M., Esposito, F., Russo, A., Picillo, M., Marcuccio, L., et al. (2012b). Resting-state brain connectivity in patients with Parkinson’s disease and freezing of gait. Park. Relat. Disord. 18, 781–787. doi: 10.1016/j.parkreldis.2012.03.018
Tzourio-Mazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., et al. (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 15, 273–289. doi: 10.1006/nimg.2001.0978
Villarreal, M. F., Huerta-Gutierrez, R., and Fregni, F. (2018). “Parkinson’s disease,” in Clinical Trials in Neurology, Vol. 138, ed. F. Fregni, (New York, NY: Humana Press).
Wang, J., Zhang, J. R., Zang, Y. F., and Wu, T. (2018). Consistent decreased activity in the putamen in Parkinson’s disease: a meta-analysis and an independent validation of resting-state fMRI. Gigascience 7, 1–13. doi: 10.1093/gigascience/giy071
Wang, M., Jiang, S., Yuan, Y., Zhang, L., Ding, J., Wang, J., et al. (2016). Alterations of functional and structural connectivity of freezing of gait in Parkinson’s disease. J. Neurol. 263, 1583–1592. doi: 10.1007/s00415-016-8174-8174
Woo, C. W., Krishnan, A., and Wager, T. D. (2014). Cluster-extent based thresholding in fMRI analyses: pitfalls and recommendations. Neuroimage 91, 412–419. doi: 10.1016/j.neuroimage.2013.12.058
Wu, T., Long, X., Zang, Y., Wang, L., Hallett, M., Li, K., et al. (2009). Regional homogeneity changes in patients with parkinson’s disease. Hum. Brain Mapp. 30, 1502–1510. doi: 10.1002/hbm.20622
Xia, M., Wang, J., and He, Y. (2013). BrainNet viewer: a network visualization tool for human brain connectomics. PLoS One 8:e68910. doi: 10.1371/journal.pone.0068910
Yan, C. G., Cheung, B., Kelly, C., Colcombe, S., Craddock, R. C., Di Martino, A., et al. (2013). A comprehensive assessment of regional variation in the impact of head micromovements on functional connectomics. Neuroimage 76, 183–201. doi: 10.1016/j.neuroimage.2013.03.004
Yan, C. G., Wang, X., Zuo, X. N., and Zang, Y. F. (2016). DPABI: data processing & analysis for (resting-state) brain imaging. Neuroinformatics 14, 339–351. doi: 10.1007/s12021-016-9299-9294
Zang, Y., Jiang, T., Lu, Y., He, Y., and Tian, L. (2004). Regional homogeneity approach to fMRI data analysis. Neuroimage 22, 394–400. doi: 10.1016/j.neuroimage.2003.12.030
Zhang, H., Yin, X., Ouyang, Z., Chen, J., Zhou, S., Zhang, C., et al. (2016). A prospective study of freezing of gait with early Parkinson disease in chinese patients. Medicine 95:e4056. doi: 10.1097/MD.0000000000004056
Keywords: resting-state fMRI, Parkinson’s disease, freezing of gait, regional homogeneity, movement function
Citation: Liu Y, Li M, Chen H, Wei X, Hu G, Yu S, Ruan X, Zhou J, Pan X, Li Z, Luo Z and Xie Y (2019) Alterations of Regional Homogeneity in Parkinson’s Disease Patients With Freezing of Gait: A Resting-State fMRI Study. Front. Aging Neurosci. 11:276. doi: 10.3389/fnagi.2019.00276
Received: 27 February 2019; Accepted: 25 September 2019;
Published: 15 October 2019.
Edited by:
Xun Chen, University of Science and Technology of China, ChinaReviewed by:
Yu-Feng Zang, Hangzhou Normal University, ChinaXiaodong Zhu, Tianjin Medical University General Hospital, China
Martin J. McKeown, The University of British Columbia, Canada
Copyright © 2019 Liu, Li, Chen, Wei, Hu, Yu, Ruan, Zhou, Pan, Li, Luo and Xie. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mengyan Li, ZG9jbG15QDE2My5jb20=; Xinhua Wei, d2VpeGluaHVhQGFsaXl1bi5jb20=; Yaoqin Xie, eXEueGllQHNpYXQuYWMuY24=