Identification of motor unit discharges from ultrasound images: Analysis of in silico and in vivo experiments

Abstract

ObjectiveUltrasound (US) images during a muscle contraction can be decoded into individual motor unit (MU) activity, i.e., trains of neural discharges from the spinal cord. However, current decoding algorithms assume a stationary mixing matrix, i.e. equal mechanical twitches at each discharge. This study aimed to investigate the accuracy of these approaches in non-ideal conditions when the mechanical twitches in response to neural discharges vary over time and are partially fused in tetanic contractions. MethodsWe performed an in silico experiment to study the decomposition accuracy for changes in simulation parameters, including the twitch waveforms, spatial territories, and motoneuron-driven activity. Then, we explored the consistency of the in silico findings with an in vivo experiment on the tibialis anterior muscle at varying contraction forces. ResultsA large population of MU spike trains across different excitatory drives, and noise levels could be identified. The identified MUs with varying twitch waveforms resulted in varying amplitudes of the estimated sources correlated with the ground truth twitch amplitudes. The identified spike trains had a wide range of firing rates, and the later recruited MUs with larger twitch amplitudes were easier to identify than those with small amplitudes. Finally, the in silico and in vivo results were consistent, and the method could identify MU spike trains in US images at least up to 40% of the maximal voluntary contraction force. ConclusionThe decoding method was accurate irrespective of the varying twitch-like shapes or the degree of twitch fusion, indicating robustness, important for neural interfacing applications.