Search

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.

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.