Motor unit recruitment strategies investigated by surface EMG variables.

During isometric contractions of increasing strength, motor units (MUs) are recruited by the central nervous system in an orderly manner starting with the smallest, with muscle fibers that usually show the lowest conduction velocity (CV). Theory predicts that the higher the velocity of propagation of the action potential, the higher the power at high frequencies of the detected surface signal. These considerations suggest that the power spectral density of the surface detected electromyogram (EMG) signal may give indications about the MU recruitment process. The purpose of this paper is to investigate the potential and limitations of spectral analysis of the surface EMG signal as a technique for the investigation of muscle force control. The study is based on a simulation approach and on an experimental investigation of the properties of surface EMG signals detected from the biceps brachii during isometric linearly increasing torque contractions. Both simulation and experimental data indicate that volume conductor properties play an important role as confounding factors that may mask any relation between EMG spectral variables and estimated CV as a size principle parameter during ramp contractions. The correlation between spectral variables and CV is thus significantly lower when the MU pool is not stable than during constant-torque isometric contractions. Our results do not support the establishment of a general relationship between spectral EMG variables and torque or recruitment strategy.     

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.     

Effect of synchronization of firings of different motor unit types on the force variability in a model of the rat medial gastrocnemius muscle

The synchronized firings of active motor units (MUs) increase the oscillations of muscle force, observed as physiological tremor. This study aimed to investigate the effects of synchronizing the firings within three types of MUs (slow—S, fast resistant to fatigue–FR, and fast fatigable–FF) on the muscle force production using a mathematical model of the rat medial gastrocnemius muscle. The model was designed based on the actual proportion and physiological properties of MUs and motoneurons innervating the muscle. The isometric muscle and MU forces were simulated by a model predicting non-synchronized firing of a pool of 57 MUs (including 8 S, 23 FR, and 26 FF) to ascertain a maximum excitatory signal when all MUs were recruited into the contraction. The mean firing frequency of each MU depended upon the twitch contraction time, whereas the recruitment order was determined according to increasing forces (the size principle). The synchronization of firings of individual MUs was simulated using four different modes and inducing the synchronization of firings within three time windows (± 2, ± 4, and ± 6 ms) for four different combinations of MUs. The synchronization was estimated using two parameters, the correlation coefficient and the cross-interval synchronization index. The four scenarios of synchronization increased the values of the root-mean-square, range, and maximum force in correlation with the increase of the time window. Greater synchronization index values resulted in higher root-mean-square, range, and maximum of force outcomes for all MU types as well as for the whole muscle output; however, the mean spectral frequency of the forces decreased, whereas the mean force remained nearly unchanged. The range of variability and the root-mean-square of forces were higher for fast MUs than for slow MUs; meanwhile, the relative values of these parameters were highest for slow MUs, indicating their important contribution to muscle tremor, especially during weak contractions.     

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.     

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.     

Muscular forearm activation in hand-grip tasks with superimposition of mechanical vibrations

The purpose of this paper is to evaluate the muscular activation of the forearm, with or without vibration stimuli at different frequencies while performing a grip tasks of 45 s at various level of exerted force. In 16 individuals, 9 females and 7 males, the surface electromyogram (EMG) of extensor carpi radialis longus and the flexor carpi ulnari muscles were assessed. At a short latency from onset EMG, RMS and the level of MU synchronization were assessed to evaluate the muscular adaptations. Whilst a trend of decay of EMG Median frequency (MDF_d) was employed as an index of muscular fatigue. Muscular tasks consists of the grip of an instrumented handle at a force level of 20%, 30%, 40%, 60% of the maximum voluntary force. Vibration was supplied by a shaker to the hand in mono-frequential waves at 20, 30, 33 and 40 Hz. In relation to EMG, RMS and MU synchronization, the muscular activation does not seem to change with the superimposition of the mechanical vibrations, on the contrary a lower MDF_d was observed at 33 Hz than in absence of vibration. This suggests an early muscular fatigue induced by vibration due to the fact that 33 Hz is a resonance frequency for the hand-arm system.     

Oxygen availability and motor unit activity in humans.

Six men were studied to determine the interrelationships among blood supply, motor unit (MU) activity and lactate concentrations during intermittent isometric contractions of the hand grip muscles. The subjects performed repeated contractions at 20% of maximal voluntary contraction (MVC) for 2 s followed by 2-s rest for 4 min with either unhindered blood circulation or arterial occlusion given between the 1st and 2nd min. The simultaneously recorded intramuscular MU spikes and surface electromyogram (EMG) data indicated that mean MU spike amplitude, firing frequency and the parameters of surface EMG power spectra (mean power frequency and root mean square amplitude) remained constant during the experiment with unhindered circulation, providing no electrophysiological signs of muscle fatigue. Significant increases in mean MU spike amplitude and frequency were, however, evident during the contractions with arterial occlusion. Similar patterns of significant changes in the surface EMG spectra parameters and venous lactate concentration were also observed, while the integrated force-time curves remained constant. These data would suggest that the metabolic state of the active muscles may have played an important role in the regulation of MU recruitment and rate coding patterns during exercise.     

Motor unit activity and surface electromyogram power spectrum during increasing force of contraction

Twelve male subjects were tested to determine the relationship between motor unit (MU) activities and surface electromyogram (EMG) power spectral parameters with contractions increasing linearly from zero to 80% of maximal voluntary contraction (MVC). Intramuscular spike and surface EMG signals recorded simultaneously from biceps brachii were analyzed by means of a computer-aided intramuscular MU spike amplitude-frequency (ISAF) histogram and an EMG frequency power spectral analysis. All measurements were made in triplicate and averaged. Results indicate that there were highly significant increases in surface EMG amplitude (71 +/- 31.3 to 505 +/- 188 microV, p less than 0.01) and mean power frequency (89 +/- 13.3 to 123 +/- 23.5 Hz, p less than 0.01) with increasing force. These changes were accompanied by progressive increases in the firing frequency of MU's initially recruited, and of newly recruited MU's with relatively larger spike amplitudes. The group data in the ISAF histograms revealed significant increases in mean spike amplitude (412 +/- 79 to 972 +/- 117 microV, p less than 0.01) and mean firing frequency (17.8 +/- 5.4 to 24.7 +/- 4.1 Hz, p less than 0.01). These data suggest that surface EMG spectral analysis can provide a sensitive measure of the relative changes in MU activity during increasing force output.     

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.     

Nonlinear surface EMG analysis to detect changes of motor unit conduction velocity and synchronization.

Amplitude and frequency content of the surface electromyographic (EMG) signal reflect central and peripheral modifications of the neuromuscular system. Classic surface EMG spectral variables applied to assess muscle functions are the centroid and median power spectral frequencies. More recently, nonlinear tools have been introduced to analyze the surface EMG; among them, the recurrence quantification analysis (RQA) was shown to be particularly promising for the detection of muscle status changes. The purpose of this work was to analyze the effect of motor unit short-term synchronization and conduction velocity (CV) on EMG spectral variables and two variables extracted by RQA, the percentage of recurrence (%Rec) and determinism (%Det). The study was performed on the basis of a simulation model, which allowed changing the degree of synchronization and mean CV of a number of motor units, and of an experimental investigation of the surface EMG signal properties detected during high-force-level isometric fatiguing contractions of the biceps brachii muscle. Simulations and experimental results were largely in agreement and show that 1) spectral variables, %Rec, and %Det are influenced by CV and degree of synchronization; 2) spectral variables are highly correlated with %Det (R = -0.95 in the simulations and -0.78 and -0.75 for the initial values and normalized slopes, respectively, in the experimental signals), and thus the information they provide on muscle properties is basically the same; and 3) variations of %Det and %Rec in response to changes in muscle properties are significantly larger than the variations of spectral variables. This study validates RQA as a means for fatigue assessment with potential advantages (such as the higher sensitivity to changes of muscle status) with respect to the classic spectral analysis.     

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.     

Assessment of low back muscle fatigue by surface EMG signal analysis: methodological aspects.

This paper focuses on methodological issues related to surface electromyographic (EMG) signal detection from the low back muscles. In particular, we analysed (1) the characteristics (in terms of propagating components) of the signals detected from these muscles; (2) the effect of electrode location on the variables extracted from surface EMG; (3) the effect of the inter-electrode distance (IED) on the same variables; (4) the possibility of assessing fatigue during high and very low force level contractions. To address these issues, we detected single differential surface EMG signals by arrays of eight electrodes from six locations on the two sides of the spine, at the levels of the first (L1), the second (L2), and the fifth (L5) lumbar vertebra. In total, 42 surface EMG channels were acquired at the same time during both high and low force, short and long duration contractions. The main results were: (1) signal quality is poor with predominance of non-travelling components; (2) as a consequence of point (1), in the majority of the cases it is not possible to reliably estimate muscle fiber conduction velocity; (3) despite the poor signal quality, it was possible to distinguish the fatigue properties of the investigated muscles and the fatigability at different contraction levels; (4) IED affects the sensitivity of surface EMG variables to electrode location and large IEDs are suggested when spectral and amplitude analysis is performed; (5) the sensitivity of surface EMG variables to changes in electrode location is on average larger than for other muscles with less complex architecture; (6) IED influences amplitude initial values and slopes, and spectral variable initial values; (7) normalized slopes for both amplitude and spectral variables are not affected by IED and, thus, are suggested for fatigue analysis at different postures or during movement, when IED may change in different conditions (in case of separated electrodes); (8) the surface EMG technique at the global level of amplitude and spectral analysis cannot be used to characterize fatigue properties of low back muscles during very low level, long duration contractions since in these cases the non-stable MU pool has a major influence on the EMG variables. These considerations clarify issues only partially investigated in past studies. The limitations indicated above are important and should be carefully discussed when presenting surface EMG results as a means for low back muscle assessment in clinical practice.     

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.     

Simulation analysis of interference EMG during fatiguing voluntary contractions. Part I: What do the intramuscular spike amplitude–frequency histograms reflect?

Decline in amplitude of EMG signals and in the rate of counts of intramuscularly recorded spikes during fatigue is often attributed to a progressive reduction of the neural drive only. As a rule, alterations in intracellular action potential (IAP) are not taken into account. To test correctness of the hypothesis, the effect of various discharge frequency patterns as well as changes in IAP shape and muscle fibre propagation velocity (MFPV) on the spike amplitude-frequency histogram of intramuscular interference EMG signals were simulated and analyzed. It was assumed that muscle was composed of four types of motor units (MUs): slow-twitch fatigue resistant, fast-twitch fatigue resistant, fast intermediate, and fast fatigable. MFPV and IAP duration at initial stage before fatigue as well as their changes differed for individual MU types. Fatigability of individual MU types in normal conditions as well as in the case of ischaemic or low oxygen conditions due to restricted blood flow was also taken into account. It was found that spike amplitude-frequency histogram is poorly sensitive to MU firing frequency, while it is highly sensitive to IAP profile lengthening. It is concluded that spike amplitude-frequency analysis can hardly provide a correct measure of MU rate-coding pattern during fatigue.     

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.     

Acoustic myography for investigating human skeletal muscle fatigue.

Sounds produced during voluntary isometric contractions of the quadriceps muscle were studied by acoustic myography (AMG) in five healthy adults. With the subject seated, isometric force, surface electromyography (EMG), and AMG were recorded over rectus femoris, and the EMG and AMG signals were integrated (IEMG and IAMG). Contractions lasting 5 s each were performed at 10, 25, 50, 60, 75, and 100% of maximum voluntary contraction (MVC) force. Fatigue was then induced by repeated voluntary contractions (10 s on, 10 s off) at 75% MVC until only 40% MVC could be sustained. After 15 min of rest, the different force levels were again tested in relation to the fresh MVC. Both before and after fatiguing activity the relationships between force and IEMG [r = 0.99 +/- 0.01 (SD), n = 10] and force and IAMG (r = 0.98 +/- 0.02) were linear. After activity, however, the slopes of the regression lines for force and IEMG increased (P less than 0.01) but those for force and IAMG remained the same (P greater than 0.05). The present results clarify the relationship between AMG and isometric force in fatigued muscle without the problem of fatigue-induced tremor, which hampered previous studies of prolonged activity. This study contributes to the validation of AMG and shows that it is a potentially useful method for noninvasive assessment of force production and fatigue. Further studies to establish the origin of AMG activity are required before AMG can be accepted for use in neuromuscular physiology or rehabilitation.     

Synchronization of human motor unit firing during voluntary contraction of muscle and the tonic vibration reflex

We constructed cross-correlograms (CCGs) of action potentials of pairs of motor units (MUs) of human soleus, triceps brachii, and the first and second dorsal interosseous muscles. During voluntary muscle contraction, a pronounced peak in the zero bin was found in 21 out of 126 pairs investigated with the aid of the CCG; this indicates that the number of coincidental firings exceeded chance. The width of the peak did not exceed 5 msec (synchronization for a brief interval, i.e., short-term synchronization). When motoneurons of the soleus muscle were activated by vibration, correlations were found in 12 out 89 pairs of MUs investigated. On the CCGs of action potentials of MU pairs in two muscles (the first and second dorsal interosseous muscles), such correlations were found in four out of 10 pairs investigated. In all of these cases, the ratio of above-change coincidences relative to the total number of MU discharges was small, from 3.0 to 6.1%. Synchronization within a brief time interval can be considered a result of simultaneous creation of EPSPs in motoneourons reached by endings of a single pre-motor nerve fiber. In some pairs of MUs, long term synchronization (clustering) occurred, ie., synchronization lasting several tens of milliseconds. The long-term synchronization can be considered a manifestation of fatigue accompanying tremor.Institute for Problems of Information Transmission, Russian Academy of Sciences, Moscow. Translated from Neirofiziologiya, Vol. 23, No. 6, pp. 691–698, November–December, 1991.     

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.