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.

Effect of changing data collection parameters on statistical motor unit number estimates

The effect of number of samples and selection of data for analysis on the calculation of surface motor unit potential (SMUP) size in the statistical method of motor unit number estimates (MUNE) was determined in 10 normal subjects and 10 with amyotrophic lateral sclerosis (ALS). We recorded 500 sequential compound muscle action potentials (CMAPs) at three different stable stimulus intensities (10–50% of maximal CMAP). Estimated mean SMUP sizes were calculated using Poisson statistical assumptions from the variance of 500 sequential CMAP obtained at each stimulus intensity. The results with the 500 data points were compared with smaller subsets from the same data set. The results using a range of 50–80% of the 500 data points were compared with the full 500. The effect of restricting analysis to data between 5–20% of the CMAP and to standard deviation limits was also assessed. No differences in mean SMUP size were found with stimulus intensity or use of different ranges of data. Consistency was improved with a greater sample number. Data within 5% of CMAP size gave both increased consistency and reduced mean SMUP size in many subjects, but excluded valid responses present at that stimulus intensity. These changes were more prominent in ALS patients in whom the presence of isolated SMUP responses was a striking difference from normal subjects. Noise, spurious data, and large SMUP limited the Poisson assumptions. When these factors are considered, consistent statistical MUNE can be calculated from a continuous sequence of data points. A 2 to 2.5 SD or 10% window are reasonable methods of limiting data for analysis. Muscle Nerve 27: 320–331, 2003

A MULTILEVEL MODEL TO ADDRESS BATCH EFFECTS IN COPY NUMBER ESTIMATION USING SNP ARRAYS

Submicroscopic changes in chromosomal DNA copy number dosage are common and have been implicated in many heritable diseases and cancers. Recent high-throughput technologies have a resolution that permits the detection of segmental changes in DNA copy number that span thousands of basepairs across the genome. Genome-wide association studies (GWAS) may simultaneously screen for copy number-phenotype and SNP-phenotype associations as part of the analytic strategy. However, genome-wide array analyses are particularly susceptible to batch effects as the logistics of preparing DNA and processing thousands of arrays often involves multiple laboratories and technicians, or changes over calendar time to the reagents and laboratory equipment. Failure to adjust for batch effects can lead to incorrect inference and requires inefficient post-hoc quality control procedures that exclude regions that are associated with batch. Our work extends previous model-based approaches for copy number estimation by explicitly modeling batch effects and using shrinkage to improve locus-specific estimates of copy number uncertainty. Key features of this approach include the use of diallelic genotype calls from experimental data to estimate batch- and locus-specific parameters of background and signal without the requirement of training data. We illustrate these ideas using a study of bipolar disease and a study of chromosome 21 trisomy. The former has batch effects that dominate much of the observed variation in quantile-normalized intensities, while the latter illustrates the robustness of our approach to datasets where as many as 25% of the samples have altered copy number. Locus-specific estimates of copy number can be plotted on the copy-number scale to investigate mosaicism and guide the choice of appropriate downstream approaches for smoothing the copy number as a function of physical position. The software is open source and implemented in the R package CRLMM available at Bioconductor (http:www.bioconductor.org).

USING THE R PACKAGE crlmm FOR GENOTYPING AND COPY NUMBER ESTIMATION

Genotyping platforms such as Affymetrix can be used to assess genotype-phenotype as well as copy number-phenotype associations at millions of markers. While genotyping algorithms are largely concordant when assessed on HapMap samples, tools to assess copy number changes are more variable and often discordant. One explanation for the discordance is that copy number estimates are susceptible to systematic differences between groups of samples that were processed at different times or by different labs. Analysis algorithms that do not adjust for batch effects are prone to spurious measures of association. The R package crlmm implements a multilevel model that adjusts for batch effects and provides allele-specific estimates of copy number. This paper illustrates a workflow for the estimation of allele-specific copy number, develops markerand study-level summaries of batch effects, and demonstrates how the marker-level estimates can be integrated with complimentary Bioconductor software for inferring regions of copy number gain or loss. All analyses are performed in the statistical environment R. A compendium for reproducing the analysis is available from the author’s website (http://www.biostat.jhsph.edu/~rscharpf/crlmmCompendium/index.html).

Preliminary hazards summary report : Bodega Bay Atomic Park unit number 1.

"December 28, 1962."

A large animal model of spinal muscular atrophy and correction of phenotype.

OBJECTIVES Spinal muscular atrophy (SMA) is caused by reduced levels of survival motor neuron (SMN) protein, which results in motoneuron loss. Therapeutic strategies to increase SMN levels including drug compounds, antisense oligonucleotides, and scAAV9 gene therapy have proved effective in mice. We wished to determine whether reduction of SMN in postnatal motoneurons resulted in SMA in a large animal model, whether SMA could be corrected after development of muscle weakness, and the response of clinically relevant biomarkers. METHODS Using intrathecal delivery of scAAV9 expressing an shRNA targeting pig SMN1, SMN was knocked down in motoneurons postnatally to SMA levels. This resulted in an SMA phenotype representing the first large animal model of SMA. Restoration of SMN was performed at different time points with scAAV9 expressing human SMN (scAAV9-SMN), and electrophysiology measurements and pathology were performed. RESULTS Knockdown of SMN in postnatal motoneurons results in overt proximal weakness, fibrillations on electromyography indicating active denervation, and reduced compound muscle action potential (CMAP) and motor unit number estimation (MUNE), as in human SMA. Neuropathology showed loss of motoneurons and motor axons. Presymptomatic delivery of scAAV9-SMN prevented SMA symptoms, indicating that all changes are SMN dependent. Delivery of scAAV9-SMN after symptom onset had a marked impact on phenotype, electrophysiological measures, and pathology. INTERPRETATION High SMN levels are critical in postnatal motoneurons, and reduction of SMN results in an SMA phenotype that is SMN dependent. Importantly, clinically relevant biomarkers including CMAP and MUNE are responsive to SMN restoration, and abrogation of phenotype can be achieved even after symptom onset.

Assessment of disease progression in motor neuron disease

Motor neuron disease (MND) is characterised by progressive deterioration of the corticospinal tract, brainstem, and anterior horn cells of the spinal cord. There is no pathognomonic test for the diagnosis of MND, and physicians rely on clinical criteria-upper and lower motor neuron signs-for diagnosis. The presentations, clinical phenotypes, and outcomes of MND are diverse and have not been combined into a marker of disease progression. No single algorithm combines the findings of functional assessments and rating scales, such as those that assess quality of life, with biological markers of disease activity and findings from imaging and neurophysiological assessments. Here, we critically appraise developments in each of these areas and discuss the potential of such measures to be included in the future assessment of disease progression in patients with MND.

Estimation of vegetation carbon stock in Portugal using land use/ land cover data

Dissertation submitted in partial fulfilment of the requirements for the Degree of Master of Science in Geospatial Technologies

Gaussian process framework for pervasive estimation of swimming velocity with body-worn IMU

Presented is an accurate swimming velocity estimation method using an inertial measurement unit (IMU) by employing a simple biomechanical constraint of motion along with Gaussian process regression to deal with sensor inherent errors. Experimental validation shows a velocity RMS error of 9.0 cm/s and high linear correlation when compared with a commercial tethered reference system. The results confirm the practicality of the presented method to estimate swimming velocity using a single low-cost, body-worn IMU.

Sequential sampling of Euschistus heros (Heteroptera: Pentatomidae) in soybean

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Gene number in an invertebrate chordate, Ciona intestinalis

Gene number can be considered a pragmatic measure of biological complexity, but reliable data is scarce. Estimates for vertebrates are 50-100,000 genes per haploid genome, whereas invertebrate estimates fall below 25,000. We wished to test the hypothesis that the origin of vertebrates coincided with extensive gene creation. A prediction is that gene number will differ sharply between invertebrate and vertebrate members of the chordate phylum. A gene number estimation method requiring limited sequence sampling of genomic DNA was developed and validated by using data for Caenorhabditis elegans. Using the method, we estimated that the invertebrate chordate Ciona intestinalis has 15,500 protein-coding genes (±3,700). This number is significantly lower than gene numbers of vertebrate chordates, but similar to those of invertebrates in distantly related phyla. The data indicate that evolution of vertebrates was accompanied by a dramatic increase in protein-coding capacity of the genome.

Simultaneous A8344G heteroplasmy and mitochondrial DNA copy number quantification in Myoclonus Epilepsy and Ragged-Red Fibers (MERRF) syndrome by a multiplex Molecular Beacon based real-time fluorescence PCR

The association of a particular mitochondrial DNA (mtDNA) mutation with different clinical phenotypes is a well-known feature of mitochondrial diseases. A simple genotype–phenotype correlation has not been found between mutation load and disease expression. Tissue and intercellular mosaicism as well as mtDNA copy number are thought to be responsible for the different clinical phenotypes. As disease expression of mitochondrial tRNA mutations is mostly in postmitotic tissues, studies to elucidate disease mechanisms need to be performed on patient material. Heteroplasmy quantitation and copy number estimation using small patient biopsy samples has not been reported before, mainly due to technical restrictions. In order to resolve this problem, we have developed a robust assay that utilizes Molecular Beacons to accurately quantify heteroplasmy levels and determine mtDNA copy number in small samples carrying the A8344G tRNALys mutation. It provides the methodological basis to investigate the role of heteroplasmy and mtDNA copy number in determining the clinical phenotypes.

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.

Transpecific microsatellites for hard pines M. Shepherd, M. Cross, T. Maguire, M. Dieters, C. Williams, R. Henry

Microsatellites are difficult to recover from large plant genomes so cross-specific utilisation is an important source of markers. Fifty micro satellites were tested for cross-specific amplification and polymorphism to two New World hard pine species, slash pine (Pinus elliottii var. elliottii) and Caribbean pine (R caribaea var. hondurensis). Twenty-nine (58%) markers amplified in both hard pine species, and 23 of these 29 were polymorphic. Soft pine (subgenus Strobus) microsatellite markers did amplify, but none were polymorphic. Pinus elliottii var. elliottii and R caribaea var. hondurensis showed mutational changes in the flanking regions and the repeat motif that were informative for Pinus spp. phylogenetic relationships. Most allele length variation could be attributed to variability in repeat unit number. There was no evidence for ascertainment bias.

Lack of the transcription factor hypoxia-inducible factor (HIF)-1α in macrophages accelerates the necrosis of Mycobacterium avium-induced granulomas

Accepted Manuscript