Motor units are the smallest functional units of our movements. The study of their activation provides a window into the mechanisms of neural control of movement in humans. The classic methods for ...motor unit investigations date to several decades ago. They are based on invasive recordings with selective needle or wire electrodes. Conversely, the noninvasive (surface) EMG has been commonly processed as an interference signal, with the extraction of its global characteristics, e.g., amplitude. These characteristics, however, are only crudely associated to the underlying motor unit activities. In the last decade, methods have been proposed for reliably extracting individual motor unit activities from the interference surface EMG signal. We describe these methods in this review, with a focus on blind source separation (BSS) and techniques used on decomposed EMG signals. For example, from the motor unit discharge timings, information can be extracted regarding the synaptic input received by the corresponding motor neurons. In reviewing these methods, we also provide examples of applications in representative conditions, such as pathological tremor. In conclusion, we provide an overview of processing methods of the surface EMG signal that allow a reliable characterization of individual motor units in vivo in humans.
Quercetin is a polyphenolic flavonoid that has reported to block the binding of adenosine to A1 receptors at central nervous system and increase calcium release from the sarcoplasmic reticulum at ...skeletal muscle. The aim of the present study was to investigate the acute effect of quercetin ingestion on motor unit activation and muscle contractile properties. High-density surface electromyography during submaximal contractions and electrically elicited contraction torque in knee extensor muscles were measured before (PRE) and 60 min after (POST) quercetin glycosides or placebo ingestions in 13 young males. Individual motor units of the vastus lateralis muscle were identified from high-density surface electromyography by the Convolution Kernel Compensation technique. Firing rates of motor units recruited at 30–50% of the maximal voluntary contraction torque (MVC) were increased from PRE to POST only with quercetin (9.0 ± 2.3 to 10.5 ± 2.0 pps,
p
= 0.034). Twitch torque during doublet stimulation was decreased from PRE to POST with placebo (77.1 ± 17.1 to 73.9 ± 17.6 Nm,
p
= 0.005), but not with quercetin (
p
> 0.05). For motor units recruited at < 10% of MVC, normalized firing rate were decreased with quercetin (1.52 ± 0.33 to 1.58 ± 0.35%MVC/pps,
p
= 0.002) but increased with placebo (1.61 ± 0.32 to 1.57 ± 0.31%MVC/pps,
p
= 0.005). These results suggest that ingested quercetin has the functional roles to: mitigate reduction in the muscle contractile properties, enhance activations of relatively higher recruitment threshold motor units, and inhibit activation of relatively lower recruitment threshold motor units.
Despite not recording directly from neural cells, the surface electromyogram (EMG) signal contains information on the neural drive to muscles, i.e, the spike trains of motor neurons. Using this ...property, myoelectric control consists of the recording of EMG signals for extracting control signals to command external devices, such as hand prostheses. In commercial control systems, the intensity of muscle activity is extracted from the EMG and used for single degrees of freedom activation (direct control). Over the past 60 years, academic research has progressed to more sophisticated approaches but, surprisingly, none of these academic achievements has been implemented in commercial systems so far. We provide an overview of both commercial and academic myoelectric control systems and we analyze their performance with respect to the characteristics of the ideal myocontroller. Classic and relatively novel academic methods are described, including techniques for simultaneous and proportional control of multiple degrees of freedom and the use of individual motor neuron spike trains for direct control. The conclusion is that the gap between industry and academia is due to the relatively small functional improvement in daily situations that academic systems offer, despite the promising laboratory results, at the expense of a substantial reduction in robustness. None of the systems so far proposed in the literature fulfills all the important criteria needed for widespread acceptance by the patients, i.e. intuitive, closed-loop, adaptive, and robust real-time (<;200 ms delay) control, minimal number of recording electrodes with low sensitivity to repositioning, minimal training, limited complexity and low consumption. Nonetheless, in recent years, important efforts have been invested in matching these criteria, with relevant steps forwards.
This paper studies a novel decomposition technique, suitable for blind separation of linear mixtures of signals comprising finite-length symbols. The observed symbols are first modeled as channel ...responses in a multiple-input-multiple-output (MIMO) model, while the channel inputs are conceptually considered sparse positive pulse trains carrying the information about the symbol arising times. Our decomposition approach compensates channel responses and aims at reconstructing the input pulse trains directly. The algorithm is derived first for the overdetermined noiseless MIMO case. A generalized scheme is then provided for the underdetermined mixtures in noisy environments. Although blind, the proposed technique approaches Bayesian optimal linear minimum mean square error estimator and is, hence, significantly noise resistant. The results of simulation tests prove it can be applied to considerably underdetermined convolutive mixtures and even to the mixtures of moderately correlated input pulse trains, with their cross-correlation up to 10% of its maximum possible value.
We describe the method for identification of motor unit (MU) firings from high-density surface electromyograms (hdEMG), recorded during repeated dynamic muscle contractions. A new convolutive data ...model for dynamic hdEMG is presented, along with the pulse-to-noise ratio (PNR) metric for assessment of MU identification accuracy and analysis of the impact of MU action potential (MUAP) changes in dynamic muscle contractions on MU identification. We tested the presented methodology on signals from biceps brachii, vastus lateralis, and rectus famoris muscles, all during different speeds of dynamic contractions. In synthetic signals with excitation levels of 10%, 30% and 50%, and MUAPs experimentally recorded from biceps brachii muscle, the presented method identified 15 ± 1, 18 ± 1, and 20 ± 1 MUs per contraction, respectively, all with average sensitivity and precision >90% and PNR >30dB. In experimental signals acquired during low force contractions of vastus lateralis and rectus femoris muscle, the method identified 9.4±1.9 and 7.8±1.4 MUs with PNR values of 35.4±3.6 and 34.1±2.7 dB. In comparison with the previously introduced Convolution Kernel Compensation method, the capability of the new method to follow dynamic MUAP changes is confirmed, also in relatively fast muscle contractions.
Key points
We propose and validate a method for accurately identifying the activity of populations of motor neurons during contractions at maximal rate of force development in humans.
The behaviour ...of the motor neuron pool during rapid voluntary contractions in humans is presented.
We show with this approach that the motor neuron recruitment speed and maximal motor unit discharge rate largely explains the individual ability in generating rapid force contractions.
The results also indicate that the synaptic inputs received by the motor neurons before force is generated dictate human potential to generate force rapidly.
This is the first characterization of the discharge behaviour of a representative sample of human motor neurons during rapid contractions.
During rapid contractions, motor neurons are recruited in a short burst and begin to discharge at high frequencies (up to >200 Hz). In the present study, we investigated the behaviour of relatively large populations of motor neurons during rapid (explosive) contractions in humans, applying a new approach to accurately identify motor neuron activity simultaneous to measuring the rate of force development. The activity of spinal motor neurons was assessed by high‐density electromyographic decomposition from the tibialis anterior muscle of 20 men during isometric explosive contractions. The speed of motor neuron recruitment and the instantaneous motor unit discharge rate were analysed as a function of the impulse (the time–force integral) and the maximal rate of force development. The peak of motor unit discharge rate occurred before force generation and discharge rates decreased thereafter. The maximal motor unit discharge rate was associated with the explosive force variables, at the whole population level (r2 = 0.71 ± 0.12; P < 0.001). Moreover, the peak motor unit discharge and maximal rate of force variables were correlated with an estimate of the supraspinal drive, which was measured as the speed of motor unit recruitment before the generation of afferent feedback (P < 0.05). We show for the first time the full association between the effective neural drive to the muscle and human maximal rate of force development. The results obtained in the present study indicate that the variability in the maximal contractile explosive force of the human tibialis anterior muscle is determined by the neural activation preceding force generation.
Key points
We propose and validate a method for accurately identifying the activity of populations of motor neurons during contractions at maximal rate of force development in humans.
The behaviour of the motor neuron pool during rapid voluntary contractions in humans is presented.
We show with this approach that the motor neuron recruitment speed and maximal motor unit discharge rate largely explains the individual ability in generating rapid force contractions.
The results also indicate that the synaptic inputs received by the motor neurons before force is generated dictate human potential to generate force rapidly.
This is the first characterization of the discharge behaviour of a representative sample of human motor neurons during rapid contractions.
Abstract This brief review discusses the methods used to estimate the neural drive to muscles from the surface electromyogram (EMG). Surface EMG has been classically used to infer the neural ...activation of muscle by associating its amplitude with the number of action potentials discharged by a population of motor neurons. Although this approach is valuable in some applications, the amplitude of the surface EMG is only a crude indicator of the neural drive to muscle. More advanced methods are now available to estimate the neural drive to muscle from the surface EMG. These approaches identify the discharge times of a few motor units by decomposing the EMG signal to determine the relative changes in neural activation. This approach is reliable in several conditions and muscles for isometric contractions of moderate force, but is limited to the few superficial units that can be identified in the recordings.
We introduce two novel methods for the estimation of muscle excitation from surface electromyograms (EMGs), the so called cumulative motor unit activity index (CAI) and robust CAI (rCAI). Both ...methods aim to remove the detected motor unit action potential (MUAP) contributions from EMG but do not assess the individual motor unit spike trains. Instead, they directly estimate the cumulative motor unit spike train (CST). We compared the methods with the spatially averaged root-mean-square (RMS) envelope of the EMG signals and with the CST, estimated by the previously introduced convolution kernel compensation (CKC) method. The tests on synthetic EMG with known muscle excitation profiles demonstrated superior accuracy of newly introduced methods. In the case of 64 EMG channels and 20-dB noise, the RMS, CAI, rCAI, and CKC estimators, calculated on 0.125-s-long signal epochs, yielded the normalized RMS error (NRMSE) of 14.5% ± 2.8%, 4.4% ± 3.2%, 4.1% ± 1.8%, and 6.3% ± 4.6%, respectively. In the experimental signals from wrist extensors and flexors, the RMS, CAI, rCAI, and CKC estimations were compared to exerted muscle force. When calculated on 0.125-s-long signal epochs, they yielded the NRMSE of 11.2% ± 3.5%, 8% ± 5.6%, 10.7% ± 6.8%, and 9.0% ± 4.9%, respectively. Therefore, the newly introduced methods exhibit accuracy that is comparable to at least 200-times slower CKC method.
This study aimed to determine whether neural drive is redistributed between muscles during a fatiguing isometric contraction, and if so, whether the initial level of common synaptic input between ...these muscles constrains this redistribution. We studied two muscle groups: triceps surae (14 participants) and quadriceps (15 participants). Participants performed a series of submaximal isometric contractions and a torque-matched contraction maintained until task failure. We used high-density surface electromyography to identify the behavior of 1,874 motor units from the soleus, gastrocnemius medialis (GM), gastrocnemius lateralis (GL), rectus femoris, vastus lateralis (VL), and vastus medialis (VM). We assessed the level of common drive between muscles in the absence of fatigue using a coherence analysis. We also assessed the redistribution of neural drive between muscles during the fatiguing contraction through the correlation between their cumulative spike trains (index of neural drive). The level of common drive between VL and VM was significantly higher than that observed for the other muscle pairs, including GL-GM. The level of common drive increased during the fatiguing contraction, but the differences between muscle pairs persisted. We also observed a strong positive correlation of neural drive between VL and VM during the fatiguing contraction (
= 0.82). This was not observed for the other muscle pairs, including GL-GM, which exhibited differential changes in neural drive. These results suggest that less common synaptic input between muscles allows for more flexible coordination strategies during a fatiguing task, i.e., differential changes in neural drive across muscles. The role of this flexibility on performance remains to be elucidated.
Redundancy of the neuromuscular system theoretically allows for a redistribution of the neural drive across muscles (i.e., between-muscle compensation) during a fatiguing contraction. Our results suggest that a high level of common input between muscles (e.g., vastus lateralis and medialis) represents a neural constraint making it less likely to redistribute the neural drive across these muscles. In this way, redistribution was only observed across muscles that share little common synaptic input (e.g., gastrocnemius lateralis and medialis).