Spatial variability in cortex-muscle coherence investigated with magnetoencephalography and high-density surface electromyography

EMG-projected MEG high-resolution source imaging of human motor execution: Brain-muscle coupling above movement frequencies

Magnetoencephalography (MEG) is a non-invasive functional imaging technique for pre-surgical mapping. However, movement-related MEG functional mapping of primary motor cortex (M1) has been challenging in presurgical patients with brain lesions and sensorimotor dysfunction due to the large numbers of trials needed to obtain adequate signal to noise. Moreover, it is not fully understood how effective the brain communication is with the muscles at frequencies above the movement frequency and its harmonics. We developed a novel Electromyography (EMG)-projected MEG source imaging technique for localizing early-stage (-100 to 0 ms) M1 activity during ~l min recordings of left and right self-paced finger movements (~1 Hz). High-resolution MEG source images were obtained by projecting M1 activity towards the skin EMG signal without trial averaging. We studied delta (1-4 Hz), theta (4-7 Hz), alpha (8-12 Hz), beta (15-30 Hz), gamma (30-90 Hz), and upper-gamma (60-90 Hz) bands in 13 healthy participants (26 datasets) and three presurgical patients with sensorimotor dysfunction. In healthy participants, EMG-projected MEG accurately localized M1 with high accuracy in delta (100.0%), theta (100.0%), and beta (76.9%) bands, but not alpha (34.6%) or gamma/upper-gamma (0.0%) bands. Except for delta, all other frequency bands were above the movement frequency and its harmonics. In three presurgical patients, M1 activity in the affected hemisphere was also accurately localized, despite highly irregular EMG movement patterns in one patient. Altogether, our EMG-projected MEG imaging approach is highly accurate and feasible for M1 mapping in presurgical patients. The results also provide insight into movement-related brain-muscle coupling above the movement frequency and its harmonics.

EMG-projected MEG High-Resolution Source Imaging of Human Motor Execution: Brain-Muscle Coupling above Movement Frequencies

Magnetoencephalography (MEG) is a non-invasive functional imaging technique for pre-surgical mapping. However, movement-related MEG functional mapping of primary motor cortex (M1) has been challenging in presurgical patients with brain lesions and sensorimotor dysfunction due to the large numbers of trails needed to obtain adequate signal to noise. Moreover, it is not fully understood how effective the brain communication is with the muscles at frequencies above the movement frequency and its harmonics. We developed a novel Electromyography (EMG)-projected MEG source imaging technique for localizing M1 during ~1 minute recordings of left and right self-paced finger movements (~1 Hz). High-resolution MEG source images were obtained by projecting M1 activity towards the skin EMG signal without trial averaging. We studied delta (1-4 Hz), theta (4-7 Hz), alpha (8-12 Hz), beta (15-30 Hz), and gamma (30-90 Hz) bands in 13 healthy participants (26 datasets) and two presurgical patients with sensorimotor dysfunction. In healthy participants, EMG-projected MEG accurately localized M1 with high accuracy in delta (100.0%), theta (100.0%), and beta (76.9%) bands, but not alpha (34.6%) and gamma (0.0%) bands. Except for delta, all other frequency bands were above the movement frequency and its harmonics. In both presurgical patients, M1 activity in the affected hemisphere was also accurately localized, despite highly irregular EMG movement patterns in one patient. Altogether, our EMG-projected MEG imaging approach is highly accurate and feasible for M1 mapping in presurgical patients. The results also provide insight into movement related brain-muscle coupling above the movement frequency and its harmonics.

Corticomuscular and intermuscular coherence are correlated after stroke: a simplified motor control?

During movement, corticomuscular coherence is a measure of central-peripheral communication, while intermuscular coherence is a measure of the amount of common central drive to the muscles. Although these two measures are modified in stroke subjects, no author has explored a correlation between them, neither in stroke subjects nor in healthy subjects. Twenty-four chronic stroke subjects and 22 healthy control subjects were included in this cohort study, and they performed 20 active elbow extension movements. The electroencephalographic and electromyographic activity of the elbow flexors and extensors were recorded. Corticomuscular and intermuscular coherence were calculated in the time-frequency domain for each limb of stroke and control subjects. Partial rank correlations were performed to study the link between these two variables. Our results showed a positive correlation between corticomuscular and intermuscular coherence only for stroke subjects, for their paretic and non-paretic limbs (P < 0.022; Rho > 0.50). These results suggest, beyond the cortical and spinal hypotheses to explain them, that stroke subjects present a form of simplification of motor control. When central-peripheral communication increases, it is less modulated and more common to the muscles involved in the active movement. This motor control simplification suggests a new way of understanding the plasticity of the neuromuscular system after stroke.

Combined effect of contraction type and intensity on corticomuscular coherence during isokinetic plantar flexions

During isometric contractions, corticomuscular coherence (CMC) may be modulated along with the contraction intensity. Furthermore, CMC may also vary between contraction types due to the contribution of spinal inhibitory mechanisms. However, the interaction between the effect of the contraction intensity and of the contraction type on CMC remains hitherto unknown. Therefore, CMC and spinal excitability modulations were compared during submaximal isometric, shortening and lengthening contractions of plantar flexor muscles at 25, 50, and 70% of the maximal soleus (SOL) EMG activity. CMC was computed in the time-frequency domain between the Cz EEG electrode signal and the SOL or medial gastrocnemius (MG) EMG signals. The results indicated that beta-band CMC was decreased in the SOL only between 25 and 50-70% contractions for both isometric and anisometric contractions, but remained similar for all contraction intensities in the MG. Spinal excitability was similar for all contraction intensities in both muscles. Meanwhile a divergence of the EEG and the EMG signals mean frequency was observed only in the SOL and only between 25 and 50-70% contractions, independently from the contraction type. Collectively, these findings confirm an effect of the contraction intensity on beta-band CMC, although it was only measured in the SOL, between low-level and high-level contraction intensities. Furthermore, the current findings provide new evidence that the observed modulations of beta-band CMC with the contraction intensity does not depend on the contraction type or on spinal excitability variations.

Temporally stable beta sensorimotor oscillations and cortico-muscular coupling underlie force steadiness

As humans, we seamlessly hold objects in our hands, and may even lose consciousness of these objects. This phenomenon raises the unsettled question of the involvement of the cerebral cortex, the core area for voluntary motor control, in dynamically maintaining steady muscle force. To address this issue, we measured magnetoencephalographic brain activity from healthy adults who maintained a steady pinch grip. Using a novel analysis approach, we uncovered fine-grained temporal modulations in the beta sensorimotor brain rhythm and its coupling with muscle activity, with respect to several aspects of muscle force (rate of increase/decrease or plateauing high/low). These modulations preceded changes in force features by ∼40 ms and possessed behavioral relevance, as less salient or absent modulation predicted a more stable force output. These findings have consequences for the existing theories regarding the functional role of cortico-muscular coupling, and suggest that steady muscle contractions are characterized by a stable rather than fluttering involvement of the sensorimotor cortex.

Design and validation of a wireless Body Sensor Network for integrated EEG and HD-sEMG acquisitions

Sensorimotor integration is the process through which the human brain plans the motor program execution according to external sources. Within this context, corticomuscular and corticokinematic coherence analyses are common methods to investigate the mechanism underlying the central control of muscle activation. This requires the synchronous acquisition of several physiological signals, including EEG and sEMG. Nevertheless, physical constraints of the current, mostly wired, technologies limit their application in dynamic and naturalistic contexts. In fact, although many efforts were made in the development of biomedical instrumentation for EEG and HD-sEMG signal acquisition, the need for an integrated wireless system is emerging. We hereby describe the design and validation of a new fully wireless body sensor network for the integrated acquisition of EEG and HD-sEMG signals. This Body Sensor Network is composed of wireless bio-signal acquisition modules, named sensor units, and a set of synchronization modules used as a general-purpose system for time-locked recordings. The system was characterized in terms of accuracy of the synchronization and quality of the collected signals. An in-depth characterization of the entire system and an end-to-end comparison of the wireless EEG sensor unit with a wired benchmark EEG device were performed. The proposed device represents an advancement of the State-of-the-Art technology allowing the integrated acquisition of EEG and HD-sEMG signals for the study of sensorimotor integration.

EEG hyperscanning in motor rehabilitation: a position paper

Studying the human brain during interpersonal interaction allows us to answer many questions related to motor control and cognition. For instance, what happens in the brain when two people walking side by side begin to change their gait and match cadences? Adapted from the neuroimaging techniques used in single-brain measurements, hyperscanning (HS) is a technique used to measure brain activity from two or more individuals simultaneously. Thus far, HS has primarily focused on healthy participants during social interactions in order to characterize inter-brain dynamics. Here, we advocate for expanding the use of this electroencephalography hyperscanning (EEG-HS) technique to rehabilitation paradigms in individuals with neurological diagnoses, namely stroke, spinal cord injury (SCI), Parkinson's disease (PD), and traumatic brain injury (TBI). We claim that EEG-HS in patient populations with impaired motor function is particularly relevant and could provide additional insight on neural dynamics, optimizing rehabilitation strategies for each individual patient. In addition, we discuss future technologies related to EEG-HS that could be developed for use in the clinic as well as technical limitations to be considered in these proposed settings.

Single channel surface electromyogram deconvolution to explore motor unit discharges

Interference surface electromyogram (EMG) reflects many bioelectric properties of active motor units (MU), which are however difficult to estimate due to the asynchronous summation of their discharges. This paper introduces a deconvolution technique to estimate the cumulative firings of MUs. Tests in simulations show that the power spectral density of the estimated MU firings has a low-frequency peak corresponding to the mean firing rate of MUs in the detection volume of the recording system, weighted by the amplitudes of MU action potentials. The peak increases in amplitude and its centroid shifts to a higher frequency when MU synchronization is simulated (mainly due to the shift of discharges of large MUs). The peak is found even at high force levels, when such a contribution does not emerge from the EMG. This result is also confirmed in preliminary applications to experimental data. Moreover, the simulated cumulative firings of MUs are estimated with a correlation above 90% (considering frequency contributions up to 150 Hz), for all force levels. The method requires a single EMG channel, thus being feasible even in applied studies using simple recording systems. It may open many potential applications, e.g., in the study of the modulation of MU firing rate induced by either fatigue or pathology and in coherency analysis. Graphical Abstract Examples of application of the deconvolution (Deconv) algorithm and comparison with the cumulative firings and the cumulated weighted firings (CWF, i.e., each firing pattern is weighted by the root mean squared amplitude of the corresponding MU action potential). Portions of data are shown on the left, the power spectral densities (PSD) on the right (Welch method applied to 3 s of data, sub-epochs of 0.5 s, mean value removed from each of them, 50% of overlap). A) Simulated signal (50% of maximal voluntary contraction, MVC) with random MU firings. B) Simulated signal (50% MVC) with a level of synchronization equal to 10%. C) Experimental data from vastus medialis at 40% MVC (data decomposed by the algorithm of Holobar and Zazula, IEEE Trans. Sig. Proc. 2007; PSD of the cumulated firings almost identical to that of CWF, as few MUs were identified).

Aging and Strength Training Influence Knee Extensor Intermuscular Coherence During Low- and High-Force Isometric Contractions

Aging is associated with reduced maximum force production and force steadiness during low-force tasks, but both can be improved by training. Intermuscular coherence measures coupling between two peripheral surface electromyography (EMG) signals in the frequency domain. It is thought to represent the presence of common input to alpha-motoneurons, but the functional meaning of intermuscular coherence, particularly regarding aging and training, remain unclear. This study investigated knee extensor intermuscular coherence in previously sedentary young (18–30 years) and older (67–73 years) subjects before and after a 14-week strength training intervention. YOUNG and OLDER groups performed maximum unilateral isometric knee extensions [100% maximum voluntary contraction (MVC)], as well as force steadiness tests at 20 and 70% MVC, pre- and post-training. Intermuscular (i.e., EMG-EMG) coherence analyses were performed for all (three) contraction intensities in vastus lateralis and medialis muscles. Pre-training coefficient of force variation (i.e., force steadiness) and MVC (i.e., maximum torque) were similar between groups. Both groups improved MVC through training, but YOUNG improved more than OLDER (42 ± 27 Nm versus 18 ± 16 Nm, P = 0.022). Force steadiness did not change during 20% MVC trials in either group, but YOUNG demonstrated increased coefficient of force variation during 70% MVC trials (1.28 ± 0.46 to 1.57 ± 0.70, P = 0.01). YOUNG demonstrated greater pre-training coherence during 20% and 70% MVC trials, particularly within the 8–14 Hz (e.g., 20%: 0.105 ± 0.119 versus 0.016 ± 0.009, P = 0.001) and 16–30 Hz (20%: 0.063 ± 0.078 versus 0.012 ± 0.007, P = 0.002) bands, but not during 100% MVC trials. Strength training led to increases in intermuscular coherence within the 40–60 Hz band during 70% MVC trials in YOUNG only, while OLDER decreased within the 8–14 Hz band during 100% MVC trials. Age-related differences in intermuscular coherence were observed between young and older individuals, even when neuromuscular performance levels were similar. The functional significance of intermuscular coherence remains unclear, since coherence within different frequency bands did not explain any of the variance in the regression models for maximum strength or force steadiness during 20 and 70% MVC trials.

Optimization of surface electrodes location for H-reflex recordings in soleus muscle

The Hoffmann reflex (H reflex) is extensively used to investigate the spinal motor neuron excitability in healthy and pathological subjects. Obtaining a representative and robust amplitude estimation of the H reflex is of marked relevance in clinical as well as in research applications. As for the motor responses, this issue is strictly related to the electrode positioning, especially for large, pinnate muscles such as the triceps surae. In this study we investigated the effect of electrode position on soleus H-reflex amplitude. A grid of 96 electrodes was used to identify maximal H reflexes (Hmax) across the whole soleus region available for surface recording. The spatial distribution of Hmax amplitude detected in monopolar and single-differential derivations was used to determine where greatest reflex responses were detected from soleus. For both derivations and for all participants, largest Hmax were detected consistently over the central soleus region, in correspondence of the muscle superficial aponeurosis. Indeed, the amplitude of Hmax provided by conventional electrodes (1cm² area, 2cm apart) located centrally was significantly greater (median: 35% for monopolar and 79% for single-differential derivations) than that obtained medially, where surface electromyograms are typically recorded from soleus. Computer simulations, used to assist in the interpretation of results, suggest the soleus pinnate architecture was the key determinant of the medio-lateral variability observed for the experimental Hmax. The presented results provide a clear indication for electrode positioning, of crucial relevance in applied studies aimed at eliciting H reflexes.

Individual difference in β-band corticomuscular coherence and its relation to force steadiness during isometric voluntary ankle dorsiflexion in healthy humans

OBJECTIVE

Magnitude of β-band coherent neural activities between the sensorimotor cortex and contracting muscle is known to vary across healthy individuals. To clarify how this variance affects actual motor function, this study examined associations between the corticomuscular coherence (CMC) and force steadiness.

METHODS

CMC was calculated between scalp electroencephalograms (EEGs) over the sensorimotor cortex and surface electromyograms (EMGs) from the tibialis anterior muscle during tonic isometric voluntary ankle dorsiflexion at 30% of maximal effort in 22 healthy individuals. We calculated the maximal peak of CMC (CMCmax), and examined its relations to some measures of force fluctuation, such as the coefficient of variation (ForceCV), the sum of the power spectral density within 1-4Hz (Forceδ-PSD), 5-14Hz (Forceα-PSD), and 15-35Hz (Forceβ-PSD) bands of force signal.

RESULTS

In all participants showing significant CMC, CMCmax was observed within the β-band. CMCmax was varied across participants (range, 0.084-0.451), and was correlated significantly and positively with ForceCV (r=0.602, p=0.003), Forceβ-PSD (r=0.637, p=0.001), Forceα-PSD (r=0.647, p=0.001), and Forceδ-PSD (r=0.518, p=0.014).

CONCLUSION

The magnitude of the CMC between EEG over the sensorimotor cortex and EMG of contracting muscle is associated with the amount of force fluctuation during tonic isometric voluntary ankle dorsiflexion in healthy humans.

SIGNIFICANCE

CMC may influence an individual's ability to stabilize their muscle force output.

The network sustaining action myoclonus: a MEG-EMG study in patients with EPM1

Background

To explore the cortical network sustaining action myoclonus and to found markers of the resulting functional impairment, we evaluated the distribution of the cortico-muscular coherence (CMC) and the frequency of coherent cortical oscillations with magnetoencephalography (MEG). All patients had EPM1 (Unverricht-Lundborg) disease known to present with prominent and disabling movement-activated myoclonus.

Methods

Using autoregressive models, we evaluated CMC on MEG sensors grouped in regions of interests (ROIs) above the main cortical areas. The movement was a repeated sustained isometric extension of the right hand and right foot. We compared the data obtained in 10 EPM1 patients with those obtained in 10 age-matched controls.

Results

As expected, CMC in beta band was significantly higher in EPM1 patients compared to controls in the ROIs exploring the sensorimotor cortex, but, it was also significantly higher in adjacent ROIs ipsilateral and contralateral to the activated limb. Moreover, the beta-CMC peak occurred at frequencies significantly slower and more stable frequencies in EPM1 patients with respect to controls. The frequency of the beta-CMC peak inversely correlated with the severity of myoclonus.

Conclusions

the high and spatially extended beta-CMC peaking in a restricted range of low-beta frequencies in EPM1 patients, suggest that action myoclonus may result not only from an enhanced local synchronization but also from a specific oscillatory activity involving an expanded neuronal pool. The significant relationship between beta-CMC peak frequency and the severity of the motor impairment can represent a useful neurophysiological marker for the patients’ evaluation and follow-up.