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.

Mechanisms that contribute to differences in motor performance between young and old adults

This paper examines the physiological mechanisms responsible for differences in the amplitude of force fluctuations between young and old adults. Because muscle force is a consequence of motor unit activity, the potential mechanisms include both motor unit properties and the behavior of motor unit populations. The force fluctuations, however, depend not only on the age of the individual but also on the muscle group performing the task, the type and intensity of the muscle contraction, and the physical activity status of the individual. Computer simulations and experimental findings performed on tasks that involved single agonist and antagonist muscles suggest that differences in force fluctuations are not attributable to motor unit twitch force, motor unit number, or nonuniform activation of the agonist muscle, but that they are influenced by the variability and common modulation of motor unit discharge in both the agonist and antagonist muscles. Because the amplitude of the force fluctuations does not vary linearly with muscle activation, these results suggest that multiple mechanisms contribute to the differences in force fluctuations between young and old adults, although the boundary conditions for each mechanism remain to be determined.

Corticospinal output and loss of force during motor fatigue

The objective of this study was to analyze central motor output changes in relation to contraction force during motor fatigue. The triple stimulation technique (TST, Magistris et al. in Brain 121(Pt 3):437-450, 1998) was used to quantify a central conduction index (CCI = amplitude ratio of central conduction response and peripheral nerve response, obtained simultaneously by the TST). The CCI removes effects of peripheral fatigue from the quantification. It allows a quantification of the percentage of the entire target muscle motor unit pool driven to discharge by a transcranial magnetic stimulus. Subjects (n = 23) performed repetitive maximal voluntary contractions (MVC) of abductor digiti minimi (duration 1 s, frequency 0.5 Hz) during 2 min. TST recordings were obtained every 15 s, using stimulation intensities sufficient to stimulate all cortical motor neurons (MNs) leading to the target muscle, and during voluntary contractions of 20% of the MVC to facilitate the responses. TST was also repetitively recorded during recovery. This basic exercise protocol was modified in a number of experiments to further characterize influences on CCI of motor fatigue (4 min exercise at 50% MVC; delayed fatigue recovery during local hemostasis, "stimulated exercise" by 20 Hz trains of 1 s duration at 0.5 Hz during 2 min). In addition, the cortical silent period was measured during the basic exercise protocol. Force fatigued to approximately 40% of MVC in all experiments and in all subjects. In all subjects, CCI decreased during exercise, but this decrease varied markedly between subjects. On average, CCI reductions preceded force reductions during exercise, and CCI recovery preceded force recovery. Exercising at 50% for 4 min reduced muscle force more markedly than CCI. Hemostasis induced by a cuff delayed muscle force recovery, but not CCI recovery. Stimulated exercise reduced force markedly, but CCI decreased only marginally. Summarized, force reduction and reduction of the CCI related poorly quantitatively and in time, and voluntary drive was particularly critical to reduce the CCI. The fatigue induced reduction of CCI may result from a central inhibitory phenomenon. Voluntary muscle activation is critical for the CCI reduction, suggesting a primarily supraspinal mechanism.

Identificação de sistemas para o estudo de controle motor.

Qualquer tarefa motora ativa se dá pela ativação de uma população de unidades motoras. Porém, devido a diversas dificuldades, tanto técnicas quanto éticas, não é possível medir a entrada sináptica dos motoneurônios em humanos. Por essas razões, o uso de modelos computacionais realistas de um núcleo de motoneurônios e as suas respectivas fibras musculares tem um importante papel no estudo do controle humano dos músculos. Entretanto, tais modelos são complexos e uma análise matemática é difícil. Neste texto é apresentada uma abordagem baseada em identificação de sistemas de um modelo realista de um núcleo de unidades motoras, com o objetivo de obter um modelo mais simples capaz de representar a transdução das entradas do núcleo de unidades motoras na força do músculo associado ao núcleo. A identificação de sistemas foi baseada em um algoritmo de mínimos quadrados ortogonal para achar um modelo NARMAX, sendo que a entrada considerada foi a condutância sináptica excitatória dendrítica total dos motoneurônios e a saída foi a força dos músculos produzida pelo núcleo de unidades motoras. O modelo identificado reproduziu o comportamento médio da saída do modelo computacional realista, mesmo para pares de sinal de entrada-saída não usados durante o processo de identificação do modelo, como sinais de força muscular modulados senoidalmente. Funções de resposta em frequência generalizada do núcleo de motoneurônios foram obtidas do modelo NARMAX, e levaram a que se inferisse que oscilações corticais na banda-beta (20 Hz) podem influenciar no controle da geração de força pela medula espinhal, comportamento do núcleo de motoneurônios até então desconhecido.

Electromyographic activity in lower limb muscles is temporally associated with the slow phase of oxygen uptake during cycling

Although the "slow" phase of pulmonary oxygen uptake (Vo2) appears to represent energetic processes in contracting muscle, electromyographic evidence tends not to support this. The present study assessed normalized integrated electromyographic (NIEMG) activity in eight muscles that act about the hip, knee and ankle during 8 min of moderate (ventilatory threshold) cycling in six male cyclists. (Vo2) was measured breath by breath during four repeated trials at each of the two intensities. Moderate and very heavy exercise followed a 4-min period of light exercise (50 W). During moderate exercise the slow (Vo2) phase was absent and NIEMG in all muscles did not increase after the first minute of exercise. During very heavy exercise, the slow phase emerged (time delay=58 ± 16 s) and increased progressively (time constant=120 ± 35 s) to an amplitude (0.83 ± 0.16 L/min) that was approximately 21% of the total (Vo2) response. This slow (Vo2) phase coincided with a significant increase in NIEMG in most muscles, and differences in NIEMG activities between the two intensities revealed "slow" muscle activation profiles that differed between muscles in terms of the onset, amplitude and shape of these profiles. This supports the hypothesis that the slow (Vo2) phase is a function of these different slow muscle activation profiles.

Oxygen consumption and muscle activation patterns during constant load exercise

The collective purpose of these two studies was to determine a link between the V02 slow component and the muscle activation patterns that occur during cycling. Six, male subjects performed an incremental cycle ergometer exercise test to determine asub-TvENT (i.e. 80% of TvENT) and supra-TvENT (TvENT + 0.75*(V02 max - TvENT) work load. These two constant work loads were subsequently performed on either three or four occasions for 8 mins each, with V02 captured on a breath-by-breath basis for every test, and EMO of eight major leg muscles collected on one occasion. EMG was collected for the first 10 s of every 30 s period, except for the very first 10 s period. The V02 data was interpolated, time aligned, averaged and smoothed for both intensities. Three models were then fitted to the V02 data to determine the kinetics responses. One of these models was mono-exponential, while the other two were biexponential. A second time delay parameter was the only difference between the two bi-exponential models. An F-test was used to determine significance between the biexponential models using the residual sum of squares term for each model. EMO was integrated to obtain one value for each 10 s period, per muscle. The EMG data was analysed by a two-way repeated measures ANOV A. A correlation was also used to determine significance between V02 and IEMG. The V02 data during the sub-TvENT intensity was best described by a mono-exponential response. In contrast, during supra-TvENT exercise the two bi-exponential models best described the V02 data. The resultant F-test revealed no significant difference between the two models and therefore demonstrated that the slow component was not delayed relative to the onset of the primary component. Furthermore, only two parameters were deemed to be significantly different based upon the two models. This is in contrast to other findings. The EMG data, for most muscles, appeared to follow the same pattern as V02 during both intensities of exercise. On most occasions, the correlation coefficient demonstrated significance. Although some muscles demonstrated the same relative increase in IEMO based upon increases in intensity and duration, it cannot be assumed that these muscles increase their contribution to V02 in a similar fashion. Larger muscles with a higher percentage of type II muscle fibres would have a larger increase in V02 over the same increase in intensity.

Bayesian algorithms with applications

Advances in algorithms for approximate sampling from a multivariable target function have led to solutions to challenging statistical inference problems that would otherwise not be considered by the applied scientist. Such sampling algorithms are particularly relevant to Bayesian statistics, since the target function is the posterior distribution of the unobservables given the observables. In this thesis we develop, adapt and apply Bayesian algorithms, whilst addressing substantive applied problems in biology and medicine as well as other applications. For an increasing number of high-impact research problems, the primary models of interest are often sufficiently complex that the likelihood function is computationally intractable. Rather than discard these models in favour of inferior alternatives, a class of Bayesian "likelihoodfree" techniques (often termed approximate Bayesian computation (ABC)) has emerged in the last few years, which avoids direct likelihood computation through repeated sampling of data from the model and comparing observed and simulated summary statistics. In Part I of this thesis we utilise sequential Monte Carlo (SMC) methodology to develop new algorithms for ABC that are more efficient in terms of the number of model simulations required and are almost black-box since very little algorithmic tuning is required. In addition, we address the issue of deriving appropriate summary statistics to use within ABC via a goodness-of-fit statistic and indirect inference. Another important problem in statistics is the design of experiments. That is, how one should select the values of the controllable variables in order to achieve some design goal. The presences of parameter and/or model uncertainty are computational obstacles when designing experiments but can lead to inefficient designs if not accounted for correctly. The Bayesian framework accommodates such uncertainties in a coherent way. If the amount of uncertainty is substantial, it can be of interest to perform adaptive designs in order to accrue information to make better decisions about future design points. This is of particular interest if the data can be collected sequentially. In a sense, the current posterior distribution becomes the new prior distribution for the next design decision. Part II of this thesis creates new algorithms for Bayesian sequential design to accommodate parameter and model uncertainty using SMC. The algorithms are substantially faster than previous approaches allowing the simulation properties of various design utilities to be investigated in a more timely manner. Furthermore the approach offers convenient estimation of Bayesian utilities and other quantities that are particularly relevant in the presence of model uncertainty. Finally, Part III of this thesis tackles a substantive medical problem. A neurological disorder known as motor neuron disease (MND) progressively causes motor neurons to no longer have the ability to innervate the muscle fibres, causing the muscles to eventually waste away. When this occurs the motor unit effectively ‘dies’. There is no cure for MND, and fatality often results from a lack of muscle strength to breathe. The prognosis for many forms of MND (particularly amyotrophic lateral sclerosis (ALS)) is particularly poor, with patients usually only surviving a small number of years after the initial onset of disease. Measuring the progress of diseases of the motor units, such as ALS, is a challenge for clinical neurologists. Motor unit number estimation (MUNE) is an attempt to directly assess underlying motor unit loss rather than indirect techniques such as muscle strength assessment, which generally is unable to detect progressions due to the body’s natural attempts at compensation. Part III of this thesis builds upon a previous Bayesian technique, which develops a sophisticated statistical model that takes into account physiological information about motor unit activation and various sources of uncertainties. More specifically, we develop a more reliable MUNE method by applying marginalisation over latent variables in order to improve the performance of a previously developed reversible jump Markov chain Monte Carlo sampler. We make other subtle changes to the model and algorithm to improve the robustness of the approach.

Is recovery driven by central or peripheral factors? A role for the brain in recovery following intermittent-sprint exercise

Prolonged intermittent-sprint exercise (i.e., team sports) induce disturbances in skeletal muscle structure and function that are associated with reduced contractile function, a cascade of inflammatory responses, perceptual soreness, and a delayed return to optimal physical performance. In this context, recovery from exercise-induced fatigue is traditionally treated from a peripheral viewpoint, with the regeneration of muscle physiology and other peripheral factors the target of recovery strategies. The direction of this research narrative on post-exercise recovery differs to the increasing emphasis on the complex interaction between both central and peripheral factors regulating exercise intensity during exercise performance. Given the role of the central nervous system (CNS) in motor-unit recruitment during exercise, it too may have an integral role in post-exercise recovery. Indeed, this hypothesis is indirectly supported by an apparent disconnect in time-course changes in physiological and biochemical markers resultant from exercise and the ensuing recovery of exercise performance. Equally, improvements in perceptual recovery, even withstanding the physiological state of recovery, may interact with both feed-forward/feed-back mechanisms to influence subsequent efforts. Considering the research interest afforded to recovery methodologies designed to hasten the return of homeostasis within the muscle, the limited focus on contributors to post-exercise recovery from CNS origins is somewhat surprising. Based on this context, the current review aims to outline the potential contributions of the brain to performance recovery after strenuous exercise.

Bayesian P-splines with a multiplicative term in EMG trace data

A method is proposed to describe force or compound muscle action potential (CMAP) trace data collected in an electromyography study for motor unit number estimation (MUNE). Experimental data was collected using incre- mental stimulation at multiple durations. However, stimulus information, vital for alternate MUNE methods, is not comparable for multiple duration data and therefore previous methods of MUNE (Ridall et al., 2006, 2007) cannot be used with any reliability. Hypothesised ring combinations of motor units are mod- elled using a multiplicative factor and Bayesian P-spline formulation. The model describes the process for force and CMAP in a meaningful way.

Statistical methods for electromyography data and associated problems

This thesis proposes three novel models which extend the statistical methodology for motor unit number estimation, a clinical neurology technique. Motor unit number estimation is important in the treatment of degenerative muscular diseases and, potentially, spinal injury. Additionally, a recent and untested statistic to enable statistical model choice is found to be a practical alternative for larger datasets. The existing methods for dose finding in dual-agent clinical trials are found to be suitable only for designs of modest dimensions. The model choice case-study is the first of its kind containing interesting results using so-called unit information prior distributions.