Synergy temporal sequences and topography in the spinal cord: evidence for a traveling wave in frog locomotion

Why Firing Rate Distributions Are Important for Understanding Spinal Central Pattern Generators

Networks in the spinal cord, which are responsible for the generation of rhythmic movements, commonly known as central pattern generators (CPGs), have remained elusive for decades. Although it is well-known that many spinal neurons are rhythmically active, little attention has been given to the distribution of firing rates across the population. Here, we argue that firing rate distributions can provide an important clue to the organization of the CPGs. The data that can be gleaned from the sparse literature indicate a firing rate distribution, which is skewed toward zero with a long tail, akin to a normal distribution on a log-scale, i.e., a “log-normal” distribution. Importantly, such a shape is difficult to unite with the widespread assumption of modules composed of recurrently connected excitatory neurons. Spinal modules with recurrent excitation has the propensity to quickly escalate their firing rate and reach the maximum, hence equalizing the spiking activity across the population. The population distribution of firing rates hence would consist of a narrow peak near the maximum. This is incompatible with experiments, that show wide distributions and a peak close to zero. A way to resolve this puzzle is to include recurrent inhibition internally in each CPG modules. Hence, we investigate the impact of recurrent inhibition in a model and find that the firing rate distributions are closer to the experimentally observed. We therefore propose that recurrent inhibition is a crucial element in motor circuits, and suggest that future models of motor circuits should include recurrent inhibition as a mandatory element.

Hierarchical control of locomotion by distinct types of spinal V2a interneurons in zebrafish

In all vertebrates, excitatory spinal interneurons execute dynamic adjustments in the timing and amplitude of locomotor movements. Currently, it is unclear whether interneurons responsible for timing control are distinct from those involved in amplitude control. Here, we show that in larval zebrafish, molecularly, morphologically and electrophysiologically distinct types of V2a neurons exhibit complementary patterns of connectivity. Stronger higher-order connections from type I neurons to other excitatory V2a and inhibitory V0d interneurons provide timing control, while stronger last-order connections from type II neurons to motor neurons provide amplitude control. Thus, timing and amplitude are coordinated by distinct interneurons distinguished not by their occupation of hierarchically-arranged anatomical layers, but rather by differences in the reliability and probability of higher-order and last-order connections that ultimately form a single anatomical layer. These findings contribute to our understanding of the origins of timing and amplitude control in the spinal cord.

V2a excitatory interneurons in the spinal cord are important for coordinating locomotion. Here the authors describe two types of V2a neuron with differences in higher order and lower order connectivity in larval zebrafish.

Differential activation of lumbar and sacral motor pools during walking at different speeds and slopes

Organization of spinal motor output has become of interest for investigating differential activation of lumbar and sacral motor pools during locomotor tasks. Motor pools are associated with functional grouping of motoneurons of the lower limb muscles. Here we examined how the spatiotemporal organization of lumbar and sacral motor pool activity during walking is orchestrated with slope of terrain and speed of progression. Ten subjects walked on an instrumented treadmill at different slopes and imposed speeds. Kinetics, kinematics, and electromyography of 16 lower limb muscles were recorded. The spinal locomotor output was assessed by decomposing the coordinated muscle activation profiles into a small set of common factors and by mapping them onto the rostrocaudal location of the motoneuron pools. Our results show that lumbar and sacral motor pool activity depend on slope and speed. Compared with level walking, sacral motor pools decrease their activity at negative slopes and increase at positive slopes, whereas lumbar motor pools increase their engagement when both positive and negative slope increase. These findings are consistent with a differential involvement of the lumbar and the sacral motor pools in relation to changes in positive and negative center of body mass mechanical power production due to slope and speed.NEW & NOTEWORTHY In this study, the spatiotemporal maps of motoneuron activity in the spinal cord were assessed during walking at different slopes and speeds. We found differential involvement of lumbar and sacral motor pools in relation to changes in positive and negative center of body mass power production due to slope and speed. The results are consistent with recent findings about the specialization of neuronal networks located at different segments of the spinal cord for performing specific locomotor tasks.

Cortical Correlates of Locomotor Muscle Synergy Activation in Humans: An Electroencephalographic Decoding Study

Muscular control during walking is believed to be simplified by the coactivation of muscles called muscle synergies. Although significant corticomuscular connectivity during walking has been reported, the level at which the cortical activity is involved in muscle activity (muscle synergy or individual muscle level) remains unclear. Here we examined cortical correlates of muscle activation during walking by brain decoding of activation of muscle synergies and individual muscles from electroencephalographic signals. We demonstrated that the activation of locomotor muscle synergies was decoded from slow cortical waves. In addition, the decoding accuracy for muscle synergies was greater than that for individual muscles and the decoding of individual muscle activation was based on muscle-synergy-related cortical information. These results indicate the cortical correlates of locomotor muscle synergy activation. These findings expand our understanding of the relationships between brain and locomotor muscle synergies and could accelerate the development of effective brain-machine interfaces for walking rehabilitation.

Adaptive hindlimb split-belt treadmill walking in rats by controlling basic muscle activation patterns via phase resetting

To investigate the adaptive locomotion mechanism in animals, a split-belt treadmill has been used, which has two parallel belts to produce left–right symmetric and asymmetric environments for walking. Spinal cats walking on the treadmill have suggested the contribution of the spinal cord and associated peripheral nervous system to the adaptive locomotion. Physiological studies have shown that phase resetting of locomotor commands involving a phase shift occurs depending on the types of sensory nerves and stimulation timing, and that muscle activation patterns during walking are represented by a linear combination of a few numbers of basic temporal patterns despite the complexity of the activation patterns. Our working hypothesis was that resetting the onset timings of basic temporal patterns based on the sensory information from the leg, especially extension of hip flexors, contributes to adaptive locomotion on the split-belt treadmill. Our hypothesis was examined by conducting forward dynamic simulations using a neuromusculoskeletal model of a rat walking on a split-belt treadmill with its hindlimbs and by comparing the simulated motions with the measured motions of rats.

Synergies in coordination: a comprehensive overview of neural, computational, and behavioral approaches

At face value, the term "synergy" provides a unifying concept within a fractured field that encompasses complementary neural, computational, and behavioral approaches. However, the term is not used synonymously by different researchers but has substantially different meanings depending on the research approach. With so many operational definitions for the one term, it becomes difficult to use as either a descriptive or explanatory concept, yet it remains pervasive and apparently indispensable. Here we provide a summary of different approaches that invoke synergies in a descriptive or explanatory context, summarizing progress, not within the one approach, but across the theoretical landscape. Bernstein's framework of flexible hierarchical control may provide a unifying framework here, since it can incorporate divergent ideas about synergies. In the current motor control literature, synergy may refer to conceptually different processes that could potentially operate in parallel, across different levels within the same hierarchical control scheme. There is evidence for the concurrent existence of synergies with different features, both "hard-wired" and "soft-wired," and task independent and task dependent. By providing a comprehensive overview of the multifaceted ideas about synergies, our goal is to move away from the compartmentalization and narrow the focus on one level and promote a broader perspective on the control and coordination of movement.

Muscle Synergies Obtained from Comprehensive Mapping of the Cortical Forelimb Representation Using Stimulus Triggered Averaging of EMG Activity

Neuromuscular control of voluntary movement may be simplified using muscle synergies similar to those found using non-negative matrix factorization. We recently identified synergies in electromyography (EMG) recordings associated with both voluntary movement and movement evoked by high-frequency long-duration intracortical microstimulation applied to the forelimb representation of the primary motor cortex (M1). The goal of this study was to use stimulus-triggered averaging (StTA) of EMG activity to investigate the synergy profiles and weighting coefficients associated with poststimulus facilitation, as synergies may be hard-wired into elemental cortical output modules and revealed by StTA. We applied StTA at low (LOW, ∼15 μA) and high intensities (HIGH, ∼110 μA) to 247 cortical locations of the M1 forelimb region in two male rhesus macaques while recording the EMG of 24 forelimb muscles. Our results show that 10-11 synergies accounted for 90% of the variation in poststimulus EMG facilitation peaks from the LOW-intensity StTA dataset while only 4-5 synergies were needed for the HIGH-intensity dataset. Synergies were similar across monkeys and current intensities. Most synergy profiles strongly activated only one or two muscles; all joints were represented and most, but not all, joint directions of motion were represented. Cortical maps of the synergy weighting coefficients suggest only a weak organization. StTA of M1 resulted in highly diverse muscle activations, suggestive of the limiting condition of requiring a synergy for each muscle to account for the patterns observed.SIGNIFICANCE STATEMENT Coordination of muscle activity and the neural origin of potential muscle synergies remains a fundamental question of neuroscience. We previously demonstrated that high-frequency long-duration intracortical microstimulation-evoked synergies were unrelated to voluntary movement synergies and were not clearly organized in the cortex. Here we present stimulus-triggered averaging facilitation-related muscle synergies, suggesting that when fundamental cortical output modules are activated, synergies approach the limit of single-muscle control. Thus, we conclude that if the CNS controls movement via linear synergies, those synergies are unlikely to be called from M1. This information is critical for understanding neural control of movement and the development of brain-machine interfaces.

Speed dependency in α-motoneuron activity and locomotor modules in human locomotion: indirect evidence for phylogenetically conserved spinal circuits

Coordinated locomotor muscle activity is generated by the spinal central pattern generators (CPGs). Vertebrate studies have demonstrated the following two characteristics of the speed control mechanisms of the spinal CPGs: (i) rostral segment activation is indispensable for achieving high-speed locomotion; and (ii) specific combinations between spinal interneuronal modules and motoneuron (MN) pools are sequentially activated with increasing speed. Here, to investigate whether similar control mechanisms exist in humans, we examined spinal neural activity during varied-speed locomotion by mapping the distribution of MN activity in the spinal cord and extracting locomotor modules, which generate basic MN activation patterns. The MN activation patterns and the locomotor modules were analysed from multi-muscle electromyographic recordings. The reconstructed MN activity patterns were divided into the following three patterns depending on the speed of locomotion: slow walking, fast walking and running. During these three activation patterns, the proportion of the activity in rostral segments to that in caudal segments increased as locomotion speed increased. Additionally, the different MN activation patterns were generated by distinct combinations of locomotor modules. These results are consistent with the speed control mechanisms observed in vertebrates, suggesting phylogenetically conserved spinal mechanisms of neural control of locomotion.

Evidence for sparse synergies in grasping actions

Converging evidence shows that hand-actions are controlled at the level of synergies and not single muscles. One intriguing aspect of synergy-based action-representation is that it may be intrinsically sparse and the same synergies can be shared across several distinct types of hand-actions. Here, adopting a normative angle, we consider three hypotheses for hand-action optimal-control: sparse-combination hypothesis (SC) – sparsity in the mapping between synergies and actions - i.e., actions implemented using a sparse combination of synergies; sparse-elements hypothesis (SE) – sparsity in synergy representation – i.e., the mapping between degrees-of-freedom (DoF) and synergies is sparse; double-sparsity hypothesis (DS) – a novel view combining both SC and SE – i.e., both the mapping between DoF and synergies and between synergies and actions are sparse, each action implementing a sparse combination of synergies (as in SC), each using a limited set of DoFs (as in SE). We evaluate these hypotheses using hand kinematic data from six human subjects performing nine different types of reach-to-grasp actions. Our results support DS, suggesting that the best action representation is based on a relatively large set of synergies, each involving a reduced number of degrees-of-freedom, and that distinct sets of synergies may be involved in distinct tasks.

Critical Points and Traveling Wave in Locomotion: Experimental Evidence and Some Theoretical Considerations

The central pattern generator (CPG) architecture for rhythm generation remains partly elusive. We compare cat and frog locomotion results, where the component unrelated to pattern formation appears as a temporal grid, and traveling wave respectively. Frog spinal cord microstimulation with N-methyl-D-Aspartate (NMDA), a CPG activator, produced a limited set of force directions, sometimes tonic, but more often alternating between directions similar to the tonic forces. The tonic forces were topographically organized, and sites evoking rhythms with different force subsets were located close to the constituent tonic force regions. Thus CPGs consist of topographically organized modules. Modularity was also identified as a limited set of muscle synergies whose combinations reconstructed the EMGs. The cat CPG was investigated using proprioceptive inputs during fictive locomotion. Critical points identified both as abrupt transitions in the effect of phasic perturbations, and burst shape transitions, had biomechanical correlates in intact locomotion. During tonic proprioceptive perturbations, discrete shifts between these critical points explained the burst durations changes, and amplitude changes occurred at one of these points. Besides confirming CPG modularity, these results suggest a fixed temporal grid of anchoring points, to shift modules onsets and offsets. Frog locomotion, reconstructed with the NMDA synergies, showed a partially overlapping synergy activation sequence. Using the early synergy output evoked by NMDA at different spinal sites, revealed a rostrocaudal topographic organization, where each synergy is preferentially evoked from a few, albeit overlapping, cord regions. Comparing the locomotor synergy sequence with this topography suggests that a rostrocaudal traveling wave would activate the synergies in the proper sequence for locomotion. This output was reproduced in a two-layer model using this topography and a traveling wave. Together our results suggest two CPG components: modules, i.e., synergies; and temporal patterning, seen as a temporal grid in the cat, and a traveling wave in the frog. Animal and limb navigation have similarities. Research relating grid cells to the theta rhythm and on segmentation during navigation may relate to our temporal grid and traveling wave results. Winfree's mathematical work, combining critical phases and a traveling wave, also appears important. We conclude suggesting tracing, and imaging experiments to investigate our CPG model.

Motor module activation sequence and topography in the spinal cord during air-stepping in human: Insights into the traveling wave in spinal locomotor circuits

Coordinated locomotor muscle activity is generated by the spinal central pattern generators (CPGs), which are modulated by peripheral and supraspinal inputs. The CPGs would consist of multiple motor modules generating basic muscle activity, which are distributed rostrocaudally along the spinal cord. To activate the motor modules in proper sequence, rostrocaudally traveling waves of activation in the spinal cord are important mechanisms in the CPGs. The traveling waves of activation have been observed in nonhuman vertebrates. However, they have not yet been confirmed during human locomotion. Although, rostrocaudal wave‐like activations in the spinal cord were observed during walking in humans in a previous study, the propagation shifted rostrally toward the upper lumbar segments at foot contact. Here, using an air stepping task to remove the foot‐contact interactions, we examined whether the traveling wave mechanism exists in the human spinal circuits based on the activation sequence of motor modules and their topography. We measured electromyographic activity of lower leg muscles during the air‐stepping task. Then, we extracted motor modules (i.e., basic patterns of sets of muscle activations: muscle synergies) from the measured muscle activities using nonnegative matrix factorization method. Next, we reconstructed motoneuron (MN) activity from each module activity based on myotomal charts. We identified four types of motor modules from muscle activities during the air‐stepping task. Each motor module represented different sets of synergistic muscle activations. MN clusters innervating each motor module were sequentially activated from the rostral to caudal region in the spinal cord, from the initial flexion to the last extension phase during air‐stepping. The rostrocaudally sequential activation of MN clusters suggests the possibility that rostrocaudally traveling waves exist in human locomotor spinal circuits. The present results advance the understanding of human locomotor control mechanisms, and provide important insights into the evolution of locomotor networks in vertebrates.

And yet it moves: Recovery of volitional control after spinal cord injury

Preclinical and clinical neurophysiological and neurorehabilitation research has generated rather surprising levels of recovery of volitional sensory-motor function in persons with chronic motor paralysis following a spinal cord injury. The key factor in this recovery is largely activity-dependent plasticity of spinal and supraspinal networks. This key factor can be triggered by neuromodulation of these networks with electrical and pharmacological interventions. This review addresses some of the systems-level physiological mechanisms that might explain the effects of electrical modulation and how repetitive training facilitates the recovery of volitional motor control. In particular, we substantiate the hypotheses that: (1) in the majority of spinal lesions, a critical number and type of neurons in the region of the injury survive, but cannot conduct action potentials, and thus are electrically non-responsive; (2) these neuronal networks within the lesioned area can be neuromodulated to a transformed state of electrical competency; (3) these two factors enable the potential for extensive activity-dependent reorganization of neuronal networks in the spinal cord and brain, and (4) propriospinal networks play a critical role in driving this activity-dependent reorganization after injury. Real-time proprioceptive input to spinal networks provides the template for reorganization of spinal networks that play a leading role in the level of coordination of motor pools required to perform a given functional task. Repetitive exposure of multi-segmental sensory-motor networks to the dynamics of task-specific sensory input as occurs with repetitive training can functionally reshape spinal and supraspinal connectivity thus re-enabling one to perform complex motor tasks, even years post injury.

Muscle synergies obtained from comprehensive mapping of the primary motor cortex forelimb representation using high-frequency, long-duration ICMS

Simplifying neuromuscular control for movement has previously been explored by extracting muscle synergies from voluntary movement electromyography (EMG) patterns. The purpose of this study was to investigate muscle synergies represented in EMG recordings associated with direct electrical stimulation of single sites in primary motor cortex (M1). We applied single-electrode high-frequency, long-duration intracortical microstimulation (HFLD-ICMS) to the forelimb region of M1 in two rhesus macaques using parameters previously found to produce forelimb movements to stable spatial end points (90-150 Hz, 90-150 μA, 1,000-ms stimulus train lengths). To develop a comprehensive representation of cortical output, stimulation was applied systematically across the full extent of M1. We recorded EMG activity from 24 forelimb muscles together with movement kinematics. Nonnegative matrix factorization (NMF) was applied to the mean stimulus-evoked EMG, and the weighting coefficients associated with each synergy were mapped to the cortical location of the stimulating electrode. Synergies were found for three data sets including 1) all stimulated sites in the cortex, 2) a subset of sites that produced stable movement end points, and 3) EMG activity associated with voluntary reaching. Two or three synergies accounted for 90% of the overall variation in voluntary movement EMG whereas four or five synergies were needed for HFLD-ICMS-evoked EMG data sets. Maps of the weighting coefficients from the full HFLD-ICMS data set show limited regional areas of higher activation for particular synergies. Our results demonstrate fundamental NMF-based muscle synergies in the collective M1 output, but whether and how the central nervous system might coordinate movements using these synergies remains unclear.NEW & NOTEWORTHY While muscle synergies have been investigated in various muscle activity sets, it is unclear whether and how synergies may be organized in the cortex. We have investigated muscle synergies resulting from high-frequency, long-duration intracortical microstimulation (HFLD-ICMS) applied throughout M1. We compared HFLD-ICMS synergies to synergies from voluntary movement. While synergies can be identified from M1 stimulation, they are not clearly related to voluntary movement synergies and do not show an orderly topographic organization across M1.

An Optogenetic Demonstration of Motor Modularity in the Mammalian Spinal Cord

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.