Influence of muscle fibre shortening on estimates of conduction velocity and spectral frequencies from surface electromyographic signals

Effects of muscle shortening on single-fiber, motor unit, and compound muscle action potentials

Even under isometric conditions, muscle contractions are associated with some degree of fiber shortening. The effects of muscle shortening on extracellular electromyographic potentials have not been characterized in detail. Moreover, the anatomical, biophysical, and detection factors influencing the muscle-shortening effects have been neither identified nor understood completely. Herein, we investigated the effects of muscle shortening on the amplitude and duration characteristics of single-fiber, motor unit, and compound muscle action potentials. We found that, at the single-fiber level, two main factors influenced the muscle-shortening effects: (1) the electrode position and distance relative to the myotendinous zone and (2) the electrode distance to the maxima of the dipole field arising from the stationary dipole created at the fiber-tendon junction. Besides, at the motor unit and muscle level, two additional factors were involved: (3) the overlapping between the propagating component of some fibers with the non-propagating component of other fibers and (4) the spatial spreading of the fiber-tendon junctions. The muscle-shortening effects depend critically on the electrode longitudinal distance to the myotendinous zone. When the electrode was placed far from the myotendinous zone, muscle shortening resulted in an enlargement and narrowing of the final (negative) phase of the potential, and this enlargement became less pronounced as the electrode approached the fiber endings. For electrode locations close to the myotendinous zone, muscle shortening caused a depression of both the main (positive) and final (negative) phases of the potential. Beyond the myotendinous zone, muscle shortening led to a decrease of the final (positive) phase. The present results provide reference information that will help to identify changes in MUPs and M waves due to muscle shortening, and thus to differentiate these changes from those caused by muscle fatigue.

M-wave potentiation after voluntary contractions of different durations and intensities in the tibialis anterior

The study was undertaken to provide insight into the mechanisms underlying the potentiation of the muscle compound action potential (M wave) after conditioning contractions. M waves were evoked in the tibialis anterior before and after isometric maximal voluntary contractions (MVC) of 1, 3, 6, 10, 30, and 60 s, and after 3-s contractions at 10, 30, 50, 70, 90, and 100% MVC. The amplitude, duration, and area of the first and second phases of the M wave, together with the median frequency (Fmedian) and muscle fiber conduction velocity (MFCV) were recorded. Furthermore, twitch force, muscle fascicle length, and pennation angle were measured at rest, before, and 1 s after the conditioning contractions. The results indicate that only the amplitude of the second phase of the M wave was significantly increased after conditioning contractions. The extent of this potentiation was similar for MVC durations ranging from 1 to 10 s and augmented progressively with contraction intensity from 30 to 70% MVC. After these conditioning contractions, the duration and area of the two M-wave phases decreased ( P < 0.05), whereas MFCV and Fmedian increased ( P < 0.05). For all of these parameters, the greatest changes occurred 1 s after the conditioning contraction. Changes in MFCV after the contractions were correlated with those in M-wave second-phase amplitude ( r2 = 0.42; P < 0.05) and Fmedian ( r2 = 0.53; P < 0.05). In contrast, fascicle length and pennation angle did not change after the conditioning contractions. It is concluded that the potentiation of the second phase of the M wave is mainly due to an increased MFCV.

Power spectral changes of the superimposed M wave during isometric voluntary contractions of increasing strength

INTRODUCTION

We examined the power spectral changes of the compound muscle action potential (M wave) evoked during isometric contractions of increasing strength.

METHODS

Surface electromyography (sEMG) of the vastus lateralis and medialis was recorded from 20 volunteers who performed 4-s step-wise isometric contractions of different intensities. A maximal M wave was elicited by a single stimulus to the femoral nerve and superimposed on the voluntary contractions. The spectral characteristics (Fmean and Fmedian) of sEMG and M-wave signals were calculated.

RESULTS

M-wave spectral indicators increased systematically with contraction intensity up to 60% MVC and then leveled off at higher forces. Over the 10-60% MVC range, the increase in spectral indicators was 3 times higher for M waves (36%) than for sEMG (12%).

CONCLUSIONS

The consistent increase in M-wave spectral characteristics with force is due to the fact that the number of motor units recruited by the superimposed supramaximal stimulus is approximately stable.

Comparison of the power spectral changes of the voluntary surface electromyogram and M wave during intermittent maximal voluntary contractions

INTRODUCTION

To compare the power spectral changes of the voluntary surface electromyogram (sEMG) and of the compound action potential (M wave) in the vastus medialis and vastus lateralis muscles during fatiguing contractions.

METHODS

Interference sEMG and force were recorded during 48 intermittent 3-s isometric maximal voluntary contractions (MVC) from 13 young, healthy subjects. M waves and twitches were evoked using supramaximal femoral nerve stimulation between the successive MVCs. Mean frequency (F mean), and median frequency were calculated from the sEMG and M waves. Muscle fiber conduction velocity (MFCV) was computed by cross-correlation.

RESULTS

The power spectral shift to lower frequencies was significantly greater for the voluntary sEMG than for the M waves (P < 0.05). Over the fatiguing protocol, the overall average decrease in MFCV (~25%) was comparable to that of sEMG F mean (~22%), but significantly greater than that of M-wave F mean (~9%) (P < 0.001). The mean decline in MFCV was highly correlated with the mean decreases in both sEMG and M-wave F mean.

CONCLUSIONS

The present findings indicated that, as fatigue progressed, central mechanisms could enhance the relative weight of the low-frequency components of the voluntary sEMG power spectrum, and/or the end-of-fiber (non-propagating) components could reduce the sensitivity of the M-wave spectrum to changes in conduction velocity.

The influence of wrist posture on the time and frequency EMG signal measures of forearm muscles

This study investigates how altering wrist posture influences the relationship between the time and frequency measures of the electromyography (EMG) signal of extensor digitorum communis (EDC) and flexor carpi ulnaris (FCU). Thirteen participants exerted handgrip force related to maximum voluntary contraction (MVC) in four tests: 20%MVC and 50%MVC in neutral wrist posture and 20%MVC in full wrist flexion and extension. EMG measurements from EDC and FCU were used to calculate normalized values of amplitude (nRMS) and mean and median frequency of the power spectrum (nMPF, nMF). During muscle shortening (wrist flexion for FCU and wrist extension for EDC) nRMS was approximately twofold higher than in neutral posture for FCU and fourfold for EDC. All measures obtained at 20%MVC in neutral posture were significantly different from 20%MVC in wrist flexion for FCU and 20%MVC in wrist extension for EDC (p<0.05). Differences between 50%MVC and 20%MVC at neutral posture (nRMS) were significant for both muscles, although in nMPF and nMF for EDC only. Muscle shortening changed the pattern of statistical significance when the time and frequency domain measures were compared, whereas muscle lengthening did not. It can be concluded that muscle shortening caused by altering wrist posture influences the relationship between the time and frequency measures in both muscles. This suggests that in studies using EMG in different wrist postures, changes in the relationship between the time and the frequency measures should be considered.

Effect of membrane properties on skeletal muscle fiber excitability: a sensitivity analysis

In this study, the sensitivity of skeletal muscle fiber excitability to changes in temperature and a range of geometrical, electrical and ionic membrane properties was examined using model simulation. A mathematical model of the propagating muscle fiber action potential (AP) was used to simulate muscle fiber APs while changing individual muscle fiber parameters in isolation to examine how they affect muscle fiber AP amplitude, shape and conduction velocity (CV). The behavior of the model was verified by comparison with previously reported experimental data from both in vivo studies conducted at physiological temperatures and in vitro and in silico studies conducted at lower temperatures. The simulation results presented demonstrate the sensitivity of AP amplitude, shape and CV and, therefore, muscle fiber excitability to small changes in a wide range of different muscle fiber parameters. Furthermore, they demonstrate the potential of computational modeling as a tool for investigating the underlying mechanisms of complex phenomena such as those which govern skeletal muscle excitation.

Spectral properties of electromyographic and mechanomyographic signals during dynamic concentric and eccentric contractions of the human biceps brachii muscle

The purpose of this study was to describe and examine the variations in recruitment patterns of motor units (MUs) in biceps brachii (BB) through a range of joint motion during dynamic eccentric and concentric contractions. Twelve healthy participants (6 females, 6 males, age=30±8.5 years) performed concentric and eccentric contractions with constant external loading at different levels. Surface electromyography (EMG) and mechanomyography (MMG) were recorded from BB. The EMGs and MMGs were decomposed into their intensities in time-frequency space using a wavelet technique. The EMG and MMG spectra were then compared using principal component analysis. Variations in total intensity, first principal component (PCI), and the angle θ formed by first component (PCI) and second component (PCII) loading scores were explained in terms of MU recruitment patterns and elbow angles. Elbow angle had a significant effect on dynamic concentric and eccentric contractions. The EMG total intensity was greater for concentric than for eccentric contractions in the present study. MMG total intensity, however, was lower during concentric than during eccentric contractions. In addition, there was no significant difference in θ between concentric and eccentric contractions for both EMG and MMG. Selective recruitment of fast MUs from BB muscle during eccentric muscle contractions was not found in the present study.

Methodology of electromyographic analysis of the trunk muscles during walking in healthy subjects: a literature review

PURPOSE

To review and discuss the literature about the use of trunk muscle electromyography - including the use of surface or fine-wire electrodes, site of application and muscle selection - during gait analysis in healthy subjects.

METHODS

The databases Pubmed, Web of Knowledge and Cochrane Library were searched. Articles were included when EMG activity of at least one trunk muscle was measured in healthy subjects during walking.

RESULTS

In the 33 selected articles 491 healthy subjects walked with different velocities on a treadmill and/or overground. The activity of the M. erector spinae, M. multifidus, M. obliquus externus and internus, M. rectus abdominus, M. trapezius, M. latissimus dorsi, M. transversus abdominus, M. iliopsoas and M. quadrates lumborum was measured. Twenty-nine studies used surface electrodes, one study fine-wire electrodes, and the other three studies used a combination. There is no consensus on the exact placement site of the electrodes.

CONCLUSION

Surface electrodes were used more often than fine-wire electrodes and the descriptions of the electrode locations were mostly vague and not consistent among the different studies. There is need for further research to make specific recommendations about the type of electrodes in combination with the optimal locations of application of these electrodes.

Effect of elbow joint angle on force–EMG relationships in human elbow flexor and extensor muscles

The purpose of this study was to examine the effect of joint angle on the relationship between force and electromyogram (EMG) amplitude and median frequency, in the biceps, brachioradialis and triceps muscles. Surface EMG were measured at eight elbow angles, during isometric flexion and extension at force levels from 10% to 100% of maximum voluntary contraction (MVC). Joint angle had a significant effect on MVC force, but not on MVC EMG amplitude in all of the muscles examined. The median frequency of the biceps and triceps EMG decreased with increasing muscle length, possibly due to relative changes in electrode position or a decrease in muscle fibre diameter. The relationship between EMG amplitude and force, normalised with respect to its maximum force at each angle, did not vary with joint angle in the biceps or brachioradialis muscles over all angles, or in the triceps between 45 degrees and 120 degrees of flexion. These results suggest that the neural excitation level to each muscle is determined by the required percentage of available force rather than the absolute force required. It is, therefore, recommended that when using surface EMG to estimate muscle excitation, force should be normalised with respect to its maximum value at each angle.

Surface EMG analysis on normal subjects based on isometric voluntary contraction

The objective of this study was to compute reference SEMG values for normal subjects of 13 parameters extracted in the time, frequency and bispectrum domain, from the Biceps Brachii (BB) muscle generated under isometric voluntary contraction (IVC). SEMG signals were recorded from 94 subjects for 5s at 10, 30, 50, 70 and 100% of maximum voluntary contraction (MVC). The Wilcoxon signed rank test was applied to detect significant differences or not at p<0.05 between force levels for each of the 13 parameters. The main findings of this study can be summarized as follows: (i) The time domain parameters turns per second and number of zero crossings per second increase significantly with force level. (ii) The power spectrum median frequency parameter decreases significantly with force level, whereas maximum power and total power increase significantly with force level. (iii) The bispectrum parameter, maximum amplitude, increases significantly with force level with the exception the transition from 30% to 50% MVC. Although, the tests for Gaussianity and linearity show no significant difference with force level, the SEMG signal exhibits a more Gaussian distribution with increase of force up to 70% MVC. The SEMG linearity test, which is a measure of how constant the bicoherence index is in the bi-frequency domain, shows that the signal's bicoherence index is less constant (hence, the signal is less linear) at 70% of MVC compared to 10, 30, 50 and 100% MVC. (iv) The time domain parameters have good correlation between them as well as, between each one of them and maximum and total power. The median frequency has a good (negative) correlation with the bispectrum peak amplitude. (v) No significant differences exist between values based on gender or age. The findings of this study can further be used for the assessment of subjects suffering with neuromuscular disorders, or in the rehabilitation laboratory for monitoring the elderly or the disabled, or in the occupational medicine laboratory.

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 considerable 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.

Electromyographic instantaneous amplitude and instantaneous mean power frequency patterns across a range of motion during a concentric isokinetic muscle action of the biceps brachii

The purpose of this study was to examine the electromyographic (EMG) instantaneous amplitude (IA) and instantaneous mean power frequency (IMPF) patterns for the biceps brachii muscle across a range of motion during maximal and submaximal concentric isokinetic muscle actions of the forearm flexors. Ten adults (mean +/- SD age = 22.0 +/- 3.4 years) performed a maximal and a submaximal [20% peak torque (PT)] concentric isokinetic forearm flexion muscle action at a velocity of 30 degrees s(-1). The surface EMG signal was detected from the biceps brachii muscle with a bipolar electrode arrangement, and the EMG IA and IMPF versus time relationships were examined for each subject using first- and second-order polynomial regression models. The results indicated that there were no consistent patterns between subjects for EMG IA or IMPF with increases in torque across the range of motion. Some of the potential nonphysiological factors that could influence the amplitude and/or frequency contents of the surface EMG signal during a dynamic muscle action include movement of the muscle fibers and innervation zone beneath the skin surface, as well as changes in muscle fiber length and the thickness of the tissue layer between the muscle and the recording electrodes. These factors may affect the EMG IA and IMPF patterns differently for each subject, thereby increasing the difficulty of drawing any general conclusions regarding the motor control strategies that increase torque across a range of motion.

A Finite Element Model for Describing the Effect of Muscle Shortening on Surface EMG

A finite-element model for the generation of single fiber action potentials in a muscle undergoing various degrees of fiber shortening is developed. The muscle is assumed fusiform with muscle fibers following a curvilinear path described by a Gaussian function. Different degrees of fiber shortening are simulated by changing the parameters of the fiber path and maintaining the volume of the muscle constant. The conductivity tensor is adapted to the muscle fiber orientation. In each point of the volume conductor, the conductivity of the muscle tissue in the direction of the fiber is larger than that in the transversal direction. Thus, the conductivity tensor changes point-by-point with fiber shortening, adapting to the fiber paths. An analytical derivation of the conductivity tensor is provided. The volume conductor is then studied with a finite-element approach using the analytically derived conductivity tensor. Representative simulations of single fiber action potentials with the muscle at different degrees of shortening are presented. It is shown that the geometrical changes in the muscle, which imply changes in the conductivity tensor, determine important variations in action potential shape, thus affecting its amplitude and frequency content. The model provides a new tool for interpreting surface EMG signal features with changes in muscle geometry, as it happens during dynamic contractions.