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

High-Density Electromyography Excitation in Front vs. Back Lat Pull-Down Prime Movers

The current study compared the spatial excitation of the primary muscles during the lat pull-down exercise with the bar passing in front (front-LPD) or behind the neck (back-LPD) using high-density electromyography. Fourteen resistance trained men performed a front-LPD or a back-LPD within a non-fatiguing set with 8-RM as the external load. The muscle excitation centroid of latissimus dorsi, middle trapezius, pectoralis major, biceps brachii, triceps brachii and posterior deltoid muscles were recorded during the ascending and the descending phase. During the descending phase, the front-LPD showed superior excitation of the latissimus dorsi (ES = 0.97) and the pectoralis major (ES = 1.17), while in the ascending phase, the back-LPD exhibited superior excitation of the latissimus dorsi (ES = 0.63), and the front-LPD showed superior excitation of the biceps brachii (ES = 0.41) and the posterior deltoid (ES = 1.77). During the descending phase, the front-LPD showed a more lateral centroid of the latissimus dorsi (ES = 0.60), the biceps brachii (ES = 0.63) and the triceps brachii (ES = 0.98), while the centroid was more medial for the middle trapezius (ES = 0.58). The centroid of the middle trapezius was also more medial in the front-LPD during the ascending phase (ES = 0.85). The pectoralis major centroid was more cranial in the front-LPD for both the descending (ES = 1.58) and the ascending phase (ES = 0.88). The front-LPD appears to provide overall greater excitation in the prime movers. However, distinct spatial excitation patterns were observed, making exercise suitable for the training routine.

Inferring position of motor units from high-density surface EMG

The spatial distribution of muscle fibre activity is of interest in guiding therapy and assessing recovery of motor function following injuries of the peripheral or central nervous system. This paper presents a new method for stable estimation of motor unit territory centres from high-density surface electromyography (HDsEMG). This completely automatic process applies principal component compression and a rotatable Gaussian surface fit to motor unit action potential (MUAP) distributions to map the spatial distribution of motor unit activity. Each estimated position corresponds to the signal centre of the motor unit territory. Two subjects were used to test the method on forearm muscles, using two different approaches. With the first dataset, motor units were identified by decomposition of intramuscular EMG and the centre position of each motor unit territory was estimated from synchronized HDsEMG data. These positions were compared to the positions of the intramuscular fine wire electrodes with depth measured using ultrasound. With the second dataset, decomposition and motor unit localization was done directly on HDsEMG data, during specific muscle contractions. From the first dataset, the mean estimated depth of the motor unit centres were 8.7, 11.6, and 9.1 mm, with standard deviations 0.5, 0.1, and 1.3 mm, and the respective depths of the fine wire electrodes were 8.4, 15.8, and 9.1 mm. The second dataset generated distinct spatial distributions of motor unit activity which were used to identify the regions of different muscles of the forearm, in a 3-dimensional and projected 2-dimensional view. In conclusion, a method is presented which estimates motor unit centre positions from HDsEMG. The study demonstrates the shifting spatial distribution of muscle fibre activity between different efforts, which could be used to assess individual muscles on a motor unit level.

The filling factor of the sEMG signal at low contraction forces in the quadriceps muscles is influenced by the thickness of the subcutaneous layer

Introduction: It has been shown that, for male subjects, the sEMG activity at low contraction forces is normally "pulsatile", i.e., formed by a few large-amplitude MUPs, coming from the most superficial motor units. The subcutaneous layer thickness, known to be greater in females than males, influences the electrode detection volume. Here, we investigated the influence of the subcutaneous layer thickness on the type of sEMG activity (pulsatile vs. continuous) at low contraction forces. Methods: Voluntary surface EMG signals were recorded from the quadriceps muscles of healthy males and females as force was gradually increased from 0% to 40% MVC. 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) The sEMG activity at low contraction forces was "continuous" in the VL, VM and RF of females, whereas this sEMG activity was "pulsatile" in the VL and VM of males. 2) The filling factor at low contraction forces was lower in males than females for the VL (p = 0.003) and VM (p = 0.002), but not for the RF (p = 0.54). 3) The subcutaneous layer was significantly thicker in females than males for the VL (p = 0.001), VM (p = 0.001), and RF (p = 0.003). 4) A significant correlation was found in the vastus muscles between the subcutaneous layer thickness and the filling factor (p < 0.05). Discussion: The present results indicate that the sEMG activity at low contraction forces in the female quadriceps muscles is "continuous" due to the thick subcutaneous layer of these muscles, which impedes an accurate assessment of the sEMG filling process.

Differences in motor unit firing properties of the vastus lateralis muscle during postural and voluntary tasks

The firing properties of the motor units are usually affected by the motor task. However, it has not been clarified whether the firing properties of the motor units of a specific muscle are different between postural and voluntary tasks. Therefore, this study investigated whether the recruitment and rate coding of the motor units differ between these two motor tasks. Thirteen healthy volunteers performed trapezoidal muscle contraction with a target value of 15% maximum electromyography (EMG) activity by voluntary left knee extension in the sitting position (voluntary task) and postural maintenance in the semi-squatting position (postural task) with a knee flexion angle of 30°. We obtained four channels of surface EMG activity during each task from left vastus lateralis muscle. We extracted the firing properties of individual motor units using the EMG decomposition algorithm. The recruitment threshold and motor unit action potential amplitude were significantly lower in the postural task than in the voluntary task, and conversely, the mean firing rate was significantly higher. These results were explained by the preferential recruitment of motor units with higher recruitment threshold and amplitude in the voluntary task, while motor units with lower recruitment threshold and higher firing rate were preferentially recruited in the postural task. Preferential activation of fatigue-resistant motor units in the postural task is a reasonable strategy as it allows for sustained postural maintenance. We provide the first evidence that motor unit firing properties are clearly different between postural and voluntary tasks, even at the same muscle activity level.

Performance Evaluation of a Wearable Tattoo Electrode Suitable for High-Resolution Surface Electromyogram Recording

OBJECTIVE

High-density surface electromyography (HD-sEMG) has been utilized extensively in neuromuscular research. Despite its potential advantages, limitations in electrode design have largely prevented widespread acceptance of the technology. Commercial electrodes have limited spatial fidelity, because of a lack of sharpness of the signal, and variable signal stability. We demonstrate here a novel tattoo electrode that addresses these issues. Our dry HD electrode grid exhibits remarkable deformability which ensures superior conformity with the skin surface, while faithfully recording signals during different levels of muscle contraction.

METHOD

We fabricated a 4 cm×3 cm tattoo HD electrode grid on a stretchable electronics membrane for sEMG applications. The grid was placed on the skin overlying the biceps brachii of healthy subjects, and was used to record signals for several hours while tracking different isometric contractions.

RESULTS

The sEMG signals were recorded successfully from all 64 electrodes across the grid. These electrodes were able to faithfully record sEMG signals during repeated contractions while maintaining a stable baseline at rest. During voluntary contractions, broad EMG frequency content was preserved, with accurate reproduction of the EMG spectrum across the full signal bandwidth.

CONCLUSION

The tattoo grid electrode can potentially be used for recording high-density sEMG from skin overlying major limb muscles. Layout programmability, good signal quality, excellent baseline stability, and easy wearability make this electrode a potentially valuable component of future HD electrode grid applications.

SIGNIFICANCE

The tattoo electrode can facilitate high fidelity recording in clinical applications such as tracking the evolution and time-course of challenging neuromuscular degenerative disorders.

Does the amplitude of biceps brachii M waves increase similarly in both limbs during staircase, electrically elicited contractions?

OBJECTIVE

Humans usually tend to control more finely muscle force production in dominant than non-dominant upper limbs. It is well established that motor unit recruitment is a key mechanism by which muscle force is controlled, and we hypothesized that a relatively smaller number of motor units may be recruited in muscles of dominant than non-dominant limbs for any given increase in synaptic input. Hence, we investigated peripheral properties of dominant and non-dominant biceps brachii through the analysis of M-wave responses to incremental electrical stimulation.

APPROACH

Current pulses at progressively greater intensities were applied in the proximal region of biceps brachii of 16 subjects while surface electromyograms were recorded with a grid of electrodes in the distal region. M-wave amplitude was averaged across channels and normalized with respect to the maximum amplitude value, separately for each stimulation intensity and limb. Amplitude-current intensity curves were interpolated to provide an equal number of stimulation levels between limbs. Differences between dominant and non-dominant arms were assessed through the average increase in M-wave amplitude for consecutive stimulation intensities (increments).

MAIN RESULTS

Wilcoxon's signed-rank test showed that increments in the M-wave amplitude were significantly smaller (p  =  0.017) in dominant than non-dominant biceps brachii.

SIGNIFICANCE

The results suggest that there was a more gradual recruitment of motor units in biceps brachii of dominant than non-dominant arms. This is in agreement with the hypothesis of a broader spectrum of motor unit recruitment thresholds in the dominant arm, which may contribute to a finer regulation of force production.

Identification of regional activation by factorization of high-density surface EMG signals: A comparison of Principal Component Analysis and Non-negative Matrix factorization

In this study, we investigated whether principal component analysis (PCA) and non-negative matrix factorization (NMF) perform similarly for the identification of regional activation within the human vastus medialis. EMG signals from 64 locations over the VM were collected from twelve participants while performing a low-force isometric knee extension. The envelope of the EMG signal of each channel was calculated by low-pass filtering (8 Hz) the monopolar EMG signal after rectification. The data matrix was factorized using PCA and NMF, and up to 5 factors were considered for each algorithm. Association between explained variance, spatial weights and temporal scores between the two algorithms were compared using Pearson correlation. For both PCA and NMF, a single factor explained approximately 70% of the variance of the signal, while two and three factors explained just over 85% or 90%. The variance explained by PCA and NMF was highly comparable (R > 0.99). Spatial weights and temporal scores extracted with non-negative reconstruction of PCA and NMF were highly associated (all p < 0.001, mean R > 0.97). Regional VM activation can be identified using high-density surface EMG and factorization algorithms. Regional activation explains up to 30% of the variance of the signal, as identified through both PCA and NMF.

Surface electromyography can quantify temporal and spatial patterns of activation of intrinsic human foot muscles

Intrinsic foot muscles (IFM) are a crucial component within the human foot. Investigating their functioning can help understand healthy and pathological behaviour of foot and ankle, fundamental for everyday activities. Recording muscle activation from IFM has been attempted with invasive techniques, mainly investigating single muscles. Here we present a novel methodology, to investigate the feasibility of recording physiological surface EMG (sEMG) non-invasively and quantify patterns of activation across the whole plantar region of the foot. sEMG were recorded with a 13 × 5 array from the sole of the foot (n = 25) during two-foot stance, two-foot tiptoe and anterior/posterior sways. Physiological features of sEMG were analysed. During anterior/posterior epochs within the sway task, sEMG patterns were analysed in terms of signal amplitude (intensity) and structure (Sample Entropy) distribution, by evaluating the centre of gravity (CoG) of each topographical map. Results suggest signals are physiological and not affected by loading. Both amplitude and sample entropy CoG coordinates were grouped in one region and overlapped, suggesting that the region with highest amplitude corresponds with the most predictable signal. Therefore, both spatial and temporal features of IFM activation may be recorded non-invasively, providing opportunity for more detailed investigation of IFM function in healthy and patient populations.

Spatial Analysis of Multichannel Surface EMG in Hemiplegic Stroke

We investigated spatial activation patterns of upper extremity muscles during isometric force generation in both intact persons and in hemispheric stroke survivors. We used a 128-channel surface electromyogram (EMG) grid to record the electrical activity of biceps brachii muscles during these contractions. EMG data were processed to develop 2-D root mean square (RMS) maps of muscle activity. Our objective was to determine whether motor impairments following stroke were associated with changes in the muscle activity maps and in the spatial distribution of muscular activation. We found that, for a given subject, spatial patterns in muscle activity maps were consistent across all measured contraction levels differing only the RMS EMG. However, the maps from opposite arms (stroke-affected versus non-affected) of stroke survivors were significantly different from each other, especially when compared with the differences observed intact participants. Our analyses revealed that chronic stroke altered the size and location of the active region in these maps. The former is potentially related to disruption of fiber and tissue structure, possibly linked to factors such as extracellular fat accumulation, connective tissue infiltration, muscle fiber atrophy, fiber shortening, and fiber loss. Changes in spatial patterns in muscle activity maps may also be linked to a shift in the location of the innervation zone or the endplate region of muscles. Furthermore, the textural analysis of EMG activity maps showed a larger pixel-to-pixel variability in stroke-affected muscles. Alterations in the muscle activity maps were also related to functional impairment (estimated using Fugl-Meyer score) and to the degree of spasticity (estimated using the modified Ashworth scale). Overall, our investigation revealed that the muscle architecture and morphology were significantly altered in the chronic stroke.

Correlation between discharge timings of pairs of motor units reveals the presence but not the proportion of common synaptic input to motor neurons

We investigated whether correlation measures derived from pairs of motor unit (MU) spike trains are reliable indicators of the degree of common synaptic input to motor neurons. Several 50-s isometric contractions of the biceps brachii muscle were performed at different target forces ranging from 10 to 30% of the maximal voluntary contraction relying on force feedback. Forty-eight pairs of MUs were examined at various force levels. Motor unit synchrony was assessed by cross-correlation analysis using three indexes: the output correlation as the peak of the cross-histogram (ρ) and the number of synchronous spikes per second (CIS) and per trigger (E). Individual analysis of MU pairs revealed that ρ, CIS, and E were most often positively associated with discharge rate (87, 85, and 76% of the MU pairs, respectively) and negatively with interspike interval variability (69, 65, and 62% of the MU pairs, respectively). Moreover, the behavior of synchronization indexes with discharge rate (and interspike interval variability) varied greatly among the MU pairs. These results were consistent with theoretical predictions, which showed that the output correlation between pairs of spike trains depends on the statistics of the input current and motor neuron intrinsic properties that differ for different motor neuron pairs. In conclusion, the synchronization between MU firing trains is necessarily caused by the (functional) common input to motor neurons, but it is not possible to infer the degree of shared common input to a pair of motor neurons on the basis of correlation measures of their output spike trains.NEW & NOTEWORTHY The strength of correlation between output spike trains is only poorly associated with the degree of common input to the population of motor neurons. The synchronization between motor unit firing trains is necessarily caused by the (functional) common input to motor neurons, but it is not possible to infer the degree of shared common input to a pair of motor neurons on the basis of correlation measures of their output spike trains.

Design of new multi-channel electrodes for surface electromyography signals for signal-processing

This paper covers the design aspects of a new multi-channel electrode for the acquisition of surface electromyography signals from a selected muscle. The new multi-channel electrode has 11 pins where the monopolar signals produced will be configured in a software either as Linear array or Laplacian configuration. The design specification of the pre-amplifier ideally was to have a voltage gain of 500 with bandpass filtering of 5 Hz-1 kHz. The final design of the pre-amplifier circuit using an INA 118 instrumentation amplifier was built and tested to give values for voltage gain of 484 with bandpass filtering of 6.8 Hz-1.02 kHz. The software configuration that gives clearer and more defined signals in terms of motor unit action potentials for future signal processing is the Laplacian rather than Linear array.

Regional activation within the vastus medialis in stimulated and voluntary contractions

This study examined the contribution of muscle fiber orientation at different knee angles to regional activation identified with high-density surface electromyography (HDsEMG). Monopolar HDsEMG signals were collected using a grid of 13 × 5 electrodes placed over the vastus medialis (VM). Intramuscular electrical stimulation was used to selectively activate two regions within VM. The distribution of EMG responses to stimulation was obtained by calculating the amplitude of the compound action potential for each channel; the position of the peak amplitude was tracked across knee angles to describe shifts of the active muscle regions under the electrodes. In a separate experiment, regional activation was investigated in 10 knee flexion-extension movements against a fixed resistance. Intramuscular stimulation of different VM regions resulted in clear differences in amplitude distribution along the columns of the electrode grid (P < 0.001); changes in knee angle resulted in consistent shifts along the rows (P < 0.01) and negligible shifts along the columns of the electrode grid. Regional VM activation was identified in dynamic movement, with distal shifts of the EMG distribution in the eccentric phase of the movement (P < 0.05) and at more flexed knee angles (P < 0.05). HDsEMG was used to describe regional activation across the VM that was not attributable to anatomic factors. Changes in muscle fiber orientation associated with knee joint angle mainly influence the amplitude distribution along the fiber direction. Future studies are needed to understand possible functional roles for regional activation within the VM in dynamic tasks.

Motor units in the human medial gastrocnemius muscle are not spatially localized or functionally grouped

KEY POINTS

Human medial gastrocnemius (MG) motor units (MUs) are thought to occupy small muscle territories or regions, with low-threshold units preferentially located distally. We used intramuscular recordings to measure the territory of muscle fibres from MG MUs and determine whether these MUs are grouped by recruitment threshold or joint action (ankle plantar flexion and knee flexion). The territory of MUs from the MG muscle varied from somewhat localized to highly distributed, with approximately half the MUs spanning at least half the length and width of the muscle. There was also no evidence of regional muscle activity based on MU recruitment thresholds or joint action. The CNS does not have the means to selectively activate regions of the MG muscle based on task requirements.

ABSTRACT

Human medial gastrocnemius (MG) motor units (MUs) are thought to occupy small muscle territories, with low-threshold units preferentially located distally. In this study, subjects (n = 8) performed ramped and sustained isometric contractions (ankle plantar flexion and knee flexion; range: ∼1-40% maximal voluntary contraction) and we measured MU territory size with spike-triggered averages from fine-wire electrodes inserted along the length (seven electrodes) or across the width (five electrodes) of the MG muscle. Of 69 MUs identified along the length of the muscle, 32 spanned at least half the muscle length (≥ 6.9 cm), 11 of which spanned all recording sites (13.6-17.9 cm). Distal fibres had smaller pennation angles (P < 0.05), which were accompanied by larger territories in MUs with fibres located distally (P < 0.05). There was no distal-to-proximal pattern of muscle activation in ramp contraction (P = 0.93). Of 36 MUs identified across the width of the muscle, 24 spanned at least half the muscle width (≥ 4.0 cm), 13 of which spanned all recording sites (8.0-10.8 cm). MUs were not localized (length or width) based on recruitment threshold or contraction type, nor was there a relationship between MU territory size and recruitment threshold (Spearman's rho = -0.20 and 0.13, P > 0.18). MUs in the human MG have larger territories than previously reported and are not localized based on recruitment threshold or joint action. This indicates that the CNS does not have the means to selectively activate regions of the MG muscle based on task requirements.

Differences in the recruitment curves obtained with monopolar and bipolar electrode configurations in the quadriceps femoris

INTRODUCTION

We sought to verify whether the stimulation intensity at which M-wave amplitude reaches a plateau actually corresponds to full motor unit activation in monopolar and bipolar configurations.

METHODS

M-waves and twitches were evoked using femoral nerve stimulation of gradually increasing intensity in 21 subjects. Recruitment curves corresponding to the amplitude of the first phase (AmpliFIRST ) and peak-to-peak amplitude (AmpliPP ) of the M-wave were obtained in the vastus lateralis, vastus medialis, and rectus femoris in monopolar and bipolar configurations.

RESULTS

In all muscles, bipolar M-waves and twitches reached plateau at a significantly lower stimulus intensity compared with monopolar M-waves (P < 0.05). The different behavior of monopolar and bipolar M-waves with stimulus intensity was found for both AmpliFIRST and AmpliPP .

CONCLUSIONS

In a bipolar configuration, the stimulus intensity at which M-waves plateau should be increased by at least 10%-20% to achieve complete motor unit recruitment. Muscle Nerve 54: 118-131, 2016.

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.

Spatial localization of electromyographic amplitude distributions associated to the activation of dorsal forearm muscles

In this study we investigated whether the spatial distribution of surface electromyographic (EMG) amplitude can be used to describe the activation of muscle portions with different biomechanical actions. Ten healthy subjects performed isometric contractions aimed to selectively activate a number of forearm muscles or muscle subportions. Monopolar electromyographic signals were collected with an electrode grid of 128 electrodes placed on the proximal, dorsal portion of the forearm. The monopolar EMG amplitude [root mean square (RMS) value] distribution was calculated for each contraction, and high-amplitude channels were identified through an automatic procedure; the position of the EMG source was estimated with the barycenter of these channels. Each of the contractions tested was associated to a specific EMG amplitude distribution, whose location in space was consistent with the expected anatomical position of the main agonist muscle (or subportion). The position of each source was significantly different from the others in at least one direction (ANOVA; transversally to the forearm: P < 0.01, F = 125.92; longitudinally: P < 0.01, F = 35.83). With such an approach, we could distinguish the spatial position of EMG distributions related to the activation of contiguous muscles [e.g., extensor carpi ulnaris (ECU) and extensor digitorum communis (EDC)], different heads of the same muscle (i.e., extensor carpi radialis (ECR) brevis and longus) and different functional compartments (i.e., EDC, middle, and ring fingers). These findings are discussed in terms of how forces along a given direction can be produced by recruiting population of motor units clustered not only in specific muscles, but also in muscle sub-portions. In addition, this study supports the use of high-density EMG systems to characterize the activation of muscle subportions with different biomechanical actions.