NeuroMotion: Open-source platform with neuromechanical and deep network modules to generate surface EMG signals during voluntary movement.

Neuromechanical studies investigate how the nervous system interacts with the musculoskeletal (MSK) system to generate volitional movements. Such studies have been supported by simulation models that provide insights into variables that cannot be measured experimentally and allow a large number of conditions to be tested before the experimental analysis. However, current simulation models of electromyography (EMG), a core physiological signal in neuromechanical analyses, remain either limited in accuracy and conditions or are computationally heavy to apply. Here, we provide a computational platform to enable future work to overcome these limitations by presenting NeuroMotion, an open-source simulator that can modularly test a variety of approaches to the full-spectrum synthesis of EMG signals during voluntary movements. We demonstrate NeuroMotion using three sample modules. The first module is an upper-limb MSK model with OpenSim API to estimate the muscle fibre lengths and muscle activations during movements. The second module is BioMime, a deep neural network-based EMG generator that receives nonstationary physiological parameter inputs, like the afore-estimated muscle fibre lengths, and efficiently outputs motor unit action potentials (MUAPs). The third module is a motor unit pool model that transforms the muscle activations into discharge timings of motor units. The discharge timings are convolved with the output of BioMime to simulate EMG signals during the movement. We first show how MUAP waveforms change during different levels of physiological parameter variations and different movements. We then show that the synthetic EMG signals during two-degree-of-freedom hand and wrist movements can be used to augment experimental data for regressing joint angles. Ridge regressors trained on the synthetic dataset were directly used to predict joint angles from experimental data. In this way, NeuroMotion was able to generate full-spectrum EMG for the first use-case of human forearm electrophysiology during voluntary hand, wrist, and forearm movements. All intermediate variables are available, which allows the user to study cause-effect relationships in the complex neuromechanical system, fast iterate algorithms before collecting experimental data, and validate algorithms that estimate non-measurable parameters in experiments. We expect this modular platform will enable validation of generative EMG models, complement experimental approaches and empower neuromechanical research.

Deep Metric Learning With Locality Sensitive Mining for Self-Correcting Source Separation of Neural Spiking Signals.

Automated source separation algorithms have become a central tool in neuroengineering and neuroscience, where they are used to decompose neurophysiological signal into its constituent spiking sources. However, in noisy or highly multivariate recordings these decomposition techniques often make a large number of errors. Such mistakes degrade online human-machine interfacing methods and require costly post-hoc manual cleaning in the offline setting. In this article we propose an automated error correction methodology using a deep metric learning (DML) framework, generating embedding spaces in which spiking events can be both identified and assigned to their respective sources. Furthermore, we investigate the relative ability of different DML techniques to preserve the intraclass semantic structure needed to identify incorrect class labels in neurophysiological time series. Motivated by this analysis, we propose locality sensitive mining, an easily implemented sampling-based augmentation to typical DML losses which substantially improves the local semantic structure of the embedding space. We demonstrate the utility of this method to generate embedding spaces which can be used to automatically identify incorrectly labeled spiking events with high accuracy.

A Particle Swarm Optimised Independence Estimator for Blind Source Separation of Neurophysiological Time Series.

The decomposition of neurophysiological recordings into their constituent neural sources is of major importance to a diverse range of neuroscientific fields and neuroengineering applications. The advent of high density electrode probes and arrays has driven a major need for novel semi-automated and automated blind source separation methodologies that take advantage of the increased spatial resolution and coverage these new devices offer. Independent component analysis (ICA) offers a principled theoretical framework for such algorithms, but implementation inefficiencies often drive poor performance in practice, particularly for sparse sources. Here we observe that the use of a single non-linear optimization function to identify spiking sources with ICA often has a detrimental effect that precludes the recovery and correct separation of all spiking sources in the signal. We go on to propose a projection-pursuit ICA algorithm designed specifically for spiking sources, which uses a particle swarm methodology to adaptively traverse a polynomial family of non-linearities approximating the asymmetric cumulants of the sources. We robustly prove state-of-the-art decomposition performance on recordings from high density intramuscular probes and demonstrate how the particle swarm quickly finds optimal contrast non-linearities across a range of neurophysiological datasets.

Conditional Generative Models for Simulation of EMG During Naturalistic Movements.

Numerical models of electromyography (EMG) signals have provided a huge contribution to our fundamental understanding of human neurophysiology and remain a central pillar of motor neuroscience and the development of human-machine interfaces. However, while modern biophysical simulations based on finite element methods (FEMs) are highly accurate, they are extremely computationally expensive and thus are generally limited to modeling static systems such as isometrically contracting limbs. As a solution to this problem, we propose to use a conditional generative model to mimic the output of an advanced numerical model. To this end, we present BioMime, a conditional generative neural network trained adversarially to generate motor unit (MU) activation potential waveforms under a wide variety of volume conductor parameters. We demonstrate the ability of such a model to predictively interpolate between a much smaller number of numerical model's outputs with a high accuracy. Consequently, the computational load is dramatically reduced, which allows the rapid simulation of EMG signals during truly dynamic and naturalistic movements.

A myoelectric digital twin for fast and realistic modelling in deep learning.

Muscle electrophysiology has emerged as a powerful tool to drive human machine interfaces, with many new recent applications outside the traditional clinical domains, such as robotics and virtual reality. However, more sophisticated, functional, and robust decoding algorithms are required to meet the fine control requirements of these applications. Deep learning has shown high potential in meeting these demands, but requires a large amount of high-quality annotated data, which is expensive and time-consuming to acquire. Data augmentation using simulations, a strategy applied in other deep learning applications, has never been attempted in electromyography due to the absence of computationally efficient models. We introduce a concept of Myoelectric Digital Twin - highly realistic and fast computational model tailored for the training of deep learning algorithms. It enables simulation of arbitrary large and perfectly annotated datasets of realistic electromyography signals, allowing new approaches to muscular signal decoding, accelerating the development of human-machine interfaces.

Online Tracking of the Phase Difference Between Neural Drives to Antagonist Muscle Pairs in Essential Tremor Patients.

Transcutaneous electrical stimulation has been applied in tremor suppression applications. Out-of-phase stimulation strategies applied above or below motor threshold result in a significant attenuation of pathological tremor. For stimulation to be properly timed, the varying phase relationship between agonist-antagonist muscle activity during tremor needs to be accurately estimated in real-time. Here we propose an online tremor phase and frequency tracking technique for the customized control of electrical stimulation, based on a phase-locked loop (PLL) system applied to the estimated neural drive to muscles. Surface electromyography signals were recorded from the wrist extensor and flexor muscle groups of 13 essential tremor patients during postural tremor. The EMG signals were pre-processed and decomposed online and offline via the convolution kernel compensation algorithm to discriminate motor unit spike trains. The summation of motor unit spike trains detected for each muscle was bandpass filtered between 3 to 10 Hz to isolate the tremor related components of the neural drive to muscles. The estimated tremorogenic neural drive was used as input to a PLL that tracked the phase differences between the two muscle groups. The online estimated phase difference was compared with the phase calculated offline using a Hilbert Transform as a ground truth. The results showed a rate of agreement of 0.88 ± 0.22 between offline and online EMG decomposition. The PLL tracked the phase difference of tremor signals in real-time with an average correlation of 0.86 ± 0.16 with the ground truth (average error of 6.40° ± 3.49°). Finally, the online decomposition and phase estimation components were integrated with an electrical stimulator and applied in closed-loop on one patient, to representatively demonstrate the working principle of the full tremor suppression system. The results of this study support the feasibility of real-time estimation of the phase of tremorogenic neural drive to muscles, providing a methodology for future tremor-suppression neuroprostheses.

Deep Learning for Robust Decomposition of High-Density Surface EMG Signals.

Blind source separation (BSS) algorithms, such as gradient convolution kernel compensation (gCKC), can efficiently and accurately decompose high-density surface electromyography (HD-sEMG) signals into constituent motor unit (MU) action potential trains. Once the separation matrix is blindly estimated on a signal interval, it is also possible to apply the same matrix to subsequent signal segments. Nonetheless, the trained separation matrices are sub-optimal in noisy conditions and require that incoming data undergo computationally expensive whitening. One unexplored alternative is to instead use the paired HD-sEMG signal and BSS output to train a model to predict MU activations within a supervised learning framework. A gated recurrent unit (GRU) network was trained to decompose both simulated and experimental unwhitened HD-sEMG signal using the output of the gCKC algorithm. The results on the experimental data were validated by comparison with the decomposition of concurrently recorded intramuscular EMG signals. The GRU network outperformed gCKC at low signal-to-noise ratios, proving superior performance in generalising to new data. Using 12 seconds of experimental data per recording, the GRU performed similarly to gCKC, at rates of agreement of 92.5% (84.5%-97.5%) and 94.9% (88.8%-100.0%) respectively for GRU and gCKC against matched intramuscular sources.

Analysis of Retinal Segmentation Changes at High Altitude With and Without Acetazolamide.

Purpose: Our aim was to assess retinal venous diameter and segmented retinal layer thickness variation in acute systemic hypoxia with and without acetazolamide and to relate these changes to high altitude headache (HAH), as a proxy for intracerebral pathophysiology. Methods: A total of 20 subjects participated in a 4-day ascent to the Margherita Hut (4,559 m) on Monte Rosa in the Italian Alps. Each participant was randomized to either oral acetazolamide 250 mg twice daily or placebo. A combination of digital imaging and optical coherence tomography was used to measure retinal vessel diameter and retinal layer thickness. Clinically-assessed HAH was recorded. Results: A total of 18 participants had usable digital and OCT images, with 12 developing HAH. Significant thickening was seen only in the two inner layers of the retina, the retinal nerve fiber layer (RNFL) and ganglion cell layer (GCL) at P = 0.012 and P = 0.010, respectively, independent of acetazolamide. There was a significant positive correlation between HAH and both retinal venous diameter (T = 4.953, P = 0.001) and retinal artery diameter (T = 2.865, P = 0.015), with both unaffected by acetazolamide (F = 0.439, P = 0.518). Conclusions: Retinal venous diameter correlates positively with HAH, adding further evidence for the proposed venous outflow limitation mechanism. The inner layers of the retina swelled disproportionately when compared to the outer layers under conditions of systemic hypoxia. Acetazolamide does not appear to influence altitudinal changes of retinal layers and vasculature.

Peer-led live research demonstrations: challenging medical student misconceptions about research.

Modern health care provision is now fundamentally evidence based, meaning competency in academic medicine is integral to medical training. The Integrated Academic Training pathway provides focussed training in this area at a postgraduate level but no such provision exists at an undergraduate level. A number of peer-led academic societies have emerged across the UK to provide education and support for undergraduates but there is little evidence about the type of peer-led interventions that are effective. We report here the findings of one such peer-led organization, the Warwick Academic Medicine Society. We found that traditional educational interventions, including didactic lectures and small-group teaching, are effective at inspiring students regarding academic medicine but poor at translating this enthusiasm into sustained involvement in research. We find this disparity to be centred on misconceptions amongst students regarding the time and skills required to meaningfully contribute to a research project. Further, we introduce the concept of the Live Research Demonstration (LRD), a novel peer-led educational intervention which aims to address these misconceptions and improve involvement of students in research. Initial pilots of the LRD concept have shown significant promise and we recommend a larger trial across multiple localities to confirm its educational benefits.