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.     

Clinical applications of high-density surface EMG: A systematic review

High density-surface EMG (HD-sEMG) is a non-invasive technique to measure electrical muscle activity with multiple (more than two) closely spaced electrodes overlying a restricted area of the skin. Besides temporal activity HD-sEMG also allows spatial EMG activity to be recorded, thus expanding the possibilities to detect new muscle characteristics. Especially muscle fiber conduction velocity (MFCV) measurements and the evaluation of single motor unit (MU) characteristics come into view. This systematic review of the literature evaluates the clinical applications of HD-sEMG. Although beyond the scope of the present review, the search yielded a large number of “non-clinical” papers demonstrating that a considarable amount of work has been done and that significant technical progress has been made concerning the feasibility and optimization of HD-sEMG techniques. Twenty-nine clinical studies and four reviews of clinical applications of HD-sEMG were considered. The clinical studies concerned muscle fatigue, motor neuron diseases (MND), neuropathies, myopathies (mainly in patients with channelopathies), spontaneous muscle activity and MU firing rates. In principle, HD-sEMG allows pathological changes at the MU level to be detected, especially changes in neurogenic disorders and channelopathies. We additionally discuss several bioengineering aspects and future clinical applications of the technique and provide recommendations for further development and implementation of HD-sEMG as a clinical diagnostic tool.     

Simulation analysis of interference EMG during fatiguing voluntary contractions. Part II--changes in amplitude and spectral characteristics.

Capabilities of amplitude and spectral methods for information extraction from interference EMG signals were assessed through simulation and preliminary experiment. Muscle was composed of 4 types of motor units (MUs). Different hypotheses on changes in firing frequency of individual MUs, intracellular action potential (IAP) and muscle fibre propagation velocity (MFPV) during fatigue were analyzed. It was found that changes in amplitude characteristics of interference signals (root mean square, RMS, or integrated rectified value, IEMG) detected by intramuscular and surface electrodes differed. RMS and IEMG of surface detected interference signals could increase even under MU firing rate reduction and without MU synchronisation. IAP profile lengthening can affect amplitude characteristics more significantly than MU firing frequency. Thus, an increase of interference EMG amplitude is unreliable to reflect changes in the neural drive. The ratio between EMG amplitude and contraction response can hardly characterise the so-called 'neuromuscular efficiency'. The recently proposed spectral fatigue indices can be used for quantification of interference EMG signals. The indices are practically insensitive to MU firing frequency. IAP profile lengthening and decrease in MFPV enhanced the index value, while recruitment of fast fatigable MUs reduced it. Sensitivity of the indices was higher than that of indices traditionally used.     

The linear electrode array: a useful tool with many applications.

In this review we describe the basic principles of operation of linear electrode arrays for the detection of surface EMG signals, together with their most relevant current applications. A linear array of electrodes is a system which detects surface EMG signals in a number of points located along a line. A spatial filter is usually placed in each point for signal detection, so that the recording of EMG signals with linear arrays corresponds to the sampling in one spatial direction of a spatially filtered version of the potential distribution over the skin. Linear arrays provide indications on motor unit (MU) anatomical properties, such as the locations of the innervation zones and tendons, and the fiber length. Such systems allow the investigation of the properties of the volume conductor and its effect on surface detected signals. Moreover, linear arrays allow to estimate muscle fiber conduction velocity with a very low standard deviation of estimation (of the order of 0.1-0.2 m/s), thus providing reliable indications on muscle fiber membrane properties and their changes in time (for example with fatigue or during treatment). Conduction velocity can be estimated from a signal epoch (global estimate) or at the single MU level. In the latter case, MU action potentials are identified from the interference EMG signals and conduction velocity is estimated for each detected potential. In this way it is possible, in certain conditions, to investigate single MU control and conduction properties with a completely non-invasive approach. Linear arrays provide valuable information on the neuromuscular system properties and appear to be promising tools for applied studies and clinical research.     

Investigation of motor unit recruitment during stimulated contractions of tibialis anterior muscle

This work investigated motor unit (MU) recruitment during transcutaneous electrical stimulation (TES) of the tibialis anterior (TA) muscle, using experimental and simulated data. Surface electromyogram (EMG) and torque were measured during electrically-elicited contractions at different current intensities, on eight healthy subjects.EMG detected during stimulation (M-wave) was simulated selecting the elicited MUs on the basis of: (a) the simulated current density distribution in the territory of each MU and (b) the excitation threshold characteristic of the MU. Exerted force was simulated by adding the contribution of each of the elicited MUs. The effects of different fat layer thickness (between 2 and 8 mm), different distributions of excitation thresholds (random excitation threshold, higher threshold for larger MUs or smaller MUs), and different MU distributions within the muscle (random distribution, larger MU deeper in the muscle, smaller MU deeper) on EMG variables and torque were tested.Increase of the current intensity led to a first rapid increase of experimental M-wave amplitude, followed by a plateau. Further increases of the stimulation current determined an increase of the exerted force, without relevant changes of the M-wave. Similar results were obtained in simulations.Rate of change of conduction velocity (CV) and leading coefficient of the second order polynomial interpolating the force vs. stimulation level curve were estimated as a function of increasing current amplitudes. Experimental data showed an increase of estimated CV with increasing levels of the stimulation current (for all subjects) and a positive leading coefficient of force vs. stimulation current curve (for five of eight subjects). Simulations matched the experimental results only when larger MUs were preferably located deeper in the TA muscle (in line with a histochemical study). Marginal effect of MU excitation thresholds was observed, suggesting that MUs closer to the stimulation electrode are recruited first during TES regardless of their excitability.     

Non-invasive characterization of single motor unit electromyographic and mechanomyographic activities in the biceps brachii muscle.

The aim of the study was to investigate amplitude and frequency content of single motor unit (MU) electromyographic (EMG) and mechanomyographic (MMG) responses. Multi-channel surface EMG and MMG signals were detected from the dominant biceps brachii muscle of 10 volunteers during isometric voluntary contractions at 20%, 50%, and 80% of the maximal voluntary contraction (MVC) force. Each contraction was performed three times in the experimental session which was repeated in three non-consecutive days. Single MU action potentials were identified from the surface EMG signals and their times of occurrence used to trigger the averaging of the MMG signal. At each contraction level, the MUs with action potentials of highest amplitude were identified. Single MU EMG and MMG amplitude and mean frequency were estimated with normalized standard error of the mean within subjects (due to repetition of the measure in different trials and experimental sessions) smaller than 15% and 7%, respectively, in all conditions. The amplitude of the action potentials of the detected MUs increased with increasing force (mean +/- SD, 244 +/- 116 microV at 20% MVC, and 1426 +/- 638 microV at 80% MVC; P < 0.001) while MU MMG amplitude increased from 20% to 50% MVC (40.5 +/- 20.9 and 150 +/- 88.4 mm/s(2), respectively; P<0.001) and did not change significantly between 50% and 80% MVC (129 +/ -82.7 mm/s(2) at 80% MVC). MU EMG mean frequency decreased with contraction level (20% MVC: 97.2 +/- 13.9 Hz; 80% MVC: 86.2 +/- 11.4 Hz; P < 0.001) while MU MMG mean frequency increased (20% MVC: 33.2 +/- 6.8 Hz; 80% MVC: 40.1 +/- 6.1 Hz; P < 0.001). EMG peak-to-peak amplitude and mean frequency of individual MUs were not correlated with the corresponding variables of MMG at any contraction level.     

Human motor unit recruitment and derecruitment during long lasting intermittent contractions.

Seven healthy subjects were investigated in cyclic ramp-and-hold long lasting isometric contractions. Wire branched electrodes were used for selective recording of single motor unit (MU) potentials from m. biceps brachii. MU behaviour was defined in terms of recruitment/derecruitment thresholds (RT and DT) and the duration of interspike intervals (ISI). A total of 63 MUs was investigated: 40 units were active from the beginning of the task performance and another 23 were recruited later. There were no changes in the recruitment pattern of MUs with fatigue development - a short first ISI followed by a very long second one and an almost constant firing rate after this transient phase. The tendency of RT to gradually decrease dominates the results. Thus, the required constant rate of force increase with fatigue development was maintained mostly by the mechanisms of space coding (i.e., decrease of RT and recruitment of additional MUs). Oppositely, the time behaviour of the DT changes was not uniform and rate coding was an essential mechanism in the adaptation of MU activity to muscle fatigue during relaxation phases. The recruitment pattern and fatigue related behaviour of the additionally recruited MUs were similar to those of MUs active from the first cycle of the motor task performance.     

Intramuscular and surface electromyogram changes during muscle fatigue

Twelve male subjects were tested to determine the effects of motor unit (MU) recruitment and firing frequency on the surface electromyogram (EMG) frequency power spectra during sustained maximal voluntary contraction (MVC) and 50% MVC of the biceps brachii muscle. Both the intramuscular MU spikes and surface EMG were recorded simultaneously and analyzed by means of a computer-aided intramuscular spike amplitude-frequency histogram and frequency power spectral analysis, respectively. Results indicated that both mean power frequency (MPF) and amplitude (rmsEMG) of the surface EMG fell significantly (P less than 0.001) together with a progressive reduction in MU spike amplitude and firing frequency during sustained MVC. During 50% MVC there was a significant decline in MPF (P less than 0.001), but this decline was accompanied by a significant increase in rmsEMG (P less than 0.001) and a progressive MU recruitment as evidenced by an increased number of MUs with relatively large spike amplitude. Our data suggest that the surface EMG amplitude could better represent the underlying MU activity during muscle fatigue and the frequency powers spectral shift may or may not reflect changes in MU recruitment and rate-coding patterns.     

Continuous, intermitted and sporadic motor unit activity in the trapezius muscle during prolonged computer work.

The Cinderella hypothesis postulates the continuous activity of specific motor units (MUs) during low-level muscle contraction. The MUs may become metabolically overloaded, with the subject developing muscle pain and strain. The hypothesis requires MUs that are active for a time long enough to actually damage muscle fibers. The aim of this study was to determine if there are continuously active MUs in the right trapezius muscle during normal computer work using a computer mouse. Fourteen healthy subjects executed an interactive computer-learning program (ErgoLight) for 30 min. Six-channel intramuscular EMG and two-channel surface EMG signals were recorded from two positions of the trapezius muscle. Decomposition was achieved with automated, multi-channel, long-term decomposition software (EMG-LODEC). In two out of the 14 subjects, three MUs were continuously active throughout the 30 min. Although the majority of the MUs were active during only part of the experimental session, an ordered on-off behavior (e.g. substitution) pattern was not observed. As long-lasting activity was verified in some subjects, the results support the Cinderella hypothesis. However, it cannot be concluded here how long the MUs could stay active. If continuous activity overloads low threshold MUs, the potential exists for selective fibre injuries in low threshold MUs of the trapezius muscle in subjects exposed to long-term computer work.     

A study of the electromyogram using a population stochastic model of skeletal muscle

The features of the electromyogram (EMG) are studied using a population model of skeletal muscle based on the differing properties and the independent activation of motor units (MUs). It is shown, both analytically and by computer simulation, that: (a) The power spectrum of the EMG is determined by the distribution of filtering and firing properties of the active MUs. (b) A tendency towards a rhythmical grouping of action potentials is to be expected from a set of asynchronous MUs firing semiregularly at similar rates; the grouped electrical activity has a phase-lead over the force output of the set of about 180°. A unified explanation of the properties of the muscle force waveform and the electromyogram, in terms of asynchronous activity of MUs, is proposed. The explanation covers the relationship and the differences between the two signals.     

Sex differences in the detection of motor unit action potentials identified using high-density surface electromyography

Sex-related disparities in force production of humans have been widely observed. Previous literature has attributed differences in peripheral traits, such as muscle size, to explain these disparities. However, less is known about potential sex-related differences in central neuromuscular traits and many comparable studies, not exploring sex-related differences, exhibit a selection-bias in the recruitment of subjects making the generalization of their findings difficult. Utilizing high-density electromyography arrays and motor unit (MU) decomposition, the aim of the current study is to compare MU yield and discharge properties of the tibialis anterior between male and female humans. Twenty-four subjects (10 females) performed two submaximal (20%) isometric dorsiflexion contractions. On average, males yielded nearly twice the amount of MUs as females. Further, females had significantly higher MU discharge rate, lower MU action potential amplitude, and lower MU action potential frequency content than males despite similar levels of torque and MU discharge variability. These findings suggest differences in central neuromuscular control of force production between sexes; however, it is unclear how lower yield counts affect the accuracy of these results.     

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.     

Impairment of human soleus motor units during ischemia

It is generally accepted that ischemia produced by limb compression affects rapidly conducting large-diameter Ia afferents in the early stage and that the motor nerve-muscle complex is blocked later. This notion, however, seems to be controversial for several reasons, so an attempt to reveal the amount of motor unit (MU) impairment during ischemia was made. Observation of human soleus muscle electromyographic (EMG) signal recorded either by bipolar needle electrode or by surface electrodes at various levels of voluntary contraction during the course of ischemia showed that low-threshold small MUs were affected first while high-threshold large MUs survived longer. The changes in EMG patterns were temporally correlated with T-reflex deterioration. It is suggested that the early loss of low-threshold MUs may play a definite role in alterations of reflexes during ischemia.     

Non-invasive detection of the single motor unit action potential by averaging the spatial potential distribution triggered on a spatially filtered motor unit action potential.

For research as well as diagnostic applications the non-invasive detection of the activity of single motor units is of interest. The most direct information is expected to be found in monopolarly recorded data. But when an array of surface electrodes is used for the monopolar recordings of the potential distribution on the skin, in most cases an additional invasive needle electrode is utilized to detect the exact points in time when a certain motor unit is firing. With this supplementary information, an averaging of the monopolar EMG tracings can be performed. In this paper, a completely non-invasive methodology is presented which replaces the invasive needle by a spatial filtering procedure. The EMG signals from the m. biceps brachii are recorded monopolarly with an electrode array. Afterwards, a spatial filtering procedure, called normal double differentiating filter, is applied to the data. The EMG signals obtained are investigated by means of an amplitude threshold to distinguish the activity of different motor units. The point of the maximum amplitude of the selected peaks then is used as trigger point to average the monopolar EMG data. The time courses of the motor unit action potential signals found after applying the described procedure show similar shapes, while two different components are to be identified: corresponding to the spread of the excitation, one is referring to stationary, the other to travelling events. These results justify the possibility to replace the needle electrode to obtain a trigger event in the future by the non-invasive spatial filtering procedure.     

Experimental and modelling investigation of learning a fast elbow flexion in the horizontal plane

Changes in the kinematic and electromyographic characteristics that occur while learning to move as fast as possible have been studied experimentally. Experimental investigation of what happens to the individual motor units (MUs) is more difficult. Access to each MU is impossible, and the recruitment and force developing properties of all individual MUs cannot be known. Thus, what is currently known about MU firing is based on experiments that have recorded relatively few MUs compared to what exists in the entire muscle. A recently developed muscle model (Raikova and Aladjov, 2002, J. Biomechanics, 35, 1123-1135) composed of MUs with different properties can be used for such investigation. The process of learning fast elbow flexion in the horizontal plane was simulated and the results were compared with experimentally measured data. Comparing the simulation results of the very first trial of a particular subject with those of the last trail (at the end of the learning process), it can be concluded that the speed of limb motion and muscle forces increase initially as a result of the more synchronous MUs activation and the increase of firing rate of active MUs. Further improvement necessitated an appreciable reduction in the motor task requirements (i.e. less muscle force and less MUs' activity) set in the computational algorithm by optimization criteria. This forced the next process-inclusion of additional MUs.     

Selective fatigue of fast motor units after electrically elicited muscle contractions.

The aim of the present study was to elucidate the electrophysiological manifestations of selective fast motor unit (MU) activation by electrical stimulation (ES) of knee extensor muscles. In six male subjects, test contraction measurement at 40% maximal voluntary contraction (MVC) was performed before and at every 5 min (5, 10, 15 and 20 min) during 20-min low intensity intermittent exercise of either ES or voluntary contractions (VC) at 10% MVC (5-s isometric contraction and 5-s rest cycles). Both isolated intramuscular MU spikes obtained from three sets of bipolar fine-wire electrodes and surface electromyogram (EMG) were simultaneously recorded and were analyzed by means of a computer-aided intramuscular spike amplitude-frequency analysis and frequency power spectral analysis, respectively. Results indicated that mean MU spike amplitude, particularly those MUs with relatively large amplitude, was significantly reduced while those MUs with small spike amplitude increased their firing rate during the 40% MVC test contraction after the ES. This was accompanied by the increased amplitude of surface EMG (rmsEMG). However, no such significant changes in the intramuscular and surface EMGs were observed after VC. These findings indicated differential MU activation patterns in terms of MU recruitment and rate coding characteristics during ES and VC, respectively. Our data strongly suggest the possibility of "an inverse size principle" of MU recruitment during ES.     

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.     

1998 ISEK Congress Keynote Lecture: The use of electromyography in applied physiology

Electromyogram (EMG) analyses (surface, intramuscular and evoked potentials) in studies of muscle function have attracted increasing attention during recent years and have been applied to assess muscle endurance capacity, anaerobic and lactate thresholds, muscle biomechanics, motor learning, neuromuscular relaxation, optimal walking and pedalling speeds, muscle soreness, neuromuscular diseases, motor unit (MU) activities (MU recruitment and rate coding), and skeletal muscle fatigue. This paper deals with the use of EMG analyses employed in the area of applied physiology and is divided into three sections: surface EMG analyses; intramuscular EMG analyses; and evoked potential analyses.     

Modelling motor units in 3D: influence on muscle contraction and joint force via a proof of concept simulation

Functional heterogeneity is a skeletal muscle’s ability to generate diverse force vectors through localised motor unit (MU) recruitment. Existing 3D macroscopic continuum-mechanical finite element (FE) muscle models neglect MU anatomy and recruit muscle volume simultaneously, making them unsuitable for studying functional heterogeneity. Here, we develop a method to incorporate MU anatomy and information in 3D models. Virtual fibres in the muscle are grouped into MUs via a novel “virtual innervation” technique, which can control the units’ size, shape, position, and overlap. The discrete MU anatomy is then mapped to the FE mesh via statistical averaging, resulting in a volumetric MU distribution. Mesh dependency is investigated using a 2D idealised model and revealed that the amount of MU overlap is inversely proportional to mesh dependency. Simultaneous recruitment of a MU’s volume implies that action potentials (AP) propagate instantaneously. A 3D idealised model is used to verify this assumption, revealing that neglecting AP propagation results in a slightly less-steady force, advanced in time by approximately 20 ms, at the tendons. Lastly, the method is applied to a 3D, anatomically realistic model of the masticatory system to demonstrate the functional heterogeneity of masseter muscles in producing bite force. We found that the MU anatomy significantly affected bite force direction compared to bite force magnitude. MU position was much more efficacious in bringing about bite force changes than MU overlap. These results highlight the relevance of MU anatomy to muscle function and joint force, particularly for muscles with complex neuromuscular architecture.

     

Neuromuscular regulation and metabolism during exercise

This article reviews current evidence regarding neuromuscular regulation and metabolism during exercise. Particular emphases are given on the relationship between motor unit (MU) activity, including single MU analysis results and spinal alpha-motoneuron excitability, and cardio-respiratory response and blood lactate during dynamic exercise. In addition, a close physiological link between muscle energy metabolism and excitation-contraction processes (failure of one will affect the extent of the other) is summarized in the light of recent nuclear magnetic resonance (NMR) studies and results of neuromuscular disorder patients.