Abstract
Interference surface electromyogram (EMG) reflects many bioelectric properties of active motor units (MU), which are however difficult to estimate due to the asynchronous summation of their discharges. This paper introduces a deconvolution technique to estimate the cumulative firings of MUs. Tests in simulations show that the power spectral density of the estimated MU firings has a low-frequency peak corresponding to the mean firing rate of MUs in the detection volume of the recording system, weighted by the amplitudes of MU action potentials. The peak increases in amplitude and its centroid shifts to a higher frequency when MU synchronization is simulated (mainly due to the shift of discharges of large MUs). The peak is found even at high force levels, when such a contribution does not emerge from the EMG. This result is also confirmed in preliminary applications to experimental data. Moreover, the simulated cumulative firings of MUs are estimated with a correlation above 90% (considering frequency contributions up to 150 Hz), for all force levels. The method requires a single EMG channel, thus being feasible even in applied studies using simple recording systems. It may open many potential applications, e.g., in the study of the modulation of MU firing rate induced by either fatigue or pathology and in coherency analysis.

Examples of application of the deconvolution (Deconv) algorithm and comparison with the cumulative firings and the cumulated weighted firings (CWF, i.e., each firing pattern is weighted by the root mean squared amplitude of the corresponding MU action potential). Portions of data are shown on the left, the power spectral densities (PSD) on the right (Welch method applied to 3 s of data, sub-epochs of 0.5 s, mean value removed from each of them, 50% of overlap). A) Simulated signal (50% of maximal voluntary contraction, MVC) with random MU firings. B) Simulated signal (50% MVC) with a level of synchronization equal to 10%. C) Experimental data from vastus medialis at 40% MVC (data decomposed by the algorithm of Holobar and Zazula, IEEE Trans. Sig. Proc. 2007; PSD of the cumulated firings almost identical to that of CWF, as few MUs were identified).







Similar content being viewed by others
References
Baker SN, Kilner JM, Pinches EM, Lemon RN (1999) The role of synchrony and oscillations in the motor output. Exp Brain Res 128(1-2):109–17
Basmajian J, De Luca CJ (1985) Muscles alive: Their function revealed by electromyography, 5th edn. Williams and Wilkins, Baltimore
Becker S, von Werder SCFA, Lassek AK, Disselhorst-Klug C (2019) Time-frequency coherence of categorized sEMG data during dynamic contractions of biceps, triceps, and brachioradialis as an approach for spasticity detection. Med Biol Eng Comput. in press
Boyd S, Vandenberghe L (2004) Convex optimization. Cambridge University Press, Cambridge
Burrus CS Iterative reweighted least squares. OpenStax CNX. Available online: http://cnx.org/contents/92b90377-2b34-49e4-b26f-7fe572db78a1@12
Chen C, Chai G, Guo W, Sheng X, Farina D, Zhu X (2019) Prediction of finger kinematics from discharge timings of motor units: implications for intuitive control of myoelectric prostheses. J Neural Eng 16 (2):026005
Contessa P, Adam A, De Luca CJ (1985) Motor unit control and force fluctuation during fatigue. J Appl Physiol 107(1):235–43
Contessa P, De Luca CJ (2013) Neural control of muscle force: indications from a simulation model. J Neurophysiol 109:1548–70
Conway BA, Halliday DM, Farmer SF, Shahani U, Maas P, Weir AI, Rosenberg JR (1995) Synchronization between motor cortex and spinal motoneuronal pool during the performance of a maintained motor task in man. J Physiol 489(Pt 3):917–24
Craven D, McGinley B, Kilmartin L, Glavin M, Jones E (2015) Compressed sensing for bioelectric signals: a review. IEEE J Biomed Health Inform 19:529–40
De Luca CJ, Erim Z (1994) Common drive of motor units in regulation of muscle force. Trends Neurosci 17:299–305
De Luca CJ, Erim Z (2002) Common drive in motor units of a synergistic muscle pair. J Neurophysiol 87 (4):2200–4
De Luca CJ, Adam A, Wotiz R, Gilmore LD, Nawab SH (2006) Decomposition of surface EMG signals. J Neurophysiol 96:1646–57
de Souza LML, Cabral HV, de Oliveira LF, Vieira TM (2018) Motor units in vastus lateralis and in different vastus medialis regions show different firing properties during low-level, isometric knee extension contraction. Hum Mov Sci 58:307–14
Enck P, Franz H, Davico E, Mastrangelo F, Mesin L, Merletti R (2010) Repeatability of innervation zone identification in the external anal sphincter muscle. Neurourol Urodyn 29(3):449–57
Farina D, Merletti R, Enoka RM (2004) The extraction of neural strategies from the surface EMG. J Appl Physiol (1985) 96(4):1486–95. Review
Farina D, Merletti R, Enoka RM (2014) The extraction of neural strategies from the surface EMG: an update. J Appl Physiol (1985) 117(11):1215–30. Review
Herda TJ, Siedlik JA, Trevino MA, Cooper MA, Weir JP (2015) Motor unit control strategies of endurance- versus resistance-trained individuals. Muscle Nerve 52:832–843
Holobar A, Zazula D (2007) Multichannel blind source separation using convolution kernel compensation. IEEE Trans Sig Proc 55:4487–96
Hu X, Suresh AK, Rymer WZ, Suresh NL (2016) Altered motor unit discharge patterns in paretic muscles of stroke survivors assessed using surface electromyography. J Neural Eng 13(4):046025
Lago PJ, Jones NB (1981) Low-frequency spectral analysis of the e.m.g. Med Biol Eng Comput 19(6):779–82
Mesin L, Cocito D (2007) A new method for the estimation of motor nerve conduction block. Clin Neurophysiol 118(4):730–40
Mesin L, Damiano L, Farina D (2007) Estimation of average muscle fiber conduction velocity from simulated surface EMG in pinnate muscles. J Neurosci Methods 160:327–34
Mesin L, Merletti R, Vieira TM (2011) Insights gained into the interpretation of surface electromyograms from the gastrocnemius muscles: A simulation study. J Biomech 44(6):1096–103
Mesin L, Dardanello D, Rainoldi A, Boccia G (2016) Motor unit firing rates and synchronisation affect the fractal dimension of simulated surface electromyogram during isometric/isotonic contraction of vastus lateralis muscle. Med Eng Phys 38(12):1530–33
Mesin L (2019) Separation of interference surface electromyogram into propagating and non-propagating components. Biomed Signal Process Control 52:238–47
Mesin L (2019) Non-propagating components of surface electromyogram reflect motor unit firing rates, submitted to IEEE Access
Myers LJ, Lowery M, O’Malley M, Vaughan CL, Heneghan C, St Clair Gibson A, Harley YX, Sreenivasan R (2003) Rectification and non-linear pre-processing of EMG signals for cortico-muscular analysis. J Neurosci Methods 124(2):157–65
Nawab SH, Chang SS, De Luca CJ (2010) High-yield decomposition of surface EMG signals. Clin Neurophysiol 121(10):1602–15
Neto OP, Christou EA (2010) Rectification of the EMG signal impairs the identification of oscillatory input to the muscle. J Neurophysiol 103(2):1093–103
Piitulainen H, Botter A, Bourguignon M, Jousmäki V, Hari R (2015) Spatial variability in cortex-muscle coherence investigated with magnetoencephalography and high-density surface electromyography. J Neurophysiol 114(5):2843–53
Van Boxtel A, Schomaker LR (1983) Motor unit firing rate during static contraction indicated by the surface EMG power spectrum. IEEE Trans Biomed Eng 30(9):601–9
Weytjens JL, van Steenberghe D (1984) The effects of motor unit synchronization on the power spectrum of the electromyogram. Biol Cybern 51(2):71–7
Acknowledgments
Experimental data were provided by LISiN (Laboratorio di Ingegneria del Sistema Neuromuscolare e della riabilitazione motoria, Turin, Italy). Decomposition of experimental EMGs was performed by A. Botter and T.M. Vieira.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Mesin, L. Single channel surface electromyogram deconvolution to explore motor unit discharges. Med Biol Eng Comput 57, 2045–2054 (2019). https://doi.org/10.1007/s11517-019-02010-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11517-019-02010-0