Different features of the cortical sensorimotor rhythms are uniquely linked to the severity of specific symptoms in Parkinson's disease

Parkinson's disease (PD) is associated with functional changes in the neural activity within the brain's sensorimotor network, which in turn are related to the characteristic motor symptoms in PD. The functional changes in PD are particularly prominent in terms of oscillatory neuronal activity in the characteristic sensorimotor alpha and beta rhythms. However, summaries in terms of alpha or beta power do not capture the full range of the complex dynamic nature of the signals from the somatosensory cortex. This raises the question of how to quantify and summarise the functional changes in such oscillatory features in a manner that captures the relevant disease- and symptom-related neural activity. We investigated the role of spontaneous cortical somatosensory activity in the electrophysiological alpha and beta bands among a cohort of early- to mid-stage PD patients (N=78) and age- and gender-matched healthy controls (N=60) using source reconstructed resting-state magnetoencephalography (MEG) recordings. We quantified the oscillatory features of the neural time series by its oscillatory alpha power, beta power, and 1/f broadband characteristics using power spectral density, and additionally by characterising "burst" properties in the signals. We examined the relationship between the signal features and disease state, age, sex, and cortical thickness. Using multiple regression, we examined the relative contribution of the oscillatory features on the clinical manifestation of motor symptoms in the PD group. Our results show that PD patients differ from healthy controls on several of the oscillatory features, showing higher beta-band power, higher burst amplitude, and steeper 1/f broadband characteristics compared to healthy controls, as well as a steeper age-related decrease in the bursts rate. While there was a high degree of correlation between some of the oscillatory features, several features also appeared functionally separated, showing independent feature-to-symptom relationships. For instance, oscillatory beta power increased with the severity of midline function symptoms, while burst rate decreased with the severity of bradykinesia. Our study shows that quantification of distinct features within the oscillatory sensorimotor neural time series in PD captures different underlying mechanisms related to disease progression and symptom severity, which in turn has a potential for a more individualised and precision-based approach to assessing functional neural changes in PD.

1 2 spontaneous activity in the cortical sensorimotor beta and alpha (8-12Hz) bands 18,19 (but see also 13,20 ). 23 Importantly, there is evidence that the beta-band changes are not in the same direction across the different 24 stages of PD. For example, there are reports of increased cortical beta-band power in the early stages of PD 21 , 25 whereas the later stages are associated with decreased beta-band power. 22 26 The beta-band power is not the only feature of the sensorimotor rhythms that is altered in PD. Several studies 27 have found a shift in the beta-band centre frequency (the frequency at which the power spectrum density 28 peaks in the beta-band) towards a lower frequency in PD patients compared to healthy controls. [23][24][25] The shift 29 towards lower beta-band centre frequency is more pronounced in PD patients with dementia 26-29 and 30 correlates with reduced cognitive ability. 25,30 Notably, the centre frequency shift is detectable already in the 31 early stages of PD 24 , and dopaminergic medication does not appear to affect the centre frequency shift. 31 The 32 changes in beta band power and centre frequency in PD could indicate that different features of the oscillatory 33 beta-band activity reflect different underlying neural functions expressed in the measured sensorimotor 34 signals. Changes in beta-band power could be functionally related to sensorimotor disturbances, and changes 35 in centre frequency could be related to cognitive function. 36 Notably, the characteristics of neuronal oscillatory activity may hold additional information of disease-related 37 changes in PD. Both beta-band power and centre frequency reflect a quantification of power spectral density 38 (PSD). While these features can provide valuable information about disease-related changes in PD, the 39 quantification of a neural time series by the PSD provides a static summary of the oscillatory activity across the 40 entire time series. PSD does not account for inherent dynamics in this activity or changes in the time series on 41 shorter time scales-as is prevalent in neural time series. The beta-band exhibits a great degree of variation 42 over time and contains characteristic high-amplitude "bursts" that last about 50-200 ms, both in the cortical 43 and sub-cortical beta-band. [32][33][34][35] Functionally, the transient bursts appear to play a pivotal role in sensorimotor 44 processing through the basal ganglia-thalamic-cortical network. For instance, the presence of a beta burst in 45 3 the sensorimotor cortex close to a tactile stimulation decreased the likelihood of tactile detection 36 , and the 46 rate of beta bursts is shown to decrease in the time leading up to a movement both in STN 37-39 and in the 47 sensorimotor cortex. 40 48 In PD, quantification of beta-band burst activity from recordings in the STN has shown that beta-burst rate and 49 duration are reduced by dopaminergic medication 41,42 and deep brain stimulation. 35 Furthermore, at the 50 cortical level, PD patients exhibit a decrease in the rate of beta burst compared to healthy controls. 17 This 51 decrease in beta burst rate scale with increased severity of bradykinesia and postural-kinetic tremor symptoms 52 but does not change as an effect of dopaminergic medication. 17 Notably, the burst rate showed a higher 53 sensitivity than PSD beta power for discriminating PD patients from healthy controls, demonstrating that the 54 choice of method for beta-band characterisation directly influences the sensitivity of subsequent analyses. 55 In sum, there is solid evidence of functional changes in the oscillatory sensorimotor neural activity in PD-56 whether assessed as PSD band power, peak shift or burst rate. The measures differ between PD patients and 57 healthy controls and scale with PD symptom severity. Sensorimotor activity measured non-invasively with 58 MEG/EEG contains rich information about the functional state of the sensorimotor system and how this 59 changes in PD. The central challenge is quantifying the measured neural signals to extract the disease's relevant 60 features from the signals, be it the spectral power, centre frequencies, or burst-like features. This is further 61 5 CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

94
The main analyses tested the effect of group (PD patients/healthy controls; Table 1), age, sex, and cortical 95 thickness in the sensorimotor cortex on the features explained in Table 2. Second, we tested for associations 96 between the clinical rating of motor symptoms in the PD group by a multiple regression analysis including the 97 sensorimotor signal features in Table 2 and age, sex, and cortical thickness to regress out the contribution 98 hereof and estimate the relative effect size of each signal feature. Figure 2 displays the PSD of the 99 sensorimotor signal. 100   For the first analysis, we wanted to examine if PSD features differed between PD patients and healthy controls 106 and how this might interact with sex, age and cortical thickness. The model parameters and statistical 107 significance of predictors in the regression analysis for all outcome measures are presented in Table 3. 108

Variable Category
Variable Explanation

Beta power
The maximum peak in the 13-30 Hz band. Estimated as the height of the Gaussian function fitted to the PSD after regressing out the 1/f regression line.

Beta centre frequency (Hz)
The dominant frequency bin in the 13-30 Hz band. Estimated as the mean of the Gaussian function fitted to the PSD after regressing out the 1/f regression line

Alpha power
The maximum peak in the 8-12 Hz band. Estimated as the height of the Gaussian function fitted to the PSD after regressing out the 1/f regression line.

Alpha centre frequency (Hz)
The dominant frequency bin in the 8-12 Hz band. Estimated as the mean of the Gaussian function fitted to the PSD after regressing out the 1/f regression line.

1/f intercept
The intercept of the log-linear regression estimated from the full PSD in the 0.5-40 Hz range.

1/f exponent
The exponent of the log-linear regression-corresponding to the slope of the loglog transformed PSD-estimated from the full PSD in the 0.5-40 Hz range

Rate
The number of burst events in the sensorimotor time series divided by the length of time series.

Duration (ms)
Duration of the burst events defined as the time the time-series is above threshold until the next time-point it drops below the threshold.

Interval (ms)
Time from the mu/beta time-series drops below threshold until the next timepoint it reaches threshold again. Amplitude The maximum amplitude of the mu/beta time-series within one burst event.
8 Table 3: regression coefficients and 95% CI for the regression models of the sensorimotor signal features 109 (Table 2)  The PSD beta power showed significant effects of group (χ 2 (1) = 5.42; p = 0.020), age (χ 2 (1) = 4.00; p=0.46) and 116 the three-way interaction between age, sex, and cortical thickness (χ 2 (1) = 8.05; p = 0.005). The relative effect 117 of group corresponded to an increase in beta power of 23.8% [CI: 3.6:48.5] for PD patients. The relative effect 118 of age corresponded to a relative change of -0.2% per year of age ( Figure 3). The direction of interaction three-119 way interaction between age, sex, and cortical thickness indicated that higher cortical thickness was associated 120 with a steeper decrease in beta power, with a larger decrease in beta power in older ages for males compared 121 to females. 122 2.1.2 Beta centre frequency 123 There were no significant effects of any predictor on beta centre frequency. 124

125
There were no significant effects of any predictor on PSD alpha power. 126 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. Higher cortical thickness was associated with a -3.6% [CI: -12.8:5.5] difference in 1/f exponent per standardised 146 unit increase in cortical thickness (Figure 3). 147 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The analysis of burst per minute showed main effects of age (χ 2 (1) = 5.10; p = 0.024) and sex (χ 2 (1) = 7.86; p = 150 0.038) as well as interactions between group and age (χ 2 (1) = 11.74; p = 0.001), age and sex (χ 2 (1) = 5.39; p = 151 0.020), age and cortical thickness (χ 2 (1) = 6.55; p = 0.010), and sex and cortical thickness (χ 2 (1) = 4.35; p = 152 0.037). 153 The age-related effect corresponded to a lower burst rate of -0.   Upper limbs bradykinesia showed significant effects of burst rate (χ 2 (1) = 9.29; p = 0.002) and alpha centre 192 frequency (χ 2 (1) = 6.14; p = 0.013) as well as a significant effect of the covariate cortical thickness (χ 2 (1) = 6.09; 193 p = 0.014). The negative effect of burst rate means that reduced burst rate was associated with increased 194 symptom rating. 195 There were no significant effects of the main predictors on symptom rating for rest tremor, rigidity, 196 postural/kinetic tremor, nor lower limb bradykinesia, though there was a significant effect of sex on rest 197 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 30, 2021. ; 14 tremor (χ 2 (1) = 7.81; p = 0.005). The standardised regression coefficients of each predictor variable on the 198 motor symptoms measured with MDS-UPDRS-III are presented in Table 4. 199 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 30, 2021. ; https://doi.org/10.1101/2021.06.27.21259592 doi: medRxiv preprint  Table 5 presents the significant correlations between the PSD summaries, burst features, age, cortical thickness 204 and clinical rating scales. There was generally high correlations between the PSD features and burst features. 205 The PSD beta power correlated with the burst rate (r=0.26, p=0.002), burst duration (r=0.20, p=0.022), burst 206 amplitude (r=0.52, p<0.001), and correlated negatively with the burst interval (r=-0.24, p=0.005). The PSD alpha 207 power also correlated with the burst rate (r=0.26 p<0.001), burst amplitude (r=0.32, p<0.001) and showed a 208 negative correlation with the bursts interval (r=-0.23, p=0.011). The PSD centre frequencies in neither the beta-209 band nor alpha-bad showed any significant correlation with burst features. 210 The 1/f broadband characteristic of the PSD showed high degree of correlation with the measures of the 211 oscillatory peaks in the PSD for both the 1/f intercept (beta power: r=0.52, p<0.001; beta centre frequency: r=-212  In this study, we aimed to explore how different features of somatosensory oscillatory neuronal activity 224 differed between PD patients and healthy controls across age and gender, and how these features relate to 225 motor symptoms in PD. In agreement with our primary hypothesis, the PD patients differed from healthy 226 controls on several oscillatory features as expected, showing higher beta-band power, higher burst amplitude, 227 a steeper broadband 1/f slope and exponent compared to healthy controls. Furthermore, the analysis of bursts 228 showed that the burst rate was reduced in PD patients compared to healthy controls, confirming previous 229 results on bursts and PD from our group 17 -here on a cohort with three times as many participants. Notably, 230 our current results show that the reduced burst rate in PD is not a static group-level difference but interacts 231 with age, gender and cortical thickness and results in a steeper reduction in burst rate in PD with age compared 232 to healthy controls. 233 Following our secondary hypothesis, we expected that different oscillatory features would reflect distinct 234 underlying functional neural properties and manifest as different motor symptoms in PD. Examining the inter-235 relationship between the features, we observed significant correlations between most oscillatory features, 236 indicating that we had some redundancy among our measures. For instance, the 1/f intercept showed a 237 significant positive correlation with the 1/f slope and beta-band PSD and a high correlation with burst 238 amplitude. We did, however, also see evidence of independent features. Using regression models to examine 239 which oscillatory features were associated with different symptoms, our results showed that the oscillatory 240 features had distinctive relationships with specific symptom scales: bradykinesia severity was related to the 241 burst rate, and alpha centre frequency (and also to cortical thickness) but midline symptoms were related to 242 burst duration and beta power. 243 While we did see correlations between the band-specific oscillatory power in the PSD and the burst rate and 244 burst amplitude-as expected since increased rate and amplitude would necessarily lead to higher PSD-it 245 was, unexpectedly, the intercept of the broadband 1/f regression that showed the highest correlation with the 246 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 30, 2021. ; https://doi.org/10.1101/2021.06.27.21259592 doi: medRxiv preprint burst features. This study adds to the growing evidence that a narrow-band focus-e.g. the beta-band in 247 isolation-could potentially miss essential aspects of the neural signals. Quantifying only the peaks in the PSD 248 misrepresents the actual oscillatory response at those frequencies as the peaks are influenced by the 249 broadband offset and 1/f decay exponent, so any unaccounted-for systematic differences in either PSD offset 250 or decay exponent can lead to a false conclusion that there is a difference in the oscillatory response. 46 The observed disease-related changes in spontaneous cortical bursts in the form of increased burst power and 261 a more rapid decrease in rate over age for PD patients could reflect inhibited projections along the thalamic-262 cortical pathways caused by disturbances in the dopamine-dependent structures projecting to the cortex. The 263 bursting properties of the cortical sensorimotor neural activity are proposed to occur due to long-range input 264 through the ascending thalamic-cortical connection to the cortex, leading to an increase in the local neural 265 excitation and resulting in a burst of synchronous activity. 34 Interestingly, we did not find significant group 266 differences in burst duration in the current study, supporting the view that the central mechanisms of the 267 cortical bursts are not primarily affected in PD. The sub-cortical beta-band activity is directly influenced by the 268 activity of dopamine-responding neurons, as there is a high density of dopamine-dependent neurons in these 269 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. We explored how age-related differences in cortical sensorimotor neural activity might interact with disease-283 related changes in PD within the current study. Age-related effects on spontaneous sensorimotor activity are 284 commonly dealt with by matching the age distributions of the patient group and the healthy control group-285 usually within a narrow age span. The analysis showed, as expected, age-related differences in the 286 sensorimotor activity-most notably in the beta PSD and for the burst duration, burst interval, and burst rate. 287 The age-related differences in burst rate differed between PD patients and healthy controls, with PD patients 288 showing a more considerable reduction of burst as a function of age than healthy controls. The steeper 289 reduction in burst rate with age in PD seems in accordance with the fact that higher age at PD onset is 290 associated with a faster disease progression and more rapid decline in motor function 54 , though a longitudinal 291 design is needed to confirm the relation between disease progression, reduction in burst rate, and age. 292 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 30, 2021. ;

21
We did not see a significant "slowing" of the beta PSD centre frequency between groups, as reported in several 293 previous studies. 43 An explanation might be that such slowing is more pronounced in PD patients with 294 dementia 26-29 and correlates with cognitive ability. 25 The PD patients in the current study were cognitively well-295 functioning and did not differ in their cognitive ability from the healthy controls. Another explanation is that we 296 focused on the activity in the sensorimotor cortex, whereas the slowing of alpha and beta PSD is usually found 297 in frontal areas and globally throughout the brain. 24,27,30 298 We included measures of cortical thickness within the same region of interest (ROI) from which we extracted 299 the function time-series, as we hypothesised that age-related effects upon the functional measures might be 300 mediated through the age-related structural changes in the cortex. However, despite the negative correlation 301 between age and cortical thickness (as expected), we did not find pervasive evidence that cortical thickness 302 affected the functional measures. Cortical thickness did show an effect on the burst rate through interaction on 303 the age-related effect. The only measure that showed a group-specific difference in cortical thickness was the 304 alpha centre frequency, where a higher cortical thickness was associated with a lower alpha PSD peak. 305 We also included sex to explore if disease-related changes in sensorimotor oscillatory activity differed between 306 males and females, as there are well-documented sex differences in the manifestation of PD. 1,55 Male sex is a 307 risk factor for developing PD, with average incidence ratios of approximately 2:1 male-female ratio across all 308 stages of the disease. 56 The disease debuts on average two years earlier in males than females and differs in 309 the initial manifestation of symptoms, with women more likely to develop tremor specific symptoms and men 310 more likely to develop rigidity. 57 We are not aware of any previous studies that explicitly included sex as a 311 factor in analysing neural oscillations in PD. 312 We found significant differences between males and females in the bursting properties in the sensorimotor 313 signal. The burst amplitude was more prominent for female PD patients than for male PD patients, and the 314 age-related differences in both burst rate and burst intervals differed between males and females. These 315 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The present study quantified the neural time series from the sensorimotor cortex based on pre-defined 324 summary measures of the PSD and burst properties. Our study is a broader approach with more sensorimotor 325 rhythm features than other studies more focused on specific signal features, and we also included more PD 326 patients and healthy controls than similar functional studies. 43 327 However, focusing the analysis on different features within a narrow ROI potentially ignores other types of 328 measurements that might be relevant to understanding the development of PD and motor symptoms; for 329 example, the long-range connectivity between the sensorimotor cortex and other cortical areas and the 330 connections between the sensorimotor cortex and the basal ganglia and thalamus (though the subcortical 331 structures are practically invisible in MEG). Treating the activity in the sensorimotor cortex as single time series 332 also means that we remove the sensitivity to spatial features of the signals, e.g. focal versus spatially blurred 333 activity in one group or the other. If the oscillatory activity extends over a larger cortical surface area, then that 334 signal will also manifest as power differences in the measured signal. 51 There are potentially many more 335 features to be uncovered, and future studies may explore how the PSD-and burst features further interact 336 with other aspects of brain activity in the global function of the brain to fully understand the interaction 337 between functional and structural changes in PD. 338 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 30, 2021. ;

23
We investigated a relatively large cohort of PD patients and healthy controls (for a neuroimaging study) to 339 make meaningful inferences about how age and sex interact with the group level difference between PD 340 patients and healthy controls, but a limitation is that our study is cross-sectional. We aim to follow this cohort 341 longitudinally to estimate the development trajectories of the sensorimotor oscillatory activity in PD compared 342 to healthy ageing. 343 There are myriads of ways to quantify brain dynamics, and finding features of neural signals that can explain 344 disease mechanism or symptoms, even if extracted along with a reduced number of dimensions, will be helpful 345 if they provide adequate information about the disease-or symptom-state. Further characterisation of the 346 association between features in the non-invasive brain signals and motor symptoms can potentially be a 347 valuable tool to aid in the diagnosis and evaluation of treatments. Understanding how features in the neural 348 time series are related to motor symptoms in PD will also help develop non-invasive neural stimulation that can 349 potentially relieve motor symptoms. 35,59 350 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 30, 2021. The study was approved by the regional ethics committee (Etikprövningsnämden Stockholm, DNR 2019-00542) 354 and followed the Declaration of Helsinki. All participants gave written informed consent before participation. 355 The PD patients were recruited from the Parkinson's Outpatient Clinic, Department of Neurology, Karolinska 356 University Hospital, Stockholm, Sweden. The healthy controls were recruited by advertising or amongst 357 spouses of PD patients. 22 participants (18 patients, 4 healthy controls) were included from a previous study 17 358 who were qualified based on the recruitment criteria of the present study and had done the same MEG and 359 MRI procedures as in the present study. All data were reanalysed following the procedure described below. One participant declined to do the MRI scanning, one participant had a scanner malfunction during MRI 368 acquisition, and eleven participants had their MRI scans cancelled and were not included in the analysis. Twp 369 PD patients and 11 healthy controls were excluded. Table 1 is a summary of the participants included in the 370 analysis. 371 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The remaining data length was not significantly different between groups (Wilcoxon rank sum test, p = 0.98). 415 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. Since the narrow-band beta power in the PSD is dependent on the broader features of the broadband 431 spectrum, we further analysed the 1/f broadband characterises of the sensorimotor activity as this could play a 432 role in the functional properties of the beta-band and has been shown to differ between healthy control and 433 PD patients. 17 We used the fitting oscillations & one over f (FOOOF) toolbox 46 to analyse the 1/f broadband 434 characteristic of the PSD (intercept and exponent) and the oscillatory peaks in the canonically defined beta 435 band (13-30 Hz) and alpha band (8-12 Hz). A log-linear regression is fitted to the PSD and subtracted before 436 fitting Gaussian functions to the peaks in the PSD. The midpoint of the Gaussians corresponds to the centre 437 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 30, 2021. ; https://doi.org/10.1101/2021.06.27.21259592 doi: medRxiv preprint 28 frequency (i.e. peak frequency) and the height represents the signal power. A new log-linear function is fitted 438 to the PSD after substracting the Gaussian function to estimate the 1/f characteristic. 439 All participants showed a discernable beta peak in the PSD. Nine PD patients and nine healthy controls did not 440 show a peak in the PSD alpha band (no difference between groups, χ 2 (1) = 0.12; p = 0.73). 441 Each symptom score was analysed by multiple regression and modelled as a function of the burst rate, median 464 burst duration, median bursts interval, median burst amplitude, PSD 1/f intercept, PSD 1/f exponent, PSD beta 465 power, PSD beta centre frequency, PSD alpha power, and PSD alpha centre frequency for each PD patient. The 466 models further included the age, sex, and cortical thickness to regress out the contribution hereof and estimate 467 the relative effect size of each signal feature. All symptom ratings and continuous predictor variables, except 468 age, were z-transformed to get the standardised effect size from the regression models. Significance testing 469 was done by removing one predictor from the model and comparing the variance explained between the full 470 model and the model with a predictor removed by log-likelihood ratio tests. 471

472
Since the features analysed in the main analysis are all extracted from the same signal, we wanted to explore 473 the mutual relation between variables. We computed the Pearson correlation between the signal features in 474 Table 2, age, cortical thickness, MoCA, and MDS-UPDRS-III subscales for all pairwise complete cases. 475

476
Parts of the data used in this analysis will be made available at a public data repository. Scripts for running the 477 analysis presented in the paper are available at www.github.com/mcvinding/PD_beta_bursts2. 478 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.  is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 30, 2021. Tables   1   Table 1 2 Group-level summary of the participants included in the analysis. Mean (standard deviation). 3  Table 2 4 Table 2: Explanation of the main outcome variables in the analysis 5 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 30, 2021. ; https://doi.org/10.1101/2021.06.27.21259592 doi: medRxiv preprint Table 3 7 Regression coefficients and 95% CI for the regression models of the sensorimotor signal features (Table 2) with  8 Group, Age, Sex, and Cortical Thickness. Values in red indicate coefficients of statistically significant factors in 9 the model comparison. LL: lower limit, UL: upper limit. 10

Variable Category
Variable Explanation

Beta power
The maximum peak in the 13-30 Hz band. Estimated as the height of the Gaussian function fitted to the PSD after regressing out the 1/f regression line.

Beta centre frequency (Hz)
The dominant frequency bin in the 13-30 Hz band. Estimated as the mean of the Gaussian function fitted to the PSD after regressing out the 1/f regression line

Alpha power
The maximum peak in the 8-12 Hz band. Estimated as the height of the Gaussian function fitted to the PSD after regressing out the 1/f regression line.

Alpha centre frequency (Hz)
The dominant frequency bin in the 8-12 Hz band. Estimated as the mean of the Gaussian function fitted to the PSD after regressing out the 1/f regression line.

1/f intercept
The intercept of the log-linear regression estimated from the full PSD in the 0.5-40 Hz range.

1/f exponent
The exponent of the log-linear regression-corresponding to the slope of the loglog transformed PSD-estimated from the full PSD in the 0.5-40 Hz range

Rate
The number of burst events in the sensorimotor time series divided by the length of time series.

Duration (ms)
Duration of the burst events defined as the time the time-series is above threshold until the next time-point it drops below the threshold.

Interval (ms)
Time from the mu/beta time-series drops below threshold until the next timepoint it reaches threshold again. Amplitude The maximum amplitude of the mu/beta time-series within one burst event.
PSD Bursts . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review)  Table 4 12 Standardised regression coefficients (95% CI) for the six regression models on motor symptoms measured with 13 the MDS-UPDRS-III. Values in bold indicate significant factors in the model comparison. 14 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted June 30, 2021. is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted June 30, 2021. ; https://doi.org/10.1101/2021.06.27.21259592 doi: medRxiv preprint  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted June 30, 2021. ; https://doi.org/10.1101/2021.06.27.21259592 doi: medRxiv preprint