An examination of innervation zone movement with increases in isometric torque production

Synaptic inputs to motor neurons underlying muscle coactivation for functionally different tasks have different spectral characteristics

The central nervous system (CNS) may produce the same endpoint trajectory or torque profile with different muscle activation patterns. What differentiates these patterns is the presence of cocontraction, which does not contribute to effective torque generation but allows to modulate joints' mechanical stiffness. Although it has been suggested that the generation of force and the modulation of stiffness rely on separate pathways, a characterization of the differences between the synaptic inputs to motor neurons (MNs) underlying these tasks is still missing. In this study, participants coactivated the same pair of upper-limb muscles, i.e., the biceps brachii and the triceps brachii, to perform two functionally different tasks: limb stiffness modulation or endpoint force generation. Spike trains of MNs were identified through decomposition of high-density electromyograms (EMGs) collected from the two muscles. Cross-correlogram showed a higher synchronization between MNs recruited to modulate stiffness, whereas cross-muscle coherence analysis revealed peaks in the β-band, which is commonly ascribed to a cortical origin. These peaks did not appear during the coactivation for force generation, thus suggesting separate cortical inputs for stiffness modulation. Moreover, a within-muscle coherence analysis identified two subsets of MNs that were selectively recruited to generate force or regulate stiffness. This study is the first to highlight different characteristics, and probable different neural origins, of the synaptic inputs driving a pair of muscles under different functional conditions. We suggest that stiffness modulation is driven by cortical inputs that project to a separate set of MNs, supporting the existence of a separate pathway underlying the control of stiffness.NEW & NOTEWORTHY The characterization of the pathways underlying force generation or stiffness modulation are still unknown. In this study, we demonstrated that the common input to motor neurons of antagonist muscles shows a high-frequency component when muscles are coactivated to modulate stiffness but not to generate force. Our results provide novel insights on the neural strategies for the recruitment of multiple muscles by identifying specific spectral characteristics of the synaptic inputs underlying functionally different tasks.

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.

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.

Imaging three-dimensional innervation zone distribution in muscles from M-wave recordings

OBJECTIVE

To localize neuromuscular junctions in skeletal muscles in vivo which is of great importance in understanding, diagnosing and managing of neuromuscular disorders.

APPROACH

A three-dimensional global innervation zone imaging technique was developed to characterize the global distribution of innervation zones, as an indication of the location and features of neuromuscular junctions, using electrically evoked high-density surface electromyogram recordings.

MAIN RESULTS

The performance of the technique was evaluated in the biceps brachii of six intact human subjects. The geometric centers of the distributions of the reconstructed innervation zones were determined with a mean distance of 9.4  ±  1.4 cm from the reference plane, situated at the medial epicondyle of the humerus. A mean depth was calculated as 1.5  ±  0.3 cm from the geometric centers to the closed points over the skin. The results are consistent with those reported in previous histology studies. It was also found that the volumes and distributions of the reconstructed innervation zones changed as the stimulation intensities increased until the supramaximal muscle response was achieved.

SIGNIFICANCE

Results have demonstrated the high performance of the proposed imaging technique in noninvasively imaging global distributions of the innervation zones in the three-dimensional muscle space in vivo, and the feasibility of its clinical applications, such as guiding botulinum toxin injections in spasticity management, or in early diagnosis of neurodegenerative progression of amyotrophic lateral sclerosis.

Re-evaluation of EMG-torque relation in chronic stroke using linear electrode array EMG recordings

The objective was to re-evaluate the controversial reports of EMG-torque relation between impaired and non-impaired sides using linear electrode array EMG recordings. Ten subjects with chronic stroke performed a series of submaximal isometric elbow flexion tasks. A 20-channel linear array was used to record surface EMG of the biceps brachii muscles from both impaired and non-impaired sides. M-wave recordings for bilateral biceps brachii muscles were also made. Distribution of the slope of the EMG-torque relations for the individual channels showed a quasi-symmetrical "M" shaped pattern. The lowest value corresponded to the innervation zone (IZ) location. The highest value from the slope curve for each side was selected for comparison to minimize the effect of electrode placement and IZ asymmetry. The slope was greater on the impaired side in 4 of 10 subjects. There were a weak correlation between slope ratio and strength ratio and a moderate to high correlation between slope ratio and M-wave ratio between two sides. These findings suggest that the EMG-torque relations are likely mediated and influenced by multiple factors. Our findings emphasize the importance of electrode placement and suggest the primary role of peripheral adaptive changes in the EMG-torque relations in chronic stroke.

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.

The influence of electromyographic recording methods and the innervation zone on the mean power frequency-torque relationships

This study examined the effects of electromyographic (EMG) recording methods and innervation zone (IZ) on the mean power frequency (MPF)-torque relationships. Nine subjects performed isometric ramp muscle actions of the leg extensors from 5% to 100% of maximal voluntary contraction with an eight channel linear electrode array over the IZ of the vastus lateralis. The slopes were calculated from the log-transformed monopolar and bipolar EMG MPF-torque relationships for each channel and subject and 95% confidence intervals (CI) were constructed around the slopes for each relationship and the composite of the slopes. Twenty-two to 55% of the subjects exhibited 95% CIs that did not include a slope of zero for the monopolar EMG MPF-torque relationships while 25-75% of the subjects exhibited 95% CIs that did not include a slope of zero for the bipolar EMG MPF-torque relationships. The composite of the slopes from the EMG MPF-torque relationships were not significantly different from zero for any method or channel, however, the method and IZ location slightly influenced the number of significant slopes on a subject-by-subject basis. The log-transform model indicated that EMG MPF-torque patterns were nonlinear regardless of recording method or distance from the IZ.

Predicting electromyographic signals under realistic conditions using a multiscale chemo-electro-mechanical finite element model

This paper presents a novel multiscale finite element-based framework for modelling electromyographic (EMG) signals. The framework combines (i) a biophysical description of the excitation-contraction coupling at the half-sarcomere level, (ii) a model of the action potential (AP) propagation along muscle fibres, (iii) a continuum-mechanical formulation of force generation and deformation of the muscle, and (iv) a model for predicting the intramuscular and surface EMG. Owing to the biophysical description of the half-sarcomere, the model inherently accounts for physiological properties of skeletal muscle. To demonstrate this, the influence of membrane fatigue on the EMG signal during sustained contractions is investigated. During a stimulation period of 500 ms at 100 Hz, the predicted EMG amplitude decreases by 40% and the AP propagation velocity decreases by 15%. Further, the model can take into account contraction-induced deformations of the muscle. This is demonstrated by simulating fixed-length contractions of an idealized geometry and a model of the human tibialis anterior muscle (TA). The model of the TA furthermore demonstrates that the proposed finite element model is capable of simulating realistic geometries, complex fibre architectures, and can include different types of heterogeneities. In addition, the TA model accounts for a distributed innervation zone, different fibre types and appeals to motor unit discharge times that are based on a biophysical description of the α motor neurons.

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.

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.

Accuracy of three different techniques for automatically estimating innervation zone location

The purpose of this study was to compare the accuracy of the estimated innervation zone (IZ) locations obtained from cross-correlation, the minimum amplitude, and maximum center frequency criteria. Eight healthy men (mean±SD age=23.0±4.3 yrs) performed isometric muscle actions of the leg extensors, and 15 separate bipolar surface electromyographic (EMG) signals were detected from the vastus lateralis. A custom software program was used to estimate the location of the IZ based on: (1) the EMG channel that demonstrated the lowest amplitude, (2) the EMG channel that showed the highest mean frequency, and (3) the EMG channel that demonstrated the lowest peak cross-correlation between the signals from adjacent channels. The IZ location estimates from the lowest amplitude and highest mean frequency criteria were accurate in only 43.75% and 7.5% of the cases, respectively. The accuracy of the cross-correlation-based method was 90%. The cross-correlation-based method was much more accurate for estimating IZ location than were the lowest amplitude and highest mean frequency criteria. Cross-correlation could potentially be used for estimating the location of the IZ without the need for visual inspection of EMG signals.

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.