Automatic Implementation of Progressive FastICA Peel-Off for High Density Surface EMG Decomposition

Myoelectric Fatigue and Motor-Unit Firing Patterns During Sinusoidal Vibration Superimposed on Low-Intensity Isometric Contraction

Vibration exercise (VE) has shown promising results for improving muscle strength and power performance when superimposed on high-level muscle contraction. However, low-level contraction may be more preferable in many rehabilitation programs due to the weakness of the patients. Unfortunately, the effects and underlying physiological mechanisms of VE superimposed on low-level contraction are unclear. This study aims to investigate the fatiguing effects and motor unit (MU) firing patterns during VE with low-level muscle contraction. Twenty-one healthy participants performed 60-s isometric contraction of the upper limb under a baseline force at 30% maximum voluntary contraction and superimposed vibration with an amplitude of 50% baseline and different frequencies at 0 Hz (control), 15, 25, 35, and 45 Hz. High-density surface electromyography (EMG) was recorded on the biceps brachii. The decay in muscle fiber conduction velocity, calculated in 3-s sliding windows, was employed as an indicator of myoelectric fatigue. MU firing patterns were obtained by decomposing the high-density EMG into MU spike trains. VE, particularly at 25 Hz, produces increased myoelectric fatigue as compared to the control condition. Besides, synchronized MU discharges are observed at the vibration frequency for 15- and 25-Hz VE and the sub-harmonics for 35- and 45-Hz VE. Furthermore, VE-induced increase in MU synchronization (as compared to control) seems to decrease with myoelectric fatigue. Significance: Our findings suggest that VE may be a suitable modality for rehabilitation programs, providing useful insights for subscribing appropriate VE training protocols.

A robust myoelectric pattern recognition framework based on individual motor unit activities against electrode array shifts

BACKGROUND AND OBJECTIVE

Electrode shift is always one of the critical factors to compromise the performance of myoelectric pattern recognition (MPR) based on surface electromyogram (SEMG). However, current studies focused on the global features of SEMG signals to mitigate this issue but it is just an oversimplified description of the human movements without incorporating microscopic neural drive information. The objective of this work is to develop a novel method for calibrating the electrode array shifts toward achieving robust MPR, leveraging individual motor unit (MU) activities obtained through advanced SEMG decomposition.

METHODS

All of the MUs from decomposition of SEMG data recorded at the original electrode array position were first initialized to train a neural network for pattern recognition. A part of decomposed MUs could be tracked and paired with MUs obtained at the original position based on spatial distribution of their MUAP waveforms, so as to determine the shift vector (describing both the orientation and distance of the shift) implicated consistently by these multiple MU pairs. Given the known shift vector, the features of the after-shift decomposed MUs were corrected accordingly and then fed into the network to finalize the MPR task. The performance of the proposed method was evaluated with data recorded by a 16 × 8 electrode array placed over the finger extensor muscles of 8 subjects performing 10 finger movement patterns.

RESULTS

The proposed method achieved a shift detection accuracy of 100 % and a pattern recognition accuracy approximating to 100 %, significantly outperforming the conventional methods with lower shift detection accuracies and lower pattern recognition accuracies (p < 0.05).

CONCLUSIONS

Our method demonstrated the feasibility of using decomposed MUAP waveforms' spatial distributions to calibrate electrode shift. This study provides a new tool to enhance the robustness of myoelectric control systems via microscopic neural drive information at an individual MU level.

2CFastICA: A Novel Method for High Density Surface EMG Decomposition Based on Kernel Constrained FastICA and Correlation Constrained FastICA

This study presents a novel high density surface electromyography (EMG) decomposition method, named as 2CFastICA, because it incorporates two key algorithms: kernel constrained FastICA and correlation constrained FastICA. The former focuses on overcoming the local convergence of FastICA without requiring the peel-off strategy used in the progressive FastICA peel-off (PFP) framework. The latter further refines the output of kernel constrained FastICA by correcting possible erroneous or missed spikes. The two constrained FastICA algorithms supplement each other to warrant the decomposition performance. The 2CFastICA method was validated using simulated surface EMG signals with different motor unit numbers and signal to noise ratios (SNRs). Two source validation was also performed by simultaneous high density surface EMG and intramuscular EMG recordings, showing a matching rate (MR) of (97.2 ± 3.5)% for 170 common motor units. In addition, a different form of two source validation was also conducted taking advantages of the high density surface EMG characteristics of patients with amyotrophic lateral sclerosis, showing a MR of (99.4 ± 0.9)% for 34 common motor units from interference and sparse datasets. Both simulation and experimental results indicate that 2CFastICA can achieve similar decomposition performance to PFP. However, the efficiency of decomposition can be greatly improved by 2CFastICA since the complex signal processing procedures associated with the peel-off strategy are not required any more. Along with this paper, we also provide the MATLAB open source code of 2CFastICA for high density surface EMG decomposition.

Decomposition strategy for surface EMG with few channels: a simulation study

Objective.In the specific use of electromyogram (EMG) driven prosthetics, the user's disability reduces the space available for the electrode array. We propose a framework for EMG decomposition adapted to the condition of a few channels (less than 30 observations), which can elevate the potential of prosthetics in terms of cost and applicability.Approach.The new framework contains a peel-off approach, a refining strategy for motor unit (MU) spike train and MU action potential and a re-subtracting strategy to adapt the framework to few channels environments. Simulated EMG signals were generated to test the framework. In addition, we quantify and analyze the effect of strategies used in the framework.Main results.The results show that the new algorithm has an average improvement of 19.97% in the number of MUs identified compared to the control algorithm. Quantitative analysis of the usage strategies shows that the re-subtracting and refining strategies can effectively improve the performance of the framework under the condition of few channels.Significance.These prove that the new framework can be applied to few channel conditions, providing a optimization space for neural interface design in cost and user adaptation.

Online prediction of sustained muscle force from individual motor unit activities using adaptive surface EMG decomposition

Decoding movement intentions from motor unit (MU) activities to represent neural drive information plays a central role in establishing neural interfaces, but there remains a great challenge for obtaining precise MU activities during sustained muscle contractions. In this paper, we presented an online muscle force prediction method driven by individual MU activities that were decomposed from prolonged surface electromyogram (SEMG) signals in real time. In the training stage of the proposed method, a set of separation vectors was initialized for decomposing MU activities. After transferring each decomposed MU activity into a twitch force train according to its action potential waveform, a neural network was designed and trained for predicting muscle force. In the subsequent online stage, a practical double-thread-parallel algorithm was developed. One frontend thread predicted the muscle force in real time utilizing the trained network and the other backend thread simultaneously updated the separation vectors. To assess the performance of the proposed method, SEMG signals were recorded from the abductor pollicis brevis muscles of eight subjects and the contraction force was simultaneously collected. With the update procedure in the backend thread, the force prediction performance of the proposed method was significantly improved in terms of lower root mean square deviation (RMSD) of around 10% and higher fitness (R2) of around 0.90, outperforming two conventional methods. This study provides a promising technique for real-time myoelectric applications in movement control and health.

Enhanced Dynamic Surface EMG Decomposition Using the Non-Negative Matrix Factorization and Three-Dimensional Motor Unit Localization

OBJECTIVE

Surface electromyography (sEMG) signal decomposition is of great importance in examining neuromuscular diseases and neuromuscular research, especially dynamic sEMG decomposition is even more technically challenging.

METHODS

A novel two-step sEMG decomposition approach was developed. The linear minimum mean square error estimation was first employed to extract estimated firing trains (EFTs) from the eigenvector matrices constructed using the non-negative matrix factorization (NMF). The firing instants of each EFT were then classified into motor units (MUs) according to their specific three-dimensional (3D) space position. The performance of the proposed approach was evaluated using simulated and experimentally recorded sEMG.

RESULTS

The simulation results demonstrated that the proposed approach can reconstruct MUAPTs with true positive rates of 89.12 ± 2.71%, 94.34 ± 1.85% and 95.45 ± 2.11% at signal-to-noise ratios of 10, 20, and 30 dB, respectively. The experimental results also demonstrated a high decomposition accuracy of 90.13 ± 1.31% in the two-source evaluation, and a high accuracy of 91.86 ± 1.14% in decompose-synthesize-decompose- compare evaluation.

CONCLUSIONS

The adoption of NMF reduces the dimension of random pattern under the restriction of non-negativity, as well as keeps the information unchanged as much as possible. The 3D space information of MUs enhances the classification accuracy by tackling the issue of relative movements between MUs and electrodes during dynamic contractions. The accuracy achieved in this study demonstrates the good performance and reliability of the proposed decomposition algorithm in dynamic surface EMG decomposition.

SIGNIFICANCE

The spatiotemporal information is applied to the dynamic surface EMG decomposition.

Single-Channel sEMG-Based Estimation of Knee Joint Angle Using a Decomposition Algorithm With a State-Space Model

Accurate human motion estimation is crucial for effective and safe human-robot interaction when using robotic devices for rehabilitation or performance enhancement. Although surface electromyography (sEMG) signals have been widely used to estimate human movements, conventional sEMG-based methods, which need sEMG signals measured from multiple relevant muscles, are usually subject to some limitations, including interference between sEMG sensors and wearable robots/environment, complicated calibration, as well as discomfort during long-term routine use. Few methods have been proposed to deal with these limitations by using single-channel sEMG (i.e., reducing the sEMG sensors as much as possible). The main challenge for developing single-channel sEMG-based estimation methods is that high estimation accuracy is difficult to be guaranteed. To address this problem, we proposed an sEMG-driven state-space model combined with an sEMG decomposition algorithm to improve the accuracy of knee joint movement estimation based on single-channel sEMG signals measured from gastrocnemius. The effectiveness of the method was evaluated via both single- and multi-speed walking experiments with seven and four healthy subjects, respectively. The results showed that the normal root-mean-squared error of the estimated knee joint angle using the method could be limited to 15%. Moreover, this method is robust with respect to variations in walking speeds. The estimation performance of this method was basically comparable to that of state-of-the-art studies using multi-channel sEMG.

High-Density Surface EMG Decomposition by Combining Iterative Convolution Kernel Compensation With an Energy-Specific Peel-off Strategy

Objective- This study aims to develop a novel framework for high-density surface electromyography (HD-sEMG) signal decomposition with superior decomposition yield and accuracy, especially for low-energy MUs. Methods- An iterative convolution kernel compensation-peel off (ICKC-P) framework is proposed, which consists of three steps: decomposition of the motor units (MUs) with relatively large energy by using the iterative convolution kernel compensation (ICKC) method and extraction of low-energy MUs with a Post-Processor and novel 'peel-off' strategy. Results- The performance of the proposed framework was evaluated by both simulated and experimental HD-sEMG signals. Our simulation results demonstrated that, with 120 simulated MUs, the proposed framework extracts more MUs compared to K-means convolutional kernel compensation (KmCKC) approach across six noise levels. And the proposed 'peel-off' strategy estimates more accurate MUAP waveforms at six noise levels than the 'peel-off' strategy proposed in the progressive FastICA peel-off (PFP) framework. For the experimental sEMG signals recorded from biceps brachii, an average of 16.1 ±3.4 MUs were identified from each contraction, while only 10.0 ± 2.8 MUs were acquired by the KmCKC method. Conclusion- The high yield and accuracy of MUs decomposed from simulated and experimental HD-sEMG signals demonstrate the superiority of the proposed framework in decomposing low-energy MUs compared to existing methods for HD-sEMG signal decomposition. Significance- The proposed framework enables us to construct a more representative motor unit pool, consequently enhancing our understanding pertaining to various neuropathological conditions and providing invaluable information for the diagnosis and treatment of neuromuscular disorders and motor neuron diseases.

Electromyography-Force Relation and Muscle Fiber Conduction Velocity Affected by Spinal Cord Injury

A surface electromyography (EMG) analysis was performed in this study to examine central neural and peripheral muscle changes after a spinal cord injury (SCI). A linear electrode array was used to record surface EMG signals from the biceps brachii (BB) in 15 SCI subjects and 14 matched healthy control subjects as they performed elbow flexor isometric contractions from 10% to 80% maximum voluntary contraction. Muscle fiber conduction velocity (MFCV) and BB EMG-force relation were examined. MFCV was found to be significantly slower in the SCI group than the control group, evident at all force levels. The BB EMG-force relation was well fit by quadratic functions in both groups. All healthy control EMG-force relations were best fit with positive quadratic coefficients. In contrast, the EMG-force relation in eight SCI subjects was best fit with negative quadratic coefficients, suggesting impaired EMG modulation at high forces. The alterations in MFCV and EMG-force relation after SCI suggest complex neuromuscular changes after SCI, including alterations in central neural drive and muscle properties.

Compressed spike-triggered averaging in iterative decomposition of surface EMG

BACKGROUND AND OBJECTIVE

Analysis of motor unit activity is important for assessing and treating diseases or injuries affecting natural movement. State-of-the-art decomposition translates high-density surface electromyography (HDsEMG) into motor unit activity. However, current decomposition methods offer far from complete separation of all motor units.

METHODS

This paper proposes a peel-off approach to automatic decomposition of HDsEMG into motor unit action potential (MUAP) trains, based on the Fast Independent Component Analysis algorithm (FastICA). The novel steps include utilizing compression by means of Principal Component Analysis and spike-triggered averaging, to estimate surface MUAP distributions with less noise, which are iteratively subtracted from the HDsEMG dataset. Furthermore, motor unit spike trains are estimated by high-dimensional density-based clustering of peaks in the FastICA source output. And finally, a new reliability measure is used to discard poor motor unit estimates by comparing the variance of the FastICA source output before and after the peel-off step. The method was validated using reconstructed synthetic data at three different signal-to-noise levels and was compared to an established deflationary FastICA approach.

RESULTS

Both algorithms had very high recall and precision, over 90%, for spikes from matching motor units, referred to as matched performance. However, the peel-off algorithm correctly identified more motor units for all noise levels. When accounting for unidentified motor units, total recall was up to 33 percentage points higher; and when accounting for duplicate estimates, total precision was up to 24 percentage points higher, compared to the state-of-the-art reference. In addition, a comparison was done using experimental data where the proposed algorithm had a matched recall of 97% and precision of 85% with respect to the reference algorithm.

CONCLUSION

These results show a substantial performance increase for decomposition of simulated HDsEMG data and serve to validate the proposed approach. This performance increase is an important step towards complete decomposition and extraction of information of motor unit activity.

Mapping Individual Motor Unit Activity to Continuous Three-DoF Wrist Torques: Perspectives for Myoelectric Control

The surface electromyography (EMG) decomposition techniques provide access to motor neuron activities and have been applied to myoelectric control schemes. However, the current decomposition-based myoelectric control mainly focuses on discrete gestures or single-DoF continuous movements. In this study, we aimed to map the motor unit discharges, which were identified from high-density surface EMG, to the three degrees of freedom (DoFs) wrist movements. The 3-DoF wrist torques and high-density surface EMG signals were recorded concurrently from eight non-disabled subjects. The experimental protocol included single-DoF movements and their various combinations. We decoded the motor unit discharges from the EMG signals using a segment-wise decomposition algorithm. Then the neural features were extracted from motor unit discharges and projected to wrist torques with a multiple linear regression model. We compared the performance of two neural features (twitch model and spike counting) and two training schemes (single-DoF and multi-DoF training). On average, 145 ± 33 motor units were identified from each subject, with a pulse-to-noise ratio of 30.8 ± 4.2 dB. Both neural features exhibited high estimation accuracy of 3-DoF wrist torques, with an average [Formula: see text] of 0.76 ± 0.12 and normalized root mean square error of 11.4 ± 3.1%. These results demonstrated the efficiency of the proposed method in continuous estimation of 3-DoF wrist torques, which has the potential to advance dexterous myoelectric control based on neural information.

Neural inputs from spinal motor neurons to lateralis vastus muscle: Comparison between sprinters and nonathletes

The adaptation of neural contractile properties has been observed in previous work. However, the neural changes on the motor unit (MU) level remain largely unknown. Voluntary movements are controlled through the precise activation of MU populations. In this work, we estimate the neural inputs from the spinal motor neurons to the muscles during isometric contractions and characterize the neural adaptation during training by comparing the MU properties decomposed from sprinters and nonathletes. Twenty subjects were recruited and divided into two groups. The high-density surface electromyography (EMG) signals were recorded from the lateralis vastus muscle during the isometric contraction of knee extension and were then decomposed into MU spike trains. Each MU’s action potentials and discharge properties were extracted for comparison across subject groups and tasks. A total of 1097 MUs were identified from all subjects. Results showed that the discharge rates and amplitudes of MUAPs from athletes were significantly higher than those from nonathletes. These results demonstrate the neural adaptations in physical training at the MU population level and indicate the great potential of EMG decomposition in physiological investigations.

Single-channel EEG signal extraction based on DWT, CEEMDAN, and ICA method

In special application scenarios, such as portable anesthesia depth monitoring, portable emotional state recognition and portable sleep monitoring, electroencephalogram (EEG) signal acquisition equipment is required to be convenient and easy to use. It is difficult to remove electrooculogram (EOG) artifacts when the number of EEG acquisition channels is small, especially when the number of observed signals is less than that of the source signals, and the overcomplete problem will arise. The independent component analysis (ICA) algorithm commonly used for artifact removal requires the number of basis vectors to be smaller than the dimension of the input data due to a set of standard orthonormal bases learned during the convergence process, so it cannot be used to solve the overcomplete problem. The empirical mode decomposition method decomposes the signal into several independent intrinsic mode functions so that the number of observed signals is more than that of the source signals, solving the overcomplete problem. However, when using this method to solve overcompleteness, the modal aliasing problem will arise, which is caused by abnormal events such as sharp signals, impulse interference, and noise. Aiming at the above problems, we propose a novel EEG artifact removal method based on discrete wavelet transform, complete empirical mode decomposition for adaptive noise (CEEMDAN) and ICA in this paper. First, the input signals are transformed by discrete wavelet (DWT), and then CEEMDAN is used to solve the overcomplete and mode aliasing problems, meeting the a priori conditions of the ICA algorithm. Finally, the components belonging to EOG artifacts are removed according to the sample entropy value of each independent component. Experiments show that this method can effectively remove EOG artifacts while solving the overcomplete and modal aliasing problems.

Segment-Wise Decomposition of Surface Electromyography to Identify Discharges Across Motor Neuron Populations

OBJECTIVE

The surface electromyography (EMG) decomposition techniques have shown promising results in neurophysiologic investigations, clinical diagnosis, and human-machine interfacing. However, current decomposition methods could only decode a limited number of motor units (MUs) because of the local convergence. The number of identified MUs remains similar even though more muscles or movements are involved, where multiple motor neuron populations are activated. The objective of this study was to develop a segment-wise decomposition strategy to increase the number of MU decoded from multiple motor neuron populations.

METHODS

The EMG signals were divided into several segments depending on the number of involved movements. The motor neurons, activated during each movement, were regarded as a population. The convolution kernel compensation (CKC) method was applied individually for each segment to decode the motor unit discharges from each motor neuron population. The MU filters were obtained in each segment and filtrated to estimate the MU spike trains (MUSTs) from the global EMG signals. The decomposition performance was validated on synthetic and experimental EMG signals.

MAIN RESULTS

From synthetic EMG signals generated by two motor neuron populations, the proposed segment-wise CKC (swCKC) decoded significantly more MUs during low and medium excitation levels, with an increased rate of 16.3% to 75.4% compared with the conventional CKC. From experimental signals recorded during ten motor tasks, 133±24 MUs with the pulse-to-noise ratio of 36.6±6.5 dB were identified for each subject by swCKC, whereas the conventional CKC identified only 43±12 MUs.

CONCLUSION AND SIGNIFICANCE

These results indicate the feasibility and superiority of the proposed swCKC to decode MU activities across motor neuron populations, extending the potential applications of EMG decomposition for neural decoding during multiple motor tasks.

Finger Force Estimation using Motor Unit Discharges Across Forearm Postures

BACKGROUND

Myoelectric-based decoding has gained popularity in upper-limb neural-machine interfaces. Motor unit (MU) firings decomposed from surface electromyographic (EMG) signals can represent motor intent, but EMG properties at different arm configurations can change due to electrode shift and differing neuromuscular states. This study investigated whether isometric fingertip force estimation using MU firings is robust to forearm rotations from a neutral to either a fully pronated or supinated posture.

METHODS

We extracted MU information from high-density EMG of the extensor digitorum communis in two ways: (1) Decomposed EMG in all three postures (MU-AllPost); and (2) Decomposed EMG in neutral posture (MU-Neu), and extracted MUs (separation matrix) were applied to other postures. Populational MU firing frequency estimated forces scaled to subjects' maximum voluntary contraction (MVC) using a regression analysis. The results were compared with the conventional EMG-amplitude method.

RESULTS

We found largely similar root-mean-square errors (RMSE) for the two MU-methods, indicating that MU decomposition was robust to postural differences. MU-methods demonstrated lower RMSE in the ring (EMG = 6.23, MU-AllPost = 5.72, MU-Neu = 5.64 %MVC) and pinky (EMG = 6.12, MU-AllPost = 4.95, MU-Neu = 5.36 %MVC) fingers, with mixed results in the middle finger (EMG = 5.47, MU-AllPost = 5.52, MU-Neu = 6.19% MVC).

CONCLUSION

Our results suggest that MU firings can be extracted reliably with little influence from forearm posture, highlighting its potential as an alternative decoding scheme for robust and continuous control of assistive devices.

A novel concatenate feature fusion RCNN architecture for sEMG-based hand gesture recognition

Hand gesture recognition tasks based on surface electromyography (sEMG) are vital in human-computer interaction, speech detection, robot control, and rehabilitation applications. However, existing models, whether traditional machine learnings (ML) or other state-of-the-arts, are limited in the number of movements. Targeting a large number of gesture classes, more data features such as temporal information should be persisted as much as possible. In the field of sEMG-based recognitions, the recurrent convolutional neural network (RCNN) is an advanced method due to the sequential characteristic of sEMG signals. However, the invariance of the pooling layer damages important temporal information. In the all convolutional neural network (ACNN), because of the feature-mixing convolution operation, a same output can be received from completely different inputs. This paper proposes a concatenate feature fusion (CFF) strategy and a novel concatenate feature fusion recurrent convolutional neural network (CFF-RCNN). In CFF-RCNN, a max-pooling layer and a 2-stride convolutional layer are concatenated together to replace the conventional simple dimensionality reduction layer. The featurewise pooling operation serves as a signal amplitude detector without using any parameter. The feature-mixing convolution operation calculates the contextual information. Complete evaluations are made on both the accuracy and convergence speed of the CFF-RCNN. Experiments are conducted using three sEMG benchmark databases named DB1, DB2 and DB4 from the NinaPro database. With more than 50 gestures, the classification accuracies of the CFF-RCNN are 88.87% on DB1, 99.51% on DB2, and 99.29% on DB4. These accuracies are the highest compared with reported accuracies of machine learnings and other state-of-the-art methods. To achieve accuracies of 86%, 99% and 98% for the RCNN, the training time are 2353.686 s, 816.173 s and 731.771 s, respectively. However, for the CFF-RCNN to reach the same accuracies, it needs only 1727.415 s, 542.245 s and 576.734 s, corresponding to a reduction of 26.61%, 33.56% and 21.19% in training time. We concluded that the CFF-RCNN is an improved method when classifying a large number of hand gestures. The CFF strategy significantly improved model performance with higher accuracy and faster convergence as compared to traditional RCNN.

Decoding Muscle Force From Motor Unit Firings Using Encoder-Decoder Networks

Appropriate interpretation of motor unit (MU) activities after surface EMG (sEMG) decomposition is a key factor to decode motor intentions in a noninvasive and physiologically meaningful way. However, there are great challenges due to the difficulty in cross-trial MU tracking and unavoidable loss of partial MU information resulting from incomplete decomposition. In light of these issues, this study presents a novel framework for interpreting MU activities and applies it to decode muscle force. The resulting MUs were clustered and classified into different categories by characterizing their spatially distributed firing waveforms. The process served as a general MU tracking method. On this basis, after transferring the MU firing trains to twitch force trains by a twitch force model, a deep network was designed to predict the normalized force. In addition, MU category distribution was examined to calibrate the actual force level, while functions of some unavailable MUs were compensated. To investigate the effectiveness of this framework, high-density sEMG signals were recorded using an 8×8 electrode array from the abductor pollicis brevis muscles of eight subjects, while thumb abduction force was measured. The proposed method outperformed three common methods ( ) yielding the lowest root mean square deviation of 6.68% ± 1.29% and the highest fitness ( R² ) of 0.94 ± 0.04 between the predicted force and the actual force. This study offers a valuable, computational solution for interpreting individual MU activities, and its effectiveness was confirmed in muscle force estimation.

First-recruited motor units adopt a faster phenotype in amyotrophic lateral sclerosis

KEY POINTS

Amyotrophic lateral sclerosis (ALS) is an incurable neurodegenerative disorder of motor neurons, carrying a short survival. High-density motor unit recordings permit analysis of motor unit size (amplitude) and firing behaviour (afterhyperpolarization duration and muscle fibre conduction velocity). Serial recordings from biceps brachii indicated that motor units fired faster and with greater amplitude as disease progressed. First-recruited motor units in the latter stages of ALS developed characteristics akin to fast-twitch motor units, possibly as a compensatory mechanism for the selective loss of this motor unit subset. This process may become maladaptive, highlighting a novel therapeutic target to reduce motor unit vulnerability.

ABSTRACT

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disorder with a median survival of 3 years. We employed serial high-density surface electromyography (HDSEMG) to characterize voluntary and ectopic patterns of motor unit (MU) firing at different stages of disease. By distinguishing MU subtypes with variable vulnerability to disease, we aimed to evaluate compensatory neuronal adaptations that accompany disease progression. Twenty patients with ALS and five patients with benign fasciculation syndrome (BFS) underwent 1-7 assessments each. HDSEMG measurements comprised 30 min of resting muscle and 1 min of light voluntary activity from biceps brachii bilaterally. MU decomposition was performed by the progressive FastICA peel-off technique. Inter-spike interval, firing pattern, MU potential area, afterhyperpolarization duration and muscle fibre conduction velocity were determined. In total, 373 MUs (ALS = 287; BFS = 86) were identified from 182 recordings. Weak ALS muscles demonstrated a lower mean inter-spike interval (82.7 ms) than strong ALS muscles (96.0 ms; P = 0.00919) and BFS muscles (95.3 ms; P = 0.0039). Mean MU potential area (area under the curve: 487.5 vs. 98.7 μV ms; P < 0.0001) and muscle fibre conduction velocity (6.2 vs. 5.1 m/s; P = 0.0292) were greater in weak ALS muscles than in BFS muscles. Purely fasciculating MUs had a greater mean MU potential area than MUs also under voluntary command (area under the curve: 679.6 vs. 232.4 μV ms; P = 0.00144). These results suggest that first-recruited MUs develop a faster phenotype in the latter stages of ALS, likely driven by the preferential loss of vulnerable fast-twitch MUs. Inhibition of this potentially maladaptive phenotypic drift may protect the longevity of the MU pool, stimulating a novel therapeutic avenue.

Properties of the surface electromyogram following traumatic spinal cord injury: a scoping review

Traumatic spinal cord injury (SCI) disrupts spinal and supraspinal pathways, and this process is reflected in changes in surface electromyography (sEMG). sEMG is an informative complement to current clinical testing and can capture the residual motor command in great detail-including in muscles below the level of injury with seemingly absent motor activities. In this comprehensive review, we sought to describe how the sEMG properties are changed after SCI. We conducted a systematic literature search followed by a narrative review focusing on sEMG analysis techniques and signal properties post-SCI. We found that early reports were mostly focused on the qualitative analysis of sEMG patterns and evolved to semi-quantitative scores and a more detailed amplitude-based quantification. Nonetheless, recent studies are still constrained to an amplitude-based analysis of the sEMG, and there are opportunities to more broadly characterize the time- and frequency-domain properties of the signal as well as to take fuller advantage of high-density EMG techniques. We recommend the incorporation of a broader range of signal properties into the neurophysiological assessment post-SCI and the development of a greater understanding of the relation between these sEMG properties and underlying physiology. Enhanced sEMG analysis could contribute to a more complete description of the effects of SCI on upper and lower motor neuron function and their interactions, and also assist in understanding the mechanisms of change following neuromodulation or exercise therapy.

Open Access Dataset, Toolbox and Benchmark Processing Results of High-Density Surface Electromyogram Recordings

We provide an open access dataset of High densitY Surface Electromyogram (HD-sEMG) Recordings (named "Hyser"), a toolbox for neural interface research, and benchmark results for pattern recognition and EMG-force applications. Data from 20 subjects were acquired twice per subject on different days following the same experimental paradigm. We acquired 256-channel HD-sEMG from forearm muscles during dexterous finger manipulations. This Hyser dataset contains five sub-datasets as: (1) pattern recognition (PR) dataset acquired during 34 commonly used hand gestures, (2) maximal voluntary muscle contraction (MVC) dataset while subjects contracted each individual finger, (3) one-degree of freedom (DoF) dataset acquired during force-varying contraction of each individual finger, (4) N-DoF dataset acquired during prescribed contractions of combinations of multiple fingers, and (5) random task dataset acquired during random contraction of combinations of fingers without any prescribed force trajectory. Dataset 1 can be used for gesture recognition studies. Datasets 2-5 also recorded individual finger forces, thus can be used for studies on proportional control of neuroprostheses. Our toolbox can be used to: (1) analyze each of the five datasets using standard benchmark methods and (2) decompose HD-sEMG signals into motor unit action potentials via independent component analysis. We expect our dataset, toolbox and benchmark analyses can provide a unique platform to promote a wide range of neural interface research and collaboration among neural rehabilitation engineers.

Muscle Force Estimation Based on Neural Drive Information From Individual Motor Units

Estimation of muscle contraction force based on the macroscopic feature of surface electromyography (SEMG) has been widely reported, but the use of microscopic neural drive information has not been thoroughly investigated. In this study, a novel method is proposed to process individual motor unit (MU) activities (firing sequences and action potential waveforms) derived from the decomposition of high density SEMG (HD-SEMG), and it is applied to muscle force estimation. In the proposed method, a supervised machine learning approach was conducted to determine the twitch force of each MU according to its action potential waveforms, which enables separate calculation of every MU's contribution to force. Thus, the muscle force was predicted through a physiologically meaningful muscle force model. In the experiment, HD-SEMG data were recorded from the abductor pollicis brevis muscles of eight healthy subjects during their performance of thumb abduction with the force increasing gradually from zero to four force levels (10%, 20%, 30%, 40% of the maximal voluntary contraction), while the true muscle force was measured simultaneously. When the proposed method was used, the root mean square difference (RMSD) of the error of the estimated force with respect to the measured force was reported to be 8.3% ± 2.8%. The proposed method also significantly outperformed the other four common methods for force estimation (RMSD: from 11.7% to 20%, ), demonstrating its effectiveness. This study offers a useful tool for exploiting the neural drive information towards muscle force estimation with improved precision. The proposed method has wide applications in precise motor control, sport and rehabilitation medicine.

Torque Estimation Using Neural Drive for a Concentric Contraction

The scope and relevance of wearable robotics spans across a number of research fields with a variety of applications. A challenge across these research areas is improving user-interface control. One established approach is using neural control interfaces derived from surface electromyography (sEMG). Although there has been some success with sEMG controlled prosthetics, the coarse nature of traditional sEMG processing has limited the development of fully functional prosthetics and wearable robotics. To solve this problem, blind source separation (BSS) techniques have been implemented to extract the user's movement intent from high-density sEMG (HDsEMG) measurements; however, current methods have only been well validated during static, low-level muscle contractions, and it is unclear how they will perform during movement. In this paper we present a neural drive based method for predicting output torque during a constant force, concentric contraction. This was achieved by modifying an existing HDsEMG decomposition algorithm to decompose 1 sec. overlapping windows. The neural drive profile was computed using both rate coding and kernel smoothing. Neither rate coding nor kernel smoothing performed as well as HDsEMG amplitude estimation, indicating that there are still significant limitations in adapting current methods to decompose dynamic contractions, and that sEMG amplitude estimation methods still remain highly reliable estimators.

Hybrid Encoder-Decoder Deep Networks for Decoding Neural Drive Information towards Precise Muscle Force Estimation

How to utilize and interpret microscopic motor unit (MU) activities after surface electromyogram (sEMG) decomposition towards accurate decoding of the neural control remains a great challenge. In this study, a novel framework of hybrid encoder-decoder deep networks is proposed to process the microscopic neural drive information and it is applied to precise muscle force estimation. After a high-density sEMG (HD-sEMG) decomposition was performed using the progressive FastICA peel-off algorithm, a muscle twitch force model was then applied to basically convert each channel's electric waveform (i.e., action potential) of each MU into a twitch force. Next, hybrid encoder-decoder deep networks were performed on every 50 ms of segment of the summation of twitch force trains from all decomposed MUs. The encoder network was designed to characterize spatial information of MU's force contribution over all channels, and the decoder network finally decoded the muscle force. This framework was validated on HD-sEMG recordings from the abductor pollicis brevis muscles of five subjects by a thumb abduction task using an 8 × 8 grid. The proposed framework yielded a mean root mean square error of 6.62% ± 1.26% and a mean coefficient of determination value of 0.95 ± 0.03 from a linear regression analysis between the estimated force and actual force over all data trials, and it outperformed three common methods with statistical significance (p < 0.001). This study offers a valuable solution for interpreting microscopic neural drive information and demonstrates its success in predicting muscle force.

A Novel Validation Approach for High-Density Surface EMG Decomposition in Motor Neuron Disease

This paper presents a novel two-source approach for validating the performance of high-density surface electromyogram (EMG) decomposition. The approach was developed taking advantage of surface EMG characteristics of amyotrophic lateral sclerosis (ALS). High-density surface EMG data from ALS patients can be divided to the sparse data set and the interference data set, with the former decomposed by expert visual inspection while the latter independently decomposed by the surface EMG decomposition algorithm. The agreement of the decomposition yields from the two data sets can be quantified for evaluating the surface EMG decomposition performance. The novel validation approach was performed for a recently developed method called automatic progressive FastICA peel-off (APFP) for high-density surface EMG decomposition. The APFP framework was used to automatically decompose high-density surface EMG signals recorded from the first dorsal interosseous muscle of ALS subjects. The common motor units independently decomposed from the interference data set and the sparse data set demonstrated an average matching rate of 99.18% ± 1.11%. The characteristics of the ALS surface EMG also facilitate a step by step illustration of the APFP framework for high-density surface EMG decomposition. The novel approach presented in this paper can supplement conventional two-source validation for accuracy assessment of decomposed motor units from experimental signals, which is essential for development of surface EMG decomposition methods.

Non-invasive Decoding of the Motoneurons: A Guided Source Separation Method Based on Convolution Kernel Compensation With Clustered Initial Points

Despite the progress in understanding of neural codes, the studies of the cortico-muscular coupling still largely rely on interferential electromyographic (EMG) signal or its rectification for the assessment of motor neuron pool behavior. This assessment is non-trivial and should be used with precaution. Direct analysis of neural codes by decomposing the EMG, also known as neural decoding, is an alternative to EMG amplitude estimation. In this study, we propose a fully-deterministic hybrid surface EMG (sEMG) decomposition approach that combines the advantages of both template-based and Blind Source Separation (BSS) decomposition approaches, a.k.a. guided source separation (GSS), to identify motor unit (MU) firing patterns. We use the single-pass density-based clustering algorithm to identify possible cluster representatives in different sEMG channels. These cluster representatives are then used as initial points of modified gradient Convolution Kernel Compensation (gCKC) algorithm. Afterwards, we use the Kalman filter to reduce the noise impact and increase convergence rate of MU filter identification by gCKC. Moreover, we designed an adaptive soft-thresholding method to identify MU firing times out of estimated MU spike trains. We tested the proposed algorithm on a set of synthetic sEMG signals with known MU firing patterns. A grid of 9 × 10 monopolar surface electrodes with 5-mm inter-electrode distances in both directions was simulated. Muscle excitation was set to 10, 30, and 50%. Colored Gaussian zero-mean noise with the signal-to-noise ratio (SNR) of 10, 20, and 30 dB, respectively, was added to 16 s long sEMG signals that were sampled at 4,096 Hz. Overall, 45 simulated signals were analyzed. Our decomposition approach was compared with gCKC algorithm. Overall, in our algorithm, the average numbers of identified MUs and Rate-of-Agreement (RoA) were 16.41 ± 4.18 MUs and 84.00 ± 0.06%, respectively, whereas the gCKC identified 12.10 ± 2.32 MUs with the average RoA of 90.78 ± 0.08%. Therefore, the proposed GSS method identified more MUs than the gCKC, with comparable performance. Its performance was dependent on the signal quality but not the signal complexity at different force levels. The proposed algorithm is a promising new offline tool in clinical neurophysiology.

Automatic Multichannel Intramuscular Electromyogram Decomposition: Progressive FastICA Peel-Off and Performance Validation

The progressive FastICA peel-off (PFP) is a recently developed blind source separation approach for high-density surface EMG decomposition. This paper explores a novel application of PFP for automatic decomposition of multi-channel intramuscular electromyogram signals. The automatic PFP (APFP) was used to decompose an open access multichannel intramuscular EMG dataset, simultaneously collected from the brachioradialis muscle using 6 to 8 fine wire or needle electrodes. Given usually a limited number of intramuscular electrodes compared with high-density surface EMG recording, a modification was made to the original APFP framework to dramatically increase the decomposition yield. A total of 131 motor units were automatically decomposed by the APFP framework from 10 multichannel intramuscular EMG signals, among which 128 motor units were also manually identified from the expert interactive EMGLAB decomposition. The average matching rate of discharge instants for all the common motor units was (98.71 ± 1.73)%. The outcomes of this study indicate that the APFP framework can also be used to automatically decompose multichannel intramuscular EMG with high accuracies, even though the number of recording channels is relatively small compared with high-density surface EMG.

Two-Source Validation of Progressive FastICA Peel-Off for Automatic Surface EMG Decomposition in Human First Dorsal Interosseous Muscle

This study aims to assess the accuracy of a novel high density surface electromyogram (SEMG) decomposition method, namely automatic progressive FastICA peel-off (APFP), for automatic decomposition of experimental electrode array SEMG signals. A two-source method was performed by simultaneous concentric needle EMG and electrode array SEMG recordings from the human first dorsal interosseous (FDI) muscle, using a protocol commonly applied in clinical EMG examination. The electrode array SEMG was automatically decomposed by the APFP while the motor unit action potential (MUAP) trains were also independently identified from the concentric needle EMG. The degree of agreement of the common motor unit (MU) discharge timings decomposed from the two different categories of EMG signals was assessed. A total of 861 and 217 MUs were identified from the 114 trials of simultaneous high density SEMG and concentric needle EMG recordings, respectively. Among them 168 common (MUs) were found with a high average matching rate of [Formula: see text] for the discharge timings. The outcomes of this study show that the APFP can reliably decompose at least a subset of MUs in the high density SEMG signals recorded from the human FDI muscle during low contraction levels using a protocol analog to clinical EMG examination.

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.