A simulation of the surface EMG for analysis of muscle activity during whole body vibratory stimulation

Many studies have accounted for whole body vibration effects in the fields of exercise physiology, sport and rehabilitation medicine. Generally, surface EMG is utilized to assess muscular activity during the treatment; however, large motion artifacts appear superimposed to the raw signal, making sEMG recording not suitable before any artifact filtering. Sharp notch filters, centered at vibration frequency and at its superior harmonics, have been used in previous studies, to remove the artifacts. [6, 10] However, to get rid of those artifacts some true EMG signal is lost. The purpose of this study was to reproduce the effect of motor-unit synchronization on a simulated surface EMG during vibratory stimulation. In addition, authors mean to evaluate the EMG power percentage in those bands in which are also typically located motion artifact components. Model characteristics were defined to take into account two main aspect: the muscle MUs discharge behavior and the triggering effects that appear during local vibratory stimulation. [7] Inter-pulse-interval, was characterized by a polimodal distribution related to the MU discharge frequency (IPI 55-80ms, σ=12ms) and to the correlation with the vibration period within the range of ±2 ms due to vibration stimulus. [1, 7] The signals were simulated using different stimulation frequencies from 30 to 70 Hz. The percentage of the total simulated EMG power within narrow bands centered at the stimulation frequency and its superior harmonics (± 1 Hz) resulted on average about 8% (± 2.85) of the total EMG power. However, the artifact in those bands may contain more than 40% of the total power of the total signal. [6] Our preliminary results suggest that the analysis of the muscular activity of muscle based on raw sEMG recordings and RMS evaluation, if not processed during vibratory stimulation may lead to a serious overestimation of muscular response.

Simulation of surface EMG for the analysis of muscle activity during whole body vibratory stimulation

This study aims to reproduce the effect of motor-unit synchronization on surface EMG recordings during vibratory stimulation to highlight vibration evoked muscle activity. The authors intended to evaluate, through numerical simulations, the changes in surface EMG spectrum in muscles undergoing whole body vibration stimulation. In some specific bands, in fact, vibration induced motion artifacts are also typically present. In addition, authors meant to compare the simulated EMGs with respect to real recordings in order to discriminate the effect of synchronization of motor units discharges with vibration frequencies from motion artifacts. Computations were performed using a model derived from previous studies and modified to consider the effect of vibratory stimulus, the motor unit synchronization and the endplates-electrodes relative position on the EMG signal. Results revealed that, in particular conditions, synchronization of MUs' discharge generates visible peaks at stimulation frequency and its harmonics. However, only a part of the total power of surface EMGs might be enclosed within artifacts related bands (±1. Hz centered at the stimulation frequency and its superior harmonics) even in case of strong synchronization of motor units discharges with the vibratory stimulus. © 2013 Elsevier Ireland Ltd.

Caracterização da fosfatase String/Cdc25 como um regulador da neurotoxicidade de Tau em Drosophila

As tauopatias, grupo onde se inclui a doença de Alzheimer (AD), são caracterizadas pela deposição intracelular de emaranhados neurofibrilares (NFTs), compostos principalmente por formas hiperfosforiladas da proteína Tau, uma proteína que se associa aos microtúbulos. Os mecanismos moleculares subjacentes à neurotoxicidade induzida por Tau não são ainda claros. Drosophila melanogaster tem sido usada para modelar diversas doenças neurodegenerativas humanas, incluindo as tauopatias. Neste trabalho foi usado o sistema visual de Drosophila como modelo para identificar os passos que podem levar à acumulação de Tau em Tauopatias. Durante o desenvolvimento do olho de Drosophila, a expressão ectópica de hTau induz um olho rugoso, em consequência da neurotoxicidade, e que pode ser utilizado para identificar modificadores do fenótipo. A fosfatase codificada por string /cdc25 (stg), um regulador universal da transição G2/M, foi previamente identificada como um supressor da neurotoxicidade associada à expressão da proteina Tau. No entanto, os mecanismos moleculares que estão na base desta interação genética nunca foram estudados, desconhecendo-se também se a atividade fosfatase de Stg/Cdc25 é essencial para modular os níveis de fosforilação de Tau. O objetivo deste projeto consistiu em elucidar os mecanismos que se encontram na base da interação Stg-Tau. Para alcançar este objectivo, usou-se uma abordagem genética e bioquímica. Os resultados obtidos sugerem que Stg é um possível modulador da neurotoxicidade de Tau.

Tensor Completion for Multidimensional Inverse Problems with Applications to Magnetic Resonance Relaxometry

This thesis deals with tensor completion for the solution of multidimensional inverse problems. We study the problem of reconstructing an approximately low rank tensor from a small number of noisy linear measurements. New recovery guarantees, numerical algorithms, non-uniform sampling strategies, and parameter selection algorithms are developed. We derive a fixed point continuation algorithm for tensor completion and prove its convergence. A restricted isometry property (RIP) based tensor recovery guarantee is proved. Probabilistic recovery guarantees are obtained for sub-Gaussian measurement operators and for measurements obtained by non-uniform sampling from a Parseval tight frame. We show how tensor completion can be used to solve multidimensional inverse problems arising in NMR relaxometry. Algorithms are developed for regularization parameter selection, including accelerated k-fold cross-validation and generalized cross-validation. These methods are validated on experimental and simulated data. We also derive condition number estimates for nonnegative least squares problems. Tensor recovery promises to significantly accelerate N-dimensional NMR relaxometry and related experiments, enabling previously impractical experiments. Our methods could also be applied to other inverse problems arising in machine learning, image processing, signal processing, computer vision, and other fields.

Smoothing a Discrete Hazard Function for the Number of Patients Colonized with Methicillin-Resistance Staphylococcus Aureus in an Intensive Care Unit

A new method for estimating the time to colonization of Methicillin-resistant Staphylococcus Aureus (MRSA) patients is developed in this paper. The time to colonization of MRSA is modelled using a Bayesian smoothing approach for the hazard function. There are two prior models discussed in this paper: the first difference prior and the second difference prior. The second difference prior model gives smoother estimates of the hazard functions and, when applied to data from an intensive care unit (ICU), clearly shows increasing hazard up to day 13, then a decreasing hazard. The results clearly demonstrate that the hazard is not constant and provide a useful quantification of the effect of length of stay on the risk of MRSA colonization which provides useful insight.

Trends in utilization of computer tomography (CT) scans for victims of motor vehicle collision (MVCs) in a level I trauma center from 1996 to 2006

Purpose. To evaluate trends in the utilization of head, abdominal, thoracic and other body regions CTs in the management of victims of MVC at a level I trauma center from 1996 to 2006.^ Method. From the trauma registry, I identified patients involved in MVC's in a level I trauma center and categorized them into three age groups of 13-18, 19-55 and ≥56. I used International Classification of Disease (ICD-9-CM) codes to find the type and number of CTs examinations performed for each patient. I plotted the mean number of CTs per patient against year of admission to find the crude estimate of change in utilization pattern for each type of CT. I used logistic regression to assess whether repetitive CTs (≥ 2) for head, abdomen, thorax and other body regions were associated with age group and year of admission for MVC patients. I adjusted the estimates for gender, ethnicity, insurance status, mechanism and severity of injury, intensive care unit admission status, patient disposition (dead or alive) and year of admission.^ Results. Utilization of head, abdominal, thoracic and other body regions CTs significantly increased over 11-year period. Utilization of head CT was greatest in the 13-18 age group, and increased from 0.58 CT/patient in 1996 to 1.37 CT/patient in 2006. Abdominal CTs were more common in the ≥56+ age group, and increased from 0.33 CT/patient in 1996 to 0.72 CT/patient in 2006. Utilization of thoracic CTs was higher in the 56+ age group, and increased from 0.01 CT/patient in 1996 to 0.42 CT/patient in 2006. Utilization of other CTs did not change materially during the study period for adolescents, adults or older adults. In the multivariable analysis, after adjustment for potential confounders, repetitive head CTs significantly increased in the 13-18 age group (95% CI: 1.29-1.87, p=<0.001) relative to the 19-55 age group. Repetitive thoracic CT use was lower in adolescents (95% CI: 0.22-0.70, p=<0.001) relative to the 19-55 age group.^ Conclusion. There has been a substantial increase in the utilization of head, abdominal, thoracic and other CTs in the management of MVC patients. Future studies need to identify if increased utilization of CTs have resulted in better health outcome for these patients. ^

Increase in particle number emissions from motor vehicles due to interruption of steady traffic flow

We assess the increase in particle number emissions from motor vehicles driving at steady speed when forced to stop and accelerate from rest. Considering the example of a signalized pedestrian crossing on a two-way single-lane urban road, we use a complex line source method to calculate the total emissions produced by a specific number and mix of light petrol cars and diesel passenger buses and show that the total emissions during a red light is significantly higher than during the time when the light remains green. Replacing two cars with one bus increased the emissions by over an order of magnitude. Considering these large differences, we conclude that the importance attached to particle number emissions in traffic management policies be reassessed in the future.

A new data source for fisheries resource assessment: genetic estimates of the effective number of spawners. Final Report to the Fisheries Research and Development Corporation.

The development of innovative methods of stock assessment is a priority for State and Commonwealth fisheries agencies. It is driven by the need to facilitate sustainable exploitation of naturally occurring fisheries resources for the current and future economic, social and environmental well being of Australia. This project was initiated in this context and took advantage of considerable recent achievements in genomics that are shaping our comprehension of the DNA of humans and animals. The basic idea behind this project was that genetic estimates of effective population size, which can be made from empirical measurements of genetic drift, were equivalent to estimates of the successful number of spawners that is an important parameter in process of fisheries stock assessment. The broad objectives of this study were to 1. Critically evaluate a variety of mathematical methods of calculating effective spawner numbers (Ne) by a. conducting comprehensive computer simulations, and by b. analysis of empirical data collected from the Moreton Bay population of tiger prawns (P. esculentus). 2. Lay the groundwork for the application of the technology in the northern prawn fishery (NPF). 3. Produce software for the calculation of Ne, and to make it widely available. The project pulled together a range of mathematical models for estimating current effective population size from diverse sources. Some of them had been recently implemented with the latest statistical methods (eg. Bayesian framework Berthier, Beaumont et al. 2002), while others had lower profiles (eg. Pudovkin, Zaykin et al. 1996; Rousset and Raymond 1995). Computer code and later software with a user-friendly interface (NeEstimator) was produced to implement the methods. This was used as a basis for simulation experiments to evaluate the performance of the methods with an individual-based model of a prawn population. Following the guidelines suggested by computer simulations, the tiger prawn population in Moreton Bay (south-east Queensland) was sampled for genetic analysis with eight microsatellite loci in three successive spring spawning seasons in 2001, 2002 and 2003. As predicted by the simulations, the estimates had non-infinite upper confidence limits, which is a major achievement for the application of the method to a naturally-occurring, short generation, highly fecund invertebrate species. The genetic estimate of the number of successful spawners was around 1000 individuals in two consecutive years. This contrasts with about 500,000 prawns participating in spawning. It is not possible to distinguish successful from non-successful spawners so we suggest a high level of protection for the entire spawning population. We interpret the difference in numbers between successful and non-successful spawners as a large variation in the number of offspring per family that survive – a large number of families have no surviving offspring, while a few have a large number. We explored various ways in which Ne can be useful in fisheries management. It can be a surrogate for spawning population size, assuming the ratio between Ne and spawning population size has been previously calculated for that species. Alternatively, it can be a surrogate for recruitment, again assuming that the ratio between Ne and recruitment has been previously determined. The number of species that can be analysed in this way, however, is likely to be small because of species-specific life history requirements that need to be satisfied for accuracy. The most universal approach would be to integrate Ne with spawning stock-recruitment models, so that these models are more accurate when applied to fisheries populations. A pathway to achieve this was established in this project, which we predict will significantly improve fisheries sustainability in the future. Regardless of the success of integrating Ne into spawning stock-recruitment models, Ne could be used as a fisheries monitoring tool. Declines in spawning stock size or increases in natural or harvest mortality would be reflected by a decline in Ne. This would be good for data-poor fisheries and provides fishery independent information, however, we suggest a species-by-species approach. Some species may be too numerous or experiencing too much migration for the method to work. During the project two important theoretical studies of the simultaneous estimation of effective population size and migration were published (Vitalis and Couvet 2001b; Wang and Whitlock 2003). These methods, combined with collection of preliminary genetic data from the tiger prawn population in southern Gulf of Carpentaria population and a computer simulation study that evaluated the effect of differing reproductive strategies on genetic estimates, suggest that this technology could make an important contribution to the stock assessment process in the northern prawn fishery (NPF). Advances in the genomics world are rapid and already a cheaper, more reliable substitute for microsatellite loci in this technology is available. Digital data from single nucleotide polymorphisms (SNPs) are likely to super cede ‘analogue’ microsatellite data, making it cheaper and easier to apply the method to species with large population sizes.

Mathematical Model of Viral Kinetics In Vitro Estimates the Number of E2-CD81 Complexes Necessary for Hepatitis C Virus Entry

Interaction between the hepatitis C virus (HCV) envelope protein E2 and the host receptor CD81 is essential for HCV entry into target cells. The number of E2-CD81 complexes necessary for HCV entry has remained difficult to estimate experimentally. Using the recently developed cell culture systems that allow persistent HCV infection in vitro, the dependence of HCV entry and kinetics on CD81 expression has been measured. We reasoned that analysis of the latter experiments using a mathematical model of viral kinetics may yield estimates of the number of E2-CD81 complexes necessary for HCV entry. Here, we constructed a mathematical model of HCV viral kinetics in vitro, in which we accounted explicitly for the dependence of HCV entry on CD81 expression. Model predictions of viral kinetics are in quantitative agreement with experimental observations. Specifically, our model predicts triphasic viral kinetics in vitro, where the first phase is characterized by cell proliferation, the second by the infection of susceptible cells and the third by the growth of cells refractory to infection. By fitting model predictions to the above data, we were able to estimate the threshold number of E2-CD81 complexes necessary for HCV entry into human hepatoma-derived cells. We found that depending on the E2-CD81 binding affinity, between 1 and 13 E2-CD81 complexes are necessary for HCV entry. With this estimate, our model captured data from independent experiments that employed different HCV clones and cells with distinct CD81 expression levels, indicating that the estimate is robust. Our study thus quantifies the molecular requirements of HCV entry and suggests guidelines for intervention strategies that target the E2-CD81 interaction. Further, our model presents a framework for quantitative analyses of cell culture studies now extensively employed to investigate HCV infection.

Neural Representations for Sensory-Motor Control, III: Learning a Body-Centered Representation of 3-D Target Position

A neural model is described of how the brain may autonomously learn a body-centered representation of 3-D target position by combining information about retinal target position, eye position, and head position in real time. Such a body-centered spatial representation enables accurate movement commands to the limbs to be generated despite changes in the spatial relationships between the eyes, head, body, and limbs through time. The model learns a vector representation--otherwise known as a parcellated distributed representation--of target vergence with respect to the two eyes, and of the horizontal and vertical spherical angles of the target with respect to a cyclopean egocenter. Such a vergence-spherical representation has been reported in the caudal midbrain and medulla of the frog, as well as in psychophysical movement studies in humans. A head-centered vergence-spherical representation of foveated target position can be generated by two stages of opponent processing that combine corollary discharges of outflow movement signals to the two eyes. Sums and differences of opponent signals define angular and vergence coordinates, respectively. The head-centered representation interacts with a binocular visual representation of non-foveated target position to learn a visuomotor representation of both foveated and non-foveated target position that is capable of commanding yoked eye movementes. This head-centered vector representation also interacts with representations of neck movement commands to learn a body-centered estimate of target position that is capable of commanding coordinated arm movements. Learning occurs during head movements made while gaze remains fixed on a foveated target. An initial estimate is stored and a VOR-mediated gating signal prevents the stored estimate from being reset during a gaze-maintaining head movement. As the head moves, new estimates arc compared with the stored estimate to compute difference vectors which act as error signals that drive the learning process, as well as control the on-line merging of multimodal information.

Neonatal pain, parenting stress and interaction, in relation to cognitive and motor development at 8 and 18 months in preterm infants

Procedural pain in the neonatal intensive care unit triggers a cascade of physiological, behavioral and hormonal disruptions which may contribute to altered neurodevelopment in infants born very preterm, who undergo prolonged hospitalization at a time of physiological immaturity and rapid brain development. The aim of this study was to examine relationships between cumulative procedural pain (number of skin-breaking procedures from birth to term, adjusted for early illness severity and overall intravenous morphine exposure), and later cognitive, motor abilities and behavior in very preterm infants at 8 and 18 months corrected chronological age (CCA), and further, to evaluate the extent to which parenting factors modulate these relationships over time. Participants were N=211 infants (n=137 born preterm 32 weeks gestational age [GA] and n=74 full-term controls) followed prospectively since birth. Infants with significant neonatal brain injury (periventricular leucomalacia, grade 3 or 4 intraventricular hemorrhage) and/or major sensori-neural impairments, were excluded. Poorer cognition and motor function were associated with higher number of skin-breaking procedures, independent of early illness severity, overall intravenous morphine, and exposure to postnatal steroids. The number of skin-breaking procedures as a marker of neonatal pain was closely related to days on mechanical ventilation. In general, greater overall exposure to intravenous morphine was associated with poorer motor development at 8 months, but not at 18 months CCA, however, specific protocols for morphine administration were not evaluated. Lower parenting stress modulated effects of neonatal pain, only on cognitive outcome at 18 months.

Assessing the accuracy of CME speed and trajectory estimates from STEREO observations through a comparison of independent methods

We have estimated the speed and direction of propagation of a number of Coronal Mass Ejections (CMEs) using single-spacecraft data from the STEREO Heliospheric Imager (HI) wide-field cameras. In general, these values are in good agreement with those predicted by Thernisien, Vourlidas, and Howard in Solar Phys. 256, 111 -aEuro parts per thousand 130 (2009) using a forward modelling method to fit CMEs imaged by the STEREO COR2 coronagraphs. The directions of the CMEs predicted by both techniques are in good agreement despite the fact that many of the CMEs under study travel in directions that cause them to fade rapidly in the HI images. The velocities estimated from both techniques are in general agreement although there are some interesting differences that may provide evidence for the influence of the ambient solar wind on the speed of CMEs. The majority of CMEs with a velocity estimated to be below 400 km s(-1) in the COR2 field of view have higher estimated velocities in the HI field of view, while, conversely, those with COR2 velocities estimated to be above 400 km s(-1) have lower estimated HI velocities. We interpret this as evidence for the deceleration of fast CMEs and the acceleration of slower CMEs by interaction with the ambient solar wind beyond the COR2 field of view. We also show that the uncertainties in our derived parameters are influenced by the range of elongations over which each CME can be tracked. In order to reduce the uncertainty in the predicted arrival time of a CME at 1 Astronomical Unit (AU) to within six hours, the CME needs to be tracked out to at least 30 degrees elongation. This is in good agreement with predictions of the accuracy of our technique based on Monte Carlo simulations.

Automatic spiral analysis for objective assessment of motor symptoms in Parkinson's disease

Objective: To develop a method for objective quantification of PD motor symptoms related to Off episodes and peak dose dyskinesias, using spiral data gathered by using a touch screen telemetry device. The aim was to objectively characterize predominant motor phenotypes (bradykinesia and dyskinesia), to help in automating the process of visual interpretation of movement anomalies in spirals as rated by movement disorder specialists. Background: A retrospective analysis was conducted on recordings from 65 patients with advanced idiopathic PD from nine different clinics in Sweden, recruited from January 2006 until August 2010. In addition to the patient group, 10 healthy elderly subjects were recruited. Upper limb movement data were collected using a touch screen telemetry device from home environments of the subjects. Measurements with the device were performed four times per day during week-long test periods. On each test occasion, the subjects were asked to trace pre-drawn Archimedean spirals, using the dominant hand. The pre-drawn spiral was shown on the screen of the device. The spiral test was repeated three times per test occasion and they were instructed to complete it within 10 seconds. The device had a sampling rate of 10Hz and measured both position and time-stamps (in milliseconds) of the pen tip. Methods: Four independent raters (FB, DH, AJ and DN) used a web interface that animated the spiral drawings and allowed them to observe different kinematic features during the drawing process and to rate task performance. Initially, a number of kinematic features were assessed including ‘impairment’, ‘speed’, ‘irregularity’ and ‘hesitation’ followed by marking the predominant motor phenotype on a 3-category scale: tremor, bradykinesia and/or choreatic dyskinesia. There were only 2 test occasions for which all the four raters either classified them as tremor or could not identify the motor phenotype. Therefore, the two main motor phenotype categories were bradykinesia and dyskinesia. ‘Impairment’ was rated on a scale from 0 (no impairment) to 10 (extremely severe) whereas ‘speed’, ‘irregularity’ and ‘hesitation’ were rated on a scale from 0 (normal) to 4 (extremely severe). The proposed data-driven method consisted of the following steps. Initially, 28 spatiotemporal features were extracted from the time series signals before being presented to a Multilayer Perceptron (MLP) classifier. The features were based on different kinematic quantities of spirals including radius, angle, speed and velocity with the aim of measuring the severity of involuntary symptoms and discriminate between PD-specific (bradykinesia) and/or treatment-induced symptoms (dyskinesia). A Principal Component Analysis was applied on the features to reduce their dimensions where 4 relevant principal components (PCs) were retained and used as inputs to the MLP classifier. Finally, the MLP classifier mapped these components to the corresponding visually assessed motor phenotype scores for automating the process of scoring the bradykinesia and dyskinesia in PD patients whilst they draw spirals using the touch screen device. For motor phenotype (bradykinesia vs. dyskinesia) classification, the stratified 10-fold cross validation technique was employed. Results: There were good agreements between the four raters when rating the individual kinematic features with intra-class correlation coefficient (ICC) of 0.88 for ‘impairment’, 0.74 for ‘speed’, 0.70 for ‘irregularity’, and moderate agreements when rating ‘hesitation’ with an ICC of 0.49. When assessing the two main motor phenotype categories (bradykinesia or dyskinesia) in animated spirals the agreements between the four raters ranged from fair to moderate. There were good correlations between mean ratings of the four raters on individual kinematic features and computed scores. The MLP classifier classified the motor phenotype that is bradykinesia or dyskinesia with an accuracy of 85% in relation to visual classifications of the four movement disorder specialists. The test-retest reliability of the four PCs across the three spiral test trials was good with Cronbach’s Alpha coefficients of 0.80, 0.82, 0.54 and 0.49, respectively. These results indicate that the computed scores are stable and consistent over time. Significant differences were found between the two groups (patients and healthy elderly subjects) in all the PCs, except for the PC3. Conclusions: The proposed method automatically assessed the severity of unwanted symptoms and could reasonably well discriminate between PD-specific and/or treatment-induced motor symptoms, in relation to visual assessments of movement disorder specialists. The objective assessments could provide a time-effect summary score that could be useful for improving decision-making during symptom evaluation of individualized treatment when the goal is to maximize functional On time for patients while minimizing their Off episodes and troublesome dyskinesias.