Use of Bayesian MUNE to show differing rate of loss of motor units in subgroups of ALS

Objectives To evaluate differences among patients with different clinical features of ALS, we used our Bayesian method of motor unit number estimation (MUNE). Methods We performed serial MUNE studies on 42 subjects who fulfilled the diagnostic criteria for ALS during the course of their illness. Subjects were classified into three subgroups according to whether they had typical ALS (with upper and lower motor neurone signs) or had predominantly upper motor neurone weakness with only minor LMN signs, or predominantly lower motor neurone weakness with only minor UMN signs. In all subjects we calculated the half life of MUs, defined as the expected time for the number of MUs to halve, in one or more of the abductor digiti minimi (ADM), abductor pollicis brevis (APB) and extensor digitorum brevis (EDB) muscles. Results The mean half life of MUs was less in subjects who had typical ALS with both upper and lower motor neurone signs than in those with predominantly upper motor neurone weakness or predominantly lower motor neurone weakness. In 18 subjects we analysed the estimated size of the MUs and demonstrated the appearance of large MUs in subjects with upper or lower motor neurone predominant weakness. We found that the appearance of large MUs was correlated with the half life of MUs. Conclusions Patients with different clinical features of ALS have different rates of loss and different sizes of MUs. Significance: These findings could indicate differences in disease pathogenesis.

Estimating the number of motor units using random sums with independently thinned terms

The problem of estimating the numbers of motor units N in a muscle is embedded in a general stochastic model using the notion of thinning from point process theory. In the paper a new moment type estimator for the numbers of motor units in a muscle is denned, which is derived using random sums with independently thinned terms. Asymptotic normality of the estimator is shown and its practical value is demonstrated with bootstrap and approximative confidence intervals for a data set from a 31-year-old healthy right-handed, female volunteer. Moreover simulation results are presented and Monte-Carlo based quantiles, means, and variances are calculated for N in{300,600,1000}.

Central and peripheral defects in motor units of the diaphragm of spinal muscular atrophy mice.

Spinal muscular atrophy (SMA) is characterized by motoneuron loss and muscle weakness. However, the structural and functional deficits that lead to the impairment of the neuromuscular system remain poorly defined. By electron microscopy, we previously found that neuromuscular junctions (NMJs) and muscle fibres of the diaphragm are among the earliest affected structures in the severe mouse SMA model. Because of certain anatomical features, i.e. its thinness and its innervation from the cervical segments of the spinal cord, the diaphragm is particularly suitable to characterize both central and peripheral events. Here we show by immunohistochemistry that, at postnatal day 3, the cervical motoneurons of SMA mice receive less stimulatory synaptic inputs. Moreover, their mitochondria become less elongated which might represent an early stage of degeneration. The NMJs of the diaphragm of SMA mice show a loss of synaptic vesicles and active zones. Moreover, the partly innervated endplates lack S100 positive perisynaptic Schwann cells (PSCs). We also demonstrate the feasibility of comparing the proteomic composition between diaphragm regions enriched and poor in NMJs. By this approach we have identified two proteins that are significantly upregulated only in the NMJ-specific regions of SMA mice. These are apoptosis inducing factor 1 (AIFM1), a mitochondrial flavoprotein that initiates apoptosis in a caspase-independent pathway, and four and a half Lim domain protein 1 (FHL1), a regulator of skeletal muscle mass that has been implicated in several myopathies.

Discharge patterns of human tensor palatini motor units during sleep onset

Study Objectives: Upper airway muscles such as genioglossus (GG) and tensor palatini (TP) reduce activity at sleep onset. In GG reduced muscle activity is primarily due to inspiratory modulated motor units becoming silent, suggesting reduced respiratory pattern generator (RPG) output. However, unlike GG, TP shows minimal respiratory modulation and presumably has few inspiratory modulated motor units and minimal input from the RPG. Thus, we investigated the mechanism by which TP reduces activity at sleep onset.

Design: The activity of TP motor units were studied during relaxed wakefulness and over the transition from wakefulness to sleep.

Setting: Sleep laboratory.

Participants: Nine young (21.4 ± 3.4 years) males were studied on a total of 11 nights.

Intervention: Sleep onset.

Measurements and Results: Two TP EMGs (thin, hooked wire electrodes), and sleep and respiratory measures were recorded. One hundred twenty-one sleep onsets were identified (13.4 ± 7.2/subject), resulting in 128 motor units (14.3 ± 13.0/subject); 29% of units were tonic, 43% inspiratory modulated (inspiratory phasic 18%, inspiratory tonic 25%), and 28% expiratory modulated (expiratory phasic 21%, expiratory tonic 7%). There was a reduction in both expiratory and inspiratory modulated units, but not tonic units, at sleep onset. Reduced TP activity was almost entirely due to de-recruitment.

Conclusions: TP showed a similar distribution of motor units as other airway muscles. However, a greater proportion of expiratory modulated motor units were active in TP and these expiratory units, along with inspiratory units, tended to become silent over sleep onset. The data suggest that both expiratory and inspiratory drive components from the RPG are reduced at sleep onset in TP.

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).

Motor unit synchronization measured by cross-correlation is not influenced by short-term strength training of a hand muscle

The purpose of the study was to quantify the strength of motor unit synchronization and coherence from pairs of concurrently active motor units before and after short-term (4–8 weeks) strength training of the left first dorsal interosseous (FDI) muscle. Five subjects (age 24.8 ± 4.3 years) performed a training protocol three times/week that consisted of six sets of ten maximal isometric index finger abductions, whereas three subjects (age 27.3 ± 6.7 years) acted as controls. Motor unit activity was recorded from pairs of intramuscular electrodes in the FDI muscle with two separate motor unit recording sessions obtained before and after strength training (trained group) or after 4 weeks of normal daily activities that did not involve training (control group). The training intervention resulted in a 54% (45.2 ± 8.3 to 69.5 ± 13.8 N, P = 0.001) increase in maximal index finger abduction force, whereas there was no change in strength in the control group. A total of 163 motor unit pairs (198 single motor units) were examined in both subject groups, with 52 motor unit pairs obtained from 10 recording sessions before training and 51 motor unit pairs from 10 recording sessions after training. Using the cross-correlation procedure, there was no change in the strength of motor unit synchronization following strength training (common input strength index; 0.71 ± 0.41 to 0.67 ± 0.43 pulses/s). Furthermore, motor unit coherence z scores at low (0–10 Hz; 3.9 ± 0.3 before to 4.4 ± 0.4 after) or high (10–30 Hz; 1.7 ± 0.1 before to 1.9 ± 0.1 after) frequencies were not influenced by strength training. These motor unit data indicate that increases in strength following several weeks of training a hand muscle are not accompanied by changes in motor unit synchronization or coherence, suggesting that these features of correlated motor unit activity are not important in the expression of muscle strength.

Properties of low-threshold motor axons in the human median nerve

This study investigated the excitability and accommodative properties of low-threshold human motor axons to test whether these motor axons have greater expression of the persistent Na(+) conductance, I(NaP). Computer-controlled threshold tracking was used to study 22 single motor units and the data were compared with compound motor potentials of various amplitudes recorded in the same experimental session. Detailed comparisons were made between the single units and compound potentials that were 40% or 5% of maximal amplitude, the former because this is the compound potential size used in most threshold tracking studies of axonal excitability, the latter because this is the compound potential most likely to be composed entirely of motor axons with low thresholds to electrical recruitment. Measurements were made of the strength-duration relationship, threshold electrotonus, current-voltage relationship, recovery cycle and latent addition. The findings did not support a difference in I(NaP). Instead they pointed to greater activity of the hyperpolarization-activated inwardly rectifying current (I(h)) as the basis for low threshold to electrical recruitment in human motor axons. Computer modelling confirmed this finding, with a doubling of the hyperpolarization-activated conductance proving the best single parameter adjustment to fit the experimental data. We suggest that the hyperpolarization-activated cyclic nucleotide-gated (HCN) channel(s) expressed on human motor axons may be active at rest and contribute to resting membrane potential.

Motor unit synchronization is reduced in anterior knee pain

Anterior knee pain (AKP) is common and has been argued to be related to poor patellofemoral joint control due to impaired coordination of the vasti muscles. However, there are conflicting data. Changes in motor unit firing may provide more definitive evidence. Synchronization of motor unit action potentials (MUAPs) in vastus medialis obliquus (VMO) and vastus lateralis (VL) may contribute to coordination in patellofemoral joint control. We hypothesized that synchronization may be reduced in AKP. Recordings of single MUAPs were made from VMO and multiunit electromyograph (EMG) recordings were made from VL. Averages of VL EMG recordings were triggered from the single MUAPs in VMO. Motor units in VL firing in association with the VMO motor units would appear as a peak in the VL EMG average. Data were compared to previous normative data. The proportion of trials in which a peak was identified in the triggered averages of VL EMG was reduced in people with AKP (38%) compared to controls (90%). Notably, although 80% of subjects had values less than controls, 20% were within normal limits. These results provide new evidence that motor unit synchronization is modified in the presence of pain and provide evidence for motor control dysfunction in AKP. Perspective: This study shows that coordination of motor units between the medial and lateral vasti muscles in people with anterior knee pain is reduced compared to people without knee pain. It confirms that motor control dysfunction is a factor in this condition and has implications for selection of rehabilitation strategies. (c) 2005 by the American Pain Society.

Motor unit synchronization between medial and lateral vasti muscles

Objective: Accurate neuromuscular control of the patellofemoral joint is important in knee joint mechanics. Strategies to coordinate the vasti muscles, such as motor unit synchronization, may simplify control of patellar tracking. This study investigated motor unit synchronization between vastus medialis (VM) and lateralis (VL). Methods: Electromyographic (EMG) recordings of single motor unit action potentials (MUAPs) were made from VM and single- and multi-unit recordings were made from VL. Synchronization was quantified from peaks in the cross-correlogram generated from single MUAP pairs in VL and VM. The proportion of motor units in VM with synchronized firing in VL was also quantified from peaks in averages of multiunit VL EMG triggered from the VM MUAP. Results: A high degree of synchronization of motor unit firing between VM and VL was identified. Results were similar for cross-correlation (similar to 45% of cases) and triggered averages (similar to 41% of cases). Conclusions: The data suggest that synchronization between VM and VL is higher than expected. Agreement between traditional cross-correlation and triggered averaging methods suggest that this new technique may provide a more clinically viable method to quantify synchronization. Significance: High synchronization between VM and VL may provide a solution to simplify control of the mechanically unstable patellofemoral joint. (c) 2005 International Federation of Clinical Neurophysiology. Published by Elsevier Ireland Ltd. All rights reserved.

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.