Abstract
Beta event-related spectral perturbation, including bilateral movement-related beta desynchronization (MRBD) and postmovement beta synchronization (PMBS), can be evoked by unilateral speed movement. A potential correlation might exist between power (de)synchronization and interhemispheric coherence during movement execution. However, during the PMBS phase, the existence of interhemispheric coupling and the effect of speed on it are largely undiscovered. To answer this question, we investigated eight healthy, right-handed volunteers using a combination of electroencephalography, transcranial magnetic stimulation, and electromyography. We explored interhemispheric (directed) coherence during isotonic right index finger abduction movements at two speeds: ballistic and self-paced. We discovered that (1) interhemispheric coherence was greater during the PMBS than during the MRBD period. Furthermore, ballistic movement induced a larger coherence during the PMBS period, but not during the MRBD period. (2) In the MRBD phase, directed coherence from the contralateral motor cortex (CM1) to the ipsilateral motor cortex (IM1) was larger, with a reverse tendency observed during the PMBS period. Additionally, in ballistic movement, directed coherence from IM1 to CM1 was stronger and positively correlated with coherence, with no effect of speed on directed coherence detected in the MRBD phase. To explore the causality of interhemispheric coherence during the PMBS period, we investigated the interhemispheric inhibition (IHI) from IM1 to CM1 at different speeds. A stronger IHI from IM1 to CM1 at PMBS peak time was demonstrated, which was enhanced during ballistic movement. Additionally, IHI was negatively correlated with PMBS, and movement speed was positively associated with interhemispheric coupling during the PMBS period and IHI from IM1 to CM1.
Significance Statement
The present study explored interhemispheric (directed) coherence during isotonic right index finger abduction movements at two speeds: ballistic and self-paced. We discovered a dominance of interhemispheric coherence during the postmovement beta synchronization (PMBS) period of ballistic movement. Furthermore, directed coherence from the contralateral motor cortex (CM1) to the ipsilateral motor cortex (IM1) was more predominant in the movement-related beta desynchronization phase, with a reverse tendency observed during the PMBS period. Advanced exploration revealed a stronger interhemispheric inhibition (IHI) from IM1 to CM1 at PMBS peak time, which was enhanced during ballistic movement. IHI was negatively correlated with PMBS, and movement speed was positively associated with interhemispheric coupling during the PMBS period and IHI.
Introduction
Unilateral movement can induce alternation of beta rhythms in bilateral hemispheres (Pfurtscheller and Da Silva, 1999; Qiu et al., 2016). Reflected by electroencephalography (EEG) signals, this phenomenon is represented by a beta event-related spectral perturbation (ERSP; Pfurtscheller and Da Silva, 1999; Qiu et al., 2016). Beta ERSP encompasses a decreased beta power change relative to the baseline, namely, movement-related beta desynchronization (MRBD), followed by an increase in beta power after the termination of movement, known as postmovement beta synchronization (PMBS; Pfurtscheller and Da Silva, 1999; Deiber et al., 2015; Qiu et al., 2016). Additionally, interhemispheric beta coherence between motor areas arises during both unimanual and bimanual rhythmic movements (Mima et al., 2000; Gross et al., 2005), serving as evidence of the functional connections between brain regions (Sklar et al., 1972; Busk and Galbraith, 1975; Deiber et al., 2015).
The association between coherence and MRBD and PMBS over different time courses has been investigated in previous studies. Leocani et al. (1997) reported an increase in coherence over sensorimotor and frontal areas during the MRBD phase. Manganotti et al. (1998) revealed active intercommunication among bilateral, mesial central, and prefrontal regions during the MRBD period. The power increase in coupling during the MRBD period indicates that the cortex is effectively activated and involved in motor control (Leocani et al., 1997; Manganotti et al., 1998). This coupling became more intense as movement complexity increased (Manganotti et al., 1998). However, prior research also indicated inconsistent results regarding the effect of the movement rate on interhemispheric coherence (Serrien and Brown, 2002; Toma et al., 2002). For example, Toma et al. (2002) suggested that motor cortical activation and coupling were greater for faster movements, while Serrien and Brown (2002) found that coupling between the primary sensorimotor cortices in the beta frequency band was reduced with increasing movement speed, and this effect was more pronounced in the antiphase than the in-phase mode. As such, understanding interhemispheric coupling under different velocity conditions during the MRBD period of movement is crucial for comprehending the physiological phenomena and functional significance of the bilateral motor cortex.
Notably, existing studies indicate that beta rhythms over bilateral sensorimotor hand areas or other frontal areas demonstrate a complete lack of bilateral coherence during the PMBS period (Leocani et al., 1997; Gross et al., 2005). Similarly, the bilateral Rolandic mu rhythms show no coupling between the hemispheres during PMBS (Van Leeuwen et al., 1978; Schoppenhorst et al., 1980; Andrew and Pfurtscheller, 1999). The disappearance of coupling corresponds to the cerebral cortex idle hypothesis (Leocani et al., 1997; Gross et al., 2005). However, increasing evidence suggests that the motor cortex does not remain in an idle state during the PMBS period but is involved in motor modulation (Cassim et al., 2001; Gaetz and Cheyne, 2006; He et al., 2020) and is associated with movement parameters (Fry et al., 2016). For instance, greater PMBS has been observed following finger extension movements performed against a heavy load compared with unloaded extensions (Stančák et al., 1997). The increase in the rate of force development is responsible for increments in the amplitude of PMBS (Fry et al., 2016). In other words, faster movement speeds strengthen PMBS (Stančák et al., 1997). Some studies have found that PMBS presence in the contralateral hemisphere or frontal areas is associated with motor inhibition (Gilbertson et al., 2005; Pogosyan et al., 2009; Little et al., 2019; He et al., 2020; Wessel, 2020). PMBS is dependent on motor parameters, which may demonstrate the inhibition mechanism from unilateral movement to the lateral motor cortex. However, several questions remain elusive: (1) the role of the ipsilateral hemisphere in movement control during the PMBS burst; (2) the interaction between hemispheres during the PMBS period, particularly during index finger abduction movements at different speeds and the influence of speed on this interaction; and (3) the specific neuromodulation mechanism manifested by this phenomenon in the sensorimotor cortex.
Aiming to address the questions above, we combined transcranial magnetic stimulation (TMS) and EEG to explore the neural mechanisms of interhemispheric functional connections during the PMBS period. Hemispheric interactions have been studied using TMS pulses applied to the bilateral primary motor cortices (M1; Schulte and Müller-Oehring, 2010). In this paradigm, a conditional stimulus (CS) was applied to M1 in one hemisphere a few milliseconds before the test stimulus (TS) on the opposite side. If the CS affects TMS-induced motor–evoked potential (MEP) amplitude, it can be inferred that the two hemispheres are functionally connected (Rothwell, 2011).
Interhemispheric inhibition (IHI) is crucial in unilateral motion (Mayston et al., 1999; Beaulé et al., 2012). Numerous studies have revealed that IHI increases during unilateral movement (Talelli et al., 2008; Vercauteren et al., 2008; Mooney et al., 2018). It has also been found that inhibition from ipsilateral to contralateral hemispheres decreases during unilateral movement execution (Chen et al., 2003; Tazoe and Perez, 2013). However, due to scarce explorations, IHI from the ipsilateral to the contralateral hemisphere during the PMBS burst after unilateral movement remains unclear. Furthermore, it has been reported that ipsilateral to contralateral IHI can help control the time course of muscle activation and the motor learning process, leading to changes in motor response speed (Meyer and Voss, 2000; Davare et al., 2007). Therefore, understanding how the ipsilateral hemisphere regulates IHI during the PMBS burst at different speeds will help elucidate the functional significance and neural mechanisms of PMBS.
In the present study, we hypothesize that (1) during the MRBD and PMBS periods, the influence of index finger abduction at different speeds on interhemispheric interaction varies and (2) the interhemispheric interaction regulated by speed and PMBS amplitude is associated with the IHI mechanism from the ipsilateral to the contralateral hemisphere.
Materials and Methods
Subjects
Eight healthy, right-handed volunteers (three females, 23.2 ± 0.47 years old) were investigated in the present study. The present study is a preliminary, small-sample investigation. The sample size was determined based on previous studies employing similar neurophysiological techniques (Stancák and Pfurtscheller, 1996; Leocani et al., 1997; Manganotti et al., 1998; Mima et al., 2000; Serrien and Brown, 2002; Toma et al., 2002; Szurhaj et al., 2003; Gross et al., 2005; Yao et al., 2007; Stevenson et al., 2011; Qiu et al., 2016; Little et al., 2019). The mean number of participants in these studies was ∼7.33. We initially recruited 10 volunteers; however, 2 withdrew due to personal circumstances. Ultimately, eight subjects completed the experiment, which exceeds the average sample size of the aforementioned studies. Additionally, a power analysis using G*Power 3.1.9.7 was conducted to determine the adequacy of our sample size. The go/no-go paradigm employed in the present study is a well-established experimental paradigm. Its robustness has been repeatedly tested in previous research conducted by our group, with relatively large effect sizes consistently observed (Zhang et al., 2020, 2024). Therefore, based on previous research (Li et al., 2022; Himmelmeier and Werheid, 2024), the effect size was set to f = 0.4 and α = 0.05 for the statistical power calculation (Cohen, 2016). The analysis yielded a statistical power of 0.68, which is considered acceptable for neuroelectrophysiological studies (Faul et al., 2007; Häger et al., 2021). All participants were confirmed to be in good physical and mental health, with no dyskinesia or mental illness. Handedness was assessed using the Edinburgh Handedness Inventory (Oldfield, 1971). Written informed consent was obtained from all participants prior to the recordings. The experimental procedures were approved by the local ethics committee at Jinan University and were conducted in accordance with the guidelines established in the Declaration of Helsinki.
EEG and EMG recordings
EEG data were acquired from 64 scalp sites (extended 10–20 system) using a cap with active Ag/AgCl electrodes (Quickcap64). Wet electrodes were used in the cap, and electrode impedance was maintained at <5 KΩ. The reference electrode was placed on the bilateral mastoid. A Neuroscan-Synamps2 amplifier was used to amplify the EEG signal. The EEG signal was sampled at 1 kHz. A bandpass filter ranging from 0.5 to 40 Hz was applied, along with a 50 Hz notch filter. During off-line data preprocessing, independent component analysis (ICA) using the Infomax algorithm from the EEGlab toolbox (Delorme and Makeig, 2004) was performed. Following visual inspection of scalp maps and time course activations for each ICA component, those clearly associated with eyeblinks were eliminated. The remaining subset of components was back-projected onto the EEG data for subsequent analysis.
The first dorsal interosseous (FDI) muscle was selected for index finger isometric abduction tasks (Yao et al., 2007). Surface electromyography (EMG) signals were recorded bilaterally from the FDI using Ag–AgCl electrodes (10 mm diameter). The EMG signals were amplified and filtered with a bandwidth of 5–2 kHz using a bioamplifier (Neurolog System, Digitimer). The signals were digitized at a rate of 5 kHz using an A/D converter [CED Micro 1401, Cambridge Electronic Design (CED)] and stored on a laboratory computer. Online visual display and off-line analysis were conducted using a custom data collection and conditional averaging software (Spike2 for Windows, Version 9.11, CED). The marks generated by the presentation software (E-Prime 3.0) were transmitted to the EEG acquisition software (Curry 8) for synchronization purposes. Simultaneously, the EMG and TMS equipment acquisition software (Spike2) operated to ensure precise timing and synchronization during the experiment.
Experimental paradigm
The subjects were seated comfortably in an armchair with both arms flexed at the elbow joint at 90° and supported by a table. During the experiment, the left hand was kept relaxed and rested on the table with the wrist in pronation.
The configuration of the experimental procedure was adapted from previous studies (Ferbert et al., 1992; Yao et al., 2007; Stevenson et al., 2011; Fry et al., 2016; Qiu et al., 2016) and is presented in Figure 1A. A fixation cross and a white circle appeared alternately on the screen. At the beginning of each trial, a fixation cross appeared in the center of the screen for a variable delay (3.25–3.75 s). Subjects remained relaxed during this period. Then, an imperative stimulus consisting of a white circle appeared for 0.20 s to denote the performance of an isometric abduction movement with the right index finger. There were two speeds for the isometric abduction movement: ballistic movement and self-paced movement. PMBS could be observed after the isometric abduction movement (Fig. 1B). When PMBS peaked, TMS pulses were applied to the primary sensorimotor areas of the left and right hemispheres. The position of the tested right hand is shown in Figure 1D. Pegs were used to limit the movement range of the index finger. The distance between the two pegs was adjusted for each subject's comfort. Once determined, the peg distance remained fixed throughout the experiment to ensure equidistant index finger abduction movement. Each speed condition included three blocks, each encompassing 20 trials. The number of trials employed in the present study is consistent with those used in previous research (Ferbert et al., 1992; Fry et al., 2016). Furthermore, we calculated the time–frequency representation (TFR), which has been demonstrated to exhibit a superior signal-to-noise ratio compared with event-related potentials (Baker et al., 1997; Morales and Bowers, 2022).
Experimental procedure. A, A schematic diagram of the single experiment task flow. Subjects were instructed to perform index finger abduction movements under different conditions while the white circle was appearing. B, The beta waveforms, which include two stages: MRBD during the execution of the motion and PMBS after movement while performing the task. C, The EMG sequence and MEP from a single task. D, A schematic diagram of the index finger abduction movement and the brain areas stimulated by TMS. MEP, motor-evoked potential; TMS, transcranial magnetic stimulation; MRBD, movement-related beta desynchronization; PMBS, postmovement beta synchronization.
Additionally, another experiment was performed for comparison. The experimental paradigm was similar to the previous one, with the only difference being the timing of the TMS pulses. In this experiment, TMS pulses were triggered as the PMBS returned to baseline, rather than during the PMBS peak. Prior to the experiments, we determined the time points when PMBS peaked and returned to the baseline. Hence, the experiment with isometric abduction movement was exclusively conducted at the beginning (Zhang et al., 2020). There were two speeds matched for ballistic movement and self-paced movement, respectively, in the index finger abduction experiment. Three blocks, each containing 30 trials in randomized order, were tested, with 15 trials for ballistic movement and the remaining trials for self-paced movement. After each trial, subjects were required to report whether they had missed the task or performed the wrong movement. Additionally, we checked the raw data during off-line analysis to identify any missed trials. For further analysis, we discarded the trials where participants either missed or performed incorrectly. If the discarded trials exceeded 20%, the block was recollected. This led to the allowable exclusion of zero to six trials per block per subject; in practice, no block needed to be recollected.
TMS
Two figure-of-eight-shaped coils (90 mm mean diameter) were connected to two Magstim 2002 magnetic stimulators (Magstim) to deliver the CS and TS, respectively. The CS was given over the ipsilateral motor cortex (IM1) and TS over the contralateral motor cortex (CM1). The handles of the coils were pointed backward and 45° away from the midline and then moved in 1 cm steps to identify the position where TMS pulses produced maximal amplitude responses (Fig. 1D). The optimal site where stimulation of slight suprathreshold intensity consistently produced the largest MEPs in the contralateral resting FDI muscles was marked with a pen for each hemisphere on the swimming cap-covered scalp. During the test, the TMS coil was held over the head by the operator's hands.
Resting motor threshold (RMT) was defined as the lowest TMS intensity that evoked an MEP amplitude above 50 μV peak-to-peak in the relaxed FDI muscle at least five times out of 10 consecutive stimuli (Rothwell et al., 1999). Additionally, we collected the TMS measurements of IHI in the following four experimental conditions: (1) PMBS induced by ballistic abduction and at peak time (ballistic, peak); (2) PMBS induced by ballistic abduction and at recovering baseline time (ballistic, baseline); (3) PMBS induced by self-paced abduction and at peak time (self-paced, peak); and (4) PMBS induced by self-paced abduction and at recovering baseline time (self-paced, baseline). Each condition consisted of three blocks with 20 trials included each. To prevent fatigue, conditions were separated by a 10 min break.
IHI
IHI from the right primary motor cortex (IM1) to the left primary motor cortex (CM1) was measured using a double-pulse TMS protocol, as previously described (Ferbert et al., 1992). The Bistim module delivered two magnetic stimuli sequentially to IM1 and CM1 via separate stimulating coils, allowing investigation of the effect of the first stimulus (CS) on the second stimulus (TS). A suprathreshold CS was applied to the optimal scalp position over IM1, followed by a suprathreshold TS delivered to CM1 10 ms later, consistently evoking IHI (Ferbert et al., 1992). The evaluation of IHI was conducted in accordance with the protocol developed by Talelli et al. (2008) and Zhang et al. (2024). The stimulation intensity for both the CS and TS was adjusted to evoke MEPs ranging from 1 to 1.5 mV in the contralateral FDI muscle at rest. The stimulation intensities of the CS and TS were maintained at consistent levels, with panel intensities set at the same value, and the adjustment range was between 110 and 120% of the RMT. Specifically, the TS intensity was set at 115 ± 5% of RMT, while the CS intensity was set at 114 ± 5% of RMT. Furthermore, when adjusting the TS intensity, it was imperative to ensure that in 10 consecutive stimulations, at least five MEP values in the right-hand FDI muscle fell between 1.5 and 2 mV to maintain consistent TS levels across different conditions. In a single-pulse mode, a suprathreshold TS was delivered exclusively to CM1 either at PMBS peak time or during recovery to the baseline. Both single-pulse and paired-pulse modes were randomly administered, with each condition involving 30 single TS MEPs and 30 CS MEPs. IHI was computed by expressing the amplitude of CS MEPs as a percentage of the amplitude of single TS MEPs (CS MEP * 100 / single TS MEP).
Data processing
EEG analysis
EEG was analyzed using the MATLAB-based FieldTrip toolbox (Matlab, 2021b), developed by the Donders Institute for Brain, Cognition, and Behavior (Oostenveld et al., 2011; http://www.ru.nl/neuroimaging/fieldtrip/), and the EEGLAB toolbox (Delorme and Makeig, 2004), developed by the Swartz Center for Computational Neuroscience (https://sccn.ucsd.edu/eeglab/). The response-locked EEG time series (−0.25 s before and 2.5 s after the stimulus) were extracted for artifact screening. Artifacts were identified as EEG signals with amplitudes exceeding M ± 5 SD. Here, M and SD were computed for each time point based on the response-locked EEG time series across all trials. The study specifically focused on channels C3 and C4 of IM1 and CM1 for analysis.
TFR
TFR was calculated across trials using a fast Fourier transform (multi)taper approach with short sliding time windows (Percival and Walden, 1993; Osipova et al., 2006). A frequency band from 5 to 40 Hz (in steps of 0.25 Hz) was analyzed, employing an adaptive time window of three cycles for each frequency (ΔT = 3 / f). The Hanning taper was employed to minimize spectral leakage and control frequency smoothing. Time windows were advanced in steps of 20 ms. EEG power TFR was computed by averaging the squared absolute values of convolutions across trials. The EEG power at each frequency and time point was normalized by subtracting the average power during the baseline period and then dividing it by the average baseline power. The resulting value was then multiplied by 100, thus expressing the event-related EEG power change relative to the baseline. The aforementioned procedure is described in detail by Nakayashiki et al. (2014) as follows:
The baseline was defined as the time window from 0 to 250 ms prior to the onset of the movement task ( Stančák et al., 1997; Fry et al., 2016; Stevenson et al., 2011). Values above 0 indicated that power at that frequency and time point was higher than the average baseline power and vice versa.
Event-related power changes
Event-related power changes in the beta band (14–30 Hz) were assessed. In both the ballistic and self-paced tasks, the duration of MRBD and PMBS typically exceeds 0.5 s (Stancák and Pfurtscheller, 1996). To visualize alterations over different time courses during the abduction tasks, we reduced the dimension to the time domain by averaging the values of beta rhythms in the time–frequency domain (as shown in Fig. 4).
The value of PMBS was defined as the average normalized power over a 200 ms window centered on the peak of the power change. Conversely, the MRBD was defined as the average normalized power over a 200 ms window centered on the trough of the power change during movement. The onset of MRBD was marked by the time at which beta power declined by 10% compared with the baseline (Zhang et al., 2008). Similarly, the end time of MRBD was defined as the point at which beta power returned to 10% below the baseline. The latency of PMBS was identified as the time at which beta power increased by >10% above the baseline, while the termination of PMBS was marked by the return of beta power to 10% above the baseline. The peak time of PMBS was defined as the time point with the maximum beta power during the PMBS burst. The movement duration was constrained within the interval between the initiation and termination of MRBD/PMBS.
It is acknowledged that bilateral PMBS appears at a fixed time in the motor cortex during index finger abduction tasks (Pfurtscheller and Neuper, 1997). The peak and recovery baseline times of PMBS can be calculated for bilateral M1. In the present study, we calculated the average PMBS peak and recovery baseline times across all trials in each speed condition when TMS pulses were triggered. Our work primarily focused on the contribution of IM1 during ipsilateral PMBS bursts. Therefore, the peak time of ipsilateral PMBS was selected to trigger TMS pulses, while the PMBS recovery baseline time was set to trigger TMS for the control condition.
Coherence between IM1 and CM1
To detect the functional coupling between bilateral M1 in the frequency domain (5–40 Hz), coherence and Granger causality were calculated. This method measures the proportion of the signal within a particular frequency band that maintains a constant phase relationship between two sources, with values ranging from 0 to 1. For EEG–EEG coherence calculation, each trial of 2.75 s was divided into 50-ms-long nonoverlapping sections. The specific procedure was as follows (Baker et al., 1997):
If the FFT derived from signal i (i = 1,2) over the lth section (l = 1, …, L) is Fi,l(f), then the cross-spectrum is as follows:
In order to process coherence values across subjects, individual coherence values were initially calculated for each subject using the previously described method. Subsequently, the values were averaged across all subjects in order to obtain group-level coherence estimates.
Causality of directed coherence
Directed coherence measures the extent to which one signal can be predicted by the past history of another (Kamiński and Blinowska, 1991). In this case, Dir-Coh can be interpreted as a measure of causality. Let X(t) be the activity of two EEG signal at time t. Estimating directed coherence is to fit a multivariable autoregressive model to the observed signals. Let the following vector:
The significance level for the directed coherence was calculated as follows (Baker et al., 2006):
To process the directed coherence values across subjects, the individual directed coherence values were initially calculated for each subject using the previously described method. Subsequently, the values were averaged across all subjects in order to obtain group-level estimates.
Statistical analysis
Repeated-measure two–way ANOVAs were performed to determine the effect of speed (ballistic, self-paced) and hemispheres (CM1, IM1) on EEG parameters (MRBD values, PMBS values, MRBD/PMBS onset time, duration, end time). Repeated-measures two–way ANOVAs were also conducted to study the effect of speed (ballistic, self-paced) and time course (MRBD period, PMBS period) on interhemispheric coherence and interhemispheric directed coherence. Additionally, repeated-measure two–way ANOVAs were performed to examine the effect of speed (ballistic, self-paced) and the time when TMS pulses were triggered (PMBS peak time, PMBS recovery baseline time) on IHI. Bonferroni’s post hoc correction was conducted to detect significant comparisons with a significance level α of 0.05. A priori comparisons were made as specified. Normal distribution was tested by the Shapiro–Wilk test (all p > 0.05). The significance level was set at p < 0.05, and group data are presented as M ± SD in the text and as SE in the figures.
We explored the correlation among speed, PMBS, coherence, directed coherence, and IHI. The simplest model, a linear model, was chosen to describe the effects of these correlations. This study is an extension of our previous research, which employed a linear model to elucidate the relationships between neural data (Zhang et al., 2020, 2024). In addition to the aforementioned studies, linear models have been demonstrated to offer a robust and straightforward methodology for elucidating correlations and causality in intricate datasets while avoiding the pitfalls of overcomplicated analysis (Gourévitch et al., 2006; Ritchie et al., 2019; Penconek, 2022). Therefore, Pearson's correlation analysis was conducted to evaluate the relationships using IBM SPSS (SPSS) with 1,000 bootstrap samples and a 95% confidence interval.
Results
Preexperiment
As shown in Figure 2, subjects were required to perform two kinds of abduction at different speeds to evoke beta event-related spectral disturbances (ballistic speed, 8.04 ± 0.93 cm/s; self-paced, 4.26 ± 0.70 cm/s; p = 0.002). A decrease in beta power was followed by a rebound synchronization after movement (Fig. 3A,B). MRBD and PMBS were consistently observed over bilateral hemispheres during abduction movement at varied speeds, with contralateral dominance at both speeds (Fig. 3A,B).
Results of index finger abduction movement at different speeds. A,B, EMG sequences after rectification. C, Average movement speeds across subjects during ballistic and self-paced abduction movements. Error bars indicate SEs. *p < 0.05.
Change in beta rhythm power over the sensorimotor cortex (CM1 and IM1 channels). A,B, Power changes in the time domain from a representative subject during ballistic and self-paced abduction movement. C, Average MRBD and PMBS values. The abscissa shows the experimental conditions, and the ordinate shows the average power change across subjects. In panels A and B, the first two rectangles mark MRBD and the second two rectangles denote PMBS. Error bars indicate SEs. *p < 0.05; MRBD, movement-related beta desynchronization; PMBS, postmovement beta synchronization; IM1, ipsilateral motor cortex; CM1, contralateral motor cortex.
A significant effect of hemisphere (CM1, IM1) was demonstrated for the MRBD value (Fig. 3C; (F(1,31) = 12.737; p = 0.003). Post hoc analysis indicated that MRBD was larger on CM1 compared with IM1 under the condition of ballistic abduction (CM1, −45.76 ± 5.16%; IM1, −37.83 ± 4.85%; p = 0.009) and self-paced abduction (CM1, −42.34 ± 4.16%; IM1, −37.02 ± 4.16%; p = 0.020). Additionally, MRBD showed no significant difference between ballistic and self-paced abduction on CM1 (p = 0.248) and IM1 (p = 0.780).
For PMBS values (Fig. 3C), a significant effect of speed (ballistic, self-paced; F(1,31) = 21.078; p < 0.001) and hemisphere (CM1, IM1; F(1,31) = 14.064; p = 0.003) was observed. Post hoc analysis revealed that PMBS was larger on CM1 compared with IM1 under ballistic abduction (CM1, 83.30 ± 12.43%; IM1, 58.64 ± 8.82%; p = 0.002) and similarly larger under self-paced abduction (CM1, 55.27 ± 9.68%; IM1, 38.24 ± 7.01%; p = 0.026). Additionally, PMBS was stronger during ballistic abduction compared with self-paced abduction on both CM1 (p < 0.001) and IM1 (p = 0.007).
As illustrated in Figure 4, the coherence between CM1 and IM1 was quantified. During the baseline period, there was no significant difference in coherence between ballistic and self-paced movements (F(1,14) = 0.149; p = 0.705). Therefore, baseline analyses were not further pursued. One-sample t tests revealed that during the MRBD period, significant beta power above zero was observed in both ballistic (Channel C3, t = 4.636; p = 0.001; channel C4, t = 4.333; p = 0.002) and self-paced (channel C3, t = 14.454; p < 0.001; Channel C4, t = 5.096; p < 0.001) movements. This verified that there was enough beta power during MRBD period to perform coherence analysis. Significant effects of speed (ballistic, self-paced; F(1,31) = 12.030; p = 0.005) and time course (MRBD period, PMBS period; F(1,31) = 14.485; p = 0.003) were observed, as well as a significant effect of their interaction (F(1,31) = 7.788; p = 0.016; Fig. 4C). Post hoc analysis elucidated that coherence was not significantly different between ballistic and self-paced movement during the MRBD period (ballistic, 0.14 ± 0.03; self-paced, 0.12 ± 0.03; p = 0.653). However, coherence under the ballistic condition was increased compared with the self-paced condition during the PMBS period (ballistic, 0.43 ± 0.08; self-paced, 0.24 ± 0.07; p < 0.001). Coherence was also larger during the PMBS period than the MRBD period under both the ballistic (p < 0.001) and self-paced (p = 0.032) conditions (Fig. 4B,C). Although more prominent coherence between bilateral hemispheres was detected during the PMBS period, the direction of information flow between bilateral hemispheres cannot be determined from these results.
Coherence between CM1 and IM1 at different speeds. A, Coherence in different time–frequency domains from a representative subject (averaged across trials) during speed movement tasks. The dashed lines represent the beta band between 14 and 30 Hz. B, Coherence in the frequency domain during different speed movements. C, Average coherence across all subjects. The abscissa shows the conditions. The ordinate shows the average coherence across subjects. Error bars indicate SEs. *p < 0.05; MRBD, movement-related beta desynchronization; PMBS, postmovement beta synchronization.
Figure 5 presents the results of directed coherence during different time courses of MRND and PMBS, which reflects the causality between the two hemispheres. During the MRBD period, a notable effect of coherence direction (CM1 to IM1, IM1 to CM1) was observed (F(1,31) = 11.220; p = 0.015; Fig. 5C). Post hoc analysis revealed that directed coherence of causality was significantly larger from CM1 to IM1 compared with IM1 to CM1 during the MRBD period induced by both ballistic abduction (CM1 to IM1, 0.087 ± 0.03; IM1 to CM1, 0.037 ± 0.03; p = 0.015) and self-paced abduction (CM1 to IM1, 0.081 ± 0.04; IM1 to CM1, 0.033 ± 0.02; p = 0.018). However, directed coherence of causality did not significantly differ between ballistic and self-paced movements for either direction (CM1 to IM1, p = 0.559; IM1 to CM1, p = 0.656). Therefore, it can be concluded that the causal direction of information interaction is predominant from CM1 to IM1 during the MRBD period, with no significant difference between the two speed conditions.
Directed coherence at different time points of MRBD and PMBS. A,B, Directed coherence across different frequency domains from a representative subject (averaged across trials) during the MRBD period. C, The average directed coherence during the MRBD period across all subjects. E,F, Directed coherence in different frequency domains from a representative subject (averaged across trials) during the PMBS period. G, The average directed coherence in the PMBS period across all subjects. The abscissa indicates the conditions, while the ordinate represents the average coherence across subjects. Dashed lines represent the threshold of directional coherence significance. Error bars denote SEs. *p < 0.05; MRBD, movement-related beta desynchronization; PMBS, postmovement beta synchronization; IM1, ipsilateral motor cortex; CM1, contralateral motor cortex.
During the PMBS period, the effects of directed coherence (CM1 to IM1, IM1 to CM1) were found to be significant (F(1,31) = 53.097; p < 0.001), as were the effects of speed (ballistic, self-paced; F(1,31) = 7.707; p = 0.024; Fig. 5G). Post hoc analysis revealed that the causal direction of coherence was greater from IM1 to CM1 compared with CM1 to IM1 during the PMBS period for both ballistic (CM1 to IM1, 0.03 ± 0.03; IM1 to CM1, 0.11 ± 0.04; p < 0.001) and self-paced movements (CM1 to IM1, 0.03 ± 0.03; IM1 to CM1, 0.08 ± 0.04; p = 0.001). Directed coherence from IM1 to CM1 was significantly increased for ballistic compared with self-paced movements (p = 0.004). However, during the PMBS period, the causal direction from CM1 to IM1 did not differ significantly between ballistic and self-paced conditions (p = 0.632). Thus, it can be concluded that the causality of information interaction is predominant from IM1 to CM1 during the PMBS period, with a greater effect observed in ballistic than self-paced movement.
Based on the results presented in Figures 4 and 5, coherence during the PMBS period is greater than in the MRBD period, with a stronger causal influence from IM1 to CM1 during PMBS, particularly induced by ballistic movement. This potentially reflects that the functional interaction between the bilateral hemispheres is not independent during the PMBS period, as supported by findings in IHI across the corpus callosum from IM1 to CM1. The results of ANOVAs are presented in Table 1.
Results of the ANOVAs
Correlation analysis during MRBD and PMBS periods
As illustrated in Figure 6A–C, no correlations were identified between speed and MRBD or hemispheric coupling. The directed coherence from CM1 to IM1 was not found to be dependent on speed during the MRBD period. In contrast, Figure 6D–F revealed positive correlations between PMBS and speed (r = 0.734; p = 0.001) and between hemispheric coupling and speed (r = 0.522; p = 0.038). Additionally, directed coherence from CM1 to IM1 was dependent on speed (r = 0.568; p = 0.022) during the PMBS period.
The scatterplots illustrating significant associations. The relationships between speed and EEG parameters are shown. During the MRBD period, no significant correlations were found between speed and MRBD (A), coherence (B), or directed coherence (C). In contrast, during the PMBS period, significant correlations were observed between speed and PMBS (D), coherence (E), and directed coherence (F). The solid and dotted lines represent the regression lines and the 95% confidence limits, respectively. MRBD, movement-related beta desynchronization; PMBS, postmovement beta synchronization.
As illustrated in Figure 7A–C, no correlation was identified between MRBD and coherence or directed coherence during the MRBD period. Conversely, Figure 7D–F illustrates a positive correlation between PMBS and coherence (r = 0.604; p = 0.013) as well as directed coherence (r = 0.631; p = 0.009) during the PMBS period. Additionally, the correlation between coherence and directed coherence is positive in both MRBD and PMBS periods (MRBD period, r = 0.770; p < 0.001; PMBS period, r = 0.694; p = 0.003).
The scatterplots with significant associations. During the MRBD period, insignificant correlations were observed between MRBD and (directed) coherence (A,B), while a significant correlation between coherence and directed coherence was found (C). In the PMBS period, significant correlations were identified between PMBS and (directed) coherence (D,E), as well as between coherence and directed coherence (F). The solid and dotted lines represent the regression lines and the 95% confidence limits, respectively. MRBD, movement-related beta desynchronization; PMBS, postmovement beta synchronization.
The results of the correlation analysis conducted during the MRBD and PMBS periods are presented in Tables 2 and 3, respectively. Since interhemispheric (directed) coherence was related to PMBS, we further assessed IHI during the PMBS period.
Results of the correlation during MRBD period
Results of the correlation during the PMBS period
Results of IHI
As illustrated in Figure 8, A and B, no notable impact of speed (ballistic, self-paced), the timing of TMS pulse initiation (PMBS peak time, PMBS recovery baseline time), or their interaction on TS MEP was discerned. There was no difference in TS MEP across the four conditions. Therefore, correction analysis for TS MEP is unnecessary and was excluded from this paper.
Single-pulse MEP under different conditions. A, Single-pulse MEP from a representative subject (averaged across trials) during various conditions. B, Average of single MEP across subjects. Error bars indicate SEs. *p < 0.05; PMBS, postmovement beta synchronization; MEP, motor-evoked potential.
As illustrated in Figure 9, the analysis revealed that there were significant effects of speed (ballistic, self-paced; F(1,31) = 16.406; p = 0.005), timing of TMS triggering (PMBS peak time, PMBS recovery baseline time; F(1,31) = 73.142; p < 0.001), and their interaction (F(1,31) = 16.415; p = 0.005) on IHI. Post hoc analysis showed that IHI during PMBS peak time was significantly strengthened compared with recovery baseline time under both ballistic movement (PMBS peak time, 0.52 ± 0.09; PMBS recovery baseline time, 0.72 ± 0.10; p < 0.001) and self-paced movement (PMBS peak time, 0.67 ± 0.07; PMBS recovery baseline time, 0.74 ± 0.11; p = 0.015). Moreover, IHI during ballistic movement was significantly more pronounced than during self-paced movement at PMBS peak time (p < 0.001). However, there was no significant difference in IHI between ballistic and self-paced movement during PMBS recovery time (p = 0.589). Therefore, we can infer that IM1 transmits more effective inhibitory regulation to CM1 during PMBS, which is enhanced under ballistic conditions compared with self-paced conditions. The results of ANOVAs of IHI are presented in Table 1.
IHI across different conditions. A, IHI from a representative subject (averaged across trials) during various conditions. B, Group data. The abscissa shows different conditions, while the ordinate displays IHI across subjects, with each gray circle representing an individual value. Error bars indicate SEs. *p < 0.05; PMBS, postmovement beta synchronization; MEP, motor-evoked potential.
To explore the association of IHI with four parameters (speed, PMBS, coherence during PMBS, and directed coherence in the PMBS period), we computed the correlation coefficients using Pearson's correlation analysis (Fig. 10). A significant negative correlation was found between speed and IHI values (r = −0.70; p = 0.003), PMBS and IHI values (r = −0.529; p = 0.035), and directed coherence from IM1 to CM1 and IHI values (r = −0.631; p = 0.009). Larger PMBS values corresponded to smaller IHI values, indicating a stronger inhibitory effect from the right hemisphere to the left hemisphere. The results of the correlation between IHI during the MRBD and PMBS periods are presented in Tables 2 and 3, respectively.
Scatterplots of significant associations. The relationships of IHI with speed (A), PMBS (B), coherence during PMBS (C), and directed coherence from IM1 to CM1 (D). The solid and dotted lines represent the regression lines and the 95% confidence limits, respectively. PMBS, postmovement beta synchronization; IM1, ipsilateral motor cortex; CM1, contralateral motor cortex; IHI, interhemispheric inhibition.
Discussion
Interhemispheric asymmetry at different speeds
The present study demonstrated that MRBD is bilaterally present during movement execution, followed by bilateral PMBS, consistent with previous research (Little et al., 2019; Pfurtscheller and Da Silva, 1999; Wessel, 2020; Zhang et al., 2020, 2024). These findings corroborate that both MRBD and PMBS exhibit contralateral dominance during unilateral movement (Pfurtscheller et al., 1997; Leocani et al., 2001; Szurhaj et al., 2003; Rothwell et al., 2021; Wischnewski et al., 2022; Long et al., 2024), highlighting the role of beta oscillations in movement initiation and cancellation (Little et al., 2019; Wessel, 2020). Additionally, Correia et al. (2022) emphasized the importance of investigating fast limb movements. Our findings extend this understanding by demonstrating that PMBS, but not MRBD, is highly velocity-dependent. Larger PMBS was observed in ballistic movements compared with self-paced movements, consistent with previous studies (Zhang et al., 2020; Correia et al., 2022).
Interhemispheric connection during the MRBD period
During the MRBD period, interhemispheric connections were evident during movement execution at different speeds. However, in the present study, the influence of speed was not observed. Coherence appears to be linked to motion complexity rather than specific motion tasks. Nevertheless, the impact of velocity on coherence remains a topic of debate. Toma et al. (2002) suggested that motor cortical activation and coupling were more pronounced during faster movements. However, Serrien and Brown (2002) found that the coupling between the primary sensorimotor cortices in the beta frequency band decreased with increasing movement speed. This effect was particularly strong in the antiphase mode compared with the in-phase mode. In the present study, the impact of speed on coherence was not detected, possibly due to differences in experimental design. The directed coherence analysis revealed that the predominant drive was from CM1 to IM1, consistent with previous studies (Amunts et al., 1996; Gross et al., 2005; Chettouf et al., 2020; Säisänen et al., 2021; Takasawa et al., 2022). This finding suggests that both hemispheres are involved in the execution of unilateral movements (Amunts et al., 1996; Gross et al., 2005; Beaulé et al., 2012; Chettouf et al., 2020; Merrick et al., 2022). Notably, interhemispheric coherence and spectral power desynchronization are distinct phenomena (Manganotti et al., 1998; Bundy et al., 2018), which may explain the lack of correlation observed between MRBD and coherence during movement execution.
Interhemispheric connection during the PMBS period
The lack of postmotion bilateral coupling has been noted in previous studies (Leocani et al., 1997; Andrew and Pfurtscheller, 1999; Chettouf et al., 2020). Paradoxically, we observed a transient increase in interhemispheric coupling after motor termination. PMBS is task-dependent and closely related to motion parameters (Barone and Rossiter, 2021; Fry et al., 2016; Kang et al., 2023; Zhang et al., 2020, 2024). Therefore, the possible reason for the contradiction might be the variation within the task design.
An increase in task complexity is associated with enhanced intercerebral coherence (Gross et al., 2005; Aliakbaryhosseinabadi et al., 2021), which may explain the observed discrepancies. In contrast to studies investigating simple, self-paced movements (Leocani et al., 1997; Andrew and Pfurtscheller, 1999), we introduced more constrained and precisely timed tasks that demand higher motor control and coordination. These more challenging tasks may prompt the brain to exceed the threshold for interhemispheric interaction compared with simpler tasks. Another critical factor in motor cortical dynamics is movement speed. Previous studies primarily focused on self-paced movements without controlling for speed (Leocani et al., 1997; Andrew and Pfurtscheller, 1999). In contrast, our study tested two distinct conditions: ballistic (∼130 ms) and self-paced (∼180 ms). This highlights the pivotal role of speed in interhemispheric interactions, with earlier studies potentially lacking coherence due to their exclusive focus on self-paced movements (Leocani et al., 1997; Andrew and Pfurtscheller, 1999; Gross et al., 2005).
Our findings revealed a correlation between movement velocity and interhemispheric connectivity. During the PMBS period, directed coherence from IM1 to CM1 is stronger in ballistic movements than in self-paced movements, indicating greater functional coordination between IM1 and CM1. The increase in PMBS in faster movements may necessitate stronger interhemispheric connectivity. Prior studies have demonstrated that contralateral PMBS reflects cortical inhibition (Gilbertson et al., 2005; Pogosyan et al., 2009; Little et al., 2019; Hervault et al., 2021; Marques et al., 2022). However, coherence during PMBS only indicates that two regions are communicating, regardless of whether this communication is inhibitory or excitatory (Manganotti et al., 1998; Bundy et al., 2018). It is possible that interhemispheric connectivity during PMBS does not directly reflect the neuroregulatory mechanisms (excitability or inhibition) between IM1 and CM1. To further investigate this phenomenon, we employed paired-pulse TMS to assess IHI across the corpus callosum. The results, as illustrated in Figure 10, demonstrate a negative correlation between IHI and PMBS, coherence during PMBS, and directed coherence from IM1 to CM1. This indicates that interhemispheric connectivity reflects an inhibitory mechanism from the right to the left hemisphere during PMBS, with the ipsilateral hemisphere playing an active role in motor control. Furthermore, a negative correlation was observed between movement speed and IHI, indicating that faster movements are associated with stronger inhibition from IM1 to CM1 (Fig. 10A). Greater inhibition at higher speeds may assist in suppressing undesired mirror movements and enhancing the precision of motor commands (Ferbert et al., 1992; Meyer et al., 1995; Correia et al., 2022). Our findings lend support to the use of PMBS as a biomarker for movement control and suggest the potential for the development of brain stimulation technologies to enhance motor-related neural circuits, thereby promoting rapid recovery through increased neural activity and plasticity.
The greater directed coherence from IM1 to CM1 in ballistic movements compared with self-paced movements during the PMBS period suggests a significant role for the ipsilateral hemisphere in movement control. This coherence may reflect an inhibitory mechanism from the right to the left hemisphere. Reinforced suppression of IM1 to CM1 in ballistic movements potentially reduces unwanted or residual right-hand movements (Hübers et al., 2008; Cabibel et al., 2020; Iwama et al., 2020). This would support motor stability of the motor system, which aligns with our previous findings on PMBS (Zhang et al., 2020, 2024).
Practical utility
The present study elucidates the excitatory and inhibitory neurophysiological mechanisms of PMBS and the role in across hemispheric information exchange. Together with previous studies from our group, the robustness of PMBS as a biomarker for the stabilization of the movement system has been repeatedly demonstrated (Zhang et al., 2020, 2024). Therefore, PMBS may have the potential to serve as a reliable neural index for brain–computer interface (BCI) paradigms. Through a real-time update module, we can monitor the temporal dynamics of PMBS and decode the state of the movement system. This constitutes the core module of a closed-loop, brain-controlled repetitive TMS (rTMS) system with the potential to enhance postapoplectic motor rehabilitation. Furthermore, using source reconstruction algorithms such as weighted minimum-norm estimation, PMBS can be localized to cortical and subcortical regions, offering a neurophysiological foundation for the optimization of target areas in rTMS and deep brain stimulation interventions (Molins et al., 2008; Chang et al., 2015; Krishnaswamy et al., 2017; Li et al., 2019).
Conclusion
The present study demonstrated that speed not only exerts control over PMBS but also influences interhemispheric functional connectivity, thereby regulating IHI during the PMBS period. The occurrence of ballistic movements has been observed to result in increased suppression from IM1 to CM1, thereby underscoring the pivotal role of the ipsilateral hemisphere in motor control. The findings in the present study highlight the potential of PMBS as a biomarker of motor state, which may pave the way for the development of a PMBS-driven closed–loop rTMS system for postapoplectic motor rehabilitation via BCI paradigms.
Data Availability
The raw and preprocessed data generated in this study have been deposited in a local database. The EEG and EMG data are available under restricted access as they contain personally identifiable information, and patients have not consented to data distribution. Access can be obtained from the corresponding author upon reasonable request. The code for this article can be accessed online at https://github.com/Morinooku/code-for-eN-TNWR-0370-24X-.git.
Footnotes
The authors declare no competing financial interests.
This research was supported by the funding from the National Natural Science Foundation of China (Grant Number, 82260364, 61773179), Gansu Provincial Science and Technology Department Youth Fund Project (Grant Number, 22JR5RA166), the Guangdong Provincial Natural Science Foundation of China (Program Number 2019A1515012175).
↵*X.Z. and S.Z. contributed equally to this work.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.