Factors affecting the stability of the spike-triggered averaged force in the human first dorsal interosseus muscle

Estimation of contractile parameters of successive twitches in unfused tetanic contractions of single motor units - A proof-of-concept study using ultrafast ultrasound imaging in vivo

During a voluntary contraction, motor units (MUs) fire a train of action potentials, causing summation of the twitch forces, resulting in fused or unfused tetanus. Twitches have been important in studying whole-muscle contractile properties and differentiation between MU types. However, there are still knowledge gaps concerning the voluntary force generation mechanisms. Current methods rely on the spike-triggered averaging technique, which cannot track changes in successive twitches' properties in response to individual neural firings. This study proposes a method that estimates successive twitches contractile parameters of single MUs during low force voluntary isometric contractions in human biceps brachii. We used a previously developed ultrafast ultrasound imaging method to estimate unfused tetanic activity signals of single MUs. A twitch decomposition model was used to decompose unfused tetanic activity signals into individual twitches. This study found that the contractile parameters varied within and across MUs. There was an association between the inter-spike interval and the contraction time (r = 0.49,p < 0.001) and the half-relaxation time (r = 0.58,p < 0.001), respectively. The method shows the proof-of-concept to study MU contractile properties of individual twitches in vivo, which can provide further insights into the force generation mechanisms of voluntary contractions and response to individual neural discharges.

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.

Limitations of the spike-triggered averaging for estimating motor unit twitch force: a theoretical analysis

Contractile properties of human motor units provide information on the force capacity and fatigability of muscles. The spike-triggered averaging technique (STA) is a conventional method used to estimate the twitch waveform of single motor units in vivo by averaging the joint force signal. Several limitations of this technique have been previously discussed in an empirical way, using simulated and experimental data. In this study, we provide a theoretical analysis of this technique in the frequency domain and describe its intrinsic limitations. By analyzing the analytical expression of STA, first we show that a certain degree of correlation between the motor unit activities prevents an accurate estimation of the twitch force, even from relatively long recordings. Second, we show that the quality of the twitch estimates by STA is highly related to the relative variability of the inter-spike intervals of motor unit action potentials. Interestingly, if this variability is extremely high, correct estimates could be obtained even for high discharge rates. However, for physiological inter-spike interval variability and discharge rate, the technique performs with relatively low estimation accuracy and high estimation variance. Finally, we show that the selection of the triggers that are most distant from the previous and next, which is often suggested, is not an effective way for improving STA estimates and in some cases can even be detrimental. These results show the intrinsic limitations of the STA technique and provide a theoretical framework for the design of new methods for the measurement of motor unit force twitch.

Extraction of individual muscle mechanical action from endpoint force

Most motor tasks require the simultaneous coordination of multiple muscles. That coordination is poorly understood in part because there is no noninvasive means of isolating a single muscle's contribution to the resultant endpoint force. The contribution of a single motor unit to isometric tasks can, however, be characterized using the spike-triggered averaging (STA) technique, applied to a single motor unit's spike train. We propose that a technique analogous to STA, which we call electromyogram (EMG)-weighted averaging (EWA), can be applied to surface EMGs to extract muscle mechanical action from the natural endpoint force fluctuations generated during steady isometric contraction. We demonstrate this technique on simultaneous recordings of fingertip force and surface EMG from the first dorsal interosseous (FDI) and extensor indicis (EI) of humans. The EWA direction was approximately the same across a wide range of fingertip force directions, and the average EWA direction was consistent with mechanical action direction of these muscles estimated from cadaveric and imaging data: the EWA directions were 193 +/- 2 degrees for the FDI and 71 +/- 5 degrees for the EI (95% confidence). EWA transient behavior also appears to capture temporal characteristics of muscle force fluctuations with peak force time and general waveform shape similar to that of the associated spike-triggered averages from single motor units. The EWA may provide a means of empirically characterizing the complex transformation between muscle force and endpoint force without the need for invasive electrode recordings or complex anatomical measurements of musculoskeletal geometry.

Afterhyperpolarization time-course and minimal discharge rate in low threshold motor units in humans

The alpha-motoneurone afterhyperpolarization (AHP) duration correlates with a number of its muscle unit properties in animal preparations. In humans, the interval death rate (IDR) analysis has been used to estimate the time course of human motoneurone AHP based on the pattern of motor unit firing. The purpose of this experiment was first, to examine the relationship between estimated AHP time course and the minimal firing rate of the motor unit and second, to examine the relationship between the AHP and motor unit contractile properties in the tibialis anterior (TA) muscle. Motor unit data were obtained from the TA muscle during low force isometric contractions lasting 600 s. Muscle unit twitch characteristics were determined using spike-triggered averaging (STA) and the motoneurone AHP time course was estimated using the IDR analysis. Minimal discharge rate and derecruitment threshold torque were determined for 2 s preceding motor unit derecruitment. The AHP time constant and minimal discharge rate were negatively correlated, whereas the derecruitment threshold torque was not associated with the AHP time constant. The estimated AHP duration, however, is considerably shorter than the mean ISI of the minimal discharge rate suggesting that synaptic noise and AHP duration are important factors in dictating the minimal discharge rate in low force voluntary contractions in humans. The AHP time constant did not vary significantly with motor unit twitch amplitude; however, significant positive relationships were found between the AHP time constant and the temporal properties of the motor unit twitch. The calculated AHP time course using the IDR analysis, therefore, is a reasonable estimate and coupled with motor unit properties attained with STA, it provides a powerful method to describe low-threshold motor units.

Variability of quadriceps femoris motor neuron discharge and muscle force in human aging

The purpose was to determine the contribution of visual feedback and the effect of aging on the variability of knee extensor (KE) muscle force and motor unit (MU) discharge. Single MUs were recorded during two types of isometric trials, (1) visual feedback provided (VIS) and then removed (NOVIS) during the trial (34 MUs from young, 32 from elderly), and (2) only NOVIS (66 MUs from young, 77 from elderly) during the trial. Recruitment threshold (RT) ranged from 0-37% MVC. Standard deviation (SD) and coefficient of variation (CV) of muscle force and MU interspike interval (ISI) was measured during steady contractions at target forces ranging from 0.3 to 54% MVC. Force drift (<0.5 Hz) was removed before analysis. VIS/NOVIS trials: the decrease in the CV of ISI from VIS to NOVIS was greater for MUs from elderly (12.5 +/- 4.1 to 9.94 +/- 2.6%) than young (10.6 +/- 3.3 to 10.3 +/- 2.8%, age group x vision interaction, P = 0.006). The change in CV of force from VIS to NOVIS was significantly greater for elderly (1.45 to 1.05%) than young (1.42 to 1.41%). NOVIS only trials: for all MUs, the average RT (6.6 +/- 7.7 % MVC), target force above RT (1.20 +/- 2.7% MVC), SD of ISI (0.012 +/- 0.005 s), and CV of ISI (11.1 +/- 3.3%) were similar for young and elderly MUs. The CV of force was similar between age groups for trials between 0 and 3% MVC (1.74 +/- 0.74%) and was greater for young subjects from 3 to 10% MVC (1.47 +/- 0.5 vs. 1.21 +/- 0.4%) and >10% MVC (1.44 +/- 0.6 vs. 1.01 +/- 0.3%). The CV of ISI was similar between age groups for MUs in 0-3, 3-10, and >10% bins of RT. Thus, the contribution of visuomotor correction to the variability of motor unit discharge and force is greater for elderly adults. The presence of visual feedback appears to be necessary to find greater discharge variability in motor units from the knee extensors of elderly adults.

Simulations of motor unit number estimation techniques

Motor unit number estimation (MUNE) is an electrodiagnostic procedure used to evaluate the number of motor axons connected to a muscle. All MUNE techniques rely on assumptions that must be fulfilled to produce a valid estimate. As there is no gold standard to compare the MUNE techniques against, we have developed a model of the relevant neuromuscular physiology and have used this model to simulate various MUNE techniques. The model allows for a quantitative analysis of candidate MUNE techniques that will hopefully contribute to consensus regarding a standard procedure for performing MUNE.

Intraneural microstimulation of motor axons in the study of human single motor units

Single motor unit activity has been studied in depth since the first intramuscular electrodes were developed more than 70 years ago. Many techniques have been combined or used in isolation since then. Intraneural motor axon microstimulation allows the detailed study of single motor units in awake human subjects in a manner most analogous to that used in reduced animal preparations. A microelectrode, inserted percutaneously into a peripheral nerve, stimulates the axon of a single alpha-motoneuron at a site remote from the contracting muscle, allowing detailed analyses of the contractile properties of a single motor unit in an otherwise quiescent muscle, that is, without interference of simultaneously active motor units or the presence of an electrode within the muscle. The methods and results obtained using this technique are described and compared to those of other studies of single motor units in human subjects. Differences have been found between human and animal motor units and between motor units of various muscles. Studying human and animal motor units using an analogous technique provides insight into the interpretation of human data when results differ from animal data, and when human motor units cannot be examined in the same way, or at a similar level of detail, as animal motor units.

The time course of the motoneurone afterhyperpolarization is related to motor unit twitch speed in human skeletal muscle

The relationship between the electrophysiological properties of motoneurones and their muscle units has been established in animal models. A functionally significant relationship exists whereby motoneurones with long post-spike afterhyperpolarizations (AHPs) innervate slow contracting muscle units. The purpose of this study was to determine whether the time course of the AHP as measured by its time constant is associated with the contractile properties of its muscle unit in humans. Using an intramuscular fine wire electrode, 46 motor units were recorded in eight subjects as they held a low force contraction of the first dorsal interosseus muscle for approximately 10 min. By applying a recently validated transform to the interspike interval histogram, the mean voltage versus time trajectory of the motoneurone AHP was determined. Spike-triggered averaging was used to extract the muscle unit twitch from the whole muscle force with strict control over force variability and motor unit discharge rate (interspike intervals between 120 and 200 ms). The AHP time constant was positively correlated to the time to half-force decay (rho = 0.36, P < 0.05) and twitch duration (rho = 0.57, P < 0.001); however, time to peak force failed to reach significance (rho = 0.27, P < 0.07). These results suggest that a similar functional relationship exists in humans between the motoneurone AHP and the muscle unit contractile properties.

Ischemia sensitivity and motoneuron afterhyperpolarization in human motor units

The purpose of this study was to examine whether motor units, grouped by speed of contraction in the human first dorsal interosseous muscle, differed in their sensitivity to ischemia and motoneuron afterhyperpolarization (AHP) time-course. Motor units were recorded while subjects held an abduction force for approximately 10 min. Subsequently, subjects abducted for 4-5 min under ischemic conditions. Motor unit twitches derived using spike-triggered averaging were allocated into "fast" or "slow" contracting groups based on twitch time to peak (TTP) force. Motor units in the "slow" group had a greater sensitivity to ischemia than the "fast" group. When upper and lower quartiles of TTP were compared, motor units with slow TTP had long AHP time-constants (as estimated by an interspike interval histogram transform). Thus, motor units grouped by speed of contraction differed in both their sensitivity to ischemia and motoneuron AHP time-course. This provides preliminary evidence that the estimated AHP time-constant may be used to deduce motor-unit type in humans.