# Functional Brain Dysconnectivity in Parkinson's Disease: A 5-Year Longitudinal Study

Sahar Yassine, Ute Gschwandtner, Manon Auffret, Sophie Achard, Marc Verin, Peter Fuhr, Mahmoud Hassan

<figure><img src="image_1.png" /></figure>

## To cite this version:

Sahar Yassine, Ute Gschwandtner, Manon Auffret, Sophie Achard, Marc Verin, et al., Functional Brain Dysconnectivity in Parkinson's Disease: A 5-Year Longitudinal Study. Movement Disorders, 2022, 37 (7), pp.1444-1453. 10.1002/mds.29026. hal-03656851

HAL Id: hal-03656851

https://hal.science/hal-03656851

Submitted on 17 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Functional Brain Dysconnectivity in Parkinson's Disease: A 5-Year Longitudinal Study

Sahar Yassine, MSc,<fnref n="1" />,<fnref n="2" /> ID Ute Gschwandtner, MD,<fnref n="3" /> Manon Auffret, PharmD, PhD,<fnref n="4" /> ID Sophie Achard, PhD,<fnref n="5" />
Marc Verin, MD, PhD,<fnref n="1" />,<fnref n="4" />,<fnref n="6" />,<fnref n="7" />* Peter Fuhr, MD,<fnref n="3" /># and Mahmoud Hassan, PhD<fnref n="8" />,<fnref n="9" />#* ID

**ABSTRACT:** **Background:** Tracking longitudinal functional brain dysconnectivity in Parkinson's disease (PD) is a key element to decoding the underlying physiopathology and understanding PD progression.

**Objectives:** The objectives of this follow-up study were to explore, for the first time, the longitudinal changes in the functional brain networks of PD patients over 5 years and to associate them with their cognitive performance and the lateralization of motor symptoms.

**Methods:** We used a 5-year longitudinal cohort of PD patients (n = 35) who completed motor and non-motor assessments and sequent resting state (RS) high-density electroencephalography (HD-EEG) recordings at three timepoints: baseline (BL), 3 years follow-up (3YFU) and 5 years follow-up (5YFU). We assessed disruptions in frequency-dependent functional networks over the course of the disease and explored their relation to clinical symptomatology.

**Results:** In contrast with HC (n = 32), PD patients showed a gradual connectivity impairment in α2 (10-13 Hz) and β (13–30 Hz) frequency bands. The deterioration

in the global cognitive assessment was strongly correlated with the disconnected networks. These disconnected networks were also associated with the lateralization of motor symptoms, revealing a dominance of the right hemisphere in terms of impaired connections in the left-affected PD patients in contrast to dominance of the left hemisphere in the right-affected PD patients.

**Conclusions:** Taken together, our findings suggest that with disease progression, dysconnectivity in the brain networks in PD can reflect the deterioration of global cognitive deficits and the lateralization of motor symptoms. RS HD-EEG may be an early biomarker of PD motor and non-motor progression. © 2022 The Authors. *Movement Disorders* published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society

**Key Words:** Parkinson's disease; functional brain networks; electroencephalography; follow-up study; movement disorders

This is an open access article under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

*Correspondence to: Marc Verin, Neurology Department, Pontchaillou University Hospital, Rue Henri Le Guilloux 35033 Rennes Cedex 9, France; E-mail: marc.verin@chu-rennes.fr
Mahmoud Hassan, Reykjavik University, Menntavegur 1, 101 Reykjavik, Iceland; E-mail: mahmoud.hassan.work@gmail.com

[Corrections added on 27 May 2022, after first online publication: The degrees for Ute Gschwandtner and Peter Fuhr have been corrected. The affiliation for Sophie Achard has been corrected.]

**Relevant conflicts of interest:** There are no financial conflicts of interest to disclose.

⁸These authors have contributed equally to this work.

**Received:** 31 January 2022; **Revised:** 19 March 2022; **Accepted:** 24 March 2022

Published online 14 April 2022 in Wiley Online Library (wileyonlinelibrary.com). DOI: 10.1002/mds.29026

<footnote n="1">¹Univ Rennes, Inserm, LTSI - U1099, Rennes, F-35000, France</footnote>

<footnote n="2">²NeuroKyma, Rennes, F-35000, France</footnote>

<footnote n="3">³Department of Neurology, Hospitals of the University of Basel, Basel, Switzerland</footnote>

<footnote n="4">⁴Comportement et noyaux gris centraux, EA 4712, CHU Rennes, Rennes, France</footnote>

<footnote n="5">⁵University Grenoble Alpes, CNRS, Inria, Grenoble INP, LJK, 38000 Grenoble, France</footnote>

<footnote n="6">⁶Movement Disorders Unit, Neurology Department, Pontchaillou University Hospital, Rennes, France</footnote>

<footnote n="7">⁷Institut des Neurosciences Cliniques de Rennes (INCR), Rennes, France</footnote>

<footnote n="8">⁸MINDig, Rennes, F-35000, France</footnote>

<footnote n="9">⁹School of Science and Engineering, Reykjavik University, Reykjavik, Iceland</footnote>
Introduction

Parkinson's disease (PD) is one of the most frequent neurodegenerative disorders in the elderly population.¹ Along with its well-known prominent motor symptoms (resting tremor, bradykinesia, rigidity, and postural imbalance), a broad spectrum of non-motor disturbances (cognitive impairment up to dementia, neuro-psychiatric disturbances, autonomic or sleep disorders, etc.) can manifest from early (even prodromal) disease stages onward diminishing the quality of life of the patients.²⁻⁴ Pathologically, the most prevalent hypotheses suggest that the degeneration of dopaminergic neurons in the nigrostriatal system associated with Lewy body inclusions are the leading features in most forms of PD at early stages.⁵,⁶ However, with the disease progression, the neuropathological process may propagate across interconnected networks reaching the cortex and causing functional alterations within and between brain regions.⁷ Still, these neuropathological insights do not sufficiently expound the heterogeneous phenotype⁸,⁹ of the patients and the progression-related changes in their brain activity.

Functional disruptions related to both motor and cognitive deficits at early and advanced PD stages were reported using mainly functional magnetic resonance imaging techniques (fMRI).¹⁰⁻¹³ Moreover, electro/magneto-encephalography (EEG/MEG) are increasingly shown as powerful, cost-effective, and non-invasive electrophysiological techniques, to explore functional brain networks in neurodegenerative disorders.¹⁴⁻¹⁸ In PD, MEG/EEG studies uncovered an increase in the cortico-cortical connectivity, characterizing PD patients from earliest disease clinical stages onward,¹⁹,²⁰ whereas decreases in the functional connectivity, mainly in the frontotemporal networks, were proclaimed to reflect the development of mild cognitive impairment and dementia conditions in advanced stages.¹⁷,²¹,²² Reductions in the weight of the functional connections associated with the cognitive phenotype of PD patients were also revealed in several EEG studies.²³,²⁴ However, an important limitation in these findings lies in the fact that these studies were all achieved in a single time point rather than longitudinally, a key point to ultimately develop biomarkers in PD.

Here, we used resting state high density (HD)-EEG recordings of 77 PD patients that underwent a follow-up study at baseline (BL), 3-year follow-up (3YFU) and 5-year follow-up (5YFU) at the Basel university hospital. Our primary goal was to explore the disruptions in the functional brain networks of those PD patients between the BL, 3YFU, and the 5YFU. We examined their progressive evolution and compared it to that of age-matched healthy control subjects. In addition, we correlated the patient-specific altered network with the

change in the clinical scores. Finally, to link the disrupted networks with the lateralized motor symptoms of the patients, we revealed different altered networks between BL and 5YFU characterizing the progression of the disease in the left-affected (LPD) and right-affected (RPD) patients.

Materials and Methods

Participants

Patients were selected from a longitudinal study cohort of patients with idiopathic PD and healthy controls (HC). Patients were recruited from the outpatient clinic of the Department of Neurology and Neurophysiology of the Hospital of the University of Basel (City of Basel, Switzerland) in the period from 2011 to 2015 based on the following selection criteria: PD according to United Kingdom Parkinson's Disease Society Brain Bank criteria,²⁵ Mini-Mental State Examination (MMSE) above 24/30, no history of vascular and/or demyelinating brain pathology, and sufficient knowledge of the German language. HC were matched in age and education. Included patients underwent neurological, cognitive, and EEG examinations at BL and follow-ups after a mean interval of 3 years and 5 years. The study was approved by the local ethics committees (Ethikkommission beider Basel, Basel; Switzerland; EK 74/09) and all patients gave written informed consent before the study inclusion. Specialists who performed the assessment of the patients were unaware of the details of this study. The main cohort included 77 patients at BL, from which 42 did not complete all the follow-ups examinations (see Table S1 for the detailed demographic and clinical measures of the main cohort). All patients were under dopaminergic medications during the examinations. To be noted, 10 patients used cholinesterase inhibitors. We first analyzed the subgroup of 35 PD patients who underwent all the examinations at BL, 3YFU, and 5YFU, and we then cross-validated the results on the remaining patients. The HC group included 32 participants at BL, from which 21 completed the 3YFU and 3 the 5YFU (see Fig. S1 for details of the subject's study flow). The 35 PD patients were also classified into two subgroups by a specialist based on the dominant side of their motor symptoms in all the three visits: LPD patients (n = 23) and RPD patients (n = 10). The division was done using the lateralized items of the Unified Parkinson's Disease Rating Scale III (UPDRS-III) score (items 20–26). Patients with symmetric symptoms in all visits (n = 2) were excluded from this subgroup division.
## Data Acquisition and Preprocessing

Resting state EEG data were recorded in the afternoons using a 256-channel EEG System (Netstation 300, EGI, Eugene, OR) during continuous 12 minutes. Patients were seated comfortably in a relaxing chair, instructed to close their eyes and relax while staying awake with minimum eye and body movements. An EEG technician present in the recording room controlled for vigilance of the patients. The sampling rate of the signals was set to 1000 Hz. The high density (HD-EEG) signals were segmented into epochs of 40 seconds each and the first epoch from each recording was discarded from the study. They were preprocessed automatically using the open-source toolbox Automagic²⁶ and the first six epochs with best quality metrics were retained for the rest of the analysis. See Appendix S1, method section, for more technical details about the preprocessing.

## Brain Networks Reconstruction

After preprocessing, the functional brain networks were estimated using the EEG source-connectivity method.²⁷ First the dynamics of cortical brain sources were reconstructed by solving the inverse problem. To do so, EEG channels and magnetic resonance imaging (MRI) template (ICBM152) were co-registered, and using the OpenMEEG toolbox,²⁸ a realistic head model was built. The weighted minimum norm estimate (wMNE)²⁹ was used to calculate the regional time series of the 210 cortical regions of interest (ROIs) of the Brainnetome atlas.³⁰ Afterward, the regional time series were filtered in different EEG frequency bands: θ (4–8 Hz), α1 (8–10 Hz), α2 (10–13 Hz), and β (13–30 Hz). For each frequency band, functional connectivity was computed between the reconstructed sources

using the phase locking value (PLV).³¹ We obtained for each participant six dynamic connectivity matrices in each frequency band. Those matrices were ultimately averaged across time and epochs to obtain a single static functional connectivity matrix per frequency band, used in the further analysis.

## Statistical Analysis

The network-based statistic (NBS)³² was used to identify the brain networks of PD patients that are significantly different between BL and 5YFU visits. This approach assumes that subnetworks of disrupted connections are more likely to indicate real alterations than isolated disconnections and has been shown to provide considerably greater statistical power than generic methods. Age, gender, education, levodopa equivalent daily dose (LEDD), and dominant side of motor symptoms were considered as confounding factors in this analysis. More details about the statistical analysis are provided in the Appendix S1, methods section.

## Network Index and Correlation Analysis

We defined a metric called network index (NI), inspired from previous work of Hassan et al,²³ as the average weight (connectivity) of the output significant networks issued from NBS:

$$NI = \frac{\sum_{i}^{N} W_{i}}{N},$$

where $W_i$ represents the connectivity value of the edge i of the significant network from NBS and N represents the total number of edges in this network. This NI was computed for each patient in each of the three visits to quantify the longitudinal disruptions in the networks as

<table><caption>TABLE 1 Demographic and clinical characteristics of the participants at BL, 3YFU, and 5YFU visits</caption><thead><tr><th rowspan="2"></th><th colspan="3">Baseline</th><th colspan="2">3YFU</th><th colspan="2">5YFU</th></tr><tr><th>PD<br>(n = 35)</th><th>HC<br>(n = 32)</th><th>P value</th><th>PD<br>(n = 35)</th><th>HC<br>(n = 21)</th><th>PD<br>(n = 35)</th><th>HC<br>(n = 3)</th></tr></thead><tbody><tr><td>Sex (M/F)</td><td>26/9</td><td>18/14</td><td>0.13</td><td>—</td><td>10/11</td><td>—</td><td>2/1</td></tr><tr><td>Age (y)</td><td>67.4 (8.2)</td><td>65.3 (5.6)</td><td>0.24</td><td>70.4 (8.2)</td><td>68.7 (5.5)</td><td>0.37</td><td>72.5 (8.2)</td><td>65.7 (4.2)</td></tr><tr><td>Education (y)</td><td>15.2 (3.2)</td><td>13.8 (2.9)</td><td>0.07</td><td>—</td><td>13.6 (3)</td><td>—</td><td>—</td><td>11 (2)</td></tr><tr><td>Disease dur. (y)</td><td>4.1 (3.7)</td><td>NA</td><td>—</td><td>7.1 (3.7)</td><td>NA</td><td>—</td><td>9.2 (3.8)</td><td>NA</td></tr><tr><td>LEDD<br>(mg/day)</td><td>555 (430)</td><td>NA</td><td>—</td><td>660 (449)</td><td>NA</td><td>—</td><td>647 (396)</td><td>NA</td></tr><tr><td>UPDRS-III</td><td>13.9 (10.2)</td><td>NA</td><td>—</td><td>18.4 (9.7)</td><td>NA</td><td>—</td><td>18 (12.4)</td><td>NA</td></tr><tr><td>MoCA</td><td>26.2 (2.4)</td><td>26.8 (2.5)</td><td>0.26</td><td>25.1 (3.8)</td><td>27.3 (2.1)</td><td>0.016*</td><td>24.5 (5.5)</td><td>23.3 (4.6)</td></tr></tbody></table>

Values are expressed as mean (standard deviation). *Indicates results are significant.

BL, baseline; 3YFU, 3 year follow-up; 5YFU, 5 year follow-up; PD, Parkinson's disease patients; HC, healthy controls; M/F, male/female; y, years; LEDD, Levodopa equivalent daily dose; UPDRS-III, Unified Parkinson's Disease Rating Scales motor ratings; MoCA, Montreal Cognitive Assessment, NA, Not Applicable.
the disease progressed in time. For the correlation analysis, the Pearson's correlation was computed to evaluate the relationship between the NI and the clinical scores of PD patients at and between different timepoints. All parameters are described in the shared GitHub containing the codes necessary to reproduce the results (https://github.com/yassinesahar/FuncDysconnectivityPD).

# Results

## Participant Characteristics

Table 1 summarized the demographic and clinical characteristics of the PD patients and HC subjects included in this study. Their complete neuropsychological and neuropsychiatric assessments are shown in Table S2. There was no significance difference in age, sex, or duration of formal education between the PD patients and the HC at BL.

## Longitudinal Alterations in Functional Connectivity Networks of PD Patients

The overall pipeline of the study is summarized in Figure 1. We followed a data-driven approach supposing that the networks at BL are different from those at 5YFU. The longitudinal change in the NI for both PD patients and HC was explored for the three timepoints. This includes the intermediate timepoint 3YFU whose networks were not used in the comparisons to confirm the overall tendency in the networks and to evaluate the potential correlations with the clinical scores.

## Decreasing Networks of PD Patients (BL > 5YFU)

At the α2 band (10–13 Hz), results revealed a statistically significant network ($t = 3.2, P = 0.021$, corrected using permutation), where the connectivity at 5YFU was significantly lower than BL. This network included 125 connections and 72 regions located principally within the right hemisphere (88.8% of the edges and 88.9% of the regions). Regions with the highest number of connections (highest degree) were mainly part of the right superior and inferior frontal gyrus, right precentral gyrus and the right precuneus, with a dominance of the frontotemporal (12%), the frontofrontal (11.2%), and frontocentral (11.2%) connections (Fig. 2A). These results are in line with previous studies that linked the decrease in frontotemporal connectivity at the α band with the severity of disease progression.¹⁷,²³ The remarkable dominance of the right hemisphere will be explored later when conducting the analysis on the lateralized subgroups. Furthermore, to better understand the longitudinal changes in the network, the NI was computed for the three visits: BL, 3YFU, and 5YFU. Results showed that the NI undergoes a significant decrease between BL and 5YFU ($P < 0.001$). This decrease was progressive and significant between BL and 3YFU ($P < 0.01$) and likewise between 3YFU and 5YFU ($P < 0.05$). All $P$ values were corrected using Bonferroni for multiple comparisons. This NI applied on the functional connectivity networks of the HC showed no significant longitudinal changes between BL and 3YFU (Fig. 2B).

Regarding the β band (13–30 Hz), the obtained network ($t = 3.7, P = 0.009$, corrected using permutation) comprised 103 connections issued from 58 regions.

<figure><img src="image_3.png" /><figcaption>**FIG. 1.** General description of the analysis. The connectivity networks of the PD patients at BL and 5YFU were compared using the Network-Based Statistics (NBS) to retrieve a subnetwork of significantly hypo (BL > 5Y)/hyper (BL < 5Y) connectivity. From this significant network, a network index (NI) was attributed to each PD patient in each of the three visits to evaluate their progression in time and to correlate their longitudinal change with the change in clinical scores. [Color figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**FIG. 2.** Dysconnectivity networks between BL and 5YFU and their corresponding highest degree regions in (A) α2, (C) β. Violin plot representing the longitudinal change of the network index of the PD patients and HC in (B) α2, (D) β. ***P < 0.001, **P < 0.01, *P < 0.05 (corrected for multiple comparisons using Bonferroni). [Color figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>

right hemisphere was also dominant in terms of number of local edges (56.3%) and regions (58.6%). Highest degree regions were largely located in the right superior and inferior frontal, precentral and occipital cortex with a dominance of the frontotemporal (12.6%), fronto-occipital (9.7%), and frontofrontal (8.7%) connections (Fig. 2C). A significant decrease in the NI was observed between BL and 5YFU (*P* < 0.001). This decrease was more severe in terms of linear slope between 3YFU and 5YFU (*P* < 0.001) than between BL and 3YFU (*P* < 0.05). As for the HC, they showed no significant change in their NI between BL and 3YFU, and because of their small sample size at 5YFU, we were not able to statistically test the decrease observed between 3YFU and 5YFU (Fig. 2D). Results of the other frequency bands, the increasing networks and the cross validation on the entire cohort are provided in the Appendix S1.

## Relationship between the NI and the Global Cognitive and Motor Scores

Results showed a global positive correlation between the NI of PD patients and their global cognitive score represented by the Montreal Cognitive Assessment score (MoCA) at both follow-up timepoints with a statistical significance at the intermediate timepoint 3YFU (*r* = 0.36, *P* < 0.05) (Fig. 3A,B). In addition, we

uncovered a subnetwork issued from the previously found dysconnectivity networks, in which the longitudinal changes in the connectivity of its edges correlates with the longitudinal changes in the MoCA scores between BL and 5YFU. In α2, a subnetwork was revealed comprising 16 connections located mostly in the right hemisphere and its corresponding NI was computed (Fig. 3C). A positive significant correlation was observed between the change in both NI and MoCA between BL and 5YFU (*r* = 0.64, *P* < 0.001). This positive correlation persists also when introducing the NI of the intermediate timepoint (3YFU) reaching a significant level mainly between 3YFU and 5YFU (*r* = 0.40, *P* < 0.05). Concerning the subnetwork of the β band, it involved 40 connections distributed between both hemispheres (Fig. 3E). The longitudinal change of its corresponding NI was highly correlated with the change in MoCA score between BL and 5YFU (*r* = 0.72, *P* < 0.001) as well as between 3YFU and 5YFU (*r* = 0.55, *P* < 0.01). To ensure that the significant correlation observed in the change between BL and 5Y in both bands was not driven by the 2 patients that showed a sharp decrease in their MoCA scores (as it might be shown in the Fig. S4), we computed the correlation after excluding these patients. Results show that the correlation remains significant in α2 (*r* = 0.38, *P* < 0.05) (Fig. 3D) and β (*r* = 0.49, *P* < 0.01) (Fig. 3F). We also found significant correlation between the
<figure><img src="image_5.png" /><figcaption>**FIG. 3.** Longitudinal change in the NI and in the MoCA score of PD patients. Correlation between the NI (issued from the hypo-connectivity networks reported above) and the MoCA score of PD patients at (A) 3YFU and (B) 5YFU. Dysconnectivity subnetworks and corresponding highest degree regions where the longitudinal change in the value of connectivity correlates with the longitudinal change in MoCA score of PD patients in (C) α2 and (E) β. Correlation between the change in NI and the change in MoCA observed between BL and 5YFU (left), 3YFU, and 5YFU (right) in (D) α2 and (F) β. [Color figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>

changes in NI and the changes in motor score (the UPDRS-III) at α2 ($r = -0.38, P < 0.05$) and β ($r = -0.41, P < 0.05$) between 3YFU and 5YFU. Figure and more details are provided in the Appendix S1.

Disrupted Networks of PD Patients with Lateralized Motor Symptoms

The above obtained networks showed a majority of disrupted connections located in the right hemisphere in α2 and β bands. To investigate this asymmetry, we divided the PD patients into two subgroups according to the dominance of their lateralized motor symptoms: LPD and RPD. No significant difference in demographic and clinical tests was observed between both

subgroups at baseline or within the same subgroup across the different visits (Table S3). To retrieve the specific dysconnectivity networks associated with the affected side of motor symptoms, we conducted the same previous analysis using NBS on the two lateralized subgroups of PD patients independently.

For the LPD patients, the α2 band revealed a disrupted network ($t = 3, P = 0.039$, corrected) between BL and 5YFU comprising predominantly regions (71%) and connections (63%) within the right hemisphere. Highest degree regions were located principally in the frontal, occipital, and central lobes of the right hemisphere. The corresponding NI of those LPD patients showed a significant decrease between BL and 3YFU ($P < 0.05$), 3YFU and 5YFU ($P < 0.05$), BL and
<figure><img src="image_6.png" /><figcaption>FIG. 4. Dysconnectivity networks with their corresponding highest degree regions (left) and the longitudinal change of the NI in both LPD and RPD patients (right). (A) Network of the LPD patients in α2, (B) network of the LPD patients in β, (C) network of the RPD patients in α2, (D) network of the RPD patients in β. **P < 0.001, *P < 0.05 (corrected for multiple comparisons using Bonferroni). [Color figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>

5YFU (*P* < 0.001). However, this same NI does not reveal any significant differences between the three timepoints when applied on the RPD patients (Fig. 4A) and between BL and 3YFU when applied on the HC (*P* = 0.65). In contrast, dysconnectivity networks of the RPD patients between BL and 5YFU revealed a dominance of the left hemisphere in terms of altered connections. In α2 (*t* = 3.2, *P* = 0.044, corrected using permutation), the disrupted network encompassed a total of 167 connections and 83 regions. Although the

majority of the edges (52.7%) were located within the left hemisphere with the highest degree regions among the left-frontal and left-central lobes, the right hemisphere comprises 40.1% of the disrupted edges mainly issued from regions in right-central, right-frontal, and right-parietal lobes. The corresponding NI presented a significant decrease between BL and 3YFU (*P* < 0.05) and between BL and 5YFU (*P* < 0.001) in RPD patients. Yet, it did not mark any significant change between timepoints in LPD patients (Fig. 4C) and
between BL and 3YFU in HC (P = 0.08). More details about these results in the βband are described in the Appendix S1.

Discussion

In this longitudinal study, we investigated the evolu-
tion of the functional brain networks over a 5-year
period in PD patients using resting state HD-EEG
recordings. We unveiled disrupted networks in both α2
and β frequency bands where the functional connectiv-
ity decreases progressively between BL, 3YFU, and
5YFU, only in PD patients and not in HC. We also
showed a positive correlation between the loss of func-
tional connectivity in time and the deterioration in the
MoCA score that is used largely to examine the global
cognitive performance of PD patients. Furthermore, we
noticed a dominance of the right hemisphere in terms
of altered connections that led us to inspect more in
depth the relationship between the motor symptoms lat-
eralization of PD patients and the disruptions in their
functional brain networks.

Our main findings revealed disrupted networks in α2 and β bands where the functional connectivity drops in function of time. Several previous cross-sectional and longitudinal studies reported a decrease in the power of both alpha and beta bands in PD patients at early and advanced stages of the disease.<sup>33-36</sup> Likewise, these frequency bands were relevant in highlighting reductions in the functional connectivity associated with the motor and cognitive deteriorations over the course of the disease.<sup>17,22,23,37</sup> Our results are in line with these previous findings. In addition, our network index (computed at the patient-level) derived from these networks has encountered a progressive decrease only in PD patients and not in healthy subjects in both bands, which could express the progression of the disease rather than the progression in age. Note that we obtained no significant differences when we performed the NBS analysis on PD patients between consecutive timepoints.

Regarding the spatial distribution of the disrupted
connections and regions in α2 and β, the right hemi-
sphere was predominant. In particular, the dominance
of the right-frontal and right-central lobes was clear in
both bands because they included most of the highest
degree regions (involving most of the altered connec-
tions). A recent meta-analysis showed that both right-
superior-frontal and right-central gyrus presented a
reduced functional connectivity in cognitively impaired
PD patients compared to HC.<sup>38</sup> Besides, reductions in
the functional connectivity of the right frontal and right
somatomotor-sensory cortex were previously reported
to reflect reductions in cortical thinning in PD
patients.<sup>39,40</sup> Moreover, we unfold here, that the
frontotemporal connections were dominant in terms of

altered connections in α2 and β. Indeed, functional dis-
turbances between the frontal and temporal lobes in the
α band were observed when comparing non-demented
to demented PD patients in several MEG study.<sup>17,22</sup> In
EEG, these patterns of connections were relevant in α2
when comparing cognitively intact to cognitively
impaired PD patients.<sup>23</sup>

We have also investigated the potential relationships
between the lateralization of motor symptoms and the
longitudinal change of the functional brain networks in
PD patients. Although previous studies have reported
changes in the cortical and subcortical areas related to
the asymmetry of the motor symptoms<sup>41,42</sup> and linked
them with the progression of both motor and cognitive
symptoms of the disease,<sup>43-45</sup> to date, no functional
connectivity study has tackled their relationship longi-
tudinally. When we conducted the same analysis based
on NBS on the two subgroups separately, we found
hypo-connectivity networks in both α2 and β bands
predominantly located in the right hemisphere and
characterizing the longitudinal progression of the dis-
ease in the LPD patients in contrast to a dominance of
the left hemisphere in the RPD patients. This contra-
laterality observed between the altered hemisphere and
the side of the symptoms may be in line with previous
studies that showed the correspondence between the
side of the motor symptoms and the dopaminergic neu-
ronal loss in the contralateral substantia nigra.<sup>46,47</sup> Fur-
thermore, we found that the α2 frequency band
revealed disturbances in the right hemisphere in RPD
patients. Notably, when we first compared the evolu-
tion of the original NI over time in LPD and RPD
patients, both subgroups showed significant decrease
between different timepoints (Fig. S6). Therefore, alter-
ations in the right hemisphere could be associated with
the progression of the disease in both subgroups. More-
over, the previous results revealing the dominance of
the right hemisphere in PD patients of the initial group
could be related not only to the outnumber of the LPD
patients over the RPD patients, (N=23 vs. N=10,
respectively) but also to a pertinent pattern of distur-
bances within the right hemisphere that could reflect a
common clinical deterioration in all patients. Another
interpretation for these patterns of disturbances could
be the fact that after more than 5 years of evolution,
and despite the worsening of motor symptoms on the
predominant side, the disease may become bilateral.<sup>47</sup>
This may also explain why the loss of connectivity in
the left hemisphere characterizing the RPD patients in β
band appears to be significant also in LPD patients
(Fig. 4D).

In this study, we used the data of the PD patients that
complete both follow-up examinations to retrieve the
disrupted networks in time, which relatively reduced
the sample size from 77 to 35. However, we cross-
validated our results on the initial cohort to confirm the
hypoconnectivity found between timepoints (Fig. S3 in Appendix S1). Unfortunately, having only three healthy subjects at 5YFU is not representative and prevents us from validating the longitudinal change in the network statistically in HC. We, therefore, only took into consideration the corresponding evolution occurring between BL and 3YFU. The relatively small sample size of RPD patients (n = 10) might also affect the statistical power of this comparison. As for the correlation analysis, we chose the global cognitive score MoCA (that covers several cognitive domains) rather than domain specific tests because our dysconnectivity analysis is more likely to reveal a large-scale network rather than a domain-specific one. For instance, the NI showed no correlation with tests that assessed only the visuospatial abilities. Finally, all patients were under dopaminergic treatments (on state) during both the EEG recordings and the neuropsychological assessments, which may mask the magnitude of their motor symptoms as well as their global cognitive performance and may, ultimately, affect our analysis. To overcome this issue, the affected side of PD patients was computed based on the lateralized items of the UPDRS-III in the three timepoints, whereas the LEDD was considered as confounding in the correlation analysis. Despite these considerations, the effect of the dopaminergic medications could still be present in the measures of functional connectivity as shown by several studies.<sup>48,49</sup>

In conclusion, the current study is the first to date to assess longitudinal changes in the functional brain networks of PD patients with the disease progression using HD-EEG. Our findings suggest that disruptions in the networks may characterize not only the evolution of the disease, but also the evolution in the cognitive performance of the patients and the lateralization of their motor symptoms. Further investigations may lead to the identification of potential neuromarkers assisting in subtyping PD patients and predicting the evolution of their motor and cognitive symptoms in time. ■

**Acknowledgments:** We thank all patients, their families, and caregivers for their participation in the study. We would like to thank also Beltrani Selina Maria for her help in EEG data management and data sharing.

## Data Availability Statement

The data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.

## References

1. Bloem BR, Okun MS, Klein C. Parkinson's disease. The Lancet 2021;397:2284-2303.

2. Poewe W, Seppi K, Tanner CM, et al. Parkinson disease. Nat Rev Dis Primer 2017;3:17013

3. Martinez-Martin P, Rodriguez-Blazquez C, Kurtis MM, Chaudhuri KR, Group, on B. of the N. V. The impact of non-motor symptoms on health-related quality of life of patients with Parkinson's disease. Mov Disord 2011;26:399-406.

4. Chaudhuri KR, Healy DG, Schapira AH. Non-motor symptoms of Parkinson's disease: diagnosis and management. Lancet Neurol 2006;5:235-245.

5. Horsager J, Andersen KB, Knudsen K, et al. Brain-first versus body-first Parkinson's disease: a multimodal imaging case-control study. Brain 2020;143:3077-3088.

6. Dickson DW, Braak H, Duda JE, et al. Neuropathological assessment of Parkinson's disease: refining the diagnostic criteria. Lancet Neurol 2009;8:1150-1157.

7. Prodoehl J, Burciu RG, Vaillancourt DE. Resting state functional magnetic resonance imaging in Parkinson's disease. Curr Neurol Neurosci Rep 2014;14:448

8. Rodriguez-Sanchez F, Rodriguez-Blazquez C, Bielza C, et al. Identifying Parkinson's disease subtypes with motor and non-motor symptoms via model-based multi-partition clustering. Sci Rep 2021;11:23645

9. Mestre TA, Eberly S, Tanner C, Grimes D, Lang AE, Oakes D, Marras C. Reproducibility of data-driven Parkinson's disease subtypes for clinical research. Parkinsonism Relat Disord 2018;56:102-106.

10. Luo C, Qiu C, Guo Z, et al. Disrupted functional brain connectivity in partial epilepsy: a resting-state fMRI study. PLOS ONE 2012;7:e28196

11. Baggio H-C, Segura B, Sala-Llonch R, et al. Cognitive impairment and resting-state network connectivity in Parkinson's disease. Hum Brain Mapp 2015;36:199-212.

12. Tuovinen N, Seppi K, de Pasquale F, et al. The reorganization of functional architecture in the early-stages of Parkinson's disease. Parkinsonism Relat Disord 2018;50:61-68.

13. Fiorenzato E, Strafella AP, Kim J, Schifano R, Weis L, Antonini A, Biundo R. Dynamic functional connectivity changes associated with dementia in Parkinson's disease. Brain 2019;142:2860-2872.

14. Stam CJ. Use of magnetoencephalography (MEG) to study functional brain networks in neurodegenerative disorders. J Neurol Sci 2010;289:128-134.

15. Stam C, Jones B, Nolte G, Breakspear M, Scheltens P. Small-world networks and functional connectivity in Alzheimer's disease. Cereb Cortex 2007;17:92-99.

16. Hata M, Kazui H, Tanaka T, et al. Functional connectivity assessed by resting state EEG correlates with cognitive decline of Alzheimer's disease – an eLORETA study. Clin Neurophysiol 2016;127:1269-1278.

17. Bosboom JLW, Stoffers D, Wolters EC, Stam CJ, Berendse HW. MEG resting state functional connectivity in Parkinson's disease related dementia. J Neural Transm 2008;116:193

18. Fingelkurts AA, Fingelkurts AA, Rytsälä H, Suominen K, Isometsä E, Kähkönen S. Impaired functional connectivity at EEG alpha and theta frequency bands in major depression. Hum Brain Mapp 2007;28:247-261.

19. Stoffers D, Bosboom JLW, Deijen JB, Wolters EC, Stam CJ, Berendse HW. Increased cortico-cortical functional connectivity in early-stage Parkinson's disease: an MEG study. Neuroimage 2008;41:212-222.

20. Sánchez-Dinorín G, Rodríguez-Violante M, Cervantes-Arriaga A, Navarro-Roa C, Ricardo-Garcell J, Rodríguez-Camacho M, Solís-Vivanco R. Frontal functional connectivity and disease duration interactively predict cognitive decline in Parkinson's disease. Clin Neurophysiol 2021;132:510-519.

21. Rucco R, Lardone A, Liparoti M, Lopez ET. Brain networks and cognitive impairment in Parkinson's disease. Brain Connect 2021. https://doi.org/10.1089/brain.2020.0985

22. Ponsen MM, Stam CJ, Bosboom JLW, Berendse HW, Hillebrand A. A three dimensional anatomical view of oscillatory resting-state activity and functional connectivity in Parkinson's disease related dementia: an MEG study using atlas-based beamforming. NeuroImage Clin. 2013;2:95-102. https://pubmed.ncbi.nlm.nih.gov/24179762/
23. Hassan M, Chaton L, Benquet P, et al. Functional connectivity disruptions correlate with cognitive phenotypes in Parkinson's disease. *NeuroImage Clin* 2017;14:591-601.

24. Utianski RL, Caviness JN, van Straaten ECW, et al. Graph theory network function in Parkinson's disease assessed with electroencephalography. *Clin Neurophysiol* 2016;127:2228-2236.

25. Gibb WR, Lees AJ. The relevance of the Lewy body to the pathogenesis of idiopathic Parkinson's disease. *J Neurol Neurosurg Psychiatry* 1988;51:745-752.

26. Pedroni A, Bahreini A, Langer N. Automagic: standardized preprocessing of big EEG data. *NeuroImage* 2019;200:460-473.

27. Hassan M, Wendling F. Electroencephalography source connectivity: aiming for high resolution of brain networks in time and space. *IEEE Signal Process Mag* 2018;35:81-96. https://ieeexplore.ieee.org/document/8351910

28. Gramfort A, Papadopoulo T, Olivi E, Clerc M. OpenMEEG: opensource software for quasistatic bioelectromagnetics. *Biomed Eng Online* 2010;9:45

29. Hämäläinen MS, Ilmoniemi RJ. Interpreting magnetic fields of the brain: minimum norm estimates. *Med Biol Eng Comput* 1994;32: 35-42.

30. Fan L, Li H, Zhuo J, et al. The human Brainnetome atlas: a new brain atlas based on connectional architecture. *Cereb Cortex* N Y NY 2016;26:3508-3526.

31. Lachaux J-P, Rodriguez E, Martinerie J, Varela FJ. Measuring phase synchrony in brain signals. *Hum Brain Mapp* 1999;8:194-208.

32. Zalesky A, Fornito A, Bullmore ET. Network-based statistic: identifying differences in brain networks. *Neuroimage* 2010;53:1197-1207.

33. Han C-X, Wang J, Yi G-S, Che Y-Q. Investigation of EEG abnormalities in the early stage of Parkinson's disease. *Cogn Neurodyn* 2013;7:351-359.

34. Caviness JN, Hentz JG, Evidente VG, et al. Both early and late cognitive dysfunction affects the electroencephalogram in Parkinson's disease. *Parkinsonism Relat Disord* 2007;13:348-354.

35. Bosboom JLW, Stoffers D, Stam CJ, van Dijk BW, Verbunt J, Berendse HW, Wolters EC. Resting state oscillatory brain dynamics in Parkinson's disease: an MEG study. *Clin Neurophysiol* 2006;117: 2521-2531.

36. Olde Dubbelink KTE, Stoffers D, Deijen JB, Twisk JWR, Stam CJ, Berendse HW. Cognitive decline in Parkinson's disease is associated with slowing of resting-state brain activity: a longitudinal study. *Neurobiol Aging* 2013;34:408-418.

37. Olde Dubbelink KTE, Stoffers D, Deijen JB, Twisk JWR, Stam CJ, Hillebrand A, Berendse HW. Resting-state functional connectivity as a marker of disease progression in Parkinson's disease: a longitudinal MEG study. *NeuroImage Clin* 2013;2:612-619.

38. Wolters AF, van de Weijer SCF, Leentjens AFG, Duits AA, Jacobs HIL, Kuijf ML. Resting-state fMRI in Parkinson's disease patients with cognitive impairment: a meta-analysis. *Parkinsonism Relat Disord* 2019;62:16-27.

39. Yau Y, Zeighami Y, Baker TE, et al. Network connectivity determines cortical thinning in early Parkinson's disease progression. *Nat Commun* 2018;9:12

40. Mak E, Su L, Williams GB, et al. Baseline and longitudinal grey matter changes in newly diagnosed Parkinson's disease: ICICLE-PD study. *Brain* 2015;138:2974-2986.

41. Heinrichs-Graham E, Santamaria PM, Gendelman HE, Wilson TW. The cortical signature of symptom laterality in Parkinson's disease. *NeuroImage Clin* 2017;14:433-440.

42. Wu T, Hou Y, Hallett M, Zhang J, Chan P. Lateralization of brain activity pattern during unilateral movement in Parkinson's disease. *Hum Brain Mapp* 2015;36:1878-1891.

43. Haaxma CA, Helmich RCG, Borm GF, Kappelle AC, Horstink MWIM, Bloem BR. Side of symptom onset affects motor dysfunction in Parkinson's disease. *Neuroscience* 2010;170:1282-1285.

44. Verreyt N, Nys GMS, Santens P, Vingerhoets G. Cognitive differences between patients with left-sided and right-sided Parkinson's disease. A review. *Neuropsychol Rev* 2011;21:405-424.

45. Fiorenzato E, Antonini A, Bisiacchi P, Weis L, Biundo R. Asymmetric dopamine transporter loss affects cognitive and motor progression in Parkinson's disease. *Mov Disord* 2021;36:2303-2313.

46. Kempster PA, Gibb WR, Stern GM, Lees AJ. Asymmetry of substantia nigra neuronal loss in Parkinson's disease and its relevance to the mechanism of levodopa related motor fluctuations. *J Neurol Neurosurg Psychiatry* 1989;52:72-76.

47. Djaldetti R, Ziv I, Melamed E. The mystery of motor asymmetry in Parkinson's disease. *Lancet Neurol* 2006;5:796-802.

48. Esposito F, Tessitore A, Giordano A, et al. Rhythm-specific modulation of the sensorimotor network in drug-naive patients with Parkinson's disease by levodopa. *Brain J Neurol* 2013;136:710-725.

49. Silberstein P, Pogosyan A, Kühn AA, et al. Cortico-cortical coupling in Parkinson's disease and its modulation by therapy. *Brain J Neurol* 2005;128:1277-1291.

## Supporting Data

Additional Supporting Information may be found in the online version of this article at the publisher's web-site.