Introduction

Beta oscillations within cortical sensorimotor areas and the basal ganglia have been proposed to play a role in movement initiation, termination, and inhibition (Aron & Poldrack, 2006; Benis et al., 2014; Jurkiewicz, Gaetz, Bostan, & Cheyne, 2006; J. R. Wessel, 2020). Altered beta band activity has been strongly linked to motor impairment in Parkinson’s disease (PD), demonstrating its relevance to proper motor performance (Brown et al., 2001; Cassidy et al., 2002; Tinkhauser et al., 2017). The beta rhythm is often interpreted as reflecting the status quo (Engel & Fries, 2010), that is, active maintenance or stabilization of current motor or cognitive output to attenuate alternatives and distractions (Espenhahn, de Berker, van Wijk, Rossiter, & Ward, 2017; Fischer et al., 2019). The basal ganglia keep cortex under inhibitory control (Bonnevie & Zaghloul, 2018) which, similarly to releasing a break in a car, must be removed to change the current motor state (Alegre et al., 2013).

Starting and stopping of movement have mostly been studied with variations of the Stop Signal Task and the Go/No-Go Task (Alegre et al., 2013; Aron & Poldrack, 2006; Ray et al., 2012) that require participants to perform simple, ballistic movements and inhibit them occasionally. Shortly before and during movement, beta oscillations are suppressed (beta suppression), reflecting a task-related active state of the motor network (Jurkiewicz et al., 2006). In contrast, beta power transiently increases above baseline levels after movement termination (beta rebound) (Fonken et al., 2016; Ray et al., 2012), indicating inhibition (Salmelin, Hämäläinen, Kajola, & Hari, 1995; Schmidt & Berke, 2017) and motor adaptation processes (Struber, Baumont, Barraud, Nougier, & Cignetti, 2021; Tan et al., 2014). Notably, whether these modulations are causally involved in motor control is still under debate (Pfurtscheller, Neuper, Brunner, & da Silva, 2005; Toledo, Manzano, Barela, & Kohn, 2016). Inhibitory processes have been associated with increased beta power or reduced suppression thereof in prefrontal cortical areas (Swann et al., 2009; Wagner, Wessel, Ghahremani, & Aron, 2018), and in the subthalamic nucleus (STN) (Bastin et al., 2014), with some studies reporting correlations with inhibitory success (Benis et al., 2014; W. Chen et al., 2020). Besides playing a critical role in stopping movement (Mosher, Mamelak, Malekmohammadi, Pouratian, & Rutishauser, 2021), the STN seems to be involved in delaying or pausing movement until sufficient evidence in favor of a motor program has accumulated (Ray et al., 2012). Recent evidence demonstrated that the STN is recruited by the prefrontal cortex via the hyperdirect pathway to implement its pausing function (W. Chen et al., 2020; Lofredi et al., 2021; Oswal et al., 2021; J. R. Wessel, Waller, & Greenlee, 2019). However, the role of cortico-subcortical beta synchronization in coordinating movements that are already ongoing remains elusive.

Communication between STN and cortex might become particularly important in tasks involving cognitive factors, such as anticipation. STN beta power has been found to index task complexity and behavioral control (Oswal et al., 2013), proactive inhibition and planning (Benis et al., 2014), non-motor decision making, and working memory (Zavala, Jang, & Zaghloul, 2017), and cue evaluation with respect to behavioral goals (Oswal, Litvak, Sauleau, & Brown, 2012). Yet, the extent to which modulation of beta oscillations in basal-ganglia cortex networks depends on expectation remains unknown.

In the present study, we address these research gaps with a paradigm that involves a rotational movement performed in a continuous fashion with occasional rapid changes in movement direction (reversals), as well as movement initiations and terminations. Accounting for the relevance of basal ganglia-cortical loops in motor control, we recorded cortical and STN oscillatory activity simultaneously in PD patients who had undergone implantation of deep brain stimulation (DBS) electrodes the day before. Patients performed the rotational movements according to visual instructions which were manipulated such that their identity and time of appearance was either predictable or unpredictable. With this design, we aimed 1) to investigate the dynamics of STN and STN-cortex beta synchronization during movement reversal compared to those of starting and stopping and 2) to assess the effect of the temporal predictability of movement instructions on the coordination of beta synchronization for starting, stopping, and reversing.

Results

Effect of predictability on behavior

Patients turned a wheel (Fig. 1b) with their index finger at their preferred speed and were prompted by visual cues to start, reverse, or stop rotational movement (Fig. 1a) while we simultaneously measured MEG and STN local field potentials (LFPs). We considered an average of 60.1 (SD = 14.8) predictable start trials, 59.3 (SD = 17.8) unpredictable start trials, 59.9 (SD = 14.2) predictable reversal trials, 58.6 (SD = 15.5) unpredictable reversal trials, 61.3 (SD = 15.3) predictable stop trials and 60.2 (SD = 17.4) unpredictable stop trials per patient for analysis. To assess whether the predictability of movement prompts had an effect on behavior, we analyzed its effect on movement speed and reaction times. Angular speed changes in start, reversal, and stop trials were similar in the predictable and the unpredictable condition when the data was aligned to action onset (Fig. 1c, Fcond(1,19) = 0.039, pcond = 0.846; see Supplementary File 1 for the complete results of the ANOVA). Aligning trials to cue onsets revealed that starting and stopping occurred slightly later in the unpredictable condition (Fig. 1d). This was reflected by a main effect of condition (Fcond(1,19) = 7.425, pcond = 0.013) and a condition*movement type interaction effect (Fcond*mov(2,18) = 5.347, pcond*mov = 0.015) on reaction times. Post-hoc t-tests revealed that reaction times to predictable start cues (M = .757, SD = 0.154) and stop cues (M = .824, SD = 0.202) were significantly shorter than reaction times to unpredictable start (M = .840, SD = 0.160) and stop (M = .889, SD = 0.233) cues (start: t = -3.469, one-sided p = 0.001, stop: t = -2.213, one-sided p = 0.020). Thus, starting and stopping were not performed at different speeds across conditions, but were initialized earlier in the predictable condition.

Paradigm and behavioral results.

(a) Patients were cued by arrows to start turning or reverse movement direction. Stop cues were presented at the end of each sequence. The timing of cues varied with the condition: in the predictable condition, the start cue was always followed by a reverse cue after 4 s and a stop cue after another 4 s (no jitter). In the unpredictable condition, there were either 0, 1, or 2 reversals (equal probability). Cue onset was jittered. CW: clockwise, CCW: counterclockwise. (b) Turning device for motor paradigm. (c) Average movement-aligned wheel speed. Red dotted lines indicate when turning began, was reversed in direction, and halted. (d) Average cue-aligned wheel speed. Red dotted lines indicate when the start, reversal and stop cue appeared, respectively. N = 20.

Modulations of spectral power associated with starting and stopping

In order to assess whether beta power modulations associated with reversals were distinct from beta suppression and rebound (research aim 1), we centered the trials on movement initiation, reversal, and termination, respectively, and assessed beta power dynamics. Besides the STN, this was done for two motor cortical regions of interest (ROIs): primary motor cortex (M1, hand knob region) and medial sensorimotor cortex (MSMC). This choice was based on the strongest movement-related modulations of beta power and coherence (Fig. 2, see Regions of interest in the Methods section for further detail). For comparison, we also present group average time-frequency spectra for the gamma frequency band.

Regions of interest.

(a, b) 3D source-reconstruction in MNI space. N = 20. White crosses mark the cortical ROIs selected for further analysis based on the strongest relative change in power (a) and the strongest absolute change in coherence (b). (c) All patients’ DBS electrodes, localized with Lead-DBS.

As expected, starting to turn the wheel was associated with a prominent beta suppression (contralateral STN: tclustersum = -2128.9, p < 0.001; ipsilateral STN: tclustersum = -2062.8, p < 0.001; contralateral M1: tclustersum = -8434.8, p < 0.001, ipsilateral M1: tclustersum = -8199.5, p < 0.001), whereas stopping resulted in a beta rebound (contralateral STN: tclustersum = 2843.0, p < 0.001; ipsilateral STN: tclustersum = 1488.8, p = 0.003; contralateral M1: tclustersum = 5958.9, p < 0.001, ipsilateral M1: tclustersum = 3834.7, p < 0.001) in both motor cortex and STN (Fig. 3a, Fig. 4a, b). The beta suppression occurred bilaterally while the beta rebound was more lateralized to the hemisphere contralateral to movement, as corroborated by a statistical analysis of the lateralization index (Fig. 4c and Supplementary File 2). Power changes in MSMC were generally similar to those in M1. Significant increases in gamma power at movement start were only observed in the contralateral STN (tclustersum = 2216.5, p = 0.003, Fig. 3a). At movement stop, there was a decrease in gamma power in contralateral STN (tclustersum = -734.8, p = 0.016, Fig. 3a) and contralateral M1 (tclustersum = -1447.4, p = 0.002, Fig. 4b). However, changes in gamma power were overall much smaller in magnitude compared to the beta suppression and rebound.

Movement-related beta power modulations in the STN.

(a) Time-frequency spectra of start, reversal, and stop trials for the STN (group average, trials averaged across predictability conditions). Time 0 marks the moment turning began, was reversed in direction, and halted (red lines). The black line in each plot represents the average wheel turning speed (scale: 0-600 deg/s). Power was baseline-corrected (baseline: -1.6-0 s). Hatched lines within black contours indicate significant changes relative to baseline. N = 20. (b) Six examples of individual patients at reversal and stop. Power was baseline-corrected (baseline: -1.6-0 s). Time 0 marks the brief pause of movement occurring during reversals, and movement stop, respectively (red lines). The black line in each plot represents each patient’s trial-average wheel turning speed (scale: 0-600 deg/s; for patient 21, the scale was adapted to 0-750 deg/s). Patient 10: contralateral, predictable; Patient 1: contralateral, unpredictable, Patient 2: contralateral, predictable; Patient 6: ipsilateral, unpredictable; Patient 21: contralateral, unpredictable; Patient 22; ipsilateral, unpredictable.

Cue-aligned beta power modulations in the STN.

Time-frequency spectra of cue-aligned start, reversal, and stop trials for the STN (group average, trials averaged across predictability conditions). Time 0 marks the appearance of the cue to start, reverse or stop turning (red lines). The black line in each plot represents the average wheel turning speed (scale: 0-600 deg/s). Power was baseline-corrected (baseline: -1.6-0 s). Hatched lines within black contours indicate significant changes relative to baseline. N = 20.

Movement-related beta power modulations in M1.

N = 20. (a) Source-localized movement-related modulation of beta power at movement start, reversal, and stop (MNI space, group average, trials averaged across predictability conditions). The hemisphere contralateral to movement is on the left. (b) Time-frequency spectra of start, reversal, and stop trials for M1. Time 0 marks the time point turning began, was reversed in direction, and halted (red lines). The black line in each plot represents the average wheel turning speed (scale: 0-600 deg/s). Power was baseline-corrected (baseline: -1.6-0 s). Hatched lines within black contours indicate significant changes relative to baseline. (c) Lateralization index for M1. LI = 0 corresponds to no lateralization; positive values refer to a contralateral lateralization and negative values to an ipsilateral lateralization. Blue: beta suppression; red: beta rebound.

Cue-aligned beta power modulations in M1.

Time-frequency spectra of cue-aligned start, reversal, and stop trials for M1. Time 0 marks the appearance of the cue to start, reverse or stop turning (red lines). The black line in each plot represents the average wheel turning speed (scale: 0-600 deg/s). Power was baseline-corrected (baseline: -1.6-0 s). Hatched lines within black contours indicate significant changes relative to baseline. N = 20.

Modulation of STN spectral power associated with reversals

Modulations of beta oscillations associated with reversals of movement direction were of particular interest to this study (research aim 1). When reversing, one first needs to stop the ongoing movement before accelerating again in the opposite direction. Stopping is followed by the beta rebound whereas starting is preceded by beta suppression. To the best of our knowledge, no study has investigated the neural signals underlying acceleration that immediately follows stopping. In the STN, reversals were associated with a brief modulation of beta power, which was weak when averaged across patients (Fig. 3a). Reversal-related beta power modulations of individual patients were variable. Some patients revealed brief increases whereas others showed decreases of STN beta power upon reaching the turning point (Fig. 3b). Reversal-related increases of beta power differed from the beta rebound, as occurring after termination of the movement sequence, with respect to amplitude and spectral content, often lacking the low-beta component of the beta rebound (Fig. 3b). These findings demonstrate distinct processing of brief pauses of action vs. a complete halt of action.

Modulation of cortical beta power associated with reversals

With respect to cortical beta power dynamics during reversals (research aim 1), we observed that reversals were associated with a brief suppression of alpha and beta power in motor cortex, particularly in M1 (contralateral M1: tclustersum = -1492.7, p < 0.001; ipsilateral M1: tclustersum = -3326.2, p < 0.001; Fig. 4a, b). The suppression occurred after the turning point had been reached (Fig. 4b) and was stronger in the hemisphere ipsilateral to movement, as demonstrated by a significant ROI*movement interaction effect on baseline-corrected beta power (FROI*movement(10,9) = 5.711, pROI*modulation = 0.008, refer to Supplementary File 3 for the full results of the ANOVA). Post-hoc t-tests confirmed that beta power was at a lower level in ipsilateral M1 (M = -0.47, SD = 0.33) compared to contralateral M1 (M = -0.21, SD = 0.23) during reversal of movement direction (t = 4.454, one-sided p < 0.001).

To test whether the ipsilateral lateralization was related to pre-event baseline (i.e., pre-reversal, pre-start, and pre-stop) levels, we re-computed the modulations using a whole recording average baseline (power averaged over all time points and movement types), thereby omitting the pre-vs. post-event contrast. Fig. 5a illustrates that movement-related power modulations were generally stronger in the hemisphere contralateral to movement, with the exception of acceleration, which was associated with bilateral suppression of beta power (compare the bilateral beta power suppression at post-start and post-reversal to the contralateral beta power modulations in all other plots). The second before reversing, beta-power was at an intermediate level in the hemisphere ipsilateral to movement (Fig. 5a). Thus, we observed a suppression relative to the pre-event baseline (Fig. 4a-b). In the contralateral hemisphere, in contrast, beta power could not be suppressed much further because it was already close to floor level prior to reversing (Fig. 5b). The lack of a pre-reversal increase of beta power is remarkable, because the second prior to reversing contained the deacceleration of the moving hand, which does not appear to involve an increase of beta power in primary motor cortex.

Pre-and post-event beta power.

N = 20. (a) Source-localized modulation of beta power before and after movement start, reversal, and stop (-1-0 and 0-1 s with respect to the movement of interest; baseline: power averaged over all time points and movement types). Plots are group-averages in MNI space, trials were averaged across predictability conditions. (b) Relative change with respect to whole recording average baseline, of ipsilateral and contralateral beta power for pre-reversal and post-reversal time windows.

Movement-related modulations of STN-cortex coherence

Beyond the local changes in beta power, we intended to investigate the dynamics of oscillatory coupling within the basal ganglia-cortex loop in the context of reversals of movement direction (research aim 1). The movement-related modulations of STN-cortex coherence were similar to modulations of power, including beta suppression (predictable start: tclustersum = -489.4, p = 0.002; unpredictable start: tclustersum = -530.1, p = 0.003), beta rebound (predictable stop: tclustersum = 802.0, p = 0.002; unpredictable stop: tclustersum = 1252.2, p < 0.001), and increases in the gamma band at movement start (predictable start: tclustersum = 120.4, p = 0.005; unpredictable start: tclustersum = 197.8, p < 0.001). Unlike motor cortical beta power, however, STN-cortex beta coherence did not decrease in the re-acceleration phase of reversals. On a qualitative level, it even increased relative to baseline. Fig. 6 depicts the movement-related changes in coherence averaged across all ROIs.

Event-related modulations of STN-cortex coherence and the effect of predictability.

N = 20 (a) Baseline-corrected group average time-frequency representations of STN-cortex coherence (averaged over ROIs) during start, reversal and stop for both the predictable and the unpredictable trials (baseline: -1.6-0 s). Time 0 marks the moment turning began, was reversed in direction, and was halted, respectively (red line). The black line in each plot represents the average wheel turning speed (scale: 0-600 deg/s). Hatched lines within black contours indicate significant changes relative to baseline. (b) Group average coherence difference between the unpredictable and predictable condition. Left: Contrast of time-frequency representations. TFRs were averaged over ROIs. Right: Contrast of source-localized, event-related coherence modulations in the beta band.

Cue-aligned modulations of STN-cortex coherence and the effect of predictability.

Baseline-corrected group average of time-frequency representations of STN-cortex coherence (averaged over ROIs) during start, reversal and stop (cue-aligned) for both the predictable and the unpredictable trials (baseline: -1.6-0 s). Time 0 marks the appearance of the cue to start, reverse or stop turning (red lines). The black line in each plot represents the average wheel turning speed (scale: 0-600 deg/s). Hatched lines within black contours indicate significant changes relative to baseline. N = 20.

Effects of predictability on beta power and STN-cortex coherence

With respect to the effect of predictability of movement instructions on beta dynamics (research aim 2), we observed that event-related modulations of beta coherence were generally stronger in the unpredictable than in the predictable condition (Fcond(1,18) = 7.210, pcond = 0.015; Supplementary File 3), suggesting that the inability to predict the next cue prompted intensified coupling between STN and cortex. Fig. 6 illustrates this general effect for coherence averaged across ROIs. When comparing region-average TFRs between the unpredictable and the predictable condition, we observed a significant difference for stopping in the beta band (tclustersum = 142.8, p = 0.023). The effect of predictability on beta coherence was most pronounced in the MSMC (Fig. 6b).

In case of power, we found an interaction between movement type and condition (Fcond*mov (2,17) = 4.394, pcond*mov = 0.029), such that beta suppression was generally stronger in the predictable (M = -.170, SD = 0.065) than in the unpredictable (M= -.154, SD = 0.070) condition across ROIs (t = -1.888, one-sided p = 0.037).

Discussion

Our study demonstrates that initiating, reversing, and stopping a continuous movement involves modulation of local and long-range beta synchronization in basal ganglia-cortex loops. Accelerating, stopping briefly, and coming to a complete halt have distinct and region-specific effects on beta oscillations in the motor system. These effects are context-dependent, with modulations of subcortico-cortical coupling intensifying when the upcoming movement instructions cannot be anticipated.

The dynamics of STN-cortex coherence

Simultaneous measurements of subthalamic and cortical oscillations in a comparably complex motor task allowed us to study the context-dependent dynamics of STN-cortex coupling. STN-M1 and STN-MSMC beta coherence decreased at movement initiation and increased after movement termination, while gamma coherence increased at movement start, corresponding to similar power changes in STN and motor cortex. A pre-movement suppression of beta coherence has been reported previously (Cassidy et al., 2002; Talakoub et al., 2016; van Wijk et al., 2017). Similarly, increases in gamma coherence have been found for the performance of ballistic movements (Alegre et al., 2013; Litvak et al., 2012). Stopping a planned movement has been found to be linked with reduced suppression of beta coherence (Alegre et al., 2013), but a post-movement increase of coherence has thus far only been described for ballistic movements (Tan et al., 2014). Considering the timing of the increase observed here, the STN’s role in movement inhibition (Benis et al., 2014; Ray et al., 2012) and the fact that frontal and prefrontal cortical areas are believed to drive subthalamic beta activity via the hyperdirect pathway (W. Chen et al., 2020; Oswal et al., 2021), it seems plausible that the increase of beta coherence reflects feedback of sensorimotor cortex to the STN in the course of post-movement processing.

STN-cortex coherence and beta suppression are modulated by predictability

In the present paradigm, patients were presented with cues that were either temporally predictable or unpredictable. We found that unpredictable movement prompts were associated with stronger modulations of STN-cortex beta coherence, likely reflecting the patients adopting a more cautious approach, paying greater attention to instructive cues. Enhanced STN-cortex synchronization might thus indicate the recruitment of additional neural resources, which might have allowed patients to maintain the same movement speed in both conditions.

The notion of beta oscillations reflecting motor and cognitive processes such as action selection, clearing, and error-monitoring has gained growing support (Fonken et al., 2016; Schmidt et al., 2019; Turner & Desmurget, 2010). Purely cognitive inhibition processes, such as inhibition of thoughts, have been found to be associated with prefrontal beta power modulations (Castiglione, Wagner, Anderson, & Aron, 2019; Schmidt et al., 2019). Further, the STN has been suggested to implement its hold your horses function, reflected by beta-band synchronization, in situations of cognitive conflict (Brittain et al., 2012). Simultaneous measurements of MEG and STN LFPs revealed that STN-cortex beta coherence increases after conflict cues in an expanded judgement task (Patai et al., 2022). Though the present paradigm did not involve any conflict as such, the context of unpredictable movement instructions possibly engaged similar cognitive processes in response to surprise/uncertainty.

With respect to power, we observed reduced beta suppression in the unpredictable condition at movement start, likely demonstrating a lower level of motor preparation. This finding aligns with MEG research that found reduced beta suppression with enhanced uncertainty in a motor task (Tzagarakis, Ince, Leuthold, & Pellizzer, 2010), and findings from an EEG study that demonstrated reduced beta suppression in response to an unpredictable sequence of rhythmic stimuli (Alegre et al., 2003). Although previous research has reported modulations of the beta rebound by cognitive factors (Fischer, Tan, Pogosyan, & Brown, 2016; Tan, Wade, & Brown, 2016; Zavala et al., 2018), we did not find an effect of predictability on the beta power rebound here.

Acceleration involves the recruitment of ipsilateral M1

As expected, we found sustained beta suppression at movement start and a strong beta power rebound at movement stop in STN, M1, and MSMC. During reversals, beta power was suppressed briefly in M1, particularly in the ipsilateral hemisphere where beta was not fully desynchronized prior to reversing. In contrast, the contralateral hemisphere revealed a floor effect: the ongoing movement resulted in more persistent beta power suppression that was only slightly intensified when reversing. Bilateral modulation of beta power, as reported during reversals, was otherwise observed during the initiation of movement, but not during ongoing movement or after movement termination, suggesting that the recruitment of ipsilateral M1 may be selective to acceleration.

Our findings are consistent with prior studies that have demonstrated a bilateral (Alegre et al., 2004; Zaepffel, Trachel, Kilavik, & Brochier, 2013) and spatially diffuse (Jurkiewicz et al., 2006) beta suppression, and more focal (Jurkiewicz et al., 2006) and predominantly contralateral (Espenhahn et al., 2017) topography of the beta rebound. Further, past research has posited a role of ipsilateral motor cortex in motor control and preparation (Jurkiewicz et al., 2006; Olson et al., 2022). Beta suppression in the ipsilateral hemisphere has been found to be related to increased corticospinal excitability, to facilitate finger movements (Rau, Plewnia, Hummel, & Gerloff, 2003), and to have a role in higher order cortical processing of fine motor programs (R. Chen, Gerloff, Hallett, & Cohen, 1997).

Brief pauses and complete stops have distinct effects on beta oscillations

We did not find evidence of a beta rebound following the short pause of movement during reversals in motor cortex. Instead, we observed a transient broadband beta power suppression in cortex, which was likely related to re-acceleration in the opposite direction. In contrast, the STN exhibited increases of high beta power in some patients, compatible with post-processing of the brief pause of movement occurring during reversals. On an observational level, the spectral patterns of these increases did not entirely match the individual stop-related beta pattern, lacking the low-beta component of the beta rebound. Thus, STN low-beta oscillations might not re-emerge when stopping briefly within a movement sequence, corroborating a dissociation of low-and high-beta oscillations, as proposed previously (Chandramouli, Iliana, & Krishna, 2019; Oswal et al., 2021; Patai et al., 2022). Given that the beta rebound has been reported to slow reaction times (Muralidharan & Aron, 2021) and to reduce corticospinal excitability (J. R. Wessel et al., 2016), and that beta power must decrease for movement to start (Heinrichs-Graham & Wilson, 2016; Khanna & Carmena, 2017), it is likely that at least the low-beta portion of the beta rebound needs to be avoided during changes of ongoing action because it would slow down re-acceleration otherwise.

In agreement with the current findings, previous research assessing STN-and cortical beta activity reported no beta rebound around the time a movement changed (Alegre et al., 2004; London et al., 2021; Oswal et al., 2016), except for one study, which did report a cortical beta rebound between two successive movements (Muralidharan & Aron, 2021). It should be noted, however, that the pauses were ∼1 to 2 s long. In our study, the beta rebound occurred only at the end of the movement sequence, when patients were already in the process of stopping and movement had already slowed. This picture emerged irrespective of whether power dynamics were analyzed in a movement-or cue-aligned fashion (see Figure 3-figure supplement 1, Figure 4-figure supplement 1 and Figure 6-figure supplement 1). A causal role of the beta rebound in stopping is therefore implausible. More likely, the rebound serves as a post-movement feedback signal reflecting task-dependent contextual information used to either confirm or update motor plans (Alegre et al., 2004; Cao & Hu, 2016). Alternatively, it might indicate the clear-out of the entire motor program (Schmidt et al., 2019).

With respect to gamma activity, we observed increases in power at movement start in the contralateral STN and decreases in power at movement termination in contralateral STN and M1. While movement-related increases in gamma power are an established finding in the literature (Litvak et al., 2012; Lofredi et al., 2018), there appears to be no consensus on its functional role during movement stopping. Previous studies using auditory stop signals reported STN gamma power increases in response to stop signals (Fischer et al., 2017; Ray et al., 2012). When assessed within a brief critical window between the stop signal and the average time of the upcoming finger tap, gamma power even correlated with stopping success, i.e. gamma was stronger when the downward movement was stopped earlier (Fischer et al., 2017). Conversely, another study using visual stop signals reported decreased STN gamma power (Alegre et al., 2013). We are unaware of studies that have assessed gamma power changes when stopping a continuous movement in response to visual cues and therefore provide first evidence for a decrease in this scenario, although we cannot rule out that focusing on different DBS contacts, or using auditory stop signals and shorter event-locked analyses windows might produce different results.

Limitations and future directions

In humans, invasive measurements of neuronal activity are only possible in patients undergoing surgery. While patients likely presented with pathologically changed oscillatory activity, we suggest that our approach at least approximates healthy brain functioning as patients were on their usual dopaminergic medication. Dopaminergic medication has been demonstrated to normalize power within the STN and globus pallidus internus, as well as STN-globus pallidus internus and STN-cortex coherence (Brown et al., 2001; Hirschmann et al., 2013). Transferring our approach to patients with different disorders, e.g. dystonia or essential tremor, or examining healthy participants solely at the cortical level, could contribute to elucidating whether the synchronization dynamics reported here are indeed independent of PD.

Due to the diversity of modulations across patients, we cannot provide a general description of how the STN responds to reversals. The diversity may result from the fact that the exact recording site varied across patients, although all recording contacts were located in the dorsolateral STN. Further, stop-processes, mediated by the hyperdirect and the indirect pathway as well as cortico-striatal go-processes may emerge in the basal ganglia close in time (Muralidharan, Aron, & Schmidt, 2022; Schmidt & Berke, 2017), potentially overlapping. The sub-populations processing these signals in the STN (Isoda & Hikosaka, 2008; Schmidt & Berke, 2017; Schmidt, Leventhal, Mallet, Chen, & Berke, 2013) might not be resolvable with macro-electrode LFP recordings.

Conclusion

In conclusion, we have revealed distinct local and long-range synchronization dynamics of motor cortex and STN during changes of ongoing action in different movement contexts. We found that stopping briefly in the course of changing movement direction and terminating a movement sequence have distinct oscillatory profiles. Moreover, movement scenarios that do not permit movement preparation and require higher levels of caution appear to involve enhanced levels of subthalamo-cortical beta synchronization, highlighting that long-range beta coherence plays an important role in coordinating movements in response to unpredictable events.

Materials and methods

Patients

23 PD patients with a mean age of 66.13 years (± 7.72 years) participated in the study (Table 1). DBS surgery was performed by the department of Functional Neurosurgery and Stereotaxy of the University Hospital Düsseldorf under full anesthesia and according to standard procedures. 21 patients were implanted with Abbott Infinity segmented leads (Abbott Laboratories, Chicago, Illinois, USA) and three patients with Medtronic SenSight electrodes (Medtronic Inc., Minneapolis, MN, USA). DBS surgery was performed in two steps, and the measurements took place in between the implantation of the electrodes and the implantation of the pulse generator. Prior to participating, all patients provided their written informed consent in agreement with the declaration of Helsinki. The study was approved by the Ethics Committee of the Medical Faculty of Heinrich Heine University Düsseldorf. The medication schedule was not changed for this experiment (Med ON state). Three participants were excluded from the analyses, two of whom were physically unable to perform the paradigm. The data of the third patient were contaminated by excessive artifacts.

Patient clinical characteristics.

Disease duration refers to the time since diagnosis. For patient 4, the time since first symptom manifestation is given. MoCa = Montreal Cognitive Assessment Test.

Recordings

Measurements took place the day after the implantation of DBS electrodes. Externalization of leads allowed us to measure LFPs from the STN in combination with MEG. For LFP recordings, we used a mastoid reference and re-referenced the signals using a bipolar montage post-measurement. MEG signals were acquired simultaneously, using a 306-channel whole-head MEG system (VectorView, MEGIN). Muscle and ocular activity were monitored via electromyography (EMG) and vertical and horizontal electrooculography (EOG), respectively. EMG surface electrodes were placed on patients’ right and left forearms, referenced to the muscle tendons at the wrist. We first recorded 5 minutes of resting-state data, followed by the motor task, which lasted for about 32 minutes in total. During the task, patients were required to turn a wheel clockwise or counterclockwise, according to visual instructions presented on a screen in front of them.

Experimental design

Patients were seated in the MEG scanner in a magnetically shielded room with a turning device (‘wheel’, Fig. 2b) placed on a table in front of them. The wheel (diameter = 14 cm, height = 6.5 cm) could be turned into both directions and had indentations, allowing comfortable placement of one index finger for turning. An MEG-compatible plastic fiber optic position sensor system (MR430 Series ZapFREE Fiber Optic Absolute Encoder System, MICRONER INC., Camarillo, CA, USA) was used to measure wheel turning. The absolute angular position was continuously measured and updated at a frequency of 1.2 kHz. Given its design, the sensing system did not introduce any magnetic interference.

Movement prompts were presented on a screen. The visual stimuli consisted of two curved arrows pointing either clockwise or counterclockwise, respectively, and a stop sign with white font on a red background (Fig. 1a). Patients were instructed to turn the wheel with their index finger following the direction of the arrows and to stop when a stop cue appeared. We did not impose requirements on turning speed or body side, so that patients could use their less affected hand and adjust the speed to their individual motor capabilities.

The experiment was conducted in two distinct blocks, where stimuli differed in their order and timing of presentation. In the predictable condition, trials consisted of a blue arrow cueing the patients to start turning clockwise, followed by a cue to change the turning direction after 4 s, and a stop cue after another 4 s. Each trial was followed by a pause lasting for 4 s. The condition was termed predictable, as the fixed timing and the fixed order of cues allowed patients to easily predict and prepare what they had to do next and when. In the unpredictable condition, the start cue was either clockwise or counterclockwise and was followed by 0, 1, or 2 reversals before the stop cue appeared. Each alternative occurred equally often (33%: go, stop; 33%: go, reverse, stop; 33%: go, reverse, reverse, stop; clockwise and counterclockwise start directions were balanced). Additionally, the intervals between the visual stimuli were unpredictable (ranging between 4-7 s), with 50% of all inter-stimulus intervals kept at 4 s, as in the predictable condition. Hence, patients could not foresee the sequence and the timing of instructions, calling for a more cautious/attentive monitoring of cues.

We recorded two blocks per condition, with 36 trials each, in a pseudo-randomized fashion. To enhance compliance, we split each block in half, allowing for a short break, and also offered breaks between blocks.

Data analysis

Data were analyzed using MATLAB R2019b (The Mathworks, Natick, Massachusetts, USA) and the toolbox FieldTrip (Oostenveld, Fries, Maris, & Schoffelen, 2011). For statistical testing, we used IBM SPSS Statistics 28 (IBM Corporation, Somers, USA).

Preprocessing

The data were visually inspected to identify and tag noisy channels and subsequently cleaned using temporal Signal Space Separation to remove artifacts originating from outside the MEG sensor array (Taulu & Simola, 2006). Then, the data were downsampled to 500 Hz. We applied a high-pass finite impulse response filter with a cut-off frequency of 1 Hz to remove low-frequency drifts and screened the data for remaining artifacts.

We used custom MATLAB-scripts for semi-automated detection of movement start, reversal and stop in the wheel data. This was achieved by applying an event-specific combination of amplitude and duration thresholds to the first temporal derivative of the rotation angle measurements of the wheel. To ensure that events were correctly marked, all events were visually inspected and manually corrected if needed. Then, we epoched MEG and LFP data with respect to the behavioral events. Trials were centered around movement events of interest, i.e., start, stop, and reversal of movement and encompassed 4 s. Movement-aligned angular speed was calculated within those time-windows and averaged over trials. Reaction times to cues were defined as the time from cue presentation until movement.

LFP channel selection

The positions of DBS electrodes were localized with the MATLAB toolbox Lead-DBS (Neudorfer et al., 2023) using the patients’ pre-operative T1-and T2-weighted MRIs (Magnetom Trio MRI scanner, Siemens, Erlangen, Germany) and postoperative CT scans (Fig. 1b). In order to select one LFP channel for each patient and hemisphere, we identified the bipolar LFP channel with the strongest beta suppression and beta rebound, as previous research has demonstrated the presence of these modulations in the dorsolateral motor STN (Benis et al., 2014; J. R. Wessel et al., 2016). Moreover, the source of subthalamic beta oscillations has been localized to the dorsal STN (Tamir et al., 2020).

Regions of interest

Similarly, we selected cortical regions of interest (ROIs) by localizing the strongest event-related modulations of beta power/beta coherence. For source localization, we first co-registered the pre-operative T1-weighted MRI scans to the MEG coordinate system. Using the segmented MRIs, we prepared forward models based on single-shell realistic head models (Nolte, 2003). Beamformer grids, specifying the position of sources, covered the entire brain and were aligned to Montreal Neurological Institute (MNI) space. Subsequently, we applied Dynamic Imaging of Coherent Sources (DICS) (Gross et al., 2001) to beta-band LFP-MEG cross-spectral densities pooled across predictability conditions. Next, we computed contrasts between post-event (0 to 2 s) and pre-event (-2 to 0 s) beta power/coherence and averaged the absolute changes across patients and events (moment start, reversal, and stop). This served to identify the regions with the strongest change in general, irrespective of sign, event, and predictability condition.

The strongest beta power modulations localized to the hand knob area of primary motor cortex (M1; Fig. 1a) and the strongest changes in STN-cortex beta coherence to medial sensorimotor cortex (MSMC; Fig. 1b). Thus, we focused our analysis on bilateral STN, M1 and MSMC. For time-frequency analysis, we represented each cortical ROI by the grid point of strongest modulation and its six nearest neighbors and extracted a time-series for each grid point, using a linearly constrained minimum variance spatial filter (Van Veen & Buckley, 1988).

Time-frequency analyses

While our main focus was the beta band, we also included other frequencies in our time-frequency analyses to get a more complete picture of power coherence changes. Specifically, we considered the frequency ranges 5 to 45 Hz and 55 to 90 Hz, omitting the 50 Hz line noise artifact, and the time range from -1.6 to 1.6 s with respect to the movement event. Fourier spectra were computed using a multi-taper approach (4 Slepian tapers for the low frequency range and 7 Slepian tapers for the high frequency range), a window size of 800 ms and a step size of 50 ms. Using the Fourier coefficients, we computed power and STN-cortex across-trial coherence for each time-frequency bin.

For illustration, we applied baseline-correction, using the mean of the pre-event time window (-1.6 to 0 s) as baseline. In case of power, we expressed changes with respect to baseline in decibel. In case of coherence, we subtracted the baseline values. Time-frequency spectra of cortical sources were averaged over neighboring grid points belonging to the same ROI.

Statistical Analysis

Repeated measures analyses of (co)variance (rmANCOVA), implemented in SPSS, were our main tool for statistical analysis. This approach provides a comprehensive, multi-factorial analysis, but requires pre-selection of brain areas (see Regions of Interest), a frequency range and a time range of interest (see Power and coherence). The dependent variable was either the event-related change in power or coherence or the hemispheric lateralization of the event-related power change (see Lateralization). The main independent variable of interest was predictability. Brain area and movement type were also included as factors due to their clear effects on power and coherence, but their main effects are not reported in the main paper. They can be found in the Supplementary Material. Because Mauchly’s test indicated violations of the sphericity assumption we report results from the multivariate test (Rasch, Friese, Hofmann, & Neumann, 2021).

The rmANCOVAs were complemented by cluster-based permutation tests for detecting significant power/coherence modulations relative to baseline. These tests are mono-factorial, but have the advantage of not requiring any preselection of time or frequency ranges while providing correction for multiple comparisons. The cluster-defining and the cluster significance threshold was set to 0.05 (two-sided test). The cluster statistic was the sum of t-values within a cluster. We performed 1000 permutations per test.

Power and coherence

To assess the effect of predictability on power, we conducted a repeated measures ANCOVA testing the influence of the factors movement (start, stop, reversal), predictability (predictable, unpredictable) and brain area (STN, M1, MSMC, ipsilateral and contralateral to the moving hand), as well as interactions between these factors, on the event-related modulation of beta power. Here, modulation refers to the difference between post-event (0 -1.6 s) and pre-event (-1.6 - 0) beta power in decibel (dB). A similar rmANCOVA was computed for event-related modulations of STN-cortex beta coherence. In this case, we considered the difference between pre- and post-event coherence. Here, the factor brain area contained of the following pairs: contralateral M1-contralateral STN, contralateral MSMC-contralateral STN, ipsilateral M1-ipsilateral STN, ipsilateral MSMC-ipsilateral STN, with the terms ipsilateral and contralateral referring to the moving hand. Because beta power is known to correlate with movement speed (Lisi & Morimoto, 2015; Lofredi et al., 2023; Pogosyan, Gaynor, Eusebio, & Brown, 2009), we added turning speed, averaged over trials and timepoints, as a covariate.

Lateralization

We compared the beta power suppression and the beta power rebound with respect to their hemispheric lateralization, using a rmANOVA with the factors brain area (STN, M1, MSMC), predictability (predictable, unpredictable) and modulation type (beta suppression, beta rebound). Lateralization was quantified by the lateralization index, defined as the difference between contralateral and ipsilateral power, normalized by power summed over both hemispheres.

Behavior

To assess whether the predictability of movement prompts had an effect on the patients’ performance in the task, we performed a rmANOVA with the factors predictability (predictable, unpredictable), movement (start, stop, reversal) and their interaction on reaction times and movement-aligned wheel turning speed, averaged over trials and time points, respectively. Epochs without movement (pre-start and post-stop) were disregarded in this analysis.

Acknowledgements

This research was funded by the Brunhilde Moll Stiftung.

The authors thank all participants for their time, cooperation, and willingness to participate. Furthermore, the authors thank Hannah Feldmann for her contributions to developing the paradigm, Dafina Sylaj for her help with patient recruitment, and Lilli Ahrenberg for her work localizing the patients’ DBS electrodes.

Behavioral effects.

(A) Effects of condition (predictable, unpredictable) and movement (start, reverse, stop) on movement-aligned speed. (B) Effects of condition (predictable, unpredictable) and movement (start, reverse, stop) on reaction times to cues.

Effects on lateralization.

Effects of modulation type (beta suppression, beta rebound), condition (predictable, unpredictable) and ROI (STN, M1, MSMC) on lateralization index.

Effects on power and coherence.

(A) Effects of condition (predictable, unpredictable), movement (start, reverse, stop) and ROI (contralateral and ipsilateral STN, M1, MSMC) on normalized power, controlling for movement speed (cov). (B) Effects of condition (predictable, unpredictable), movement (start, reverse, stop) and ROI (contralateral STN-M1, contralateral STN-MSMC, ipsilateral STN-M1, ipsilateral STN-MSMC) on coherence modulation, controlling for movement speed (cov).