Background: Individuals who have experienced spinal cord injury (SCI) may exhibit various muscle-related neurophysiological adaptations, including alterations in motor unit (MU) size and firing behavior. However, due to the technical challenges of in vivo measurement, our understanding of the alterations in the electrophysiological parameters of these MUs remains limited. This study proposed an integrated approach using high-density electromyography (HD-EMG) decomposition and motor neuron (MN) modelling to estimate the intrinsic properties of MUs in vivo and investigated alterations of these properties in persons with SCI.

Methods: HD-EMG signals were recorded during submaximal isometric dorsiflexion and plantar flexion tasks on tibialis anterior (TA), soleus, and gastrocnemius medialis muscles from twenty-six participants with SCI and eighteen non-disabled controls. The HD-EMG signals were subsequently decomposed into MN spike trains and the common synaptic input to the MN pool was estimated. A simplified leaky integrate-and-fire neuron model was then used to simulate MN spiking trains, with soma size and inert period as tunning parameters, which are crucial for MU recruitment and firing patterns, respectively. These parameters were estimated by fitting the instantaneous discharge frequencies of decomposed and simulated spike trains via a genetic algorithm.

Results: The results showed a prolonged inert period in the TA of the persons with SCI. This finding suggested that the MUs in the TA have a slower recovery period before becoming excitable again, which may result in a lower firing rate of MUs in the TA muscle. No significant differences were observed in the soleus and gastrocnemius medialis muscles between the SCI and control groups for either the soma size or inert period parameters.

Conclusions: The simplified leaky integrate-and-fire model exhibited robustness in estimating MN parameters in vivo, offering valuable insights into personalized MU behavior monitoring. To the best knowledge of authors, this is the first study to combine HD-EMG and MU modeling to investigate MU electrophysiological changes in persons with SCI in vivo. This novel approach offers a comprehensive understanding of MU properties adaptations following neurological disorders and informs the development of novel rehabilitation strategies.