The position of innervation zones in the biceps brachii investigated by surface electromyography

The probability density function of the surface electromyogram and its dependence on contraction force in the vastus lateralis

INTRODUCTION

The probability density function (PDF) of the surface electromyogram (sEMG) depends on contraction force. This dependence, however, has so far been investigated by having the subject generate force at a few fixed percentages of MVC. Here, we examined how the shape of the sEMG PDF changes with contraction force when this force was gradually increased from zero.

METHODS

Voluntary surface EMG signals were recorded from the vastus lateralis of healthy subjects as force was increased in a continuous manner vs. in a step-wise fashion. The sEMG filling process was examined by measuring the EMG filling factor, computed from the non-central moments of the rectified sEMG signal.

RESULTS

(1) In 84% of the subjects, as contraction force increased from 0 to 10% MVC, the sEMG PDF shape oscillated back and forth between the semi-degenerate and the Gaussian distribution. (2) The PDF-force relation varied greatly among subjects for forces between 0 and ~ 10% MVC, but this variability was largely reduced for forces above 10% MVC. (3) The pooled analysis showed that, as contraction force gradually increased, the sEMG PDF evolved rapidly from the semi-degenerate towards the Laplacian distribution from 0 to 5% MVC, and then more slowly from the Laplacian towards the Gaussian distribution for higher forces.

CONCLUSIONS

The study demonstrated that the dependence of the sEMG PDF shape on contraction force can only be reliably assessed by gradually increasing force from zero, and not by performing a few constant-force contractions. The study also showed that the PDF-force relation differed greatly among individuals for contraction forces below 10% MVC, but this variability was largely reduced when force increased above 10% MVC.

The contribution of the tendon electrode to M-wave characteristics in the biceps brachii, vastus lateralis and tibialis anterior

In some compound muscle action potentials (M waves) recorded using the belly-tendon configuration, the tendon electrode makes a noticeable contribution to the M wave. However, this finding has only been demonstrated in some hand and foot muscles. Here, we assessed the contribution of the tendon potential to the amplitude of the vastus lateralis, biceps brachii and tibialis anterior M waves, and we also examined the role of this tendon potential in the shoulder-like feature appearing in most M waves. M waves were recorded separately at the belly and tendon locations of the vastus lateralis, biceps brachii and tibialis anterior from 38 participants by placing the reference electrode at a distant (contralateral) site. The amplitude of the M waves and the latency of their peaks and shoulders were measured. In the vastus lateralis, the tendon potential was markedly smaller in amplitude (∼75%) compared to the belly M wave (P = 0.001), whereas for the biceps brachii and tibialis anterior, the tendon and belly potentials had comparable amplitudes. In the vastus lateralis, the tendon potential showed a small positive peak coinciding in latency with the shoulder of the belly-tendon M wave, whilst in the biceps brachii and tibialis anterior, the tendon potential showed a clear negative peak which coincided in latency with the shoulder. The tendon potential makes a significant contribution to the belly-tendon M waves of the biceps brachii and tibialis anterior muscles, but little contribution to the vastus lateralis M waves. The shoulder observed in the belly-tendon M wave of the vastus lateralis is caused by the belly potential, the shoulder in the biceps brachii M wave is generated by the tendon potential, whereas the shoulder in the tibialis anterior M wave is caused by both the tendon and belly potentials. NEW FINDINGS: What is the central question of this study? Does a tendon electrode make a noticeable contribution to the belly-tendon M wave in the vastus lateralis, biceps brachii and tibialis anterior muscles? What is the main finding and its importance? Because the patellar tendon potential is small in amplitude, it hardly influences the amplitude and shape of the belly-tendon M wave of the vastus lateralis. However, for the biceps brachii and tibialis anterior muscles, the potentials at the tendon sites show a large amplitude, and thus have a great impact on the corresponding belly-tendon M waves.

M-wave changes caused by brief voluntary and stimulated isometric contractions

INTRODUCTION

Under isometric conditions, the increase in muscle force is accompanied by a reduction in the fibers' length. The effects of muscle shortening on the compound muscle action potential (M wave) have so far been investigated only by computer simulation. This study was undertaken to assess experimentally the M-wave changes caused by brief voluntary and stimulated isometric contractions.

METHODS

Two different methods of inducing muscle shortening under isometric condition were adopted: (1) applying a brief (1 s) tetanic contraction and (2) performing brief voluntary contractions of different intensities. In both methods, supramaximal stimulation was applied to the brachial plexus and femoral nerves to evoke M waves. In the first method, electrical stimulation (20 Hz) was delivered with the muscle at rest, whereas in the second, stimulation was applied while participants performed 5-s stepwise isometric contractions at 10, 20, 30, 40, 50, 60, 70, and 100% MVC. The amplitude and duration of the first and second M-wave phases were computed.

RESULTS

The main findings were: (1) on application of tetanic stimulation, the amplitude of the M-wave first phase decreased (~ 10%, P < 0.05), that of the second phase increased (~ 50%, P < 0.05), and the M-wave duration decreased (~ 20%, P < 0.05) across the first five M waves of the tetanic train and then plateaued for the subsequent responses; (2) when superimposing a single electrical stimulus on muscle contractions of increasing forces, the amplitude of the M-wave first phase decreased (~ 20%, P < 0.05), that of the second phase increased (~ 30%, P < 0.05), and M-wave duration decreased (~ 30%, P < 0.05) as force was raised from 0 to 60-70% MVC force.

CONCLUSIONS

The present results will help to identify the adjustments in the M-wave profile caused by muscle shortening and also contribute to differentiate these adjustments from those caused by muscle fatigue and/or changes in Na+-K+ pump activity.

Distribution of innervation zone and muscle fiber conduction velocity in the biceps brachii muscle

The spatial distributions of muscle innervation zone (IZ) and muscle fiber conduction velocity (CV) were examined in nine healthy young male participants. High-density surface electromyography (EMG) was collected from the biceps brachii muscle when subjects performed isometric elbow flexions at 20% to 80% of the maximal voluntary contraction (MVC). A total of 9498 samples of IZs were identified and CVs were calculated using the Radon transform. The center and width of IZ sample distribution were compared within four different force levels and six medial to lateral electrode column positions using repeated measures ANOVA and multiple comparison tests. Significant shifts of IZ center were observed in the medial columns (Columns 5, 6, and 7) compared with the lateral columns (Columns 3 and 4) (p < 0.05). Similarly, significant differences in the IZ width were found in Column 7 and 8 compared to Column 3 (p < 0.05). In contrast, muscle CV was unaffected by column position. Instead, muscle CV was faster at 40% and 80% MVC compared to 20% MVC (p < 0.05). The findings of this study add further insights into the physiological properties of the biceps brachii muscle.

Sarcolemmal Excitability, M-Wave Changes, and Conduction Velocity During a Sustained Low-Force Contraction

This study was undertaken to investigate whether sarcolemmal excitability is impaired during a sustained low-force contraction [10% maximal voluntary contraction (MVC)] by assessing muscle conduction velocity and also by analyzing separately the first and second phases of the muscle compound action potential (M wave). Twenty-one participants sustained an isometric knee extension of 10% MVC for 3min. M waves were evoked by supramaximal single shocks to the femoral nerve given at 10-s intervals. The amplitude, duration, and area of the first and second M-wave phases were computed. Muscle fiber conduction velocity, voluntary surface electromyographic (EMG), perceived effort, MVC force, peak twitch force, and temperature were also recorded. The main findings were: (1) During the sustained contraction, conduction velocity remained unchanged. (2) The amplitude of the M-wave first phase decreased for the first ~30s (−7%, p<0.05) and stabilized thereafter, whereas the second phase amplitude increased for the initial ~30s (+7%, p<0.05), before stabilizing. (3) Both duration and area decreased steeply during the first ~30s, and then more gradually for the rest of the contraction. (4) During the sustained contraction, perceived effort increased fivefold, whereas knee extension EMG increased by ~10%. (5) Maximal voluntary force and peak twitch force decreased (respectively, −9% and −10%, p<0.05) after the low-force contraction. Collectively, the present results indicate that sarcolemmal excitability is well preserved during a sustained 10% MVC task. A depression of the M-wave first phase during a low-force contraction can occur even in the absence of changes in membrane excitability. The development of fatigue during a low-force contraction can occur without alteration of membrane excitability.

Neurophysiological Factors Affecting Muscle Innervation Zone Estimation Using Surface EMG: A Simulation Study

Surface electromyography (EMG) recorded by a linear or 2-dimensional electrode array can be used to estimate the location of muscle innervation zones (IZ). There are various neurophysiological factors that may influence surface EMG and thus potentially compromise muscle IZ estimation. The objective of this study was to evaluate how surface-EMG-based IZ estimation might be affected by different factors, including varying degrees of motor unit (MU) synchronization in the case of single or double IZs. The study was performed by implementing a model simulating surface EMG activity. Three different MU synchronization conditions were simulated, namely no synchronization, medium level synchronization, and complete synchronization analog to M wave. Surface EMG signals recorded by a 2-dimensional electrode array were simulated from a muscle with single and double IZs, respectively. For each situation, the IZ was estimated from surface EMG and compared with the one used in the model for performance evaluation. For the muscle with only one IZ, the estimated IZ location from surface EMG was consistent with the one used in the model for all the three MU synchronization conditions. For the muscle with double IZs, at least one IZ was appropriately estimated from interference surface EMG when there was no MU synchronization. However, the estimated IZ was different from either of the two IZ locations used in the model for the other two MU synchronization conditions. For muscles with a single IZ, MU synchronization has little effect on IZ estimation from electrode array surface EMG. However, caution is required for multiple IZ muscles since MU synchronization might lead to false IZ estimation.

Assessing redistribution of muscle innervation zones after spinal cord injuries

This study aimed to examine the redistribution of neuromuscular junctions or innervation zones (IZs) after spinal cord injuries (SCI). Fifteen able-bodied subjects and 15 subjects with SCI (American Spinal Injury Association Impairment Scale A to D), participated in the study. Surface electromyography (EMG) signals were collected from the biceps brachii muscle by a customized linear electrode array when subjects generated maximal isometric voluntary contractions. The Radon transform was applied to detect the IZ locations in the multiple channel surface EMG signals which were differentiated between consecutive channels. The distribution of IZs was compared between the SCI and control groups using the student-t test. Statistical analysis disclosed a significantly wider range of IZs in the SCI group compared with the control group (SCI: 3.83 ± 1.32 IED, control: 2.83 ± 0. 87 IED, IED: inter-electrode distance, p < 0.05). No remarkable shifts of the center of the distribution were observed between the two groups (SCI: 9.23 ± 2.35 IED, control: 8.53 ± 2.33 IED, p = 0.42). Changes of IZ distribution in the paralyzed muscles could be associated with the complex neuromuscular reorganization after the SCI.

Muscle fibre conduction velocity varies in opposite directions after short- vs. long-duration muscle contractions

INTRODUCTION

The effects of muscle contractions on muscle fibre conduction velocity have normally been investigated for contractions of a given duration and intensity, with most studies being focused on the decline on conduction velocity during/after prolonged contractions. Herein, we perform a systematic analysis of the changes in conduction velocity after voluntary contractions of different durations and intensities.

METHODS

Conduction velocity was estimated in the vastus lateralis before and after knee extensor isometric maximal voluntary contractions (MVCs) of 1, 3, 6, 10, 30 and 60 s, and after brief (3 s) contractions at 10, 30, 50, 70, and 90% of MVC force. Measurements were made during the 10-min period following each contraction.

RESULTS

(1) Conduction velocity was increased immediately after (1 s) the MVCs of brief (≤ 10 s) duration (12 ± 2%, P < 0.05), and then returned rapidly (within 15 s) to control levels; (2) the extent of the increase in conduction velocity was similar after the 3-s, 6-s, and 10-s MVCs (P > 0.05); (3) the magnitude of the increase in conduction velocity after a brief contraction augmented with the intensity of the contraction (increases of 4.6, 7.7, 11.4, 14.8, and 15.2% for contractions at 10, 30, 50, 70, and 90% of MVC force, respectively); (4) conduction velocity was not decreased immediately after the 30-s MVC (P > 0.05); and (5) conduction velocity did not reach its minimum 1 s after the long (≥ 30 s) MVCs.

CONCLUSIONS

Brief (≤ 10 s) muscle contractions induce a short-term increase in conduction velocity, lasting 15 s, while long (≥ 30 s) contractions produce a long-term decrease in conduction velocity, lasting more than 2 min.

The Effect of Spinal Manipulation on the Electrophysiological and Metabolic Properties of the Tibialis Anterior Muscle

There is growing evidence showing that spinal manipulation increases muscle strength in healthy individuals as well as in people with some musculoskeletal and neurological disorders. However, the underlying mechanism by which spinal manipulation changes muscle strength is less clear. This study aimed to assess the effects of a single spinal manipulation session on the electrophysiological and metabolic properties of the tibialis anterior (TA) muscle. Maximum voluntary contractions (MVC) of the ankle dorsiflexors, high-density electromyography (HDsEMG), intramuscular EMG, and near-infrared spectroscopy (NIRS) were recorded from the TA muscle in 25 participants with low level recurring spinal dysfunction using a randomized controlled crossover design. The following outcomes: motor unit discharge rate (MUDR), strength (force at MVC), muscle conduction velocity (CV), relative changes in oxy- and deoxyhemoglobin were assessed pre and post a spinal manipulation intervention and passive movement control. Repeated measures ANOVA was used to assess within and between-group differences. Following the spinal manipulation intervention, there was a significant increase in MVC (p = 0.02; avg 18.87 ± 28.35%) and a significant increase in CV in both the isometric steady-state (10% of MVC) contractions (p < 0.01; avg 22.11 ± 11.69%) and during the isometric ramp (10% of MVC) contractions (p < 0.01; avg 4.52 ± 4.58%) compared to the control intervention. There were no other significant findings. The observed TA strength and CV increase, without changes in MUDR, suggests that the strength changes observed following spinal manipulation are, in part, due to increased recruitment of larger, higher threshold motor units. Further research needs to investigate the longer term and potential functional effects of spinal manipulation in various patients who may benefit from improved muscle function and greater motor unit recruitment.

Innervation zone distribution of the biceps brachii muscle examined using voluntary and electrically-evoked high-density surface EMG

BACKGROUND

High density surface electromyography (EMG) can be used to estimate muscle innervation zones (IZ). The objective of this study was to compare the differences in the distribution of the biceps brachii (BB) IZ derived from voluntary contractions (VC) and electrical stimulation (ES) of the musculocutaneous nerve.

METHODS

Surface EMG signals were recorded from the medial and lateral BB with two 64-channel high density electrode matrices in eight healthy men. The surface EMG was recorded at different percentages of the maximal voluntary contraction (MVC) force (20-100% MVC) and at different percentages of the current needed to elicit a maximal M-wave (20-100% Imax). The IZs of the medial and lateral BB were identified from the EMG signals and expressed as a row number within a given medial-lateral column.

RESULTS

ES current intensity had no significant effect on the group mean IZ location (p > 0.05). However, The IZ during VC was located more proximally with increasing force (p < 0.05), likely due to muscle shortening. The position of the IZ varied slightly (by up to ~ 8 mm) in a medial-lateral direction under both contraction types, but this spatial effect was not significant. The IZ during ES and weak VC (20, 40% MVC) was similar (p > 0.05), but was more proximal in the latter than the former during 60-100% MVC (p < 0.05).

CONCLUSION

ES can be used to detect spatial differences in IZ location free of the confounding effects of muscle shortening and recruitment order of different sized motor units. The method may prove beneficial for locating the IZ in patients who lack voluntary control of their musculature.

Global Innervation Zone Identification With High-Density Surface Electromyography

OBJECTIVE

The aim of this study is to compare the performance of three strategies in determining the global innervation zone (IZ) distribution.

METHODS

High-density surface electromyography was recorded from the biceps brachii muscle of seven healthy subjects under isometric voluntary contractions at 20%, 50%, and 100% of the maximal voluntary contraction and supramaximal musculocutaneous nerve stimulations. IZs were detected: first, by visual identification in a column-specific manner (IZ-1D); second, based on decomposed bipolar mapping of motor unit action potentials (IZ-2D); and third, by source imaging in the three-dimensional muscle space (IZ-3D).

RESULTS

All three IZ detection approaches have exhibited excellent trial-to-trial repeatability. Consistent IZ results were found in the axial direction of the arm across all three approaches, yet a difference was observed in the mediolateral direction.

CONCLUSIONS

Among all three approaches, IZ-3D is capable of providing the most comprehensive information regarding the global IZ distribution, while maintaining high consistency with IZ-1D and IZ-2D results.

SIGNIFICANCE

IZ-3D approach can be a potential tool for global IZ imaging, which is critical to the clinical diagnosis and treatment of neuromuscular disorders.

Motor unit innervation zone localization based on robust linear regression analysis

With the aim of developing a flexible and reliable procedure for superficial muscle innervation zone (IZ) localization, we proposed a method to estimate IZ location using surface electromyogram (EMG) based on robust linear regression. Regression lines were used to model the bidirectional propagation pattern of a single motor unit action potential (MUAP) and visualize the trajectory of the MUAP propagation. IZ localization was performed by identifying the origin of the bidirectional MUAP propagation. Robust linear regression and MUAP peak detection, combined with propagation phase reversal identification, may provide an efficient way to estimate IZ location. Our method offers high resolution in locating IZs based on simulation studies and experimental tests. Furthermore, our method is flexible and may also be applied using a relatively small number of EMG channels. A comparative study of the proposed method with the cross-correlation method for IZ localization was conducted. The results obtained with simulated MUAPs and measured spontaneous MUAPs in the biceps brachii muscle in six subjects (four males and two females, 57 ± 10 years old) with amyotrophic lateral sclerosis (ALS). Our method achieved estimation performance comparable to that obtained by using the cross-correlation method but with higher resolution. This study provides an accurate and practical method to estimate IZ location.

Sarcolemmal membrane excitability during repeated intermittent maximal voluntary contractions

NEW FINDINGS

What is the central question of this study? Is impaired membrane excitability reflected by an increase or by a decrease in M-wave amplitude? What is the main finding and its importance? The magnitude of the M-wave first and second phases changed in completely different ways during intermittent maximal voluntary contractions, leading to the counterintuitive conclusion that it is an increase (and not a decrease) of the M-wave first phase that reflects impaired membrane excitability.

ABSTRACT

The study was undertaken to investigate separately the changes in the first and second phases of the muscle compound action potential (M-wave) during 4 min of intermittent maximal voluntary contractions (MVCs) of the quadriceps. M-waves were evoked by supramaximal single electrical stimulation to the femoral nerve delivered in the resting periods between 48 successive MVCs of 3 s. The amplitude, duration and area of the M-wave first and second phases were measured separately, together with muscle conduction velocity and MVC force. During the intermittent MVCs, the amplitude of the M-wave first phase increased uninterruptedly for the first 3 min (12-16%, P < 0.05) and stabilized thereafter, whereas the second phase initially increased for 55-75 s (11-22%, P < 0.05), but decreased subsequently. The enlargement of the first phase occurred in parallel with an increase in its duration, and concomitantly with a decline in conduction velocity (maximal cross-correlations, 0.89-0.97; time lag, 0 s). Also, a significant temporal association was found between the amplitude of the first phase and MVC force (time lag, 0 s; maximal cross-correlations, 0.85-0.97). Conversely, there was no temporal association between the second phase amplitude and conduction velocity or MVC force (time lag, 73-117 s; maximal cross-correlations, 0.65-0.77). It is concluded that the enlargement of the M-wave first phase is the electrical manifestation of impaired muscle membrane excitability. The results highlight the importance of independently analysing the first and second phases, as only the first phase can be used reliably to detect changes in membrane excitability, while the second might be affected by muscle architecture.

Muscle fiber conduction velocity of the vastus medilais and lateralis muscle after eccentric exercise induced-muscle damage

Change in muscle fiber conduction velocity (MFCV) has been reported after eccentric exercise induces muscle fiber damage, most likely due to a change in membrane permeability of the injured fiber. The extent of damage to the muscle fiber depends on the morphological and architectural characteristics of the muscle fibers. Morphological and architectural characteristics of the VMO muscle fibers are different from VL muscle. Thus, it is expected that eccentric exercise of quadriceps muscle results in a non-uniform fiber damage within the VMO and VL muscle and, as a consequence, non-uniform changes in membrane excitability and conduction velocity. The aim of the study was to investigate MFCV of the VMO and VL muscles before and 24 h after eccentric exercise. Multichannel surface EMG signals were concurrently recorded from the right VMO and VL muscles of 15 healthy men during sustained isometric contractions at 50% of the maximal force. Maximal voluntary force significantly reduced after eccentric exercise with respect to the pre-exercise condition (P < 0.0001). MFCV decreased over time during the sustained contractions at faster rates when assessed 24 h after exercise (VMO = -26.1; VL = -20.1) with respect to the pre-exercise condition (VMO = -9.1; VL = -13.7, P < 0.0001). Moreover, VMO showed a greater rate of reduction in MFCV over sustained contraction (26.1 ± 10.7%) in comparison with VL muscle (20.1 ± 8.5%, P < 0.025) 24 h after eccentric exercise. The result indicates that eccentric exercise contributes to a larger reduction in MFCV within the VMO muscle as compared to the VL muscle. This may abolish the ability of VMO to counteract the lateral pull of the VL muscle during knee extension, thereby leaving the knee complex more vulnerable to injury.

End-of-Fiber Signals Strongly Influence the First and Second Phases of the M Wave in the Vastus Lateralis: Implications for the Study of Muscle Excitability

It has been recurrently observed that, for compound muscle action potentials (M wave) recorded over the innervation zone of the vastus lateralis, the descending portion of the first phase generally shows an “inflection” or “shoulder.” We sought to clarify the electrical origin of this shoulder-like feature and examine its implications. M waves evoked by maximal single shocks to the femoral nerve were recorded in monopolar and bipolar configurations from 126 individuals using classical (10-mm recording diameter, 20-mm inter-electrode distance) electrodes and from eight individuals using small electrodes arranged in a linear array. The changes of the M-wave waveform at different positions along the muscle fibers' direction were examined. The shoulder was identified more frequently in monopolar (97%) than in bipolar (46%) M waves. The shoulder of M waves recorded at different distances from the innervation zone had the same latency. Furthermore, the shoulder of the M wave recorded over the innervation zone coincided in latency with the positive peak of that recorded beyond the muscle. The positive phase of the M wave detected 20 mm away from the innervation zone was essentially composed of non-propagating components. The shoulder-like feature in monopolar and bipolar M waves results from the termination of action potentials at the superficial aponeurosis of the vastus lateralis. We conclude that, only the amplitude of the first phase, and not the second, of M waves recorded monopolarly and/or bipolarly in close proximity to the innervation zone can be used reliably to monitor possible changes in muscle membrane excitability.

Speedup computation of HD-sEMG signals using a motor unit-specific electrical source model

Nowadays, bio-reliable modeling of muscle contraction is becoming more accurate and complex. This increasing complexity induces a significant increase in computation time which prevents the possibility of using this model in certain applications and studies. Accordingly, the aim of this work is to significantly reduce the computation time of high-density surface electromyogram (HD-sEMG) generation. This will be done through a new model of motor unit (MU)-specific electrical source based on the fibers composing the MU. In order to assess the efficiency of this approach, we computed the normalized root mean square error (NRMSE) between several simulations on single generated MU action potential (MUAP) using the usual fiber electrical sources and the MU-specific electrical source. This NRMSE was computed for five different simulation sets wherein hundreds of MUAPs are generated and summed into HD-sEMG signals. The obtained results display less than 2% error on the generated signals compared to the same signals generated with fiber electrical sources. Moreover, the computation time of the HD-sEMG signal generation model is reduced to about 90% compared to the fiber electrical source model. Using this model with MU electrical sources, we can simulate HD-sEMG signals of a physiological muscle (hundreds of MU) in less than an hour on a classical workstation. Graphical Abstract Overview of the simulation of HD-sEMG signals using the fiber scale and the MU scale. Upscaling the electrical source to the MU scale reduces the computation time by 90% inducing only small deviation of the same simulated HD-sEMG signals.

Tracking motor units longitudinally across experimental sessions with high-density surface electromyography

KEY POINTS

Classic motor unit (MU) recording and analysis methods do not allow the same MUs to be tracked across different experimental sessions, and therefore, there is limited experimental evidence on the adjustments in MU properties following training or during the progression of neuromuscular disorders. We propose a new processing method to track the same MUs across experimental sessions (separated by weeks) by using high-density surface electromyography. The application of the proposed method in two experiments showed that individual MUs can be identified reliably in measurements separated by weeks and that changes in properties of the tracked MUs across experimental sessions can be identified with high sensitivity. These results indicate that the behaviour and properties of the same MUs can be monitored across multiple testing sessions. The proposed method opens new possibilities in the understanding of adjustments in motor unit properties due to training interventions or the progression of pathologies.

ABSTRACT

A new method is proposed for tracking individual motor units (MUs) across multiple experimental sessions on different days. The technique is based on a novel decomposition approach for high-density surface electromyography and was tested with two experimental studies for reliability and sensitivity. Experiment I (reliability): ten participants performed isometric knee extensions at 10, 30, 50 and 70% of their maximum voluntary contraction (MVC) force in three sessions, each separated by 1 week. Experiment II (sensitivity): seven participants performed 2 weeks of endurance training (cycling) and were tested pre-post intervention during isometric knee extensions at 10 and 30% MVC. The reliability (Experiment I) and sensitivity (Experiment II) of the measured MU properties were compared for the MUs tracked across sessions, with respect to all MUs identified in each session. In Experiment I, on average 38.3% and 40.1% of the identified MUs could be tracked across two sessions (1 and 2 weeks apart), for the vastus medialis and vastus lateralis, respectively. Moreover, the properties of the tracked MUs were more reliable across sessions than those of the full set of identified MUs (intra-class correlation coefficients ranged between 0.63-0.99 and 0.39-0.95, respectively). In Experiment II, ∼40% of the MUs could be tracked before and after the training intervention and training-induced changes in MU conduction velocity had an effect size of 2.1 (tracked MUs) and 1.5 (group of all identified motor units). These results show the possibility of monitoring MU properties longitudinally to document the effect of interventions or the progression of neuromuscular disorders.

A more precise, repeatable and diagnostic alternative to surface electromyography - an appraisal of the clinical utility of acoustic myography

Acoustic myography (AMG) enables a detailed and accurate measurement of those muscles involved in a particular movement and is independent of electrical signals between the nerve and muscle, measuring solely muscle contractions, unlike surface electromyography (sEMG). With modern amplifiers and digital sound recording systems, measurements during physical activity both inside and outside a laboratory setting are now possible and accurate. Muscle sound gives a representation of the work of each muscle group during a complex movement, and under certain forms of movement even reveals both concentric and eccentric activity, something that sEMG is incapable of. Recent findings suggest that AMG has a number of advantages over sEMG, being simple to use, accurate and repeatable as well as being intuitive to interpret. The AMG signal comprises three physiological parameters, namely efficiency/coordination (E-score), spatial summation (S-score) and temporal summation (T-score). It is concluded that modern AMG units have the potential to accurately assess patients with neuromuscular and musculoskeletal complaints in hospital clinics, home monitoring situations as well as sports settings.

Analysis of the sEMG/force relationship using HD-sEMG technique and data fusion: A simulation study

The relationship between the surface Electromyogram (sEMG) signal and the force of an individual muscle is still ambiguous due to the complexity of experimental evaluation. However, understanding this relationship should be useful for the assessment of neuromuscular system in healthy and pathological contexts. In this study, we present a global investigation of the factors governing the shape of this relationship. Accordingly, we conducted a focused sensitivity analysis of the sEMG/force relationship form with respect to neural, functional and physiological parameters variation. For this purpose, we used a fast generation cylindrical model for the simulation of an 8×8 High Density-sEMG (HD-sEMG) grid and a twitch based force model for the muscle force generation. The HD-sEMG signals as well as the corresponding force signals were simulated in isometric non-fatiguing conditions and were based on the Biceps Brachii (BB) muscle properties. A total of 10 isometric constant contractions of 5s were simulated for each configuration of parameters. The Root Mean Squared (RMS) value was computed in order to quantify the sEMG amplitude. Then, an image segmentation method was used for data fusion of the 8×8 RMS maps. In addition, a comparative study between recent modeling propositions and the model proposed in this study is presented. The evaluation was made by computing the Normalized Root Mean Squared Error (NRMSE) of their fitting to the simulated relationship functions. Our results indicated that the relationship between the RMS (mV) and muscle force (N) can be modeled using a 3rd degree polynomial equation. Moreover, it appears that the obtained coefficients are patient-specific and dependent on physiological, anatomical and neural parameters.

Detection of Multiple Innervation Zones from Multi-Channel Surface EMG Recordings with Low Signal-to-Noise Ratio Using Graph-Cut Segmentation

Knowledge of the location of muscle Innervation Zones (IZs) is important in many applications, e.g. for minimizing the quantity of injected botulinum toxin for the treatment of spasticity or for deciding on the type of episiotomy during child delivery. Surface EMG (sEMG) can be noninvasively recorded to assess physiological and morphological characteristics of contracting muscles. However, it is not often possible to record signals of high quality. Moreover, muscles could have multiple IZs, which should all be identified. We designed a fully-automatic algorithm based on the enhanced image Graph-Cut segmentation and morphological image processing methods to identify up to five IZs in 60-ms intervals of very-low to moderate quality sEMG signal detected with multi-channel electrodes (20 bipolar channels with Inter Electrode Distance (IED) of 5 mm). An anisotropic multilayered cylinder model was used to simulate 750 sEMG signals with signal-to-noise ratio ranging from -5 to 15 dB (using Gaussian noise) and in each 60-ms signal frame, 1 to 5 IZs were included. The micro- and macro- averaged performance indices were then reported for the proposed IZ detection algorithm. In the micro-averaging procedure, the number of True Positives, False Positives and False Negatives in each frame were summed up to generate cumulative measures. In the macro-averaging, on the other hand, precision and recall were calculated for each frame and their averages are used to determine F1-score. Overall, the micro (macro)-averaged sensitivity, precision and F1-score of the algorithm for IZ channel identification were 82.7% (87.5%), 92.9% (94.0%) and 87.5% (90.6%), respectively. For the correctly identified IZ locations, the average bias error was of 0.02±0.10 IED ratio. Also, the average absolute conduction velocity estimation error was 0.41±0.40 m/s for such frames. The sensitivity analysis including increasing IED and reducing interpolation coefficient for time samples was performed. Meanwhile, the effect of adding power-line interference and using other image interpolation methods on the deterioration of the performance of the proposed algorithm was investigated. The average running time of the proposed algorithm on each 60-ms sEMG frame was 25.5±8.9 (s) on an Intel dual-core 1.83 GHz CPU with 2 GB of RAM. The proposed algorithm correctly and precisely identified multiple IZs in each signal epoch in a wide range of signal quality and is thus a promising new offline tool for electrophysiological studies.

Guidelines to electrode positioning for human and animal electrical impedance myography research

The positioning of electrodes in electrical impedance myography (EIM) is critical for accurately assessing disease progression and effectiveness of treatment. In human and animal trials for neuromuscular disorders, inconsistent electrode positioning adds errors to the muscle impedance. Despite its importance, how the reproducibility of resistance and reactance, the two parameters that define EIM, are affected by changes in electrode positioning remains unknown. In this paper, we present a novel approach founded on biophysical principles to study the reproducibility of resistance and reactance to electrode misplacements. The analytical framework presented allows the user to quantify a priori the effect on the muscle resistance and reactance using only one parameter: the uncertainty placing the electrodes. We also provide quantitative data on the precision needed to position the electrodes and the minimum muscle length needed to achieve a pre-specified EIM reproducibility. The results reported here are confirmed with finite element model simulations and measurements on five healthy subjects. Ultimately, our data can serve as normative values to enhance the reliability of EIM as a biomarker and facilitate comparability of future human and animal studies.

Flexible Inkjet-Printed Multielectrode Arrays for Neuromuscular Cartography

Flexible Poly(3,4-ethylenedioxythiophene)-poly(styrenesulfonate) (

PEDOT

PSS) conductive-polymer multielectrode arrays (MEAs) are fabricated without etching or aggressive lift-off processes, only by additive solution processes. Inkjet printing technology has several advantages, such as a customized design and a rapid realization time, adaptability to different patients and to different applications. In particular, inkjet printing technology, as additive and "contactless" technology, can be easily inserted into various technological fabrication steps on different substrates at low cost. In vivo electrochemical impedance spectroscopy measurements show the time stability of such MEAs. An equivalent circuit model is established for such flexible cutaneous MEAs. It is shown that the charge transfer resistance remains the same, even two months after fabrication. Surface electromyography and electrocardiography measurements show that the

PEDOT

PSS MEAs record electrophysiological activity signals that are comparable to those obtained with unitary Ag/AgCl commercial electrodes. Additionally, such MEAs offer parallel and simultaneous recordings on multiple locations at high surface density. It also proves its suitability to reconstruct an innervation zone map and opens new perspectives for a better control of amputee's myoelectric prostheses. The employment of additive technologies such as inkjet printing suggests that the integration of multifunctional sensors can improve the performances of ultraflexible brain-computer interfaces.

High-density surface electromyography provides reliable estimates of motor unit behavior

OBJECTIVE

To assess the intra- and inter-session reliability of estimates of motor unit behavior and muscle fiber properties derived from high-density surface electromyography (HDEMG).

METHODS

Ten healthy subjects performed submaximal isometric knee extensions during three recording sessions (separate days) at 10%, 30%, 50% and 70% of their maximum voluntary effort. The discharge timings of motor units of the vastus lateralis and medialis muscles were automatically identified from HDEMG by a decomposition algorithm. We characterized the number of detected motor units, their discharge rates, the coefficient of variation of their inter-spike intervals (CoVisi), the action potential conduction velocity and peak-to-peak amplitude. Reliability was assessed for each motor unit characteristics by intra-class correlation coefficient (ICC). Additionally, a pulse-to-noise ratio (PNR) was calculated, to verify the accuracy of the decomposition.

RESULTS

Good to excellent reliability within and between sessions was found for all motor unit characteristics at all force levels (ICCs>0.8), with the exception of CoVisi that presented poor reliability (ICC<0.6). PNR was high and similar for both muscles with values ranging between 45.1 and 47.6dB (accuracy>95%).

CONCLUSION

Motor unit features can be assessed non-invasively and reliably within and across sessions over a wide range of force levels.

SIGNIFICANCE

These results suggest that it is possible to characterize motor units in longitudinal intervention studies.

Motor Neuron Pools of Synergistic Thigh Muscles Share Most of Their Synaptic Input

Neural control of synergist muscles is not well understood. Presumably, each muscle in a synergistic group receives some unique neural drive and some drive that is also shared in common with other muscles in the group. In this investigation, we sought to characterize the strength, frequency spectrum, and force dependence of the neural drive to the human vastus lateralis and vastus medialis muscles during the production of isometric knee extension forces at 10 and 30% of maximum voluntary effort. High-density surface electromyography recordings were decomposed into motor unit action potentials to examine the neural drive to each muscle. Motor unit coherence analysis was used to characterize the total neural drive to each muscle and the drive shared between muscles. Using a novel approach based on partial coherence analysis, we were also able to study specifically the neural drive unique to each muscle (not shared). The results showed that the majority of neural drive to the vasti muscles was a cross-muscle drive characterized by a force-dependent strength and bandwidth. Muscle-specific neural drive was at low frequencies (<5 Hz) and relatively weak. Frequencies of neural drive associated with afferent feedback (6-12 Hz) and with descending cortical input (∼20 Hz) were almost entirely shared by the two muscles, whereas low-frequency (<5 Hz) drive comprised shared (primary) and muscle-specific (secondary) components. This study is the first to directly investigate the extent of shared versus independent control of synergist muscles at the motor neuron level.

SIGNIFICANCE STATEMENT

Precisely how the nervous system coordinates the activity of synergist muscles is not well understood. One possibility is that muscles of a synergy share a common neural drive. In this study, we directly compared the relative strength of shared versus independent neural drive to synergistically activated thigh muscles in humans. The results of this analysis support the notion that synergistically activated muscles share most of their neural drive. Scientifically, this study addressed an important gap in our current understanding of how neural drive is delivered to synergist muscles. We have also demonstrated the feasibility of a novel approach to the study of muscle synergies based on partial coherence analysis of motor unit activity.

Spatial variability in cortex-muscle coherence investigated with magnetoencephalography and high-density surface electromyography

Cortex-muscle coherence (CMC) reflects coupling between magnetoencephalography (MEG) and surface electromyography (sEMG), being strongest during isometric contraction but absent, for unknown reasons, in some individuals. We used a novel nonmagnetic high-density sEMG (HD-sEMG) electrode grid (36 mm × 12 mm; 60 electrodes separated by 3 mm) to study effects of sEMG recording site, electrode derivation, and rectification on the strength of CMC. Monopolar sEMG from right thenar and 306-channel whole-scalp MEG were recorded from 14 subjects during 4-min isometric thumb abduction. CMC was computed for 60 monopolar, 55 bipolar, and 32 Laplacian HD-sEMG derivations, and two derivations were computed to mimic "macroscopic" monopolar and bipolar sEMG (electrode diameter 9 mm; interelectrode distance 21 mm). With unrectified sEMG, 12 subjects showed statistically significant CMC in 91-95% of the HD-sEMG channels, with maximum coherence at ∼25 Hz. CMC was about a fifth stronger for monopolar than bipolar and Laplacian derivations. Monopolar derivations resulted in most uniform CMC distributions across the thenar and in tightest cortical source clusters in the left rolandic hand area. CMC was 19-27% stronger for HD-sEMG than for "macroscopic" monopolar or bipolar derivations. EMG rectification reduced the CMC peak by a quarter, resulted in a more uniformly distributed CMC across the thenar, and provided more tightly clustered cortical sources than unrectifed sEMGs. Moreover, it revealed CMC at ∼12 Hz. We conclude that HD-sEMG, especially with monopolar derivation, can facilitate detection of CMC and that individual muscle anatomy cannot explain the high interindividual CMC variability.

Prolonged passive static stretching-induced innervation zone shift in biceps brachii

The purpose of this study was to examine the influence of a bout of repeated and prolonged passive static stretching on the innervation zone (IZ) location of the human biceps brachii muscle. Eleven men performed 12 sets of 100-s passive stretches on their biceps brachii. Before (Pre) and immediately after (Post) the stretching intervention, isometric strength was tested during the maximal voluntary contractions (MVCs) of the forearm flexors. The subjects also performed several separate isometric forearm flexion muscle actions at 30%, 50%, and 70% of their predetermined MVCs for examining the locations of the IZ at different contraction intensities. The IZ was identified through multi-channel surface electromyographic (EMG) recordings from a linear electrode array. The stretching intervention induced an average of 10% isometric strength loss for the forearm flexors (mean ± SD: Pre-MVC vs. Post-MVC = 332.12 ± 59.40 N vs. 299.53 ± 70.51 N; p < 0.001). In addition, the average IZ shift was nearly 4.5 mm in average in the proximal direction. However, this shift was not specific to the contraction intensity. We believe that the IZ shift was caused by the elongation of the entire muscle-tendon unit in the proximal direction. Therefore, caution should be taken when using surface EMG technique to examine possible changes in the EMG variables after a stretching protocol, as these variables can be contaminated by the shift of the IZ.

Investigation of innervation zone shift with continuous dynamic muscle contraction

Innervation zone (IZ) has been identified as the origin of action potential propagation in isometric contraction. However, IZ shifts with changes in muscle length during muscle activity. The IZ shift has been estimated using raw EMG signals. This study aimed to investigate the movement of IZ location during continuous dynamic muscle contraction, using a computer program. Subjects flexed their elbow joint as repetitive dynamic muscle contractions. EMG signals were recorded from the biceps brachii muscle using an eight-channel surface electrode array. Approximately 100 peaks from EMG signals were detected for each channel and summed to estimate the IZ location. For each subject, the estimated IZ locations were subtracted from the IZ location during isometric contractions with the elbow flexed at 90°. The results showed that the IZ moved significantly with elbow joint movement from 45° to 135°. However, IZ movement was biased with only a 3.9 mm IZ shift on average when the elbow angle was acute but a 16 mm IZ shift on average when it was obtuse. The movement of IZ location during continuous dynamic muscle contraction can be investigated using this signal processing procedure without subjective judgment.

Spatial distribution of surface action potentials generated by individual motor units in the human biceps brachii muscle

This study analyses the spatial distribution of individual motor unit potentials (MUPs) over the skin surface and the influence of motor unit depth and recording configuration on this distribution. Multichannel surface (13×5 electrode grid) and intramuscular (wire electrodes inserted with needles of lengths 15 and 25mm) electromyographic (EMG) signals were concurrently recorded with monopolar derivations from the biceps brachii muscle of 10 healthy subjects during 60-s isometric contractions at 20% of the maximum torque. Multichannel monopolar MUPs of the target motor unit were obtained by spike-triggered averaging of the surface EMG. Amplitude and frequency characteristics of monopolar and bipolar MUPs were calculated for locations along the fibers' direction (longitudinal), and along the direction perpendicular (transverse) to the fibers. In the longitudinal direction, monopolar and bipolar MUPs exhibited marked amplitude changes that extended for 16-32mm and 16-24mm over the innervation and tendon zones, respectively. The variation of monopolar and bipolar MUP characteristics was not symmetrical about the innervation zone. Motor unit depth had a considerable influence on the relative longitudinal variation of amplitude for monopolar MUPs, but not for bipolar MUPs. The transverse extension of bipolar MUPs ranged between 24 and 32mm, whereas that of monopolar MUPs ranged between 72 and 96mm. The mean power spectral frequency of surface MUPs was highly dependent on the transverse electrode location but not on depth. This study provides a basis for the interpretation of the contribution of individual motor units to the interference surface EMG signal.

Wrist torque estimation during simultaneous and continuously changing movements: surface vs. untargeted intramuscular EMG

In this paper, the predictive capability of surface and untargeted intramuscular electromyography (EMG) was compared with respect to wrist-joint torque to quantify which type of measurement better represents joint torque during multiple degrees-of-freedom (DoF) movements for possible application in prosthetic control. Ten able-bodied subjects participated in the study. Surface and intramuscular EMG was recorded concurrently from the right forearm. The subjects were instructed to track continuous contraction profiles using single and combined DoF in two trials. The association between torque and EMG was assessed using an artificial neural network. Results showed a significant difference between the two types of EMG (P < 0.007) for all performance metrics: coefficient of determination (R(2)), Pearson correlation coefficient (PCC), and root mean square error (RMSE). The performance of surface EMG (R(2) = 0.93 ± 0.03; PCC = 0.98 ± 0.01; RMSE = 8.7 ± 2.1%) was found to be superior compared with intramuscular EMG (R(2) = 0.80 ± 0.07; PCC = 0.93 ± 0.03; RMSE = 14.5 ± 2.9%). The higher values of PCC compared with R(2) indicate that both methods are able to track the torque profile well but have some trouble (particularly intramuscular EMG) in estimating the exact amplitude. The possible cause for the difference, thus the low performance of intramuscular EMG, may be attributed to the very high selectivity of the recordings used in this study.

Quantifying the effects of electrode distance from the innervation zone on the electromyographic amplitude versus torque relationships

The present study applied a log-transformation model to compare the electromyographic (EMG) amplitude versus torque relationships from monopolar EMG signals up to 35 mm proximal and distal from the innervation zone (IZ). Seven men (age = 23 ± 2 year; mass = 82 ± 10 kg) and two women (age = 21 ± 1 year; mass = 62 ± 8 kg) performed isometric ramp contractions of the right leg extensors with an eight-channel linear electrode array positioned over the vastus lateralis with the IZ located between channels 4 and 5. Linear regression models were fit to the log-transformed monopolar EMG(RMS)-torque relationships with the b terms (slope) and the a terms (Y-intercept) calculated for each channel and subject. The b terms for channels 4, 5, and 6 were higher (P ≤ 0.05) than the more distal channels 7 and 8 (P < 0.05). In contrast, there were no differences (P > 0.05) among the a terms of the eight channels. Thus, the shapes of the monopolar EMG(RMS)-torque relationships were altered as a function of distance between the IZ and recording area, which may be helpful for clinicians and researchers who infer changes in motor control strategies based on the shapes of the EMG(RMS)-torque relationships.

Distribution of innervation zones in the human biceps brachii

Two types of multicontact surface electrodes were used to define the distribution of innervation zones (motor end-plate regions) in human skeletal muscles: a linear array electrode and a two-dimensional grid electrode. These electrodes detected motor unit action potentials (MUAPs) propagating along the muscle fibers during a voluntary contraction and marked the source of the propagation, which was assumed to indicate the position of the innervation zones. The linear array electrode defined the absolute position of the innervation zones in the muscle, but it could not clarify the detailed distribution, within the innervation zones. On the other hand, the grid electrode clarified the detailed configuration of innervation zones under the pickup area of the electrode, but the absolute position of the sources was inaccurate. The records obtained with these two electrodes were combined to draw a detailed map of the innervation zones over the entire surface of the muscle. In seven healthy normal adults, the biceps brachii was studied. The innervation zones were positioned at the middle length of the muscle in a band 30-60 mm wide. The detailed configuration of the innervation zones varied between the subjects.

A novel approach to surface electromyography: an exploratory study of electrode-pair selection based on signal characteristics

A 3 × 4 electrode array was placed over each of seven muscles and surface electromyography (sEMG) data were collected during isometric contractions. For each array, nine bipolar electrode pairs were formed off-line and sEMG parameters were calculated and evaluated based on repeatability across trials and comparison to an anatomically placed electrode pair. The use of time-domain parameters for the selection of an electrode pair from within a grid-like array may improve upon existing electrode placement methodologies.

Effect of innervation zones in estimating biceps brachii force-EMG relationship during isometric contraction

Measuring muscle forces in vivo is invasive and consequently indirect methods e.g., electromyography (EMG) are used in estimating muscular force production. The aim of the present paper was to examine what kind of effect the disruption of the physiological signal caused by the innervation zone has in predicting the force/torque output from surface EMG. Twelve men (age 26 (SD ±3)years; height 179 (±6)cm; body mass 73 (±6)kg) volunteered as subjects. They were asked to perform maximal voluntary isometric contraction (MVC) in elbow flexion, and submaximal contractions at 10%, 20%, 30%, 40%, 50% and 75% of the recorded MVC. EMG was measured from biceps brachii muscle with an electrode grid of 5 columns×13 rows. Force-EMG relationships were determined from individual channels and as the global mean value. The relationship was deemed inconsistent if EMG value did not increase in successive force levels. Root mean squared errors were calculated for 3rd order polynomial fits. All subjects had at least one (4-52) inconsistent channel. Two subjects had inconsistent relationship calculated from the global mean. The mean root mean squared error calculated using leave one out method for the fits of the individual channels (0.33±0.17) was higher (P<0.001) than the error for the global mean fit (0.16±0.08). It seems that the disruption of the physiological signal caused by the innervation zone affects the consistency of the force-EMG relationship on single bipolar channel level. Multichannel EMG recordings used for predicting force overcame this disruption.

Mapping the motor point in the human tibialis anterior muscle

OBJECTIVE

Percutaneous electrical stimulation of the motor point permits selective activation of a muscle. However, the definition and number of motor points reported for a given muscle varies. Our goal was to address these problems.

METHODS

The area, location and number of motor points in human tibialis anterior were examined, using isometric dorsiflexion torque responses to electrical stimuli. Three methods were used: lowest electrical threshold, maximum muscle response, and approximate motor point.

RESULTS

A single motor point was identified in 39/40 subjects regardless of method. The area of the site of lowest electrical threshold was smaller (median, 35 mm(2)) than those using the maximum muscle response (144 mm(2)) and approximate motor point (132 mm(2)). There was substantial, but not significant, between-subject variation in motor point location. Fifty three percent of motor points would have been missed if located only by reference to anatomical landmarks.

CONCLUSIONS

These results suggested that the motor point's location cannot be determined a priori and that the identification method will affect both area and location.

SIGNIFICANCE

If it is important to maximally activate a single muscle in isolation, the motor point is best represented by the site producing a maximal but isolated muscle response at the lowest stimulation intensity.

Reliability of surface EMG matrix in locating the innervation zone of upper trapezius muscle

The identification of the motor unit (MU) innervation zone (IZ) using surface electromyographic (sEMG) signals detected on the skin with a linear array or a matrix of electrodes has been recently proposed in the literature. However, an analysis of the reliability of this procedure and, therefore, of the suitability of the sEMG signals for this purpose has not been reported. The purpose of this work is to describe the intra and inter-rater reliability and the suitability of surface EMG in locating the innervation zone of the upper trapezius muscle. Two operators were trained on electrode matrix positioning and sEMG signal analysis. Ten healthy subjects, instructed to perform a series of isometric contractions of the upper trapezius muscle participated in the study. The two operators collected sEMG signals and then independently estimated the IZ location through visual analysis. Results showed an almost perfect agreement for intra-rater and inter-rater reliability. The constancy of IZ location could be affected by the factors reflecting the population of active MUs and their IZs, including: the contraction intensity, the acquisition period analyzed, the contraction repetition. In almost all cases the IZ location shift due to these factors did not exceed 4mm. Results generalization to other muscles should be made with caution.

Muscle fibre conduction velocity during a 30-s Wingate anaerobic test

Ten male volunteers (age 29.2 ± 5.2 years, mean ± SD) were recruited to test the hypothesis that muscle fibre conduction velocity (MFCV) would decrease with power output during a 30-s Wingate test on a mechanically-braked cycle ergometer. Prior to the main test, the optimal pre-fixed load corresponding to the highest power output was selected following a random series of six 10-s sprints. Surface electromyographic (EMG) signals were detected from the right vastus lateralis with linear adhesive arrays of eight electrodes. Power output decreased significantly from 6-s until the end of the test (860.9 ± 207.8 vs. 360.9 ± 11.4 W, respectively) and was correlated with MFCV (R=0.543, P<0.01), which also declined significantly by 26.8 ± 11% (P<0.05). There was a tendency for the mean frequency of the EMG power spectrum (MNF) to decrease, but average rectified values (ARV) remained unchanged throughout the test. The parallel decline of MFCV with power output suggests changes in fibre membrane properties. The unaltered ARV, together with the declined MFCV, would indicate either a decrease in discharge rate, de-recruitment of fatigued motor units or elongation of still present motor unit action potentials.