Use of surface EMG for studying motor unit recruitment during isometric linear force ramp.

The aim of this work was to demonstrate the rank order of motor unit (MU) recruitment by surface EMG based on a Laplacian detection technique and to document the MU features at their recruitment threshold. Surface EMG signals were recorded on the biceps brachii of 10 healthy subjects during linear force ramps. When achievable, the signals were decomposed into MU action potential (MUAP) trains. MU inter-pulse interval (IPI), conduction velocity (MUCV) and amplitude were estimated on the first 12 MUAPs of each detectable train in order to characterize the MU features at their firing onset. A strong correlation was found between MU recruitment threshold and IPI, MUCV, and amplitude, showing that the size principle can be demonstrated by a fully non-invasive EMG technique. However, signal decomposition was not possible on seven subjects due to the effects of the volume conductor when the skinfold thickness was too large. When requirements for an optimal detection of MUAP trains are met, surface EMG may be used to improve our understanding of MU activity.     

Automatic decomposition of multichannel intramuscular EMG signals

We describe an automatic algorithm for decomposing multichannel EMG signals into their component motor unit action potential (MUAP) trains, including signals from widely separated recording sites in which MUAPs exhibit appreciable interchannel offset and jitter. The algorithm has two phases. In the clustering phase, the distinct, recurring MUAPs in each channel are identified, the ones that correspond to the same motor units are determined by their temporal relationships, and multichannel templates are computed. In the identification stage, the MUAP discharges in the signal are identified using matched filtering and superimposition resolution techniques. The algorithm looks for the MUAPs with the largest single channel components first, using matches in one channel to guide the search in other channels, and using information from the other channels to confirm or refute each identification. For validation, the algorithm was used to decompose 10 real 6-to-8-channel EMG signals containing activity from up to 25 motor units. Comparison with expert manual decomposition showed that the algorithm identified more than 75% of the total 176 MUAP trains with an accuracy greater than 95%. The algorithm is fast, robust, and shows promise to be accurate enough to be a useful tool for decomposing multichannel signals. It is freely available at http://emglab.stanford.edu.     

Computer simulation study of the shape of motor unit action potential

Motor unit action potentials (MUAPs) of brachial biceps were simulated. A simulated MUAP was obtained as a sum of single fibre action potentials (SFAPs) from all the muscle fibres of a motor unit (MU). The influence of the following factors on MUAP shape for different kinds of recording electrode was studied: fibre density, neuromuscular jitter, temporal dispersion and electrode displacements. The simulation confirms that typical MUAPs recorded with needle electrodes from muscles of low fibre density such as brachial biceps are usually triphasic. Increased fibre density produces MUAPs of more complex shape and higher amplitude. Normal neuromuscular jitter is responsible for the variability of shape of subsequent potentials from the same MU as well as for electromyographic shimmer. Pathologic (increased) jitter makes the shapes of subsequent potentials unrecognizable. The influence of temporal dispersion is interconnected with other factors but rather of minor importance. The simulation shows how big changes in MUAP shape can be expected due to electrode displacements during single experiment or during estimation of MU territory.     

Using two-dimensional spatial information in decomposition of surface EMG signals.

Recently, high-density surface EMG electrode grids and multi-channel amplifiers became available for non-invasive recording of human motor units (MUs). We present a way to decompose surface EMG signals into MU firing patterns, whereby we concentrate on the importance of two-dimensional spatial differences between the MU action potentials (MUAPs). Our method is exemplified with high-density EMG data from the vastus lateralis muscle of a single subject. Bipolar and Laplacian spatial filtering was applied to the monopolar raw signals. From the single recording in this subject six different simultaneously active MUs could be distinguished using the spatial differences between MUAPs in the direction perpendicular to the muscle fiber direction. After spike-triggered averaging, 125-channel two-dimensional MUAP templates were obtained. Template-matching allowed tracking of all MU firings. The impact of spatial information was measured by using subsets of the MUAP templates, either in parallel or perpendicular to the muscle fiber direction. The use of one-dimensional spatial information perpendicular to the muscle fiber direction was superior to the use of a linear array electrode in the longitudinal direction. However, to detect the firing events of the MUs with a high accuracy, as needed for instance for estimation of firing synchrony, two-dimensional information from the complete grid electrode appears essential.     

Motor unit tracking with high-density surface EMG

Following (tracking) individual motor units over time can provide important new insights, both into the relationships among various motor unit (MU) morphological and functional properties and into how these properties are influenced by neuromuscular disorders or interventions. The present study aimed to determine whether high-density surface EMG (HD-sEMG) recordings, which use an array of surface electrodes over a muscle, can increase the yield of MU tracking studies in terms of the number of MUs that can be tracked. For that purpose, four HD-sEMG recording sessions were performed on the thenar muscles of ten healthy subjects. Decomposition of the recorded composite responses yielded a study total of 2849 motor unit action potentials (MUAPs). MUAPs that were found in both of the first two sessions, performed on the same day, were defined as trackable MUAPs. Our results show that 22 (median value; range, 13–34) MUAPs per nerve were trackable, which represented approximately 5% of the total MU population. Of these trackable MUAPs, 16 (11–26) could also be found in one or both of the third and fourth sessions, which were performed between 1 and 13 weeks after the initial studies. Nine (4–18) MUAPs were found in all four sessions. Many of the characteristic MUAP shapes matched well between sessions, even when these sessions were several weeks apart. However, some MUAPs seem very sensitive to changes in arm position or in the muscle’s morphology (e.g., to changes in muscle fiber length due to variable degrees of thumb flexion or extension), particularly those from larger and/or superficial MUs. Standardization is, therefore, essential to detect even small MUAP changes, as may occur with pathology or interventions. If this is accomplished, MU tracking with HD-sEMG may prove to be a powerful tool for a promising type of neurophysiological investigation.     

Short term bed-rest reduces conduction velocity of individual motor units in leg muscles

Space permanence simulations such as prolonged bed-rest can mimic some of the physiological modifications in the human body and provide study conditions that are more accessible than during space flight. A short term bed-rest experiment was organized to simulate the effects of weightlessness for studying the adaptation to this condition. Eight healthy young volunteers were studied before and immediately after the 14 day periods of strict bed-rest.Surface EMG signals were detected with linear electrode arrays from vastus medialis, vastus lateralis and tibialis anterior muscle during isometric voluntary contractions at 20% MVC. Motor unit action potentials (MUAPs) of individual motor units were extracted from the interference EMG signals with a partial decomposition algorithm and averaged.MUAP templates generated by the same motor unit could be retrieved before and after bed-rest period. Muscle fiber conduction velocity (CV) was estimated from each averaged MUAP template and from the global EMG signal. Both global and single MU conduction velocity was observed to decrease by about 10% after the bed-rest period (p < 0.05). Amplitude and power spectral parameters did not significantly change after the bed-rest period.It is concluded that a short term bed-rest reduces the CV of individual motor units without a significant effect on muscle force or on other electrophysiological parameters.     

Single motor unit and fiber action potentials during fatigue

Muscle fatigue is defined as a loss of tension development during constant stimulation. Although the relationship is not well documented, muscle fatigue has been inferred from electromyogram (EMG) signals. The purpose of this study was to determine the relationship between the amplitude and duration of single motor unit action potentials (MUAPs) and the loss of tension development (fatigue) in the medial gastrocnemius muscles of cats. Single motor units were fatigued by continuous stimulation at 10 or 80 Hz or with trains of 40-Hz stimuli. When motor units were stimulated at 10 Hz and with trains at 40 Hz (low frequency), tension declined and remained depressed during recovery. The changes in the MUAP correlated poorly with changes in tension. During and after stimulation at 80 Hz (high frequency), changes in the amplitude and duration of MUAPs correlated highly with changes in tension development. Since the EMG signal is dependent on a summation and cancellation of individual MUAPs, the EMG provides a reasonable estimate of high-frequency fatigue but an unreliable measure of low-frequency fatigue.     

Automatic localisation of innervation zones: A simulation study of the external anal sphincter

Traumas of the innervation zone (IZ) of the external anal sphincter (EAS), e.g. during delivery, can promote the development of faecal incontinence. Recently developed probes allow high-resolution detection of EMG signals from the EAS. The analysis of pelvic floor muscles by surface EMG (in particular, the estimation of the location of the IZ) has potential applications in the diagnosis and investigation of the mechanisms of incontinence.An automatic method (based on matched filter approach) for the estimation of the IZ distribution of EAS from surface EMG is discussed and tested using an analytical model of generation of EMG signals from sphincter muscles. Simulations are performed varying length of the fibres, thickness of the mucosa, position of the motor units, and force level. Different distributions of IZs are simulated.The performance of the proposed method in the estimation of the IZ distribution is affected by surface MUAP amplitude (as the estimation made by visual inspection), by mucosa thickness (performance decreases when fibre length is higher) and by different MU distributions. However, in general the method is able to identify the position of two IZ locations and can measure asymmetry of the IZ distribution. This strengthens the potential applications of high density surface EMG in the prevention and investigation of incontinence.     

Modeling the vestibular evoked myogenic potential

Measuring the vestibular evoked myogenic potential (VEMP) promises to become a routine method for assessing vestibular function, although the technique is not yet standardized. To overcome the problem that the VEMP amplitude depends not only on the inhibition triggered by the acoustic stimulation of the vestibular end organs in the inner ear, but also on the tone of the muscle from which the potential is recorded, the VEMP is often normalized by dividing through a measure of the electromyogram (EMG) activity. The underlying idea is that VEMP amplitude and EMG activity are proportional. But this would imply that the muscle tone is irrelevant for a successful VEMP recording, contradicting experimental evidence. Here, an analytical model is presented that allows to resolve the contradiction. The EMG is modeled as the sum of motor unit action potentials (MUAPs). A brief inhibition can be characterized by its equivalent rectangular duration (ERD), irrespective of the actual time course of the inhibition. The VEMP resembles a polarity-inverted MUAP under such circumstances. Its amplitude is proportional to both the ERD and the MUAP rate. The EMG activity, by contrast, is proportional to the square root of the MUAP rate. Thus, the normalized VEMP still depends on the muscle tone. To avoid confounding effects of the muscle tone, the standard deviation of the EMG could be considered. But the inhibition effect on the standard deviation is small so that the measuring time would have to be much longer than usual today.     

Detection of individual motor units of the puborectalis muscle by non-invasive EMG electrode arrays

PURPOSE: The purpose of the study was to demonstrate that anatomical features of individual motor units of the puborectalis muscle can be detected with non-invasive electromyography (EMG) and to evaluate differences in electrophysiological properties of the puborectalis muscles in a small group of healthy and pathologic subjects. METHODS: Multichannel EMG was recorded by means of a flexible probe applied on the gloved index finger and carrying an array of eight equally spaced (1.15 mm) electrodes. A multichannel EMG amplifier provided seven outputs corresponding to the pairs of adjacent electrodes. Tests were performed in three different positions (dorsal, left and right) over the puborectalis muscle on 20 subjects (nine healthy, seven constipated and four incontinent patients). Motor unit action potentials (MUAPs) generated at the innervation zone of a MU and propagating along the muscle fibers generated repetitive characteristic patterns on the seven output channels allowing identification of anatomical features of the motor units. RESULTS: MUAPs were observed travelling in either one or both directions with the array in dorsal position, and mainly in ventral-to-dorsal direction in either lateral position. MUAP amplitude was lower in constipated and incontinent patients with respect to healthy subjects. The conduction velocity estimated on the identified MUAPs was lower for constipated patients with respect to healthy subjects suggesting different mechanical properties of the active motor units. CONCLUSIONS: This technique allows the extraction of relevant information about the anatomical features (innervation zone position and overlapping of motor unit branches) of the puborectalis muscle and its electrophysiological properties and maybe can be applied as an novel methodology for assessing the anorectal function in patients.     

Determination of feature relevance for the grouping of motor unit action potentials through a generative mixture model

The study of motor unit action potential (MUAP) activity from electromyographic signals is important for neurological investigations aiming to understand the state of the neuromuscular system. In this context, the identification and clustering of MUAPs that exhibit common characteristics, and the assessment of which data features are most relevant for the definition of such cluster structure, are central issues. In this paper, we propose the application of an unsupervised feature relevance determination (FRD) method to the analysis of experimental MUAPs. This method is embedded in a constrained mixture of distributions model that simultaneously performs data clustering and visualization. The experimental results of the analysis of a data set consisting of MUAPs measured from the First Dorsal Interosseous, a hand muscle, indicate that the features corresponding to the hyperpolarization period in the physiological process of generating muscle fibre action potentials are consistently estimated to be the most relevant. Moreover, the MUAP cluster structure of the data is shown to be only partially attributable to inter-subject differences, with the hyperpolarization period providing the best discrimination of the data by subject.     

Fuzzy MUAP recognition in HSR-EMG detection basing on morphological features

The idea of ‘besides the MU properties and depending on the recording techniques, MUAPs can have unique pattern’ was adopted. The aim of this work was to recognise whether a Laplacian-detected MUAP is isolated or overlapped basing on novel morphological features using fuzzy classifier. Training data set was constructed to elaborate and test the ‘if-then’ fuzzy rules using signals provided by three muscles: the abductor pollicis brevis (APB), the first dorsal interosseous (FDI) and the biceps brachii (BB) muscles of 11 healthy subjects. The proposed fuzzy classier recognized automatically the isolated MUAPs with a performance of 95.03% which was improved to 97.8% by adjusting the certainty grades of rules using genetic algorithms (GA). Synthetic signals were used as reference to further evaluate the performance of the elaborated classifier. The recognition of the isolated MUAPs depends largely on noise level and is acceptable down to the signal to noise ratio of 20 dB with a detection probability of 0.96. The recognition of overlapped MUAPs depends slightly on the noise level with a detection probability of about 0.8. The corresponding misrecognition is caused principally by the synchronisation and the small overlapping degree.     

The application of independent component analysis to the multi-channel surface electromyographic signals for separation of motor unit action potential trains: part I-measuring techniques.

The purpose of this study is to examine whether or not the application of independent component analysis (ICA) is useful for separation of motor unit action potential trains (MUAPTs) from the multi-channel surface EMG (sEMG) signals. In this study, the eight-channel sEMG signals were recorded from tibialis anterior muscles during isometric dorsi-flexions at 5%, 10%, 15% and 20% maximal voluntary contraction. Recording MUAP waveforms with little time delay mounted between the channels were obtained by vertical sEMG channel arrangements to muscle fibers. The independent components estimated by FastICA were compared with the sEMG signals and the principal components calculated by principal component analysis (PCA). From our results, it was shown that FastICA could separate groups of similar MUAP waveforms of the sEMG signals separated into each independent component while PCA could not sufficiently separate the groups into the principal components. A greater reduction of interferences between different MUAP waveforms was demonstrated by the use of FastICA. Therefore, it is suggested that FastICA could provide much better discrimination of the properties of MUAPTs for sEMG signal decomposition, i.e. waveforms, discharge intervals, etc., than not only PCA but also the original sEMG signals.     

Effect of aging on properties of motor unit action potentials in the rat medial gastrocnemius muscle

The purpose of this study was to investigate whether age-related changes in motor unit (MU) contractile properties are reflected in parameters of motor unit action potentials (MUAPs). MUs of the medial gastrocnemius muscle were functionally isolated in anaesthetized Wistar rats. A control group of young animals (5–10 mo) was compared to two groups of old rats (24–25 mo and 28–30 mo). The basic contractile properties of MUs as well as the amplitude, total duration, peak-to-peak time, and number of turns within MUAPs were measured. Effects of aging were mainly observed for fast fatigable MUs (a prolongation of MUAPs and increased number of turns). The MUAP amplitude did not change significantly with aging in either MU type, but it correlated to the twitch or tetanic forces, which tended to increase with age, especially for slow MUs. We concluded that the prolongation of MUAPs and the greater incidence of signal turns was probably a result of a decrease in muscle fiber conduction velocity and/or an increase in their dispersion, and enlargement of MU territories – presumably caused by axonal sprouting of surviving motoneurons. The latter might also be responsible for the observed age-related tendency for a increase in MUAP amplitudes in slow MUs.     

Rigorous performance assessment of the algorithms for resolving motor unit action potential superpositions

It is necessary to decompose the intra-muscular EMG signal to extract motor unit action potential (MUAP) waveforms and firing times. Some algorithms were proposed in the literature to resolve superimposed MUAPs, including Peel-Off (PO), branch and bound (BB), genetic algorithm (GA), and particle swarm optimization (PSO). This study aimed to compare these algorithms in terms of overall accuracy and running time. Two sets of two-to-five MUAP templates (set1: a wide range of energies, and set2: a high degree of similarity) were used. Such templates were time-shifted, and white Gaussian noise was added. A total of 1000 superpositions were simulated for each template and were resolved using PO (also, POI: interpolated PO), BB, GA, and PSO algorithms. The generalized estimating equation was used to identify which method significantly outperformed, while the overall rank product was used for overall ranking. The rankings were PSO, BB, GA, PO, and POI in the first, and BB, PSO, GA, PO, POI in the second set. The overall ranking was BB, PSO, GA, PO, and POI in the entire dataset. Although the BB algorithm is generally fast, there are cases where the BB algorithm is too slow and it is thus not suitable for real-time applications.     

Separation of propagating and non propagating components in surface EMG

Surface electromyogram (EMG) detected by electrode arrays along the muscle fibre direction can be approximated by the sum of propagating and non propagating components. A technique to separate propagating and non propagating components in surface EMG signals is developed. The first step is an adaptive filter, which allows obtaining an estimation of the delay between signals detected at different channels and a first estimate of propagating and non propagating components; the second step is used to optimise the estimation of the two components. The method is applicable to signals with one propagating and one non propagating component. It was optimised on simulated signals, and then applied to single motor unit action potentials (MUAP) and to electrically elicited EMG (M-waves).

The new method was first tested on phenomenological signals constituted by the sum of a propagating and a non propagating signal and then applied to simulated and experimental EMG signals. Simulated signals were generated by a cylindrical, layered volume conductor model. Experimental signals were monopolar surface EMG signals collected from the abductor pollicis brevis muscle and M-waves recorded during transcutaneous electrical stimulation of the biceps muscle. The technique may find different applications: in single motor unit (MU) studies (a) for decreasing the variability and bias of CV estimates due to the presence of the non propagating components, (b) for estimating automatically the length of the muscle fibres from only three detected channels and (c) for removal of the stimulation artifact M-waves.     

The need for computer-aided electromyography

The limits of quantitative manual electromyography (EMG) are discussed. The role of computers in EMG laboratories is at present to imitate or replace the physician performing the test, extending his memory, to gather valuable information which cannot be obtained in a conventional way or, on the contrary, to delete the redundant information. The difficulties in standardizing the parameters of a single motor unit action potential (MUAP) are mainly related to the complexity of the EMG signal and its variability, particularly in pathological states. A computer-aided quantification of interference pattern is presented. The novel methods of examination applied in computerized EMG laboratories are discussed. The scope and limits of computer-aided EMG should be taken into action potential (MUAP) are mainly related to the complexicity of the EMG be accepted are listed.     

Motor unit number estimation as a complementary test to routine electromyography in the diagnosis of amyotrophic lateral sclerosis

Electromyographic (EMG) abnormalities that reveal denervation and reinnervation caused by lower motor neuron degeneration do not reflect the number of motor units that determines muscle strength. Consequently, motor unit activity potential (MUAP) parameters do not reflect muscle dysfunction.The aim of the study was to compare the value of motor unit number estimation (MUNE) and MUAP parameters as indicators of clinical muscle dysfunction in patients with amyotrophic lateral sclerosis (ALS), and to analyze the role of MUNE as a supplement to the EMG criteria for the diagnosis of ALS.In 25 patients with ALS, MUNE by the multipoint incremental method in the abductor digiti minimi (ADM) and quantitative EMG in the first dorsal interosseous (FDI) were obtained. The Medical Research Council (MRC) scale was used to evaluate clinical muscle dysfunction. A strong correlation between the number of motor units evaluated by MUNE and ADM clinical function by the MRC scale was found (P < 0.001). An increased value of surface-detected single motor action potential was associated with a decreased MRC score for ADM (P < 0.1). No relation was found between MUAP parameters in FDI and MRC scores. Our data support the value of the MUNE method for the detection of motor unit loss in ALS, and it could be postulated that MUNE studies may be considered complementary tests for ALS in a future revision of ALS criteria.     

Measuring the accuracies of motor unit firing times and action potential waveforms derived from surface electromyographic decomposition

This study included spike trigger averaging (STA) procedures to examine the acceptability of the Precision Decomposition (PD) III derived motor unit action potential (MUAP) trains that met the >90% accuracy criteria from the reconstruct-and-test. MUs met the >90% accuracy criteria from the reconstruct-and-test with STA procedures then applied. Y-intercepts and slopes were calculated for the firing rate- and MUAP amplitude-recruitment threshold relationships. Gaussian noise (1% of the SD of the mean interspike interval) was added to the firing times with the changes in MUAPs quantified. A total of 455 MUs were decomposed with 155 MUs removed as a result of the reconstruct-and-test. Five additional MUs were excluded via the STA criteria. The MUAP waveforms deteriorated with the inclusion of Gaussian noise. There were differences in the derived action potentials amplitudes of higher-threshold MUs between the PD III algorithm and the STA procedure. There was excellent agreement among the slopes and y-intercepts between the relationships that included or excluded MUs that did not meet the STA criteria. There was good agreement between the MUAP amplitude-recruitment threshold relationships derived from the PD III and STA procedure. The addition of the STA procedures did not alter the MU-derived relationships.     

Comparative analysis of motor unit action potentials of the medial gastrocnemius muscle in cats and rats

Properties of motor unit action potentials (MUAPs) were compared for medial gastrocnemius (MG) motor units (MUs) in cats and rats. The experiments on functionally isolated MUs were performed under general anaesthesia, under comparable conditions (surgery, stimulating protocol and recording methods) for both species investigated. The proportions of motor units and contractile properties of the sample used in the study were consistent with previous studies performed on the MG muscle in both animal species, so comparisons of action potentials of individual types of MUs were acknowledged as fully reliable. The most prominent differences concerning MUAPs were observed in total duration and peak-to-peak times which for all MU types were about twice longer in cat MUs, in comparison to the rat MUs. The considerable disproportions were observed between the MUAP amplitudes of FF (fast fatigable), FR (fast resistant to fatigue) and S (slow) MUs in each species (the highest amplitudes were measured for FF and the lowest for S MUs), but there were no significant differences between cat and rat when respective types of MUs were compared. The shapes of MUAPs were commonly characterized by biphasic waveforms composed of two or three turns in all types of units, and no interspecies differences were revealed. Several factors influencing MUAP parameters were discussed indicating most of all importance of variable length of cat and rat muscle fibres and ambiguous influence of motor unit size, thickness of muscle fibres and their density around the recording electrode in the MG muscle of both species.