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

Persistent inward currents in tibialis anterior motoneurons can be reliably estimated within the same session

The response of spinal motoneurons to synaptic input greatly depends on the activation of persistent inward currents (PICs), the contribution of which can be estimated through the paired motor unit technique. Yet, the intra-session test-retest reliability of this measurement remains to be fully established. Twenty males performed isometric triangular dorsiflexion contractions to 20 and 50 % of maximal torque at baseline and after a 15-min resting period. High-density electromyographic signals (HD-EMG) of the tibialis anterior were recorded with a 64-electrode matrix. HD-EMG signals were decomposed, and motor units tracked across time points to estimate the contribution of PICs to motoneuron firing through quantification of motor unit recruitment-derecruitment hysteresis (ΔF). A good intraclass correlation coefficient (ICC = 0.75 [0.63, 0.83]) and a large repeated measures correlation coefficient (rrm = 0.65 [0.49, 0.77]; p < 0.001) were found between ΔF values obtained at both time points for 20 % MVC ramps. For 50 % MVC ramps, a good ICC (0.77 [0.65, 0.85]) and a very large repeated measures correlation coefficient (rrm = 0.73 [0.63, 0.80]; p < 0.001) were observed. Our data suggest that ΔF scores can be reliably investigated in tibialis anterior motor units during both low- and moderate-intensity contractions within a single experimental session.

Changes in motor unit behaviour across repeated bouts of eccentric exercise

Unaccustomed eccentric exercise (EE) is protective against muscle damage following a subsequent bout of similar exercise. One hypothesis suggests the existence of an alteration in motor unit (MU) behaviour during the second bout, which might contribute to the adaptive response. Accordingly, the present study investigated MU changes during repeated bouts of EE. During two bouts of exercise where maximal lengthening dorsiflexion (10 repetitions × 10 sets) was performed 3 weeks apart, maximal voluntary isometric torque (MVIC) and MU behaviour (quantified using high-density electromyography; HDsEMG) were measured at baseline, during (after set 5), and post-EE. The HDsEMG signals were decomposed into individual MU discharge timings, and a subset were tracked across each time point. MVIC was reduced similarly in both bouts post-EE (Δ27 vs. 23%, P = 0.144), with a comparable amount of total work performed (∼1,300 J; P = 0.905). In total, 1,754 MUs were identified and the decline in MVIC was accompanied by a stepwise increase in discharge rate (∼13%; P < 0.001). A decrease in relative recruitment was found immediately after EE in Bout 1 versus baseline (∼16%; P < 0.01), along with reductions in derecruitment thresholds immediately after EE in Bout 2. The coefficient of variation of inter-spike intervals was lower in Bout 2 (∼15%; P < 0.001). Our data provide new information regarding a change in MU behaviour during the performance of a repeated bout of EE. Importantly, such changes in MU behaviour might contribute, at least in part, to the repeated bout phenomenon.

Can non-invasive motor unit analysis reveal distinct neural strategies of force production in young with uncomplicated type 1 diabetes?

PURPOSE

to investigate the early consequences of type 1 diabetes (T1D) on the neural strategies of muscle force production.

METHODS

motor unit (MU) activity was recorded from the vastus lateralis muscle with High-Density surface Electromyography during isometric knee extension at 20 and 40% of maximum voluntary contraction (MVC) in 8 T1D (4 males, 4 females, 30.5 ± 3.6 years) and 8 matched control (4 males, 4 females, 27.3 ± 5.9 years) participants. Muscle biopsies were also collected from vastus lateralis for fiber type analysis, including myosin heavy chain (MyHC) isoform content via protein and mRNA expression.

RESULTS

MVC was comparable between groups as well as MU conduction velocity, action potentials' amplitude and proportions of MyHC protein isoforms. Nonetheless, MU discharge rate, relative derecruitment thresholds and mRNA expression of MyHC isoform I were lower in T1D.

CONCLUSIONS

young people with uncomplicated T1D present a different neural control of muscle force production. Furthermore, differences are detectable non-invasively in absence of any functional manifestation (i.e., force production and fiber type distribution). These novel findings suggest that T1D has early consequences on the neuromuscular system and highlights the necessity of a better characterization of neural control in this population.

Acute neuromuscular and perceptual responses to blood flow restriction exercise in adults with severe haemophilia: A pilot study

INTRODUCTION

No previous studies have implemented a standard blood flow restriction (BFR) training session in people with severe haemophilia (PwH), where this type of training has been contraindicated.

AIMS

The purpose of this study was to evaluate the tolerability, adverse events, and neuromuscular and perceptual responses to an acute session of low load (LL) knee extensions with BFR in PwH under prophylaxis.

METHODS

Eight PwH performed one LL-BFR session with 40% arterial occlusion pressure (AOP). Perceptual responses and adverse effects were assessed, together with high-density surface electromyography of vastus medialis (VM) and lateralis (VL).

RESULTS

Significant normalized root mean square differences were found within each set, but not between sets. Spatial distribution (centroid displacement (p > .05), modified entropy (VM, set two, cycles three and five, p = .032) and coefficient of variation (VM, set two, cycles four and five lower than cycle three (p = .049; p = .036)) showed changes within each set. Median frequency showed a slight increase during cycle four of set four (p = .030). Rate of perceived exertion slightly increased with each set while tolerability slightly decreased in the last set and fear of training with BFR generally decreased after the session.

CONCLUSIONS

In PwH, a LL-BFR session at 40% AOP is safe and feasible. Our results suggest that potential muscle impairments may blunt neuromuscular adaptations induced by BFR.

Adaptive EMG decomposition in dynamic conditions based on online learning metrics with tunable hyperparameters

Objective. Developing neural decoders robust to non-stationary conditions is essential to ensure their long-term accuracy and stability. This is particularly important when decoding the neural drive to muscles during dynamic contractions, which pose significant challenges for stationary decoders.Approach. We propose a novel adaptive electromyography (EMG) decomposition algorithm that builds on blind source separation methods by leveraging the Kullback-Leibler divergence and kurtosis of the signals as metrics for online learning. The proposed approach provides a theoretical framework to tune the adaptation hyperparameters and compensate for non-stationarities in the mixing matrix, such as due to dynamic contractions, and to identify the underlying motor neuron (MN) discharges. The adaptation is performed in real-time (∼22 ms of computational time per 100 ms batches).Main results. The hyperparameters of the proposed adaptation captured anatomical differences between recording locations (forearm vs wrist) and generalised across subjects. Once optimised, the proposed adaptation algorithm significantly improved all decomposition performance metrics with respect to the absence of adaptation in a wide range of motion of the wrist (80∘). The rate of agreement, sensitivity, and precision were⩾90%in⩾80%of the cases in both simulated and experimentally recorded data, according to a two-source validation approach.Significance. The findings demonstrate the suitability of the proposed online learning metrics and hyperparameter optimisation to compensate the induced modulation and accurately decode MN discharges in dynamic conditions. Moreover, the study proposes an experimental validation method for EMG decomposition in dynamic tasks.

Neural Filtering of Physiological Tremor Oscillations to Spinal Motor Neurons Mediates Short-Term Acquisition of a Skill Learning Task

The acquisition of a motor skill involves adaptations of spinal and supraspinal pathways to alpha motoneurons. In this study, we estimated the shared synaptic contributions of these pathways to understand the neural mechanisms underlying the short-term acquisition of a new force-matching task. High-density surface electromyography (HDsEMG) was acquired from the first dorsal interosseous (FDI; 7 males and 6 females) and tibialis anterior (TA; 7 males and 4 females) during 15 trials of an isometric force-matching task. For two selected trials (pre- and post-skill acquisition), we decomposed the HDsEMG into motor unit spike trains, tracked motor units between trials, and calculated the mean discharge rate and the coefficient of variation of interspike interval (COVISI). We also quantified the post/pre ratio of motor units' coherence within delta, alpha, and beta bands. Force-matching improvements were accompanied by increased mean discharge rate and decreased COVISI for both muscles. Moreover, the area under the curve within alpha band decreased by ∼22% (TA) and ∼13% (FDI), with no delta or beta bands changes. These reductions correlated significantly with increased coupling between force/neural drive and target oscillations. These results suggest that short-term force-matching skill acquisition is mediated by attenuation of physiological tremor oscillations in the shared synaptic inputs. Supported by simulations, a plausible mechanism for alpha band reductions may involve spinal interneuron phase-cancelling descending oscillations. Therefore, during skill learning, the central nervous system acts as a matched filter, adjusting synaptic weights of shared inputs to suppress neural components unrelated to the specific task.

Muscle contractile properties directly influence shared synaptic inputs to spinal motor neurons

Alpha band oscillations in shared synaptic inputs to the alpha motor neuron pool can be considered an involuntary source of noise that hinders precise voluntary force production. This study investigated the impact of changing muscle length on the shared synaptic oscillations to spinal motor neurons, particularly in the physiological tremor band. Fourteen healthy individuals performed low-level dorsiflexion contractions at ankle joint angles of 90° and 130°, while high-density surface electromyography (HDsEMG) was recorded from the tibialis anterior (TA). We decomposed the HDsEMG into motor units spike trains and calculated the motor units' coherence within the delta (1-5 Hz), alpha (5-15 Hz), and beta (15-35 Hz) bands. Additionally, force steadiness and force spectral power within the tremor band were quantified. Results showed no significant differences in force steadiness between 90° and 130°. In contrast, alpha band oscillations in both synaptic inputs and force output decreased as the length of the TA was moved from shorter (90°) to longer (130°), with no changes in delta and beta bands. In a second set of experiments (10 participants), evoked twitches were recorded with the ankle joint at 90° and 130°, revealing longer twitch durations in the longer TA muscle length condition compared to the shorter. These experimental results, supported by a simple computational simulation, suggest that increasing muscle length enhances the muscle's low-pass filtering properties, influencing the oscillations generated by the Ia afferent feedback loop. Therefore, this study provides valuable insights into the interplay between muscle biomechanics and neural oscillations. KEY POINTS: We investigated whether changes in muscle length, achieved by changing joint position, could influence common synaptic oscillations to spinal motor neurons, particularly in the tremor band (5-15 Hz). Our results demonstrate that changing muscle length from shorter to longer induces reductions in the magnitude of alpha band oscillations in common synaptic inputs. Importantly, these reductions were reflected in the oscillations of muscle force output within the alpha band. Longer twitch durations were observed in the longer muscle length condition compared to the shorter, suggesting that increasing muscle length enhances the muscle's low-pass filtering properties. Changes in the peripheral contractile properties of motor units due to changes in muscle length significantly influence the transmission of shared synaptic inputs into muscle force output. These findings prove the interplay between muscle mechanics and neural adaptations.

More Variability in Tibialis Anterior Function during the Adduction of the Foot than Dorsiflexion of the Ankle

INTRODUCTION

The aim of the study was to compare maximal force, force steadiness, and the discharge characteristics of motor units in the tibialis anterior (TA) muscle during submaximal isometric contractions for ankle dorsiflexion and adduction of the foot.

METHODS

Nineteen active young adults performed maximal and submaximal isometric dorsiflexion and adduction contractions at five target forces (5%, 10%, 20%, 40%, and 60% maximal voluntary contraction [MVC]). The activity of motor units in TA was recorded by high-density EMG.

RESULTS

The maximal force was similar between dorsiflexion and adduction, despite EMG amplitude for TA being greater ( P < 0.05) during dorsiflexion than adduction. Τhe coefficient of variation (CV) for force (force steadiness) during dorsiflexion was always less ( P < 0.05) than during adduction, except of 5% MVC force. No differences were observed for mean discharge rate; however, the regression between the changes in discharge rate relative to the change of force was significant for dorsiflexion ( R2 = 0.25, P < 0.05) but not for adduction. Discharge variability, however, was usually less during dorsiflexion. The CV for interspike interval was less ( P < 0.05) at 10%, 20%, and 40% MVC but greater at 60% MVC during dorsiflexion than adduction. Similarly, the SD values of the filtered cumulative spike train of the motor units in TA were less ( P < 0.05) at 5%, 10%, 20%, and 40% MVC during dorsiflexion than adduction.

CONCLUSIONS

Although the mean discharge rate of motor units in TA was similar during foot adduction and ankle dorsiflexion, discharge variability was less during dorsiflexion resulting in less accurate performance of the steady adduction contractions. The neural drive to bifunctional muscles differs during their accessory function, which must be considered for training and rehabilitation interventions.

Advancements and Challenges in Non-Invasive Sensor Technologies for Swallowing Assessment: A Review

Dysphagia is a pervasive health issue that impacts diverse demographic groups worldwide, particularly the elderly, stroke survivors, and those suffering from neurological disorders. This condition poses substantial health risks, including malnutrition, respiratory complications, and increased mortality. Additionally, it exacerbates economic burdens by extending hospital stays and escalating healthcare costs. Given that this disorder is frequently underestimated in vulnerable populations, there is an urgent need for enhanced diagnostic and therapeutic strategies. Traditional diagnostic tools such as the videofluoroscopic swallowing study (VFSS) and flexible endoscopic evaluation of swallowing (FEES) require interpretation by clinical experts and may lead to complications. In contrast, non-invasive sensors offer a more comfortable and convenient approach for assessing swallowing function. This review systematically examines recent advancements in non-invasive swallowing function detection devices, focusing on the validation of the device designs and their implementation in clinical practice. Moreover, this review discusses the swallowing process and the associated biomechanics, providing a theoretical foundation for the technologies discussed. It is hoped that this comprehensive overview will facilitate a paradigm shift in swallowing assessments, steering the development of technologies towards more accessible and accurate diagnostic tools, thereby improving patient care and treatment outcomes.

Prediction of 1-year change in knee extension strength by neuromuscular properties in older adults

Improving muscle strength and preventing muscle weakness are important for older adults. The change in strength can be effectively explained by skeletal muscle mass and neural factors. Neural factors are important for older adults because the variation of neural components is greater in older than in young adults, and any decline in strength cannot solely be explained by a decrease in skeletal muscle mass. The purpose of the present study was to investigate whether skeletal muscle mass or motor unit firing properties could explain the change in muscle strength after 1 year. Thirty-eight older adults (75.0 ± 4.7 years, 156.6 ± 7.7 cm, 55.5 ± 9.4 kg, 26 women) performed maximum voluntary knee extension and their skeletal muscle mass was measured using a bioimpedance device. During a submaximal contraction task, high-density surface electromyography was recorded and the signals were decomposed into individual motor unit firing. As an index of motor unit firing properties, the slope and y-intercept (MU intercept) were calculated from the regression line between recruitment thresholds and firing rates in each participant. After 1 year, their maximum knee extension torque was evaluated again. A stepwise multiple regression linear model with sex and age as covariates indicated that MU intercept was a significant explanation with a negative association for the 1-year change in muscle strength (β =  - 0.493, p = 0.004), but not skeletal muscle mass (p = 0.364). The results suggest that neural components might be predictors of increasing and decreasing muscle strength rather than skeletal muscle mass.

Horizontal foot orientation affects the distribution of neural drive between gastrocnemii during plantarflexion, without changing neural excitability

The distribution of activation among muscles from the same anatomical group can be affected by the mechanical constraints of the task, such as limb orientation. For example, the distribution of activation between the gastrocnemius medialis (GM) and lateralis (GL) muscles during submaximal plantarflexion depends on the orientation of the foot in the horizontal plane. The neural mechanisms behind these modulations are not known. The overall aim of this study was to determine whether the excitability of the two gastrocnemius muscles is differentially affected by changes in foot orientation. Nineteen males performed isometric plantarflexions with their foot internally (toes-in) or externally (toes-out) rotated. GM and GL motor unit discharge characteristics were estimated from high-density surface electromyography to estimate neural drive. GM and GL corticospinal excitability and intracortical activity were assessed using transcranial magnetic stimulation through motor-evoked potentials. The efficacy of synaptic transmission between Ia-afferent fibers and α-motoneurons of the GM and GL was evaluated through the Hoffmann reflex. We observed a differential change in neural drive between GM (toes-out > toes-in) and GL (toes-out < toes-in). However, there was no foot orientation-related modulation in corticospinal excitability of the GM or GL, either at the cortical level or through modulation of the efficacy of Ia-α-motoneuron transmission. These results demonstrate that change in the motor pathway excitability is not the mechanism controlling the different distribution of neural drive between GM and GL with foot orientation.NEW & NOTEWORTHY Horizontal foot orientation affects the distribution of neural drive between the gastrocnemii during plantarflexion. There is no foot orientation-related modulation in the corticospinal excitability of the gastrocnemii, either at the cortical level or through modulation of the efficacy of Ia-α-motoneuron transmission. Change in motor pathway excitability is not the mechanism controlling the different distribution of neural drive between gastrocnemius medialis and lateralis with foot orientation.

Adaptive HD-sEMG decomposition: towards robust real-time decoding of neural drive

Objective. Neural interfacing via decomposition of high-density surface electromyography (HD-sEMG) should be robust to signal non-stationarities incurred by changes in joint pose and contraction intensity.Approach. We present an adaptive real-time motor unit decoding algorithm and test it on HD-sEMG collected from the extensor carpi radialis brevis during isometric contractions over a range of wrist angles and contraction intensities. The performance of the algorithm was verified using high-confidence benchmark decompositions derived from concurrently recorded intramuscular electromyography.Main results. In trials where contraction conditions between the initialization and testing data differed, the adaptive decoding algorithm maintained significantly higher decoding accuracies when compared to static decoding methods.Significance. Using "gold standard" verification techniques, we demonstrate the limitations of filter re-use decoding methods and show the necessity of parameter adaptation to achieve robust neural decoding.

Lower limb suspension induces threshold-specific alterations of motor units properties that are reversed by active recovery

PURPOSE

This study aimed to non-invasively test the hypothesis that (a) short-term lower limb unloading would induce changes in the neural control of force production (based on motor units (MUs) properties) in the vastus lateralis muscle and (b) possible changes are reversed by active recovery (AR).

METHODS

Ten young males underwent 10 days of unilateral lower limb suspension (ULLS) followed by 21 days of AR. During ULLS, participants walked exclusively on crutches with the dominant leg suspended in a slightly flexed position (15°-20°) and with the contralateral foot raised by an elevated shoe. The AR was based on resistance exercise (leg press and leg extension) and executed at 70% of each participant's 1 repetition maximum, 3 times/week. Maximal voluntary isometric contraction (MVC) of knee extensors and MUs properties of the vastus lateralis muscle were measured at baseline, after ULLS, and after AR. MUs were identified using high-density electromyography during trapezoidal isometric contractions at 10%, 25%, and 50% of the current MVC, and individual MUs were tracked across the 3 data collection points.

RESULTS

We identified 1428 unique MUs, and 270 of them (18.9%) were accurately tracked. After ULLS, MVC decreased by 29.77%, MUs absolute recruitment/derecruitment thresholds were reduced at all contraction intensities (with changes between the 2 variables strongly correlated), while discharge rate was reduced at 10% and 25% but not at 50% MVC. Impaired MVC and MUs properties fully recovered to baseline levels after AR. Similar changes were observed in the pool of total as well as tracked MUs.

CONCLUSION

Our novel results demonstrate, non-invasively, that 10 days of ULLS affected neural control predominantly by altering the discharge rate of lower-threshold but not of higher-threshold MUs, suggesting a preferential impact of disuse on motoneurons with a lower depolarization threshold. However, after 21 days of AR, the impaired MUs properties were fully restored to baseline levels, highlighting the plasticity of the components involved in neural control.

Time-of-day effects on motor unit firing and muscle contractile properties in humans

Intrinsic factors related to neuromuscular function are time-of-day dependent, but diurnal rhythms in neural and muscular components of the human neuromuscular system remain unclear. The present study aimed to investigate the time-of-day effects on neural excitability and muscle contractile properties by assessing the firing properties of tracked motor units and electrically evoked twitch muscle contraction. In 15 young adults (22.9 ± 4.7 yr), neuromuscular function was measured in the morning (10:00), at noon (13:30), in the evening (17:00), and at night (20:30). Four measurements were completed within 24 h. The measurements consisted of maximal voluntary contraction (MVC) strength of knee extension, recording of high-density surface electromyography (HDsEMG) from the vastus lateralis during ramp-up contraction to 50% of MVC, and evoked twitch torque of knee extensors by electrical stimulation. Recorded HDsEMG signals were decomposed to individual motor unit firing behaviors and the same motor units were tracked among the times of day, and recruitment thresholds and firing rates were calculated. The number of detected and tracked motor units was 127. Motor unit firing rates significantly increased from morning to noon, evening, and night (P < 0.01), but there were no significant differences in recruitment thresholds among the times of day (P > 0.05). Also, there were no significant effects of time of day on evoked twitch torque (P > 0.05). Changes in the motor unit firing rate and evoked twitch torque were not significantly correlated (P > 0.05). These findings suggest that neural excitability may be affected by the time of day, but it did not accompany changes in peripheral contractile properties in a diurnal manner.NEW & NOTEWORTHY We investigated the variations of tracked motor unit firing properties and electrically evoked twitch contraction during the day within 24 h. The variation of motor unit firing rate was observed, and tracked motor unit firing rate increased at noon, in the evening, and at night compared with that in the morning. The variation in motor unit firing rate was independent of changes in twitch contraction. Motor unit firing rate may be affected by diurnal rhythms.

Inferring position of motor units from high-density surface EMG

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

Tutorial: Analysis of central and peripheral motor unit properties from decomposed High-Density surface EMG signals with openhdemg

High-Density surface Electromyography (HD-sEMG) is the most established technique for the non-invasive analysis of single motor unit (MU) activity in humans. It provides the possibility to study the central properties (e.g., discharge rate) of large populations of MUs by analysis of their firing pattern. Additionally, by spike-triggered averaging, peripheral properties such as MUs conduction velocity can be estimated over adjacent regions of the muscles and single MUs can be tracked across different recording sessions. In this tutorial, we guide the reader through the investigation of MUs properties from decomposed HD-sEMG recordings by providing both the theoretical knowledge and practical tools necessary to perform the analyses. The practical application of this tutorial is based on openhdemg, a free and open-source community-based framework for the automated analysis of MUs properties built on Python 3 and composed of different modules for HD-sEMG data handling, visualisation, editing, and analysis. openhdemg is interfaceable with most of the available recording software, equipment or decomposition techniques, and all the built-in functions are easily adaptable to different experimental needs. The framework also includes a graphical user interface which enables users with limited coding skills to perform a robust and reliable analysis of MUs properties without coding.

Alterations of Motor Unit Characteristics Associated With Muscle Fatigue

This study aims to characterize motor unit (MU) features associated with muscle fatigue, using high-density surface electromyography (HD-sEMG). The same MUs recruited before/after, and during muscle fatigue were identified for analysis. The surface location of the innervation zones (IZs) of the MUs was identified from the HD-sEMG bipolar motor unit action potential (MUAP) map. The depth of the MU was also identified from the decay pattern of the MUAP along the muscle fiber transverse direction. Both the surface IZ location and the MU depth information were utilized to ensure the same MU was examined during the contraction before/after muscle fatigue. The MUAP similarity, defined as the correlation coefficient between MUAP morphology, was adopted to reveal the alterations in MU characteristics under the condition of fatigue. The biomarkers of the same MUs were compared before/after fatigue (task 1) at 5%, 10%, and 15% maximal voluntary contraction (MVC) and in the process of continuous fatigue (task 2) at 20% MVC. Our results indicate that the MUAP morphology similarity of the same MUs was 0.91 ± 0.06 (task 1) and 0.93 ± 0.04 (task 2). The results showed that MUAP morphology maintained good stability before/after, and during muscle fatigue. The findings of this study may advance our understanding of the mechanism of MU neuromuscular fatigue.

Motor unit firing rate modulation is more impaired during flexion synergy-driven contractions of the biceps brachii in chronic stroke

Following a hemiparetic stroke, individuals exhibit altered motor unit firing patterns during voluntary muscle contractions, including impairments in firing rate modulation and recruitment. These individuals also exhibit abnormal muscle coactivation through multi-joint synergies (e.g., flexion synergy). Here, we investigate whether motor unit firing activity during flexion synergy-driven contractions of the paretic biceps brachii differs from that of voluntary contractions and use these differences to predict changes in descending motor commands. To accomplish this, we characterized motor unit firing patterns of the biceps brachii in individuals with chronic hemiparetic stroke during voluntary isometric elbow flexion contractions in the paretic and non-paretic limbs, as well as during contractions driven by voluntary effort and by flexion synergy expression in the paretic limb. We observed significant reductions in motor unit firing rate modulation from the non-paretic to paretic limb (non-paretic - paretic: 0.14 pps/%MVT, 95% CI: [0.09 0.19]) that were further reduced during synergy-driven contractions (voluntary paretic - synergy driven: 0.19 pps/%MVT, 95% CI: [0.14 0.25]). Moreover, using recently developed metrics, we evaluated how a stroke-induced reliance on indirect motor pathways alters the inputs that motor units receive and revealed progressive increases in neuromodulatory and inhibitory drive to the motor pool in the paretic limb, with the changes greatest during synergy-driven contractions. These findings suggest that an interplay between heightened neuromodulatory drive and alterations in inhibitory command structure may account for the observed motor unit impairments, further illuminating underlying neural mechanisms involved in the flexion synergy and its impact on motor unit firing patterns post-stroke.

Characterization of Motor Unit Firing and Twitch Properties for Decoding Musculoskeletal Force in the Human Ankle Joint In Vivo

Understanding how motor units (MUs) contribute to skeletal mechanical force is crucial for unraveling the underlying mechanism of human movement. Alterations in MU firing, contractile and force-generating properties emerge in response to physical training, aging or injury. However, how changes in MU firing and twitch properties dictate skeletal muscle force generation in healthy and impaired individuals remains an open question. In this work, we present a MU-specific approach to identify firing and twitch properties of MU samples and employ them to decode musculoskeletal function in vivo. First, MU firing events were decomposed offline from high-density electromyography (HD-EMG) of six lower leg muscles involved in ankle plantar-dorsi flexion. We characterized their twitch responses based on the statistical distributions of their firing properties and employed them to compute MU-specific activation dynamics. Subsequently, we decoded ankle joint moments by linking our framework to a subject-specific musculoskeletal model. We validated our approach at different ankle positions and levels of activation and compared it with traditional EMG-driven models. Our proposed MU-specific formulation achieves higher generalization across conditions than the EMG-driven models, with significantly lower coefficients of variation in torque predictions. Furthermore, our approach shows distinct neural strategies across a large repertoire of contractile conditions in different muscles. Our proposed approach may open new avenues for characterizing the relationship between MU firing and twitch properties and their influence on force capacity. This can facilitate the development of targeted rehabilitation strategies tailored to individuals with specific neuromuscular conditions.

Two motor neuron synergies, invariant across ankle joint angles, activate the triceps surae during plantarflexion

Recent studies have suggested that the nervous system generates movements by controlling groups of motor neurons (synergies) that do not always align with muscle anatomy. In this study, we determined whether these synergies are robust across tasks with different mechanical constraints. We identified motor neuron synergies using principal component analysis (PCA) and cross-correlations between smoothed discharge rates of motor neurons. In part 1, we used simulations to validate these methods. The results suggested that PCA can accurately identify the number of common inputs and their distribution across active motor neurons. Moreover, the results confirmed that cross-correlation can separate pairs of motor neurons that receive common inputs from those that do not receive common inputs. In part 2, 16 individuals performed plantarflexion at three ankle angles while we recorded EMG signals from the gastrocnemius lateralis (GL) and medialis (GM) and the soleus (SOL) with grids of surface electrodes. The PCA revealed two motor neuron synergies. These motor neuron synergies were relatively stable, with no significant differences in the distribution of motor neuron weights across ankle angles (P = 0.62). When the cross-correlation was calculated for pairs of motor units tracked across ankle angles, we observed that only 13.0% of pairs of motor units from GL and GM exhibited significant correlations of their smoothed discharge rates across angles, confirming the low level of common inputs between these muscles. Overall, these results highlight the modularity of movement control at the motor neuron level, suggesting a sensible reduction of computational resources for movement control. KEY POINTS: The CNS might generate movements by activating groups of motor neurons (synergies) with common inputs. We show here that two main sources of common inputs drive the motor neurons innervating the triceps surae muscles during isometric ankle plantarflexions. We report that the distribution of these common inputs is globally invariant despite changing the mechanical constraints of the tasks, i.e. the ankle angle. These results suggest the functional relevance of the modular organization of the CNS to control movements.

Extracting Individual Muscle Drive and Activity From High-Density Surface Electromyography Signals Based on the Center of Gravity of Motor Unit

Neural interfacing has played an essential role in advancing our understanding of fundamental movement neurophysiology and the development of human-machine interface. However, direct neural interfaces from brain and nerve recording are currently limited in clinical areas for their invasiveness and high selectivity. Here, we applied the surface electromyogram (EMG) in studying the neural control of movement and proposed a new non-invasive way of extracting neural drive to individual muscles. Sixteen subjects performed isometric contractions to complete six hand tasks. High-density surface EMG signals (256 channels in total) recorded from the forearm muscles were decomposed into motor unit firing trains. The location of each decomposed motor unit was represented by its center of gravity and was put into clustering for distinct muscle regions. All the motor units in the same cluster served as a muscle-specific motor pool from which individual muscle drive could be extracted directly. Moreover, we cross-validated the self-clustered muscle regions by magnetic resonance imaging (MRI) recorded from the subjects' forearms. All motor units that fall within the MRI region are considered correctly clustered. We achieved a clustering accuracy of 95.72% ± 4.01% for all subjects. We provided a new framework for collecting experimental muscle-specific drives and generalized the way of surface electrode placement without prior knowledge of the targeting muscle architecture.

Sensorimotor integration is affected by acute whole-body vibration: a coherence study

Introduction: Several whole-body vibration (WBV) effects on performance have been related to potential changes in the neural drive, motor unit firing rate, and sensorimotor integration. In the present paper, motor unit coherence analysis was performed to detect the source of neural modulation based on the frequency domain. Methods: Thirteen men [25 ± 2.1 years; Body Mass Index (BMI) = 23.9 ± 1.3 kg m2; maximal voluntary force (MVF): 324.36 ± 41.26 N] performed sustained contractions of the Tibialis Anterior (TA) at 10%MVF before and after acute WBV. The vibrating stimulus was applied barefoot through a platform to target the TA. High-Density surface Electromyography (HDsEMG) was used to record the myoelectrical activity of TA to evaluate coherence from motor unit cumulative spike-trains (CSTs). Results: Mean coherence showed a significant decrease in the alpha and low-beta bandwidths (alpha: from 0.143 ± 0.129 to 0.132 ± 0.129, p = 0.035; low-beta: from 0.117 ± 0.039 to 0.086 ± 0.03, p = 0.0001), whereas no significant changes were found in the other ones (p > 0.05). The discharge rate (DR) and the Force Covariance (CovF%) were not significantly affected by acute WBV exposure (p > 0.05). Discussion: According to the significant effects found in alpha and low-beta bandwidths, which reflect sensorimotor integration parameters, accompanied by no differences in the DR and CovF%, the present results underlined that possible neural mechanisms at the base of the previously reported performance enhancements following acute WBV are likely based on sensorimotor integration rather than direct neural drive modulation.

Acute adaptation of central and peripheral motor unit features to exercise-induced fatigue differs with concentric and eccentric loading

NEW FINDINGS

What is the central question of this study? Conflicting evidence exists on motor unit (MU) firing rate in response to exercise-induced fatigue, possibly due to the contraction modality used: Do MU properties adapt similarly following concentric and eccentric loading? What is the main finding and its importance? MU firing rate increased following eccentric loading only despite a decline in absolute force. Force steadiness deteriorated following both loading methods. Central and peripheral MU features are altered in a contraction type-dependant manner, which is an important consideration for training interventions.

ABSTRACT

Force output of muscle is partly mediated by the adjustment of motor unit (MU) firing rate (FR). Disparities in MU features in response to fatigue may be influenced by contraction type, as concentric (CON) and eccentric (ECC) contractions demand variable amounts of neural input, which alters the response to fatigue. This study aimed to determine the effects of fatigue following CON and ECC loading on MU features of the vastus lateralis (VL). High-density surface (HD-sEMG) and intramuscular (iEMG) electromyography were used to record MU potentials (MUPs) from bilateral VLs of 12 young volunteers (six females) during sustained isometric contractions at 25% and 40% of the maximum voluntary contraction (MVC), before and after completing CON and ECC weighted stepping exercise. Multi-level mixed effects linear regression models were performed with significance assumed as P < 0.05. MVC decreased in both CON and ECC legs post-exercise (P < 0.0001), as did force steadiness at both 25% and 40% MVC (P < 0.004). MU FR increased in ECC at both contraction levels (P < 0.001) but did not change in CON. FR variability increased in both legs at 25% and 40% MVC following fatigue (P < 0.01). From iEMG measures at 25% MVC, MUP shape did not change (P > 0.1) but neuromuscular junction transmission instability increased in both legs (P < 0.04), and markers of fibre membrane excitability increased following CON only (P = 0.018). These data demonstrate that central and peripheral MU features are altered following exercise-induced fatigue and differ according to exercise modality. This is important when considering interventional strategies targeting MU function.

Adjustments in the motor unit discharge behavior following neuromuscular electrical stimulation compared to voluntary contractions

Introduction: The application of neuromuscular electrical stimulation superimposed on voluntary muscle contractions (NMES+) has demonstrated a considerable potential to enhance or restore muscle function in both healthy and individuals with neurological or orthopedic disorders. Improvements in muscle strength and power have been commonly associated with specific neural adaptations. In this study, we investigated changes in the discharge characteristics of the tibialis anterior motor units, following three acute exercises consisting of NMES+, passive NMES and voluntary isometric contractions alone. Methods: Seventeen young participants participated in the study. High-density surface electromyography was used to record myoelectric activity in the tibialis anterior muscle during trapezoidal force trajectories involving isometric contractions of ankle dorsi flexors with target forces set at 35, 50% and 70% of maximal voluntary isometric contraction (MVIC). From decomposition of the electromyographic signal, motor unit discharge rate, recruitment and derecruitment thresholds were extracted and the input-output gain of the motoneuron pool was estimated. Results: Global discharge rate increased following the isometric condition compared to baseline at 35% MVIC while it increased after all experimental conditions at 50% MVIC target force. Interestingly, at 70% MVIC target force, only NMES + led to greater discharge rate compared to baseline. Recruitment threshold decreased after the isometric condition, although only at 50% MVIC. Input-output gain of the motoneurons of the tibialis anterior muscle was unaltered after the experimental conditions. Discussion: These results indicated that acute exercise involving NMES + induces an increase in motor unit discharge rate, particularly when higher forces are required. This reflects an enhanced neural drive to the muscle and might be strongly related to the distinctive motor fiber recruitment characterizing NMES+.

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.

Facilitation of sensory transmission to motoneurons during cortical or sensory-evoked primary afferent depolarization (PAD) in humans

Sensory and corticospinal tract (CST) pathways activate spinal GABAergic interneurons that have axoaxonic connections onto proprioceptive (Ia) afferents that cause long-lasting depolarizations (termed primary afferent depolarization, PAD). In rodents, sensory-evoked PAD is produced by GABAA receptors at nodes of Ranvier in Ia afferents, rather than at presynaptic terminals, and facilitates spike propagation to motoneurons by preventing branch-point failures, rather than causing presynaptic inhibition. We examined in 40 human participants whether putative activation of Ia-PAD by sensory or CST pathways can also facilitate Ia afferent activation of motoneurons via the H-reflex. H-reflexes in several leg muscles were facilitated by prior conditioning from low-threshold proprioceptive, cutaneous or CST pathways, with a similar long-lasting time course (∼200 ms) to phasic PAD measured in rodent Ia afferents. Long trains of cutaneous or proprioceptive afferent conditioning produced longer-lasting facilitation of the H-reflex for up to 2 min, consistent with tonic PAD in rodent Ia afferents mediated by nodal α5-GABAA receptors for similar stimulation trains. Facilitation of H-reflexes by this conditioning was likely not mediated by direct facilitation of the motoneurons because isolated stimulation of sensory or CST pathways did not alone facilitate the tonic firing rate of motor units. Furthermore, cutaneous conditioning increased the firing probability of single motor units (motoneurons) during the H-reflex without increasing their firing rate at this time, indicating that the underlying excitatory postsynaptic potential was more probable, but not larger. These results are consistent with sensory and CST pathways activating nodal GABAA receptors that reduce intermittent failure of action potentials propagating into Ia afferent branches. KEY POINTS: Controlled execution of posture and movement requires continually adjusted feedback from peripheral sensory pathways, especially those that carry proprioceptive information about body position, movement and effort. It was previously thought that the flow of proprioceptive feedback from Ia afferents was only reduced by GABAergic neurons in the spinal cord that sent axoaxonic projections to the terminal endings of sensory axons (termed GABAaxo neurons). Based on new findings in rodents, we provide complementary evidence in humans to suggest that sensory and corticospinal pathways known to activate GABAaxo neurons that project to dorsal parts of the Ia afferent also increase the flow of proprioceptive feedback to motoneurons in the spinal cord. These findings support a new role for spinal GABAaxo neurons in facilitating afferent feedback to the spinal cord during voluntary or reflexive movements.

Decoding firings of a large population of human motor units from high-density surface electromyogram in response to transcranial magnetic stimulation

We describe a novel application of methodology for high-density surface electromyography (HDsEMG) decomposition to identify motor unit (MU) firings in response to transcranial magnetic stimulation (TMS). The method is based on the MU filter estimation from HDsEMG decomposition with convolution kernel compensation during voluntary isometric contractions and its application to contractions elicited by TMS. First, we simulated synthetic HDsEMG signals during voluntary contractions followed by simulated motor evoked potentials (MEPs) recruiting an increasing proportion of the motor pool. The estimation of MU filters from voluntary contractions and their application to elicited contractions resulted in high (>90%) precision and sensitivity of MU firings during MEPs. Subsequently, we conducted three experiments in humans. From HDsEMG recordings in first dorsal interosseous and tibialis anterior muscles, we demonstrated an increase in the number of identified MUs during MEPs evoked with increasing stimulation intensity, low variability in the MU firing latency and a proportion of MEP energy accounted for by decomposition similar to voluntary contractions. A negative relationship between the MU recruitment threshold and the number of identified MU firings was exhibited during the MEP recruitment curve, suggesting orderly MU recruitment. During isometric dorsiflexion we also showed a negative association between voluntary MU firing rate and the number of firings of the identified MUs during MEPs, suggesting a decrease in the probability of MU firing during MEPs with increased background MU firing rate. We demonstrate accurate identification of a large population of MU firings in a broad recruitment range in response to TMS via non-invasive HDsEMG recordings. KEY POINTS: Transcranial magnetic stimulation (TMS) of the scalp produces multiple descending volleys, exciting motor pools in a diffuse manner. The characteristics of a motor pool response to TMS have been previously investigated with intramuscular electromyography (EMG), but this is limited in its capacity to detect many motor units (MUs) that constitute a motor evoked potential (MEP) in response to TMS. By simulating synthetic signals with known MU firing patterns, and recording high-density EMG signals from two human muscles, we show the feasibility of identifying firings of many MUs that comprise a MEP. We demonstrate the identification of firings of a large population of MUs in the broad recruitment range, up to maximal MEP amplitude, with fewer required stimuli compared to intramuscular EMG recordings. The methodology demonstrates an emerging possibility to study responses to TMS on a level of individual MUs in a non-invasive manner.

Intrinsic motor neuron excitability is increased after resistance training in older adults

This study investigated the effects of high-intensity resistance training on estimates of the motor neuron persistent inward current (PIC) in older adults. Seventeen participants (68.5 ± 2.8 yr) completed a 2-wk nonexercise control period followed by 6 wk of resistance training. Surface electromyographic signals were collected with two 32-channel electrodes placed over soleus to investigate motor unit discharge rates. Paired motor unit analysis was used to calculate delta frequency (ΔF) as an estimate of PIC amplitudes during 1) triangular-shaped contractions to 20% of maximum torque capacity and 2) trapezoidal- and triangular-shaped contractions to 20% and 40% of maximum torque capacity, respectively, to understand their ability to modulate PICs as contraction intensity increases. Maximal strength and functional capacity tests were also assessed. For the 20% triangular-shaped contractions, ΔF [0.58-0.87 peaks per second (pps); P ≤ 0.015] and peak discharge rates (0.78-0.99 pps; P ≤ 0.005) increased after training, indicating increased PIC amplitude. PIC modulation also improved after training. During the control period, mean ΔF differences between 20% trapezoidal-shaped and 40% triangular-shaped contractions were 0.09-0.18 pps (P = 0.448 and 0.109, respectively), which increased to 0.44 pps (P < 0.001) after training. Also, changes in ΔF showed moderate to very large correlations (r = 0.39-0.82) with changes in peak discharge rates and broad measures of motor function. Our findings indicate that increased motor neuron excitability is a potential mechanism underpinning training-induced improvements in motor neuron discharge rate, strength, and motor function in older adults. This increased excitability is likely mediated by enhanced PIC amplitudes, which are larger at higher contraction intensities.NEW & NOTEWORTHY Resistance training elicited important alterations in soleus intrinsic motor neuronal excitability, likely mediated by enhanced persistent inward current (PIC) amplitude, in older adults. Estimates of PICs increased after the training period, accompanied by an enhanced ability to increase PIC amplitudes at higher contraction intensities. Our data also suggest that changes in PIC contribution to self-sustained discharging may contribute to increases in motor neuron discharge rates, maximal strength, and functional capacity in older adults after resistance training.

Consensus for experimental design in electromyography (CEDE) project: Single motor unit matrix

The analysis of single motor unit (SMU) activity provides the foundation from which information about the neural strategies underlying the control of muscle force can be identified, due to the one-to-one association between the action potentials generated by an alpha motor neuron and those received by the innervated muscle fibers. Such a powerful assessment has been conventionally performed with invasive electrodes (i.e., intramuscular electromyography (EMG)), however, recent advances in signal processing techniques have enabled the identification of single motor unit (SMU) activity in high-density surface electromyography (HDsEMG) recordings. This matrix, developed by the Consensus for Experimental Design in Electromyography (CEDE) project, provides recommendations for the recording and analysis of SMU activity with both invasive (needle and fine-wire EMG) and non-invasive (HDsEMG) SMU identification methods, summarizing their advantages and disadvantages when used during different testing conditions. Recommendations for the analysis and reporting of discharge rate and peripheral (i.e., muscle fiber conduction velocity) SMU properties are also provided. The results of the Delphi process to reach consensus are contained in an appendix. This matrix is intended to help researchers to collect, report, and interpret SMU data in the context of both research and clinical applications.

A geometric approach to quantifying the neuromodulatory effects of persistent inward currents on individual motor unit discharge patterns

Objective.All motor commands flow through motoneurons, which entrain control of their innervated muscle fibers, forming a motor unit (MU). Owing to the high fidelity of action potentials within MUs, their discharge profiles detail the organization of ionotropic excitatory/inhibitory as well as metabotropic neuromodulatory commands to motoneurons. Neuromodulatory inputs (e.g. norepinephrine, serotonin) enhance motoneuron excitability and facilitate persistent inward currents (PICs). PICs introduce quantifiable properties in MU discharge profiles by augmenting depolarizing currents upon activation (i.e. PIC amplification) and facilitating discharge at lower levels of excitatory input than required for recruitment (i.e. PIC prolongation).Approach. Here, we introduce a novel geometric approach to estimate neuromodulatory and inhibitory contributions to MU discharge by exploiting discharge non-linearities introduced by PIC amplification during time-varying linear tasks. In specific, we quantify the deviation from linear discharge ('brace height') and the rate of change in discharge (i.e. acceleration slope, attenuation slope, angle). We further characterize these metrics on a simulated motoneuron pool with known excitatory, inhibitory, and neuromodulatory inputs and on human MUs (number of MUs; Tibialis Anterior: 1448, Medial Gastrocnemius: 2100, Soleus: 1062, First Dorsal Interosseus: 2296).Main results. In the simulated motor pool, we found brace height and attenuation slope to consistently indicate changes in neuromodulation and the pattern of inhibition (excitation-inhibition coupling), respectively, whereas the paired MU analysis (ΔF) was dependent on both neuromodulation and inhibition pattern. Furthermore, we provide estimates of these metrics in human MUs and show comparable variability in ΔFand brace height measures for MUs matched across multiple trials.Significance. Spanning both datasets, we found brace height quantification to provide an intuitive method for achieving graded estimates of neuromodulatory and inhibitory drive to individual MUs. This complements common techniques and provides an avenue for decoupling changes in the level of neuromodulatory and pattern of inhibitory motor commands.

The acute effects of whole-body vibration on motor unit recruitment and discharge properties

Introduction: several studies have reported improved neuromuscular parameters in response to whole-body vibration (WBV). This is likely achieved by modulation of the central nervous system (CNS). Reduced recruitment threshold (RT), which is the % of Maximal Voluntary Force (%MVF) at which a given Motor Unit (MU) is recruited, may be responsible for the force/power improvements observed in several studies. Methods: 14 men (25 ± 2.3 years; BMI = 23.3 ± 1.5 kg m² MVF: 319.82 ± 45.74 N) performed trapezoidal isometric contractions of the tibialis anterior (TA) at 35-50-70 %MVF before and after three conditions: WBV, STAND (standing posture), and CNT (no intervention). The vibration was applied through a platform for targeting the TA. High-density surface electromyography (HDsEMG) recordings and analysis were used to detect changes in the RT and Discharge Rate (DR) of the MUs. Results: Mean motor unit recruitment threshold (MURT) reached 32.04 ± 3.28 %MVF before and 31.2 ± 3.72 %MVF after WBV, with no significant differences between conditions (p > 0.05). Additionally, no significant changes were found in the mean motor unit discharge rate (before WBV: 21.11 ± 2.94 pps; after WBV: 21.19 ± 2.17 pps). Discussion: The present study showed no significant changes in motor unit properties at the base of neuromuscular changes documented in previous studies. Further investigations are needed to understand motor unit responses to different vibration protocols and the chronic effect of vibration exposure on motor control strategies.

Common synaptic input between motor units from the lateral and medial posterior soleus compartments does not differ from that within each compartment

The human soleus muscle is anatomically divided into four separate anatomical compartments. The functional role of this compartmentalization remains unclear. Here, we tested the hypothesis that the common synaptic input to motor units between the medial and lateral posterior compartments is less than within each compartment. Fourteen male participants performed three different heel-raise tasks that were considered to place a different mechanical demand on the medial and lateral soleus compartments. High-density electromyography (EMG) signals from the medial and lateral soleus compartments and the medial gastrocnemius of the right leg were decomposed into individual motor unit spike trains. The coherence between cumulative spike trains of the motor units was estimated. The coherence analysis was also repeated for motor units that were matched across all three tasks. Furthermore, we calculated the ratio of significant correlations between the spike trains of pairs of motor units. We observed that the coherence between motor units of the two soleus compartments was similar as the coherence between motor units within each compartment, regardless of the task. The correlation analysis performed on pairs of motor units confirmed these results. We conclude that the level of common synaptic input between the motor units innervating the medial and lateral posterior soleus compartment is not different than the common synaptic input between motor units innervating each of these compartments, which contrasts with findings from previous studies on finger muscles. This suggests that there is no independent neural control for the individual posterior soleus compartments.NEW & NOTEWORTHY The human soleus muscle is anatomically subdivided into four compartments. The functional role for this compartmentalization remains unknown. Here, we showed that, contrary to previous findings in finger muscles, the common synaptic input between motor units innervating the medial and lateral posterior soleus compartment was similar as that between motor units within the individual compartments. We suggest that the contradictory findings with other compartmentalized muscles may be explained by differences in muscle-tendon anatomy and function.

Firing behavior of single motor units of the tibialis anterior in human walking as non-invasively revealed by HDsEMG decomposition

Objective.Investigation of the firing behavior of motor units (MUs) provides essential neuromuscular control information because MUs are the smallest organizational component of the neuromuscular system. The MUs activated during human infants' leg movements and rodent locomotion, mainly controlled by the spinal central pattern generator (CPG), show highly synchronous firing. In addition to spinal CPGs, the cerebral cortex is involved in neuromuscular control during walking in human adults. Based on the difference in the neural control mechanisms of locomotion between rodent, human infants and adults, MU firing behavior during adult walking probably has some different features from the other populations. However, so far, the firing activity of MUs in human adult walking has been largely unknown due to technical issues.Approach.Recent technical advances allow noninvasive investigation of MU firing by high-density surface electromyogram (HDsEMG) decomposition. We investigated the MU firing behavior of the tibialis anterior (TA) muscle during walking at a slow speed by HDsEMG decomposition.Main results.We found recruitment threshold modulation of MU between walking and steady isometric contractions. Doublet firings, and gait phase-specific firings were also observed during walking. We also found high MU synchronization during walking over a wide range of frequencies, probably including cortical and spinal CPG-related components. The amount of MU synchronization was modulated between the gait phases and motor tasks. These results suggest that the central nervous system flexibly controls MU firing to generate appropriate force of TA during human walking.Significance.This study revealed the MU behavior during walking at a slow speed and demonstrated the feasibility of noninvasive investigation of MUs during dynamic locomotor tasks, which will open new frontiers for the study of neuromuscular systems in the fields of neuroscience and biomedical engineering.

Neural inputs from spinal motor neurons to lateralis vastus muscle: Comparison between sprinters and nonathletes

The adaptation of neural contractile properties has been observed in previous work. However, the neural changes on the motor unit (MU) level remain largely unknown. Voluntary movements are controlled through the precise activation of MU populations. In this work, we estimate the neural inputs from the spinal motor neurons to the muscles during isometric contractions and characterize the neural adaptation during training by comparing the MU properties decomposed from sprinters and nonathletes. Twenty subjects were recruited and divided into two groups. The high-density surface electromyography (EMG) signals were recorded from the lateralis vastus muscle during the isometric contraction of knee extension and were then decomposed into MU spike trains. Each MU’s action potentials and discharge properties were extracted for comparison across subject groups and tasks. A total of 1097 MUs were identified from all subjects. Results showed that the discharge rates and amplitudes of MUAPs from athletes were significantly higher than those from nonathletes. These results demonstrate the neural adaptations in physical training at the MU population level and indicate the great potential of EMG decomposition in physiological investigations.

Interfacing Motor Units in Nonhuman Primates Identifies a Principal Neural Component for Force Control Constrained by the Size Principle

Motor units convert the last neural code of movement into muscle forces. The classic view of motor unit control is that the CNS sends common synaptic inputs to motoneuron pools and that motoneurons respond in an orderly fashion dictated by the size principle. This view, however, is in contrast with the large number of dimensions observed in motor cortex, which may allow individual and flexible control of motor units. Evidence for flexible control of motor units may be obtained by tracking motor units longitudinally during tasks with some level of behavioral variability. Here we identified and tracked populations of motor units in the brachioradialis muscle of two macaque monkeys during 10 sessions spanning >1 month with a broad range of rate of force development (1.8-38.6 N · m · s-1). We found a very stable recruitment order and discharge characteristics of the motor units over sessions and contraction trials. The small deviations from orderly recruitment were fully predicted by the motor unit recruitment intervals, so that small shifts in recruitment thresholds happened only during contractions at a high rate of force development. Moreover, we also found that one component explained more than ∼50% of the motor unit discharge rate variance, and that the remaining components represented a time-shifted version of the first. In conclusion, our results show that the recruitment of motoneurons is determined by the interplay of the size principle and common input and that this recruitment scheme is not violated over time or by the speed of the contractions.SIGNIFICANCE STATEMENT With a new noninvasive high-density electromyographic framework, we show the activity of motor unit ensembles in macaques during voluntary contractions. The discharge characteristics of brachioradialis motor units revealed a relatively fixed recruitment order and discharge characteristics across days and rate of force developments. These results were further confirmed through invasive axonal stimulation and recordings of intramuscular electromyographic activity from 16 arm muscles. The study shows for the first time the feasibility of longitudinal noninvasive motor unit interfacing and tracking of the same motor units in nonhuman primates.

Cervical spine and muscle adaptation after spaceflight and relationship to herniation risk: protocol from 'Cervical in Space' trial

Background

Astronauts have a higher risk of cervical intervertebral disc herniation. Several mechanisms have been attributed as causative factors for this increased risk. However, most of the previous studies have examined potential causal factors for lumbar intervertebral disc herniation only. Hence, we aim to conduct a study to identify the various changes in the cervical spine that lead to an increased risk of cervical disc herniation after spaceflight.

Methods

A cohort study with astronauts will be conducted. The data collection will involve four main components: a) Magnetic resonance imaging (MRI); b) cervical 3D kinematics; c) an Integrated Protocol consisting of maximal and submaximal voluntary contractions of the neck muscles, endurance testing of the neck muscles, neck muscle fatigue testing and questionnaires; and d) dual energy X-ray absorptiometry (DXA) examination. Measurements will be conducted at several time points before and after astronauts visit the International Space Station. The main outcomes of interest are adaptations in the cervical discs, muscles and bones.

Discussion

Astronauts are at higher risk of cervical disc herniation, but contributing factors remain unclear. The results of this study will inform future preventive measures for astronauts and will also contribute to the understanding of intervertebral disc herniation risk in the cervical spine for people on Earth. In addition, we anticipate deeper insight into the aetiology of neck pain with this research project.

Trial registration

German Clinical Trials Register, DRKS00026777. Registered on 08 October 2021.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12891-022-05684-0.

Training induced improvements in knee extensor force accuracy are associated with reduced vastus lateralis motor unit firing variability

NEW FINDINGS

What is the central question of this study? We aimed to determine levels of bilateral knee extensor force accuracy and any subsequent alterations to central and/or peripheral motor unit features, following 4 weeks of unilateral force accuracy training. What is the main finding and its importance? In the trained limb only, knee extensor force tracking accuracy improved with reduced motor unit firing rate variability in the vastus lateralis, and no change to neuromuscular junction transmission instability. Interventional strategies to improve force accuracy may be directed to older/clinical populations where such improvements may aid performance of daily living activities.

ABSTRACT

Background Muscle force output during sustained submaximal isometric contractions fluctuates around an average value and is partly influenced by variation in motor unit (MU) firing rates. MU firing rate (FR) variability seemingly reduces following exercise training interventions, however, much less is known with respect to peripheral MU properties. We therefore investigated whether targeted force accuracy training could lead to improved muscle functional capacity and control, in addition to determining any alterations of individual MU features. Methods Ten healthy participants (7 females, 3 males, 27±6 years, 170±8 cm, 69±16kg) underwent a 4-week supervised, unilateral knee extensor force accuracy training intervention. The coefficient of variation for force (FORCE^CoV ) and sinusoidal wave force tracking accuracy (FORCE^Sinu ) were determined at 25% maximal voluntary contraction (MVC) pre- and post-training. Intramuscular electromyography was utilised to record individual MU potentials from the vastus lateralis (VL) muscles at 25% MVC during sustained contractions, pre- and post-training. Results Knee extensor muscle strength remained unchanged following training, with no improvements in unilateral leg-balance. FORCE^CoV and FORCE^Sinu significantly improved in only the trained knee extensors by ∼13% (p = 0.01) and ∼30% (p<0.0001) respectively. MU FR variability significantly reduced in the trained VL by ∼16% (n = 8; p = 0.001), with no further alterations to MU FR or neuromuscular junction transmission instability. Conclusion Our results suggest muscle force control and tracking accuracy is a trainable characteristic in the knee extensors, which is likely explained by the reduction in MU FR variability which was apparent in the trained limb only. This article is protected by copyright. All rights reserved.

Neuromuscular and structural tendon adaptations after 6 weeks of either concentric or eccentric exercise in individuals with non-insertional Achilles tendinopathy: protocol for a randomised controlled trial

INTRODUCTION

There is limited evidence on the neural strategies employed by the central nervous system to control muscle force in the presence of non-insertional Achilles tendinopathy (NIAT). Additionally, the neuromuscular mechanisms by which exercise may help to resolve tendon pain remain unclear.

OBJECTIVE

This study aims to first establish changes in the gastrocnemius-soleus motor unit firing properties after applying a training protocol of 6 weeks based on either controlled eccentric or concentric contractions in individuals with NIAT. Second, we want to determine changes in the level of pain and function and mechanical and structural properties of the Achilles tendon after applying the same training protocol. Additionally, we want to compare these variables at baseline between individuals with NIAT and asymptomatic controls.

METHODS AND ANALYSIS

A total of 26 individuals with chronic (>3 months) NIAT and 13 healthy controls will participate in the study. Individuals with NIAT will be randomised to perform eccentric or concentric training for 6 weeks. Motor unit firing properties of the medial gastrocnemius, lateral gastrocnemius and soleus muscles will be assessed using high-density surface electromyography, as well as Achilles tendon length, cross-sectional area, thickness and stiffness using B-mode ultrasonography and shear wave elastography. Moreover, participants will complete a battery of questionnaires to document their level of pain and function.

ETHICS AND DISSEMINATION

Ethical approval (ERN-20-0604A) for the study was obtained from the Science, Technology, Engineering and Mathematics Ethical Review Committee of the University of Birmingham. The results of the study will be published in peer-review journals.

TRIAL REGISTRATION NUMBER

ISRCTN46462385.

Neural decoding from surface high-density EMG signals: influence of anatomy and synchronization on the number of identified motor units

OBJECTIVE

High-density surface electromyography (HD-sEMG) allows the reliable identification of individual motor unit (MU) action potentials. Despite the accuracy in decomposition, there is a large variability in the number of identified MUs across individuals and exerted forces. Here we present a systematic investigation of the anatomical and neural factors that determine this variability.

APPROACH

We investigated factors of influence on HD-sEMG decomposition, such as synchronization of MU discharges, distribution of MU territories, muscle-electrode distance (MED - subcutaneous adipose tissue thickness), maximum anatomical cross-sectional area (ACSA_max), and fiber CSA. For this purpose, we recorded HD-sEMG signals, ultrasound and, magnetic resonance images, and took a muscle biopsy from the biceps brachii muscle from 30 male participants drawn from two groups to ensure variability within the factors - untrained-controls (UT=14) and strength-trained individuals (ST=16). Participants performed isometric ramp contractions with elbow flexors (at 15, 35, 50 and 70% maximum voluntary torque - MVT). We assessed the correlation between the number of accurately detected MUs by HD-sEMG decomposition and each measured parameter, for each target force level. Multiple regression analysis was then applied.

MAIN RESULTS

ST subjects showed lower MED (UT=5.1±1.4 mm; ST=3.8±0.8 mm) and a greater number of identified motor units (UT:21.3±10.2 vs ST:29.2±11.8 MUs/subject across all force levels). The entire cohort showed a negative correlation between MED and the number of identified MUs at low forces (r= -0.6, p=0.002 at 15%MVT). Moreover, the number of identified MUs was positively correlated to the distribution of MU territories (r=0.56, p=0.01) and ACSA_max(r=0.48, p=0.03) at 15%MVT. By accounting for all anatomical parameters, we were able to partly predict the number of decomposed MUs at low but not at high forces.

SIGNIFICANCE

Our results confirmed the influence of subcutaneous tissue on the quality of HD-sEMG signals and demonstrated that MU spatial distribution and ACSA_maxare also relevant parameters of influence for current decomposition algorithms.

Detecting motor unit abnormalities in amyotrophic lateral sclerosis using high-density surface EMG

OBJECTIVE

The purpose of this study was to detect specific motor unit (MU) abnormalities in people with amyotrophic lateral sclerosis (ALS) compared to controls using high-density surface electromyography (HD-SEMG).

METHODS

Sixteen people with ALS and 16 control subjects. The participants performed ramp up and sustained contractions at 30% of their maximal voluntary contraction. HD-SEMG signals were recorded in the vastus lateralis muscle and decomposed into individual MU firing behavior using a convolution blind source separation method.

RESULTS

In total, 339 MUs were detected (people with ALS; n = 93, control subjects; n = 246). People with ALS showed significantly higher mean firing rate, recruitment threshold, coefficient of variation of the MU firing rate, MU firing rate at recruitment, and motoneurons excitability than those of control subjects (p < 0.001). The number of MU, MU firing rate, recruitment threshold, and MU firing rate at recruitment were significantly correlated with disease severity (p < 0.001). Multivariable analysis revealed that an increased MU firing rate at recruitment was independently associated with ALS.

CONCLUSIONS

These results suggest increased excitability at recruitment, which is consistent with neurodegeneration results in a compensatory increase in MU activity.

SIGNIFICANCE

Abnormal MU firing behavior provides an important physiological index for understanding the pathophysiology of ALS.

Physical and electrophysiological motor unit characteristics are revealed with simultaneous high-density electromyography and ultrafast ultrasound imaging

Electromyography and ultrasonography provide complementary information about electrophysiological and physical (i.e. anatomical and mechanical) muscle properties. In this study, we propose a method to assess the electrical and physical properties of single motor units (MUs) by combining High-Density surface Electromyography (HDsEMG) and ultrafast ultrasonography (US). Individual MU firings extracted from HDsEMG were used to identify the corresponding region of muscle tissue displacement in US videos. The time evolution of the tissue velocity in the identified region was regarded as the MU tissue displacement velocity. The method was tested in simulated conditions and applied to experimental signals to study the local association between the amplitude distribution of single MU action potentials and the identified displacement area. We were able to identify the location of simulated MUs in the muscle cross-section within a 2 mm error and to reconstruct the simulated MU displacement velocity (cc > 0.85). Multiple regression analysis of 180 experimental MUs detected during isometric contractions of the biceps brachii revealed a significant association between the identified location of MU displacement areas and the centroid of the EMG amplitude distribution. The proposed approach has the potential to enable non-invasive assessment of the electrical, anatomical, and mechanical properties of single MUs in voluntary contractions.

Neuromotor control associates with muscle weakness observed with McArdle sign of multiple sclerosis

OBJECTIVE

Multiple Sclerosis (MS) is often accompanied by myelopathy, which may be associated with progressive worsening. A specific finding of MS-associated myelopathy is McArdle sign, wherein neck flexion is associated with prominent increased limb weakness relative to that detected with neck extension. In this study, we characterized neuromotor control properties of finger extensors in association with the McArdle sign.

METHODS

A custom-built device was utilized to monitor torque production of the wrist extensors with simultaneous recording of surface electromyography of the extensor digitorum. The electromyography was decomposed and analyzed via both linear and nominal regressions.

RESULTS

Linear regressions demonstrated a strong difference between groups for MS from healthy controls and other myelopathies for motor unit action potential amplitude and average firing rate (p < 0.001). Further, linear regression demonstrated good correlations of neuromotor variables to mechanical torque output (0.24 ≤ R²  ≤ 0.76). Nominal regression distinguished MS from healthy controls with an AUC of 0.87, specificity of 0.97, and sensitivity of 0.64. Nominal regression of MS from other myelopathies demonstrated an AUC of 0.88, specificity of 0.85, and sensitivity of 0.79.

INTERPRETATION

These data demonstrate the neuromotor control factors that largely determine muscle force production change with the observation of McArdle sign; these neuromotor control factors can differentiate MS from both healthy controls and other myelopathy conditions.

Endurance-exercise training adaptations in spinal motoneurones: potential functional relevance to locomotor output and assessment in humans

It is clear from non-human animal work that spinal motoneurones undergo endurance training (chronic) and locomotor (acute) related changes in their electrical properties and thus their ability to fire action potentials in response to synaptic input. The functional implications of these changes, however, are speculative. In humans, data suggests that similar chronic and acute changes in motoneurone excitability may occur, though the work is limited due to technical constraints. To examine the potential influence of chronic changes in human motoneurone excitability on the acute changes that occur during locomotor output, we must develop more sophisticated recording techniques or adapt our current methods. In this review, we briefly discuss chronic and acute changes in motoneurone excitability arising from non-human and human work. We then discuss the potential interaction effects of chronic and acute changes in motoneurone excitability and the potential impact on locomotor output. Finally, we discuss the use of high-density surface electromyogram recordings to examine human motor unit firing patterns and thus, indirectly, motoneurone excitability. The assessment of single motor units from high-density recording is mainly limited to tonic motor outputs and minimally dynamic motor output such as postural sway. Adapting this technology for use during locomotor outputs would allow us to gain a better understanding of the potential functional implications of endurance training-induced changes in human motoneurone excitability on motor output.

The Effect of Resistance Training on Motor Unit Firing Properties: A Systematic Review and Meta-Analysis

While neural changes are thought to be responsible for early increases in strength following resistance training (RT), the exact changes in motor unit (MU) firing properties remain unclear. This review aims to synthesize the available evidence on the effect of RT on MU firing properties. MEDLINE (OVID interface), EMBASE (OVID interface), Web of Science (all databases), Cochrane Library, EBSCO CINAHL Plus, PubMed, and EBSCO SportDiscus were searched from inception until June 2021. Randomized controlled trials and non-randomized studies of interventions that compared RT to no intervention (control) were included. Two reviewers independently extracted data from each trial, assessed the risk of bias and rated the cumulative quality of evidence. Motor unit discharge rate (MUDR), motor unit recruitment threshold (MURT), motor unit discharge rate variability (MUDRV), MU discharge rate at recruitment vs. recruitment threshold relationship, and MU discharge rate vs. recruitment threshold relationship were assessed. Seven trials including 167 participants met the inclusion criteria. Meta-analysis (four studies) revealed that MUDR did not change significantly (P = 0.43), but with considerable heterogeneity likely to be present (I 2 = 91). Low to moderate evidence supports changes in MUDRV, MUDR at recruitment vs. recruitment threshold relationship, and the MUDR vs. recruitment threshold relationship. Overall, this systematic review revealed that there is a lack of high-quality evidence for the effect of RT on MU firing properties. Heterogeneity across studies undermines the quality of the evidence for multiple outcomes and affects the conclusions that can be drawn.

The role of the neural stimulus in regulating skeletal muscle hypertrophy

Resistance training is frequently performed with the goal of stimulating muscle hypertrophy. Due to the key roles motor unit recruitment and mechanical tension play to induce muscle growth, when programming, the manipulation of the training variables is oriented to provoke the correct stimulus. Although it is known that the nervous system is responsible for the control of motor units and active muscle force, muscle hypertrophy researchers and trainers tend to only focus on the adaptations of the musculotendinous unit and not in the nervous system behaviour. To better guide resistance exercise prescription for muscle hypertrophy and aiming to delve into the mechanisms that maximize this goal, this review provides evidence-based considerations for possible effects of neural behaviour on muscle growth when programming resistance training, and future neurophysiological measurement that should be tested when training to increase muscle mass. Combined information from the neural and muscular structures will allow to understand the exact adaptations of the muscle in response to a given input (neural drive to the muscle). Changes at different levels of the nervous system will affect the control of motor units and mechanical forces during resistance training, thus impacting the potential hypertrophic adaptations. Additionally, this article addresses how neural adaptations and fatigue accumulation that occur when resistance training may influence the hypertrophic response and propose neurophysiological assessments that may improve our understanding of resistance training variables that impact on muscular adaptations.

A computational approach for generating continuous estimates of motor unit discharge rates and visualizing population discharge characteristics

OBJECTIVE

Successive improvements in high density surface electromyography and decomposition techniques have facilitated an increasing yield in decomposed motor unit (MU) spike times. Though these advancements enhance the generalizability of findings and promote the application of MU discharge characteristics to inform the neural control of motor output, limitations remain. Specifically, 1) common approaches for generating smooth estimates of MU discharge rates introduce artifacts in quantification, which may bias findings, and 2) discharge characteristics of large MU populations are often difficult to visualize.

APPROACH

In the present study, we propose support vector regression (SVR) as an improved approach for generating smooth continuous estimates of discharge rate and compare the fit characteristics of SVR to traditionally used methods, including Hanning window filtering and polynomial regression. Furthermore, we introduce ensembles as a method to visualize the discharge characteristics of large MU populations. We define ensembles as the average discharge profile of a subpopulation of MUs, composed of a time normalized ensemble average of all units within this subpopulation. Analysis was conducted with MUs decomposed from the tibialis anterior (N = 2128), medial gastrocnemius (N = 2673), and soleus (N = 1190) during isometric plantarflexion and dorsiflexion contractions.

MAIN RESULT

Compared to traditional approaches, we found SVR to alleviate commonly observed inaccuracies and produce significantly less absolute fit error in the initial phase of MU discharge and throughout the entire duration of discharge. Additionally, we found the visualization of MU populations as ensembles to intuitively represent population discharge characteristics with appropriate accuracy for visualization.

SIGNIFICANCE

The results and methods outlined here provide an improved method for generating estimates of MU discharge rate with SVR and present a unique approach to visualizing MU populations with ensembles. In combination, the use of SVR and generation of ensembles represent an efficient method for rendering population discharge characteristics.

Skilled independent control of individual motor units via a non-invasive neuromuscular-machine interface

Objective.Brain-machine interfaces (BMIs) have the potential to augment human functions and restore independence in people with disabilities, yet a compromise between non-invasiveness and performance limits their relevance.Approach.Here, we hypothesized that a non-invasive neuromuscular-machine interface providing real-time neurofeedback of individual motor units within a muscle could enable independent motor unit control to an extent suitable for high-performance BMI applications.Main results.Over 6 days of training, eight participants progressively learned to skillfully and independently control three biceps brachii motor units to complete a 2D center-out task. We show that neurofeedback enabled motor unit activity that largely violated recruitment constraints observed during ramp-and-hold isometric contractions thought to limit individual motor unit controllability. Finally, participants demonstrated the suitability of individual motor units for powering general applications through a spelling task.Significance.These results illustrate the flexibility of the sensorimotor system and highlight individual motor units as a promising source of control for BMI applications.

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.

Identification of Motor Unit Twitch Properties in the Intact Human In Vivo

Restoring natural motor function in neurologically injured individuals is challenging, largely due to the lack of personalization in current neurorehabilitation technologies. Signal-driven neuro-musculoskeletal models may offer a novel paradigm for devising novel closed-loop rehabilitation strategies according to an individual's physiology. However, current modelling techniques are constrained to bipolar electromyography (EMG), thereby lacking the resolution necessary to extract the activity of individual motor units (MUs) in vivo. In this work, we decoded MU spike trains from high-density (HD)-EMG to obtain relevant neural properties across multiple isometric plantar-dorsiflexion tasks. Then, we sampled MU statistical distributions and used them to reproduce MU specific activation profiles. Results showed bimodal distributions which may correspond to slow and fast MU populations. The estimated activation profiles showed a high degree of similarity to the reference torque (R²>0.8) across the recorded muscles. This suggests that the estimation of MU twitch properties is a crucial step for the translation of neural information into muscle force.Clinical Relevance- This work has multiple implications for understanding the underlying mechanism of motor impairment and for developing closed-loop strategies for modulating alpha motor circuitries in neurologically injured individuals.