Low-frequency spectral analysis of the e.m.g

Single channel surface electromyogram deconvolution to explore motor unit discharges

Interference surface electromyogram (EMG) reflects many bioelectric properties of active motor units (MU), which are however difficult to estimate due to the asynchronous summation of their discharges. This paper introduces a deconvolution technique to estimate the cumulative firings of MUs. Tests in simulations show that the power spectral density of the estimated MU firings has a low-frequency peak corresponding to the mean firing rate of MUs in the detection volume of the recording system, weighted by the amplitudes of MU action potentials. The peak increases in amplitude and its centroid shifts to a higher frequency when MU synchronization is simulated (mainly due to the shift of discharges of large MUs). The peak is found even at high force levels, when such a contribution does not emerge from the EMG. This result is also confirmed in preliminary applications to experimental data. Moreover, the simulated cumulative firings of MUs are estimated with a correlation above 90% (considering frequency contributions up to 150 Hz), for all force levels. The method requires a single EMG channel, thus being feasible even in applied studies using simple recording systems. It may open many potential applications, e.g., in the study of the modulation of MU firing rate induced by either fatigue or pathology and in coherency analysis. Graphical Abstract Examples of application of the deconvolution (Deconv) algorithm and comparison with the cumulative firings and the cumulated weighted firings (CWF, i.e., each firing pattern is weighted by the root mean squared amplitude of the corresponding MU action potential). Portions of data are shown on the left, the power spectral densities (PSD) on the right (Welch method applied to 3 s of data, sub-epochs of 0.5 s, mean value removed from each of them, 50% of overlap). A) Simulated signal (50% of maximal voluntary contraction, MVC) with random MU firings. B) Simulated signal (50% MVC) with a level of synchronization equal to 10%. C) Experimental data from vastus medialis at 40% MVC (data decomposed by the algorithm of Holobar and Zazula, IEEE Trans. Sig. Proc. 2007; PSD of the cumulated firings almost identical to that of CWF, as few MUs were identified).

Increased intensity and reduced frequency of EMG signals from feline self-reinnervated ankle extensors during walking do not normalize excessive lengthening

Kinematics of cat level walking recover after elimination of length-dependent sensory feedback from the major ankle extensor muscles induced by self-reinnervation. Little is known, however, about changes in locomotor myoelectric activity of self-reinnervated muscles. We examined the myoelectric activity of self-reinnervated muscles and intact synergists to determine the extent to which patterns of muscle activity change as almost normal walking is restored following muscle self-reinnervation. Nerves to soleus (SO) and lateral gastrocnemius (LG) of six adult cats were surgically transected and repaired. Intramuscular myoelectric signals of SO, LG, medial gastrocnemius (MG), and plantaris (PL), muscle fascicle length of SO and MG, and hindlimb mechanics were recorded during level and slope (±27°) walking before and after (10-12 wk postsurgery) self-reinnervation of LG and SO. Mean myoelectric signal intensity and frequency were determined using wavelet analysis. Following SO and LG self-reinnervation, mean myoelectric signal intensity increased and frequency decreased in most conditions for SO and LG as well as for intact synergist MG (P < 0.05). Greater elongation of SO muscle-tendon unit during downslope and unchanged magnitudes of ankle extensor moment during the stance phase in all walking conditions suggested a functional deficiency of ankle extensors after self-reinnervation. Possible effects of morphological reorganization of motor units of ankle extensors and altered sensory and central inputs on the changes in myoelectric activity of self-reinnervated SO and LG are discussed.

Task-dependent activity of motor unit populations in feline ankle extensor muscles

Understanding the functional significance of the morphological diversity of mammalian skeletal muscles is limited by technical difficulties of estimating the contribution of motor units with different properties to unconstrained motor behaviours. Recently developed wavelet and principal components analysis of intramuscular myoelectric signals has linked signals with lower and higher frequency contents to the use of slower and faster motor unit populations. In this study we estimated the relative contributions of lower and higher frequency signals of cat ankle extensors (soleus, medial and lateral gastrocnemii, plantaris) during level, downslope and upslope walking and the paw-shake response. This was done using the first two myoelectric signal principal components (PCI, PCII), explaining over 90% of the signal, and an angle θ, a function of PCI/PCII, indicating the relative contribution of slower and faster motor unit populations. Mean myoelectric frequencies in all walking conditions were lowest for slow soleus (234 Hz) and highest for fast gastrocnemii (307 and 330 Hz) muscles. Motor unit populations within and across the studied muscles that demonstrated lower myoelectric frequency (suggesting slower populations) were recruited during tasks and movement phases with lower mechanical demands on the ankle extensors--during downslope and level walking and in early walking stance and paw-shake phases. With increasing mechanical demands (upslope walking, mid-phase of paw-shake cycles), motor unit populations generating higher frequency signals (suggesting faster populations) contributed progressively more. We conclude that the myoelectric frequency contents within and between feline ankle extensors vary across studied motor behaviours, with patterns that are generally consistent with muscle fibre-type composition.

Decreased motor unit discharge rate in the potentiated human tibialis anterior muscle

AIM

The purpose of this study was to examine the influence of post-activation potentiation (PAP), the transient increase in low-frequency isometric force observed after muscle activity, on motor unit discharge rates measured during submaximal contractions.

METHODS

A quadrifilar needle electrode was inserted into the tibialis anterior muscle to determine discharge rate of individual motor units while monopolar electrodes were used to monitor the root-mean-square (RMS) and mean power frequency (MPF) of the surface EMG signal. Control (unpotentiated) and experimental (potentiated) measures were obtained during a 5 s voluntary contraction at 50% of maximal. In between these measures, subjects performed a 10 s maximal voluntary contraction (MVC) to induce PAP.

RESULTS

All subjects data are reported as means ± SEM (n = 10). Compared to baseline values measured prior to the MVC, isometric twitch force measured immediately after the MVC was increased by 260 ± 16% (day 3). Motor unit discharge rate in the potentiated tibialis anterior muscle decreased by approx. 10%, from 20.3 ± 0.8 (before) to 18.3 ± 0.99 pps (P = 0.01) (after). Moreover, the MPF was decreased by approx. 9% (from 58.1 ± 2.84 to 53.6 ± 2.85 Hz; P = 0.01) in the potentiated tibialis anterior. On the other hand, consistent with the absence of fatigue during the MVC, the RMS signal was not altered in the potentiated tibialis anterior (0.29 ± 0.03 vs. 0.33 ± 0.04 mV; P = 0.07).

CONCLUSION

Motor unit discharge rates determined during a brief, submaximal contraction were decreased in the potentiated human tibialis anterior muscle.

Spatio-temporal representation of multichannel EMG firing patterns and its clinical applications

Analyzing motor unit (MU) activity is essential for studying the neurological dysfunction of upper motor neuron disorders (UMND). This study employs multichannel surface electromyographic (EMG) signals, as recorded from the upper arm during elbow flexion and extension, to analyze the temporal changes and spatial distribution of the dominant firing rate. To estimate the dominant firing rate, the autoregressive (AR) spectrum analysis method is utilized to detect the peaks and poles of the AR model, of the surface EMG spectrum below 40 Hz. The temporal changes in firing rates are also observed by using the spectrogram representation of low-frequency EMG spectra. The EMG spectrogram facilitates examination of the time-varying characteristics of firing rates and recruitment of MUs from surface EMG signal. The low-frequency spectra of multichannel EMG are then represented in a polar form to visualize the spatial distribution of firing patterns across muscles. Via spatio-temporal representation techniques, this study provides a viable approach of observing both the spatial and temporal patterns of MU activities in normal subjects and patients with UMND, including cerebrovascular disease and Parkinson's disease.

Single motor unit myoelectric signal analysis with nonstationary data

The information content of the myoelectric signal (MES) is commonly revealed by statistical measures in the time or frequency domain. Empirical analyses of the MES from a single motor unit have generally assumed that features are invariant with time. Theoretical and experimental work has been done to demonstrate how nonstationary behavior in the discharge statistics of a motor neuron may affect estimates of features extracted from the motor unit's contribution to the MES. Specifically, it has been shown that nonstationary behavior can markedly influence estimates of features describing motor neuron firing behavior and consequently, the low-frequency portion of the MES power spectral density. These results may help to explain the discrepancies in the literature which report empirical models of motor neuron firing statistics.

Changes in the myoelectric signal (MES) power spectra during dynamic contractions

Power spectra of surface myoelectric signals (MES) have been estimated during dynamic contractions of the biceps brachii. The effect of applied torque, velocity of contraction and muscle length have been studied with respect to median frequency and total spectral power. The results of the investigation carried out thus far can be summarised as follows: (1) The median frequency is observed to vary with the applied torque and the muscle length but not with the velocity of contraction. (2) The expected dependence of the MES power with the applied torque, velocity and muscle length was confirmed. These results suggest that the power spectrum could be used to reveal changes in the electrophysiological characteristics of a muscle during dynamic contractions as is the case with static contractions.

The fatigability of two agonistic muscles in human isometric voluntary submaximal contraction: an EMG study. II. Motor unit firing rate and recruitment

The recruitment and firing rate of biceps brachii (BB) and brachioradialis (BR) motor units (MUs) were studied in the course of fatiguing isometric contractions at 20%-30% of maximal voluntary contraction (MVC). MU recruitment generally occurred throughout the maintained contraction and was similar for BB and BR muscles. Newly recruited MUs started to discharge in the form of bursts, the duration of which increased until a continuous rhythmical firing was achieved. Within each burst, the first interval between two consecutive discharges was usually the shortest. MU threshold was lowered just after the limit time of the maintained contraction. The MU's firing rate either increased or remained stable as a function of the elapsed time. It is concluded that (1) in fatiguing isometric contractions at 20%-30% MVC contractile failure is mainly compensated for by MU recruitment and a lowered MU threshold and (2) differences between in surface changes in the electromyogram of BB and BR muscles cannot easily be explained by related differences in MU firing rate and recruitment.

The motor unit firing rate and the power spectrum of EMG in humans

The EMG power spectrum is influenced by many factors such as the conduction velocity of the muscle fiber, the action potential of the motor unit, the number of motor units firing near the electrode, and the recording conditions. Model studies of the relation between motor unit firing rate and power spectrum of EMG have produced conflicting results. To examine this relation in vivo the brachial biceps muscle was examined in 14 controls at a force of 10% of maximum. The motor unit firing intervals were obtained from 164 motor units, sampled with a single fiber electrode. The EMG was sampled at 10 sites in each muscle with a concentric electrode and the power spectrum was obtained using fast Fourier transformation. The mean power frequency of the interference pattern as well as the relative power at 1400 Hz both decreased with increasing motor unit firing intervals between subjects. The study thus indicates that the amount of high frequencies in the power spectrum is greater in a subject with a high firing rate of the motor units than in a subject with a low firing rate.

Spectral analysis of the surface electromyogram as a tool for studying rate modulation: a comparison between theory, simulation, and experiment

Theoretical work suggests that if the interpulse intervals ( IPIs ) of motor unit action potential trains ( MUAPTs ) are independently and normally distributed, then spectral analysis of the electromyogram could be a useful tool for studying rate modulation by virtue of the presence of a peak in the power spectrum at the average firing frequency of all active motor units. It is shown in this paper that IPIs need not be normally distributed, specifically that the results are very much the same if the IPIs are distributed according to a Gamma probability density function ( PDF ). Simulation of the electromyogram based on this theory proved the applicability of the method. Experimental results obtained for the masseter, biceps brachii and first dorsal interosseus (FDI) muscles, however, were in disagreement with both theory and simulation except for the biceps muscle at force levels up to 20% of the maximal force and for the masseter and FDI muscles in 1 out of 5 subjects. This indicates that the models for MUAPTs hitherto used might not be generally correct. Apart from this discrepancy, our results reveal differences between masseter and FDI muscles on the one hand and the biceps brachii on the other, which indicate that motor unit synchronisation is much more pronounced in the latter muscle.

Influence of motor unit firing statistics on the median frequency of the EMG power spectrum

Changes in the EMG power spectrum during static fatiguing contractions are often attributed to changes in muscle fibre action potential conduction velocity. Mathematical models of the EMG power spectrum, which have been empirically confirmed, predict that under certain conditions a distinct maximum occurs in the low-frequency part of the spectrum, indicating the dominant firing rate of the motor units. The present study investigated the influence of this firing rate peak on the spectral changes during a static fatiguing contraction at 50% of maximum EMG amplitude in the frontalis and corrugator supercilii muscles. An exponential decrease of the median frequency (MF) of the EMG power spectrum was observed when the firing rate peak was absent. When the firing rate peak was present, an exaggerated decrease of MF in the beginning of the contraction was found, which was associated with an increase in firing rate peak magnitude. In later stages of the contraction, a partial recovery of MF occurred, concomitant with a decrease in firing rate peak magnitude. The influence of the firing rate peak on MF was also investigated during nonfatiguing contractions of the frontalis muscle at 20, 40, 60, and 80% of maximum EMG amplitude. A curvilinear relationship between MF and contraction strength was found, whether firing rate peaks were present or absent. The presence of firing rate peaks, however, was associated with a decrease in MF which was inversely related to contraction strength, due to the inverse relationship between firing rate peak magnitude and contraction strength.