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.

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.

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.

Generative topographic mapping applied to clustering and visualization of motor unit action potentials

The identification and visualization of clusters formed by motor unit action potentials (MUAPs) is an essential step in investigations seeking to explain the control of the neuromuscular system. This work introduces the generative topographic mapping (GTM), a novel machine learning tool, for clustering of MUAPs, and also it extends the GTM technique to provide a way of visualizing MUAPs. The performance of GTM was compared to that of three other clustering methods: the self-organizing map (SOM), a Gaussian mixture model (GMM), and the neural-gas network (NGN). The results, based on the study of experimental MUAPs, showed that the rate of success of both GTM and SOM outperformed that of GMM and NGN, and also that GTM may in practice be used as a principled alternative to the SOM in the study of MUAPs. A visualization tool, which we called GTM grid, was devised for visualization of MUAPs lying in a high-dimensional space. The visualization provided by the GTM grid was compared to that obtained from principal component analysis (PCA). (c) 2005 Elsevier Ireland Ltd. All rights reserved.

Motor-unit Coherence During Isometric Contractions Is Greater in a Hand Muscle of Older Adults

The purpose of this study was to quantify the strength of motor-unit coherence from the first dorsal interosseus muscle in young and old adults using data obtained in a previous study, where no differences in motor-unit synchronization between the two groups were observed.

Motor-unit coherence and its relation with synchrony are influenced by training

The purpose of the study was to quantify the strength of motor-unit coherence from the left and right first dorsal interosseous muscles in untrained, skill-trained (musicians), and strength-trained (weightlifters) individuals who had long-term specialized use of their hand muscles. The strength of motor-unit coherence was quantified from a total of 394 motor-unit pairs in 13 subjects using data from a previous study in which differences were found in the strength of motor-unit synchronization depending on training status. In the present study, we found that the strength of motor-unit coherence was significantly greater in the left compared with the right hand of untrained right-handed subjects with the largest differences observed between 21 and 24 Hz. The strength of motor-unit coherence was lower in both hands of skill-trained subjects (21–27 Hz) and the right (skilled) hand of untrained subjects (21–24 Hz), whereas the largest motor-unit coherence was observed in both hands of strength-trained subjects (3–9 and 21–27 Hz). A strong curvilinear association was observed between motor-unit synchronization and the integral of coherence at 10–30 Hz in all motor-unit pairs (r² = 0.77), and was most pronounced in strength-trained subjects (r² = 0.90). Furthermore, this association was accentuated when using synchronization data with broad peaks (>11 ms), suggesting that the 10- to 30-Hz coherence is due to oscillatory activity in indirect branched common inputs. The altered coherence with training may be due to an interaction between cortical inhibition and the number of direct common inputs to motor neurons in skill- or strength-trained hands.