Motor unit pool organization examined via spike-triggered averaging of the surface electromyogram

Developmental changes in motor unit activity patterns: child-adult comparison using discrete motor unit analysis

Using global surface electromyography (sEMG) and the sEMG threshold it has been suggested that children activate their type-II motor unit (MU) to a lesser extent compared with adults. However, when age-related differences in discrete MU activation are examined using sEMG decomposition this phenomenon is not observed. Furthermore, findings from these studies are inconsistent and conflicting. Therefore, the purpose of this study was to examine differences in discrete MU activation of the vastus lateralis (VL) between boys and men during moderate-intensity knee extensions. Seventeen boys and 20 men completed two laboratory sessions. Following a habituation session, maximal voluntary isometric knee extension (MVIC) torque was determined before completing trapezoidal contractions at 70% MVIC. sEMG of the VL was captured and mathematically decomposed into individual MU action potential trains. Motor unit action potential amplitude (MUAPamp), recruitment threshold (RT), and MU firing rates (MUFR) were calculated. We observed that MUAPamp-RT slope was steeper in men compared with boys (p < 0.05) even after accounting for fat thickness and quadriceps muscle depth. The mean MUFR and y-intercept of the MUFR-RT relationship were significantly (p < 0.001) lower in boys than in men. The slope of the MUFR-RT relationship tended to be steeper in men, but the differences did not reach statistical significance (p = 0.056). Overall, our results suggest that neural strategies used to produce torque are different among boys and men. Such differences may be related, in part, to boys' lower MUFR and lesser ability to activate their higher-threshold MUs. Although, other factors (e.g., muscle composition) likely also play a role.

Influence of Blood Flow Restriction on Neuromuscular Function and Fatigue During Forearm Flexion in Men

ABSTRACT

Montgomery, TR Jr, Olmos, A, Sears, KN, Succi, PJ, Hammer, SM, Bergstrom, HC, Hill, EC, Trevino, MA, and Dinyer-McNeely, TK. Influence of blood flow restriction on neuromuscular function and fatigue during forearm flexion in men. J Strength Cond Res 38(7): e349-e358, 2024-To determine the effects of blood flow restriction (BFR) on the mean firing rate (MFR) and motor unit action potential amplitude (MUAPAMP) vs. recruitment threshold (RT) relationships during fatiguing isometric elbow flexions. Ten men (24.5 ± 4.0 years) performed isometric trapezoidal contractions at 50% maximum voluntary contraction to task failure with or without BFR, on 2 separate days. For BFR, a cuff was inflated to 60% of the pressure required for full brachial artery occlusion at rest. During both visits, surface electromyography was recorded from the biceps brachii of the dominant limb and the signal was decomposed. A paired-samples t test was used to determine the number of repetitions completed between BFR and CON. ANOVAs (repetition [first, last] × condition [BFR, CON]) were used to determine differences in MFR vs. RT and MUAPAMP vs. RT relationships. Subjects completed more repetitions during CON (12 ± 4) than BFR (9 ± 2; p = 0.012). There was no significant interaction (p > 0.05) between the slopes and y-intercepts during the repetition × condition interaction for MUAPAMP vs. MFR. However, there was a main effect of repetition for the slopes of the MUAPAMP vs. RT (p = 0.041) but not the y-intercept (p = 0.964). Post hoc analysis (collapsed across condition) indicated that the slopes of the MUAPAMP vs. RT during the first repetition was less than the last repetition (first: 0.022 ± 0.003 mv/%MVC; last: 0.028 ± 0.004 mv/%MVC; p = 0.041). Blood flow restriction resulted in the same amount of higher threshold MU recruitment in approximately 75% of the repetitions. Furthermore, there was no change in MFR for either condition, even when taken to task failure. Thus, BFR training may create similar MU responses with less total work completed than training without BFR.

Antagonist Activation Measurement in Triceps Surae Using High-Density and Bipolar Surface EMG in Chronic Hemiparesis

After a stroke, antagonist muscle activation during agonist command impedes movement. This study compared measurements of antagonist muscle activation using surface bipolar EMG in the gastrocnemius medialis (GM) and high-density (HD) EMG in the GM and soleus (SO) during isometric submaximal and maximal dorsiflexion efforts, with knee flexed and extended, in 12 subjects with chronic hemiparesis. The coefficients of antagonist activation (CAN) of GM and SO were calculated according to the ratio of the RMS amplitude during dorsiflexion effort to the maximal agonist effort for the same muscle. Bipolar CAN (BipCAN) was compared to CAN from channel-specific (CsCAN) and overall (OvCAN) normalizations of HD-EMG. The location of the CAN centroid was explored in GM, and CAN was compared between the medial and lateral portions of SO. Between-EMG system differences in GM were observed in maximal efforts only, between BipCAN and CsCAN with lower values in BipCAN (p < 0.001), and between BipCAN and OvCAN with lower values in OvCAN (p < 0.05). The CAN centroid is located mid-height and medially in GM, while the CAN was similar in medial and lateral SO. In chronic hemiparesis, the estimates of GM hyperactivity differ between bipolar and HD-EMGs, with channel-specific and overall normalizations yielding, respectively, higher and lower CAN values than bipolar EMG. HD-EMG would be the way to develop personalized rehabilitation programs based on individual antagonist activations.

Force and Neuromuscular Responses During a Handgrip Hold to Failure Anchored to a Moderate Perceptual Intensity in Males

OBJECTIVES

The current study investigated performance fatigability (PF) and time course of changes in force, electromyographic amplitude (EMG AMP) and frequency (EMG MPF), and neuromuscular efficiency (NME) during a sustained, isometric, handgrip hold to failure (HTF) using the rating of perceived exertion (RPE)-Clamp Model.

METHODS

Twelve males performed a handgrip HTF anchored to RPE=5. The time to task failure (Tlim), force (N), EMG AMP and MPF, and NME (normalized force/ normalized EMG AMP) were recorded. Analyses included a paired samples t-test for PF at an alpha of p<0.05, 1-way repeated measures ANOVA across time and post-hoc t-tests (p<0.0025) for force, EMG AMP and MPF, and NME responses.

RESULTS

The PF (pre- to post- maximal force % decline) was 38.2±11.5%. There were decreases in responses, relative to 0% Tlim, from 40% to 100% Tlim (force), at 30%, 60%, and 100% Tlim (EMG AMP), from 10% to 100% Tlim(EMP MPF), and from 50% to 65%, and 80% to 100% Tlim (NME) (p<0.0025).

CONCLUSIONS

The RPE-Clamp Model in this study demonstrated that pacing strategies may be influenced by the integration of anticipatory, feedforward, and feedback mechanisms, and provided insights into the relationship between neuromuscular and perceptual responses, and actual force generating capacity.

Motor unit firing rates during slow and fast contractions in boys and men

BACKGROUND

Motor unit (MU) activation during maximal contractions is lower in children compared with adults. Among adults, discrete MU activation differs, depending on the rate of contraction. We investigated the effect of contraction rate on discrete MU activation in boys and men.

METHODS

Following a habituation session, 14 boys and 20 men completed two experimental sessions for knee extension and wrist flexion, in random order. Maximal voluntary isometric torque (MVIC) was determined before completing trapezoidal isometric contractions (70%MVIC) at low (10%MVIC/s) and high (35%MVIC/s) contraction rates. Surface electromyography was captured from the vastus lateralis (VL) and flexor carpi radialis (FCR) and decomposed into individual MU action potential (MUAP) trains.

RESULTS

In both groups and muscles, the initial MU firing rate (MUFR) was greater (p < 0.05) at high compared with low contraction rates. The increase in initial MUFR at the fast contraction in the VL was greater in men than boys (p < 0.05). Mean MUFR was significantly lower during fast contractions only in the FCR (p < 0.05). In both groups and muscles, the rate of decay of MUFR with increasing MUAP amplitude was less steep (p < 0.05) during fast compared with slow contractions.

CONCLUSION

In both groups and muscles, initial MUFRs, as well as MUFRs of large MUs were higher during fast compared with slow contractions. However, in the VL, the increase in initial MUFR was greater in men compared with boys. This suggests that in large muscles, men may rely more on increasing MUFR to generate torque at faster rates compared with boys.

Spectral characterization of human leg EMG signals from an open access dataset for the development of computational models

Large-scale neuromusculoskeletal models have been used for predicting mechanisms underlying neuromuscular functions in humans. Simulations of such models provide several types of signals of practical interest, such as surface electromyographic signals (EMG), which are compared with experimental data for interpretations of neurophysiological phenomena under study. Specifically, realistic characterization of spectral properties of simulated EMG signals is important for achieving powerful inferences, whereas considerations should be taken for myoelectric signals of different muscles. In this study, we characterized spectral properties of surface interference pattern EMG signals and motor unit action potentials (MUAP) acquired from three plantar flexor muscles: Soleus (SO), Medial Gastrocnemius (MG), and Lateral Gastrocnemius (LG); and one dorsiflexor muscle: Tibialis Anterior (TA). Surface EMG signals were acquired from 20 participants using the same convention for electrode placement. Specifically, interference pattern EMG signals were obtained during isometric constant force contractions at 5%, 10% and 20% of maximum voluntary contraction (MVC), whereas surface MUAPs were decomposed from surface EMG signals obtained at low contraction forces. We compared the spectrum median frequency (MDF) estimated from interference pattern EMG signals across muscles and contraction intensities. Additionally, we compared MDF and durations of MUAPs between muscles. Our results showed that MDF of interference pattern EMG signals acquired from TA were higher compared to SO, MG, and LG for all contraction intensities i.e., 5%, 10%, and 20% MVC. Consistently, MUAPs acquired from TA also had higher MDF values and shorter durations compared to the other leg muscles. We provide herein a dataset with the surface MUAPs waveforms and interference pattern EMG signals obtained for this study, which should be useful for implementing and validating the simulation of myoelectrical signals of leg muscles. Importantly, these results indicate that spectral properties of myoelectrical signals should be considered for improving EMG modeling in large-scale neuromusculoskeletal models.

Motor unit firing rates increase in prepubescent youth following linear periodization resistance exercise training

PURPOSE

The purpose was to examine the effects of 8-weeks (3 days/week) of linear periodization resistance exercise training (RET) on neuromuscular function in prepubescent youth.

METHODS

Twenty-five healthy prepubescent youth (11 males, 14 females, age = 9.1 ± 0.8 years) completed the RET (n = 17) or served as controls (CON, n = 8). Isometric maximal voluntary contractions (MVCs) and trapezoidal submaximal contractions at 35 and 60% MVC of the right leg extensors were performed with surface electromyography (EMG) recorded from the leg extensors [vastus lateralis (VL), rectus femoris, and vastus medialis] and flexors (biceps femoris and semitendinosus). EMG amplitude of the leg extensors and flexors were calculated during the MVCs. Motor unit (MU) action potential trains were decomposed from the surface EMG of the VL for the 35 and 60% MVCs. MU firing rates and action potential amplitudes were regressed against recruitment threshold with the y-intercepts and slopes calculated for each contraction. Total leg extensor muscle cross-sectional area (CSA) was collected using ultrasound images. ANOVA models were used to examine potential differences.

RESULTS

Isometric strength increased post-RET (P = 0.006) with no changes in leg extensor and flexor EMG amplitude. Furthermore, there were no changes in total CSA or the MU action potential amplitude vs. recruitment threshold relationships. However, there were increases in the firing rates of the higher-threshold MUs post-RET as indicated with greater y-intercepts (P = 0.003) from the 60% MVC and less negative slope (P = 0.004) of the firing rates vs. recruitment threshold relationships at 35% MVC.

CONCLUSIONS

MU adaptations contribute to strength increases following RET in prepubescent youth.

Motor unit activity and synaptic inputs to motoneurons in the caudal part of the injured spinal cord

Spinal cord injury (SCI) disrupts neuronal function below the lesion epicenter, causing disuse muscle atrophy. We investigated motor unit (MU) activity and synaptic inputs to motoneurons in the caudal region of the injured spinal cord. Participants with C4-C7 cervical injuries were studied. The extensor digitorum communis (EDC) muscle, which is mainly innervated by C8, was assessed for disuse muscle atrophy. Using advanced electromyography and signal-processing techniques, we examined the concurrent activation of a substantial population of MUs during force-tracking tasks. We found that in participants with SCI (n = 9), both MU discharge rates and the amplitudes of MU action potentials were significantly lower than in controls (n = 9). After SCI, MUs were recruited in a limited force range as the strength of muscle contractions increased, implying a disruption in the orderly MU recruitment pattern. Coherence analysis revealed reduced synaptic inputs to motoneurons in the delta band (0.5-5 Hz) for participants with SCI, suggesting diminished common synaptic inputs to the EDC muscle. In addition, participants with SCI exhibited greater muscle force variability. Using principal component analysis on low-frequency MU discharge rates, we found that the first common component (FCC) captured the most discharge variability in participants with SCI. The coefficients of variation (CV) of the FCC correlated with force signal CVs, suggesting force variability mainly results from common synaptic inputs to the EDC muscle after SCI. These results advance our understanding of the neurophysiology of disuse muscle atrophy in human SCI, paving the way for therapeutic interventions to restore muscle function.NEW & NOTEWORTHY This study analyzed motor unit (MU) function below the lesion epicenter in patients with spinal cord injury (SCI). We found reduced MU discharge rates and action potential amplitudes in participants with SCI compared with controls. The strength of common synaptic inputs to motoneurons was reduced in patients with SCI, with increased force variability primarily due to low-frequency oscillations of common inputs. This study enhances understanding of neurophysiological and behavioral changes in disuse muscle atrophy post-SCI.

Modified motor unit properties in residual muscle following transtibial amputation

Objective.Neural signals in residual muscles of amputated limbs are frequently decoded to control powered prostheses. Yet myoelectric controllers assume muscle activities of residual muscles are similar to that of intact muscles. This study sought to understand potential changes to motor unit (MU) properties after limb amputation.Approach.Six people with unilateral transtibial amputation were recruited. Surface electromyogram (EMG) of residual and intacttibialis anterior(TA) andgastrocnemius(GA) muscles were recorded while subjects traced profiles targeting up to 20% and 35% of maximum activation for each muscle (isometric for intact limbs). EMG was decomposed into groups of MU spike trains. MU recruitment thresholds, action potential amplitudes (MU size), and firing rates were correlated to model Henneman's size principle, the onion-skin phenomenon, and rate-size associations. Organization (correlation) and modulation (rates of change) of relations were compared between intact and residual muscles.Main results.The residual TA exhibited significantly lower correlation and flatter slopes in the size principle and onion-skin, and each outcome covaried between the MU relations. The residual GA was unaffected for most subjects. Subjects trained prior with myoelectric prostheses had minimally affected slopes in the TA. Rate-size association correlations were preserved, but both residual muscles exhibited flatter decay rates.Significance.We showed peripheral neuromuscular damage also leads to spinal-level functional reorganizations. Our findings suggest models of MU recruitment and discharge patterns for residual muscle EMG generation need reparameterization to account for disturbances observed. In the future, tracking MU pool adaptations may also provide a biomarker of neuromuscular control to aid training with myoelectric prostheses.

Ankle Torque Estimation With Motor Unit Discharges in Residual Muscles Following Lower-Limb Amputation

There has been increased interest in using residual muscle activity for neural control of powered lower-limb prostheses. However, only surface electromyography (EMG)-based decoders have been investigated. This study aims to investigate the potential of using motor unit (MU)-based decoding methods as an alternative to EMG-based intent recognition for ankle torque estimation. Eight people without amputation (NON) and seven people with amputation (AMP) participated in the experiments. Subjects conducted isometric dorsi- and plantarflexion with their intact limb by tracing desired muscle activity of the tibialis anterior (TA) and gastrocnemius (GA) while ankle torque was recorded. To match phantom limb and intact limb activity, AMP mirrored muscle activation with their residual TA and GA. We compared neuromuscular decoders (linear regression) for ankle joint torque estimation based on 1) EMG amplitude (aEMG), 2) MU firing frequencies representing neural drive (ND), and 3) MU firings convolved with modeled twitch forces (MUDrive). In addition, sensitivity analysis and dimensionality reduction of optimization were performed on the MUDrive method to further improve its practical value. Our results suggest MUDrive significantly outperforms (lower root-mean-square error) EMG and ND methods in muscles of NON, as well as both intact and residual muscles of AMP. Reducing the number of optimized MUDrive parameters degraded performance. Even so, optimization computational time was reduced and MUDrive still outperformed aEMG. Our outcomes indicate integrating MU discharges with modeled biomechanical outputs may provide a more accurate torque control signal than direct EMG control of assistive, lower-limb devices, such as exoskeletons and powered prostheses.

A comparison of techniques for verifying the accuracy of precision decomposition-derived relationships between motor unit firing rates and recruitment thresholds from surface EMG signals

Approaches for validating motor unit firing times following surface electromyographic (EMG) signal decomposition with the precision decomposition III (PDIII) algorithm have not been agreed upon. Two approaches have been common: (1) "reconstruct-and-test" and (2) spike-triggered averaging (STA). We sought to compare motor unit results following the application of these approaches. Surface EMG signals were recorded from the vastus lateralis of 13 young males performing trapezoidal, isometric knee extensions at 50% and 80% of maximum voluntary contraction (MVC) force. The PDIII algorithm was used to quantify motor unit firing rates. Motor units were excluded using eight combinations of the reconstruct-and-test approach with accuracy thresholds of 0, 90, 91, and 92% with and without STA. The mean firing rate versus recruitment threshold relationship was minimally affected by STA. At 80% MVC, slopes acquired at the 0% accuracy threshold were significantly greater (i.e., less negative) than when 91% (p = .010) and 92% (p = .030) accuracy thresholds were applied. The application of STA has minimal influence on surface EMG signal decomposition results. Stringent reconstruct-and-test accuracy thresholds influence motor unit-derived relationships at high forces, perhaps explained through the increased presence of large motor unit action potentials. Investigators using the PDIII algorithm can expect negligible changes in motor unit-derived linear regression relationships with the application of secondary validation procedures.

Sex-related differences in motor unit behavior are influenced by myosin heavy chain during high- but not moderate-intensity contractions

AIMS

Motor unit recruitment and firing rate patterns of the vastus lateralis (VL) have not been compared between sexes during moderate- and high-intensity contraction intensities. Additionally, the influence of fiber composition on potential sex-related differences remains unquantified.

METHODS

Eleven males and 11 females performed 40% and 70% maximal voluntary contractions (MVCs). Surface electromyographic (EMG) signals recorded from the VL were decomposed. Recruitment thresholds (RTs), MU action potential amplitudes (MUAPAMP ), initial firing rates (IFRs), mean firing rates (MFRs), and normalized EMG amplitude (N-EMGRMS ) at steady torque were analyzed. Y-intercepts and slopes were calculated for MUAPAMP , IFR, and MFR versus RT relationships. Type I myosin heavy chain isoform (MHC) was determined with muscle biopsies.

RESULTS

There were no sex-related differences in MU characteristics at 40% MVC. At 70% MVC, males exhibited greater slopes (p = 0.002) for the MUAPAMP , whereas females displayed greater slopes (p = 0.001-0.007) for the IFR and MFR versus RT relationships. N-EMGRMS at 70% MVC was greater for females (p < 0.001). Type I %MHC was greater for females (p = 0.006), and was correlated (p = 0.018-0.031) with the slopes for the MUAPAMP , IFR, and MFR versus RT relationships at 70% MVC (r = -0.599-0.585).

CONCLUSION

Both sexes exhibited an inverse relationship between MU firing rates and recruitment thresholds. However, the sex-related differences in MU recruitment and firing rate patterns and N-EMGRMS at 70% MVC were likely due to greater type I% MHC and smaller twitch forces of the higher threshold MUs for the females. Evidence is provided that muscle fiber composition may explain divergent MU behavior between sexes.

Transformer-based hand gesture recognition from instantaneous to fused neural decomposition of high-density EMG signals

Designing efficient and labor-saving prosthetic hands requires powerful hand gesture recognition algorithms that can achieve high accuracy with limited complexity and latency. In this context, the paper proposes a Compact Transformer-based Hand Gesture Recognition framework referred to as [Formula: see text], which employs a vision transformer network to conduct hand gesture recognition using high-density surface EMG (HD-sEMG) signals. Taking advantage of the attention mechanism, which is incorporated into the transformer architectures, our proposed [Formula: see text] framework overcomes major constraints associated with most of the existing deep learning models such as model complexity; requiring feature engineering; inability to consider both temporal and spatial information of HD-sEMG signals, and requiring a large number of training samples. The attention mechanism in the proposed model identifies similarities among different data segments with a greater capacity for parallel computations and addresses the memory limitation problems while dealing with inputs of large sequence lengths. [Formula: see text] can be trained from scratch without any need for transfer learning and can simultaneously extract both temporal and spatial features of HD-sEMG data. Additionally, the [Formula: see text] framework can perform instantaneous recognition using sEMG image spatially composed from HD-sEMG signals. A variant of the [Formula: see text] is also designed to incorporate microscopic neural drive information in the form of Motor Unit Spike Trains (MUSTs) extracted from HD-sEMG signals using Blind Source Separation (BSS). This variant is combined with its baseline version via a hybrid architecture to evaluate potentials of fusing macroscopic and microscopic neural drive information. The utilized HD-sEMG dataset involves 128 electrodes that collect the signals related to 65 isometric hand gestures of 20 subjects. The proposed [Formula: see text] framework is applied to 31.25, 62.5, 125, 250 ms window sizes of the above-mentioned dataset utilizing 32, 64, 128 electrode channels. Our results are obtained via 5-fold cross-validation by first applying the proposed framework on the dataset of each subject separately and then, averaging the accuracies among all the subjects. The average accuracy over all the participants using 32 electrodes and a window size of 31.25 ms is 86.23%, which gradually increases till reaching 91.98% for 128 electrodes and a window size of 250 ms. The [Formula: see text] achieves accuracy of 89.13% for instantaneous recognition based on a single frame of HD-sEMG image. The proposed model is statistically compared with a 3D Convolutional Neural Network (CNN) and two different variants of Support Vector Machine (SVM) and Linear Discriminant Analysis (LDA) models. The accuracy results for each of the above-mentioned models are paired with their precision, recall, F1 score, required memory, and train/test times. The results corroborate effectiveness of the proposed [Formula: see text] framework compared to its counterparts.

Differences in psychological readiness for return to sport after anterior cruciate ligament injury is evident in thigh musculature motor unit characteristics

Background

Following anterior cruciate ligament (ACL) injury, many athletes that undergo surgery and 6-9 months of rehabilitation struggle to return to sport. Evidence suggests that psychological factors contribute to this failure to return-to-sport.

Objective

Determine the motor control relationship between thigh musculature motor unit characteristics and psychological readiness to return to sport between ACL-injured and healthy controls.

Study design

A longitudinal cohort study.

Methods

Athletes longitudinally completed the ACL Return to Sport after Injury (ACL-RSI) survey and isometric strength measures with a measurement of electromyography (EMG) of the vastus lateralis, vastus medialis, biceps femoris, and semitendinosus. A score cut-off of 61 on the ACL-RSI was used to divide ACL-injured groups. EMG was decomposed to provide each identified motor unit's characteristics (amplitude, average firing rate, etc).

Results

Data demonstrated increased average firing rate for hamstrings (p<0.001), decreased average firing rate for vastus lateralis (p<0.001) and decreased motor unit size for both the quadriceps and hamstrings at return-to-sport post-ACL reconstruction compared with sex-matched and age-matched healthy controls (p<0.001). Furthermore, there were marked differences in disparate ACL-RSI scores between ACL-injured athletes.

Conclusions

At return to sport, ACL-injured athletes have major alterations of thigh musculature motor control, with smaller motor units used by those with low ACL-RSI scores. This study uniquely demonstrates objective thigh muscle motor unit characteristics that coincide with subjective reports of psychological readiness. This information will be important to address psychomotor complexes of injury for future rehabilitation protocols.

Arthrogenic muscle inhibition manifests in thigh musculature motor unit characteristics after anterior cruciate ligament injury

Joint trauma induces a presynaptic reflex inhibition termed arthrogenic muscle inhibition (AMI) that prevents complete activation of muscles. Reduced motor unit (MU) output is a hypothesised mechanism for persistent strength deficits. The objective of this study was to determine MU characteristics of thigh musculature and determine how they change with anterior cruciate ligament (ACL) injury compared to healthy controls. A randomised protocol of knee flexion/extension isometric contractions (10-50% maximal voluntary isometric contraction) was performed for each leg with surface EMG 5-pin array electrodes placed on the vastus medialis, vastus lateralis, semitendinosus and biceps femoris. Longitudinal assessments for average rate coding, recruitment thresholds and MU action potentials were acquired at 6-month intervals. With exception of the vastus medialis, all thigh musculature of ACL-injured demonstrated smaller MU action potential peak-to-peak amplitude. For average rate coding, ACL-injured demonstrated lower coding rates than Controls for the quadriceps (p < .05) and higher rates than Controls for the hamstrings (p < .05). These MU characteristics were different from Controls after ACL reconstruction up to 12 months post-surgery, yet maximal strength increased during this time frame. As thigh MU characteristics are known across phases of ACL rehabilitation, future studies can assess these patterns of motor control and their potential to determine risk of re-injury. Further, future rehabilitation can target specific intervention programmes to restore motor control.HighlightsMotor unit strategies of arthrogenic muscle inhibition are characterised for the first time via decomposed EMG.Motor unit deficits of thigh musculature persist throughout all phases of ACL rehabilitation, even after return-to-sport.After ACL injury, motor unit sizes at similar recruitment thresholds were smaller than those of healthy controls.

The effects of resistance training to near failure on strength, hypertrophy, and motor unit adaptations in previously trained adults

Limited research exists examining how resistance training to failure affects applied outcomes and single motor unit characteristics in previously trained individuals. Herein, resistance-trained adults (24 ± 3 years old, self-reported resistance training experience was 6 ± 4 years, 11 men and 8 women) were randomly assigned to either a low-repetitions-in-reserve (RIR; i.e., training near failure, n = 10) or high-RIR (i.e., not training near failure, n = 9) group. All participants implemented progressive overload during 5 weeks where low-RIR performed squat, bench press, and deadlift twice weekly and were instructed to end each training set with 0-1 RIR. high-RIR performed identical training except for being instructed to maintain 4-6 RIR after each set. During week 6, participants performed a reduced volume-load. The following were assessed prior to and following the intervention: (i) vastus lateralis (VL) muscle cross-sectional area (mCSA) at multiple sites; (ii) squat, bench press, and deadlift one-repetition maximums (1RMs); and (iii) maximal isometric knee extensor torque and VL motor unit firing rates during an 80% maximal voluntary contraction. Although RIR was lower in the low- versus high-RIR group during the intervention (p < 0.001), total training volume did not significantly differ between groups (p = 0.222). There were main effects of time for squat, bench press, and deadlift 1RMs (all p-values < 0.05), but no significant condition × time interactions existed for these or proximal/middle/distal VL mCSA data. There were significant interactions for the slope and y-intercept of the motor unit mean firing rate versus recruitment threshold relationship. Post hoc analyses indicated low-RIR group slope values decreased and y-intercept values increased after training suggesting low-RIR training increased lower-threshold motor unit firing rates. This study provides insight into how resistance training in proximity to failure affects strength, hypertrophy, and single motor unit characteristics, and may inform those who aim to program for resistance-trained individuals.

Asynchronous axonal firing patterns evoked via continuous subthreshold kilohertz stimulation

Objective.Transcutaneous electrical stimulation of peripheral nerves is a common technique to assist or rehabilitate impaired muscle activation. However, conventional stimulation paradigms activate nerve fibers synchronously with action potentials time-locked with stimulation pulses. Such synchronous activation limits fine control of muscle force due to synchronized force twitches. Accordingly, we developed a subthreshold high-frequency stimulation waveform with the goal of activating axons asynchronously.Approach.We evaluated our waveform experimentally and through model simulations. During the experiment, we delivered continuous subthreshold pulses at frequencies of 16.67, 12.5, or 10 kHz transcutaneously to the median and ulnar nerves. We obtained high-density electromyographic (EMG) signals and fingertip forces to quantify the axonal activation patterns. We used a conventional 30 Hz stimulation waveform and the associated voluntary muscle activation for comparison. We modeled stimulation of biophysically realistic myelinated mammalian axons using a simplified volume conductor model to solve for extracellular electric potentials. We compared the firing properties under kHz and conventional 30 Hz stimulation.Main results.EMG activity evoked by kHz stimulation showed high entropy values similar to voluntary EMG activity, indicating asynchronous axon firing activity. In contrast, we observed low entropy values in EMG evoked by conventional 30 Hz stimulation. The muscle forces evoked by kHz stimulation also showed more stable force profiles across repeated trials compared with 30 Hz stimulation. Our simulation results provide direct evidence of asynchronous firing patterns across a population of axons in response to kHz frequency stimulation, while 30 Hz stimulation elicited synchronized time-locked responses across the population.Significance.We demonstrate that the continuous subthreshold high-frequency stimulation waveform can elicit asynchronous axon firing patterns, which can lead to finer control of muscle forces.

Motor Unit Number Estimation in Spastic Biceps Brachii Muscles of Chronic Stroke Survivors Before and After BoNT Injection

OBJECTIVE

The study aims to characterize the motor unit (MU) loss in spastic biceps brachii muscle (BBM) of chronic stroke survivors before and after botulinum neurotoxin (BoNT) injection.

METHODS

High-density weighted average (HDWA) motor unit number estimation (MUNE) was employed to estimate the number of functioning motor units of BBMs of eight chronic stroke survivors 1-week before (1st visit) and 3-week after (2nd visit) BoNT injection based on the surface electromyography (sEMG) signals recorded during voluntary contraction and supramaximal electrical stimulation.

RESULT

Significant lower MUNE was estimated from the spastic BBMs compared to the non-spastic MUNEs during two visits. A surprisingly higher MUNE was obtained from the spastic side during the 2nd visit after BoNT injection.

CONCLUSIONS

The HDWA MUNE technique can be employed to characterize the motor unit loss in spastic muscle caused by upper motor neuro lesions at contraction level up to 30% MVC, but may fail to detect the MU loss caused by the chemodenervation effect of BoNT due to the non-uniform denervation of small and large size MUs.

SIGNIFICANCE

This study presents the first effort to evaluate the applicability of HDWA MUNE technique to characterize the MU loss in the spastic muscle following stroke and the subsequent BoNT injection for the treatment of post-stroke spasticity. The finding of this study suggests that HDWA MUNE can be a sensitive approach to detect the MU loss in spastic muscles after stroke, but the large inter-subject MUNE variability after the BoNT injection should be interpreted with caution.

Arthrogenic muscle inhibition after anterior cruciate ligament injury: Injured and uninjured limb recovery over time

Introduction

It is well documented that marked weakness of the quadriceps is present after knee joint injury. This joint trauma induces a presynaptic reflex inhibition of musculature surrounding the joint, termed arthrogenic muscle inhibition (AMI). The extent to which anterior cruciate ligament (ACL) injury affects thigh musculature motor unit activity, which may affect restoration of thigh muscle strength after injury, is undetermined.

Methods

A randomized protocol of knee flexion and extension isometric contractions (10%-50% maximal voluntary isometric contraction) were performed for each leg on 54 subjects with electromyography array electrodes placed on the vastus medialis, vastus lateralis, semitendinosus, and biceps femoris. Longitudinal assessments for motor unit recruitment and average firing rate were acquired at 6-month intervals for 1 year post ACL injury.

Results

The ACL-injured population demonstrated smaller quadriceps and hamstrings motor unit size (assessed via motor unit action potential peak-to-peak amplitude) and altered firing rate activity in both injured and uninjured limbs compared to healthy controls. Motor unit activity remained altered compared to healthy controls at 12 months post ACL reconstruction (ACLR).

Discussion

Motor unit activity was altered after ACLR up to 12 months post-surgery. Further research is warranted to optimize rehabilitation interventions that adequately address altered motor unit activity and improve safety and success with return to sport after ACLR. In the interim, evidence based clinical reasoning with a focus on development of muscular strength and power capacity should be the impetus behind rehabilitation programming to address motor control deficits.

Motor unit distribution and recruitment in spastic and non-spastic bilateral biceps brachii muscles of chronic stroke survivors

Objective.This study aims to characterize the motor units (MUs) distribution and recruitment pattern in the spastic and non-spastic bilateral biceps brachii muscles (BBMs) of chronic stroke survivors.Approach.High-density surface electromyography (HD-sEMG) signals were collected from both spastic and non-spastic BBMs of fourteen chronic stroke subjects during isometric elbow flexion at 10%, 30%, 50% and 100% maximal voluntary contractions (MVCs). By combining HD-sEMG decomposition and bioelectrical source imaging, MU innervation zones (MUIZs) of the decomposed MUs were first localized in the 3D space of spastic and non-spastic BBMs. The MU depth defined as the distance between the localized MUIZ and its normal projection on the skin surface was then normalized to the arm radius of each subject and averaged at given contraction level. The averaged MU depth at different contraction levels on a specific arm side (intra-side) and the bilateral depths under a specific contraction level (inter-side) were compared.Main results.The average depth of decomposed MUs increased with the contraction force and significant differences observed between 10% vs 50% (p< 0.0001), 10% vs 100% (p< 0.0001) and 30% vs 100% MVC (p= 0.0017) on the non-spastic side, indicating that larger MUs with higher recruitment threshold locate in deeper muscle regions. In contrast, no force-related difference in MU depth was observed on the spastic side, suggesting a disruption of orderly recruitment of MUs with increase of force level, or the MU denervation and the subsequent collateral reinnervation secondary to upper motor neuron lesions. Inter-side comparison demonstrated significant MU depth difference at 10% (p= 0.0048) and 100% force effort (p= 0.0026).Significance.This study represents the first effort to non-invasively characterize the MU distribution inside spastic and non-spastic bilateral BBM of chronic stroke patients by combining HD-sEMG recording, EMG signal decomposition and bioelectrical source imaging. The findings of this study advances our understanding regarding the neurophysiology of human muscles and the neuromuscular alterations following stroke. It may also offer important MU depth information for botulinum toxin injection in clinical post-stroke spasticity management.

Resistance exercise training and the motor unit

Resistance exercise training (RET) is a key modality to enhance sports performance, injury prevention and rehabilitation, and improving overall health via increases in muscular strength. Yet, the contribution of neural mechanisms to increases in muscular strength are highly debated. This is particularly true for the involvement of the motor unit, which is the link between neural (activation) and mechanical (muscle fiber twitch forces) mechanisms. A plethora of literature that examines the effects of RET on skeletal muscle speculate the role of motor units, such as increases in firing rates partially explains muscular strength gains. Results, however, are mixed regarding changes in firing rates in studies that utilize single motor unit recordings. The lack of clarity could be related to vast or subtle differences in RET programs, methods to record motor units, muscles tested, types of contractions and intensities used to record motor units, etc. Yet to be discussed, mixed findings could be the result of non-uniform MU behavior that is not typically accounted for in RET research. The purpose of this narration is to discuss the effects of acute resistance exercise training studies on MU behavior and to provide guidance for future research.

The Time-Course Changes in Knee Flexion Range of Motion, Muscle Strength, and Rate of Force Development After Static Stretching

Previous studies have shown that longer-duration static stretching (SS) interventions can cause a decrease in muscle strength, especially explosive muscle strength. Furthermore, force steadiness is an important aspect of muscle force control, which should also be considered. However, the time course of the changes in these variables after an SS intervention remains unclear. Nevertheless, this information is essential for athletes and coaches to establish optimal warm-up routines. The aim of this study was to investigate the time course of changes in knee flexion range of motion (ROM), maximal voluntary isometric contraction (MVIC), rate of force development (RFD), and force steadiness (at 5 and 20% of MVIC) after three 60-s SS interventions. Study participants were sedentary healthy adult volunteers (n = 20) who performed three 60-s SS interventions of the knee extensors, where these variables were measured before and after SS intervention at three different periods, i.e., immediately after, 10 min, and 20 min the SS intervention (crossover design). The results showed an increase in ROM at all time points (d = 0.86-1.01). MVIC was decreased immediately after the SS intervention (d = -0.30), but MVIC showed a recovery trend for both 10 min (d = -0.17) and 20 min (d = -0.20) after the SS intervention. However, there were significant impairments in RFD at 100 m (p = 0.014, F = 6.37, η_p ² = 0.101) and 200 m (p < 0.01, F = 28.0, η_p ² = 0.33) up to 20 min after the SS intervention. Similarly, there were significant impairments in force steadiness of 5% (p < 0.01, F = 16.2, η_p ² = 0.221) and 20% MVIC (p < 0.01, F = 16.0, η_p ² = 0.219) at 20 min after the SS intervention. Therefore, it is concluded that three 60-s SS interventions could increase knee flexion ROM but impair explosive muscle strength and muscle control function until 20 min after the SS intervention.

The Effects of the Biceps Brachii and Brachioradialis on Elbow Flexor Muscle Strength and Spasticity in Stroke Patients

Objective

Muscle weakness and spasticity are common consequences of stroke, leading to a decrease in physical activity. The effective implementation of precision rehabilitation requires detailed rehabilitation evaluation. We aimed to analyze the surface electromyography (sEMG) signal features of elbow flexor muscle (biceps brachii and brachioradialis) spasticity in maximum voluntary isometric contraction (MVIC) and fast passive extension (FPE) in stroke patients and to explore the main muscle groups that affect the active movement and spasticity of the elbow flexor muscles to provide an objective reference for optimizing stroke rehabilitation.

Methods

Fifteen patients with elbow flexor spasticity after stroke were enrolled in this study. sEMG signals of the paretic and nonparetic elbow flexor muscles (biceps and brachioradialis) were detected during MVIC and FPE, and root mean square (RMS) values were calculated. The RMS values (mean and peak) of the biceps and brachioradialis were compared between the paretic and nonparetic sides. Additionally, the correlation between the manual muscle test (MMT) score and the RMS values (mean and peak) of the paretic elbow flexors during MVIC was analyzed, and the correlation between the modified Ashworth scale (MAS) score and the RMS values (mean and peak) of the paretic elbow flexors during FPE was analyzed.

Results

During MVIC exercise, the RMS values (mean and peak) of the biceps and brachioradialis on the paretic side were significantly lower than those on the nonparetic side (p < 0.01), and the RMS values (mean and peak) of the bilateral biceps were significantly higher than those of the brachioradialis (p < 0.01). The MMT score was positively correlated with the mean and peak RMS values of the paretic biceps and brachioradialis (r = 0.89, r = 0.91, r = 0.82, r = 0.85; p < 0.001). During FPE exercise, the RMS values (mean and peak) of the biceps and brachioradialis on the paretic side were significantly higher than those on the nonparetic side (p < 0.01), and the RMS values (mean and peak) of the brachioradialis on the paretic side were significantly higher than those of the biceps (p < 0.01). TheMAS score was positively correlated with the mean RMS of the paretic biceps and brachioradialis (r = 0.62, p = 0.021; r = 0.74, p = 0.004), and the MAS score was positively correlated with the peak RMS of the paretic brachioradialis (r = 0.59, p = 0.029) but had no significant correlation with the peak RMS of the paretic biceps (r = 0.49, p > 0.05).

Conclusions

The results confirm that the biceps is a vital muscle in active elbow flexion and that the brachioradialis plays an important role in elbow flexor spasticity, suggesting that the biceps should be the focus of muscle strength training of the elbow flexors and that the role of the brachioradialis should not be ignored in the treatment of elbow flexor spasticity. This study also confirmed the application value of sEMG in the objective assessment of individual muscle strength and spasticity in stroke patients.

Method of analysis influences interpretations of sex-related differences in firing rates during prolonged submaximal isometric contractions

OBJECTIVE

This study examined motor unit (MU) firing rates during a prolonged isometric contraction of the vastus lateralis (VL) for females and males.

METHODS

Surface electromyographic (sEMG) signals were recorded from the VL for eleven females and twelve males during a 45-second isometric trapezoid muscle actions at 40% of maximal voluntary contraction (MVC). For each MU, mean firing rate (MFR) was calculated for the initial and final 10-second epochs of the steady torque segment and regressed against recruitment threshold (RT, expressed as %MVC), as well as time at recruitment (T_REC, seconds). MFR was also averaged for each subject.

RESULTS

Significant differences existed across epochs for the y-intercepts (P=0.009) of the MFR vs. T_REC relationship, as well as the grouped MFR analysis (P<0.001); no differences were observed between epochs for the MFR vs. RT relationship. Significant differences existed between sexes for the grouped MFR analysis (P=0.049), but no differences were observed for the MFR vs. T_REC or MFR vs. RT relationships.

CONCLUSION

Analysis method may impact interpretation of firing rate behavior; increases in MU firing rates across a prolonged isometric contraction were observed in the MFR vs. T_REC relationship and the grouped MFR analysis.

Decoding Muscle Force From Motor Unit Firings Using Encoder-Decoder Networks

Appropriate interpretation of motor unit (MU) activities after surface EMG (sEMG) decomposition is a key factor to decode motor intentions in a noninvasive and physiologically meaningful way. However, there are great challenges due to the difficulty in cross-trial MU tracking and unavoidable loss of partial MU information resulting from incomplete decomposition. In light of these issues, this study presents a novel framework for interpreting MU activities and applies it to decode muscle force. The resulting MUs were clustered and classified into different categories by characterizing their spatially distributed firing waveforms. The process served as a general MU tracking method. On this basis, after transferring the MU firing trains to twitch force trains by a twitch force model, a deep network was designed to predict the normalized force. In addition, MU category distribution was examined to calibrate the actual force level, while functions of some unavailable MUs were compensated. To investigate the effectiveness of this framework, high-density sEMG signals were recorded using an 8×8 electrode array from the abductor pollicis brevis muscles of eight subjects, while thumb abduction force was measured. The proposed method outperformed three common methods ( ) yielding the lowest root mean square deviation of 6.68% ± 1.29% and the highest fitness ( R² ) of 0.94 ± 0.04 between the predicted force and the actual force. This study offers a valuable, computational solution for interpreting individual MU activities, and its effectiveness was confirmed in muscle force estimation.

An examination of motor unit firing rates during steady torque of maximal efforts with either an explosive or slower rate of torque development

NEW FINDINGS

What is the central question of this study? The aim was to explore agonist and antagonist motor unit firing rates during maximal efforts performed with either an explosive or a slower rate of torque development. What is the main finding and its importance? The antagonist muscle presented a motor unit firing rate relationship similar to the agonist muscle. Additionally, the motor units of both muscles exhibited higher firing rates during explosive maximal contractions than during maximal contractions performed at a slower rate of torque development. These results could prove useful to future research analysing the effects of age, disease, resistance training and/or fatigue-related alterations to motor unit firing rates.

ABSTRACT

The primary purpose of the present study was to examine motor unit (MU) firing rates in agonist and antagonist muscles during periods of steady, maximal efforts using explosive and slower rates of torque development. A secondary purpose was to analyse the MU firing rate versus action potential amplitude relationships of the agonist and antagonist muscles during maximal efforts. Thirteen subjects (mean ± SD; age, 21.2 ± 3.6 years; mass 81.1 ± 21.3 kg; and stature, 177.1±9.9 cm) performed two maximal isometric trapezoid muscle actions of the elbow flexors that included either an explosive or a slower, linearly increasing rate (ramp) of torque development. Surface EMG signals of the biceps brachii (BB) and triceps brachii (TB) muscles were collected and decomposed into their constituent MU action potential trains. The MU firing rate versus action potential amplitude relationships of the BB (agonist) and TB (antagonist) muscles were analysed. Moderate to strong relationships (|r| ≥ 0.65) were present for the explosive and ramp contractions in the agonist and antagonist muscles. Firing rates of smaller and larger MUs were higher during the explosive [mean ± SD; agonist = 18.1 ± 6.9 pulses per second (pps), antagonist = 22.0±3.9 pps] than the ramp (agonist = 14.0 ± 5.1 pps, antagonist = 18.3 ± 4.4 pps) contractions for the agonist (P = 0.013) and antagonist muscles (P = 0.007). The antagonist muscle exhibits a similar MU firing rate versus action potential amplitude relationship to the agonist muscle at maximal efforts. Future research should investigate the effects of short-term resistance training on antagonist firing rates and the involvement of peripheral feedback on firing rates during maximal efforts performed at various rates of torque development.

The reliability of the slopes and y-intercepts of the motor unit firing times and action potential waveforms versus recruitment threshold relationships derived from surface electromyography signal decomposition

PURPOSE

Changes in motor unit (MU) activity pre- to post-interventions is of high interest. However, there is minimal information regarding day-to-day changes or the reliability of measuring MU activity. This study examined the reliability of relationships calculated via the MU action potential (AP) trains derived from surface electromyography signal decomposition. A comparison between reliability statistics was made between MUAP trains verified with only the reconstruct-and-test versus verification including reconstruct-and-test with spike trigger average (STA) procedures.

METHODS

Twenty-one individuals performed isometric muscle actions at 50% of maximal voluntary contraction of the first dorsal interosseous on two separate visits. MUs included for reliability analyses initially met the > 90% accuracy from the reconstruct-and-test. STA was applied as an additional exclusionary procedure. Linear regressions were applied to the firing rate and AP amplitude versus recruitment threshold relationships with and without MUs that met the STA criteria. Reliability statistics were also performed on relationships that met a strict range of recruitment thresholds. Reliability was established with intraclass correlation coefficients (ICC) along with other traditional parameters.

RESULTS

The firing rate versus recruitment threshold relationships were reliable (ICC > 0.56) and improved (ICC > 0.84) when recruitment ranges were controlled. The slopes of the MUAP amplitude versus recruitment threshold relationships were reliable (ICC > 0.78) while the y-intercepts were reliable (ICC > 0.81) once corrections were made to combat negative scores.

CONCLUSIONS

Electromyographic signal decomposition without the secondary STA verification procedures can be used to detect intervention-related changes in neural drive with confidence when recorded MU recruitment thresholds are similar across days.

Adaptive Real-time Decomposition of Electromyogram During Sustained Muscle Activation: A Simulation Study

OBJECTIVE

Real-time decomposition of electromyogram (EMG) into constituent motor unit (MU) activity has shown promising applications in neurophysiology and human-machine interactions. Existing decomposition methods could not accommodate stochastic variations in EMG signals such as drifts of action potential amplitudes and MU recruitment-derecruitment (rotation) patterns during long-term recordings. The objective of this study was to develop an adaptive real-time decomposition approach suitable for prolonged muscle activation.

METHODS

We developed a parallel-double-thread computation algorithm. The backend thread initiated and periodically refined and updated the MU information (separation matrix) using independent component analysis and convolution kernel compensation. The frontend thread performed the real-time decomposition. We evaluated our algorithm on synthesized high-density EMG signals, in which MUs were recruited-derecruited sporadically and MU action potentials amplitude drifted over time. Different signal-to-noise levels were also simulated.

RESULTS

Compared with the decomposition without the adaptive processes, periodically fine-tuned and updated separation matrix increased identifiable MU number by 3-4 fold over 30-minute of recordings. The increased MU number was more prominent at higher signal-to-noise ratios. The decomposition accuracy also increased by up to 10% with greater improvement observed at higher muscle contraction levels.

CONCLUSION

The adaptive algorithm can maintain the decomposition performance over time, allows us to continuously track the same MUs during sustained activation, and, at the same time, can add newly recruited MU information to existing separation matrix.

SIGNIFICANCE

Our approach showed robust performance over time, which has the potential to longitudinally evaluate MU firing and recruitment properties and improve neural decoding performance for neural-machine interactions.

An examination of a potential organized motor unit firing rate and recruitment scheme of an antagonist muscle during isometric contractions

The primary purpose of the present study is to determine if an organized control scheme exists for the antagonist muscle during steady isometric torque. A secondary focus is to better understand how firing rates of the antagonist muscle change from a moderate- to higher-contraction intensity. Fourteen subjects performed two submaximal isometric trapezoid muscle actions of the forearm flexors that included a linearly increasing, steady force at both 40% and 70% maximum voluntary contraction, and linearly decreasing segments. Surface electromyographic signals of the biceps and triceps brachii were collected and decomposed into constituent motor unit action potential trains. Motor unit firing rate versus recruitment threshold, motor unit action potential amplitude versus recruitment threshold, and motor unit firing rate versus action potential amplitude relationships of the biceps brachii (agonist) and triceps brachii (antagonist) muscles were analyzed. Moderate- to-strong relationships (|r| ≥ 0.69) were present for the agonist and antagonist muscles for each relationship with no differences between muscles (P = 0.716, 0.428, 0.182). The y-intercepts of the motor unit firing rate versus recruitment threshold relationship of the antagonist did not increase from 40% to 70% maximal voluntary contractions (P = 0.96), unlike for the agonist (P = 0.009). The antagonist muscle exhibits a similar motor unit control scheme to the agonist. Unlike the agonist, however, the firing rates of the antagonist did not increase with increasing intensity. Future research should investigate how antagonist firing rates adapt to resistance training and changes in antagonist firing rates in the absence of peripheral feedback.NEW & NOTEWORTHY This is the first study to explore a potential motor unit control scheme and quantify changes in firing rates with increasing intensity of an antagonist muscle during isometric contractions. We demonstrate that the antagonist muscle possesses an organized motor unit firing rate and recruitment scheme similar to the agonist muscle during isometric forearm flexion, but unlike the agonist muscle, there was no significant increase in firing rates from a moderate- to higher-intensity isometric contraction.

Concurrent Prediction of Finger Forces Based on Source Separation and Classification of Neuron Discharge Information

A reliable neural-machine interface is essential for humans to intuitively interact with advanced robotic hands in an unconstrained environment. Existing neural decoding approaches utilize either discrete hand gesture-based pattern recognition or continuous force decoding with one finger at a time. We developed a neural decoding technique that allowed continuous and concurrent prediction of forces of different fingers based on spinal motoneuron firing information. High-density skin-surface electromyogram (HD-EMG) signals of finger extensor muscle were recorded, while human participants produced isometric flexion forces in a dexterous manner (i.e. produced varying forces using either a single finger or multiple fingers concurrently). Motoneuron firing information was extracted from the EMG signals using a blind source separation technique, and each identified neuron was further classified to be associated with a given finger. The forces of individual fingers were then predicted concurrently by utilizing the corresponding motoneuron pool firing frequency of individual fingers. Compared with conventional approaches, our technique led to better prediction performances, i.e. a higher correlation ([Formula: see text] versus [Formula: see text]), a lower prediction error ([Formula: see text]% MVC versus [Formula: see text]% MVC), and a higher accuracy in finger state (rest/active) prediction ([Formula: see text]% versus [Formula: see text]%). Our decoding method demonstrated the possibility of classifying motoneurons for different fingers, which significantly alleviated the cross-talk issue of EMG recordings from neighboring hand muscles, and allowed the decoding of finger forces individually and concurrently. The outcomes offered a robust neural-machine interface that could allow users to intuitively control robotic hands in a dexterous manner.

Motor unit action potential amplitude during low torque fatiguing contractions versus high torque non-fatiguing contractions: a multilevel analysis

PURPOSE

The ability to maintain an absolute, submaximal torque level during fatiguing contractions is controlled, in part, by the recruitment of larger motor units. These motor units are commonly identified based on greater action potential peak-to-peak amplitude values. It is unclear, however, if motor unit action potential (MUAP) amplitude values during low torque, fatiguing contractions reach similar levels as those observed during non-fatigued, high torque contractions. To establish a clearer understanding of motor unit control during fatigue, we compared MUAP amplitude during 50 and 80% maximum voluntary contraction (MVC) torque contractions and at the beginning, middle, and end of a 30% MVC fatigue protocol.

METHODS

Eleven untrained men (mean age = 24 years) performed isometric contractions at 50 and 80% MVC, followed by repeated contractions at 30% MVC. Surface electromyographic (EMG) signals were detected from the vastus lateralis and decomposed to quantify the peak-to-peak amplitude of individual MUAPs. A two-level multilevel model was estimated, allowing examination of simultaneous measures of MUAP amplitude within participants and controlling for the dependence between measures within participants.

RESULTS

Results from the multilevel analyses suggested that there were not statistically significant differences in MUAP amplitude between 80% MVC and end fatigue. Separate repeated-measures analyses of variance indicated that there were not statistically significant mean differences in greatest MUAP or surface EMG amplitude between 80% MVC and end fatigue.

CONCLUSIONS

MUAP and surface EMG amplitude values during a 30% MVC fatiguing protocol appear to be comparable to those observed during a non-fatigued 80% MVC condition.

Improving Botulinum Toxin Efficiency in Treating Post-Stroke Spasticity Using 3D Innervation Zone Imaging

Spasticity is a common post-stroke syndrome that imposes significant adverse impacts on patients and caregivers. This study aims to improve the efficiency of botulinum toxin (BoNT) in managing spasticity, by utilizing a three-dimensional innervation zone imaging (3DIZI) technique based on high-density surface electromyography (HD-sEMG) recordings. Stroke subjects were randomly assigned to two groups: the control group ([Formula: see text]) which received standard ultrasound-guided injections, and the experimental group ([Formula: see text]) which received 3DIZI-guided injections. The amount of BoNT given was consistent for all subjects. The Modified Ashworth Scale (MAS), compound muscle action potential (CMAP) and muscle activation volume (MAV) from bilateral biceps brachii muscles were obtained at the baseline, 3 weeks, and 3 months after injection. Intra-group and inter-group comparisons of MAS, CMAP amplitude and MAV were performed. An overall improvement in MAS of spastic elbow flexors was observed during the 3-week visit ([Formula: see text]), yet no statistically significant difference found with intra-group or inter-group analysis. Compared to the baseline, a significant reduction of CMAP amplitude and MAV were observed in the spastic biceps muscles of both groups at 3-week post-injection, and returned to approximate baseline value at 12-week post injection. A significantly higher reduction was found in CMAP amplitude ([Formula: see text]% versus [Formula: see text]%, [Formula: see text]) and MAV ([Formula: see text]% versus [Formula: see text]%, [Formula: see text]) in the experimental group compared to the control group. The study has demonstrated preliminary evidence that precisely directing BoNT to the innervation zones (IZs) localized by 3DIZI leads to a significantly higher treatment efficiency improvement in spasticity management. Results have also shown the feasibility of developing a personalized BoNT injection technique for the optimization of clinical treatment for post-stroke spasticity using proposed 3DIZI technique.

Neural Drive is Greater for a High-Intensity Contraction Than for Moderate-Intensity Contractions Performed to Fatigue

Miller, JD, Lippman, JD, Trevino, MA, and Herda, TJ. Larger motor units are recruited for high-intensity contractions than for fatiguing moderate-intensity contractions. J Strength Cond Res 34(11): 3013-3021, 2020-The purpose of this study was to investigate whether moderate-intensity contractions performed to fatigue activate the motor unit (MU) pool to the same extent as a higher-intensity contraction. Subjects (7 men, 2 women, age = 22.78 ± 4.15 years, height = 173.78 ± 14.19 cm, mass = 87.39 ± 21.19 kg) performed 3 isometric maximum voluntary contractions (MVCs), an isometric trapezoidal contraction at 90% MVC (REP90), and repetitive isometric trapezoidal contractions at 50% MVC performed to failure with the first (REP1) and final repetition (REPL) used for analysis. Surface EMG was recorded from the vastus lateralis. Action potentials were extracted into firing events of single MUs with recruitment thresholds (RTs), MU action potential amplitudes (MUAPAMP), and mean firing rates (MFRs) recorded. Linear MFR and MUAPAMP vs. RT and exponential MFR vs. MUAPAMP relationships were calculated for each subject. The level of significance was set at p ≤ 0.05. B terms for the MFR vs. MUAPAMP relationships (p = 0.001, REPL = -4.77 ± 1.82 pps·mV, REP90 = -2.63 ± 1.00 pps·mV) and predicted MFRs for MUs recruited at 40% MVC (p < 0.001, REPL = 11.14 ± 3.48 pps, REP90 = 18.38 ± 2.60 pps) were greater for REP90 than REPL indicating firing rates were greater during REP90. In addition, larger mean (p = 0.038, REPL = 0.178 ± 0.0668 mV, REP90 = 0.263 ± 0.128 mV) and maximum (p = 0.008, REPL = 0.320 ± 0.127 mV, Rep90 = 0.520 ± 0.234 mV) MUAPAMPS were recorded during REP90 than REPL. Larger MUs were recruited and similar sized MUs maintained greater firing rates during a high-intensity contraction in comparison to a moderate-intensity contraction performed at fatigue. Individuals seeking maximized activation of the MU pool should use high-intensity resistance training paradigms rather than moderate-intensity to fatigue.

Does strict validation criteria for individual motor units alter population-based regression models of the motor unit pool?

The purpose of this study was to determine if the implementation of a strict validation procedure, designed to limit the inclusion of inaccuracies from the decomposition of surface electromyographic (sEMG) signals, affects population-based motor unit (MU) analyses. Four sEMG signals were obtained from the vastus lateralis of 59 participants during isometric contractions at different relative intensities [30%, 70%, and 100% of maximal voluntary contraction (MVC)], and its individual motor unit potential trains (MUPTs) were extracted. The MUPTs were then excluded (ISIval) based on the coefficient of variation and histogram of the interspike intervals (ISI), the absence of additional clusters that reveals missed or additional firings, and more. MU population-based regression models (i.e., modeling the entire motor unit pool) were performed between motor unit potential size (MUP_SIZE), mean firing rate (MFR), and recruitment threshold (RT%) separately for DSDC_Only (includes all MUPTs without the additional validation performed) and ISIval data at each contraction intensity. The only significant difference in regression coefficients between DSDC_Only and ISIval was for the intercepts of the MUP_SIZE/MFR at 100% MVC. The validation had no other significant effect on any of the other regression coefficients for each of the contraction intensities. Our findings suggest that even though the decomposition of surface signals leads to some inaccuracies, these errors have limited effects on the regression models used to estimate the behavior of the whole pool. Therefore, we propose that motor unit population-based regression models may be robust enough to overcome decomposition-induced errors at the individual MU level.

Differences in the firing rate versus recruitment threshold relationships of the vastus lateralis in children ages 7-10 years and adults

Motor unit (MU) firing rates of the vastus lateralis in children and adults were examined. Seven healthy adult males (mean ± SD, age = 21 ± 2.6 yrs) and six healthy male children (mean ± SD age = 8.8 ± 1.7) volunteered. Surface electromyography (EMG) signals were recorded from 20% and 60% maximal voluntary contractions (MVC). Surface EMG signals were decomposed into firing events of individual MUs and slopes and y-intercepts were calculated for the mean firing rate (MFR, pps) at steady torque vs. recruitment thresholds (RT) relationships for each subject. Muscle cross-sectional area (mCSA) was measured, via ultrasonography, with specific torque calculated (MVC/mCSA). Adults possessed greater mCSA (p = .002; children = 11.5 ± 2.1 cm²; adults = 31.80 ± 12.15 cm²) and greater specific torque (p = .018; children = 4.63 ± 1.4 Nm/cm²; adults = 7.1 ± 1.8 Nm/cm²) compared to children. The y-intercepts were significantly (p < .001) greater during the 60% (28.91 ± 4.56 pps) than the 20% MVC (23.5 ± 4.9 pps) collapsed across groups while the children had significantly (p = .036) lower y-intercepts (23.9 ± 5.4 pps) than the adults (28.2 ± 4.8 pps) when collapsed across intensities. Slopes of the MFR vs RT relationships were greater for the 60% (-0.342 ± 0.127 pps/%MVC) contraction than the 20% (-0.50 ± 0.159 pps/%MVC) MVC when collapsed across groups. Adults had greater firing rates regardless of recruitment threshold than children. This may be due to lower recruitment potential and overall excitation to the motoneuron pool of children as indicated with differences in specific torque and/or differences in antagonist co-activation.

Measuring the accuracies of motor unit firing times and action potential waveforms derived from surface electromyographic decomposition

This study included spike trigger averaging (STA) procedures to examine the acceptability of the Precision Decomposition (PD) III derived motor unit action potential (MUAP) trains that met the >90% accuracy criteria from the reconstruct-and-test. MUs met the >90% accuracy criteria from the reconstruct-and-test with STA procedures then applied. Y-intercepts and slopes were calculated for the firing rate- and MUAP amplitude-recruitment threshold relationships. Gaussian noise (1% of the SD of the mean interspike interval) was added to the firing times with the changes in MUAPs quantified. A total of 455 MUs were decomposed with 155 MUs removed as a result of the reconstruct-and-test. Five additional MUs were excluded via the STA criteria. The MUAP waveforms deteriorated with the inclusion of Gaussian noise. There were differences in the derived action potentials amplitudes of higher-threshold MUs between the PD III algorithm and the STA procedure. There was excellent agreement among the slopes and y-intercepts between the relationships that included or excluded MUs that did not meet the STA criteria. There was good agreement between the MUAP amplitude-recruitment threshold relationships derived from the PD III and STA procedure. The addition of the STA procedures did not alter the MU-derived relationships.

Sex-related differences in motor unit firing rates and action potential amplitudes of the first dorsal interosseous during high-, but not low-intensity contractions

Despite ample evidence that females are weaker and possess smaller muscle cross-sectional areas (CSAs) compared to males, it remains unclear if there are sex-related differences in the properties of motor units (MU). Eleven males (age 22 ± 3 years) and 12 females (age 21 ± 1 years) performed isometric trapezoid muscle actions at 10% and 70% of maximal voluntary contraction (MVC). Surface electromyography signals were recorded and decomposed into MU action potential (AP) waveforms and firing instances. Average MUAP amplitudes (MUAP_AMPS), mean firing rates (MFRs), initial firing rates (IFRs), and recruitment thresholds (RT) were calculated for the 10% MVC, while MUAP_AMPS, IFRs, and MFRs were regressed against RT for the 70% MVC. Ultrasonography was used to measure CSA of the first dorsal interosseous (FDI). Males had greater CSAs (p < 0.001; males 2.34 ± 0.28 cm², females 1.82 ± 0.18 cm²) and MVC strength (p < 0.001; males 25.9 ± 5.5 N, females 16.44 ± 2.5 N). No differences existed for MUAP_AMPS, IFRs, MFRs, or RTs (p > 0.05) during the 10% MVC. For the 70% MVC, the y-intercepts from the MUAP_AMPS vs. RT relationships were greater (p < 0.05) for the males (males - 0.19 ± 0.53 mV; females - 0.78 ± 0.75 mV), while the inverse was true for the MFR vs. RT relationships (males 31.55 ± 6.92 pps, females 38.65 ± 6.71 pps) with no differences (p > 0.05) in the slopes. Therefore, smaller CSAs and weaker MVCs are likely the result of smaller higher-threshold MUs for females.

Coupling Robot-Aided Assessment and Surface Electromyography (sEMG) to Evaluate the Effect of Muscle Fatigue on Wrist Position Sense in the Flexion-Extension Plane

Proprioception is a crucial sensory modality involved in the control and regulation of coordinated movements and in motor learning. However, the extent to which proprioceptive acuity is influenced by local muscle fatigue is obscured by methodological differences in proprioceptive and fatiguing protocols. In this study, we used high resolution kinematic measurements provided by a robotic device, as well as both frequency and time domain analysis of signals captured via surface electromyography (sEMG) to examine the effects of local muscle fatigue on wrist proprioceptive acuity in 16 physically and neurologically healthy young adults. To this end, participants performed a flexion/extension ipsilateral joint position matching test (JPM), after which a high-resistive robotic task was used to induce muscle fatigue of the flexor carpi radialis (FCR) muscle. The JPM test was then repeated in order to analyze potential changes in proprioceptive acuity. Results indicated that the fatigue protocol had a significant effect on movements performed in flexion direction, with participants exhibiting a tendency to undershoot the target before the fatigue protocol (−1.218°), but overshooting after the fatigue protocol (0.587°). In contrast, in the extension direction error bias values were similar before and after the fatigue protocol as expected (pre = −1.852°, post = −1.237°) and reflected a tendency to undershoot the target. Moreover, statistical analysis indicated that movement variability was not influenced by the fatigue protocol or movement direction. In sum, results of the present study demonstrate that an individual’s estimation of wrist joint displacement (i.e., error bias), but not precision (i.e., variability), is affected by muscular fatigue in a sample of neurologically and physically healthy adults.

Motor unit firing rates of the first dorsal interosseous differ between male and female children aged 8-10 years

The purpose of this study was to examine possible differences in motor unit action potential amplitudes (MUAP_AMPS) and firing rates of the first dorsal interosseous (FDI) in male and female children aged 8-10 years. Eight male (mean ± SD, age = 8.8 ± 0.7 yrs; BMI = 16.5 ± 1.3 kg/m²) and eight female (age = 9.3 ± 0.9 yrs; BMI = 16.1 ± 1.5 kg/m²) children volunteered to complete isometric trapezoidal muscle actions of the first dorsal interosseous at 50% of maximal voluntary contraction (MVC). Electromyographic signals were decomposed to yield MUAP_AMPS and mean firing rates (MFR) at the targeted force. An exponential model was fitted to the MUAP_AMPS vs. recruitment threshold (RT) while linear models were fitted to the MFRs vs. RT relationships for each subject. Ultrasonography determined the muscle cross-sectional area (CSA) of the FDI. Independent samples t-tests were used to examine possible differences between the male and female children for MVC strength, CSA, and the coefficients from the MU relationships. There were no differences in MVC strength, CSA, or the MUAP_AMP vs. RT relationships between the male and female children (P < 0.05). Males, however, had greater MFRs of lower-threshold MUs as evident by significantly larger y-intercepts (P = 0.019) and more negative slopes (P = 0.004) from the MFR vs. RT relationships. Despite no differences in muscle strength, CSA, and MUAP_AMPS, differences in firing rates existed between male and female children aged 8-10 years. Neural mechanisms may primarily contribute to sex-related differences in firing rates.