Factors influencing the estimates of correlation between motor unit activities in humans

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.

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

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

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.

The neural control of movement: a century of in vivo motor unit recordings is the legacy of Adrian and Bronk

In two papers dated 1928 to 1929 in The Journal of Physiology, Edgar Adrian and Detlev Bronk described recordings from motor nerve and muscle fibres. The recordings from motor nerve fibres required progressive dissection of the nerve until a few fibres remained, from which isolated single fibre activity could be detected. The muscle fibre recordings were performed in humans during voluntary contractions with an intramuscular electrode - the concentric needle electrode - that they describe for the first time in the second paper. They recognised that muscle fibres would respond to each impulse sent by the innervating motor neurone and that therefore muscle fibre recordings provided information on the times of activation of the motor nerve fibres which were as accurate as a direct record from the nerve. These observations and the description of the concentric needle electrode opened the era of motor unit recordings in humans, which have continued for almost a century and have provided a comprehensive view of the neural control of movement at the motor unit level. Despite important advances in technology, many of the principles of motor unit behaviour that would be investigated in the subsequent decades were canvassed in the two papers by Adrian and Bronk. For example, they described the concomitant motor neurones' recruitment and rate coding for force modulation, synchronisation of motor unit discharges, and the dependence of discharge rate on motor unit recruitment threshold. Here, we summarise their observations and discuss the impact of their work. We highlight the advent of the concentric needle, and its subsequent influence on motor control research.

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.

Correlation networks of spinal motor neurons that innervate lower limb muscles during a multi-joint isometric task

Movements are reportedly controlled through the combination of synergies that generate specific motor outputs by imposing an activation pattern on a group of muscles. To date, the smallest unit of analysis of these synergies has been the muscle through the measurement of its activation. However, the muscle is not the lowest neural level of movement control. In this human study (n = 10), we used a purely data-driven method grounded on graph theory to extract networks of motor neurons based on their correlated activity during an isometric multi-joint task. Specifically, high-density surface electromyography recordings from six lower limb muscles were decomposed into motor neurons spiking activity. We analysed these activities by identifying their common low-frequency components, from which networks of correlated activity to the motor neurons were derived and interpreted as networks of common synaptic inputs. The vast majority of the identified motor neurons shared common inputs with other motor neuron(s). In addition, groups of motor neurons were partly decoupled from their innervated muscle, such that motor neurons innervating the same muscle did not necessarily receive common inputs. Conversely, some motor neurons from different muscles-including distant muscles-received common inputs. The study supports the theory that movements are produced through the control of small numbers of groups of motor neurons via common inputs and that there is a partial mismatch between these groups of motor neurons and muscle anatomy. We provide a new neural framework for a deeper understanding of the structure of common inputs to motor neurons. KEY POINTS: A central and unresolved question is how spinal motor neurons are controlled to generate movement. We decoded the spiking activities of dozens of spinal motor neurons innervating six muscles during a multi-joint task, and we used a purely data-driven method grounded on graph theory to extract networks of motor neurons based on their correlated activity (considered as common input). The vast majority of the identified motor neurons shared common inputs with other motor neuron(s). Groups of motor neurons were partly decoupled from their innervated muscle, such that motor neurons innervating the same muscle did not necessarily receive common inputs. Conversely, some motor neurons from different muscles, including distant muscles, received common inputs. The study supports the theory that movement is produced through the control of groups of motor neurons via common inputs and that there is a partial mismatch between these groups of motor neurons and muscle anatomy.

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.

Common synaptic input, synergies and size principle: Control of spinal motor neurons for movement generation

Understanding how movement is controlled by the CNS remains a major challenge, with ongoing debate about basic features underlying this control. In current established views, the concepts of motor neuron recruitment order, common synaptic input to motor neurons and muscle synergies are usually addressed separately and therefore seen as independent features of motor control. In this review, we analyse the body of literature in a broader perspective and we identify a unified approach to explain apparently divergent observations at different scales of motor control. Specifically, we propose a new conceptual framework of the neural control of movement, which merges the concept of common input to motor neurons and modular control, together with the constraints imposed by recruitment order. This framework is based on the following assumptions: (1) motor neurons are grouped into functional groups (clusters) based on the common inputs they receive; (2) clusters may significantly differ from the classical definition of motor neuron pools, such that they may span across muscles and/or involve only a portion of a muscle; (3) clusters represent functional modules used by the CNS to reduce the dimensionality of the control; and (4) selective volitional control of single motor neurons within a cluster receiving common inputs cannot be achieved. Here, we discuss this framework and its underlying theoretical and experimental evidence.

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.

Handedness is associated with less common input to spinal motor neurons innervating different hand muscles

Whether the neural control of manual behaviours differs between the dominant and non-dominant hand is poorly understood. This study aimed to determine whether the level of common synaptic input to motor neurons innervating the same or different muscles differs between the dominant and the non-dominant hand. Seventeen participants performed two motor tasks with distinct mechanical requirements: an isometric pinch and an isometric rotation of a pinched dial. Each task was performed at 30% of maximum effort and was repeated with the dominant and non-dominant hand. Motor units were identified from two intrinsic (flexor digitorum interosseous and thenar) and one extrinsic muscle (flexor digitorum superficialis) from high-density surface electromyography recordings. Two complementary approaches were used to estimate common synaptic inputs. First, we calculated the coherence between groups of motor neurons from the same and from different muscles. Then, we estimated the common input for all pairs of motor neurons by correlating the low-frequency oscillations of their discharge rate. Both analyses led to the same conclusion, indicating less common synaptic input between motor neurons innervating different muscles in the dominant hand than in the non-dominant hand, which was only observed during the isometric rotation task. No between-side differences in common input were observed between motor neurons of the same muscle. This lower level of common input could confer higher flexibility in the recruitment of motor units, and therefore, in mechanical outputs. Whether this difference between the dominant and non-dominant arm is the cause or the consequence of handedness remains to be determined.

Motor Unit Discharge Patterns in Response to Focal Tendon Vibration of the Lower Limb in Cats and Humans

High-frequency vibration of the tendon provides potent activation of Ia afferents time-locked to the stimulation frequency and provides excitatory ionotropic activation of homonymous motor pools. In cats, the evoked motor unit discharge is constrained to discharge at integer multiples of the vibration frequency, resulting in a probability of discharge that is highly punctuated. Here we quantify the robustness of this punctuated response in the cat and evaluate whether it is present in the human. Soleus electromyography (EMG) was collected from eight cats using 64 channel electrodes during three modes of motoneuron activation. First, tendon vibration parameters were modified. Second, secondary reflex inputs are applied concurrently with tendon vibration. Third, the state of the spinal cord was altered through pharmacological or surgical manipulations. Analogous surface high-density EMG was collected from the lower leg of six humans during both vibration evoked and matched volitional contractions. Array EMG signals from both the cat and human were decomposed into corresponding motor unit action potential spike trains, and the punctuation in discharge was quantified. In the cat, regardless of vibration parameters, secondary synaptic drive, and state of spinal circuitry, focal tendon vibration evoked punctuated motor unit discharge. However, in the human lower limb, the vibration-evoked contractions do not produce punctuated motor unit discharge.

Online tracking of the phase difference between neural drives to antagonist muscle pairs in essential tremor patients

Transcutaneous electrical stimulation has been applied in tremor suppression applications. Out-of-phase stimulation strategies applied above or below motor threshold result in a significant attenuation of pathological tremor. For stimulation to be properly timed, the varying phase relationship between agonist-antagonist muscle activity during tremor needs to be accurately estimated in real-time. Here we propose an online tremor phase and frequency tracking technique for the customized control of electrical stimulation, based on a phase-locked loop (PLL) system applied to the estimated neural drive to muscles. Surface electromyography signals were recorded from the wrist extensor and flexor muscle groups of 13 essential tremor patients during postural tremor. The EMG signals were pre-processed and decomposed online and offline via the convolution kernel compensation algorithm to discriminate motor unit spike trains. The summation of motor unit spike trains detected for each muscle was bandpass filtered between 3 to 10 Hz to isolate the tremor related components of the neural drive to muscles. The estimated tremorogenic neural drive was used as input to a PLL that tracked the phase differences between the two muscle groups. The online estimated phase difference was compared with the phase calculated offline using a Hilbert Transform as a ground truth. The results showed a rate of agreement of 0.88 ± 0.22 between offline and online EMG decomposition. The PLL tracked the phase difference of tremor signals in real-time with an average correlation of 0.86 ± 0.16 with the ground truth (average error of 6.40° ± 3.49°). Finally, the online decomposition and phase estimation components were integrated with an electrical stimulator and applied in closed-loop on one patient, to representatively demonstrate the working principle of the full tremor suppression system. The results of this study support the feasibility of real-time estimation of the phase of tremorogenic neural drive to muscles, providing a methodology for future tremor-suppression neuroprostheses.

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.

Influence of central and peripheral motor unit properties on isometric muscle force entropy: A computer simulation study

Approximate entropy of isometric force is a popular measure to characterize behavioral changes across muscle contraction conditions. The degree to which force entropy characterizes the randomness of the motor control strategy, however, is not known. In this study, we used a computational model to investigate the correlation between approximate entropy of the synaptic input to a motor neuron pool, the neural drive to muscle (cumulative spike train; CST), and the force. This comparison was made across several simulation conditions, that included different synaptic command signal bandwidths, motor neuron pool sizes, and muscle contractile properties. The results indicated that although force entropy to some degree reflects the entropy of the synaptic command to motor neurons, it is biased by changes in motor unit properties. As a consequence, there was a low correlation between approximate entropy of force and the motor neuron input signal across all simulation conditions (r² = 0.13). Therefore, force entropy should only be used to compare motor control strategies across conditions where motor neuron properties can be assumed to be maintained. Instead, we recommend that the entropy of the descending motor commands should be estimated from CSTs comprising spike trains of multiple motor units.

Feedforward modulation of gamma motor neuron activity can improve motor command accuracy

Objective. Coactivation of gamma and alpha motor neuron activity ensures that muscle spindle responsiveness is maintained during muscle contractions. However, some evidence suggests that the activity of gamma motor neurons is phase-advanced with respect to that of alpha motor neurons during manual control tasks. We hypothesized that this might be a deliberate control strategy to maximize movement accuracy.Approach. Using a computational model of the neural activation of a muscle and its type Ia sensory feedback to the motor neurons, we systematically investigated the impact of the phase difference between oscillatory descending input to alpha and dynamic gamma motor neurons. Specifically, the amplification of the alpha motor neuron drive to the muscle was investigated as a function of the frequency of the synaptic input (1-9 Hz individually or superimposed) and the alpha-gamma phase difference (0-2π).Main results. Simulation results showed that when the phase advance of the dynamic gamma drive resulted in delays between muscle velocity and type Ia afferent feedback similar to those previously observed experimentally, low-frequency components (1 and 2 Hz) of the motor neuron synaptic input were amplified (gain up to 1.7). On the other hand, synaptic input at higher frequencies was unaffected.Significance. This finding suggests that by imposing a phase advance of the input to dynamic gamma motor neurons, components of the neural drive usually associated with voluntary control are amplified. In this way, our study suggests that this neural strategy increases the control-to-neural-noise ratio of the motor output during movement.

Individual differences in the neural strategies to control the lateral and medial head of the quadriceps during a mechanically constrained task

We observed that the distribution of the strength of neural drive between the vastus lateralis and vastus medialis during a single-joint isometric task varied across participants. Also, we observed that the proportion of neural drive that was shared within and between these muscles also varied across participants. These results provide evidence that the neural strategies to control the vastus lateralis and vastus medialis muscles widely vary across individuals, even during a mechanically constrained task.

Interfacing With Alpha Motor Neurons in Spinal Cord Injury Patients Receiving Trans-spinal Electrical Stimulation

Trans-spinal direct current stimulation (tsDCS) provides a non-invasive, clinically viable approach to potentially restore physiological neuromuscular function after neurological impairment, e.g., spinal cord injury (SCI). Use of tsDCS has been hampered by the inability of delivering stimulation patterns based on the activity of neural targets responsible to motor function, i.e., α-motor neurons (α-MNs). State of the art modeling and experimental techniques do not provide information about how individual α-MNs respond to electrical fields. This is a major element hindering the development of neuro-modulative technologies highly tailored to an individual patient. For the first time, we propose the use of a signal-based approach to infer tsDCS effects on large α-MNs pools in four incomplete SCI individuals. We employ leg muscles spatial sampling and deconvolution of high-density fiber electrical activity to decode accurate α-MNs discharges across multiple lumbosacral segments during isometric plantar flexion sub-maximal contractions. This is done before, immediately after and 30 min after sub-threshold cathodal stimulation. We deliver sham tsDCS as a control measure. First, we propose a new algorithm for removing compromised information from decomposed α-MNs spike trains, thereby enabling robust decomposition and frequency-domain analysis. Second, we propose the analysis of α-MNs spike trains coherence (i.e., frequency-domain) as an indicator of spinal response to tsDCS. Results showed that α-MNs spike trains coherence analysis sensibly varied across stimulation phases. Coherence analyses results suggested that the common synaptic input to α-MNs pools decreased immediately after cathodal tsDCS with a persistent effect after 30 min. Our proposed non-invasive decoding of individual α-MNs behavior may open up new avenues for the design of real-time closed-loop control applications including both transcutaneous and epidural spinal electrical stimulation where stimulation parameters are adjusted on-the-fly.

Fatigue-related modulation of low-frequency common drive to motor units

PURPOSE

This study investigated fatigue-related modulation of common neural inputs to motor units (MUs) under 5 Hz, which determines force precision control.

METHODS

Twenty-seven adults performed a sequence of fatiguing contractions. The participants were assessed with a static isometric index abduction at 20% maximal voluntary contraction in the pre-test and post-test. Discharge characteristics of MUs of the first dorsal interosseous muscle were analyzed with decomposed EMG signals.

RESULTS

Along with increases in the mean (58.40 ± 11.76 ms → 62.55 ± 10.83 ms, P = 0.029) and coefficient of variation (0.204 ± .014 → 0.215 ± 0.017, P = 0.002) in inter-spike intervals, the fatiguing contraction caused reductions in the mean frequency (16.84 ± 3.31 Hz → 15.59 ± 3.21 Hz, P = 0.027) and spectral dispersions (67.54 ± 4.49 → 62.64 ± 6.76 Hz, P = 0.007) of common neural drive, as estimated with smoothed cumulative motor unit spike trains (SCMUSTs). Stabilogram diffusion analysis of SCMUSTs revealed significant fatigue-related reductions in the long-term effective diffusion coefficient (1.91 ± 0.77 Hz²/s → 1.61 ± 0.61 Hz²/s, P = 0.020) and long-term scaling exponent (0.480 ± 0.013 Hz²/s → 0.471 ± 0.017 Hz²/s, P = 0.014). After fatiguing contraction, mutual information of force fluctuations and SCMUSTs was augmented roughly by 12.95% (P = 0.041).

CONCLUSIONS

Muscular fatigue could compress and shift the low-frequency common drive to MUs toward lower spectral bands, thereby enhancing transmission of twitch forces through the muscle-tendon complex with a low-pass filter property. The fatigue-induced changes involve increased closed-loop control of the common modulation of MU discharge rates.

Neuro-Musculoskeletal Mapping for Man-Machine Interfacing

We propose a myoelectric control method based on neural data regression and musculoskeletal modeling. This paradigm uses the timings of motor neuron discharges decoded by high-density surface electromyogram (HD-EMG) decomposition to estimate muscle excitations. The muscle excitations are then mapped into the kinematics of the wrist joint using forward dynamics. The offline tracking performance of the proposed method was superior to that of state-of-the-art myoelectric regression methods based on artificial neural networks in two amputees and in four out of six intact-bodied subjects. In addition to joint kinematics, the proposed data-driven model-based approach also estimated several biomechanical variables in a full feed-forward manner that could potentially be useful in supporting the rehabilitation and training process. These results indicate that using a full forward dynamics musculoskeletal model directly driven by motor neuron activity is a promising approach in rehabilitation and prosthetics to model the series of transformations from muscle excitation to resulting joint function.

Exogenous neuromodulation of spinal neurons induces beta-band coherence during self-sustained discharge of hind limb motor unit populations

The spontaneous or self-sustained discharge of spinal motoneurons can be observed in both animals and humans. Although the origins of this self-sustained discharge are not fully known, it can be generated by activation of persistent inward currents intrinsic to the motoneuron. If self-sustained discharge is generated exclusively through this intrinsic mechanism, the discharge of individual motor units will be relatively independent of one another. Alternatively, if increased activation of premotor circuits underlies this prolonged discharge of spinal motoneurons, we would expect correlated activity among motoneurons. Our aim is to assess potential synaptic drive by quantifying coherence during self-sustained discharge of spinal motoneurons. Electromyographic activity was collected from 20 decerebrate animals using a 64-channel electrode grid placed on the isolated soleus muscle before and following intrathecal administration of methoxamine, a selective α₁-noradrenergic agonist. Sustained muscle activity was recorded and decomposed into the discharge times of ~10-30 concurrently active individual motor units. Consistent with previous reports, the self-sustained discharge of motor units occurred at low mean discharge rates with low-interspike variability. Before methoxamine administration, significant low-frequency coherence (<2 Hz) was observed, while minimal coherence was observed within higher frequency bands. Following intrathecal administration of methoxamine, increases in motor unit discharge rates and strong coherence in both the low-frequency and 15- to 30-Hz beta bands were observed. These data demonstrate beta-band coherence among motor units can be observed through noncortical mechanisms and that neuromodulation of spinal/brainstem neurons greatly influences coherent discharge within spinal motor pools.NEW & NOTEWORTHY The correlated discharge of spinal motoneurons is often used to describe the input to the motor pool. We demonstrate spinal/brainstem neurons devoid of cortical input can generate correlated motor unit discharge in the 15- to 30-Hz beta band, which is amplified through neuromodulation. Activity in the beta band is often ascribed to cortical drive in humans; however, these data demonstrate the capability of the mammalian segmental motor system to generate and modulate this coherent state of motor unit discharge.

The effect of visuospatial resolution on discharge variability among motor units and force-discharge relation

Although force steadiness varies with visuospatial information, accountable motor unit (MU) behaviors are not fully understood. This study investigated the modulation of MU discharges and force-discharge relation due to variations in the spatial resolution of visual feedback, with a particular focus on discharge variability among MUs. Fourteen young adults produced isometric force at 10% of maximal voluntary contraction (MVC) through index abduction, under the conditions of force trajectory displayed with low visual gain (LVG) and high visual gain (HVG). Together with smaller and more complex force fluctuations, HVG resulted in greater variabilities of the mean interspike interval and discharge irregularity among MUs than LVG did. Estimated via smoothening of a cumulative spike train of all MUs, global discharge rate was tuned to visual gain, with a more complex global discharge rate and a lower force-discharge relation in the HVG condition. These higher discharge variabilities were linked to larger variance of the common drive received by MUs for regulation of muscle force with higher visuospatial information. In summary, higher visuospatial information improves force steadiness with more complex force fluctuations, underlying joint effects of low-pass filter property of the musculotendon complex and central modulation of discharge variability among MUs.

Voluntary and tremorogenic inputs to motor neuron pools of agonist/antagonist muscles in essential tremor patients

Pathological tremor is an oscillation of body parts at 3–10 Hz, determined by the output of spinal motor neurons (MNs), which receive synaptic inputs from supraspinal centers and muscle afferents. The behavior of spinal MNs during tremor is not well understood, especially in relation to the activation of the multiple muscles involved. Recent studies on patients with essential tremor have shown that antagonist MN pools receive shared input at the tremor frequency. In this study, we investigated the synaptic inputs related to tremor and voluntary movement, and their coordination across antagonist muscles. We analyzed the spike trains of motor units (MUs) identified from high-density surface electromyography from the forearm extensor and flexor muscles in 15 patients with essential tremor during postural tremor. The shared synaptic input was quantified by coherence and phase difference analysis of the spike trains. All pairs of spike trains in each muscle showed coherence peaks at the voluntary drive frequency (1–3 Hz, 0.2 ± 0.2, mean ± SD) and tremor frequency (3–10 Hz, 0.6 ± 0.3) and were synchronized with small phase differences (3.3 ± 25.2° and 3.9 ± 22.0° for the voluntary drive and tremor frequencies, respectively). The coherence between MN spike trains of antagonist muscle groups at the tremor frequency was significantly smaller than intramuscular coherence. We predominantly observed in-phase activation of MUs between agonist/antagonist muscles at the voluntary frequency band (0.6 ± 48.8°) and out-of-phase activation at the tremor frequency band (126.9 ± 75.6°). Thus MNs innervating agonist/antagonist muscles concurrently receive synaptic inputs with different phase shifts in the voluntary and tremor frequency bands.

NEW & NOTEWORTHY Although the mechanical characteristics of tremor have been widely studied, the activation of the affected muscles is still poorly understood. We analyzed the behavior of motor units of pairs of antagonistic wrist muscle groups in patients with essential tremor and studied their activity at voluntary movement- and tremor-related frequencies. We found that the phase relation between inputs to antagonistic muscles is different at the voluntary and tremor frequency bands.

The relative strength of common synaptic input to motor neurons is not a determinant of the maximal rate of force development in humans

Correlation between motor unit discharge times, often referred to as motor unit synchronization, is determined by common synaptic input to motor neurons. Although it has been largely speculated that synchronization should influence the rate of force development, the association between the degree of motor unit synchronization and rapid force generation has not been determined. In this study, we examined this association with both simulations and experimental motor unit recordings. The analysis of experimental motor unit discharges from the tibialis anterior muscle of 20 healthy individuals during rapid isometric contractions revealed that the average motor unit discharge rate was associated with the rate of force development. Moreover, the extent of motor unit synchronization was entirely determined by the average motor unit discharge rate ( R > 0.7, P < 0.0001). The simulation model demonstrated that the relative proportion of common synaptic input received by motor neurons, which determines motor unit synchronization, does not influence the rate of force development ( R = 0.03, P > 0.05). Nonetheless, the estimates of correlation between motor unit spike trains were significantly correlated with the rate of force generation ( R > 0.8, P < 0.0001). These results indicate that the average motor unit discharge rate, but not the degree of motor unit synchronization, contributes to most of the variance of human contractile speed among individuals. In addition, estimates of correlation between motor unit discharge times depend strongly on the number of identified motor units and therefore are not indicative of the strength of common input.

NEW & NOTEWORTHY It is commonly assumed that motor unit synchronization has an impact on the rate of force development of a muscle. Here we present computer simulations and experimental data of human tibialis anterior motor units during rapid contractions that show that motor unit synchronization is not a determinant of the rate of force production. This conclusion clarifies the neural determinants of rapid force generation.

Multiscale modeling of the neuromuscular system: Coupling neurophysiology and skeletal muscle mechanics

Mathematical models and computer simulations have the great potential to substantially increase our understanding of the biophysical behavior of the neuromuscular system. This, however, requires detailed multiscale, and multiphysics models. Once validated, such models allow systematic in silico investigations that are not necessarily feasible within experiments and, therefore, have the ability to provide valuable insights into the complex interrelations within the healthy system and for pathological conditions. Most of the existing models focus on individual parts of the neuromuscular system and do not consider the neuromuscular system as an integrated physiological system. Hence, the aim of this advanced review is to facilitate the prospective development of detailed biophysical models of the entire neuromuscular system. For this purpose, this review is subdivided into three parts. The first part introduces the key anatomical and physiological aspects of the healthy neuromuscular system necessary for modeling the neuromuscular system. The second part provides an overview on state-of-the-art modeling approaches representing all major components of the neuromuscular system on different time and length scales. Within the last part, a specific multiscale neuromuscular system model is introduced. The integrated system model combines existing models of the motor neuron pool, of the sensory system and of a multiscale model describing the mechanical behavior of skeletal muscles. Since many sub-models are based on strictly biophysical modeling approaches, it closely represents the underlying physiological system and thus could be employed as starting point for further improvements and future developments. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

Amplitude cancellation influences the association between frequency components in the neural drive to muscle and the rectified EMG signal

The rectified surface EMG signal is commonly used as an estimator of the neural drive to muscles and therefore to infer sources of synaptic input to motor neurons. Loss of EMG amplitude due to the overlap of motor unit action potentials (amplitude cancellation), however, may distort the spectrum of the rectified EMG and thereby its correlation with the neural drive. In this study, we investigated the impact of amplitude cancelation on this correlation using analytical derivations and a computational model of motor neuron activity, force, and the EMG signal. First, we demonstrated analytically that an ideal rectified EMG signal without amplitude cancellation (EMG_nc) is superior to the actual rectified EMG signal as estimator of the neural drive to muscle. This observation was confirmed by the simulations, as the average coefficient of determination (r²) between the neural drive in the 1–30 Hz band and EMG_nc (0.59±0.08) was matched by the correlation between the rectified EMG and the neural drive only when the level of amplitude cancellation was low (<40%) at low contraction levels (<5% of maximum voluntary contraction force; MVC). This correlation, however, decreased linearly with amplitude cancellation (r = -0.83) to values of r² <0.2 at amplitude cancellation levels >60% (contraction levels >15% MVC). Moreover, the simulations showed that a stronger (i.e. more variable) neural drive implied a stronger correlation between the rectified EMG and the neural drive and that amplitude cancellation distorted this correlation mainly for low-frequency components (<5 Hz) of the neural drive. In conclusion, the results indicate that amplitude cancellation distorts the spectrum of the rectified EMG signal. This implies that valid use of the rectified EMG as an estimator of the neural drive requires low contraction levels and/or strong common synaptic input to the motor neurons.

The rectified surface EMG signal is commonly used to analyze the neural activation of muscles. However, since this signal is most often exposed to so-called amplitude cancellation (loss of EMG amplitude due to overlap of positive and negative phases of different motor unit action potentials), the frequency content of the rectified EMG may not fully reflect that of the neural drive to the muscle. In this study we prove this notion analytically and demonstrate, using simulations, that the rectified EMG signal accurately reflects the neural drive to the muscle only in a limited set of conditions. Specifically, these conditions include low contraction levels and/or high variability of the neural drive. In other conditions, the rectified EMG signal from a muscle is a poor predictor of its neural input. This finding has potentially large implications for the way neural drive to muscles and neural connectivity (e.g. across muscles or between the brain and a muscle) should be analyzed.

Coherence of the Surface EMG and Common Synaptic Input to Motor Neurons

Coherence between electromyographic (EMG) signals is often used to infer the common synaptic input to populations of motor neurons. This analysis, however, may be limited due to the filtering effect of the motor unit action potential waveforms. This study investigated the ability of surface EMG–EMG coherence to predict common synaptic input to motor neurons. Surface and intramuscular EMG were recorded from two locations of the tibialis anterior muscle during steady ankle dorsiflexions at 5 and 10% of the maximal force in 10 healthy individuals. The intramuscular EMG signals were decomposed to identify single motor unit spike trains. For each trial, the strength of the common input in different frequency bands was estimated from the coherence between two cumulative spike trains, generated from sets of single motor unit spike trains (reference measure). These coherence values were compared with those obtained from the coherence between the surface EMG signals (raw, rectified, and high-passed filtered at 250 Hz before rectification) using linear regression. Overall, the high-pass filtering of the EMG prior to rectification did not substantially change the results with respect to rectification only. For both signals, the correlation of EMG coherence with motor unit coherence was strong at 5% MVC (r² > 0.8; p < 0.01), but only for frequencies > 5 Hz. At 10% MVC, the correlation between EMG and motor unit coherence was only significant for frequencies > 15 Hz (r² > 0.8; p < 0.01). However, when using raw EMG for coherence analysis, the only significant relation with motor unit coherence was observed for the bandwidth 5–15 Hz (r² > 0.68; p = 0.04). In all cases, there was no association between motor unit and EMG coherence for frequencies < 5 Hz (r² ≤ 0.2; p ≥ 0.51). In addition, a substantial error in the best linear fit between motor unit and EMG coherence was always present. In conclusion, high-frequency (>5 Hz) common synaptic inputs to motor neurons can partly be estimated from the rectified surface EMG at low-level steady contractions. The results, however, suggest that this association is weakened with increasing contraction intensity and that input at lower frequencies during steady isometric contractions cannot be detected accurately by surface EMG coherence.

Variations in Static Force Control and Motor Unit Behavior with Error Amplification Feedback in the Elderly

Error amplification (EA) feedback is a promising approach to advance visuomotor skill. As error detection and visuomotor processing at short time scales decline with age, this study examined whether older adults could benefit from EA feedback that included higher-frequency information to guide a force-tracking task. Fourteen young and 14 older adults performed low-level static isometric force-tracking with visual guidance of typical visual feedback and EA feedback containing augmented high-frequency errors. Stabilogram diffusion analysis was used to characterize force fluctuation dynamics. Also, the discharge behaviors of motor units and pooled motor unit coherence were assessed following the decomposition of multi-channel surface electromyography (EMG). EA produced different behavioral and neurophysiological impacts on young and older adults. Older adults exhibited inferior task accuracy with EA feedback than with typical visual feedback, but not young adults. Although stabilogram diffusion analysis revealed that EA led to a significant decrease in critical time points for both groups, EA potentiated the critical point of force fluctuations [Formula: see text], short-term effective diffusion coefficients (Ds), and short-term exponent scaling only for the older adults. Moreover, in older adults, EA added to the size of discharge variability of motor units and discharge regularity of cumulative discharge rate, but suppressed the pooled motor unit coherence in the 13-35 Hz band. Virtual EA alters the strategic balance between open-loop and closed-loop controls for force-tracking. Contrary to expectations, the prevailing use of closed-loop control with EA that contained high-frequency error information enhanced the motor unit discharge variability and undermined the force steadiness in the older group, concerning declines in physiological complexity in the neurobehavioral system and the common drive to the motoneuronal pool against force destabilization.

Origins of Common Neural Inputs to Different Compartments of the Extensor Digitorum Communis Muscle

The extensor digitorum communis (EDC) is a multi-compartment muscle that allows dexterous extension of the four digits. However, the level of common input shared across different compartments of this muscle is not well understood. We seek to systematically characterize the common and independent neural input, originated from different levels of the central nervous system, to the different compartments. A motor unit (MU) coherence analysis was used to capture the different sources of common and independent input, by quantifying the coherence of MU discharge between different compartments. The MU activities were obtained from decomposition of surface electromyogram recordings. Our results showed that the MU coherence across different muscle compartments accounted for only a small proportion (<20%) of the total input in the alpha (5–12 Hz) and beta (15–30 Hz) bands, but was a major driver (>60%) in the delta (1–4 Hz) band. Additionally, cross-compartment coherence between the middle and ring-little fingers tended to be higher as compared with other finger combinations. Overall, the common input shared across different fingers are found to be at low to moderate levels, in comparison with the total input, which allows dexterous control of individual digits with some degree of coordinated control of multiple digits.

Robust estimation of average twitch contraction forces of populations of motor units in humans

The characteristics of motor unit force twitch profiles provide important information for the understanding of the muscle force generation. The twitch force is commonly estimated with the spike-triggered averaging technique, which, despite the many limitations, has been important for clarifying central issues in force generation. In this study, we propose a new technique for the estimation of the average twitch profile of populations of motor units with uniform contractile properties. The method encompasses a model-based deconvolution of the force signal using the identified discharge times of a population of motor units. The proposed technique was validated using simulations and tested on signals recorded during voluntary activation. The results of the simulations showed that the proposed method provides accurate estimates (relative error <25%) of the main parameters of the average twitch force when the number of identified motor units is between 5% and 15% of the total number of active motor units. It is discussed that current detection and decomposition methods of multi-channel surface EMG signals allow decoding this relative sample of the active motor unit pool. However, even when this condition is not met, our results show that the estimates provided by the new method are anyway always superior to those obtained by the spike triggered average approach, especially for high motor unit synchronization levels and when a relatively small number of triggers is available. In conclusion, we present a new method that overcome the main limitations of the spike-triggered average for the study of contractile properties of individual motor units. The method provides a new reliable tool for the investigation of the determinants of muscle force.

Altered Motor Unit Discharge Coherence in Paretic Muscles of Stroke Survivors

After a cerebral stroke, a series of changes at the supraspinal and spinal nervous system can alter the control of muscle activation, leading to persistent motor impairment. However, the relative contribution of these different levels of the nervous system to impaired muscle activation is not well understood. The coherence of motor unit (MU) spike trains is considered to partly reflect activities of higher level control, with different frequency band representing different levels of control. Accordingly, the objective of this study was to quantify the different sources of contribution to altered muscle activation. We examined the coherence of MU spike trains decomposed from surface electromyogram (sEMG) of the first dorsal interosseous muscle on both paretic and contralateral sides of 14 hemispheric stroke survivors. sEMG was obtained over a range of force contraction levels at 40, 50, and 60% of maximum voluntary contraction. Our results showed that MU coherence increased significantly in delta (1–4 Hz), alpha (8–12 Hz), and beta (15–30 Hz) bands on the affected side compared with the contralateral side, but was maintained at the same level in the gamma (30–60 Hz) band. In addition, no significant alteration was observed across medium–high force levels (40–60%). These results indicated that the common synaptic input to motor neurons increased on the paretic side, and the increased common input can originate from changes at multiple levels, including spinal and supraspinal levels following a stroke. All these changes can contribute to impaired activation of affected muscles in stroke survivors. Our findings also provide evidence regarding the different origins of impaired muscle activation poststroke.

Nonlinear Frequency-Domain Analysis of the Transformation of Cortical Inputs by a Motoneuron Pool-Muscle Complex

Corticomotor coherence in the beta and/or gamma bands has been described in different motor tasks, but the role of descending brain oscillations on force control has been elusive. Large-scale computational models of a motoneuron pool and the muscle it innervates have been used as tools to advance the knowledge of how neural elements may influence force control. Here, we present a frequency domain analysis of a NARX model fitted to a large-scale neuromuscular model by the means of generalized frequency response functions (GFRF). The results of such procedures indicated that the computational neuromuscular model was capable of transforming an oscillatory synaptic input (e.g., at 20 Hz) into a constant mean muscle force output. The nonlinearity uncovered by the GFRFs of the NARX model was responsible for the demodulation of an oscillatory input (e.g., a beta band oscillation coming from the brain and forming the input to the motoneuron pool). This suggests a manner by which brain rhythms descending as command signals to the spinal cord and acting on a motoneuron pool can regulate a maintained muscle force. In addition to the scientific aspects of these results, they provide new interpretations that may further neural engineering applications associated with quantitative neurological diagnoses and robotic systems for artificial limbs.

Paradigm Shifts in Voluntary Force Control and Motor Unit Behaviors with the Manipulated Size of Visual Error Perception

The detection of error information is an essential prerequisite of a feedback-based movement. This study investigated the differential behavior and neurophysiological mechanisms of a cyclic force-tracking task using error-reducing and error-enhancing feedback. The discharge patterns of a relatively large number of motor units (MUs) were assessed with custom-designed multi-channel surface electromyography following mathematical decomposition of the experimentally-measured signals. Force characteristics, force-discharge relation, and phase-locking cortical activities in the contralateral motor cortex to individual MUs were contrasted among the low (LSF), normal (NSF), and high scaling factor (HSF) conditions, in which the sizes of online execution errors were displayed with various amplification ratios. Along with a spectral shift of the force output toward a lower band, force output with a more phase-lead became less irregular, and tracking accuracy was worse in the LSF condition than in the HSF condition. The coherent discharge of high phasic (HP) MUs with the target signal was greater, and inter-spike intervals were larger, in the LSF condition than in the HSF condition. Force-tracking in the LSF condition manifested with stronger phase-locked EEG activity in the contralateral motor cortex to discharge of the (HP) MUs (LSF > NSF, HSF). The coherent discharge of the (HP) MUs during the cyclic force-tracking predominated the force-discharge relation, which increased inversely to the error scaling factor. In conclusion, the size of visualized error gates motor unit discharge, force-discharge relation, and the relative influences of the feedback and feedforward processes on force control. A smaller visualized error size favors voluntary force control using a feedforward process, in relation to a selective central modulation that enhance the coherent discharge of (HP) MUs.

The potential for understanding the synaptic organization of human motor commands via the firing patterns of motoneurons

Motoneurons are unique in being the only neurons in the CNS whose firing patterns can be easily recorded in human subjects. This is because of the one-to-one relationship between the motoneuron and muscle cell behavior. It has long been appreciated that the connection of motoneurons to their muscle fibers allows their action potentials to be amplified and recorded, but only recently has it become possible to simultaneously record the firing pattern of many motoneurons via array electrodes placed on the skin. These firing patterns contain detailed information about the synaptic organization of motor commands to the motoneurons. This review focuses on parameters in these firing patterns that are directly linked to specific features of this organization. It is now well established that motor commands consist of three components, excitation, inhibition, and neuromodulation; the importance of the third component has become increasingly evident. Firing parameters linked to each of the three components are discussed, along with consideration of potential limitations in their utility for understanding the underlying organization of motor commands. Future work based on realistic computer simulations of motoneurons may allow quantitative "reverse engineering" of human motoneuron firing patterns to provide good estimates of the relative amplitudes and temporal patterns of all three components of motor commands.

Alterations in Neural Control of Constant Isometric Contraction with the Size of Error Feedback

Discharge patterns from a population of motor units (MUs) were estimated with multi-channel surface electromyogram and signal processing techniques to investigate parametric differences in low-frequency force fluctuations, MU discharges, and force-discharge relation during static force-tracking with varying sizes of execution error presented via visual feedback. Fourteen healthy adults produced isometric force at 10% of maximal voluntary contraction through index abduction under three visual conditions that scaled execution errors with different amplification factors. Error-augmentation feedback that used a high amplification factor (HAF) to potentiate visualized error size resulted in higher sample entropy, mean frequency, ratio of high-frequency components, and spectral dispersion of force fluctuations than those of error-reducing feedback using a low amplification factor (LAF). In the HAF condition, MUs with relatively high recruitment thresholds in the dorsal interosseous muscle exhibited a larger coefficient of variation for inter-spike intervals and a greater spectral peak of the pooled MU coherence at 13-35 Hz than did those in the LAF condition. Manipulation of the size of error feedback altered the force-discharge relation, which was characterized with non-linear approaches such as mutual information and cross sample entropy. The association of force fluctuations and global discharge trace decreased with increasing error amplification factor. Our findings provide direct neurophysiological evidence that favors motor training using error-augmentation feedback. Amplification of the visualized error size of visual feedback could enrich force gradation strategies during static force-tracking, pertaining to selective increases in the discharge variability of higher-threshold MUs that receive greater common oscillatory inputs in the β-band.

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

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

Theoretical Model and Experimental Validation of the estimated proportions of common and independent input to motor neurons

Motor neurons in the spinal cord receive synaptic input that comprises common and independent components. The part of synaptic input that is common to all motor neurons is the one regulating the production of force. Therefore, its quantification is important to assess the strategy used by Central Nervous System (CNS) to control and regulate movements, especially in physiological conditions such as fatigue. In this study we present and validate a method to estimate the ratio between strengths of common and independent inputs to motor neurons and we apply this method to investigate its changes during fatigue. By means of coherence analysis we estimated the level of correlation between motor unit spike trains at the beginning and at the end of fatiguing contractions of the Tibialis Anterior muscle at three different force targets. Combining theoretical modeling and experimental data we estimated the strength of the common synaptic input with respect to the independent one. We observed a consistent increase in the proportion of the shared input to motor neurons during fatigue. This may be interpreted as a strategy used by the CNS to counteract the occurrence of fatigue and the concurrent decrease of generated force.

Accurate and representative decoding of the neural drive to muscles in humans with multi-channel intramuscular thin-film electrodes

Intramuscular electrodes developed over the past 80 years can record the concurrent activity of only a few motor units active during a muscle contraction. We designed, produced and tested a novel multi-channel intramuscular wire electrode that allows in vivo concurrent recordings of a substantially greater number of motor units than with conventional methods. The electrode has been extensively tested in deep and superficial human muscles. The performed tests indicate the applicability of the proposed technology in a variety of conditions. The electrode represents an important novel technology that opens new avenues in the study of the neural control of muscles in humans. We describe the design, fabrication and testing of a novel multi-channel thin-film electrode for detection of the output of motoneurones in vivo and in humans, through muscle signals. The structure includes a linear array of 16 detection sites that can sample intramuscular electromyographic activity from the entire muscle cross-section. The structure was tested in two superficial muscles (the abductor digiti minimi (ADM) and the tibialis anterior (TA)) and a deep muscle (the genioglossus (GG)) during contractions at various forces. Moreover, surface electromyogram (EMG) signals were concurrently detected from the TA muscle with a grid of 64 electrodes. Surface and intramuscular signals were decomposed into the constituent motor unit (MU) action potential trains. With the intramuscular electrode, up to 31 MUs were identified from the ADM muscle during an isometric contraction at 15% of the maximal force (MVC) and 50 MUs were identified for a 30% MVC contraction of TA. The new electrode detects different sources from a surface EMG system, as only one MU spike train was found to be common in the decomposition of the intramuscular and surface signals acquired from the TA. The system also allowed access to the GG muscle, which cannot be analysed with surface EMG, with successful identification of MU activity. With respect to classic detection systems, the presented thin-film structure enables recording from large populations of active MUs of deep and superficial muscles and thus can provide a faithful representation of the neural drive sent to a muscle.

The human motor neuron pools receive a dominant slow-varying common synaptic input

KEY POINTS

Motor neurons in a pool receive both common and independent synaptic inputs, although the proportion and role of their common synaptic input is debated. Classic correlation techniques between motor unit spike trains do not measure the absolute proportion of common input and have limitations as a result of the non-linearity of motor neurons. We propose a method that for the first time allows an accurate quantification of the absolute proportion of low frequency common synaptic input (<5 Hz) to motor neurons in humans. We applied the proposed method to three human muscles and determined experimentally that they receive a similar large amount (>60%) of common input, irrespective of their different functional and control properties. These results increase our knowledge about the role of common and independent input to motor neurons in force control.

ABSTRACT

Motor neurons receive both common and independent synaptic inputs. This observation is classically based on the presence of a significant correlation between pairs of motor unit spike trains. The functional significance of different relative proportions of common input across muscles, individuals and conditions is still debated. One of the limitations in our understanding of correlated input to motor neurons is that it has not been possible so far to quantify the absolute proportion of common input with respect to the total synaptic input received by the motor neurons. Indeed, correlation measures of pairs of output spike trains only allow for relative comparisons. In the present study, we report for the first time an approach for measuring the proportion of common input in the low frequency bandwidth (<5 Hz) to a motor neuron pool in humans. This estimate is based on a phenomenological model and the theoretical fitting of the experimental values of coherence between the permutations of groups of motor unit spike trains. We demonstrate the validity of this theoretical estimate with several simulations. Moreover, we applied this method to three human muscles: the abductor digiti minimi, tibialis anterior and vastus medialis. Despite these muscles having different functional roles and control properties, as confirmed by the results of the present study, we estimate that their motor pools receive a similar and large (>60%) proportion of common low frequency oscillations with respect to their total synaptic input. These results suggest that the central nervous system provides a large amount of common input to motor neuron pools, in a similar way to that for muscles with different functional and control properties.

Preferential distribution of nociceptive input to motoneurons with muscle units in the cranial portion of the upper trapezius muscle

Pain is associated with changes in the neural drive to muscles. For the upper trapezius muscle, surface electromyography (EMG) recordings have indicated that acute noxious stimulation in either the cranial or the caudal region of the muscle leads to a relative decrease in muscle activity in the cranial region. It is, however, not known if this adaption reflects different recruitment thresholds of the upper trapezius motor units in the cranial and caudal region or a nonuniform nociceptive input to the motor units of both regions. This study investigated these potential mechanisms by direct motor unit identification. Motor unit activity was investigated with high-density surface EMG signals recorded from the upper trapezius muscle of 12 healthy volunteers during baseline, control (intramuscular injection of isotonic saline), and painful (hypertonic saline) conditions. The EMG was decomposed into individual motor unit spike trains. Motor unit discharge rates decreased significantly from control to pain conditions by 4.0 ± 3.6 pulses/s (pps) in the cranial region but not in the caudal region (1.4 ± 2.8 pps; not significant). These changes were compatible with variations in the synaptic input to the motoneurons of the two regions. These adjustments were observed, irrespective of the location of noxious stimulation. These results strongly indicate that the nociceptive synaptic input is distributed in a nonuniform way across regions of the upper trapezius muscle.

Principles of Motor Unit Physiology Evolve With Advances in Technology

Movements are generated by the coordinated activation of motor units. Recent technological advances have made it possible to identify the concurrent activity of several tens of motor units, in contrast with much smaller samples available in classic studies. We discuss how these advances in technology have enabled the development of a population perspective of how the central nervous system controls motor unit activity and thereby the forces exerted by muscles.

Fast Oscillatory Commands from the Motor Cortex Can Be Decoded by the Spinal Cord for Force Control

Oscillations in the beta and gamma bands (13-30 Hz; 35-70 Hz) have often been observed in motor cortical outputs that reach the spinal cord, acting on motoneurons and interneurons. However, the frequencies of these oscillations are above the muscle force frequency range. A current view is that the transformation of the motoneuron pool inputs into force is linear. For this reason possible roles for these oscillations are unclear, since if this transformation is linear, the high frequencies in the motoneuron inputs (e.g., 20 Hz from pyramidal tract neurons) would be filtered out by the muscle and have no effect on force control. A biologically inspired mathematical model of the neuromuscular system was used to investigate the impact of high-frequency cortical oscillatory activity on force control. The model simulation results evidenced that a typical motoneuron pool has a nonlinear behavior that enables the decoding of a high-frequency oscillatory input. An input at a single frequency (e.g., beta band) leads to an increase in the steady-state force generated by the muscle. When the input oscillation was amplitude modulated at a given low frequency, the force oscillated at this frequency. In both cases, the mechanism relies on the recruitment and derecruitment of motor units in response to the oscillatory descending drive. Therefore, the results from this study suggest a potential role in force control for cortical oscillations at frequencies at or above the beta band, despite the low-pass behavior of the muscles.

SIGNIFICANCE STATEMENT

The role of cortical oscillations in motor control has been a long-standing question, one view being that they are an epiphenomenon. Fast oscillations are known to reach the spinal cord, and hence they have been thought to affect muscle behavior. However, experimental limitations have hampered further advances to explain how they could influence muscle force. An approach for such a challenge was adopted in the present research: to study the problem through computer simulations of an advanced biologically compatible mathematical model. Using such a model, we found that the well-known mechanism of recruitment and derecruitment of the spinal cord motoneurons can allow the muscle to respond to cortical oscillations, suggesting that these oscillations are not epiphenomena in motor control.

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

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

SIGNIFICANCE STATEMENT

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

High-density surface EMG decomposition allows for recording of motor unit discharge from proximal and distal flexion synergy muscles simultaneously in individuals with stroke

Analysis of motor unit discharge can provide insight into the neural control of movement in healthy and pathological states, but it is typically completed in one muscle at a time. For some research investigations, it would be advantageous to study motor unit discharge from multiple muscles simultaneously. One such example is investigation of the flexion synergy, an abnormal muscle co-activation pattern in post-stroke individuals in which activation of shoulder abductors is involuntarily coupled with that of elbow and finger flexors. However, limitations in available technology have hindered the ability to efficiently extract motor unit discharge from multiple muscles simultaneously. In this study, we propose the use of high-density surface EMG decomposition from proximal and distal flexion synergy muscles (deltoid, biceps, wrist/finger flexors) in combination with an isometric joint torque recording device in individuals with chronic stroke. This innovative approach provides the ability to efficiently analyze both motor units and joint torques that have been simultaneously recorded from the shoulder, elbow, and fingers. In preliminary experiments, 3 stroke and 5 control participants generated shoulder abduction, elbow flexion, and finger flexion torques at 10, 20, 30 and 40% of maximum torque. Motor unit spike trains could be extracted from all muscles at each torque level. Mean motor unit firing rates were significantly lower in the stroke group than in the control group for all three muscles. Within the stroke group, wrist/finger flexor motor units had the lowest coefficient of variation. Additionally, modulation of mean firing rates across torque levels was significantly impaired in all three paretic muscles. The implications of these findings and overall impact of this approach are discussed.

The proportion of common synaptic input to motor neurons increases with an increase in net excitatory input

α-Motor neurons receive synaptic inputs from spinal and supraspinal centers that comprise components either common to the motor neuron pool or independent. The input shared by motor neurons--common input--determines force control. The aim of the study was to investigate the changes in the strength of common synaptic input delivered to motor neurons with changes in force and with fatigue, two conditions that underlie an increase in the net excitatory drive to the motor neurons. High-density surface electromyogram (EMG) signals were recorded from the tibialis anterior muscle during contractions at 20, 50, and 75% of the maximal voluntary contraction force (in 3 sessions separated by at least 2 days), all sustained until task failure. EMG signal decomposition identified the activity of a total of 1,245 motor units. The coherence values between cumulative motor unit spike trains increased with increasing force, especially for low frequencies. This increase in coherence was not observed when comparing two subsets of motor units having different recruitment thresholds, but detected at the same force level. Moreover, the coherence values for frequencies <5 Hz increased at task failure with respect to the beginning of the contractions for all force levels. In conclusion, the results indicated that the relative strength of common synaptic input to motor neurons increases with respect to independent input when the net excitatory drive to motor neurons increases as a consequence of a change in force and fatigue.

The phase difference between neural drives to antagonist muscles in essential tremor is associated with the relative strength of supraspinal and afferent input

The pathophysiology of essential tremor (ET), the most common movement disorder, is not fully understood. We investigated which factors determine the variability in the phase difference between neural drives to antagonist muscles, a long-standing observation yet unexplained. We used a computational model to simulate the effects of different levels of voluntary and tremulous synaptic input to antagonistic motoneuron pools on the tremor. We compared these simulations to data from 11 human ET patients. In both analyses, the neural drive to muscle was represented as the pooled spike trains of several motor units, which provides an accurate representation of the common synaptic input to motoneurons. The simulations showed that, for each voluntary input level, the phase difference between neural drives to antagonist muscles is determined by the relative strength of the supraspinal tremor input to the motoneuron pools. In addition, when the supraspinal tremor input to one muscle was weak or absent, Ia afferents provided significant common tremor input due to passive stretch. The simulations predicted that without a voluntary drive (rest tremor) the neural drives would be more likely in phase, while a concurrent voluntary input (postural tremor) would lead more frequently to an out-of-phase pattern. The experimental results matched these predictions, showing a significant change in phase difference between postural and rest tremor. They also indicated that the common tremor input is always shared by the antagonistic motoneuron pools, in agreement with the simulations. Our results highlight that the interplay between supraspinal input and spinal afferents is relevant for tremor generation.

The optimal neural strategy for a stable motor task requires a compromise between level of muscle cocontraction and synaptic gain of afferent feedback

Increasing joint stiffness by cocontraction of antagonist muscles and compensatory reflexes are neural strategies to minimize the impact of unexpected perturbations on movement. Combining these strategies, however, may compromise steadiness, as elements of the afferent input to motor pools innervating antagonist muscles are inherently negatively correlated. Consequently, a high afferent gain and active contractions of both muscles may imply negatively correlated neural drives to the muscles and thus an unstable limb position. This hypothesis was systematically explored with a novel computational model of the peripheral nervous system and the mechanics of one limb. Two populations of motor neurons received synaptic input from descending drive, spinal interneurons, and afferent feedback. Muscle force, simulated based on motor unit activity, determined limb movement that gave rise to afferent feedback from muscle spindles and Golgi tendon organs. The results indicated that optimal steadiness was achieved with low synaptic gain of the afferent feedback. High afferent gains during cocontraction implied increased levels of common drive in the motor neuron outputs, which were negatively correlated across the two populations, constraining instability of the limb. Increasing the force acting on the joint and the afferent gain both effectively minimized the impact of an external perturbation, and suboptimal adjustment of the afferent gain could be compensated by muscle cocontraction. These observations show that selection of the strategy for a given contraction implies a compromise between steadiness and effectiveness of compensations to perturbations. This indicates that a task-dependent selection of neural strategy for steadiness is necessary when acting in different environments.