Skip to main content

ORIGINAL RESEARCH article

Front. Aging Neurosci., 29 March 2022
Sec. Parkinson’s Disease and Aging-related Movement Disorders
This article is part of the Research Topic Translational Advances in Alzheimer's, Parkinson's, and other Dementia: Molecular Mechanisms, Biomarkers, Diagnosis, and Therapies, Volume II View all 58 articles

Aberrant Volume-Wise and Voxel-Wise Concordance Among Dynamic Intrinsic Brain Activity Indices in Parkinson’s Disease: A Resting-State fMRI Study

\r\nYuan Tian,Yuan Tian1,2Hai-Bo ChenHai-Bo Chen1Xin-Xin MaXin-Xin Ma1Shu-Hua LiShu-Hua Li1Chun-Mei LiChun-Mei Li3Shao-Hui Wu,Shao-Hui Wu1,2Feng-Zhi Liu,Feng-Zhi Liu1,2Yu Du,Yu Du1,2Kai Li*Kai Li1*Wen Su,*Wen Su1,2*
  • 1Department of Neurology, Beijing Hospital, National Center of Gerontology, Institute of Geriatric Medicine, Chinese Academy of Medical Sciences, Beijing, China
  • 2Graduate School of Peking Union Medical College, Dongcheng, Beijing, China
  • 3Department of Radiology, Beijing Hospital, National Center of Gerontology, Institute of Geriatric Medicine, Chinese Academy of Medical Sciences, Beijing, China

Researches using resting-state functional magnetic resonance imaging (rs-fMRI) have applied different regional measurements to study the intrinsic brain activity (IBA) of patients with Parkinson’s disease (PD). Most previous studies have only examined the static characteristics of IBA in patients with PD, neglecting the dynamic features. We sought to explore the concordance between the dynamics of different rs-fMRI regional indices. This study included 31 healthy controls (HCs) and 57 PD patients to calculate the volume-wise (across voxels) and voxel-wise (across periods) concordance using a sliding time window approach. This allowed us to compare the concordance of dynamic alterations in frequently used metrics such as degree centrality (DC), global signal connectivity (GSC), voxel-mirrored heterotopic connectivity (VMHC), the amplitude of low-frequency fluctuations (ALFF), and regional homogeneity (ReHo). We analyzed the changes of concordance indices in the PD patients and investigated the relationship between aberrant concordance values and clinical/neuropsychological assessments in the PD patients. We found that, compared with the HCs, the PD patients had lower volume concordance in the whole brain and lower voxel-wise concordance in the posterior cerebellar lobe, cerebellar tonsils, superior temporal gyrus, and supplementary motor region. We also found negative correlations between these concordance alterations and patients’ age. The exploratory results contribute to a better understanding of IBA alterations and pathophysiological mechanisms in PD.

Introduction

Parkinson’s disease (PD) is characterized by the decrease of dopaminergic neurons in the nigrostriatal pathway, resulting in a neurodegenerative disorder with progressive movement disorders (Kalia and Lang, 2015). Moreover, non-motor symptoms such as emotional disorders and cognitive impairment are also common in PD patients (Rana et al., 2015; Raza et al., 2019). Many neuroimaging studies have demonstrated that several regions of the brain are also impacted in PD patients (Weingarten et al., 2015). Resting-state functional magnetic resonance imaging (rs-fMRI) has also been widely used to study the IBA in PD, including both brain functional connectivity (FC; Tessitore et al., 2019) and regional cerebral neuronal activity (Yue et al., 2020), which has produced inconsistent results. Some meta-analyses of neuroimaging data have attempted to identify abnormalities in the consistency of IBA patterns in PD, such as studying the amplitude of low-frequency fluctuations (ALFF; Pan et al., 2017; Wang et al., 2018) or regional homogeneity (ReHo; Pan et al., 2016), and functional connectivity network (Tahmasian et al., 2017). There are also relevant studies on the degree centrality (DC; Guo et al., 2020; Chen et al., 2021), and voxel-minored homotopic connectivity (VMHC; Luo et al., 2015). However, most current rs-fMRI research in PD patients has focused on static characteristics of human brain activity and neglected the dynamic properties of IBA in the temporal dimension (Hutchison et al., 2013). In recent years, the sliding window approach has demonstrated the altered dynamic functional connectivity (dFC) of IBA in individuals with common neuropsychiatric disorders such as Alzheimer’s disease (de Vos et al., 2018), major depressive disorder (Hosenfeld et al., 2015), and primary generalized epilepsy (Liu et al., 2017).

Many studies have assessed dFC in PD and demonstrated that dynamical changes in IBA were associated with motor (Kim et al., 2017) and cognitive impairments (Fiorenzato et al., 2019). However, it is not enough to focus only on the dynamic functional connectivity features, because much information is found in the fluctuations of regional neural activity (Yu et al., 2019). Different rs-fMRI indicators represent different IBA aspects. The degree of integration is represented by the concordance between the different dynamic regional indices. Yan et al. investigated the volume-wise and voxel-wise concordance between these regional measures and found they demonstrated stable inter-individual differences that were correlated with age (Yan et al., 2017). Fu et al. found consistency in the changes in dynamic ALFF and dynamic FC over time, and this consistency differed between patients with schizophrenia and HCs (Fu et al., 2018). Whether the consistency between regional indicators is abnormal in PD patients remains a direction to be explored. Previous studies on dynamic indicators in PD have been performed with the dynamical amplitude of low-frequency fluctuation (dALFF) analysis (Zhang et al., 2019). Though the PD patients exhibit abnormalities in dynamic and static brain activity, the concordance between different dynamic regional indices remains unclear. In this study, we compared the differences in volume-wise and voxel-wise concordance among IBA measures between PD patients and HCs to explore the pathophysiological mechanisms of PD.

Materials and Methods

Participants

We recruited 67 individuals with PD and 34 healthy controls. Both groups were similar in age and gender and had no history of neurological or mental disorders. In this study, all the PD patients were diagnosed based on the diagnostic criteria established by the United Kingdom Parkinson’s Disease Society [17]. The Helsinki Declaration was closely followed at every step of this study. Before participating in the trial, all participants were required to provide written informed consent. This research was conducted according to all relevant regulations and laws and the study was authorized by the hospital’s ethical review committee. All of the individuals were right-handed. We have ruled out conditions such as deep brain stimulation, head trauma, moderate to severe tremors in the head, alcohol/drug abuse, dementia, and other mental and neurological issues. When all PD medications were stopped for 12 h, all patients with Parkinson’s disease were administered an MRI and non-motor/motor function tests in an “off” state: several demographic and clinical questions were posed to each PD participant. These tests included the Mini-Mental State Examination (MMSE), the Hamilton Depression Rating Scale (HAMD), the Hamilton Anxiety Rating Scale (HAMA), the Non-Motor Symptoms Questionnaire (NMSQ), the Unified Parkinson’s Disease Rating Scale (UPDRS), and the Hoehn-Yahr staging. Additionally, the MMSE was used to examine the participants in the control group.

MRI Data Acquisition

MRI scans were performed on all patients at Beijing Hospital using a 3.0-T scanner (Achieva TX; Philips Medical Systems, Best, Netherlands). During the scan, the participants were instructed to remain motionless, relax, close their eyes, and maintain their alertness. The use of foam cushion and headphones helped to keep the scanning noise and head movements to a minimum. The following settings were used to generate high-resolution T1-weighted images (three-dimensional turbo field echo): field of view (FOV) = 240 × 240 mm, echo time (TE) = 3.0 ms, repetition time (TR) = 7.4 ms, flip angle = 8°, slice thickness = 1.2 mm, matrix size = 256 × 256, voxel dimensions = 0.94 × 0.94 × 1.20 mm, 140 slices. The functional pictures were generated using the following parameters: FOV = 240 × 240 mm, flip angle = 90, TR = 3,000 ms, TE = 35 ms, voxel dimensions 3.75 × 3.75 × 4.00 mm, matrix size = 64 × 64, slices = 33, time points = 210, slice thickness = 4 mm.

Resting-State Functional Magnetic Resonance Imaging Data Preprocessing

Data preprocessing for the rs-fMRI was performed using RESTPlus V 1.2 (Jia et al., 2019) and the SPM12 program1 in MATLAB (MathWorks, Inc., Natick, MA, United States). Ten volumes were initially eliminated to stabilize the magnetization and adjust it to the subject’s magnetic fields. Participants with head motion above 2 mm in displacement or 2° in rotation were then excluded from the study (Supplementary Table 1). The remaining 200 volumes were adjusted according to the slice time and readjusted to account for the head motion. Individual T1WI pictures were co-registered to the mean functional images, and then segmented into the cerebrospinal fluid (CSF), white matter, and gray matter. The functional volume of each individual was segmented and then normalized to the Montreal Neurological Institute (MNI) space using the Diffeomorphic Anatomical Registration Through Exponentiated Lie Algebra (DARTEL) toolbox for two-dimensional anatomical registration. They were then resampled to 3 mm isotropic voxels (Ashburner, 2007), which were used to determine the functional volumes of each individual. The rs-fMRI indices were extracted after functional volumes were smoothed with a 6-mm full width at half maximum (FWHM) Gaussian kernel, and the linear trend of the time course was eliminated. The volumes were then filtered with a bandpass filter between 0.01 and 0.08 Hz. The nuisance covariate regression included the Friston 24-head motion parameters, white matter and CSF signals (Power et al., 2012).

Calculation of Dynamic Resting-State Functional Magnetic Resonance Imaging Measures

The dynamic rs-fMRI measures and their concordance were computed using DPABI V6.0 (Yan et al., 2016). Sliding time window analysis was used to examine the temporal dynamics of the five rs-fMRI measurements listed below. The preprocessed functional data from each participant was subjected to hamming windows (window size = 30 TR; window step = 1 TR) in order to produce a series of BOLD signal windows. For each window, the dynamics of these rs-fMRI measurements were then calculated (Supplementary Figures 1, 2).

(1) Amplitude of low-frequency fluctuation (ALFF): We used a Fast Fourier Transform to create the power spectrum by transforming the BOLD time course to a frequency domain (FFT). In this study, ALFF was defined as the mean power spectrum in a particular low-frequency band (0.01–0.1 Hz) (Zang et al., 2007), while fALFF was defined as the ratio of the power spectrum in the low-frequency band (0.01–0.1 Hz) to the entire power spectrum in the frequency range (Zou et al., 2008). Because ALFF is sensitive when evaluating aberrant brain activity in individuals with neuropsychiatric illnesses, we only used ALFF in the concordance analysis (de Vos et al., 2018).

(2) Regional homogeneity (ReHo): Regional brain activity was quantified using the ReHo method, which presumes that neural activity is more likely to occur in a cluster than in a single voxel, and is thus more consistent. Kendall’s coordination coefficient (Kendall’s W) was obtained when the time series of a particular voxel was compared with the time series of its nearest neighbors (Zang et al., 2004).

(3) Voxel-mirrored homotopic connectivity (VMHC): To enhance the relationship between symmetric voxels, we converted the functional pictures to a symmetric space from a functional image space. We then averaged all of the normalized T1 pictures to obtain a mean T1 image. This was then averaged with its left and right mirror counterparts to create a group-specific symmetrical template using the mean T1 image as a guide. The non-linear registration of each subject was then adjusted to the standard space using the symmetric template and the functional data of each subject was changed to the symmetric space using the symmetric template. The Pearson’s correlation coefficient between the time series of the symmetric inter-hemispheric counterpart and the time series of each voxel was defined as the Pearson’s correlation coefficient between any pair of symmetric inter-hemispheric voxels (Zuo et al., 2010).

(4) Degree centrality (DC): DC is defined as the number of edges connecting a node (binary graphs) or the total of their weights (weighted graphs). The Pearson correlation coefficients were between the BOLD time courses of all pairs of voxels, after which a gray matter functional connectivity matrix was constructed for each participant. The DC was calculated as the total number of positive functional connections larger than the threshold of 0.25 between a specific voxel and all other voxels in the gray matter surrounding that voxel (Zuo et al., 2012).

(5) Global signal connectivity (GSC): The gray matter correlation coefficient (GSC) was determined as the Pearson’s correlation coefficient between the mean time course and the BOLD time course of each voxel throughout the whole gray matter. The Fisher Z transformation was then applied to the resulting GSC maps for further analysis (Hahamy et al., 2014).

Volume-Wise and Voxel-Wise Concordance

Since Kendall’s W is a non-parametric statistic that does not require distributional assumptions and is not impacted by scale discrepancies between these rs-fMRI indices, it was used to quantify the volume and voxel consistency among them (Yan et al., 2017). We calculated two kinds of consistency among the five abovementioned dynamic indices: (1) the volume-wise consistency index is an index reflecting global-level consistency and used Kendall’s W to calculate rs-fMRI indices in each subject’s brain in the scan and across all time windows; (2) the voxel-wise concordance index used Kendall’s W to calculate each voxel coordination across time windows and provided an interactive voxel concordance map (Yan et al., 2017).

Statistical Analysis

A two-sample t-test was used to assess the demographic information and neuropsychological scale scores of each group. The quantitative data were showed as a mean ± standard deviation. The chi-square test was used to compare participant gender characteristics across groups, while the inter-group volume-wise concordance was also analyzed using a two-sample t-test. Clinical data were analyzed using SPSS (SPSS Version 21.0. Inc., Chicago, IL, United States). P < 0.05 was considered significant. The normalized voxel-wise concordance maps of each rs-fMRI index were then compared between the two groups using ANCOVA analysis using DPABI V6.0. Age and gray matter volume (GMV) were included as covariates and subsequently adjusted by regression throughout the statistical analysis to reduce the impact of confounding factors. The voxel-level and cluster-level P < 0.05 significance levels were used to correct multiple comparisons using a Gaussian random field correction (GRF) with a voxel-level P < 0.001 and cluster-level P < 0.05 (Chen et al., 2018). The concordance indices were extracted from areas with significant between group differences and the correlations between the extracted concordance indices values and HAMA, HAMD, NMSQ, MMSE, UPSR scores, H-Y staging, and disease duration were assessed using Spearman correlation analyses.

Results

Demographic and Clinical Characteristics

This research assessed a total of 57 individuals with PD and 31 HCs. Four PD patients and two controls were excluded from the study because they were left-handed. Five PD patients and one healthy control were omitted from the study due to excessive head motion. One PD disease patient was excluded from the study due to low picture quality. Table 1 summarizes the demographic and clinical data. No statistically significant differences were found in the MMSE scores, age, and gender between the PD and the HCs (p > 0.05). Table 1 summarizes demographic and clinical data. There were no differences in the age, gender, or MMSE scores between PD patients and HCs (p > 0.05).

TABLE 1
www.frontiersin.org

Table 1. Clinical and demographic characteristics of the Parkinson’s disease (PD) patients and the healthy controls (HCs).

Volume-Wise Concordance Alterations in Parkinson’s Disease

We observed a significant degree of volume-wise concordance across the whole brain in the HCs and PD patients (mean Kendall’s W ≥ =0.4), indicating a high degree of spatial distribution consistency across these rs-fMRI measurements. The mean value of volume consistency significantly differed between the two groups (P = 0.007), indicating a statistically significant difference between the two groups. However, there was no statistically significant difference in the standard deviation (SD) of volume-wise concordance (P = 0.299) (Figure 1 and Table 2). We also identified statistically significant negative associations between volume-wise concordance and age/UPDRS scores in individuals with PD (r = −0.439, −0.309; P = 0.001, 0.019) (Figure 2).

FIGURE 1
www.frontiersin.org

Figure 1. Comparison of volume-wise concordance indices between PD patients and HCs. (A) Volume-wise concordance time series of the mean value in PD patients and HCs. (B) Comparison of mean volume-wise concordance indices. (C) Comparison of the SD of volume-wise concordance indices. PD, Parkinson’s disease; HCs, healthy controls. SD, Standard deviation. **P < 0.01; ns, no significance.

TABLE 2
www.frontiersin.org

Table 2. Comparison of volume-wise concordance indices between PD patients and HCs.

FIGURE 2
www.frontiersin.org

Figure 2. Negative associations between volume-wise concordance and age (A) UPDRS (B) scores in individuals with PD (r = –0.439, –0.309; P = 0.001, 0.019).

Voxel-Wise Concordance Alterations in Parkinson’s Disease

As displayed in Figure 3, inter-group comparisons indicated that PD patients exhibited worse voxel-wise concordance than the HCs in the following regions: the cerebellum posterior lobe (CPL), cerebellar tonsil (CTO), superior temporal gyrus (STG), and supplementary motor area (SMA) (GRF, voxel level, p < 0.001; cluster level, p < 0.05) (Figure 3 and Table 3). Additionally, the Spearman correlation analyses demonstrated a connection between age/MMSA scores and voxel-wise concordance values in the left CPL between PD patients and HCs (r = −0.282, 0.331; P = 0.033, 0.012) (Figure 4A). The concordance value in the CTO is negatively correlated with UPDRS scores (r = −0.269; P = 0.043) (Figure 4B), and there is a negative association between age, UPDRS scores, and the concordance values in STG (r = −0.278, −0.304; P = 0.036, 0.021) (Figures 4C–E).

FIGURE 3
www.frontiersin.org

Figure 3. Brain regions of voxel-wise concordance differences between PD patients and HCs: the cerebellum posterior lobe (CPL) and cerebellar tonsil (CTO) (A), superior temporal gyrus (STG) (B), and supplementary motor area (SMA) (C).

TABLE 3
www.frontiersin.org

Table 3. Brain regions of voxel-wise concordance differences between PD patients and HCs.

FIGURE 4
www.frontiersin.org

Figure 4. Connection between age/MMSA scores and voxel-wise concordance value in the left CPL in PD patients from HCs (r = –0.282,0.331; P = 0.033, 0.012). (A,B) Negative correlation between concordance value in the CTO and UPDRS scores (r = –0.269; P = 0.043). (C) Negative association between age, UPDRS scores, and the concordance value in STG (r = –0.278, –0.304; P = 0.036, 0.021) (D,E).

Discussion

Concordance abnormalities in rs-fMRI brain activity data were significantly altered in PD patients based on a dynamic analytic technique. PD patients had worse voxel concordance in certain cerebral regions and lower volume concordances throughout the whole brain than did the healthy controls. Furthermore, age was related to changes in volume concordance, and motor symptoms and age were associated with decreases in voxel concordance regions. This demonstrates that coordination in Parkinson’s patients is clinically meaningful.

Previous research has applied various static and dynamic indicators to explore the characteristics of spontaneous brain activity in PD patients, which provided helpful information about certain brain regions (Li et al., 2018; Wang et al., 2018; Filippi et al., 2019). However, certain regions demonstrated various degrees of dysfunction, and some measures were the only way of detecting functional impairments in certain brain regions. The distinctions and overlaps between these rs-fMRI measures remain unknown. Resting-state indicators are constantly being developed, leading researchers to examine the consistency between multiple indices. Therefore, based on the dynamic nature of IBA, we further investigated the concordance of different rs- fMRI indices. Different indices reflect different aspects of IBA, which means that their concordance can be measured by reflecting the integration of different functional levels (Yan et al., 2017). The volume-wise concordance demonstrates the consistency of the spatial distribution of these rs-fMRI values. As previously reported (Yan et al., 2017), a high degree of spatial distribution consistency was observed with high volume-wise concordance in HCs (mean Kendall’s W > 0.4). We also observed volume-wise concordance in PD patients (mean Kendall’s W≈0.4), which was significantly lower than HCs. This indicates a less consistent spatial distribution, which could be due to the inconsistent pattern of spatial variation of rs-fMRI measurements in PD. Moreover, the inter-voxel agreement reflects the concordance of temporal dynamics between these rs-fMRI measures. In this study, the frontal lobe, temporal gyrus, and cerebellum of PD patients displayed decreased voxel-wise concordance. The reason for voxel-wise concordance disruptions could be the inconsistent pattern of temporal dynamic changes in rs-fMRI measurements. For detecting regional functional abnormalities, voxel-wise concordance may be the novel symbols based on the single rs-fMRI measurements. For example, concordance studies have been performed on schizophrenia (Zhu et al., 2018), subjective cognitive decline (SCD; Yang et al., 2021), and Alzheimer’s disease (AD; Chen et al., 2021), which have identified abnormalities in spontaneous dynamic brain activity and integration.

The SMA is located on the medial side of the cerebral hemisphere, in front of the primary motor cortex. It is closely related to motor planning, learning, and cognitive functions (Nachev et al., 2008). The SMA is connected to the basal ganglia through the thalamus, forming the striatal-thalamocortical loop (STC loop). On the other hand, the SMA connects to the cerebellum through the thalamus to form the cerebellum-thalamocortical loop (CTC loop). Studies have demonstrated that the primary pathogenesis of PD is damage to the STC loop, which disrupts the intracranial control circuit and causes a relatively active CTC loop to compensate for functional deficits caused by this damage (Palmer et al., 2010). PD patients have impaired SMA function, which decreases connectivity between the SMA and other functional brain areas and results in motor impairments in PD patients (Yu et al., 2007; Prodoehl et al., 2010; Wu et al., 2010).

Parkinson’s disease patients struggle in rhythm recognition assignments. The right STG plays a crucial role in this process and increased activity in the right STG region suggests inhibited temporal auditory perception (Vikene et al., 2019). The medial temporal lobe is an important center of emotional integration and is associated with motor control, self-perception, and emotional regulation (Paus et al., 1996). While PD has typically been characterized by motor impairment, it can also produce several different non-motor symptoms (Reichmann et al., 2016). We found a significant correlation between the UPDRS score and the voxel-wise concordance value of STG. This could represent enhanced temporal lobe activity as a compensatory mechanism for impaired motor symptoms.

The cerebellum coordinates somatic balance, muscle tone, and casual movements, and is a key node in the CTC loop, controlling and influencing somatic movements (Yang et al., 2013; Hanakawa et al., 2017). Prior studies have found that neuronal activity in cerebellar hemispheres increases in patients with PD compared with HCs (Yu et al., 2007), and that deep stimulation of the subthalamic nuclei inhibits glucose metabolism in the cerebellum (Hilker et al., 2004). Therefore, increased neuronal activity in the cerebellum of PD patients could be functional compensation for the hypofunction of the basal ganglia (Lewis et al., 2007). Several studies have demonstrated that Parkinson’s tremor is closely related to the CTC loop (Hacker et al., 2012; Helmich et al., 2012; Zhang et al., 2015). Additionally, pathological and compensatory effects are believed to result from the cerebellum in PD patients. Certain clinical symptoms could be caused by pathological cerebellar alterations resulting from degeneration of aberrant drives from the basal ganglia, dopaminergic system, and dopaminergic therapy. This compensatory effect could help maintain both functional motor and non-motor systems (Wu and Hallett, 2013). Additionally, recruitment of the CTC circuit is positively related with the severity of symptoms or development in PD of mild to moderate stages (Sen et al., 2010). Neuronal activity in the cerebellar hemispheres is increased in patients with PD compared to HCs (Yu et al., 2007), and its deep stimulation of the subthalamic nuclei clusters aligns glucose metabolism in the cerebellar brain (Hilker et al., 2004). Pathological damages increase in severity as the disease progresses, while the compensatory effect may weaken or fail altogether (Jankovic, 2005). This study demonstrates that concordance indexes in the posterior cerebellar lobes and cerebellar tonsil decreased in PD patients and that concordance values are correlated with age, motor, and non-motor symptoms. This indicates that neuronal activity increased, which compensates for decreases in the STC loop function.

Limitation

This study has some limitations. First, the size of this study sample was insufficient; therefore, a study with a larger number of participants is necessary. For the concordance studies, we only used five common rs-fMRI measures since these data-driven measures are less likely to be influenced by variables such as component selection in independent component analysis and seed specification in seed-based correlation analysis. Accounting for additional rs-fMRI measurements could allow for better IBA characteristics and a more nuanced understanding of the physiological mechanisms responsible for these activities.

Conclusion

In conclusion, this study demonstrates that PD patients have altered patterns of coherence obtained from rs-fMRI across several routine IBA measurements. We believe that concordance measures can provide novel insights into IBA mechanisms in PD patients.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by Beijing Hospital. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

WS, H-BC, and KL conceived and designed the experiments. YT and KL analyzed the fMRI data. C-ML was responsible for the fMRI scan and helped perform fMRI data analyses. WS, S-HL, and H-BC recruited the subjects. X-XM, F-ZL, S-HW, and YD collected the demographic, clinical, and neuropsychological information of the subjects. YT, KL, and WS wrote the manuscript. All authors read and approved the final manuscript.

Funding

This study was funded by the National Key Research and Development Program of China (No. 2020YFC2006404, 2020YFC2006402, and 2020YFC2006400), the project “National Major Multidisciplinary Cooperative Diagnosis and Treatment Capacity Building Project” from the National Health Commission of the People’s Republic of China, Beijing Hospital grant (BJ-2021-182), the National Key R&D Program of China (No. 2017YFC1310200), and the 12th Five-year Plan for National Science and Technology Supporting Program (Nos. 2012BAI10B03 and 2012BAI10B04).

Conflict of Interest

The 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.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We are grateful for the cooperation of all participants.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2022.814893/full#supplementary-material

Footnotes

  1. ^ www.fifil.ion.ucl.ac.uk/spm

References

Ashburner, J. (2007). A fast diffeomorphic image registration algorithm. Neuroimage 38, 95–113. doi: 10.1016/j.neuroimage.2007.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Lu, B., and Yan, C. G. (2018). Reproducibility of R-fMRI metrics on the impact of different strategies for multiple comparison correction and sample sizes. Hum. Brain Mapp. 39, 300–318. doi: 10.1002/hbm.23843

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Onur, O. A., Richter, N., Fassbender, R., Gramespacher, H., Befahr, Q., et al. (2021). Concordance of intrinsic brain connectivity measures is disrupted in Alzheimer’s Disease. Brain Connect. [Epub ahead of print]. doi: 10.1089/brain.2020.0918

PubMed Abstract | CrossRef Full Text | Google Scholar

de Vos, F., Koini, M., Schouten, T. M., Seiler, S., van der Grond, J., Lechner, A., et al. (2018). A comprehensive analysis of resting state fMRI measures to classify individual patients with Alzheimer’s disease. Neuroimage 167, 62–72. doi: 10.1016/j.neuroimage.2017.11.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Filippi, M., Spinelli, E. G., Cividini, C., and Agosta, F. (2019). Resting state dynamic functional connectivity in neurodegenerative conditions : a review of magnetic resonance imaging findings. Front. Neurosci. 13:657. doi: 10.3389/fnins.2019.00657

PubMed Abstract | CrossRef Full Text | Google Scholar

Fiorenzato, E., Strafella, A. P., Kim, J., Schifano, R., Weis, L., Antonini, A., et al. (2019). Dynamic functional connectivity changes associated with dementia in Parkinson’s disease. Brain 142, 2860–2872. doi: 10.1093/brain/awz192

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, Z., Tu, Y., Di, X., Du, Y., Pearlson, G. D., Turner, J. A., et al. (2018). Characterizing dynamic amplitude of low-frequency fluctuation and its relationship with dynamic functional connectivity: an application to schizophrenia. Neuroimage 180, 619–631. doi: 10.1016/j.neuroimage.2017.09.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, M. R., Ren, Y., Yu, H. M., Yang, H. G., Cao, C. H., Li, Y. M., et al. (2020). Alterations in degree centrality and functional connectivity in Parkinson’s disease patients with freezing of gait: a resting-state functional magnetic resonance imaging study. Front. Neurosci. 14:582079. doi: 10.3389/fnins.2020.582079

PubMed Abstract | CrossRef Full Text | Google Scholar

Hacker, C. D., Perlmutter, J. S., Criswell, S. R., Ances, B. M., and Snyder, A. Z. (2012). Resting state functional connectivity of the striatum in Parkinson’s disease. Brain 135, 3699–3711. doi: 10.1093/brain/aws281

PubMed Abstract | CrossRef Full Text | Google Scholar

Hahamy, A., Calhoun, V., Pearlson, G., Harel, M., Stern, N., Attar, F., et al. (2014). Save the global: global signal connectivity as a tool for studying clinical populations with functional magnetic resonance imaging. Brain Connect. 4, 395–403. doi: 10.1089/brain.2014.0244

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanakawa, T., Goldfine, A. M., and Hallett, M. (2017). A common function of basal ganglia-cortical circuits subserving speed in both motor and cognitive domains. eNeuro 4:ENEURO.0200-17.2017. doi: 10.1523/ENEURO.0200-17.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Helmich, R. C., Hallett, M., Deuschl, G., Toni, I., and Bloem, B. R. (2012). Cerebral causes and consequences of parkinsonian resting tremor: A tale of two circuits? Brain 135, 3206–3226. doi: 10.1093/brain/aws023

PubMed Abstract | CrossRef Full Text | Google Scholar

Hilker, R., Voges, J., Weisenbach, S., Kalbe, E., Burghaus, L., Ghaemi, M., et al. (2004). Subthalamic nucleus stimulation restores glucose metabolism in associative and limbic cortices and in cerebellum: evidence from a FDG-PET Study in Advanced Parkinson’s Disease. J. Cereb. Blood Flow Metab. 24, 7–16. doi: 10.1097/01.WCB.0000092831.44769.09

CrossRef Full Text | Google Scholar

Hosenfeld, B., Bos, E. H., Wardenaar, K. J., Conradi, H. J., van der Maas, H. L. J., Visser, I., et al. (2015). Major depressive disorder as a nonlinear dynamic system: bimodality in the frequency distribution of depressive symptoms over time. BMC Psychiatry 15:222. doi: 10.1186/s12888-015-0596-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Hutchison, R. M., Womelsdorf, T., Allen, E. A., Bandettini, P. A., Calhoun, V. D., Corbetta, M., et al. (2013). Dynamic functional connectivity: promise, issues, and interpretations. Neuroimage 80, 360–378. doi: 10.1016/j.neuroimage.2013.05.079

PubMed Abstract | CrossRef Full Text | Google Scholar

Jankovic, J. (2005). Progression of Parkinson Disease. Arch. Neurol. 62, 351–352. doi: 10.1001/archneur.62.3.351

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, X., Li, Y., Li, K., Liang, P., and Fu, X. (2019). Precuneus dysfunction in Parkinson’s disease with mild cognitive impairment. Front. Aging Neurosci. 11:427. doi: 10.3389/fnagi.2018.00427

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalia, L. V., and Lang, A. E. (2015). Parkinson’s disease. Lancet 386, 896–912. doi: 10.1016/S0140-6736(14)61393-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J., Criaud, M., Cho, S. S., Díez-Cirarda, M., Mihaescu, A., Coakeley, S., et al. (2017). Abnormal intrinsic brain functional network dynamics in Parkinson’s disease. Brain 140, 2955–2967. doi: 10.1093/brain/awx233

PubMed Abstract | CrossRef Full Text | Google Scholar

Lewis, M. M., Slagle, C. G., Smith, A. B., Truong, Y., Bai, P., McKeown, M. J., et al. (2007). Task specific influences of Parkinson’s disease on the striato-thalamo-cortical and cerebello-thalamo-cortical motor circuitries. Neuroscience 147, 224–235. doi: 10.1016/j.neuroscience.2007.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, K., Su, W., Li, S. H., Jin, Y., and Chen, H. B. (2018). Resting State fMRI: a valuable tool for studying cognitive dysfunction in PD. Parkinsons Dis. 2018:6278649. doi: 10.1155/2018/6278649

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, F., Wang, Y., Li, M., Wang, W., Li, R., Zhang, Z., et al. (2017). Dynamic functional network connectivity in idiopathic generalized epilepsy with generalized tonic – Clonic Seizure. Hum. Brain Mapp. 38, 957–973. doi: 10.1002/hbm.23430

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, C., Guo, X., Song, W., Zhao, B., Cao, B., Yang, J., et al. (2015). Decreased resting-state interhemispheric functional connectivity in Parkinson’s Disease. Biomed Res. Int. 2015:692684. doi: 10.1155/2015/692684

PubMed Abstract | CrossRef Full Text | Google Scholar

Nachev, P., Kennard, C., and Husain, M. (2008). Functional role of the supplementary and pre-supplementary motor areas. Nat. Rev. Neurosci. 9, 856–869. doi: 10.1038/nrn2478

PubMed Abstract | CrossRef Full Text | Google Scholar

Palmer, S. J., Li, J., Wang, Z. J., and McKeown, M. J. (2010). Joint amplitude and connectivity compensatory mechanisms in Parkinson’s disease. Neuroscience 166, 1110–1118. doi: 10.1016/j.neuroscience.2010.01.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, P., Zhan, H., Xia, M., Zhang, Y., Guan, D., and Xu, Y. (2016). 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, P., Zhang, Y., Liu, Y., Zhang, H., Guan, D., and Xu, Y. (2017). Abnormalities of regional brain function in Parkinson’s disease: a meta-analysis of resting state functional magnetic resonance imaging studies. Sci. Rep. 7:40469. doi: 10.1038/srep40469

PubMed Abstract | CrossRef Full Text | Google Scholar

Paus, T., Tomaiuolo, F., Otaky, N., Petrides, M., Atlas, J., Morris, R., et al. (1996). Human cingulate and paracingulate sulci: pattern, variability, asymmetryand probabilistic map. Cereb. Cortex 6, 207–214. doi: 10.1093/cercor/6.2.207

PubMed Abstract | CrossRef Full Text | Google Scholar

Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L., and Petersen, S. E. (2012). NeuroImage Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59, 2142–2154. doi: 10.1016/j.neuroimage.2011.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Prodoehl, J., Spraker, M., Corcos, D., Comella, C., and Vaillancourt, D. (2010). Blood oxygenation level-dependent activation in basal ganglia nuclei relates to specific symptoms in de novo Parkinson’s disease. Mov. Disord. 25, 2035–2043. doi: 10.1002/mds.23360

PubMed Abstract | CrossRef Full Text | Google Scholar

Rana, A. Q., Ahmed, U. S., Chaudry, Z. M., and Vasan, S. (2015). Parkinson’s disease: a review of non-motor symptoms. Expert Rev. Neurother. 15, 549–562. doi: 10.1586/14737175.2015.1038244

PubMed Abstract | CrossRef Full Text | Google Scholar

Raza, C., Anjum, R., and Shakeel, N. U. A. (2019). Parkinson’s disease: mechanisms, translational models and management strategies. Life Sci. 226, 77–90. doi: 10.1016/j.lfs.2019.03.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Reichmann, H., Brandt, M. D., and Klingelhoefer, L. (2016). The nonmotor features of Parkinson’s disease: pathophysiology and management advances. Curr. Opin. Neurol. 29, 467–473. doi: 10.1097/WCO.0000000000000348

PubMed Abstract | CrossRef Full Text | Google Scholar

Sen, S., Kawaguchi, A., Truong, Y., Lewis, M. M., and Huang, X. (2010). Dynamic changes in cerebello-thalamo-cortical motor circuitry during progression of Parkinson’s disease. Neuroscience 166, 712–719. doi: 10.1016/j.neuroscience.2009.12.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Tahmasian, M., Eickhoff, S. B., Giehl, K., Schwartz, F., Herz, D. M., Drzezga, A., et al. (2017). Resting-state functional reorganization in Parkinson’s disease: an activation likelihood estimation meta-analysis. Cortex 92, 119–138. doi: 10.1016/j.cortex.2017.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Tessitore, A., Cirillo, M., and De Micco, R. (2019). Functional connectivity signatures of Parkinson’s disease. J. Parkinsons Dis. 9, 637–652. doi: 10.3233/JPD-191592

PubMed Abstract | CrossRef Full Text | Google Scholar

Vikene, K., Skeie, G. O., and Specht, K. (2019). Compensatory task-specific hypersensitivity in bilateral planum temporale and right superior temporal gyrus during auditory rhythm and omission processing in Parkinson’s disease. Sci. Rep. 9:12623. doi: 10.1038/s41598-019-48791-0

PubMed Abstract | CrossRef Full Text | Google Scholar

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:giy071. doi: 10.1093/gigascience/giy071

PubMed Abstract | CrossRef Full Text | Google Scholar

Weingarten, C. P., Sundman, M. H., Hickey, P., and Chen, N. K. (2015). Neuroimaging of Parkinson’s disease: expanding views. Neurosci. Biobehav. Rev. 59, 16–52. doi: 10.1016/j.neubiorev.2015.09.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, T., and Hallett, M. (2013). The cerebellum in Parkinson’s disease. Brain 136, 696–709. doi: 10.1093/brain/aws360

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, T., Wang, L., Hallett, M., Li, K., and Chan, P. (2010). Neural correlates of bimanual anti-phase and in-phase movements in Parkinson’s disease. Brain 133, 2394–2409. doi: 10.1093/brain/awq151

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, C., Wang, X., Zuo, X., and Zang, Y. (2016). DPABI : data processing & analysis for (Resting-State) brain imaging. Neuroinformatics 14, 339–351. doi: 10.1007/s12021-016-9299-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, C. G., Yang, Z., Colcombe, S. J., Zuo, X. N., and Milham, M. P. (2017). Concordance among indices of intrinsic brain function: Insights from inter-individual variation and temporal dynamics. Sci. Bull. 62, 1572–1584. doi: 10.1016/j.scib.2017.09.015

CrossRef Full Text | Google Scholar

Yang, H., Zhou, X. J., Zhang, M.-M., Zheng, X. N., Zhao, Y. L., and Wang, J. (2013). Changes in spontaneous brain activity in early Parkinson’s disease. Neurosci. Lett. 549, 24–28. doi: 10.1016/j.neulet.2013.05.080

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Zha, X., Zhang, X., Ke, J., Hu, S., Wang, X., et al. (2021). Dynamics and concordance abnormalities among indices of intrinsic brain activity in individuals with subjective cognitive decline: a temporal dynamics resting-state functional magnetic resonance imaging analysis. Front. Aging Neurosci. 12:584863. doi: 10.3389/fnagi.2020.584863

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, H., Sternad, D., Corcos, D. M., and Vaillancourt, D. E. (2007). Role of hyperactive cerebellum and motor cortex in Parkinson’s disease. Neuroimage 35, 222–233. doi: 10.1016/j.neuroimage.2006.11.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, Y., Li, Z., Lin, Y., Yu, J., Peng, G., Zhang, K., et al. (2019). Depression affects intrinsic brain activity in patients with mild cognitive impairment. Front. Neurosci. 13:1333. doi: 10.3389/fnins.2019.01333

PubMed Abstract | CrossRef Full Text | Google Scholar

Yue, Y., Jiang, Y., Shen, T., Pu, J., Lai, H. Y., and Zhang, B. (2020). ALFF and ReHo mapping reveals different functional patterns in early- and Late-Onset Parkinson’s Disease. Front. Neurosci. 14:141. doi: 10.3389/fnins.2020.00141

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zang, Y. F., Yong, H., Chao-Zhe, Z., Qing-Jiu, C., Man-Qiu, S., Meng, L., et al. (2007). Altered baseline brain activity in children with ADHD revealed by resting-state functional MRI. Brain Dev. 29, 83–91. doi: 10.1016/j.braindev.2006.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Dou, B., Wang, J., Xu, K., Zhang, H., Sami, M. U., et al. (2019). Dynamic alterations of spontaneous neural activity in Parkinson’s disease: a resting-state fMRI Study. Front. Neurol. 10:1052. doi: 10.3389/fneur.2019.01052

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Wei, L., Hu, X., Xie, B., Zhang, Y., Wu, G. R., et al. (2015). Akinetic-rigid and tremor-dominant Parkinson’s disease patients show different patterns of intrinsic brain activity. Parkinsonism Relat. Disord. 21, 23–30. doi: 10.1016/j.parkreldis.2014.10.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, J., Zhu, D.-M., Qian, Y., Li, X., and Yu, Y. (2018). Altered spatial and temporal concordance among intrinsic brain activity measures in schizophrenia. J. Psychiatr. Res. 106, 91–98. doi: 10.1016/j.jpsychires.2018.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, Q. H., Zhu, C. Z., Yang, Y., Zuo, X. N., Long, X. Y., Cao, Q. J., et al. (2008). An improved approach to detection of amplitude of low-frequency fluctuation (ALFF) for resting-state fMRI: fractional ALFF. J. Neurosci. Methods 172, 137–141. doi: 10.1016/j.jneumeth.2008.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Zuo, X. N., Ehmke, R., Mennes, M., Imperati, D., Castellanos, F. X., Sporns, O., et al. (2012). Network centrality in the human functional connectome. Cereb. Cortex 22, 1862–1875. doi: 10.1093/cercor/bhr269

PubMed Abstract | CrossRef Full Text | Google Scholar

Zuo, X. N., Kelly, C., Di Martino, A., Mennes, M., Margulies, D. S., Bangaru, S., et al. (2010). Growing together and growing apart: regional and sex differences in the lifespan developmental trajectories of functional homotopy. J. Neurosci. 30, 15034–15043. doi: 10.1523/JNEUROSCI.2612-10.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Parkinson’s disease, resting-state fMRI, intrinsic brain activity, volume-wise concordance, voxel-wise concordance

Citation: Tian Y, Chen H-B, Ma X-X, Li S-H, Li C-M, Wu S-H, Liu F-Z, Du Y, Li K and Su W (2022) Aberrant Volume-Wise and Voxel-Wise Concordance Among Dynamic Intrinsic Brain Activity Indices in Parkinson’s Disease: A Resting-State fMRI Study. Front. Aging Neurosci. 14:814893. doi: 10.3389/fnagi.2022.814893

Received: 14 November 2021; Accepted: 07 March 2022;
Published: 29 March 2022.

Edited by:

Kuangyu Shi, University of Bern, Switzerland

Reviewed by:

David Eidelberg, Northwell Health, United States
Mahdieh Ghasemi, University of Neyshabur, Iran

Copyright © 2022 Tian, Chen, Ma, Li, Li, Wu, Liu, Du, Li and Su. 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: Kai Li, likai_bjh@hotmail.com; Wen Su, suwendy@126.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.