Marginal reversible jump Markov chain Monte Carlo with application to motor unit number estimation

Motor unit number estimation (MUNE) is a method which aims to provide a quantitative indicator of progression of diseases that lead to loss of motor units, such as motor neurone disease. However the development of a reliable, repeatable and fast real-time MUNE method has proved elusive hitherto. Ridall et al. (2007) implement a reversible jump Markov chain Monte Carlo (RJMCMC) algorithm to produce a posterior distribution for the number of motor units using a Bayesian hierarchical model that takes into account biological information about motor unit activation. However we find that the approach can be unreliable for some datasets since it can suffer from poor cross-dimensional mixing. Here we focus on improved inference by marginalising over latent variables to create the likelihood. In particular we explore how this can improve the RJMCMC mixing and investigate alternative approaches that utilise the likelihood (e.g. DIC (Spiegelhalter et al., 2002)). For this model the marginalisation is over latent variables which, for a larger number of motor units, is an intractable summation over all combinations of a set of latent binary variables whose joint sample space increases exponentially with the number of motor units. We provide a tractable and accurate approximation for this quantity and also investigate simulation approaches incorporated into RJMCMC using results of Andrieu and Roberts (2009).

Decreased motor unit firing rate in the potentiated tibialis anterior in humans

With repeated activity, force production, rate of force production, and relaxation time are impaired. These are characteristics ofa fatigued muscle (Vandenboom, 2004). However, brief bouts of near maximal to maximal activity results in the increased ability of the muscle to generate force, termed post activation potentiation (P AP)(V andervoort et aI., 1983). The purpose of the present study was to characterize motor unit firing rate (MUFR) in the unfatigued, potentiated tibialis anterior (TA). Using a quadrifilar needle electrode, MUFR was measured during a 5s 50% MVC in which the TA was either potentiated or unpotentiated; monopolar electrodes measured surface parameters. A lOs MVC was used to potentiate the muscle. Firing rate decreased significantly from 20.15±2.9Opps to 18.27±2.99pps, while mean power frequency decreased significantly from 60. 13±7.75 Hz to 53.62±8.56 Hz. No change in root mean square (RMS) was observed. Therefore, in the present study, MUFR decreases in response to a potentiated TA.

Motor unit synchronization of the vasti muscles in closed and open chain tasks

Objectives: To investigate motor unit synchronization between medial and lateral vasti and whether such synchronization differs in closed and open chain tasks. Design: Electromyographic recordings of single motor unit action potentials were made from the vastus medialis obliquus (VMO) and multiunit recordings from vastus lateralis during isometric contractions at 30 degrees of knee flexion in closed and open chain conditions. Setting: Laboratory. Participants: Five volunteers with no history of knee pain (age, 30 +/- 3.32y). Interventions: Not applicable. Main Outcome Measure: The degree of synchronization between motor unit firing was evaluated by identifying peaks in the electromyographic averages of the vastus lateralis, triggered from motor unit action potentials in the VMO, and the proportion of power in the power spectral density of the triggered average at the firing frequency of the reference motor unit. The proportion of cases in which there was significant power and peaks in the triggered averages was calculated. Results: The proportion of trials with peaks in the triggered averages of the vastus lateralis electromyographic activity was greater than 61.5% in all tasks, and there was a significantly greater proportion of cases where power in the spectrum was greater than 7.5% (P = .01) for the closed chain condition. Conclusions: There was a high proportion of synchronized motor units between the 2 muscles during isometric contractions, with evidence for greater common drive between the VMO and vastus lateralis in closed chain tasks. This has implications for rehabilitation because it suggests that closed chain tasks may generate better coordination between the vasti muscles.

Muscle fiber and motor unit behavior in the longest human skeletal muscle

The sartorius muscle is the longest muscle in the human body. It is strap-like, up to 600 mm in length, and contains five to seven neurovascular compartments, each with a neuromuscular endplate zone. Some of its fibers terminate intrafascicularly, whereas others may run the full length of the muscle. To assess the location and timing of activation within motor units of this long muscle, we recorded electromyographic potentials from multiple intramuscular electrodes along sartorius muscle during steady voluntary contraction and analyzed their activity with spike-triggered averaging from a needle electrode inserted near the proximal end of the muscle. Approximately 30% of sartorius motor units included muscle fibers that ran the full length of the muscle, conducting action potentials at 3.9 +/- 0.1 m/s. Most motor units were innervated within a single muscle endplate zone that was not necessarily near the midpoint of the fiber. As a consequence, action potentials reached the distal end of a unit as late as 100 ms after initiation at an endplate zone. Thus, contractile activity is not synchronized along the length of single sartorius fibers. We postulate that lateral transmission of force from fiber to endomysium and a wide distribution of motor unit endplates along the muscle are critical for the efficient transmission of force from sarcomere to tendon and for the prevention of muscle injury caused by overextension of inactive regions of muscle fibers.

The stimulus-response curve and motor unit variability in normal subjects and subjects with amyotrophic lateral sclerosis

The behavior and stability of motor units (MUs) in response to electrical stimulation of different intensities can be assessed with the stimulus-response curve, which is a graphical representation of the size of the compound muscle action potential (CMAP) in relation to stimulus intensity. To examine MU characteristics across the whole stimulus range, the variability of CMAP responses to electrical stimulation, and the differences that occur between normal and disease states, the curve was studied in 11 normal subjects and 16 subjects with amyotrophic lateral sclerosis (ALS). In normal subjects, the curve showed a gradual increase in CMAP size with increasing stimulus intensity, although one or two discrete steps were sometimes observed in the upper half of the curve, indicating the activation of large MUs at higher intensities. In ALS subjects, large discrete steps, due to loss of MUs and collateral sprouting, were frequently present. Variability of the CMAP responses was greater than baseline variability, indicating variability of MU responses, and at certain levels this variability was up to 100 mu Vms. The stimulus-response curve shows differences between normal and ALS subjects and provides information on MU activation and variability throughout the curve.

Motor unit number estimation - A Bayesian approach

All muscle contractions are dependent on the functioning of motor units. In diseases such as amyotrophic lateral sclerosis (ALS), progressive loss of motor units leads to gradual paralysis. A major difficulty in the search for a treatment for these diseases has been the lack of a reliable measure of disease progression. One possible measure would be an estimate of the number of surviving motor units. Despite over 30 years of motor unit number estimation (MUNE), all proposed methods have been met with practical and theoretical objections. Our aim is to develop a method of MUNE that overcomes these objections. We record the compound muscle action potential (CMAP) from a selected muscle in response to a graded electrical stimulation applied to the nerve. As the stimulus increases, the threshold of each motor unit is exceeded, and the size of the CMAP increases until a maximum response is obtained. However, the threshold potential required to excite an axon is not a precise value but fluctuates over a small range leading to probabilistic activation of motor units in response to a given stimulus. When the threshold ranges of motor units overlap, there may be alternation where the number of motor units that fire in response to the stimulus is variable. This means that increments in the value of the CMAP correspond to the firing of different combinations of motor units. At a fixed stimulus, variability in the CMAP, measured as variance, can be used to conduct MUNE using the "statistical" or the "Poisson" method. However, this method relies on the assumptions that the numbers of motor units that are firing probabilistically have the Poisson distribution and that all single motor unit action potentials (MUAP) have a fixed and identical size. These assumptions are not necessarily correct. We propose to develop a Bayesian statistical methodology to analyze electrophysiological data to provide an estimate of motor unit numbers. Our method of MUNE incorporates the variability of the threshold, the variability between and within single MUAPs, and baseline variability. Our model not only gives the most probable number of motor units but also provides information about both the population of units and individual units. We use Markov chain Monte Carlo to obtain information about the characteristics of individual motor units and about the population of motor units and the Bayesian information criterion for MUNE. We test our method of MUNE on three subjects. Our method provides a reproducible estimate for a patient with stable but severe ALS. In a serial study, we demonstrate a decline in the number of motor unit numbers with a patient with rapidly advancing disease. Finally, with our last patient, we show that our method has the capacity to estimate a larger number of motor units.

The relationship between Bayesian motor unit number estimation and histological measurements of motor neurons in wild-type and SOD1 G93A mice

Objective: To assess the relationship between Bayesian MUNE and histological motor neuron counts in wild-type mice and in an animal model of ALS. Methods: We performed Bayesian MUNE paired with histological counts of motor neurons in the lumbar spinal cord of wild-type mice and transgenic SOD1 G93A mice that show progressive weakness over time. We evaluated the number of acetylcholine endplates that were innervated by a presynaptic nerve. Results: In wild-type mice, the motor unit number in the gastrocnemius muscle estimated by Bayesian MUNE was approximately half the number of motor neurons in the region of the spinal cord that contains the cell bodies of the motor neurons supplying the hindlimb crural flexor muscles. In SOD1 G93A mice, motor neuron numbers declined over time. This was associated with motor endplate denervation at the end-stage of disease. Conclusion: The number of motor neurons in the spinal cord of wild-type mice is proportional to the number of motor units estimated by Bayesian MUNE. In SOD1 G93A mice, there is a lower number of estimated motor units compared to the number of spinal cord motor neurons at the end-stage of disease, and this is associated with disruption of the neuromuscular junction. Significance: Our finding that the Bayesian MUNE method gives estimates of motor unit numbers that are proportional to the numbers of motor neurons in the spinal cord supports the clinical use of Bayesian MUNE in monitoring motor unit loss in ALS patients. © 2012 International Federation of Clinical Neurophysiology.

Quantitative studies of lower motor neuron degeneration in amyotrophic lateral sclerosis : Evidence for exponential decay of motor unit numbers and greatest rate of loss at the site of onset

Objective: To use our Bayesian method of motor unit number estimation (MUNE) to evaluate lower motor neuron degeneration in ALS. Methods: In subjects with ALS we performed serial MUNE studies. We examined the repeatability of the test and then determined whether the loss of MUs was fitted by an exponential or Weibull distribution. Results: The decline in motor unit (MU) numbers was well-fitted by an exponential decay curve. We calculated the half life of MUs in the abductor digiti minimi (ADM), abductor pollicis brevis (APB) and/or extensor digitorum brevis (EDB) muscles. The mean half life of the MUs of ADM muscle was greater than those of the APB or EDB muscles. The half-life of MUs was less in the ADM muscle of subjects with upper limb than in those with lower limb onset. Conclusions: The rate of loss of lower motor neurons in ALS is exponential, the motor units of the APB decay more quickly than those of the ADM muscle and the rate of loss of motor units is greater at the site of onset of disease. Significance: This shows that the Bayesian MUNE method is useful in following the course and exploring the clinical features of ALS. 2012 International Federation of Clinical Neurophysiology.

Determining feature relevance for the grouping of motor unit action potentials through generative topographic mapping

The study of motor unit action potential (MUAP) activity from electrornyographic signals is an important stage on neurological investigations that aim to understand the state of the neuromuscular system. In this context, the identification and clustering of MUAPs that exhibit common characteristics, and the assessment of which data features are most relevant for the definition of such cluster structure are central issues. In this paper, we propose the application of an unsupervised Feature Relevance Determination (FRD) method to the analysis of experimental MUAPs obtained from healthy human subjects. In contrast to approaches that require the knowledge of a priori information from the data, this FRD method is embedded on a constrained mixture model, known as Generative Topographic Mapping, which simultaneously performs clustering and visualization of MUAPs. The experimental results of the analysis of a data set consisting of MUAPs measured from the surface of the First Dorsal Interosseous, a hand muscle, indicate that the MUAP features corresponding to the hyperpolarization period in the physisiological process of generation of muscle fibre action potentials are consistently estimated as the most relevant and, therefore, as those that should be paid preferential attention for the interpretation of the MUAP groupings.