Adaptive real-time identification of motor unit discharges from non-stationary high-density surface electromyographic signals

Objective. Estimation of the discharge pattern of motor units by electromyography (EMG) decomposition has been applied for neurophysiologic investigations, clinical diagnosis, and human-machine interfacing. However, most of the methods for EMG decomposition are currently applied offline. Here, we propose an approach for high-density surface EMG decomposition in real-time. Methods. A real-time decomposition scheme including two sessions, offline training and online decomposition, is proposed based on the convolutional kernel compensation algorithm. The estimation parameters, separation vectors and the thresholds for spike extraction, are first computed during offline training, and then they are directly applied to estimate motor unit spike trains (MUSTs) during the online decomposition. The estimation parameters are updated with the identification of new discharges to adapt to non-stationary conditions. The decomposition accuracy was validated on simulated EMG signals by convolving synthetic MUSTs with motor unit action potentials (MUAPs). Moreover, the accuracy of the online decomposition was assessed from experimental signals recorded from forearm muscles using a signal-based performance metrics (pulse-to-noise ratio, PNR). Main results. The proposed algorithm yielded a high decomposition accuracy and robustness to non-stationary conditions. The accuracy of MUSTs identified from simulated EMG signals was > 80% for most conditions. From experimental EMG signals, on average, 12±2 MUSTs were identified from each electrode grid with PNR of 25.0±1.8 dB, corresponding to an estimated decomposition accuracy > 75%. Conclusion and Significance. These results indicate the feasibility of real-time identification of motor unit activities non-invasively during variable force contractions, extending the potential applications of high-density EMG as a neural interface

Semi-Automated Identification of Motor Units Concurrently Recorded in High-Density Surface and Intramuscular Electromyography

: An increasing focus on extending automated surface electromyography (EMG) decomposition algorithms to operate under non-stationary conditions requires rigorous and robust validation. However, relevant benchmarks derived manually from iEMG are laborsome to obtain and this is further exacerbated by the need to consider multiple contraction conditions. This work demonstrates a semi-automatic technique for extracting motor units (MUs) whose activities are present in concurrently recorded high-density surface EMG (HD-sEMG) and intramuscular EMG (iEMG) during isometric contractions. We leverage existing automatic surface decomposition algorithms for initial identification of active MUs. Resulting spike times are then used to identify (trigger) the sources that are concurrently detectable in iEMG. We demonstrate this technique on recordings targeting the extensor carpi radialis brevis in five human subjects. This dataset consists of 117 trials across different force levels and wrist angles, from which the presented method yielded a set of 367 high-confidence decompositions. Thus, our approach effectively alleviates the overhead of manual decomposition as it efficiently produces reliable benchmarks under different conditions.Clinical Relevance- We present an efficient method for obtaining high-quality in-vivo decomposition particularly useful in the verification of new surface decomposition approaches

Identification of Motor Unit Twitch Properties in the Intact Human In Vivo

Restoring natural motor function in neurologically injured individuals is challenging, largely due to the lack of personalization in current neurorehabilitation technologies. Signal-driven neuro-musculoskeletal models may offer a novel paradigm for devising novel closed-loop rehabilitation strategies according to an individual's physiology. However, current modelling techniques are constrained to bipolar electromyography (EMG), thereby lacking the resolution necessary to extract the activity of individual motor units (MUs) in vivo. In this work, we decoded MU spike trains from high-density (HD)-EMG to obtain relevant neural properties across multiple isometric plantar-dorsiflexion tasks. Then, we sampled MU statistical distributions and used them to reproduce MU specific activation profiles. Results showed bimodal distributions which may correspond to slow and fast MU populations. The estimated activation profiles showed a high degree of similarity to the reference torque (R2>0.8) across the recorded muscles. This suggests that the estimation of MU twitch properties is a crucial step for the translation of neural information into muscle force.Clinical Relevance- This work has multiple implications for understanding the underlying mechanism of motor impairment and for developing closed-loop strategies for modulating alpha motor circuitries in neurologically injured individuals

HYDRA-HGR: A Hybrid Transformer-based Architecture for Fusion of Macroscopic and Microscopic Neural Drive Information

Development of advance surface Electromyogram (sEMG)-based Human-Machine Interface (HMI) systems is of paramount importance to pave the way towards emergence of futuristic Cyber-Physical-Human (CPH) worlds. In this context, the main focus of recent literature was on development of different Deep Neural Network (DNN)-based architectures that perform Hand Gesture Recognition (HGR) at a macroscopic level (i.e., directly from sEMG signals). At the same time, advancements in acquisition of High-Density sEMG signals (HD-sEMG) have resulted in a surge of significant interest on sEMG decomposition techniques to extract microscopic neural drive information. However, due to complexities of sEMG decomposition and added computational overhead, HGR at microscopic level is less explored than its aforementioned DNN-based counterparts. In this regard, we propose the HYDRA-HGR framework, which is a hybrid model that simultaneously extracts a set of temporal and spatial features through its two independent Vision Transformer (ViT)-based parallel architectures (the so called Macro and Micro paths). The Macro Path is trained directly on the pre-processed HD-sEMG signals, while the Micro path is fed with the p-to-p values of the extracted Motor Unit Action Potentials (MUAPs) of each source. Extracted features at macroscopic and microscopic levels are then coupled via a Fully Connected (FC) fusion layer. We evaluate the proposed hybrid HYDRA-HGR framework through a recently released HD-sEMG dataset, and show that it significantly outperforms its stand-alone counterparts. The proposed HYDRA-HGR framework achieves average accuracy of 94.86% for the 250 ms window size, which is 5.52% and 8.22% higher than that of the Macro and Micro paths, respectively

Multichannel surface EMG decomposition based on measurement correlation and LMMSE

A method based on measurement correlation (MC) and linear minimum mean square error (LMMSE) for multichannel surface electromyography (sEMG) signal decomposition was developed in this study. This MC-LMMSE method gradually and iteratively increases the correlation between an optimized vector and a reconstructed matrix that is correlated with the measurement matrix. The performance of the proposed MC-LMMSE method was evaluated with both simulated and experimental sEMG signals. Simulation results show that the MC-LMMSE method can successfully reconstruct up to 53 innervation pulse trains with a true positive rate greater than 95%. The performance of the MC-LMMSE method was also evaluated using experimental sEMG signals collected with a 64-channel electrode array from the first dorsal interosseous muscles of three subjects at different contraction levels. A maximum of 16 motor units were successfully extracted from these multichannel experimental sEMG signals. The performance of the MC-LMMSE method was further evaluated with multichannel experimental sEMG data by using the “two sources” method. The large population of common MUs extracted from the two independent subgroups of sEMG signals demonstrates the reliability of the MC-LMMSE method in multichannel sEMG decomposition

Deep learning for robust decomposition of high-density surface EMG signals

Blind source separation (BSS) algorithms, such as gradient convolution kernel compensation (gCKC), can efficiently and accurately decompose high-density surface electromyography (HD-sEMG) signals into constituent motor unit (MU) action potential trains. Once the separation matrix is blindly estimated on a signal interval, it is also possible to apply the same matrix to subsequent signal segments. Nonetheless, the trained separation matrices are sub-optimal in noisy conditions and require that incoming data undergo computationally expensive whitening. One unexplored alternative is to instead use the paired HD-sEMG signal and BSS output to train a model to predict MU activations within a supervised learning framework. A gated recurrent unit (GRU) network was trained to decompose both simulated and experimental unwhitened HD-sEMG signal using the output of the gCKC algorithm. The results on the experimental data were validated by comparison with the decomposition of concurrently recorded intramuscular EMG signals. The GRU network outperformed gCKC at low signal-to-noise ratios, proving superior performance in generalising to new data. Using 12 seconds of experimental data per recording, the GRU performed similarly to gCKC, at rates of agreement of 92.5% (84.5%-97.5%) and 94.9% (88.8%-100.0%) respectively for GRU and gCKC against matched intramuscular sources

A real-time and convex model for the estimation of muscle force from surface electromyographic signals in the upper and lower limbs

Surface electromyography (sEMG) is a signal consisting of different motor unit action potential trains and records from the surface of the muscles. One of the applications of sEMG is the estimation of muscle force. We proposed a new real-time convex and interpretable model for solving the sEMG-force estimation. We validated it on the upper limb during isometric voluntary flexions-extensions at 30%, 50%, and 70% Maximum Voluntary Contraction in five subjects, and lower limbs during standing tasks in thirty-three volunteers, without a history of neuromuscular disorders. Moreover, the performance of the proposed method was statistically compared with that of the state-of-the-art (13 methods, including linear-in-the-parameter models, Artificial Neural Networks and Supported Vector Machines, and non-linear models). The envelope of the sEMG signals was estimated, and the representative envelope of each muscle was used in our analysis. The convex form of an exponential EMG-force model was derived, and each muscle's coefficient was estimated using the Least Square method. The goodness-of-fit indices, the residual signal analysis (bias and Bland-Altman plot), and the running time analysis were provided. For the entire model, 30% of the data was used for estimation, while the remaining 20% and 50% were used for validation and testing, respectively. The average R-square (%) of the proposed method was 96.77 +/- 1.67 [94.38, 98.06] for the test sets of the upper limb and 91.08 +/- 6.84 [62.22, 96.62] for the lower-limb dataset (MEAN +/- SD [min, max]). The proposed method was not significantly different from the recorded force signal (p-value = 0.610); that was not the case for the other tested models. The proposed method significantly outperformed the other methods (adj. p-value < 0.05). The average running time of each 250 ms signal of the training and testing of the proposed method was 25.7 +/- 4.0 [22.3, 40.8] and 11.0 +/- 2.9 [4.7, 17.8] in microseconds for the entire dataset. The proposed convex model is thus a promising method for estimating the force from the joints of the upper and lower limbs, with applications in load sharing, robotics, rehabilitation, and prosthesis control for the upper and lower limbs

A real-time and convex model for the estimation of muscle force from surface electromyographic signals in the upper and lower limbs

Surface electromyography (sEMG) is a signal consisting of different motor unit action potential trains and records from the surface of the muscles. One of the applications of sEMG is the estimation of muscle force. We proposed a new real-time convex and interpretable model for solving the sEMG—force estimation. We validated it on the upper limb during isometric voluntary flexions-extensions at 30%, 50%, and 70% Maximum Voluntary Contraction in five subjects, and lower limbs during standing tasks in thirty-three volunteers, without a history of neuromuscular disorders. Moreover, the performance of the proposed method was statistically compared with that of the state-of-the-art (13 methods, including linear-in-the-parameter models, Artificial Neural Networks and Supported Vector Machines, and non-linear models). The envelope of the sEMG signals was estimated, and the representative envelope of each muscle was used in our analysis. The convex form of an exponential EMG-force model was derived, and each muscle’s coefficient was estimated using the Least Square method. The goodness-of-fit indices, the residual signal analysis (bias and Bland-Altman plot), and the running time analysis were provided. For the entire model, 30% of the data was used for estimation, while the remaining 20% and 50% were used for validation and testing, respectively. The average R-square (%) of the proposed method was 96.77 ± 1.67 [94.38, 98.06] for the test sets of the upper limb and 91.08 ± 6.84 [62.22, 96.62] for the lower-limb dataset (MEAN ± SD [min, max]). The proposed method was not significantly different from the recorded force signal (p-value = 0.610); that was not the case for the other tested models. The proposed method significantly outperformed the other methods (adj. p-value &lt; 0.05). The average running time of each 250 ms signal of the training and testing of the proposed method was 25.7 ± 4.0 [22.3, 40.8] and 11.0 ± 2.9 [4.7, 17.8] in microseconds for the entire dataset. The proposed convex model is thus a promising method for estimating the force from the joints of the upper and lower limbs, with applications in load sharing, robotics, rehabilitation, and prosthesis control for the upper and lower limbs

Effect of Resistance Training and Fish Protein Intake on Motor Unit Firing Pattern and Motor Function of Elderly

We investigated the effect of resistance training and fish protein intake on the motor unit firing pattern and motor function in elderly. Fifty healthy elderly males and females (69.2 ± 4.7 years) underwent 6 weeks of intervention. We applied the leg-press exercise as resistance training and fish protein including Alaska pollack protein (APP) as nutritional supplementation. Subjects were divided into four groups: fish protein intake without resistance training (APP-CN, n = 13), placebo intake without resistance training (PLA-CN, n = 12), fish protein intake with resistance training (APP-RT, n = 12), and placebo intake with resistance training (PLA-RT, n = 13). Motor unit firing rates were calculated from multi-channel surface electromyography by the Convolution Kernel. For the chair-stand test, while significant increases were observed at 6 weeks compared with 0 week in all groups (p &lt; 0.05), significant increases from 0 to 3 weeks and 6 weeks were observed in APP-RT (18.2 ± 1.9 at 0 week to 19.8 ± 2.2 at 3 weeks and 21.2 ± 1.9 at 6 weeks) (p &lt; 0.05). Increase and/or decrease in the motor unit firing rate were mainly noted within motor units with a low-recruitment threshold in APP-RT and PLA-RT at 3 and 6 weeks (12.3 pps at 0 week to 13.6 pps at 3 weeks and 12.1 pps at 6 weeks for APP-RT and 12.9 pps at 0 week to 13.9 pps at 3 weeks and 14.1 pps at 6 weeks for PLA-RT at 50% of MVC) (p &lt; 0.05), but not in APP-CN or PLA-CN (p &gt; 0.05). Time courses of changes in the results of the chair-stand test and motor unit firing rate were different between APP-RT and PLA-RT. These findings suggest that, in the elderly, the effect of resistance training on the motor unit firing rate is observed in motor units with a low-recruitment threshold, and additional fish protein intake modifies these adaptations in motor unit firing patterns and the motor function following resistance training

Tremor severity in Parkinson’s disease and cortical changes of areas controlling movement sequencing: a preliminary study.

. There remains much to learn about the changes in cortical anatomy that are associated with tremor severity in Parkinson’s disease (PD). For this reason, we used a combination of structural neuroimaging to measure cortical thickness and neurophysiological studies to analyze whether PD tremor was associated with cortex integrity. Magnetic resonance imaging and neurophysiological assessment were performed in 13 nondemented PD patients (9 women, 69.2%) with a clearly tremor-dominant phenotype. Cortical reconstruction and volumetric segmentation was performed with the Freesurfer image analysis software. Assessment of tremor was performed by means of high-density surface electromyography (hdEMG) and inertial measurement units (IMUs). Individual motor unit discharge patterns were identified from surface hdEMG and tremor metrics quantifying motor unit synchronization from IMUs were defined. Increased motor unit synchronization (i.e., more severe tremor) was associated with cortical changes (i.e., atrophy) in dorsal premotor cortices, left posterior parietal cortex, left lateral orbitofrontal cortex, cingulate cortex bilaterally, left posterior and transverse temporal cortex, and left occipital lobe, as well as reduced left middle temporal volume. Given that the majority of these areas are involved in controlling movement sequencing, our results support Albert’s classic hypothesis that PD tremor may be the result of an involuntary activation of a program of motor behavior used in the genesis of rapid voluntary alternating movements.pre-print670 K