Next Article in Journal
Father–Child Interactions in Preschool Children with ASD: A Systematic Review
Next Article in Special Issue
The Directionality of Fronto-Posterior Brain Connectivity Is Associated with the Degree of Individual Autistic Traits
Previous Article in Journal
A Therapeutic Matrix: Virtual Reality as a Clinical Tool for Spinal Cord Injury-Induced Neuropathic Pain
Previous Article in Special Issue
Changes in Default Mode Network Connectivity in Resting-State fMRI in People with Mild Dementia Receiving Cognitive Stimulation Therapy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reduced Effective Connectivity in the Motor Cortex in Parkinson’s Disease

1
Department of Neuroscience, Section of Rehabilitation, University of Padova, Via Gustiniani 3, 35128 Padova, Italy
2
Department of Information Engineering, University of Padova, Via Gradenigo 6/A, 35131 Padova, Italy
3
Parkinson and Movement Disorders Unit, Study Centre on Neurodegeneration (CESNE), Department of Neuroscience, University of Padova, Via Giustiniani 5, 35121 Padova, Italy
4
Padova Neuroscience Center, University of Padova, Via Orus, 35128 Padova, Italy
*
Author to whom correspondence should be addressed.
Brain Sci. 2021, 11(9), 1200; https://doi.org/10.3390/brainsci11091200
Submission received: 22 July 2021 / Revised: 6 September 2021 / Accepted: 9 September 2021 / Published: 12 September 2021
(This article belongs to the Special Issue Neural Networks and Connectivity among Brain Regions)

Abstract

:
Fast rhythms excess is a hallmark of Parkinson’s Disease (PD). To implement innovative, non-pharmacological, neurostimulation interventions to restore cortical-cortical interactions, we need to understand the neurophysiological mechanisms underlying these phenomena. Here, we investigated effective connectivity on source-level resting-state electroencephalography (EEG) signals in 15 PD participants and 10 healthy controls. First, we fitted multivariate auto-regressive models to the EEG source waveforms. Second, we estimated causal connections using Granger Causality, which provide information on connections’ strength and directionality. Lastly, we sought significant differences connectivity patterns between the two populations characterizing the network graph features—i.e., global efficiency and node strength. Causal brain networks in PD show overall poorer and weaker connections compared to controls quantified as a reduction of global efficiency. Motor areas appear almost isolated, with a strongly impoverished information flow particularly from parietal and occipital cortices. This striking isolation of motor areas may reflect an impaired sensory-motor integration in PD. The identification of defective nodes/edges in PD network may be a biomarker of disease and a potential target for future interventional trials.

1. Introduction

Parkinson’s disease (PD) main clinical features are motor symptoms, which have a causal factor in dopamine depletion in the basal ganglia circuit [1]. This degenerative process translates in abnormally fast oscillatory activity expressed in the cortico-basal loop, with frequencies ranging up to 200–250 Hz [2]. At scalp level, this activity may be only recorded as beta frequency, due to the impedance and attenuation processes to which the travelling signal undergoes while passing through brain tissue, cerebrospinal fluid and bone. Pathological firing of the basal ganglia entrains oscillations within the supplementary motor area, lateral premotor cortex, as well as frontal areas in various electroencephalographic (EEG) frequency bands [3], with beta frequency band being the most affected.
More recently EEG modulation and connectivity profiling have been proposed to investigate PD progression stages [2]. Brain connectivity explores physical connections (structural connectivity), statistical dependencies (functional connectivity), and casual interactions (effective connectivity) between brain areas. In particular, effective connectivity studies the direct influence that one cortical area may exerts over another and describes the direction of the interaction.
A computational approach to estimate effective connectivity is Granger Causality (GC) [4] measured by MultiVariate AutoRegressive (MVAR) models. This interconnected representation of the brain can be conveniently treated as a network—i.e., in Mathematics a graph. Graph analysis constitutes a simple model to decode network structures, capturing their complexity and providing measures to characterize their topological features, functional organization, and/or the global or local properties of each node.
A network perspective may provide an intervention target on defective nodes: neurostimulation techniques, whether invasive or non-invasive, may restore cortico-cortical interactions [5,6,7]. Abnormalities in oscillatory activity reflect alterations in connectivity in the cortico-striato-thalamo-cortical circuit [8]. Impaired connectivity has also been described between cortical areas—e.g., prefrontal and motor cortex [9]. Most of these data stem from magnetic resonance imaging (MRI) studies, using either functional or structural connectivity methods [2]. Granger Causality applied to functional Magnetic Resonance Imaging (fMRI) has been deployed to study cortico-thalamo-striato-cortical connections in PD [10,11], confirming the weaker interaction between the caudate and the motor cortex.
Cortico-cortical interactions and their directionality in PD are less known [12]. MRI data focused almost exclusively on motor areas and demonstrated an increased causal connectivity from the left premotor cortex (PMC) to right primary motor cortex (M1). Applying GC to EEG data exploits a non-invasive, low-cost, and bed-side functional brain imaging technique [5,6,7].
The aim of the present study is to identify causal connectivity among cortical areas using source-waveforms in a cohort of PD individuals. Our hypothesis is that we may able to confirm MRI derived connectivity alterations in PD, with a particular involvement of the motor areas, with EEG-derived data.

2. Materials and Methods

2.1. Participants

Fifteen people with PD (9 M; age in the interval [60–83] years; mean disease duration 6.1 ± 4.6 years; mean L-dopa dose 450 ± 185.2 mg; subjects tested in ON-phase) were recruited from Movement Disorders Unit at the Neurology clinic of the Padova University Hospital (Table 1 and Table 2). All showed a bradykinetic phenotype. Inclusion criteria were: diagnosis of idiopathic PD within the last five years; stable dose of antiparkinsonian therapy for at least four weeks; off-medication motor Hoehn and Yahr 1–2 at enrolment. Exclusion criteria were: psychiatric disorder; benzodiazepine and neuroleptic treatment; Mini Mental State Examination (MMSE) < 23; contraindications to neurostimulation. Subjects were part of a randomized controlled trial to study the effect of neurostimulation and motor rehabilitation [13]. EEG data used for the present analysis were collected at baseline to avoid any interference with intervention. Control data from ten healthy volunteers (3 M; age in the interval [51–71] years) were obtained using the same EEG system. All participants provided written informed consent. The Ethics Committee of Padova University Hospital approved the protocol (protocol n. 3507/AO/15).

2.2. EEG Data

Five minutes of eyes-open resting state EEG signal (32-channels system; BrainAmp 32MRplus, BrainProducts GmbH, Munich, Germany) were acquired using an analogic anti-aliasing band pass-filter at 0.1–1000 Hz and converted from analog to digital using a sampling rate of 5 KHz. The reference was between Fz/Cz and ground anterior to Fz. The data were processed in Matlab R2018b (MathWorks, Natick, MA) using personalized scripts based on EEGLAB toolbox (http://www.sccn.ucsd.edu/eeglab (accessed on 1 September 2020) [14] and in Brainstorm (https://neuroimage.usc.edu/brainstorm (accessed on 1 September 2020) [15].
EEG data were band-pass filtered from 1 to 30 Hz using a even-order linear phase FIR filter avoiding phase-distortion and down-sampled using a sampling rate of 500 Hz by using Brainstorm software. Artefacts—i.e., eye movements, cardiac activity, and scalp muscle contraction—were visually identified and removed using independent component analysis implemented in EEGLAB (Infomax). The EEG data were further processed using a common average reference. Source-reconstructed time-series were obtained by using Brainstorm software [15]. First, a head model was created using a symmetric boundary element method based on the anatomy derived from the ICBM152 brain [16]. Time-series of neuronal activity were reconstructed using minimum norm algorithm. Sources were constrained to the cortex and the reconstructed time-series were projected onto regions of interest (ROIs) as defined by the Brodmann atlas [17], where time-series for voxels within a ROI were averaged and used for connectivity analysis. Eight ROIs were defined merging Brodmann areas (BA) in frontal, central, parietal and occipital regions, grouped separately for left and right hemisphere, defined as follow: central areas—BA 1, 2, 3, 4, 6; frontal areas—BA 8, 9, 10, 11, 44, 45, 46; parietal areas—BA 5, 7, 30, 39, 40, 43; occipital areas—BA 17, 18, 19. Time-series were divided into non-overlapping epochs of 2 s.

MVAR Model and Time-Domain Formulation of Granger Causality Index

Analysis was performed for each epoch using the MultiVariate Granger Causality Matlab® Toolbox (MVGC) [18]. MVAR models represent a parametric method for discrete-time signal description. They extend to multiple time-series the AutoRegressive (AR) approach, which expresses the current value of a time-series as a weighted sum of its previous p values, plus a noise term. Under the assumption of having m stochastic processes y i ( n ) ; i = 1…, m; n = p + 1, …, N (in our study, m = 8 as the number of ROIs, N = 1000 , i.e., the number of samples for each epoch), MVAR models allow to take into account and reconstruct the evolution of each of them in terms of its past history, the past history of all the other signals and a contribution from white noise. Our MVAR model, thus, consists of m linear equations, which can be conveniently written in matrix form as:
y ( n ) = k = 1 p A k · y ( n k ) + u ( n )
where y ( n ) is the vector of the n-th sample of the m processes, which has a deterministic component accompanied by the contribution of u ( n ) that represents the m-dimension vector of input white noise, with zero mean and diagonal covariance matrix Σ u . The term p and the A k matrices are respectively the model order and the matrices of model coefficients. The model order defines the number of previous samples, from all the processes, that contribute to the expression of the current samples y ( n ) ; the weights of these samples are defined by the corresponding elements of the A k matrices, namely A k ( i , j ) , and represents the influence of y j ( n k ) on y i ( n ) . The elements of the A k ( i , j ) matrices, together with the variances of the input noise appearing on the main diagonal of the Σ u matrix, represent the model parameters (their number is equal to 64 × p + 8 ) to be identified on EEG data of the eight ROIs.
According to Granger definition, a cause-effect relationship exists between two time-series y i ( n ) and y j ( n ) , with ij, if the variance of the prediction error of y i ( n ) , estimated with a MVAR model including all the m time-series of y ( n ) (unrestricted model), is lower than the one resulting from a MVAR model from which y j ( n ) is excluded (restricted model). Hence, the Granger causality from y j ( n ) to y i ( n ) can be measured with Geweke’s Granger Causality index, defined as:
G C i j = G C y j y i / y = ln σ i ˜ 2 σ i 2
where σ i ˜ 2 and σ i 2 are the variance of predicted error u i ( n ) for restricted and unrestricted model regression, respectively.
For each subject of each group the analysis was firstly carried out epoch-wise. At first, the order of MVAR model was selected: Akaike parsimony criterion (AIC) was used for this purpose [19]. AIC optimum order was evaluated for each epoch with an upper bound of p = 20. The order that was most frequently selected (p = 9, i.e., 18 ms) was then fixed for all epoch/subjects. Equation (2) was evaluated for all the possible pairs of the eight ROIs defined in this study and the 8 × 8 matrix of GC indices represents causal strength among them—i.e., the (i, j) element of the GC matrix represents the strength of the influence that ROI j exerts on ROI i. F-test with Bonferroni correction was applied to check the statistical significance of each index (p < 0.000892): if the null hypothesis, that error variances were not significantly different, was accepted, the value of the corresponding GC coefficient was forced to zero for the current epoch.
The overall estimate of Granger connection magnitudes at a single subject level was obtained by averaging the GC matrices across all epochs. The percentage of epochs recognized as significant throughout the whole recording was also monitored for each subject. In order to graphically represent only the most relevant connections in the network, a threshold was computed using a surrogate data strategy: on the same dataset, the null distribution was determined using phase randomization—i.e., the dataset was transformed in frequency domain via fast Fourier Transform (FFT), randomly shuffled in order to change phase information and then back transformed in time domain via inverse FFT. This procedure was iterated 100 times, thus obtaining 100 realizations under the null hypothesis, and for each realization GC was computed. The threshold was fixed at the 95th percentile.

2.3. Graph Representation

Once the relevant connections between regions were estimated, the dynamics of global organization of the network was described using graph theory measures [20]. We mathematically represented the network as a graph—i.e., a set of vertices (or nodes) that represent our variables (ROIs) and a set of edges representing effective connections among them. The obtained graph was directed—i.e., the edges had a directionality—and weighted—i.e., the edges had variable weights measured by the magnitude of the Granger coefficients. The graph was exhaustively characterized by its adjacency matrix which shows the existing links among nodes. In our case of weighted directed graph, the adjacency matrix is non-symmetrical, thus allowing the differentiation between the incoming and outgoing connections of a node.
To obtain statistical inference regarding groups’ differences in non-threshold connectivity matrices (GC values), a Mann-Whitney-Wilcoxon test (p < 0.05) was performed. The significant connections, corresponding to p < 0.05, were plotted onto a MRI template using the BrainNet Viewer toolbox [21].

2.4. Graph Analysis

Graph analysis describes the network model properties by quantifying topologies of the network representation. Starting from the matrix, different topological measures can be derived for graph description [22]. To characterize the brain networks, parameters related to global efficiency and node strength were calculated by using an open-source toolbox (Brain Connectivity Toolbox, https://sites.google.com/site/bctnet/Home; accessed date: 1 December 2020). Global efficiency was computed as the inverse of the shortest weighted path length and evaluates the ability of the brain to rapidly combine specialized information from distributed brain regions. Structural networks usually are similarly organized and share a high global efficiency, whereas functional network have weaker connections between modules and consequently a weaker global efficiency [23,24]. Node strength was computed as the sum of weights of links connected to the node: the in-strength is the sum of inward link weights and the out-strength is the sum of outward link weights. In brain networks, this basic feature is not distributed homogeneously across the nodes since some nodes have a very high strength, marking them as putative network hubs [25].
A Mann-Whitney-Wilcoxon test (p < 0.05) was applied to statistically compare global efficiency and node strength between the two groups.

3. Results

Average results from GC analysis are summarized in Figure 1 for healthy controls (left panels) and PD participants (right panels): matrices of GC indices representing causal strength among the eight ROIs (upper panels), matrices of the percentage of epochs recognized as significant (medium panels), brain networks showing the connections significant in at least 79% of epochs (lower panels), where 79% is the threshold selected on surrogate data. Controls showed a higher number of connections; almost all of them, regardless of their origin, were directed to bilateral central regions (the motor cortex: BA 1, 2, 3, 4, 6) (Figure 1). The PD participants’ graph was overall poorer; only the connections from left parietal to left central node exceeded the threshold.
Comparison between the average GC matrices in the two groups (Figure 2) indicated significantly stronger connections (p < 0.05) in controls vs. PD participants for the edges from the left (L) frontal to the right (R) central node; from the R frontal to the R central and R parietal nodes; form L central to R central node; from the R occipital to R and L central, and to R frontal nodes; from L occipital to L central and R frontal nodes; from the L parietal to the R central node; and from the R parietal to the R central node.
Comparing Figure 1e with Figure 2 only 4 connections are not represented in Figure 1e (i.e., from R frontal node to R parietal node; from L occipital node to R frontal node, from R central node to R frontal node, from R occipital node to R frontal node). It has to be noted that the different threshold could affect the comparison: the 79% used to select significant connections within the group may be more stringent than the 0.05 p-value threshold used to select significant connections between the two groups.
Global efficiency, a graph global measure expressed by a single value per network, was higher in controls (0.035) compared to PD participants (0.017) (Figure 3a), although the difference was not statistically significant. Node strength, a node-level measure based on the number of connections that a given node has with any other node of the network, showed higher values both for nodes with in-coming and out-coming information flow in controls (Figure 3b,c). The significant larger information flow to/from the motor nodes was clearly visible in both groups: for the in-coming information flow over R central node (p = 0.025), for the out-coming over R central (p = 0.0407) and R occipital (p = 0.0024) nodes, for the overall information flow over R frontal (p = 0.0327), central (p = 0.0043) and occipital (p = 0.0368) nodes (Figure 3b). The larger information flow to/from the motor nodes was clearly visible in both groups, especially for the in-coming flow (Figure 3b).

4. Discussion

Causal brain networks in PD, computed at a source level from EEG data, show overall poorer and weaker connections, calculated as a reduction of global efficiency, although the difference was not statistically significant, most likely due to the intra-subjects variability. Measuring cortical areas reciprocal interactions, we observed a striking isolation of motor areas in PD, especially for the information flow from the parietal and occipital areas to the motor cortices.
Our results highlight a critical parietal and occipital disconnection in PD and are in line with recent MRI data showing resting-state changes in dynamic functional connectivity and reduced number of transition between states in more severely compromised PD [26]. Other imaging studies suggest alterations of these areas: a Positron Emission Tomography (PET) investigation revealed that PD participants showed significantly lower relative cerebral blood flow in the bilateral occipital and posterior parietal cortex [27]. The parieto-occipital network oversees visual and sensory integration, processing the information to be incorporated into the motor program and action. The strict interplay between motor and visual system in the complex process of motor control appears impaired in PD [28,29] and is one of the core neuropsychological features of PD disease progression.
Impaired sensory-motor integration, reflected by the weakening of connections between parietal and central areas, is also confirmed by clinical signs: altered proprioception—i.e., the capability of recognizing body segment position in space—is considered a potential contributor to some peculiar PD clinical pictures—i.e., Pisa syndrome—and becomes even more evident as in the end stages of disease.
We observed a defective incoming information flow at the right motor area in PD. An interesting hypothesis rests on the differential function of dominant and non-dominant limb [30,31,32]. In reaching movements, the non-dominant arm appears better adapted for achieving accurate final positions and the dominant one for specifying initial trajectory features, such as movement direction and peak acceleration. According to this view, the visuo-spatial integration inflow to the right motor cortex, even at rest, may signal the more detailed visuo-spatial control exerted on the left arm, which reaches targets with more accuracy. In fact, people with PD show hypometria in reach-to-grasp tasks [33], interpreted as a deficit of sequential implementation of complex movements. Recent work reported also a difference in body representation during action in PD, which in fact rests, according to the authors, on complex cognitive integrative functions [34].
We noticed a reduction of network global efficiency in PD. The efficiency of a network is a measure of how efficiently it exchanges information across the whole set of edges. This finding suggest that overall information flow is defective in PD. Similar results have been observed in network topology analysis based on MRI data [35] and correlate with cognitive [36] and motor symptoms [37]. These observations confirm the reliability of GC computation based on EEG source analysis reconstruction; they may pave the way for a more widespread use of EEG data, which are part of clinical practice and low-cost, for this kind of analysis. Supporting this thesis, in our previous study [38], resting state relative power of the same PD group was compared to the relative power of a different group of healthy participants in each frequency band. All PD participants showed a statistically significant increase of beta rhythm over sensorimotor areas and a substantial decrease in delta range compared with healthy subjects (p < 0.05). These results are in line with emerging theories suggesting PD is characterized by excessive synchronization in the beta frequency band. The relationship between beta rhythm and clinical state is less clear but some studies have shown strong correlations with rigidity and bradykinesia at rest [39,40]. The average Montreal Cognitive Assessment (MoCA) score in PD participants was 24.79 [38], suggesting a slight degree of cognitive deficit. We did not analyze a single subject correlation between EEG power spectra and cognitive performance at MoCA. In fact, a progressive increase of slow (theta and delta) rhythms is a well known feature in PD. On this ground, the detection of beta band excess may be considered an even stronger finding considering the likely relative increase in this cohort of slow rhythms.
The abnormal oscillatory activity may be modulated by Deep Brain Stimulation (DBS) to stop the extrapyramidal tremor [41]. Based on the same concept, extradural motor cortex stimulation was evaluated in PD with the objective to entrain the pathological circuit at the cortical site [42]. Lastly, non-invasive brain stimulation techniques, namely repetitive transcranial magnetic stimulation (rTMS) and transcranial direct current/alternating stimulation (tDCS/tACS), have been deployed with the same rationale and to extend the efficacy of cognitive rehabilitation [13,43,44]. Given the results that non-invasive stimulation can reduce abnormal brain activity, results from the EEG connectivity could be potentially used to direct the stimulation and optimize restoration of connectivity.
Since recent work [45,46] provided direct evidence that scalp EEG can sense both cortical and subcortical signals, the promising results of this preliminary study support a future data acquisition with high-density EEG systems (≥64 channels), exploiting individual head models based on magnetic resonance imaging.
The choice of the quantification of causal interactions through the time domain formulation of Granger causality in the source-waveforms domain enabled us to overcome some interpretation caveats and common pitfalls that may arise when performing effective connectivity on scalp EEG [47]. Among them, (1) the large distance between EEG sensors and neural sources, and the spatial blurring effect of the skull on the electrical potential EEG distribution —i.e., volume conduction, which causes spurious correlation values between sensors; (2) the differences in signal-to-noise (SNR) ratio—i.e., impedances are equated as much as possible prior to data recordings, but in practice they cannot be equal for all sensors and stationary during the entire recording between signals, which can drive the asymmetries in the causality estimation; (3) common inputs (not included in the computational model), which can lead to spurious inference. Granger causality in the source domain allows to mitigate the adverse effect of volume conduction, also by discarding the interactions at 0 phase and to putatively provide all the common inputs/generators of EEG scalp data. In our model, signals are a combination of the signal-of-interest and a superimposed noise to minimise the effect of variable SNR.
Eventually, to evaluate how much the system is fault tolerant and how much the communication is efficient, the causal interactions were modelled as the weighted edges of a graph: strength and global efficiency of the nodes—i.e., brain regions—of the graph were estimated and compared between PD participants and controls.
The use of threshold is a typical approach for eliminating the weakest connections. This methodology is still not universally accepted, because it may influence the number of links included in the graph analysis and could affect the indices extracted from the network. We used an automatic method based on surrogate and uncorrelated data that enables us to control for spurious findings.
Our study has some limitations. First, the accuracy of EEG source-waveforms reconstruction depends on a sufficient sampling of the surface potential field and increases as the sampling density and coverage increases [48]. Second, despite the method for source reconstruction is computationally cheap and easy to implement, it is based on ad-hoc assumptions and constraints: different assumptions embedded in the source localization method used to solve the inverse problem might lead to different results. Results might be also influenced by the choice of the brain parcellation and the method used to deal with the dipole orientation of all solution points in the same region of interest [49].
In summary, the use of MVAR models and a consolidated index such as the Granger Causality index in the time domain allowed to assess alteration in brain connectivity in PD participants. Future studies might finalize the analysis by using extended version of Granger causality, e.g., to account for the frequency features of the connections and minimize the volume conduction artefact—e.g., partial Granger causality [50], conditional Granger causality [51,52], directed transfer function [53] or partial directed coherence and their derivatives [54,55].

5. Conclusions

Our results provide evidence of an overall resting-state cortico-cortical reduced connectivity in PD: motor areas appear almost isolated, with a strongly impoverished information flow particularly from parietal and occipital cortices. The identification of defective nodes/edges of a network may provide a therapeutic target for future interventional trials.

Author Contributions

Conceptualization and supervision: A.A., S.M., G.M.T., A.D.F.; Data collection: E.F.; A.D.F.; Data analyses: E.F., M.R., J.R.; Original draft preparation: E.F., M.R.; Funding acquisition: A.D.F., S.M.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the ITALIAN MINISTRY FOR FOREIGN AFFAIRS grant number PGR-01045 - SoftAct (call ”Progetti di Grande Rilevanza Internazionale”).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Ethics Committee of Padova University Hospital (protocol n. 3507/AO/15, 10 September 2015).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent was obtained from the subjects to publish this paper.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to patients’ privacy.

Acknowledgments

The authors thank Leonora Castiglia and Michele Tonellato for data collection.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
PDParkinson’s Disease
EEGElectroEncephaloGraphy
MVARMultiVariate AutoRegressive
MRIMagnetic Resonance Imaging
fMRIfunctional Magnetic Resonance Imaging
PMCPreMotor Cortex
M1Primary Motor Cortex
GCGranger Causality
MMSEMini Mental State Examination
ROIRegion Of Interest
BABrodmann Areas
MVGCMultiVariate Granger Causality
AICAkaike parsimony criterion
FFTFast Fourier Transform
tDCS/tACStranscranial Direct Current/Alternating Stimulation
rTMSrepetitive Transcranial Magnetic Stimulation

References

  1. Sharman, M.; Valabregue, R.; Perlbarg, V.; Marrakchi-Kacem, L.; Vidailhet, M.; Benali, H.; Brice, A.; Lehéricy, S. Parkinson’s disease patients show reduced cortical-subcortical sensorimotor connectivity. Mov. Disord. 2013, 28, 447–454. [Google Scholar] [CrossRef]
  2. Wong, J.K.; Middlebrooks, E.H.; Grewal, S.S.; Almeida, L.; Hess, C.W.; Okun, M.S. A Comprehensive Review of Brain Connectomics and Imaging to Improve Deep Brain Stimulation Outcomes. Mov. Disord. 2020, 35, 741–751. [Google Scholar] [CrossRef] [PubMed]
  3. Cassidy, M.; Mazzone, P.; Oliviero, A.; Insola, A.; Tonali, P.; Di Lazzaro, V.; Brown, P. Movement-related changes in synchronization in the human basal ganglia. Brain 2002, 125, 1235–1246. [Google Scholar] [CrossRef]
  4. Granger, C.W. Investigating causal relations by econometric models and cross-spectral methods. Econom. J. Econom. Soc. 1969, 37, 424–438. [Google Scholar] [CrossRef]
  5. Damborská, A.; Honzírková, E.; Barteček, R.; Hořínková, J.; Fedorová, S.; Ondruš, Š.; Michel, C.M.; Rubega, M. Altered directed functional connectivity of the right amygdala in depression: High-density EEG study. Sci. Rep. 2020, 10, 4398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Del Felice, A.; Magalini, A.; Masiero, S. Brain StimulSlow-oscillatory Transcranial Direct Current Stimulation Modulates Memory in Temporal Lobe Epilepsy by Altering Sleep Spindle Generators: A Possible Rehabilitation Tool. Brain Stimul. 2015, 8, 567–573. [Google Scholar] [CrossRef]
  7. Kühn, A.A.; Kempf, F.; Brücke, C.; Doyle, L.G.; Martinez-Torres, I.; Pogosyan, A.; Trottenberg, T.; Kupsch, A.; Schneider, G.H.; Hariz, M.I.; et al. High-frequency stimulation of the subthalamic nucleus suppresses oscillatory beta activity in patients with Parkinson’s disease in parallel with improvement in motor performance. J. Neurosci. 2008, 28, 6165–6173. [Google Scholar] [CrossRef] [Green Version]
  8. Hammond, C.; Bergman, H.; Brown, P. Pathological synchronization in Parkinson’s disease: Networks, models and treatments. Trends Neurosci. 2007, 30, 357–364. [Google Scholar] [CrossRef]
  9. Nettersheim, F.S.; Loehrer, P.A.; Weber, I.; Jung, F.; Dembek, T.A.; Pelzer, E.A.; Dafsari, H.S.; Huber, C.A.; Tittgemeyer, M.; Timmermann, L. Dopamine substitution alters effective connectivity of cortical prefrontal, premotor, and motor regions during complex bimanual finger movements in Parkinson’s disease. Neuroimage 2019, 190, 118–132. [Google Scholar] [CrossRef]
  10. Florin, E.; Pfeifer, J.; Visser-Vandewalle, V.; Schnitzler, A.; Timmermann, L. Parkinson subtype-specific Granger-causal coupling and coherence frequency in the subthalamic area. Neuroscience 2016, 332, 170–180. [Google Scholar] [CrossRef]
  11. Ghasemi, M.; Mahloojifar, A. Disorganization of Equilibrium Directional Interactions in the Brain Motor Network of Parkinson’s disease: New Insight of Resting State Analysis Using Granger Causality and Graphical Approach. J. Med. Signals Sens. 2013, 3, 69–78. [Google Scholar] [CrossRef]
  12. Hao, L.; Sheng, Z.; Ruijun, W.; Kun, H.Z.; Peng, Z.; Yu, H. Altered Granger causality connectivity within motor-related regions of patients with Parkinson’s disease: A resting-state fMRI study. Neuroradiology 2020, 62, 63–69. [Google Scholar] [CrossRef]
  13. Del Felice, A.; Castiglia, L.; Formaggio, E.; Cattelan, M.; Scarpa, B.; Manganotti, P.; Tenconi, E.; Masiero, S. Personalized transcranial alternating current stimulation (tACS) and physical therapy to treat motor and cognitive symptoms in Parkinson’s disease: A randomized cross-over trial. Neuroimage Clin. 2019, 22, 101768. [Google Scholar] [CrossRef]
  14. Delorme, A.; Makeig, S. EEGLAB: An open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J. Neurosci. Methods 2004, 134, 9–21. [Google Scholar] [CrossRef] [Green Version]
  15. Tadel, F.; Baillet, S.; Mosher, J.C.; Pantazis, D.; Leahy, R.M. Brainstorm: A user-friendly application for MEG/EEG analysis. Comput. Intell. Neurosci. 2011, 2011, 879716. [Google Scholar] [CrossRef] [PubMed]
  16. Mazziotta, J.; Toga, A.; Evans, A.; Fox, P.; Lancaster, J.; Zilles, K.; Woods, R.; Paus, T.; Simpson, G.; Pike, B.; et al. A four-dimensional probabilistic atlas of the human brain. J. Am. Med. Inform. Assoc. 2001, 8, 401–430. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Zilles, K.; Amunts, K.; Brodmann, K. Centenary of Brodmann’s map–conception and fate. Nat. Rev. Neurosci. 2010, 11, 139–145. [Google Scholar] [CrossRef]
  18. Barnett, L.; Seth, A.K. The MVGC multivariate Granger causality toolbox: A new approach to Granger-causal inference. J. Neurosci. Methods 2014, 223, 50–68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–723. [Google Scholar] [CrossRef]
  20. Sporns, O.; Tononi, G.; Edelman, G.M. Theoretical neuroanatomy: Relating anatomical and functional connectivity in graphs and cortical connection matrices. Cereb. Cortex 2000, 10, 127–141. [Google Scholar] [CrossRef] [Green Version]
  21. Xia, M.; Wang, J.; He, Y. BrainNet Viewer: A network visualization tool for human brain connectomics. PLoS ONE 2013, 8, e68910. [Google Scholar] [CrossRef] [Green Version]
  22. Bullmore, E.; Sporns, O. Complex brain networks: Graph theoretical analysis of structural and functional systems. Nat. Rev. Neurosci. 2009, 10, 186–198. [Google Scholar] [CrossRef] [PubMed]
  23. Carboni, M.; Rubega, M.; Iannotti, G.R.; De Stefano, P.; Toscano, G.; Tourbier, S.; Pittau, F.; Hagmann, P.; Momjian, S.; Schaller, K.; et al. The network integration of epileptic activity in relation to surgical outcome. Clin. Neurophysiol. 2019, 130, 2193–2202. [Google Scholar] [CrossRef] [PubMed]
  24. Honey, C.J.; Kötter, R.; Breakspear, M.; Sporns, O. Network structure of cerebral cortex shapes functional connectivity on multiple time scales. Proc. Natl. Acad. Sci. USA 2007, 104, 10240–10245. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Fornito, A.; Zalesky, A.; Bullmore, E. Fundamentals of Brain Network Analysis; Academic Press: Cambridge, MA, USA, 2016. [Google Scholar]
  26. Fiorenzato, E.; Strafella, A.P.; Kim, J.; Schifano, R.; Weis, L.; Antonini, A.; Biundo, R. BrainDynamic functional connectivity changes associated with dementia in Parkinson’s disease. Brain 2019, 142, 2860–2872. [Google Scholar] [CrossRef] [Green Version]
  27. Abe, Y.; Kachi, T.; Kato, T.; Arahata, Y.; Yamada, T.; Washimi, Y.; Iwai, K.; Ito, K.; Yanagisawa, N.; Sobue, G. Occipital hypoperfusion in Parkinson’s disease without dementia: Correlation to impaired cortical visual processing. J. Neurol. Neurosurg. Psychiatry 2003, 74, 419–422. [Google Scholar] [CrossRef] [Green Version]
  28. Kloeters, S.; Hartmann, C.J.; Pundmann, V.D.; Schnitzler, A.; Südmeyer, M.; Lange, J. Impaired perception of human movements in Parkinson’s disease. Behav. Brain Res. 2017, 317, 88–94. [Google Scholar] [CrossRef]
  29. Tropini, G.; Chiang, J.; Wang, Z.J.; Ty, E.; McKeown, M.J. Altered directional connectivity in Parkinson’s disease during performance of a visually guided task. Neuroimage 2011, 56, 2144–2156. [Google Scholar] [CrossRef]
  30. Formaggio, E.; Masiero, S.; Bosco, A.; Izzi, F.; Piccione, F.; Del Felice, A. Quantitative EEG Evaluation During Robot-Assisted Foot Movement. IEEE Trans. Neural Syst. Rehabil. Eng. 2017, 25, 1633–1640. [Google Scholar] [CrossRef] [PubMed]
  31. Molteni, F.; Formaggio, E.; Bosco, A.; Guanziroli, E.; Piccione, F.; Masiero, S.; Del Felice, A. Connectivity Modulation After Exoskeleton-Assisted Gait in Chronic Hemiplegic Stroke Survivors: A Pilot Study. Am. J. Phys. Med. Rehabil. 2020, 99, 694–700. [Google Scholar] [CrossRef]
  32. Schaefer, S.Y.; Haaland, K.Y.; Sainburg, R.L. Ipsilesional motor deficits following stroke reflect hemispheric specializations for movement control. Brain 2007, 130, 2146–2158. [Google Scholar] [CrossRef] [PubMed]
  33. Parma, V.; Zanatto, D.; Straulino, E.; Scaravilli, T.; Castiello, U. Kinematics of the Reach-to-Grasp Movement in Vascular Parkinsonism: A Comparison with Idiopathic Parkinson’s Disease Patients. Front. Neurol. 2014, 5, 75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Scarpina, F.; Magnani, F.G.; Tagini, S.; Priano, L.; Mauro, A.; Sedda, A. Mental representation of the body in action in Parkinson’s disease. Exp. Brain Res. 2019, 237, 2505–2521. [Google Scholar] [CrossRef]
  35. Kim, J.; Criaud, M.; Cho, S.S.; Díez-Cirarda, M.; Mihaescu, A.; Coakeley, S.; Ghadery, C.; Valli, M.; Jacobs, M.F.; Houle, S.; et al. Abnormal intrinsic brain functional network dynamics in Parkinson’s disease. Brain 2017, 140, 2955–2967. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Klobušiaková, P.; Mareček, R.; Fousek, J.; Výtvarová, E.; Rektorová, I. Connectivity Between Brain Networks Dynamically Reflects Cognitive Status of Parkinson’s Disease: A Longitudinal Study. J. Alzheimers Dis. 2019, 67, 971–984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Wang, L.; Wang, M.; Si, Q.; Yuan, Y.; Ma, K.; Gan, C.; Zhang, K. Altered brain structural topological properties in Parkinson’s disease with levodopa-induced dyskinesias. Park. Relat. Disord. 2019, 67, 36–41. [Google Scholar] [CrossRef]
  38. Formaggio, E.; Tonellato, M.; Antonini, A.; Castiglia, L.; Gallo, L.; Manganotti, P.; Masiero, S.; Del Felice, A. Oscillatory EEG-TMS Reactivity in Parkinson Disease. J. Clin. Neurophysiol. Off. Publ. Am. Electroencephalogr. Soc. 2021. [Google Scholar] [CrossRef]
  39. Chen, C.C.; Hsu, Y.T.; Chan, H.L.; Chiou, S.M.; Tu, P.H.; Lee, S.T.; Tsai, C.H.; Lu, C.S.; Brown, P. Complexity of subthalamic 13–35 Hz oscillatory activity directly correlates with clinical impairment in patients with Parkinson’s disease. Exp. Neurol. 2010, 224, 234–240. [Google Scholar] [CrossRef]
  40. Little, S.; Pogosyan, A.; Kuhn, A.; Brown, P. Beta band stability over time correlates with Parkinsonian rigidity and bradykinesia. Exp. Neurol. 2012, 236, 383–388. [Google Scholar] [CrossRef] [Green Version]
  41. Benabid, A.L.; Pollak, P.; Louveau, A.; Henry, S.; de Rougemont, J. Combined (thalamotomy and stimulation) stereotactic surgery of the VIM thalamic nucleus for bilateral Parkinson disease. Appl. Neurophysiol. 1987, 50, 344–346. [Google Scholar] [CrossRef]
  42. Cilia, R.; Landi, A.; Vergani, F.; Sganzerla, E.; Pezzoli, G.; Antonini, A. Extradural motor cortex stimulation in Parkinson’s disease. Mov. Disord. 2007, 22, 111–114. [Google Scholar] [CrossRef] [PubMed]
  43. Biundo, R.; Weis, L.; Antonini, A. tDCS effect on cognitive performance in Parkinson’s disease. Mov. Disord. 2016, 31, 1253–1254. [Google Scholar] [CrossRef] [PubMed]
  44. Delgado-Alvarado, M.; Marano, M.; Santurtún, A.; Urtiaga-Gallano, A.; Tordesillas-Gutierrez, D.; Infante, J. Nonpharmacological, nonsurgical treatments for freezing of gait in Parkinson’s disease: A systematic review. Mov. Disord. 2020, 35, 204–214. [Google Scholar] [CrossRef]
  45. Pizzo, F.; Roehri, N.; Villalon, S.M.; Trébuchon, A.; Chen, S.; Lagarde, S.; Carron, R.; Gavaret, M.; Giusiano, B.; McGonigal, A.; et al. Author Correction: Deep brain activities can be detected with magnetoencephalography. Nat. Commun. 2021, 12, 2566. [Google Scholar] [CrossRef] [PubMed]
  46. Seeber, M.; Cantonas, L.M.; Hoevels, M.; Sesia, T.; Visser-Vandewalle, V.; Michel, C.M. Subcortical electrophysiological activity is detectable with high-density EEG source imaging. Nat. Commun. 2019, 10, 753. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Bastos, A.M.; Schoffelen, J.M. A tutorial review of functional connectivity analysis methods and their interpretational pitfalls. Front. Syst. Neurosci. 2016, 9, 175. [Google Scholar] [CrossRef] [Green Version]
  48. Song, J.; Davey, C.; Poulsen, C.; Luu, P.; Turovets, S.; Anderson, E.; Li, K.; Tucker, D. EEG source localization: Sensor density and head surface coverage. J. Neurosci. Methods 2015, 256, 9–21. [Google Scholar] [CrossRef] [Green Version]
  49. Rubega, M.; Carboni, M.; Seeber, M.; Pascucci, D.; Tourbier, S.; Toscano, G.; Van Mierlo, P.; Hagmann, P.; Plomp, G.; Vulliemoz, S.; et al. Estimating EEG Source Dipole Orientation Based on Singular-value Decomposition for Connectivity Analysis. Brain Topogr. 2019, 32, 704–719. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Guo, S.; Seth, A.K.; Kendrick, K.M.; Zhou, C.; Feng, J. Partial Granger causality–eliminating exogenous inputs and latent variables. J. Neurosci. Methods 2008, 172, 79–93. [Google Scholar] [CrossRef] [Green Version]
  51. Chen, Y.; Bressler, S.L.; Ding, M. Frequency decomposition of conditional Granger causality and application to multivariate neural field potential data. J. Neurosci. Methods 2006, 150, 228–237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Wen, X.; Rangarajan, G.; Ding, M. Multivariate Granger causality: An estimation framework based on factorization of the spectral density matrix. Philos. Trans. A Math. Phys. Eng. Sci. 2013, 371, 20110610. [Google Scholar] [CrossRef] [PubMed]
  53. Blinowska, K.J. Review of the methods of determination of directed connectivity from multichannel data. Med. Biol. Eng. Comput. 2011, 49, 521–529. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Baccalá, L.A.; Sameshima, K. Partial directed coherence: A new concept in neural structure determination. Biol. Cybern. 2001, 84, 463–474. [Google Scholar] [CrossRef]
  55. Pascucci, D.; Rubega, M.; Plomp, G. Modeling time-varying brain networks with a self-tuning optimized Kalman filter. PLoS Comput. Biol. 2020, 16, e1007566. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a,b) Grand-average GC matrices for (a) controls and (b) PD participants. Element (i,j) - where i indicates the row and j the column - represents the weight of the connection from ROI j to ROI i. CR and cL denote central areas of right (R) and left (L) hemisphere, fR and fL frontal areas, pR and pL parietal areas and oR and oL occipital areas. The diagonal elements of the matrix are set to zero, since they represent connections from a ROI to itself (loops). (c,d) Grand-average percentage matrices for (c) controls (d) and PD participants: values in the matrix are reported as the result of the ratio, e.g., 0.77 = 77 100 which corresponds to 77 % . Element (i,j) represents the percentage of epochs where the connection from ROI j to ROI i was recognized as significant throughout the whole recording. (e,f) Connectivity networks for (e) controls and (f) PD participants. Only the most relevant connections are depicted, i.e., the ones recognized as significant in at least 79% of epochs. Circles stand for the eight ROIs and arrows for directed connections. Colorbars represent the magnitude of the brain connections.
Figure 1. (a,b) Grand-average GC matrices for (a) controls and (b) PD participants. Element (i,j) - where i indicates the row and j the column - represents the weight of the connection from ROI j to ROI i. CR and cL denote central areas of right (R) and left (L) hemisphere, fR and fL frontal areas, pR and pL parietal areas and oR and oL occipital areas. The diagonal elements of the matrix are set to zero, since they represent connections from a ROI to itself (loops). (c,d) Grand-average percentage matrices for (c) controls (d) and PD participants: values in the matrix are reported as the result of the ratio, e.g., 0.77 = 77 100 which corresponds to 77 % . Element (i,j) represents the percentage of epochs where the connection from ROI j to ROI i was recognized as significant throughout the whole recording. (e,f) Connectivity networks for (e) controls and (f) PD participants. Only the most relevant connections are depicted, i.e., the ones recognized as significant in at least 79% of epochs. Circles stand for the eight ROIs and arrows for directed connections. Colorbars represent the magnitude of the brain connections.
Brainsci 11 01200 g001
Figure 2. Statistically significantly different edges (p < 0.05) in controls vs. PD participants as calculated from the GC matrices. Edges values are derived from the Mann-Whitney-Wilcoxon test, which compares the non-threshold connectivity matrices between the two groups. Only the p-values of the significant connections (p < 0.05) were plotted onto the MRI template. Colorbar represents p values.
Figure 2. Statistically significantly different edges (p < 0.05) in controls vs. PD participants as calculated from the GC matrices. Edges values are derived from the Mann-Whitney-Wilcoxon test, which compares the non-threshold connectivity matrices between the two groups. Only the p-values of the significant connections (p < 0.05) were plotted onto the MRI template. Colorbar represents p values.
Brainsci 11 01200 g002
Figure 3. (a) Global efficiency and (b) in-going, (c) out-going (d) and overall strength for each node, in controls and PD participants. Indices were computed after applying the threshold at GC matrices for each subject. (ad) represent the grand average. The error bar in (a) represents the standard deviation. The size of the node is proportional to the value of the index.
Figure 3. (a) Global efficiency and (b) in-going, (c) out-going (d) and overall strength for each node, in controls and PD participants. Indices were computed after applying the threshold at GC matrices for each subject. (ad) represent the grand average. The error bar in (a) represents the standard deviation. The size of the node is proportional to the value of the index.
Brainsci 11 01200 g003
Table 1. Demographic characteristics of included participants.
Table 1. Demographic characteristics of included participants.
Subject (S)AgeSexEducation
172 yMale16 y
264 yFemale5 y
373 yMale13 y
479 yMale5 y
580 yFemale13 y
669 yMale8 y
761 yFemale11 y
875 yMale13 y
971 yMale17 y
1063 yMale10 y
1168 yMale17 y
1260 yMale17 y
1365 yFemale5 y
1466 yFemale17 y
1583 yFemale5 y
Table 2. Clinical characteristics of included participants.
Table 2. Clinical characteristics of included participants.
SDurationL-Dopa DoseUPDRS IIIBradykinesiaTremorAxial Symptoms
12 y800 mg36411
22 y450 mg28601
31 y200 mg30201
410 y300 mg31420
53 y300 mg37500
69 y750 mg35102
711 y600 mg36000
818 y600 mg29301
96 y400 mg35411
107 y450 mg32202
112 y300 mg27102
122 y200 mg38211
138 y400 mg32321
146 y600 mg35300
154 y400 mg36301
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Formaggio, E.; Rubega, M.; Rupil, J.; Antonini, A.; Masiero, S.; Toffolo, G.M.; Del Felice, A. Reduced Effective Connectivity in the Motor Cortex in Parkinson’s Disease. Brain Sci. 2021, 11, 1200. https://doi.org/10.3390/brainsci11091200

AMA Style

Formaggio E, Rubega M, Rupil J, Antonini A, Masiero S, Toffolo GM, Del Felice A. Reduced Effective Connectivity in the Motor Cortex in Parkinson’s Disease. Brain Sciences. 2021; 11(9):1200. https://doi.org/10.3390/brainsci11091200

Chicago/Turabian Style

Formaggio, Emanuela, Maria Rubega, Jessica Rupil, Angelo Antonini, Stefano Masiero, Gianna Maria Toffolo, and Alessandra Del Felice. 2021. "Reduced Effective Connectivity in the Motor Cortex in Parkinson’s Disease" Brain Sciences 11, no. 9: 1200. https://doi.org/10.3390/brainsci11091200

APA Style

Formaggio, E., Rubega, M., Rupil, J., Antonini, A., Masiero, S., Toffolo, G. M., & Del Felice, A. (2021). Reduced Effective Connectivity in the Motor Cortex in Parkinson’s Disease. Brain Sciences, 11(9), 1200. https://doi.org/10.3390/brainsci11091200

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop