Distinguishing synchronous and time-varying synergies using point process interval statistics: motor primitives in frog and rat

Muscle synergies are associated with intermuscular coherence and cortico-synergy coherence in an isometric upper limb task

To elucidate the underlying physiological mechanisms of muscle synergies, we investigated long-range functional connectivity by cortico-muscular (CMC), intermuscular (IMC) and cortico-synergy (CSC) coherence. Fourteen healthy participants executed an isometric upper limb task in synergy-tuned directions. Cortical activity was recorded using 32-channel electroencephalography (EEG) and muscle activity using 16-channel electromyography (EMG). Using non-negative matrix factorisation (NMF), we calculated muscle synergies from two different tasks. A preliminary multidirectional task was used to identify synergy-preferred directions (PDs). A subsequent coherence task, consisting of generating forces isometrically in the synergy PDs, was used to assess the functional connectivity properties of synergies. Overall, we were able to identify four different synergies from the multidirectional task. A significant alpha band IMC was consistently present in all extracted synergies. Moreover, IMC alpha band was higher between muscles with higher weights within a synergy. Interestingly, CSC alpha band was also significantly higher across muscles with higher weights within a synergy. In contrast, no significant CMC was found between the motor cortex area and synergy muscles. The presence of a shared input onto synergistic muscles within a synergy supports the idea of neurally derived muscle synergies that build human movement. Our findings suggest cortical modulation of some of the synergies and the consequential existence of shared input between muscles within cortically modulated synergies.

Using different matrix factorization approaches to identify muscle synergy in stroke survivors

Over the past several decades, many scholars have investigated muscle synergy as a promising tool for evaluating motor function. However, it is challenging to obtain favorable robustness using the general muscle synergy identification algorithms, namely non-negative matrix factorization (NMF), independent component analysis (ICA), and factor analysis (FA). Some scholars have proposed improved muscle synergy identification algorithms to overcome the shortcomings of these approaches, such as singular value decomposition NMF (SVD-NMF), sparse NMF (S-NMF), and multivariate curve resolution-alternating least squares (MCR-ALS). However, performance comparisons of these algorithms are seldom conducted. In this study, experimental electromyography (EMG) data collected from healthy individuals and stroke survivors were applied to assess the repeatability and intra-subject consistency of NMF, SVD-NMF, S-NMF, ICA, FA, and MCR-ALS. MCR-ALS presented higher repeatability and intra-subject consistencies than the other algorithms. More synergies and lower intra-subject consistencies were observed in stroke survivors than in healthy individuals. Thus, MCR-ALS is considered a favorable muscle synergy identification algorithm for patients with neural system disorders.

Evaluation of Synergy-Based Hand Gesture Recognition Method Against Force Variation for Robust Myoelectric Control

The non-stationary characteristics of surface electromyography (sEMG) and possible adverse variations in real-world conditions make it still an open challenge to realize robust myoelectric control (MEC) for multifunctional prostheses. Variable muscle contraction level is one of the handicaps that may degrade the performance of MEC. In this study, we proposed a force-invariant intent recognition method based on muscle synergy analysis (MSA) in the setting of three self-defined force levels (low, medium, and high). Specifically, a fast matrix factorization algorithm based on alternating non-negativity constrained least squares (NMF/ANLS) was chosen to extract task-specific synergies associated with each of six hand gestures in the training stage; while for the testing samples, we used the non-negative least square (NNLS) method to estimate neural commands for movement classification. The performance of proposed method was compared with conventional pattern recognition (PR) method consisting of LDA (linear discrimination analysis) classifier and representative features in three offline evaluation scenarios. Statistical tests on ten able-bodied subjects revealed no significant difference in intra-force-level (p = 0.353) and multi-force-level (p = 0.695) accuracy; But the synergy-based method performed significantly better than conventional PR-based method under inter-force-level conditions (p < 0.05). Similar results were observed for nine amputee subjects though there was a drop in the classification accuracy. This study was the first to concurrently demonstrate the robustness and predictive power of task-specific synergies under variant force levels and explore their potential for reliable intent recognition against force variation. Although the online performance is yet to be demonstrated, the proposed method is characterized by simple training procedure and acceptable computational efficiency, which would potentially provide an alternative approach for the development of clinically viable prostheses and rehabilitation robots driven by sEMG.

A Novel Muscle Synergy Extraction Method Used for Motor Function Evaluation of Stroke Patients: A Pilot Study

In this paper, we present a novel muscle synergy extraction method based on multivariate curve resolution–alternating least squares (MCR-ALS) to overcome the limitation of the nonnegative matrix factorization (NMF) method for extracting non-sparse muscle synergy, and we study its potential application for evaluating motor function of stroke survivors. Nonnegative matrix factorization (NMF) is the most widely used method for muscle synergy extraction. However, NMF is susceptible to components’ sparseness and usually provides inferior reliability, which significantly limits the promotion of muscle synergy. In this study, MCR-ALS was employed to extract muscle synergy from electromyography (EMG) data. Its performance was compared with two other matrix factorization algorithms, NMF and self-modeling mixture analysis (SMMA). Simulated data sets were utilized to explore the influences of the sparseness and noise on the extracted synergies. As a result, the synergies estimated by MCR-ALS were the most similar to true synergies as compared with SMMA and NMF. MCR-ALS was used to analyze the muscle synergy characteristics of upper limb movements performed by healthy (n = 11) and stroke (n = 5) subjects. The repeatability and intra-subject consistency were used to evaluate the performance of MCR-ALS. As a result, MCR-ALS provided much higher repeatability and intra-subject consistency as compared with NMF, which were important for the reliability of the motor function evaluation. The stroke subjects had lower intra-subject consistency and seemingly had more synergies as compared with the healthy subjects. Thus, MCR-ALS is a promising muscle synergy analysis method for motor function evaluation of stroke patients.

A three-dimensional musculoskeletal model of the dog

The domestic dog is interesting to investigate because of the wide range of body size, body mass, and physique in the many breeds. In the last several years, the number of clinical and biomechanical studies on dog locomotion has increased. However, the relationship between body structure and joint load during locomotion, as well as between joint load and degenerative diseases of the locomotor system (e.g. dysplasia), are not sufficiently understood. Collecting this data through in vivo measurements/records of joint forces and loads on deep/small muscles is complex, invasive, and sometimes unethical. The use of detailed musculoskeletal models may help fill the knowledge gap. We describe here the methods we used to create a detailed musculoskeletal model with 84 degrees of freedom and 134 muscles. Our model has three key-features: three-dimensionality, scalability, and modularity. We tested the validity of the model by identifying forelimb muscle synergies of a walking Beagle. We used inverse dynamics and static optimization to estimate muscle activations based on experimental data. We identified three muscle synergy groups by using hierarchical clustering. The activation patterns predicted from the model exhibit good agreement with experimental data for most of the forelimb muscles. We expect that our model will speed up the analysis of how body size, physique, agility, and disease influence neuronal control and joint loading in dog locomotion.

Non-negative matrix factorisation is the most appropriate method for extraction of muscle synergies in walking and running

Muscle synergies provide a simple description of a complex motor control mechanism. Synergies are extracted from muscle activation patterns using factorisation methods. Despite the availability of several factorisation methods in the literature, the most appropriate method for muscle synergy extraction is currently unknown. In this study, we compared four muscle synergy extraction methods: non-negative matrix factorisation, principal component analysis, independent component analysis, and factor analysis. Probability distribution of muscle activation patterns were compared with the probability distribution of synergy excitation primitives obtained from the four factorisation methods. Muscle synergies extracted using non-negative matrix factorisation best matched the probability distribution of muscle activation patterns across different walking and running speeds. Non-negative matrix factorisation also best tracked changes in muscle activation patterns compared to the other factorisation methods. Our results suggest that non-negative matrix factorisation is the best factorisation method for identifying muscle synergies in dynamic tasks with different levels of muscle contraction.

Neuromechanical Model of Rat Hindlimb Walking with Two-Layer CPGs

This work demonstrates a neuromechanical model of rat hindlimb locomotion undergoing nominal walking with perturbations. In the animal, two types of responses to perturbations are observed: resetting and non-resetting deletions. This suggests that the animal locomotor system contains a memory-like organization. To model this phenomenon, we built a synthetic nervous system that uses separate rhythm generator and pattern formation layers to activate antagonistic muscle pairs about each joint in the sagittal plane. Our model replicates the resetting and non-resetting deletions observed in the animal. In addition, in the intact (i.e., fully afferented) rat walking simulation, we observe slower recovery after perturbation, which is different from the deafferented animal experiment. These results demonstrate that our model is a biologically feasible description of some of the neural circuits in the mammalian spinal cord that control locomotion, and the difference between our simulation and fictive motion shows the importance of sensory feedback on motor output. This model also demonstrates how the pattern formation network can activate muscle synergies in a coordinated way to produce stable walking, which motivates the use of more complex synergies activating more muscles in the legs for three-dimensional limb motion.

Making the Coupled Gaussian Process Dynamical Model Modular and Scalable with Variational Approximations

We describe a sparse, variational posterior approximation to the Coupled Gaussian Process Dynamical Model (CGPDM), which is a latent space coupled dynamical model in discrete time. The purpose of the approximation is threefold: first, to reduce training time of the model; second, to enable modular re-use of learned dynamics; and, third, to store these learned dynamics compactly. Our target applications here are human movement primitive (MP) models, where an MP is a reusable spatiotemporal component, or "module" of a human full-body movement. Besides re-usability of learned MPs, compactness is crucial, to allow for the storage of a large library of movements. We first derive the variational approximation, illustrate it on toy data, test its predictions against a range of other MP models and finally compare movements produced by the model against human perceptual expectations. We show that the variational CGPDM outperforms several other MP models on movement trajectory prediction. Furthermore, human observers find its movements nearly indistinguishable from replays of natural movement recordings for a very compact parameterization of the approximation.

A Systematic Review on Muscle Synergies: From Building Blocks of Motor Behavior to a Neurorehabilitation Tool

The central nervous system (CNS) is believed to utilize specific predefined modules, called muscle synergies (MS), to accomplish a motor task. Yet questions persist about how the CNS combines these primitives in different ways to suit the task conditions. The MS hypothesis has been a subject of debate as to whether they originate from neural origins or nonneural constraints. In this review article, we present three aspects related to the MS hypothesis: (1) the experimental and computational evidence in support of the existence of MS, (2) algorithmic approaches for extracting them from surface electromyography (EMG) signals, and (3) the possible role of MS as a neurorehabilitation tool. We note that recent advances in computational neuroscience have utilized the MS hypothesis in motor control and learning. Prospective advances in clinical, medical, and engineering sciences and in fields such as robotics and rehabilitation stand to benefit from a more thorough understanding of MS.

Model selection for the extraction of movement primitives

A wide range of blind source separation methods have been used in motor control research for the extraction of movement primitives from EMG and kinematic data. Popular examples are principal component analysis (PCA), independent component analysis (ICA), anechoic demixing, and the time-varying synergy model (d'Avella and Tresch, 2002). However, choosing the parameters of these models, or indeed choosing the type of model, is often done in a heuristic fashion, driven by result expectations as much as by the data. We propose an objective criterion which allows to select the model type, number of primitives and the temporal smoothness prior. Our approach is based on a Laplace approximation to the posterior distribution of the parameters of a given blind source separation model, re-formulated as a Bayesian generative model. We first validate our criterion on ground truth data, showing that it performs at least as good as traditional model selection criteria [Bayesian information criterion, BIC (Schwarz, 1978) and the Akaike Information Criterion (AIC) (Akaike, 1974)]. Then, we analyze human gait data, finding that an anechoic mixture model with a temporal smoothness constraint on the sources can best account for the data.