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

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

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

Spatially repeatable components from ultrafast ultrasound are associated with motor unit activity in human isometric contractions. J Neural Eng 2023;20:046016. [PMID: 37437598 DOI: 10.1088/1741-2552/ace6fc]

Objective.Ultrafast ultrasound (UUS) imaging has been used to detect intramuscular mechanical dynamics associated with single motor units (MUs). Detecting MUs from ultrasound sequences requires decomposing a velocity field into components, each consisting of an image and a signal. These components can be associated with putative MU activity or spurious movements (noise). The differentiation between putative MUs and noise has been accomplished by comparing the signals with MU firings obtained from needle electromyography (EMG). Here, we examined whether the repeatability of the images over brief time intervals can serve as a criterion for distinguishing putative MUs from noise in low-force isometric contractions.Approach.UUS images and high-density surface EMG (HDsEMG) were recorded simultaneously from 99 MUs in the biceps brachii of five healthy subjects. The MUs identified through HDsEMG decomposition were used as a reference to assess the outcomes of the ultrasound-based components. For each contraction, velocity sequences from the same eight-second ultrasound recording were separated into consecutive two-second epochs and decomposed. To evaluate the repeatability of components' images across epochs, we calculated the Jaccard similarity coefficient (JSC). JSC compares the similarity between two images providing values between 0 and 1. Finally, the association between the components and the MUs from HDsEMG was assessed.Main results.All the MU-matched components had JSC > 0.38, indicating they were repeatable and accounted for about one-third of the HDsEMG-detected MUs (1.8 ± 1.6 matches over 4.9 ± 1.8 MUs). The repeatable components (JSC > 0.38) represented 14% of the total components (6.5 ± 3.3 components). These findings align with our hypothesis that intra-sequence repeatability can differentiate putative MUs from noise and can be used for data reduction.Significance.This study provides the foundation for developing stand-alone methods to identify MU in UUS sequences and towards real-time imaging of MUs. These methods are relevant for studying muscle neuromechanics and designing novel neural interfaces.

Skeletal muscle models composed of motor units: A review

The mathematical muscle models should include several aspects of muscle structure and physiology. First, muscle force is the sum of forces of multiple motor units (MUs), which have different contractile properties and play different roles in generating muscle force. Second, whole muscle activity is an effect of net excitatory inputs to a pool of motoneurons innervating the muscle, which have different excitability, influencing MU recruitment. In this review, we compare various methods for modeling MU twitch and tetanic forces and then discuss muscle models composed of different MU types and number. We first present four different analytical functions used for twitch modeling and show limitations related to the number of twitch describing parameters. We also show that a nonlinear summation of twitches should be considered in modeling tetanic contractions. We then compare different muscle models, most of which are variations of Fuglevand's model, adopting a common drive hypothesis and the size principle. We pay attention to integrating previously developed models into a consensus model based on physiological data from in vivo experiments on the rat medial gastrocnemius muscle and its respective motoneurons. Finally, we discuss the shortcomings of existing models and potential applications for studying MU synchronization, potentiation, and fatigue.

Estimating the neural spike train from an unfused tetanic signal of low-threshold motor units using convolutive blind source separation

BACKGROUND

Individual motor units have been imaged using ultrafast ultrasound based on separating ultrasound images into motor unit twitches (unfused tetanus) evoked by the motoneuronal spike train. Currently, the spike train is estimated from the unfused tetanic signal using a Haar wavelet method (HWM). Although this ultrasound technique has great potential to provide comprehensive access to the neural drive to muscles for a large population of motor units simultaneously, the method has a limited identification rate of the active motor units. The estimation of spikes partly explains the limitation. Since the HWM may be sensitive to noise and unfused tetanic signals often are noisy, we must consider alternative methods with at least similar performance and robust against noise, among other factors.

RESULTS

This study aimed to estimate spike trains from simulated and experimental unfused tetani using a convolutive blind source separation (CBSS) algorithm and compare it against HWM. We evaluated the parameters of CBSS using simulations and compared the performance of CBSS against the HWM using simulated and experimental unfused tetanic signals from voluntary contractions of humans and evoked contraction of rats. We found that CBSS had a higher performance than HWM with respect to the simulated firings than HWM (97.5 ± 2.7 vs 96.9 ± 3.3, p < 0.001). In addition, we found that the estimated spike trains from CBSS and HWM highly agreed with the experimental spike trains (98.0% and 96.4%).

CONCLUSIONS

This result implies that CBSS can be used to estimate the spike train of an unfused tetanic signal and can be used directly within the current ultrasound-based motor unit identification pipeline. Extending this approach to decomposing ultrasound images into spike trains directly is promising. However, it remains to be investigated in future studies where spatial information is inevitable as a discriminating factor.