- Open Access
Oscillation-specific nodal alterations in early to middle stages Parkinson’s disease
Translational Neurodegeneration volume 8, Article number: 36 (2019)
Different oscillations of brain networks could carry different dimensions of brain integration. We aimed to investigate oscillation-specific nodal alterations in patients with Parkinson’s disease (PD) across early stage to middle stage by using graph theory-based analysis.
Eighty-eight PD patients including 39 PD patients in the early stage (EPD) and 49 patients in the middle stage (MPD) and 36 controls were recruited in the present study. Graph theory-based network analyses from three oscillation frequencies (slow-5: 0.01–0.027 Hz; slow-4: 0.027–0.073 Hz; slow-3: 0.073–0.198 Hz) were analyzed. Nodal metrics (e.g. nodal degree centrality, betweenness centrality and nodal efficiency) were calculated.
Our results showed that (1) a divergent effect of oscillation frequencies on nodal metrics, especially on nodal degree centrality and nodal efficiency, that the anteroventral neocortex and subcortex had high nodal metrics within low oscillation frequencies while the posterolateral neocortex had high values within the relative high oscillation frequency was observed, which visually showed that network was perturbed in PD; (2) PD patients in early stage relatively preserved nodal properties while MPD patients showed widespread abnormalities, which was consistently detected within all three oscillation frequencies; (3) the involvement of basal ganglia could be specifically observed within slow-5 oscillation frequency in MPD patients; (4) logistic regression and receiver operating characteristic curve analyses demonstrated that some of those oscillation-specific nodal alterations had the ability to well discriminate PD patients from controls or MPD from EPD patients at the individual level; (5) occipital disruption within high frequency (slow-3) made a significant influence on motor impairment which was dominated by akinesia and rigidity.
Coupling various oscillations could provide potentially useful information for large-scale network and progressive oscillation-specific nodal alterations were observed in PD patients across early to middle stages.
Parkinson’s disease (PD) is a chronic and progressive movement disorder characterized by heterogenous motor symptoms including tremor, akinesia and rigidity [1, 2]. Damage to substantia nigra pars compacta resulting in the depletion of dopamine is frequently considered as an important pathological hallmark of PD [3, 4], which functionally leads to disruption of the basal ganglia, a trigger of clinical motor symptoms . Clinical magnetic resonance imaging (MRI) has detected robust abnormalities in substantia nigra [6, 7], basal ganglia [8,9,10], neocortex [11,12,13] and cerebellar cortex [14, 15] across the various parkinsonian statuses. Nevertheless, the parkinsonian network abnormalities are not fully understood, so more integrated and comprehensive approaches to identify possible pathogenesis are highly desired.
Neurons and networks are endowed with complex dynamics, including their intrinsic abilities to resonate and oscillate at multiple frequencies, and these different oscillations of brain networks could carry information about different dimensions of brain integration . Although resting functional MRI (fMRI) provides dynamical brain network information with a broad power spectrum, oscillatory coupling is usually examined within a single frequency band [8,9,10, 17,18,19,20]. As exceptions, there have been studies that decomposed fMRI oscillation into multiple distinct frequency bands. However, these studies revealed oscillation-specific functional abnormalities in PD patients by measuring local spontaneous brain function [21,22,23]. Because human brain works as a large-scale network connected with intricate edges [24, 25], local measurements could not thoroughly clarify the potential PD pathogenesis [21,22,23] neither the measurement merely focusing on the alterations of striatum-based functional connectivity . Therefore, restoring human brain to a large-scale network with multiple oscillation frequencies could make insightful contributions to a better understanding of PD.
Graph theory-based analysis models the human brain as a complex large-scale network and provides a powerful mathematical framework to characterize topological organization of the human brain network [24, 27, 28]. By employing this method, a number of studies found that the integration and segregation of the global network topology were disturbed in PD patients [7, 17, 18, 29, 30]. In the estimation of each node function in the constructed large-scale network, nodal abnormalities in the sensorimotor and temporal-occipital regions were observed in PD patients [17, 18]; in particular, one of the studies reported that the nodal centrality was associated with disease stages . However, although these studies clarified PD pathogenesis in a network perspective, the BOLD signal oscillation was still analyzed within a single frequency band, overlooking the distinctive information provided by multiple frequency bands [21, 22, 26, 31,32,33]. Therefore, the influence of different frequency bands on nodal properties in the large-scale network in PD patients and aging population remains unclear. More specifically, the potential progressive alterations of oscillation-specific nodal properties across different disease stages and their relationships with motor function deficits are still unexplored in PD. Since it was suggested that specific firing patterns of the motor network mediate the heterogenous motor symptoms in PD , the relationships between oscillation-specific nodal alterations and specific motor symptoms, i.e. akinesia/rigidity and/or tremor, are unknown.
To address these questions, we first constructed large-scale network matrixes within various oscillation frequencies. Next, graph theory-based analysis was performed to measure the nodal properties within each oscillation frequency. Furthermore, the effects of different oscillation frequencies on nodal properties were investigated. We also explored oscillation-specific progressive nodal abnormalities in PD patients across early to middle stages. Finally, clinical correlation analyses were conducted between altered nodal properties and motor severity, taking into account which component of motor symptoms (akinesia/rigidity and tremor) contributed dominantly.
All PD patients and control subjects signed informed consent forms and all procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. PD diagnosis was made by a senior neurologist according to United Kindom Parkinson’s Disease Society Brain Bank criteria . For each patient, demographic information including age, gender, education, disease history and clinical assessments including the Unified Parkinson’s Disease Rating Scale (UPDRS) score, the Mini-Mental State Examination (MMSE) score, Hoehn-Yahr stage and disease duration were obtained. In detail, the UPDRS motor tremor score (sum of items 20 and 21) and the UPDRS motor akinesia/rigidity score (sum of items 22–27 and 31) as representatives of different motor impairments were calculated [8, 35]. For normal controls, above demographic and clinical information including age, gender, disease history, UPDRS motor score and MMSE score were collected. Clinical data and fMRI data were obtained after overnight withdrawal of treatment (at least 12 h) for PD patients taking anti-parkinsonian drugs.
Subjects were excluded depending on their disease history and medication status: (1) with schizophrenia, n = 1; (2) with blindness, n = 1; (3) with metal artifact for MRI scanning, n = 8; (4) with potentially cognitive impairment according to previous suggestions weighted by Chinese education (MMSE score ≤ 17 for illiterate subjects, ≤ 20 for grade-school literate, and ≤ 23 for junior high school and higher education literate) [36, 37], n = 3; (5) with excessive head motion (greater than 2 mm in transformation and 2° in rotation), n = 6; (6) with more than 1/3 bad time points removed after scrubbing, n = 4; (7) without drug withdrawal for PD patients, n = 1.
In total, 88 PD patients with Hoehn-Yahr stage ranging from 1 to 2.5 and 36 controls were included in the study. Patients with Hoehn-Yahr stage ranging from 1 to 1.5 were defined as early stage PD (EPD) while patients with Hoehn-Yahr stage ranging from 2 to 2.5 were defined as middle stage PD (MPD). Due to the advanced PD patients were not the target cohort of current research project, the small sample size of PD patients with Hoehn-Yahr stage ranging from 3 to 5 made it impractical to be included in this study.
MRI scanning was performed on a 3.0 Tesla system (GE Medical Systems, Discovery 750) equipped with an eight-channel head coil. During MRI scanning, the head of each subject was stabilized with restraining foam pads. Earplugs were provided to reduce the noise during scanning. The fMRI scanning was performed in darkness, and the participants were explicitly instructed to relax, close their eyes and not fall asleep during the fMRI acquisition. fMRI images were acquired using a Gradient Recalled Echo/Echo Planar Imaging sequence: repetition time = 2000 ms; echo time = 30 ms; flip angle =77 degrees; field of view = 240 × 240 mm2; matrix = 64 × 64; slice thickness = 4 mm; slice gap = 0 mm; 38 interleaved slices. A total of 205 volumes were acquired from each subject. Structural T1 images were acquired using a Fast Spoiled Gradient Recalled sequence: repetition time = 7.336 ms; echo time = 3.036 ms; inversion time = 450 ms; flip angle = 11 degrees; field of view = 260 × 260 mm2; matrix = 256 × 256; slice thickness = 1.2 mm; 196 continuous sagittal slices.
fMRI data preprocessing
fMRI data preprocessing was performed using the Data Processing & Analysis for (Resting-State) Brain Imaging, DPABI (http://rfmri.org/dpabi) . The first 10 volumes of rsfMRI data were discarded due to the consideration of instability of the initial MRI signal, thus 195 time points were implemented into following procedures: slice timing, realignment, nuisance covariates (Friston’s 24 head motion parameters, white matter and cerebrospinal fluid signal) regression, spatial normalization with resampling to 3 × 3 × 3 mm3 through structure images, smoothing with a Gaussian kernel of 6 × 6 × 6 mm3 full width at half maximum, detrending and scrubbing. Head motion parameter from each subject was collected for further regression analysis. To investigate functional network oscillations with different frequencies, we further divided the full frequency range (0–0.25 Hz) into five oscillation frequencies: slow-6 (0–0.01 Hz), slow-5 (0.01–0.027 Hz), slow-4 (0.027–0.073 Hz), slow-3 (0.073–0.198 Hz) and slow-2 (0.198–0.25 Hz) [16, 32]. Because signal from slow-6 and slow-2 mainly reflect low frequency drift, white matter signals, and high-frequency physiological noises, respectively [32, 39], we did not construct functional network within these two oscillations.
Whole brain network is composed of nodes and edges between nodes [24, 27]. Nodes represent brain regions, i.e., a collection of voxels. Edges represent the interregional statistical coherences of blood oxygen level-dependent signals (functional connectivity). First of all, to define the nodes, as basal ganglia dysfunction is cardinally involved in PD [5, 8,9,10], the atlas with 112 subcortical and cortical regions of interest generated from a probabilistic atlas of Harvard-Oxford Structural Atlas that defines regions based on standard anatomical boundaries (probability threshold = 25%)  was used in the present study. Each region of interest was representing a node of the network. Then, we extracted the mean time course of each node and interregional resting-state functional connectivity was evaluated by calculating the Pearson correlation between the time courses of each node pair. Fisher’s r-to-z transformation was applied to improve data distributions for parametric statistical analysis. Finally, network matrixes (112 × 112) within three corresponding oscillations (slow-5, slow-4 and slow-3) for each subject were generated for following graph theory-based network analysis.
Graph theory-based network analysis
Before calculating network metrics, each constructed matrix was thresholded into a binarized matrix with a fixed sparsity value which defined as the total number of edges in a network divided by the maximum possible number of edges. Because the topological property computation has a strong dependency on network sparsity, a range of sparsity from 5 to 50% with an interval of 5% instead of a single threshold was selected, in which the constructed network has prominent small-world properties . Therefore, by setting a sparsity-specific threshold, the networks to be analyzed from each group had the same number of edges and the potential discrepancies in the overall functional connectivity would be minimized. Finally, to investigate group differences in these networks, we calculated the area under the curve (AUC) for each network metric (property), which provides a summarized scalar for topological characterization of brain networks independent of single threshold selection [17, 18, 41].
For the constructed functional networks, we calculated topological nodal metrics for each subject [24, 27, 28, 42]. The metrics included nodal degree centrality, nodal efficiency and nodal betweenness centrality. Briefly, degree centrality for a given node reveals its information communication ability in the functional network; nodal betweenness centrality reflects its effect on information flow between other nodes; and nodal efficiency characterizes the efficiency of parallel information transfer of that node in the network.
To detect group differences in AUC values for each nodal metric, general linear model was used with age, gender and head motion as covariates of no interests. For multiple comparisons of nodal properties, we used a false-positive correction of each one (112 nodes in total), p = 1/112 (1/N) = 0.009, where N is the number of comparisons, which implies that we expected less than one false positive per analysis on average [41, 43]. Of note, when the comparisons of nodal metrics in general linear model mentioned above were performed among three groups (EPD, MPD and controls), the p value was adjusted by using Bonferroni correction. For example, the threshold of ‘p = 0.009’ from the comparison of each nodal metric among three groups was first adjusted on multiple-group level with Bonferroni correction and then on multiple-node level (1/112).
All these network analyses were conducted using the GRETNA toolbox (http://www.nitrc.org/projects/gretna/)  and IBM SPSS 19.0 and the results were visualized using BrainNet Viewer (http://www.nitrc.org/projects/bnv/) .
Statistics analysis for demographic and clinical information
The normal distribution of data was confirmed using the one-sample Kolmogorov–Smirnov test. Differences in the age, education, disease duration, UPDRS motor score (including akinesia/rigidity score and tremor score) and Hoehn-Yahr stage distribution between whole PD group and normal controls or PD groups or among EPD, LPD and controls were compared with independent t test or analysis of variance. Gender distribution between/among groups was compared with Pearson chi-square. Due to the non-normality of data distribution, the differences of MMSE score and head motion between the two PD groups or among three groups were compared with Mann-Whitney U test or Kruskal-Wallis test appropriately. The AUC value from each nodal metric showing significant differences in intergroup comparisons was used to perform partial correlation analysis with clinical motor score (i.e. UPDRS motor score) first. Age, gender and head motion were regressed out. Then, considering the heterogeneity of motor impairment in PD, motor subscale, like akinesia/rigidity score and tremor score, was separately input into partial correlation analysis with additionally regressing out the other motor subscale. These analyses were conducted by using IBM SPSS 19.0.
(1) Since no existing study had investigated the effect of oscillation frequencies on large-scale networks, previous documents detected that on a local measurement of brain function (amplitude of low-frequency fluctuation, ALFF; or fractional ALFF, fALFF) [22, 23, 31,32,33]. Therefore, in the present study, we replicated their results to confirm the fluctuations of regional measurement on different oscillations in our database, which could help compare the oscillatory effect on large-scale network. In the present study, we first used the unfiltered preprocessed data to calculate the fALFF which has an improved sensitivity and specificity in detecting spontaneous brain activities compared with ALFF [45, 46]. For a timeseries, ALFF is calculated as the sum of amplitudes within a specific low frequency range. fALFF is the ALFF of a given frequency band expressed as a fraction of the sum of amplitudes across the entire frequency range detectable in a given signal. Consequently, the fALFF from each oscillation frequency (slow-5, slow-4 and slow-3) was computed. Finally, z transformation was applied to improve data distributions for parametric statistical analysis. All these steps were completed in DAPBI (http://rfmri.org/dpabi) . Multiple comparisons correction was performed using false discovery rate (FDR) p < 0.05 with an extending cluster size > 10.
(2) Given that the constructions of large-scale functional networks from each oscillation frequency could provide sensitive and comprehensive measurements, as a control, we repeated the calculation of functional network within the commonly used frequency band (0.01–0.1 Hz) and compared the nodal properties among the groups.
(3) Though all of the subjects recruited in the present study were cognitively normal and the intergroup difference of MMSE did not survive after education regression, PD patients (in particular for MPD patients) had a lower MMSE score than normal controls. Therefore, we reworked graph theory-based network analysis and partial correlation analysis by regressing out MMSE score as an extra covariate.
(4) To validate the robustness of network properties at individual level, we first performed logistic regression to identify the most important variables that having ability to discriminate PD patients from controls or between EPD and MPD. To fully identify all the relevant variables, the initial valuables input into logistic regression model were composed of two parts: i. those nodal properties with statistical intergroup difference (p < 0.009); ii. nodal properties contralateral to them with a trend to be statistically significant (p < 0.05) because they might have similar discriminative ability, otherwise these variables with low importance could be removed by logistic regression. Then, receiver operating characteristic (ROC) curve was used to plot the compositive score (probability) from the logistic regression model. Some indices, e.g. sensitivity, specificity and AUC value, were calculated, and Youden index was used to determine the best discriminative result.
Demographic and clinical information
Though potential cognitively impaired subjects were excluded in the present study, compared with normal controls, PD patients had decreased MMSE score (p = 0.001) and education (p = 0.002). In detail, we observed that both MMSE score and education score in MPD patients were lower than that in EPD patients (p = 0.004 and p = 0.005, respectively) and normal controls (p < 0.001 for each). Of note, such difference did not survive after education was regressed out as a covariate of no interest (Table 1). Among EPD, MPD and normal controls, MPD patients were older than EPD patients (p = 0.005). Finally, except for the significant differences of disease severity (e.g. UPDRS motor score, akinesia/rigidity score, tremor score and Hoehn-Yahr stage) between/among groups, no other significant difference was observed. Table 1 showed the specified information.
The effect of oscillation frequency on network properties
Figure 1 visualized the alterations of nodal metrics among three oscillation frequencies and Additional file 1 showed the detailed statistical information (exact p values and t values). For degree centrality, in both PD patients and normal controls, higher values were observed in the frontal lobe (excluding frontal pole), inferior temporal gyrus, temporal fusiform cortex and subcortex (e.g. basal ganglia, thalamus etc.) and lower values were observed in the regions mainly located in parietal, middle and superior temporal and occipital lobes within low oscillation frequencies (slow-5 and slow-4) than within high oscillation frequency (slow-3). Moreover, compared with normal controls, substantial perturbations of network properties calculated in each oscillation frequency were visually observable in PD patients. Similar alterations could be observed in nodal efficiency. For nodal betweenness centrality, the nodal alterations among three oscillation frequencies distributed more sparsely compared with degree centrality and nodal efficiency though a similar trend could still be found. In summary, the anteroventral neocortex (excluding frontal pole) and deep grey matter had high values within the low oscillation frequencies (slow-5 and slow-4) while the posterolateral neocortex had high values within the relatively high oscillation frequency (slow-3), and PD patients suffered from more observable widespread network perturbations across different oscillation frequencies compared with normal controls.
Oscillation-specific nodal alterations in parkinsonian status
For nodal degree centrality (Fig. 2 and Table 2), compared with normal controls, significantly increased degree centrality in basal ganglia (bilateral putamen and right thalamus) and decreased degree centrality in left occipital pole were observed within the slow-5 oscillation frequency in PD patients. Within slow-4 oscillation frequency, PD patients had a significantly decreased degree centrality in left inferior frontal gyrus and right occipital pole and increased degree centrality in left accumbens and right frontal orbital cortex. There was a single region (left frontal medial cortex) showing increased degree centrality in PD patients within slow-3 oscillation frequency. For nodal efficiency (Fig. 3 and Table 3), within slow-5 oscillation frequency, significantly increased efficiency was observed in bilateral putamen, right brain-stem and right angular gyrus while reduced efficiency was observed in left occipital pole in PD patients compared with controls. Within slow-4 oscillation frequency, significantly increased efficiency in left accumbens and frontal orbital cortex and decreased efficiency in left inferior frontal gyrus were observed in PD patients compared with controls. No significant difference was found in nodal efficiency within slow-3 oscillation frequency and in betweenness centrality within all three oscillation frequencies.
Further, among three oscillation frequencies, the only node showing both lower degree centrality and efficiency in EPD than those in controls was left inferior frontal gyrus within slow-4. Severely, widespread alterations of degree centrality and efficiency located in the subcortex (e.g. putamen, pallidum, accumbens and thalamus,), frontal lobe (e.g. frontal orbital cortex), temporal lobe (e.g. middle temporal gyrus), parietal lobe (e.g. angular gyrus) and occipital lobe (e.g. occipital pole and lateral occipital cortex superior division) were investigated in MPD patients.
In brief, nodal alterations were more widespread in MPD patients than in EPD patients. More interestingly, oscillation-specific alterations could be observed in PD patients in different stages, e.g. enhanced nodal properties in basal ganglia and thalamus were specifically observed within slow-5 oscillation frequency in MPD patients.
Clinical relationships between altered oscillation-specific network and motor severity
Among those intergroup comparisons within each oscillation frequency, we performed partial correlation analyses to obtain clinically related nodal alterations. Though several weak correlations were observed and occipital disruption was detected in the three oscillation frequencies, only the correlations between nodal properties in occipital lobe within slow-3 oscillation frequency and motor severity survived after multiple comparisons correction (Bonferroni approach). Specifically, the nodal degree centrality of bilateral lateral occipital cortices superior division (r = − 0.3266, p = 0.0023 for left side, and r = − 0.3741, p = 0.0004 for right side) and bilateral occipital poles (r = − 0.3708, p = 0.0005 for left side, and r = − 0.4328, p < 0.0001 for right side) had significantly negative correlations with UPDRS motor scores. Further, by splitting UPDRS motor score into tremor and akinesia/rigidity scores, we detected that the above correlations were only related to akinesia/rigidity but not tremor (Fig. 4). Similar correlations were shown in nodal efficiency within slow-3 oscillation frequency (Fig. 5). Of note, the correlation between left occipital pole and akinesia/rigidity score did not survive after Bonferroni correction. No other significant clinical correlation was observed in other nodal metrics.
(1) To a large extent, the effect of oscillation frequencies on fALFF was different from that on nodal properties of large-scale network (see Fig. 1 and Additional file 2). Compared with fALFF within slow-4, fALFF within slow-5 was significantly decreased in subcortex (e.g. basal ganglia, thalamus etc.) and dorsolateral cortex while it was significantly increased in ventral cortex and cerebellum, which was largely consistent with previous studies [22, 31,32,33]. And the fALFF perturbation was visually more widespread in PD patients than that in controls among different oscillation frequencies. While comparing fALFF in low frequencies (slow-5 and slow-4) with high frequency (slow-3), subcortex showed lower fALFF and cortex showed higher fALFF within slow-5 and slow-4 than that within slow-3, which was also largely consistent with one recent investigation .
(2) For the nodal degree centrality, among the groups, all the findings observed in the network constructed in the commonly used frequency band (0.01–0.1 Hz) could be detected but assigned to different networks calculated from corresponding oscillation frequencies. More interestingly, the result of bilateral putaminal dysfunction did not survive in the network with 0.01–0.1 Hz in the PD (MPD) patients (Additional file 3). For the nodal efficiency, enhanced brain-stem was additionally observed in the networks with oscillation frequency of slow-5, while other findings kept similar between two methods (Additional file 4). These findings reflects that oscillation-specific network construction is sensitive in detecting nodal abnormalities and could provide distinctive information of nodal alterations in PD.
(3) After regressing out MMSE score, the results of nodal alterations were largely similar to those obtained without MMSE regression except for some individual fluctuation(s) and these statistical data were shown in Additional files 5 and 6. For correlation analysis, no significant change was observed.
(4) In the discrimination between PD patients and controls, the variables finally entered into the logistic regression model included degree centrality in the left inferior frontal gyrus in the oscillation frequency of slow-4 (X1) and left frontal medial cortex in the oscillation frequency of slow-3 (X2), and nodal efficiency in the right angular gyrus (X3) and left putamen (X4) in the oscillation frequency of slow-5. The sensitivity, specificity, AUC value and p value were 0.830, 0.639, 0.820 and < 0.001, respectively (Fig. 6a). The logistic regression formula was:
In the discrimination between EPD and MPD patients, the variables selected by the logistic regression model were the degree centrality in right occipital pole in the oscillation frequency of slow-5 (X1), right accumbens in the oscillation frequency of slow-3 (X2) and right occipital pole in the oscillation frequency of slow-3 (X3), and the nodal efficiency in the right angular gyrus in the oscillation frequency of slow-4 (X4) and lateral occipital cortex superior division in the frequency of slow-3 (X5). The sensitivity, specificity, AUC value and p value were 0.878, 0.744, 0.869 and p < 0.001, respectively (Fig. 6b). The logistic regression formula was:
In the current study, we aimed to investigate oscillation-specific nodal alterations in PD patients across early stage to middle stage. Here, we had three main findings. First, different from the fALFF, the anteroventral neocortex and subcortex had high nodal properties within low oscillation frequencies while the posterolateral neocortex had high values within the relatively high oscillation frequency, where observable perturbation of nodal properties among three oscillation frequencies was detected in PD. Second, oscillation-specific progressive alterations of nodal properties could be observed in PD patients across early to middle stages, which had the ability to well discriminate PD patients from controls or between EPD and MPD (both AUC values > 0.8) at the individual level. Third, occipital dysfunction within slow-3 oscillation frequency had significant correlations with motor severity which was dominated by akinesia/rigidity.
The coupling of two or more oscillations could provide enhanced combinatorial opportunities for storing complex temporal patterns and optimizing synaptic weights . Although the relationship between the oscillation-specific power distribution and physiological functions is not yet fully understood, it was considered that neuronal properties and cytoarchitectonic complexity may contribute to oscillatory power . In the present study, despite the organization hierarchy of large-scale network was different from that of fALFF, we observed a divergent effect of specific oscillations on brain intrinsic function that the anteroventral neocortex and subcortex had high nodal properties within low oscillation frequencies while the posterolateral neocortex had high values within the relatively high oscillation frequency in both PD patients and controls, which uncovered the intrinsic hierarchies of large-scale network in human brain. Since fALFF reflects the intensity of spontaneous brain activity locally , nodal properties from large-scale network have the advantage to evaluate the large-scale network integration [24, 27]. The different effects of oscillations on the different network measurements, e.g. in comparison with slow-4 oscillation frequency, basal ganglia where long time scales (low frequency) of neuron spiking activity, could be observed  showed increased nodal properties in the graph-theory analysis but decreased values in the fALFF measurement within slow-5 oscillation frequency, probably indicated that human brain network is composed of high-dimension intrinsic self-organization varying from different oscillation frequencies and algorithmic hierarchies . More than that, by visual cross-sectional comparisons, we observed that both distributions of altered nodal properties and fALFF in PD patients were more extensive and significant than that in controls, which showed that network was perturbed in parkinsonian status.
Through analyzing PD patients across early to middle stages, we observed that there were relatively preserved nodal properties in PD patients in the early stage while widespread abnormalities occurred in the middle stage, which was consistently seen within all three oscillation frequencies. Luo et al.  systematically investigated the widespread distribution of abnormal network topology and found functional disruption in temporal-occipital, sensorimotor and orbitofrontal regions in PD patients with Hoehn-Yahr stages = 1 and 2. Though a recent study showed a correlation between nodal property and Hoehn-Yahr stage  in PD, the early alteration and the potential progression of network were not detected. Pathologically, neocortex has been involved following an upward topological sequence with disease progression in symptomatic PD patients . Compatible with them, in the present study, nodal dysfunction of neocortex including orbitofrontal, temporal-occipital and parietal regions were observed in MPD patients. Likewise, growing neuroimaging evidence demonstrated structural atrophy and decreased metabolism in inferior frontal gyrus [50, 51], decreased function in occipital lobe [8, 22, 23], and increased function in frontal medial cortex, frontal orbital cortex, middle temporal gyrus and angular gyrus [19, 20, 26] in PD patients. In brief, with disease exacerbation, the perturbation of large-scale nodal properties became conspicuous, which indicated a potential neural mechanism of PD.
Importantly, we found that some of these alterations were intrinsically oscillation-specific, e.g. the cardinal involved region of PD, basal ganglia and thalamus , showing increased degree centrality and nodal efficiency, could be specifically detected within slow-5 oscillation frequency, which was significantly affected in MPD patients. Consistent with a pilot experiment in our previous study using another independent PD cohort , the dysfunction of basal ganglia and thalamus was not significantly observed in early stage PD but enhanced eigenvector centrality in these regions was detected in whole PD group. High nodal degree centrality and efficiency in basal ganglia and thalamus require hyper metabolism to support . Metabolic analyses consistently reported the hyper metabolism and high perfusion in these regions [53,54,55,56,57], and study also confirmed the specificity of hyper metabolism in these regions in differentiating PD from atypical parkinsonism . These findings have been considered as reflecting the release of the basal ganglia from nigral dopaminergic inhibition . Taken together, overconnected basal ganglia and thalamus observed in the large-scale network provided new information to support their roles in the physiological mechanism of PD, and such observation specifically seen within low oscillation frequency (slow-5) indicated that the dysfunction of basal ganglia and thalamus had a specific rhythm of network organization that was perturbed in parkinsonian status.
Though accumbens was not commonly introduced in neuroimaging studies of PD [17, 18, 29], we did observe enhanced nodal function in accumbens for the first time. It was known that as an important node of ventral striatum, accumbens has projections from dopaminergic neurons of mesolimbic pathway . Functionally, no significant devoid of dopaminergic innervation occurs in parkinsonian accumbens , but its structure would be enlarged under chronic exposure of levodopa treatment , which was accordance with our enhanced accumbens function observed in PD patients in the later disease stage, a majority of whom were undertaking levodopa medication. Therefore, future studies through recruiting drug-naïve PD patients and putting insight into the accumbens alterations may help deepen our finding.
Visual impairment is becoming recognized in PD patients [11, 61], which has a correlation with gait impairment  and akinesia and rigidity . In the present study, disrupted occipital nodal degree centrality and efficiency were observed in MPD patients in all three oscillation frequencies. It has been described that PD patients show a higher dependence on visual information for motor control , which was confirmed by the evidence that visual cueing can improve walking in PD . Though the mechanism of visual-motor loop in PD has not been fully studied , the significant correlations between impaired occipital nodes and motor severity dominated by akinesia/rigidity not tremor specifically within slow-3 oscillation frequency further supported the critical role of visual modulation in parkinsonian akinesia and rigidity. Thus, our result suggested that the oscillation-specific occipital nodes whose disruption accompanied with motor exacerbation might be a critical target to modulate motor controls.
Our study had several limitations. First of all, for the 59 PD patients who were under medication, lasting anti-parkinsonian treatment might influence the network organization though clinical assessments and image scanning were carried out after withholding anti-parkinsonian medicine for overnight. Future investigations to detect the immediate influence of levodopa on oscillation-specific network in PD are needed, which is also in progress as an extension of current study. Second, since it was the first study to detect oscillation-specific alterations of large-scale network in parkinsonian brains across early to middle stages, replicated studies are expecting to confirm our findings including the consideration of PD subtypes, e.g. PD with depression, freezing of gait and rapid eye movement sleep behavior disorder. Third, it would be helpful to fully understand PD pathogenesis by exploring the network alterations in PD patients in the advanced stage, which was not done in our study for the limited population. Finally, to date the commonly used scanning protocols studying low oscillations BOLD signal were lasting 290 s to 480 s (390 s in the present study) [21, 22, 31, 32], which might be influencing the stability of network analysis in the low oscillations, therefore, improved protocol with longer scanning time would be necessary to validate current findings in the low oscillations.
By coupling various oscillation frequencies, the intrinsic hierarchy of functional large-scale network was explored and progressive oscillation-specific perturbations of nodal properties were observed in PD. Clinical motor impairment, dominated by akinesia and rigidity, was significantly linearly influenced by the occipital disruption within slow-3 frequency.
Availability of data and materials
The materials used and/or analyzed during the current study are available from the corresponding author on reasonable request.
The area under curve
Early stage PD
Fractional amplitude of low-frequency fluctuation
Mini-mental state examination
Middle stage PD
- ROC curve:
Receiver operating characteristic curve
Unified Parkinson’s disease rating scale
Lees AJ, Hardy J, Revesz T. Parkinson's disease. LANCET. 2009;373:2055–66.
Zhang ZX, Roman GC, Hong Z, Wu CB, Qu QM, Huang JB, et al. Parkinson's disease in China: prevalence in Beijing, Xian, and Shanghai. LANCET. 2005;365:595–7.
Damier P, Hirsch EC, Agid Y, Graybiel AM. The substantia nigra of the human brain. II. Patterns of loss of dopamine-containing neurons in Parkinson's disease. BRAIN. 1999;122(Pt 8):1437–48.
Fearnley JM, Lees AJ. Ageing and Parkinson's disease: substantia nigra regional selectivity. BRAIN. 1991;114(Pt 5):2283–301.
Obeso JA, Rodriguez-Oroz MC, Rodriguez M, Lanciego JL, Artieda J, Gonzalo N, et al. Pathophysiology of the basal ganglia in Parkinson's disease. Trends Neurosci. 2000;23:S8–S19.
Guan X, Xuan M, Gu Q, Xu X, Huang P, Wang N, et al. Influence of regional iron on the motor impairments of Parkinson's disease: a quantitative susceptibility mapping study. J Magn Reson Imaging. 2017;45:1335–42.
Guan X, Zhang Y, Wei H, Guo T, Zeng Q, Zhou C, et al. Iron-related nigral degeneration influences functional topology mediated by striatal dysfunction in Parkinson's disease. Neurobiol Aging. 2019;75:83–97.
Guan X, Zeng Q, Guo T, Wang J, Xuan M, Gu Q, et al. Disrupted functional connectivity of basal ganglia across tremor-dominant and Akinetic/rigid-dominant Parkinson's disease. Front Aging Neurosci. 2017;9:360.
Rolinski M, Griffanti L, Piccini P, Roussakis AA, Szewczyk-Krolikowski K, Menke RA, et al. Basal ganglia dysfunction in idiopathic REM sleep behaviour disorder parallels that in early Parkinson's disease. BRAIN. 2016;139:2224–34.
Szewczyk-Krolikowski K, Menke RA, Rolinski M, Duff E, Salimi-Khorshidi G, Filippini N, et al. Functional connectivity in the basal ganglia network differentiates PD patients from controls. NEUROLOGY. 2014;83:208–14.
Weil RS, Schrag AE, Warren JD, Crutch SJ, Lees AJ, Morris HR. Visual dysfunction in Parkinson's disease. BRAIN 2016 139(11):2827-2843.
Cerasa A, Novellino F, Quattrone A. Connectivity changes in Parkinson's disease. Curr Neurol Neurosci Rep. 2016;16:91.
Xu X, Guan X, Guo T, Zeng Q, Ye R, Wang J, et al. Brain atrophy and reorganization of structural network in Parkinson's disease with Hemiparkinsonism. Front Hum Neurosci. 2018;12.
Wu T, Hallett M. The cerebellum in Parkinson's disease. BRAIN. 2013;136:696–709.
Zeng Q, Guan X, Law Yan Lun JCF, Shen Z, Guo T, Xuan M, et al. Longitudinal alterations of local spontaneous brain activity in Parkinson’s disease. Neurosci Bull. 2017;33:501–9.
Buzsaki G, Draguhn A. Neuronal oscillations in cortical networks. SCIENCE. 2004;304:1926–9.
Suo X, Lei D, Li N, Cheng L, Chen F, Wang M, et al. Functional brain Connectome and its relation to Hoehn and Yahr stage in Parkinson disease. RADIOLOGY. 2017;285:904–13.
Luo CY, Guo XY, Song W, Chen Q, Cao B, Yang J, et al. Functional connectome assessed using graph theory in drug-naive Parkinson's disease. J Neurol. 2015;262:1557–67.
Choe IH, Yeo S, Chung KC, Kim SH, Lim S. Decreased and increased cerebral regional homogeneity in early Parkinson's disease. Brain Res. 2013;1527:230–7.
Wu T, Long X, Zang Y, Wang L, Hallett M, Li K, et al. Regional homogeneity changes in patients with Parkinson's disease. Hum Brain Mapp. 2009;30:1502–10.
Song X, Hu X, Zhou S, Xu Y, Zhang Y, Yuan Y, et al. Association of specific frequency bands of functional MRI signal oscillations with motor symptoms and depression in Parkinson's disease. Sci Rep. 2015;5:16376.
Hou Y, Wu X, Hallett M, Chan P, Wu T. Frequency-dependent neural activity in Parkinson's disease. Hum Brain Mapp. 2014;35:5815–33.
Zhang J, Wei L, Hu X, Zhang Y, Zhou D, Li C, et al. Specific frequency band of amplitude low-frequency fl uctuation predicts Parkinson's disease. Behav Brain Res. 2013;252:18–23.
Rubinov M, Sporns O. Complex network measures of brain connectivity: uses and interpretations. NEUROIMAGE. 2010;52:1059–69.
Watts DJ, Strogatz SH. Collective dynamics of 'small-world' networks. NATURE. 1998;393:440–2.
Kwak Y, Peltier S, Bohnen NI, Muller ML, Dayalu P, Seidler RD. Altered resting state cortico-striatal connectivity in mild to moderate stage Parkinson's disease. Front Syst Neurosci. 2010;4:143.
Liao X, Vasilakos AV, He Y. Small-world human brain networks: perspectives and challenges. Neurosci Biobehav Rev. 2017;77:286–300.
Wang J, Wang X, Xia M, Liao X, Evans A, He Y. GRETNA: a graph theoretical network analysis toolbox for imaging connectomics. Front Hum Neurosci. 2015;9:386.
Skidmore F, Korenkevych D, Liu Y, He G, Bullmore E, Pardalos PM. Connectivity brain networks based on wavelet correlation analysis in Parkinson fMRI data. Neurosci Lett. 2011;499:47–51.
Utianski RL, Caviness JN, van Straaten EC, Beach TG, Dugger BN, Shill HA, et al. Graph theory network function in Parkinson's disease assessed with electroencephalography. Clin Neurophysiol. 2016;127:2228–36.
Han Y, Wang J, Zhao Z, Min B, Lu J, Li K, et al. Frequency-dependent changes in the amplitude of low-frequency fluctuations in amnestic mild cognitive impairment: a resting-state fMRI study. NEUROIMAGE. 2011;55:287–95.
Zuo XN, Di Martino A, Kelly C, Shehzad ZE, Gee DG, Klein DF, et al. The oscillating brain: complex and reliable. NEUROIMAGE. 2010;49:1432–45.
Hong JY, Muller-Oehring EM, Pfefferbaum A, Sullivan EV, Kwon D, Schulte T. Aberrant blood-oxygen-level-dependent signal oscillations across frequency bands characterize the alcoholic brain. Addict Biol. 2017;23(2):824–35.
Hughes AJ, Daniel SE, Kilford L, Lees AJ. Accuracy of clinical diagnosis of idiopathic Parkinson's disease: a clinico-pathological study of 100 cases. J Neurol Neurosurg Psychiatry. 1992;55:181–4.
Kang GA, Bronstein JM, Masterman DL, Redelings M, Crum JA, Ritz B. Clinical characteristics in early Parkinson's disease in a Central California population-based study. Mov Disord. 2005;20:1133–42.
Zhang MY, Katzman R, Salmon D, Jin H, Cai GJ, Wang ZY, et al. The prevalence of dementia and Alzheimer's disease in Shanghai, China: impact of age, gender, and education. Ann Neurol. 1990;27:428–37.
Katzman R, Zhang MY, Ouang-Ya-Qu WZY, Liu WT, Yu E, et al. A Chinese version of the mini-mental state examination; impact of illiteracy in a Shanghai dementia survey. J Clin Epidemiol. 1988;41:971–8.
Yan CG, Wang XD, Zuo XN, Zang YF. DPABI: data processing & analysis for (resting-state) brain imaging. NEUROINFORMATICS. 2016;14:339–51.
Biswal B, Yetkin FZ, Haughton VM, Hyde JS. Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med. 1995;34:537–41.
Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, Johansen-Berg H, et al. Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage. 2004;23(Suppl 1):S208–19.
Yang X, Liu J, Meng Y, Xia M, Cui Z, Wu X, et al. Network analysis reveals disrupted functional brain circuitry in drug-naive social anxiety disorder. Neuroimage 2017 190:213-223.
Achard S, Bullmore E. Efficiency and cost of economical brain functional networks. PLoS Comput Biol. 2007;3:e17.
Lynall ME, Bassett DS, Kerwin R, McKenna PJ, Kitzbichler M, Muller U, et al. Functional connectivity and brain networks in schizophrenia. J Neurosci. 2010;30:9477–87.
Xia M, Wang J, He Y. BrainNet viewer: a network visualization tool for human brain connectomics. PLoS One. 2013;8:e68910.
Zou QH, Zhu CZ, Yang Y, Zuo XN, Long XY, Cao QJ, et al. An improved approach to detection of amplitude of low-frequency fluctuation (ALFF) for resting-state fMRI: fractional ALFF. J Neurosci Methods. 2008;172:137–41.
Zang YF, He Y, Zhu CZ, Cao QJ, Sui MQ, Liang M, et al. Altered baseline brain activity in children with ADHD revealed by resting-state functional MRI. Brain and Development. 2007;29:83–91.
Baria AT, Baliki MN, Parrish T, Apkarian AV. Anatomical and functional assemblies of brain BOLD oscillations. J Neurosci. 2011;31:7910–9.
Ruskin DN, Bergstrom DA, Kaneoke Y, Patel BN, Twery MJ, Walters JR. Multisecond oscillations in firing rate in the basal ganglia: robust modulation by dopamine receptor activation and anesthesia. J Neurophysiol. 1999;81:2046–55.
Braak H, Del TK, Rub U, de Vos RA, Jansen SE, Braak E. Staging of brain pathology related to sporadic Parkinson's disease. Neurobiol Aging. 2003;24:197–211.
Guimaraes RP, Arci SM, Dagher A, Campos LS, Azevedo P, Piovesana LG, et al. Pattern of reduced functional connectivity and structural abnormalities in Parkinson's disease: an exploratory study. Front Neurol. 2016;7:243.
Maier F, Williamson KL, Tahmasian M, Rochhausen L, Ellereit AL, Prigatano GP, et al. Behavioural and neuroimaging correlates of impaired self-awareness of hypo- and hyperkinesia in Parkinson's disease. CORTEX. 2016;82:35–47.
van den Heuvel MP, Sporns O. Network hubs in the human brain. Trends Cogn Sci. 2013;17:683–96.
Eckert T, Tang C, Eidelberg D. Assessment of the progression of Parkinson's disease: a metabolic network approach. Lancet Neurol. 2007;6:926–32.
Eidelberg D, Moeller JR, Dhawan V, Spetsieris P, Takikawa S, Ishikawa T, et al. The metabolic topography of parkinsonism. J Cereb Blood Flow Metab. 1994;14:783–801.
Ma Y, Huang C, Dyke JP, Pan H, Alsop D, Feigin A, et al. Parkinson's disease spatial covariance pattern: noninvasive quantification with perfusion MRI. J Cereb Blood Flow Metab. 2010;30:505–9.
Tang CC, Poston KL, Eckert T, Feigin A, Frucht S, Gudesblatt M, et al. Differential diagnosis of parkinsonism: a metabolic imaging study using pattern analysis. Lancet Neurol. 2010;9:149–58.
Ko JH, Spetsieris PG, Eidelberg D. Network structure and function in Parkinson's disease. Cereb Cortex. 2017:1–15.
Ikemoto S. Brain reward circuitry beyond the mesolimbic dopamine system: a neurobiological theory. Neurosci Biobehav Rev. 2010;35:129–50.
Villalba RM, Smith Y. Differential structural plasticity of corticostriatal and thalamostriatal axo-spinous synapses in MPTP-treated Parkinsonian monkeys. J Comp Neurol. 2011;519:989–1005.
Funamizu Y, Nishijima H, Ueno T, Ueno S, Mizukami H, Yagihashi S, et al. Morphological dendritic spine changes of medium spiny neurons in the nucleus accumbens in 6-hydroxydopamine-lesioned rats treated with levodopa. Neurosci Res. 2017;121:49–53.
Arrigo A, Calamuneri A, Milardi D, Mormina E, Rania L, Postorino E, et al. Visual system involvement in patients with newly diagnosed Parkinson disease. RADIOLOGY. 2017;285:885–95.
Uc EY, Rizzo M, Anderson SW, Qian S, Rodnitzky RL, Dawson JD. Visual dysfunction in Parkinson disease without dementia. NEUROLOGY. 2005;65:1907–13.
Bronstein AM, Hood JD, Gresty MA, Panagi C. Visual control of balance in cerebellar and parkinsonian syndromes. BRAIN. 1990;113(Pt 3):767–79.
Lewis GN, Byblow WD, Walt SE. Stride length regulation in Parkinson's disease: the use of extrinsic, visual cues. BRAIN. 2000;123(Pt 10):2077–90.
The authors would like to thank other members (e.g. Mr. Victor Han and Mr. Steven Cao) from the Liu’s lab (EECS and Helen Willis Neuroscience Institute, University of California, Berkeley) for their intellectual contributions including language editing. The authors would like to thank all participants included in current research project. Xiaojun Guan would like to personally acknowledge the project of Zhejiang University Academic Award for Outstanding Doctoral Candidates.
This work was supported by the 13th Five-year Plan for National Key Research and Development Program of China (Grant No. 2016YFC1306600), the Fundamental Research Funds for the Central Universities of China (Grant No. 2017XZZX001–01), the 12th Five-year Plan for National Science and Technology Supporting Program of China (Grant No. 2012BAI10B04), the National Natural Science Foundation of China (Grant Nos. 81571654, 81371519 and 81701647), the Cooperative Project by Ministry of Health and Provincial Department (Grant No. 2016149022) and the Projects of Medical and Health Technology Development Program in Zhejiang Province (Grant No. 2015KYB174).
Ethics approval and consent to participate
All PD patients and control subjects signed informed consent forms and all procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
Consent for publication
The authors declare that they have no competing of interests.
Additional file 1.
Statistical data for the distribution of nodal properties among different oscillatory frequencies in PD patients and controls.
Additional file 2.
The effect of oscillatory frequencies on fractional amplitude of low-frequency fluctuation (fALFF) in patients with Parkinson’s disease and normal controls (p < 0.05 FDR corrected with an extending cluster > 10).
Additional file 3.
Oscillation-specific alterations of degree centrality between/among groups in the network constructed from the commonly used frequency (0.01–0.1 Hz).
Additional file 4.
Oscillation-specific alterations of nodal efficiency between/among groups in the network constructed from the commonly used frequency (0.01–0.1 Hz).
Additional file 5.
Oscillation-specific alterations of degree centrality between/among groups after MMSE regression.
Additional file 6.
Oscillation-specific alterations of nodal efficiency between/among groups after MMSE regression.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Guan, X., Guo, T., Zeng, Q. et al. Oscillation-specific nodal alterations in early to middle stages Parkinson’s disease. Transl Neurodegener 8, 36 (2019). https://doi.org/10.1186/s40035-019-0177-5
- Parkinson’s disease
- Functional magnetic resonance imaging
- Oscillation frequency
- Graph theory analysis
- Akinesia and rigidity