Is the 'reverse onion skin' phenomenon more prevalent than we thought during intramuscular myoelectric recordings from low to maximal force outputs?

Characterising sex-related differences in lower- and higher-threshold motor unit behaviour through high-density surface electromyography

Emerging questions in neuromuscular physiology revolve around whether males and females share similar neural control in diverse tasks across a broad range of intensities. In order to explore these features, high-density electromyography was used to record the myoelectrical activity of biceps brachii during trapezoidal isometric contractions at 35% and 70% of maximal voluntary force (MVF) on 11 male and 13 female participants. Identified motor units were then classified as lower-threshold (recruited at ≤30%MVF) and higher-threshold (recruited at >30%MVF). The discharge rate, interspike interval variability, recruitment and derecruitment thresholds, and estimates of neural drive to motor neurons were assessed. Female lower-threshold motor units showed higher neural drive (P < 0.001), accompanied by higher discharge rate at recruitment (P = 0.006), plateau (P = 0.001) and derecruitment (P = 0.001). On the other hand, male higher-threshold motor units showed greater neural drive (P = 0.04), accompanied by higher discharge rate at recruitment (P = 0.005), plateau (P = 0.04) and derecruitment (P = 0.01). Motor unit discharge rate normalised by the recruitment threshold was significantly higher in female lower-threshold motor units (P < 0.001), while no differences were observed in higher-threshold motor units. Recruitment and derecruitment thresholds are higher in males across all intensities (P < 0.01). However, males and females have similar activation and deactivation strategies, as evidenced by similar recruitment-to-derecruitment ratios (P > 0.05). This study encompasses a broad intensity range to analyse motor unit sex-related differences, highlighting higher neural drive and discharge rates in female lower-threshold motor units, elevated recruitment and derecruitment thresholds in males, and convergences in activation and deactivation strategies. HIGHLIGHTS: What is the central question of the study? Do male and female motor units behave similarly in low- and high-intensity contractions? What is the main finding and its importance? Female motor units show higher discharge rates in low-intensity tasks and lower discharge rates in high-intensity tasks, with no differences in recruitment behaviour. A broader inter-spike interval variability was also observed in females. These findings underline that there are sex-specific differences concern the firing strategies based on task intensity.

Sex differences in laterality of motor unit firing behavior of the first dorsal interosseous muscle in strength-matched healthy young males and females

PURPOSE

The purpose of this study was to compare laterality in motor unit firing behavior between females and males.

METHODS

Twenty-seven subjects (14 females) were recruited for this study. The participants performed ramp up and hold isometric index finger abduction at 10, 30, and 60% of their maximum voluntary contraction (MVC). High-density surface electromyography (HD-sEMG) signals were recorded in the first dorsal interosseous (FDI) muscle and decomposed into individual motor unit (MU) firing behavior using a convolution blind source separation method.

RESULTS

In total, 769 MUs were detected (females, n = 318 and males, n = 451). Females had a significantly higher discharge rate than males at each relative torque level (10%: male dominant hand, 13.4 ± 2.7 pps vs. female dominant hand, 16.3 ± 3.4 pps; 30%: male dominant hand, 16.1 ± 3.9 pps vs. female dominant hand, 20.0 ± 5.0 pps; and 60%: male dominant hand, 19.3 ± 3.8 vs. female dominant hand, 25.3 ± 4.8 pps; p < 0.0001). The recruitment threshold was also significantly higher in females than in males at 30 and 60% MVC. Furthermore, males exhibited asymmetrical discharge rates at 30 and 60% MVC and recruitment thresholds at 30 and 60% MVC, whereas no asymmetry was observed in females.

CONCLUSION

In the FDI muscle, compared to males, females exhibited different neuromuscular strategies with higher discharge rates and recruitment thresholds and no asymmetrical MU firing behavior. Notably, the findings that sex differences in neuromuscular activity also occur in healthy individuals provide important information for understanding the pathogenesis of various diseases.

Onion skin is not a universal firing pattern for spinal motoneurons: simulation study

Muscle force is modulated by sequential recruitment and firing rates of motor units (MUs). However, discrepancies exist in the literature regarding the relationship between MU firing rates and their recruitment, presenting two contrasting firing-recruitment schemes. The first firing scheme, known as "onion skin," exhibits low-threshold MUs firing faster than high-threshold MUs, forming separate layers akin to an onion. This contradicts the other firing scheme, known as "reverse onion skin" or "afterhyperpolarization (AHP)," with low-threshold MUs firing slower than high-threshold MUs. To study this apparent dichotomy, we used a high-fidelity computational model that prioritizes physiological fidelity and heterogeneity, allowing versatility in the recruitment of different motoneuron types. Our simulations indicate that these two schemes are not mutually exclusive but rather coexist. The likelihood of observing each scheme depends on factors such as the motoneuron pool activation level, synaptic input activation rates, and MU type. The onion skin scheme does not universally govern the encoding rates of MUs but tends to emerge in unsaturated motoneurons (cells firing < their fusion frequency that generates peak force), whereas the AHP scheme prevails in saturated MUs (cells firing at their fusion frequency), which is highly probable for slow (S)-type MUs. When unsaturated, fast fatigable (FF)-type MUs always show the onion skin scheme, whereas S-type MUs do not show either one. Fast fatigue-resistant (FR)-type MUs are generally similar but show weaker onion skin behaviors than FF-type MUs. Our results offer an explanation for the longstanding dichotomy regarding MU firing patterns, shedding light on the factors influencing the firing-recruitment schemes.NEW & NOTEWORTHY The literature reports two contrasting schemes, namely the onion skin and the afterhyperpolarization (AHP) regarding the relationship between motor units (MUs) firing rates and recruitment order. Previous studies have examined these schemes phenomenologically, imposing one scheme on the firing-recruitment relationship. Here, we used a high-fidelity computational model that prioritizes biological fidelity and heterogeneity to investigate motoneuron firing schemes without bias toward either scheme. Our objective findings offer an explanation for the longstanding dichotomy on MU firing patterns.

Association between physical fitness tests and neuromuscular properties

PURPOSE

While various fitness tests have been developed to assess physical performances, it is unclear how these tests are affected by differences, such as, in morphological and neural factors. This study was aimed to investigate associations between individual differences in physical fitness tests and neuromuscular properties.

METHODS

One hundred and thirty-three young adults participated in various general physical fitness tests and neuromuscular measurements. The appendicular skeletal muscle mass (ASM) was estimated by bioelectrical impedance analysis. Echo intensity (EI) was evaluated from the vastus lateralis. During submaximal knee extension force, high-density surface electromyography of the vastus lateralis was recorded and individual motor unit firings were detected. Y-intercept (i-MU) and slope (s-MU) from the regression line between the recruitment threshold and motor unit firing rate were calculated.

RESULTS

Stepwise multiple regression analyses revealed that knee extension strength could be explained (adjusted R2 = 0.712) by ASM (β = 0.723), i-MU (0.317), EI (- 0.177), and s-MU (0.210). Five-sec stepping could be explained by ASM (adjusted R2 = 0.212). Grip strength, side-stepping, and standing broad jump could be explained by ASM and echo intensity (adjusted R2 = 0.686, 0.354, and 0.627, respectively). Squat jump could be explained by EI (adjusted R2 = 0.640). Counter-movement jump could be explained by EI and s-MU (adjusted R2 = 0.631). On the other hand, i-MU and s-MU could be explained by five-sec stepping and counter-movement jump, respectively, but the coefficients of determination were low (adjusted R2 = 0.100 and 0.045).

CONCLUSION

Generally developed physical fitness tests were mainly explained by morphological factors, but were weakly affected by neural factors involved in performance.

Spatial Dependence of Log-Transformed Electromyography-Force Relation: Model-Based Sensitivity Analysis and Experimental Study of Biceps Brachii

This study investigated electromyography (EMG)-force relations using both simulated and experimental approaches. A motor neuron pool model was first implemented to simulate EMG-force signals, focusing on three different conditions that test the effects of small or large motor units located more or less superficially in the muscle. It was found that the patterns of the EMG-force relations varied significantly across the simulated conditions, quantified by the slope (b) of the log-transformed EMG-force relation. b was significantly higher for large motor units, which were preferentially located superficially rather than for random depth or deep depth conditions (p < 0.001). The log-transformed EMG-force relations in the biceps brachii muscles of nine healthy subjects were examined using a high-density surface EMG. The slope (b) distribution of the relation across the electrode array showed a spatial dependence; b in the proximal region was significantly larger than the distal region, whereas b was not different between the lateral and medial regions. The findings of this study provide evidence that the log-transformed EMG-force relations are sensitive to different motor unit spatial distributions. The slope (b) of this relation may prove to be a useful adjunct measure in the investigation of muscle or motor unit changes associated with disease, injury, or aging.

Neuromuscular recruitment strategies of the vastus lateralis according to sex

AIM

Despite males typically exhibiting greater muscle strength and fatigability than females, it remains unclear if there are sex-based differences in neuromuscular recruitment strategies e.g. recruitment and modulation of motor unit firing rate (MU FR) at normalized forces and during progressive increases in force.

METHODS

The study includes 29 healthy male and 31 healthy female participants (18-35 years). Intramuscular electromyography (iEMG) was used to record individual motor unit potentials (MUPs) and near-fibre MUPs from the vastus lateralis (VL) during 10% and 25% maximum isometric voluntary contractions (MVC), and spike-triggered averaging was used to obtain motor unit number estimates (MUNE) of the VL.

RESULTS

Males exhibited greater muscle strength (P < .001) and size (P < .001) than females, with no difference in force steadiness at 10% or 25% MVC. Females had 8.4% and 6.5% higher FR at 10% and 25% MVC, respectively (both P < .03), while the MUP area was 33% smaller in females at 10% MVC (P < .02) and 26% smaller at 25% MVC (P = .062). However, both sexes showed similar increases in MU size and FR when moving from low- to mid-level contractions. There were no sex differences in any near-fibre MUP parameters or in MUNE.

CONCLUSION

In the vastus lateralis, females produce muscle force via different neuromuscular recruitment strategies to males which is characterized by smaller MUs discharging at higher rates. However, similar strategies are employed to increase force production from low- to mid-level contractions. These findings of similar proportional increases between sexes support the use of mixed sex cohorts in studies of this nature.

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.

Neurophysiological Factors Affecting Muscle Innervation Zone Estimation Using Surface EMG: A Simulation Study

Surface electromyography (EMG) recorded by a linear or 2-dimensional electrode array can be used to estimate the location of muscle innervation zones (IZ). There are various neurophysiological factors that may influence surface EMG and thus potentially compromise muscle IZ estimation. The objective of this study was to evaluate how surface-EMG-based IZ estimation might be affected by different factors, including varying degrees of motor unit (MU) synchronization in the case of single or double IZs. The study was performed by implementing a model simulating surface EMG activity. Three different MU synchronization conditions were simulated, namely no synchronization, medium level synchronization, and complete synchronization analog to M wave. Surface EMG signals recorded by a 2-dimensional electrode array were simulated from a muscle with single and double IZs, respectively. For each situation, the IZ was estimated from surface EMG and compared with the one used in the model for performance evaluation. For the muscle with only one IZ, the estimated IZ location from surface EMG was consistent with the one used in the model for all the three MU synchronization conditions. For the muscle with double IZs, at least one IZ was appropriately estimated from interference surface EMG when there was no MU synchronization. However, the estimated IZ was different from either of the two IZ locations used in the model for the other two MU synchronization conditions. For muscles with a single IZ, MU synchronization has little effect on IZ estimation from electrode array surface EMG. However, caution is required for multiple IZ muscles since MU synchronization might lead to false IZ estimation.