Abstract
The role of beta band activity in cortico-basal ganglia interactions during motor control has been studied extensively in resting-state and for simple movements, such as button pressing. However, little is known about how beta oscillations change and interact in more complex situations involving rapid changes of movement in various contexts.
To close this knowledge gap, we combined magnetoencephalography (MEG) and local field potential recordings from the subthalamic nucleus (STN) in Parkinson’s disease patients to study beta dynamics during initiation, stopping, and rapid reversal of rotational movements. The action prompts were manipulated to be predictable vs. unpredictable.
We observed movement-related beta suppression at motor sequence start, and a beta rebound after motor sequence stop in STN power, motor cortical power, and STN-cortex coherence. Despite involving a brief stop of movement, no clear rebound was observed during reversals of turning direction. On the cortical level, beta power decreased bilaterally following reversals, but more so in the hemisphere ipsilateral to movement, due to a floor effect on the contralateral side. In the STN, power modulations varied across patients, with patients revealing brief increases or decreases of high-beta power. Importantly, cue predictability affected these modulations. Event-related changes of STN-cortex beta coherence were generally stronger in the unpredictable than in the predictable condition.
In summary, this study reveals the influence of movement context on beta oscillations in basal ganglia-cortex loops when humans change ongoing movements according to external cues. We find that movement scenarios requiring higher levels of caution involve enhanced modulations of subthalamo-cortical beta synchronization. Further, our results confirm that beta oscillations reflect the start and end of motor sequences better than movement changes within a sequence.
Significance Statement
Beta synchrony between motor cortex and the subthalamic nucleus is intensified when instructional cues within a continuous motor sequence become less predictable, calling for more cautious behavior.
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.
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.
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.
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.
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.
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.
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.
References
- Alpha and beta oscillatory changes during stimulus-induced movement paradigms: Effect of stimulus predictabilityNeuroreport 14:381–385https://doi.org/10.1097/00001756-200303030-00017
- Alpha and beta oscillatory activity during a sequence of two movementsClinical Neurophysiology 115:124–130https://doi.org/10.1016/S1388-2457(03)00311-0
- The subthalamic nucleus is involved in successful inhibition in the stop-signal task: A local field potential study in Parkinson’s diseaseExperimental Neurology 239:1–12https://doi.org/10.1016/j.expneurol.2012.08.027
- Cortical and subcortical contributions to stop signal response inhibition: Role of the subthalamic nucleusThe Journal of Neuroscience 26:2424–2433https://doi.org/10.1523/jneurosci.4682-05.2006
- Inhibitory control and error monitoring by human subthalamic neuronsTranslational Psychiatry 4https://doi.org/10.1038/tp.2014.73
- Subthalamic nucleus activity dissociates proactive and reactive inhibition in patients with Parkinson’s diseaseNeuroimage 91:273–281https://doi.org/10.1016/j.neuroimage.2013.10.070
- The subthalamic nucleus: Unravelling new roles and mechanisms in the control of actionThe Neuroscientist 25:48–64https://doi.org/10.1177/1073858418763594
- A role for the subthalamic nucleus in response inhibition during conflictThe Journal of Neuroscience 32:13396–13401https://doi.org/10.1523/jneurosci.2259-12.2012
- Dopamine dependency of oscillations between subthalamic nucleus and pallidum in Parkinson’s diseaseThe Journal of Neuroscience 21:1033–1038https://doi.org/10.1523/jneurosci.21-03-01033.2001
- Beta rebound in visuomotor adaptation: Still the status quo?The Journal of Neuroscience 36:6365–6367https://doi.org/10.1523/JNEUROSCI.1007-16.2016
- Movement-related changes in synchronization in the human basal gangliaBrain 125:1235–1246https://doi.org/10.1093/brain/awf135
- Preventing a thought from coming to mind elicits increased right frontal beta just as stopping action doesCerebral Cortex 29:2160–2172https://doi.org/10.1093/cercor/bhz017
- Frequency shifts and depth dependence of premotor beta band activity during perceptual cecision-makingThe Journal of Neuroscience 39:1420–1435https://doi.org/10.1523/JNEUROSCI.1066-18.2018
- Involvement of the ipsilateral motor cortex in finger movements of different complexitiesAnnals of Neurology 41:247–254https://doi.org/10.1002/ana.410410216
- Prefrontal-subthalamic hyperdirect pathway modulates movement inhibition in humansNeuron 106:579–588https://doi.org/10.1016/j.neuron.2020.02.012
- Beta-band oscillation -signalling the status quo?Current Opinion in Neurobiology 20:156–165https://doi.org/10.1016/j.conb.2010.02.015
- Movement-related beta oscillations show high intra-individual reliabilityNeuroimage 147:175–185https://doi.org/10.1016/j.neuroimage.2016.12.025
- Beta synchrony in the cortico-basal ganglia network during regulation of force control on and off dopamineNeurobiology of Disease 127:253–263https://doi.org/10.1016/j.nbd.2019.03.004
- Subthalamic nucleus gamma activity increases not only during movement but also during movement inhibitioneLife 6https://doi.org/10.7554/eLife.23947
- High post-movement parietal low-beta power during rhythmic tapping facilitates performance in a stop taskEuropean Journal of Neuroscience 44:2202–2213https://doi.org/10.1111/ejn.13328
- Frontal and motor cortex contributions to response inhibition: Evidence from electrocorticographyJournal of Neurophysiology 115:2224–2236https://doi.org/10.1152/jn.00708.2015
- Dynamic imaging of coherent sources: Studying neural interactions in the human brainProceedings of the National Academy of Sciences 98:694–699https://doi.org/10.1073/pnas.98.2.694
- Is an absolute level of cortical beta suppression required for proper movement? Magnetoencephalographic evidence from healthy agingNeuroimage 134:514–521https://doi.org/10.1016/j.neuroimage.2016.04.032
- Differential modulation of STN-cortical and cortico-muscular coherence by movement and levodopa in Parkinson’s diseaseNeuroimage 68:203–213https://doi.org/10.1016/j.neuroimage.2012.11.036
- Role for subthalamic nucleus neurons in switching from automatic to controlled eye movementThe Journal of Neuroscience 28:7209–7218https://doi.org/10.1523/JNEUROSCI.0487-08.2008
- Post-movement beta rebound is generated in motor cortex: Evidence from neuromagnetic recordingsNeuroimage 32:1281–1289https://doi.org/10.1016/j.neuroimage.2006.06.005
- Beta band oscillations in motor cortex reflect neural population signals that delay movement onseteLife 6https://doi.org/10.7554/eLife.24573
- EEG single-trial detection of gait speed changes during treadmill walkPLOS ONE 10https://doi.org/10.1371/journal.pone.0125479
- Movement-Related Changes in Local and Long-Range Synchronization in Parkinson’s Disease Revealed by Simultaneous Magnetoencephalography and Intracranial RecordingsThe Journal of Neuroscience 32https://doi.org/10.1523/JNEUROSCI.0767-12.2012
- Subthalamic stimulation impairs stopping of ongoing movementsBrain 144:44–52https://doi.org/10.1093/brain/awaa341
- Dopamine-dependent scaling of subthalamic gamma bursts with movement velocity in patients with Parkinson’s diseaseeLife 7https://doi.org/10.7554/eLife.31895
- Pallidal beta activity is linked to stimulation-induced slowness in dystoniaMovement Disorders 38:894–899https://doi.org/10.1002/mds.29347
- Distinct population code for movement kinematics and changes of ongoing movements in human subthalamic nucleuseLife 10https://doi.org/10.7554/eLife.64893
- Distinct roles of dorsal and ventral subthalamic neurons in action selection and cancellationNeuron 109:869–881https://doi.org/10.1016/j.neuron.2020.12.025
- Behavioral induction of a high beta state in sensorimotor cortex leads to movement slowingJournal of Cognitive Neuroscience 33:1311–1328https://doi.org/10.1162/jocn_a_01717
- Transient beta modulates decision thresholds during human action-stoppingNeuroimage 254https://doi.org/10.1016/j.neuroimage.2022.119145
- Lead-DBS v3.0: Mapping deep brain stimulation effects to local anatomy and global networksNeuroimage 268https://doi.org/10.1016/j.neuroimage.2023.119862
- The magnetic lead field theorem in the quasi-static approximation and its use for magnetoencephalography forward calculation in realistic volume conductorsPhysics in Medicine and Biology 48:3637–3652https://doi.org/10.1088/0031-9155/48/22/002
- Cortical and Subthalamic Nucleus Spectral Changes During Limb Movements in Parkinson’s Disease Patients with and Without DystoniaMovement Disorders 37:1683–1692https://doi.org/10.1002/mds.29057
- FieldTrip: Open source software for advanced analysis of MEG, EEG, and invasive electrophysiological dataComputational Intelligence and Neuroscience 2011https://doi.org/10.1155/2011/156869
- Deep brain stimulation modulates synchrony within spatially and spectrally distinct resting state networks in Parkinson’s diseaseBrain 139:1482–1496https://doi.org/10.1093/brain/aww048
- Neural signatures of hyperdirect pathway activity in Parkinson’s diseaseNature Communications 12https://doi.org/10.1038/s41467-021-25366-0
- Cognitive factors modulate activity within the human subthalamic nucleus during voluntary movement in Parkinson’s diseaseThe Journal of Neuroscience 33:15815–15826https://doi.org/10.1523/JNEUROSCI.1790-13.2013
- Beta reactivity, prospective facilitation of executive processing, and its dependence on dopaminergic therapy in Parkinson’s diseaseThe Journal of Neuroscience 32:9909–9916https://doi.org/10.1523/JNEUROSCI.0275-12.2012
- Conflict detection in a sequential decision task is associated with increased cortico-subthalamic coherence and prolonged subthalamic oscillatory response in the β bandThe Journal of Neuroscience 42https://doi.org/10.1523/JNEUROSCI.0572-21.2022
- Beta rebound after different types of motor imagery in manNeuroscience Letters 378:156–159https://doi.org/10.1016/j.neulet.2004.12.034
- Boosting cortical activity at beta-band frequencies slows movement in humansCurrent Biology 19:1637–1641https://doi.org/10.1016/j.cub.2009.07.074
- Quantitative Methoden (5 ed. Vol. 2)Heidelberg: Springer
- Event-related desynchronization and excitability of the ipsilateral motor cortex during simple self-paced finger movementsClinical Neurophysiology 114:1819–1826https://doi.org/10.1016/S1388-2457(03)00174-3
- The role of the subthalamic nucleus in response inhibition: Evidence from local field potential recordings in the human subthalamic nucleusNeuroimage 60:271–278https://doi.org/10.1016/j.neuroimage.2011.12.035
- Functional segregation of movement-related rhythmic activity in the human brainNeuroimage 2:237–243https://doi.org/10.1006/nimg.1995.1031
- A pause-then-cancel model of stopping: Evidence from basal ganglia neurophysiologyPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 372https://doi.org/10.1098/rstb.2016.0202
- Beta oscillations in working memory, executive control of movement and thought, and sensorimotor functionThe Journal of Neuroscience 39:8231–8238https://doi.org/10.1523/JNEUROSCI.1163-19.2019
- Canceling actions involves a race between basal ganglia pathwaysNat Neurosci 16:1118–1124https://doi.org/10.1038/nn.3456
- Brain oscillatory correlates of visuomotor adaptive learningNeuroimage 245https://doi.org/10.1016/j.neuroimage.2021.118645
- Intracranial EEG reveals a time-and frequency-specific role for the right inferior frontal gyrus and primary motor cortex in stopping initiated responsesThe Journal of Neuroscience 29https://doi.org/10.1523/JNEUROSCI.3359-09.2009
- Time-course of coherence in the human basal ganglia during voluntary movementsScientific Reports 6https://doi.org/10.1038/srep34930
- Eight cylindrical contact lead recordings in the subthalamic region localize beta oscillations source to the dorsal STNNeurobiology of Disease 146https://doi.org/10.1016/j.nbd.2020.105090
- Post-movement beta activity in sensorimotor cortex indexes confidence in the estimations from internal modelsThe Journal of Neuroscience 36https://doi.org/10.1523/JNEUROSCI.3204-15.2016
- Human subthalamic nucleus in movement error detection and its evaluation during visuomotor adaptationThe Journal of Neuroscience 34https://doi.org/10.1523/JNEUROSCI.3414-14.2014
- Spatiotemporal signal space separation method for rejecting nearby interference in MEG measurementsPhysics in Medicine and Biology 51:1759–1768https://doi.org/10.1088/0031-9155/51/7/008
- Beta burst dynamics in Parkinson’s disease OFF and ON dopaminergic medicationBrain : a journal of neurology 140:2968–2981https://doi.org/10.1093/brain/awx252
- Cortical correlates of response time slowing in older adults: ERP and ERD/ERS analyses during passive ankle movementClinical Neurophysiology 127:655–663https://doi.org/10.1016/j.clinph.2015.05.003
- Basal ganglia contributions to motor control: A vigorous tutorCurrent Opinion in Neurobiology 20:704–716https://doi.org/10.1016/j.conb.2010.08.022
- Beta-band activity during motor planning reflects response uncertaintyThe Journal of Neuroscience 30:11270–11277https://doi.org/10.1523/jneurosci.6026-09.2010
- Beamforming: A versatile approach to spatial filteringIEEE ASSP Magazine 5:4–24https://doi.org/10.1109/53.665
- Low-beta cortico-pallidal coherence decreases during movement and correlates with overall reaction timeNeuroimage 159:1–8https://doi.org/10.1016/j.neuroimage.2017.07.024
- Establishing a right frontal beta signature for stopping action in scalp eeg: Implications for testing inhibitory control in other task contextsJournal of Cognitive Neuroscience 30:107–118https://doi.org/10.1162/jocn_a_01183
- β-Bursts Reveal the Trial-to-Trial Dynamics of Movement Initiation and CancellationThe Journal of Neuroscience 40https://doi.org/10.1523/JNEUROSCI.1887-19.2019
- Stop-related subthalamic beta activity indexes global motor suppression in Parkinson’s diseaseMovement Disorders 31:1846–1853https://doi.org/10.1002/mds.26732
- Non-selective inhibition of inappropriate motor-tendencies during response-conflict by a fronto-subthalamic mechanismeLife 8https://doi.org/10.7554/eLife.42959
- Modulations of EEG beta power during planning and execution of grasping movementsPLOS ONE 8https://doi.org/10.1371/journal.pone.0060060
- Cognitive control involves theta power within trials and beta power across trials in the prefrontal-subthalamic networkBrain 141:3361–3376https://doi.org/10.1093/brain/awy266
- Human subthalamic nucleus activity during non-motor decision makingeLife 6https://doi.org/10.7554/eLife.31007
Article and author information
Author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
Copyright
© 2024, Winkler et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.