Search

The central nervous system (CNS) may produce coordinated motor outputs via the combination of motor modules representable as muscle synergies. Identification of muscle synergies has hitherto relied on applying factorization algorithms to multimuscle electromyographic data (EMGs) recorded during motor behaviors. Recent studies have attempted to validate the neural basis of the muscle synergies identified by independently retrieving the muscle synergies through CNS manipulations and analytic techniques such as spike-triggered averaging of EMGs. Experimental data have demonstrated the pivotal role of the spinal premotor interneurons in the synergies’ organization and the presence of motor cortical loci whose stimulations offer access to the synergies, but whether the motor cortex is also involved in organizing the synergies has remained unsettled. We argue that one difficulty inherent in current approaches to probing the synergies’ neural basis is that the EMG generative model based on linear combination of synergies and the decomposition algorithms used for synergy identification are not grounded on enough prior knowledge from neurophysiology. Progress may be facilitated by constraining or updating the model and algorithms with knowledge derived directly from CNS manipulations or recordings. An investigative framework based on evaluating the relevance of neurophysiologically constrained models of muscle synergies to natural motor behaviors will allow a more sophisticated understanding of motor modularity, which will help the community move forward from the current debate on the neural versus nonneural origin of muscle synergies.

Recent studies have investigated muscle synergies as biomarkers for stroke, but it remains controversial if muscle synergies and clinical observation convey the same information on motor impairment. We aim to identify whether muscle synergies and clinical scales convey the same information or not. Post-stroke patients were administered an upper limb treatment. Before (T0) and after (T1) treatment, we assessed motor performance with clinical scales and motor output with EMG-derived muscle synergies. We implemented an exploratory factor analysis (EFA) and a confirmatory factor analysis (CFA) to identify the underlying relationships among all variables, at T0 and T1, and a general linear regression model to infer any relationships between the similarity between the affected and unaffected synergies (Median-sp) and clinical outcomes at T0. Clinical variables improved with rehabilitation whereas muscle-synergy parameters did not show any significant change. EFA and CFA showed that clinical variables and muscle-synergy parameters (except Median-sp) were grouped into different factors. Regression model showed that Median-sp could be well predicted by clinical scales. The information underlying clinical scales and muscle synergies are therefore different. However, clinical scales well predicted the similarity between the affected and unaffected synergies. Our results may have implications on personalizing rehabilitation protocols.

The Electromyography-based Pattern-Recognition (EMG-PR) framework has been investigated for almost three decades towards developing an intuitive myoelectric prosthesis. To utilize the knowledge of the underlying neurophysiological processes of natural movements, the concept of muscle synergy has been applied in prosthesis control and proved to be of great potential recently. For a muscle-synergy-based myoelectric system, the variation of muscle contraction force is also a confounding factor. This study evaluates the robustness of muscle synergies under a variant force level for forearm movements. Six channels of forearm surface EMG were recorded from six healthy subjects when they performed 4 movements (hand open, hand close, wrist flexion, and wrist extension) using low, moderate, and high force, respectively. Muscle synergies were extracted from the EMG using the alternating nonnegativity constrained least squares and active set (NNLS) algorithm. Three analytic strategies were adopted to examine whether muscle synergies were conserved under different force levels. Our results consistently showed that there exists fixed, robust muscle synergies across force levels. This outcome would provide valuable insights to the implementation of muscle- synergy-based assistive technology for the upper extremity.

Chronic stroke survivors often suffer from gait impairment resistant to intervention. Recent rehabilitation strategies based on gait training with powered exoskeletons appear promising, but whether chronic survivors may benefit from them remains controversial. We evaluated the potential of exoskeletal gait training in restoring normal motor outputs in chronic survivors (N = 10) by recording electromyographic signals (EMGs, 28 muscles both legs) as they adapted to exoskeletal perturbations, and examined whether any EMG alterations after adaptation were underpinned by closer-to-normal muscle synergies. A unilateral ankle-foot orthosis that produced dorsiflexor torque on the paretic leg during swing was tested. Over a single session, subjects walked overground without exoskeleton (FREE), then with the unpowered exoskeleton (OFF), and finally with the powered exoskeleton (ON). Muscle synergies were identified from EMGs using non-negative matrix factorization. During adaptation to OFF, some paretic-side synergies became more dissimilar to their nonparetic-side counterparts. During adaptation to ON, in half of the subjects some paretic-side synergies became closer to their nonparetic references relative to their similarity at FREE as these paretic-side synergies became sparser in muscle components. Across subjects, level of inter-side similarity increase correlated negatively with the degree of gait temporal asymmetry at FREE. Our results demonstrate the possibility that for some survivors, exoskeletal training may promote closer-to-normal muscle synergies. But to fully achieve this, the active force must trigger adaptive processes that offset any undesired synergy changes arising from adaptation to the device’s mechanical properties while also fostering the reemergence of the normal synergies.

Complex motor commands for human locomotion are generated through the combination of motor modules representable as muscle synergies. Recent data have argued that muscle synergies are inborn or determined early in life, but development of the neuro-musculoskeletal system and acquisition of new skills may demand fine-tuning or reshaping of the early synergies. We seek to understand how locomotor synergies change during development and training by studying the synergies for running in preschoolers and diverse adults from sedentary subjects to elite marathoners, totaling 63 subjects assessed over 100 sessions. During development, synergies are fractionated into units with fewer muscles. As adults train to run, specific synergies coalesce to become merged synergies. Presences of specific synergy-merging patterns correlate with enhanced or reduced running efficiency. Fractionation and merging of muscle synergies may be a mechanism for modifying early motor modules (Nature) to accommodate the changing limb biomechanics and influences from sensorimotor training (Nurture)., Motor commands for human locomotion are generated by combination of muscle synergies. In humans, muscle synergies for running exhibit considerable plasticity during child-to-adult development and adult training to meet the constantly changing biomechanical and efficiency demands.

Identification of runner’s performance level is critical to coaching, performance enhancement and injury prevention. Machine learning techniques have been developed to measure biomechanical parameters with body-worn inertial measurement unit (IMU) sensors. However, a robust method to classify runners is still unavailable. In this paper, we developed two models to classify running performance and predict biomechanical parameters of 30 subjects. We named the models RunNet-CNN and RunNet-MLP based on their architectures: convolutional neural network (CNN) and multilayer perceptron (MLP), respectively. In addition, we examined two validation approaches, subject-wise (leave-one-subject-out) and record-wise. RunNet-MLP classified runner’s performance levels with an overall accuracy of 97.1%. Our results also showed that RunNet-CNN outperformed RunNet-MLP and gradient boosting decision tree in predicting biomechanical parameters. RunNet-CNN showed good agreement (R2 > 0.9) with the ground-truth reference on biomechanical parameters. The prediction accuracy for the record-wise method was better than the subject-wise method regardless of biomechanical parameters or models. Our findings showed the viability of using IMUs to produce reliable prediction of runners’ performance levels and biomechanical parameters.

The past decades have witnessed remarkable progress in neural technologies such as functional electrical stimulation (FES) and their applications in neurorehabilitation and neuromodulation. These advances are powered by new neuroscientific understandings of the organization and compositionality of neuromuscular control illuminating how muscle groups may be activated together as discrete units known as muscle synergies. These parallel developments have promoted novel approaches to clinical rehabilitation for neurological disorders that are insurmountable to current treatments. One such breakthrough is the evolution of FES as a therapeutic tool in poststroke rehabilitation with an interventional strategy particularly inspired by the concept that muscles in humans may be purposefully coordinated through neuromotor modules represented as muscle synergies. This paper will review recent advances in multichannel FES technology, its potential applications in poststroke rehabilitation, new findings that support the neurological basis of the muscle synergies for generating natural motor tasks, and the application of the muscle-synergy concept in poststroke assessment and rehabilitation of motor impairment. Finally, we will recommend future directions of development in relation to assistive FES and synergy-driven adaptive training for poststroke rehabilitation.

Humans respond to mechanical perturbations that affect their gait by changing their motor control strategy. Previous work indicates that adaptation during gait is context dependent, and perturbations altering long-term stability are compensated for even at the cost of higher energy expenditure. However, it is unclear if gait adaptation is driven by unilateral or bilateral mechanisms, and what the roles of feedback and feedforward control are in the generation of compensatory responses. Here, we used a robot-based adaptation paradigm to investigate if feedback/feedforward and unilateral/bilateral contributions to locomotor adaptation are also context dependent in healthy adults. A robot was used to induce two opposite unilateral mechanical perturbations affecting the step length over multiple gait cycles. Electromyographic signals were collected and analyzed to determine how muscle synergies change in response to perturbations. The results unraveled different unilateral modulation dynamics of the muscle-synergy activations during adaptation, characterized by the combination of a slow-progressive feedforward process and a fast-reactive feedback-driven process. The relative unilateral contributions of the two processes to motor-output adjustments, however, depended on which perturbation was delivered. Overall, these observations provide evidence that, in humans, both descending and afferent drives project onto the same spinal interneuronal networks that encode locomotor muscle synergies. Enterprise Ireland European Commission Horizon 2020 Science Foundation Ireland Insight Research Centre German Federal Ministry of Research and Education Bavarian Ministry of Economic Affairs European Space Agency German Federal Ministry of Economic Affairs and Energy National Institutes of Health (NIH) Lemelson Foundation Medtronic Hong Kong Research Grants Council Chinese University of Hong Kong Guangdong Provincial Hospital of Traditional Chinese Medicine Golden King Shun Hing Institute of Advanced Engineering American Heart Association Department of Defense Michael J Fox Foundation National Science Foundation Peabody Foundation Barrett Technology BioSensics Veristride Emerge Diagnostics MC10 Mitsui Chemicals Shimmer Research SynPhNe

Objective: Motor variability – performance variations across task repetitions – has been assumed to be undesirable. But recent studies argue that variability facilitates early motor learning by allowing exploratory search of reward-generating motion, and that variability's structure may be modulated by neural circuits for furthering learning. What are the neural sources of learning-relevant motor variability and its modulation in humans of different ages? Methods: Elderlies and young adults played a 3-session virtual bowling while multi-muscle electromyographic signals were collected. We quantified trial-to-trial variability of muscle synergies – neuromotor control modules – and of their activations. Results: In elderlies, bowling-score gain correlated with change of activation timing variability of specific synergies, but in young adults, with variability changes of synergy-activation magnitude, and of the synergies themselves. Conclusions: Variability modulation of specific muscle synergies and their activations contribute to early motor learning. Elderly and young individuals may rely on different aspects of motor variability to drive learning.

Background Neuromuscular electrical stimulation (NMES) is extensively used in stroke motor rehabilitation. How it promotes motor recovery remains only partially understood. NMES could change muscular properties, produce altered sensory inputs, and modulate fluctuations of cortical activities; but the potential contribution from cortico-muscular couplings during NMES synchronized with dynamic movement has rarely been discussed. Method We investigated cortico-muscular interactions during passive, active, and NMES rhythmic pedaling in healthy subjects and chronic stroke survivors. EEG (128 channels), EMG (4 unilateral lower limb muscles) and movement parameters were measured during 3 sessions of constant-speed pedaling. Sensory-level NMES (20 mA) was applied to the muscles, and cyclic stimulation patterns were synchronized with the EMG during pedaling cycles. Adaptive mixture independent component analysis was utilized to determine the movement-related electro-cortical sources and the source dipole clusters. A directed cortico-muscular coupling analysis was conducted between representative source clusters and the EMGs using generalized partial directed coherence (GPDC). The bidirectional GPDC was compared across muscles and pedaling sessions for post-stroke and healthy subjects. Results Directed cortico-muscular coupling of NMES cycling was more similar to that of active pedaling than to that of passive pedaling for the tested muscles. For healthy subjects, sensory-level NMES could modulate GPDC of both ascending and descending pathways. Whereas for stroke survivors, NMES could modulate GPDC of only the ascending pathways. Conclusions By clarifying how NMES influences neuromuscular control during pedaling in healthy and post-stroke subjects, our results indicate the potential limitation of sensory-level NMES in promoting sensorimotor recovery in chronic stroke survivors.

Background: How to promote true recovery from poststroke upper limb motor impairment has remained an urgent public health problem. Acupuncture has the potential to facilitate poststroke recovery. Abdominal acupuncture, based on the recently discovered acupoint system on the anterior abdominal wall, appears attractive because it induces less pain, and allows concurrent limb rehabilitative training during treatment. However, its clinical efficacy has not been systematically demonstrated, and its neurophysiological mechanism has remained obscure. Methods: First-onset stroke survivors (0.5-3 months post-stroke) will be randomly divided into 3 groups (N=22 in each), respectively receiving (1) abdominal acupuncture, (2) abdominal acupuncture with sham needles, and (3) no acupuncture. All subjects will concurrently receive basic treatment, including upper limb rehabilitative training and measures for secondary stroke prevention. Clinical scores reflective of motor functions and impairment (Wolf motor function test, Fugl-Meyer assessment, Brunnstrom staging), evaluation of daily life ability, surface electromyography, and motor-imagery functional magnetic resonance imaging will be collected as outcome measures before and after intervention. Upper-limb muscle synergies will be identified from the collected surface electromyography. Discussion: The study will use abdominal acupuncture to improve recovery from motor dysfunction of the upper limb after stroke, to observe the effects of abdominal acupuncture on post-stroke upper limb motor functions, and to analyze the relationship between changes in upper-limb functions and measurements from both multi-muscle surface electromyographic data and brain activations during motor imagery from functional magnetic resonance imaging, so as to explore possible mechanisms of neuroplasticity associated with abdominal acupuncture. Trial registration: This trial was registered with the ClinicalTrials.gov (ID: NCT03712085) on 7th July 2018, and last updated on 16th Oct 2018.

A unified approach to nonnegative matrix factorization based on the theory of generalized linear models is proposed. This approach embeds a variety of statistical models, including the exponential family, within a single theoretical framework and provides a unified view of such factorizations from the perspective of quasi-likelihood. Using this framework, a family of algorithms for handling signal-dependent noise is developed and its convergence proved using the expectation-maximization algorithm. In addition, a measure to evaluate the goodness of fit of the resulting factorization is described. The proposed methods allow modeling of nonlinear effects using appropriate link functions and are illustrated using an application in biomedical signal processing.

OBJECTIVE To promote clinical applications of muscle-synergy-based neurorehabilitation techniques, this study aims to clarify any potential modulations of both the muscular compositions and temporal activations of forearm muscle synergies for multiple movements under variant force levels and arm positions. APPROACH Two groups of healthy subjects participated in this study. Electromyography (EMG) signals were collected when they performed four hand and wrist movements under variant constraints-three different force levels for one group and five arm positions for the other. Muscle synergies were extracted from the EMGs, and their robustness across variant force levels and arm positions was separately assessed by evaluating their across-condition structure similarity, cross-validation, and cluster analysis. The synergies' activation coefficients across the variant constraints were also compared, and the coefficients were used to discriminate the different force levels and the arm positions, respectively. MAIN RESULTS Overall, the muscle synergies were relatively fixed across variant constraints, but they were more robust to variant forces than to changing arm positions. The activations of muscle synergies depended largely on the level of contraction force and could discriminate the force levels very well, but the coefficients corresponding to different arm positions discriminated the positions with lower accuracy. Similar results were found for all types of forearm movement analyzed. SIGNIFICANCE With our experiment and subject-specific analysis, only slight modulation of the muscular compositions of forearm muscle synergies was found under variant force and arm position constraints. Our results may shed valuable insights to future design of both muscle-synergy-based assistive robots and motor-function assessments.

Neural Plasticity is an interdisciplinary journal dedicated to the publication of articles related to all aspects of neural plasticity, with special emphasis on its functional significance as reflected in behavior and in psychopathology.

Constant-force isometric muscle training is useful for increasing the maximal strength , rehabilitation and work-fatigue assessment. Earlier studies have shown that muscle fatigue characteristics can be used for evaluating muscle endurance limit.To predict muscle endurance time during isometric task using frequency spectrum characteristics of surface electromyography signals along with analysis of frequency spectrum shape and scale during fatigue accumulation.Thirteen subjects performed isometric lateral raise at 60% MVC of deltoid (lateral) till endurance limit. Time windowed sEMG frequency spectrum was modelled using 2-parameter distributions namely Gamma and Weibull for spectrum analysis and endurance prediction.Gamma distribution provided better spectrum fitting (P 0.001) than Weibull distribution. Spectrum Distribution demonstrated no change in shape but shifted towards lower frequency with increase of magnitude at characteristic mode frequency. Support Vector Regression based algorithm was developed for endurance time estimation using features derived from fitted frequency spectrum. Time taken till endurance limit for acquired dataset 38.53 ± 17.33 s (Mean ± Standard Deviation) was predicted with error of 0.029 ± 4.19 s . R-square: 0.956, training and test sets RMSE was calculated as 3.96 and 4.29 s respectively. The application of the algorithm suggested that model required 70% of sEMG signal from maximum time of endurance for high prediction accuracy.Endurance Limit prediction algorithm was developed for quantification of endurance time for optimizing isometric training and rehabilitation. Our method could help personalize and change conventional training method of same weight and duration for all subjects with optimized training parameters, based upon individual sEMG activity.

Motor modules are neural entities hypothesized to be building blocks of movement construction. How motor modules are underpinned by neural circuits has remained obscured. As a first step towards dissecting these circuits, we optogenetically evoked motor outputs from the lumbosacral spinal cord of two strains of transgenic mice – the Chat, with channelrhodopsin (ChR2) expressed in motoneurons, and the Thy1, expressed in putatively excitatory neurons. Motor output was represented as a spatial field of isometric ankle force. We found that Thy1 force fields were more complex and diverse in structure than Chat fields: the Thy1 fields comprised mostly non-parallel vectors while the Chat fields, mostly parallel vectors. In both, most fields elicited by co-stimulation of two laser beams were well explained by linear combination of the separately-evoked fields. We interpreted the Thy1 force fields as representations of spinal motor modules. Our comparison of the Chat and Thy1 fields allowed us to conclude, with reasonable certainty, that the structure of neuromotor modules originates from excitatory spinal interneurons. Our results not only demonstrate, for the first time using optogenetics, how the spinal modules follow linearity in their combinations, but also provide a reference against which future optogenetic studies of modularity can be compared.

There exists a gap from experimental data to the understanding of neural control of movements. This research topic was dedicated to promote computational modeling approach that can facilitate data interpretation (Niu et al.; Ranjbaran and Galiana; Pearson et al.; Sharif Razavian et al.; Malik et al.), elucidate control theories (Ueyama; Ota et al.; Takemura et al.), shed light on systemic mechanisms (Buhrmann and DiPaolo; Pearson et al.; Li et al.), suggest testable hypothesis (Loeb and Tsianos; Jiang et al.), and aid design of rehabilitation or therapeutic strategies (Zitella et al.). The 14 articles reflected these different aspects of computational modeling in bridging this gap between functions of neural circuits and observable behaviors. This research topic demonstrated that computational modeling is playing a more and more prominent role in sensorimotor control studies.

Previous studies using intact and spinalized animals have suggested that coordinated movements can be generated by appropriate combinations of muscle synergies controlled by the central nervous system (CNS). However, which CNS regions are responsible for expressing muscle synergies remains an open question. We address whether the brain stem and spinal cord are involved in expressing muscle synergies used for executing a range of natural movements. We analyzed the electromyographic (EMG) data recorded from frog leg muscles before and after transection at different levels of the neuraxis—rostral midbrain (brain stem preparations), rostral medulla (medullary preparations), and the spinal-medullary junction (spinal preparations). Brain stem frogs could jump, swim, kick, and step, while medullary frogs could perform only a partial repertoire of movements. In spinal frogs, cutaneous reflexes could be elicited. Systematic EMG analysis found two different synergy types: 1) synergies shared between pre- and posttransection states and 2) synergies specific to individual states. Almost all synergies found in natural movements persisted after transection at rostral midbrain or medulla but not at the spinal-medullary junction for swim and step. Some pretransection- and posttransection-specific synergies for a certain behavior appeared as shared synergies for other motor behaviors of the same animal. These results suggest that the medulla and spinal cord are sufficient for the expression of most muscle synergies in frog behaviors. Overall, this study provides further evidence supporting the idea that motor behaviors may be constructed by muscle synergies organized within the brain stem and spinal cord and activated by descending commands from supraspinal areas.

It has been suggested that the motor system may circumvent the difficulty of controlling many degrees of freedom in the musculoskeletal apparatus by generating motor outputs through a combination of discrete muscle synergies. How a discretely organized motor system compensates for diverse perturbations has remained elusive. Here, we investigate whether motor responses observed after an inertial-load perturbation can be generated by altering the recruitment of synergies normally used for constructing unperturbed movements. Electromyographic (EMG, 13 muscles) data were collected from the bullfrog hindlimb during natural behaviors before, during, and after the same limb was loaded by a weight attached to the calf. Kinematic analysis reveals the absence of aftereffect on load removal, suggesting that load-related EMG changes were results of immediate motor pattern adjustments. We then extracted synergies from EMGs using the nonnegative matrix factorization algorithm and developed a procedure for assessing the extent of synergy sharing across different loading conditions. Most synergies extracted were found to be activated in all loaded and unloaded conditions. However, for certain synergies, the amplitude, duration, and/or onset time of their activation bursts were up- or down-modulated during loading. Behavioral parameterizations reveal that load-related modulation of synergy activations depended on the behavioral variety (e.g., kick direction and amplitude) and the movement phase performed. Our results suggest that muscle synergies are robust across different dynamic conditions and immediate motor adjustments can be accomplished by modulating synergy activations. An appendix describes the novel procedure we developed, useful for discovering shared and specific features from multiple data sets.

We review experiments supporting the hypothesis that the vertebrate motor system produces movements by combining a small number of units of motor output. Using a variety of approaches such as microstimulation of the spinal cord, NMDA iontophoresis, and an examination of natural behaviors in intact and deafferented animals we have provided evidence for a modular organization of the spinal cord. A module is a functional unit in the spinal cord that generates a specific motor output by imposing a specific pattern of muscle activation. Such an organization might help to simplify the production of movements by reducing the degrees of freedom that need to be specified.

The non-negative matrix factorization algorithm (NMF) decomposes a data matrix into a set of non-negative basis vectors, each scaled by a coefficient. In its original formulation, the NMF assumes the data samples and dimensions to be independently distributed, making it a less-than-ideal algorithm for the analysis of time series data with temporal correlations. Here, we seek to derive an NMF that accounts for temporal dependencies in the data by explicitly incorporating a very simple temporal constraint for the coefficients into the NMF update rules. We applied the modified algorithm to 2 multi-dimensional electromyographic data sets collected from the human upper-limb to identify muscle synergies. We found that because it reduced the number of free parameters in the model, our modified NMF made it possible to use the Akaike Information Criterion to objectively identify a model order (i.e., the number of muscle synergies composing the data) that is more functionally interpretable, and closer to the numbers previously determined using ad hoc measures.

Several recent studies have used matrix factorization algorithms to assess the hypothesis that behaviors might be produced through the combination of a small number of muscle synergies. Although generally agreeing in their basic conclusions, these studies have used a range of different algorithms, making their interpretation and integration difficult. We therefore compared the performance of these different algorithms on both simulated and experimental data sets. We focused on the ability of these algorithms to identify the set of synergies underlying a data set. All data sets consisted of nonnegative values, reflecting the nonnegative data of muscle activation patterns. We found that the performance of principal component analysis (PCA) was generally lower than that of the other algorithms in identifying muscle synergies. Factor analysis (FA) with varimax rotation was better than PCA, and was generally at the same levels as independent component analysis (ICA) and nonnegative matrix factorization (NMF). ICA performed very well on data sets corrupted by constant variance Gaussian noise, but was impaired on data sets with signal-dependent noise and when synergy activation coefficients were correlated. Nonnegative matrix factorization (NMF) performed similarly to ICA and FA on data sets with signal-dependent noise and was generally robust across data sets. The best algorithms were ICA applied to the subspace defined by PCA (ICAPCA) and a version of probabilistic ICA with nonnegativity constraints (pICA). We also evaluated some commonly used criteria to identify the number of synergies underlying a data set, finding that only likelihood ratios based on factor analysis identified the correct number of synergies for data sets with signal-dependent noise in some cases. We then proposed an ad hoc procedure, finding that it was able to identify the correct number in a larger number of cases. Finally, we applied these methods to an experimentally obtained data set. The best performing algorithms (FA, ICA, NMF, ICAPCA, pICA) identified synergies very similar to one another. Based on these results, we discuss guidelines for using factorization algorithms to analyze muscle activation patterns. More generally, the ability of several algorithms to identify the correct muscle synergies and activation coefficients in simulated data, combined with their consistency when applied to physiological data sets, suggests that the muscle synergies found by a particular algorithm are not an artifact of that algorithm, but reflect basic aspects of the organization of muscle activation patterns underlying behaviors.

Nonnegative matrix factorization (NMF) by the multiplicative updates algorithm is a powerful machine learning method for decomposing a high-dimensional nonnegative matrix V into two nonnegative matrices, W and H, where [Formula: see text]. It has been successfully applied in the analysis and interpretation of large-scale data arising in neuroscience, computational biology, and natural language processing, among other areas. A distinctive feature of NMF is its nonnegativity constraints that allow only additive linear combinations of the data, thus enabling it to learn parts that have distinct physical representations in reality. In this letter, we describe an information-theoretic approach to NMF for signal-dependent noise based on the generalized inverse gaussian model. Specifically, we propose three novel algorithms in this setting, each based on multiplicative updates, and prove monotonicity of updates using the EM algorithm. In addition, we develop algorithm-specific measures to evaluate their goodness of fit on data. Our methods are demonstrated using experimental data from electromyography studies, as well as simulated data in the extraction of muscle synergies, and compared with existing algorithms for signal-dependent noise.

Support against gravity is a promising strategy to enable upper extremity rehabilitation in neurological patients. While such support facilitates reaching in subjects with severe impairments, its effect on the reorganization of muscle activations has not been investigated. In particular, it remains unclear what amount of weight compensation should be employed to achieve maximal functional recovery. To provide insights into these issues, we studied muscle synergies during the execution of 3D reaching movements under six different levels of arm-weight support. Our preliminary results suggest that gravity compensation decreases the level of muscle synergy activations required for task performance. However, whether weight support could induce changes in the muscle-synergy structures awaits further confirmation.

The primary motor cortex (M1) supports motor skill learning, yet little is known about the genes that contribute to motor cortical plasticity. Such knowledge could identify candidate molecules whose targeting might enable a new understanding of motor cortical functions, and provide new drug targets for the treatment of diseases which impair motor function, such as ischemic stroke. Here, we assess changes in the motor-cortical transcriptome across different stages of motor skill acquisition. Adult rats were trained on a gradually acquired appetitive reach and grasp task that required different strategies for successful pellet retrieval, or a sham version of the task in which the rats received pellet reward without needing to develop the reach and grasp skill. Tissue was harvested from the forelimb motor-cortical area either before training commenced, prior to the initial rise in task performance, or at peak performance. Differential classes of gene expression were observed at the time point immediately preceding motor task improvement. Functional clustering revealed that gene expression changes were related to the synapse, development, intracellular signaling, and the fibroblast growth factor (FGF) family, with many modulated genes known to regulate synaptic plasticity, synaptogenesis, and cytoskeletal dynamics. The modulated expression of synaptic genes likely reflects ongoing network reorganization from commencement of training till the point of task improvement, suggesting that motor performance improves only after sufficient modifications in the cortical circuitry have accumulated. The regulated FGF-related genes may together contribute to M1 remodeling through their roles in synaptic growth and maturation., McGovern Institute for Brain Research at MIT, National Institutes of Health (U.S.) ((NIH grant 1-RC1-NS068103-01), National Institutes of Health (U.S.) (NIH grant R01-MH084966), Roberto Rocca Education Program (Fellowship), Massachusetts Institute of Technology. Undergraduate Research Opportunities Program (Fellowship), Italy. Ministero dell'istruzione, dell'università e della ricerca (MIUR grant RBIN04H5AS), Italy. Ministero dell'istruzione, dell'università e della ricerca (MIUR grant RBLA03FLJC), Italy. Ministero dell'istruzione, dell'università e della ricerca (FIRB n. RBAP10L8TY)

The Brain and Cognitive Sciences (BCS) Department of Massachusetts Institute of Technology (MIT) is a place where science and engineering intersect. Housed under the roof of the recently built BCS Complex (MIT Building 46) (Figures 1 and 2), the largest neuroscience facility in the world, BCS investigators work closely with engineers both to elucidate neural mechanisms behind our cognitive activities and to apply their neuroscientific knowledge to improve treatments for diseases of the brain. Here, we discuss Edward Boyden and Emery Brown, two BCS faculty members who employ genetic engineering techniques and signalprocessing methods, respectively, to further enhance their understanding of the mysterious brain.

Production of voluntary movements relies critically on the functional integration of several motor cortical areas, such as the primary motor cortex, and the spinal circuitries. Surprisingly, after almost 40 years of research, how the motor cortices specify descending neural signals destined for the downstream interneurons and motoneurons has remained elusive. In light of the many recent experimental demonstrations that the motor system may coordinate muscle activations through a linear combination of muscle synergies, we hypothesize that the motor cortices may function to select and activate fixed muscle synergies specified by the spinal or brainstem networks. To test this hypothesis, we recorded electromyograms (EMGs) from 12–16 upper arm and shoulder muscles from both the unaffected and the stroke-affected arms of stroke patients having moderate-to-severe unilateral ischemic lesions in the frontal motor cortical areas. Analyses of EMGs using a nonnegative matrix factorization algorithm revealed that in seven of eight patients the muscular compositions of the synergies for both the unaffected and the affected arms were strikingly similar to each other despite differences in motor performance between the arms, and differences in cerebral lesion sizes and locations between patients. This robustness of muscle synergies that we observed supports the notion that descending cortical signals represent neuronal drives that select, activate, and flexibly combine muscle synergies specified by networks in the spinal cord and/or brainstem. Our conclusion also suggests an approach to stroke rehabilitation by focusing on those synergies with altered activations after stroke.

Previous studies have suggested that the motor system may simplify control by combining a small number of muscle synergies represented as activation profiles across a set of muscles. The role of sensory feedback in the activation and organization of synergies has remained an open question. Here, we assess to what extent the motor system relies on centrally organized synergies activated by spinal and/or supraspinal commands to generate motor outputs by analyzing electromyographic (EMG) signals collected from 13 hindlimb muscles of the bullfrog during swimming and jumping, before and after deafferentation. We first established that, for both behaviors, the intact and deafferented data sets possess low and similar dimensionalities. Subsequently, we used a novel reformulation of the non-negative matrix factorization algorithm to simultaneously search for synergies shared by, and synergies specific to, the intact and deafferented data sets. Most muscle synergies were identified as shared synergies, suggesting that EMGs of locomotor behaviors are generated primarily by centrally organized synergies. Both the amplitude and temporal patterns of the activation coefficients of most shared synergies, however, were altered by deafferentation, suggesting that sensory inflow modulates activation of those centrally organized synergies. For most synergies, effects of deafferentation on the activation coefficients were not consistent across frogs, indicating substantial interanimal variability of feedback actions. We speculate that sensory feedback might adapt recruitment of muscle synergies to behavioral constraints, and the few synergies specific to the intact or deafferented states might represent afferent-specific modules or feedback reorganization of spinal neuronal networks.

We developed non-negative factorization algorithms based on statistical distributions which are members of the exponential family, and using multiplicative update rules. We compared in detail the performance of algorithms derived using two particular exponential family distributions, assuming either constant variance noise (Gaussian) or signal dependent noise (gamma). These algorithms were compared on both simulated data sets and on muscle activation patterns collected from behaving animals. We found that on muscle activation patterns, which are expected to be corrupted by signal dependent noise, the factorizations identified by the algorithm assuming gamma distributed data were more robust than those identified by the algorithm assuming Gaussian distributed data.

Background Compensating for the effect of gravity by providing arm-weight support (WS) is a technique often utilized in the rehabilitation of patients with neurological conditions such as stroke to facilitate the performance of arm movements during therapy. Although it has been shown that, in healthy subjects as well as in stroke survivors, the use of arm WS during the performance of reaching movements leads to a general reduction, as expected, in the level of activation of upper limb muscles, the effects of different levels of WS on the characteristics of the kinematics of motion and of the activity of upper limb muscles have not been thoroughly investigated before. Methods In this study, we systematically assessed the characteristics of the kinematics of motion and of the activity of 14 upper limb muscles in a group of 9 healthy subjects who performed 3-D arm reaching movements while provided with different levels of arm WS. We studied the hand trajectory and the trunk, shoulder, and elbow joint angular displacement trajectories for different levels of arm WS. Besides, we analyzed the amplitude of the surface electromyographic (EMG) data collected from upper limb muscles and investigated patterns of coordination via the analysis of muscle synergies. Results The characteristics of the kinematics of motion varied across WS conditions but did not show distinct trends with the level of arm WS. The level of activation of upper limb muscles generally decreased, as expected, with the increase in arm WS. The same eight muscle synergies were identified in all WS conditions. Their level of activation depended on the provided level of arm WS. Conclusions The analysis of muscle synergies allowed us to identify a modular organization underlying the generation of arm reaching movements that appears to be invariant to the level of arm WS. The results of this study provide a normative dataset for the assessment of the effects of the level of arm WS on muscle synergies in stroke survivors and other patients who could benefit from upper limb rehabilitation with arm WS., National Institutes of Health (U.S.) (NIH grant # RC1-NS068103-01), Swiss National Science Foundation (grant PBEZ3_137336), National Institutes of Health (U.S.) (NIH grant # R24HD050821)

The experimental findings herein reported are aimed at gaining a perspective on the complex neural events that follow lesions of the motor cortical areas. Cortical damage, whether by trauma or stroke, interferes with the flow of descending signals to the modular interneuronal structures of the spinal cord. These spinal modules subserve normal motor behaviors by activating groups of muscles as individual units (muscle synergies). Damage to the motor cortical areas disrupts the orchestration of the modules, resulting in abnormal movements. To gain insights into this complex process, we recorded myoelectric signals from multiple upper-limb muscles in subjects with cortical lesions. We used a factorization algorithm to identify the muscle synergies. Our factorization analysis revealed, in a quantitative way, three distinct patterns of muscle coordination—including preservation, merging, and fractionation of muscle synergies—that reflect the multiple neural responses that occur after cortical damage. These patterns varied as a function of both the severity of functional impairment and the temporal distance from stroke onset. We think these muscle-synergy patterns can be used as physiological markers of the status of any patient with stroke or trauma, thereby guiding the development of different rehabilitation approaches, as well as future physiological experiments for a further understanding of postinjury mechanisms of motor control and recovery.