Key points: Previous studies have indicated that several weeks of strength training is sufficient to elicit significant adaptations in the neural drive sent to the muscles. There are few data, however, on the changes elicited by strength training in the recruitment and rate coding of motor units during voluntary contractions. We show for the first time that the discharge characteristics of motor units in the tibialis anterior muscle tracked across the intervention are changed by 4 weeks of strength training with isometric voluntary contractions. The specific adaptations included significant increases in motor unit discharge rate, decreases in the recruitment‐threshold force of motor units and a similar input–output gain of the motor neurons. The findings suggest that the adaptations in motor unit function may be attributable to changes in synaptic input to the motor neuron pool or to adaptations in intrinsic motor neuron properties. The strength of a muscle typically begins to increase after only a few sessions of strength training. This increase is usually attributed to changes in the neural drive to muscle as a result of adaptations at the cortical or spinal level. We investigated the change in the discharge characteristics of large populations of longitudinally tracked motor units in tibialis anterior before and after 4 weeks of strength training the ankle‐dorsiflexor muscles with isometric contractions. The adaptations exhibited by 14 individuals were compared with 14 control subjects. High‐density electromyogram grids with 128 electrodes recorded the myoelectric activity during isometric ramp contractions to the target forces of 35%, 50% and 70% of maximal voluntary force. The motor unit recruitment and derecruitment thresholds, discharge rate, interspike intervals and estimates of synaptic inputs to motor neurons were assessed. The normalized recruitment‐threshold forces of the motor units were decreased after strength training (P < 0.05). Moreover, discharge rate increased by 3.3 ± 2.5 pps (average across subjects and motor units) during the plateau phase of the submaximal isometric contractions (P < 0.001). Discharge rates at recruitment and derecruitment were not modified by training (P < 0.05). The association between force and motor unit discharge rate during the ramp‐phase of the contractions was also not altered by training (P < 0.05). These results demonstrate for the first time that the increase in muscle force after 4 weeks of strength training is the result of an increase in motor neuron output from the spinal cord to the muscle.

Key points: We propose and validate a method for accurately identifying the activity of populations of motor neurons during contractions at maximal rate of force development in humans. The behaviour of the motor neuron pool during rapid voluntary contractions in humans is presented. We show with this approach that the motor neuron recruitment speed and maximal motor unit discharge rate largely explains the individual ability in generating rapid force contractions. The results also indicate that the synaptic inputs received by the motor neurons before force is generated dictate human potential to generate force rapidly. This is the first characterization of the discharge behaviour of a representative sample of human motor neurons during rapid contractions. During rapid contractions, motor neurons are recruited in a short burst and begin to discharge at high frequencies (up to >200 Hz). In the present study, we investigated the behaviour of relatively large populations of motor neurons during rapid (explosive) contractions in humans, applying a new approach to accurately identify motor neuron activity simultaneous to measuring the rate of force development. The activity of spinal motor neurons was assessed by high‐density electromyographic decomposition from the tibialis anterior muscle of 20 men during isometric explosive contractions. The speed of motor neuron recruitment and the instantaneous motor unit discharge rate were analysed as a function of the impulse (the time–force integral) and the maximal rate of force development. The peak of motor unit discharge rate occurred before force generation and discharge rates decreased thereafter. The maximal motor unit discharge rate was associated with the explosive force variables, at the whole population level (r² = 0.71 ± 0.12; P < 0.001). Moreover, the peak motor unit discharge and maximal rate of force variables were correlated with an estimate of the supraspinal drive, which was measured as the speed of motor unit recruitment before the generation of afferent feedback (P < 0.05). We show for the first time the full association between the effective neural drive to the muscle and human maximal rate of force development. The results obtained in the present study indicate that the variability in the maximal contractile explosive force of the human tibialis anterior muscle is determined by the neural activation preceding force generation.

It has been hypothesized that the coordinated activation of muscles is controlled by the central nervous system by means of a small alphabet of control signals (also referred to as activation signals) and motor modules (synergies). We analyzed the locomotion of 10 patients recently affected by stroke (maximum of 20 wk) and compared it with that of healthy controls. The aim was to assess whether the walking of subacute stroke patients is based on the same motor modules and/or activation signals as healthy subjects. The activity of muscles of the lower and upper limb and the trunk was measured and used for extracting motor modules. Four modules were sufficient to explain the majority of variance in muscle activation in both controls and patients. Modules from the affected side of stroke patients were different from those of healthy controls and from the unaffected side of stroke patients. However, the activation signals were similar between groups and between the affected and unaffected side of stroke patients, and were characterized by impulses at specific time instants within the gait cycle, underlying an impulsive controller of gait. In conclusion, motor modules observed in healthy subjects during locomotion are different from those used by subacute stroke patients, despite similar impulsive activation signals. We suggest that this pattern is consistent with a neuronal network in which the timing of activity generated by central pattern generators is directed to the motoneurons via a premotor network that distributes the activity in a task-dependent manner determined by sensory and descending control information.

The surface interference EMG signal provides some information on the neural drive to muscles. However, the association between neural drive to muscle and muscle activation has long been debated with controversial indications due to the unavailability of motor unit population data. In this study, we clarify the potential and limitations of interference EMG analysis to infer motor unit recruitment strategies with an experimental investigation of several concurrently active motor units and of the associated features of the surface EMG. For this purpose, we recorded high-density surface EMG signals during linearly increasing force contractions of the tibialis anterior muscle, up to 70% of maximal force. The recruitment threshold (RT), conduction velocity (MUCV), median frequency (MDF_MU), and amplitude (RMS_MU) of action potentials of 587 motor units from 13 individuals were assessed and associated with features of the interference EMG. MUCV was positively associated with RT ( R² = 0.64 ± 0.14), whereas MDF_MU and RMS_MU showed a weaker relation with RT ( R² = 0.11 ± 0.11 and 0.39 ± 0.24, respectively). Moreover, the changes in average conduction velocity estimated from the interference EMG predicted well the changes in MUCV ( R² = 0.71), with a strong association to ankle dorsiflexion force ( R² = 0.81 ± 0.12). Conversely, both the average EMG MDF and RMS were poorly associated with motor unit recruitment. These results clarify the limitations of EMG spectral and amplitude analysis in inferring the neural strategies of muscle control and indicate that, conversely, the average conduction velocity could provide relevant information on these strategies. NEW & NOTEWORTHY The surface EMG provides information on the neural drive to muscles. However, the associations between EMG features and neural drive have been long debated due to unavailability of motor unit population data. Here, by using novel highly accurate decomposition of the EMG, we related motor unit population behavior to a wide range of voluntary forces. The results fully clarify the potential and limitation of the surface EMG to provide estimates of the neural drive to muscles.

Amplitude and frequency content of the surface electromyographic (EMG) signal reflect central and peripheral modifications of the neuromuscular system. Classic surface EMG spectral variables applied to assess muscle functions are the centroid and median power spectral frequencies. More recently, nonlinear tools have been introduced to analyze the surface EMG; among them, the recurrence quantification analysis (RQA) was shown to be particularly promising for the detection of muscle status changes. The purpose of this work was to analyze the effect of motor unit short-term synchronization and conduction velocity (CV) on EMG spectral variables and two variables extracted by RQA, the percentage of recurrence (%Rec) and determinism (%Det). The study was performed on the basis of a simulation model, which allowed changing the degree of synchronization and mean CV of a number of motor units, and of an experimental investigation of the surface EMG signal properties detected during high-force-level isometric fatiguing contractions of the biceps brachii muscle. Simulations and experimental results were largely in agreement and show that 1) spectral variables, %Rec, and %Det are influenced by CV and degree of synchronization; 2) spectral variables are highly correlated with %Det ( R = −0.95 in the simulations and −0.78 and −0.75 for the initial values and normalized slopes, respectively, in the experimental signals), and thus the information they provide on muscle properties is basically the same; and 3) variations of %Det and %Rec in response to changes in muscle properties are significantly larger than the variations of spectral variables. This study validates RQA as a means for fatigue assessment with potential advantages (such as the higher sensitivity to changes of muscle status) with respect to the classic spectral analysis.