Modelling spinal circuitry involved in locomotor pattern generation: insights from the effects of afferent stimulation

Operation regimes of spinal circuits controlling locomotion and the role of supraspinal drives and sensory feedback

Locomotion in mammals is directly controlled by the spinal neuronal network, operating under the control of supraspinal signals and somatosensory feedback that interact with each other. However, the functional architecture of the spinal locomotor network, its operation regimes, and the role of supraspinal and sensory feedback in different locomotor behaviors, including at different speeds, remain unclear. We developed a computational model of spinal locomotor circuits receiving supraspinal drives and limb sensory feedback that could reproduce multiple experimental data obtained in intact and spinal-transected cats during tied-belt and split-belt treadmill locomotion. We provide evidence that the spinal locomotor network operates in different regimes depending on locomotor speed. In an intact system, at slow speeds (<0.4 m/s), the spinal network operates in a non-oscillating state-machine regime and requires sensory feedback or external inputs for phase transitions. Removing sensory feedback related to limb extension prevents locomotor oscillations at slow speeds. With increasing speed and supraspinal drives, the spinal network switches to a flexor-driven oscillatory regime and then to a classical half-center regime. Following spinal transection, the model predicts that the spinal network can only operate in the state-machine regime. Our results suggest that the spinal network operates in different regimes for slow exploratory and fast escape locomotor behaviors, making use of different control mechanisms.

EFFECTS OF SPINAL TRANSECTION AND LOCOMOTOR SPEED ON MUSCLE SYNERGIES OF THE CAT HINDLIMB

It was suggested that during locomotion, the nervous system controls movement by activating groups of muscles, or muscle synergies. Analysis of muscle synergies can reveal the organization of spinal locomotor networks and how it depends on the state of the nervous system, such as before and after spinal cord injury, and on different locomotor conditions, including a change in speed. The goal of this study was to investigate the effects of spinal transection and locomotor speed on hindlimb muscle synergies and their time-dependent activity patterns in adult cats. EMG activities of 15 hindlimb muscles were recorded in 9 adult cats of either sex during tied-belt treadmill locomotion at speeds of 0.4, 0.7, and 1.0 m/s before and after recovery from a low thoracic spinal transection. We determined EMG burst groups using cluster analysis of EMG burst onset and offset times and muscle synergies using non-negative matrix factorization. We found five major EMG burst groups and five muscle synergies in each of six experimental conditions (2 states × 3 speeds). In each case, the synergies accounted for at least 90% of muscle EMG variance. Both spinal transection and locomotion speed modified subgroups of EMG burst groups and the composition and activation patterns of selected synergies. However, these changes did not modify the general organization of muscle synergies. Based on the obtained results, we propose an organization for a pattern formation network of a two-level central pattern generator that can be tested in neuromechanical simulations of spinal circuits controlling cat locomotion.

Regeneration and Plasticity Induced by Epidural Stimulation in a Rodent Model of Spinal Cord Injury

Traumatic spinal cord injury is a major cause of disability for which there are currently no fully effective treatments. Recent studies using epidural electrical stimulation have shown significant advances in motor rehabilitation, even when applied during chronic phases of the disease. The present study aimed to investigate the effectiveness of epidural electric stimulation in the motor recovery of rats with spinal cord injury. Furthermore, we aimed to elucidate the neurophysiological mechanisms underlying motor recovery. First, we improved upon the impact spinal cord injury model to cause severe and permanent motor deficits lasting up to 2 months. Next, we developed and tested an implantable epidural spinal cord stimulator device for rats containing an electrode and an implantable generator. Finally, we evaluated the efficacy of epidural electrical stimulation on motor recovery after spinal cord injury in Wistar rats. A total of 60 animals were divided into the following groups: (i) severe injury with epidural electrical stimulation (injury + stim, n = 15), (ii) severe injury without stimulation (group injury, n = 15), (iii) sham implantation without battery (sham, n = 15), and (iv) a control group, without surgical intervention (control, n = 15). All animals underwent weekly evaluations using the Basso, Beattie, Bresnahan (BBB) locomotor rating scale index, inclined plane, and OpenField test starting one week before the lesion and continuing for eight weeks. After this period, the animals were sacrificed and their spinal cords were explanted and prepared for histological analysis (hematoxylin-eosin) and immunohistochemistry for NeuN, β-III-tubulin, synaptophysin, and Caspase 3. Finally, NeuN-positive neuronal nuclei were quantified through stereology; fluorescence signal intensities for β-tubulin, synaptophyin, and Caspase 3 were quantified using an epifluorescence microscope. The injury + stim group showed significant improvement on the BBB scale compared with the injured group after the 5th week (p < 0.05). Stereological analysis showed a significantly higher average count of neural cells in the injury + stim group in relation to the injury group (1783 ± 2 vs. 897 ± 3, p < 0.001). Additionally, fluorescence signal intensity for synaptophysin was significantly higher in the injury + stim group in relation to the injury group (1294 ± 46 vs. 1198 ± 23, p < 0.01); no statistically significant difference was found in β-III-tubulin signal intensity. Finally, Caspase 3 signal intensity was significantly lower in the stim group (727 ± 123) compared with the injury group (1225 ± 87 p < 0.05), approaching levels observed in the sham and control groups. Our data suggest a regenerative and protective effect of epidural electrical stimulation in rats subjected to impact-induced traumatic spinal cord injury.

Operation regimes of spinal circuits controlling locomotion and role of supraspinal drives and sensory feedback

Locomotion in mammals is directly controlled by the spinal neuronal network, operating under the control of supraspinal signals and somatosensory feedback that interact with each other. However, the functional architecture of the spinal locomotor network, its operation regimes, and the role of supraspinal and sensory feedback in different locomotor behaviors, including at different speeds, remain unclear. We developed a computational model of spinal locomotor circuits receiving supraspinal drives and limb sensory feedback that could reproduce multiple experimental data obtained in intact and spinal-transected cats during tied-belt and split-belt treadmill locomotion. We provide evidence that the spinal locomotor network operates in different regimes depending on locomotor speed. In an intact system, at slow speeds (< 0.4 m/s), the spinal network operates in a non-oscillating state-machine regime and requires sensory feedback or external inputs for phase transitions. Removing sensory feedback related to limb extension prevents locomotor oscillations at slow speeds. With increasing speed and supraspinal drives, the spinal network switches to a flexor-driven oscillatory regime and then to a classical half-center regime. Following spinal transection, the model predicts that the spinal network can only operate in the state-machine regime. Our results suggest that the spinal network operates in different regimes for slow exploratory and fast escape locomotor behaviors, making use of different control mechanisms.

Rare phenomena of central rhythm and pattern generation in a case of complete spinal cord injury

Lumbar central pattern generators (CPGs) control the basic rhythm and coordinate muscle activation underlying hindlimb locomotion in quadrupedal mammals. The existence and function of CPGs in humans have remained controversial. Here, we investigated a case of a male individual with complete thoracic spinal cord injury who presented with a rare form of self-sustained rhythmic spinal myoclonus in the legs and rhythmic activities induced by epidural electrical stimulation (EES). Analysis of muscle activation patterns suggested that the myoclonus tapped into spinal circuits that generate muscle spasms, rather than reflecting locomotor CPG activity as previously thought. The EES-induced patterns were fundamentally different in that they included flexor-extensor and left-right alternations, hallmarks of locomotor CPGs, and showed spontaneous errors in rhythmicity. These motor deletions, with preserved cycle frequency and period when rhythmic activity resumed, were previously reported only in animal studies and suggest a separation between rhythm generation and pattern formation. Spinal myoclonus and the EES-induced activity demonstrate that the human lumbar spinal cord contains distinct mechanisms for generating rhythmic multi-muscle patterns.

Spinal cords: Symphonies of interneurons across species

Vertebrate movement is orchestrated by spinal inter- and motor neurons that, together with sensory and cognitive input, produce dynamic motor behaviors. These behaviors vary from the simple undulatory swimming of fish and larval aquatic species to the highly coordinated running, reaching and grasping of mice, humans and other mammals. This variation raises the fundamental question of how spinal circuits have changed in register with motor behavior. In simple, undulatory fish, exemplified by the lamprey, two broad classes of interneurons shape motor neuron output: ipsilateral-projecting excitatory neurons, and commissural-projecting inhibitory neurons. An additional class of ipsilateral inhibitory neurons is required to generate escape swim behavior in larval zebrafish and tadpoles. In limbed vertebrates, a more complex spinal neuron composition is observed. In this review, we provide evidence that movement elaboration correlates with an increase and specialization of these three basic interneuron types into molecularly, anatomically, and functionally distinct subpopulations. We summarize recent work linking neuron types to movement-pattern generation across fish, amphibians, reptiles, birds and mammals.

Combining Reflexes and External Sensory Information in a Neuromusculoskeletal Model to Control a Quadruped Robot

This article examines the importance of integrating locomotion and cognitive information for achieving dynamic locomotion from a viewpoint combining biology and ecological psychology. We present a mammalian neuromusculoskeletal model from external sensory information processing to muscle activation, which includes: 1) a visual-attention control mechanism for controlling attention to external inputs; 2) object recognition representing the primary motor cortex; 3) a motor control model that determines motor commands traveling down the corticospinal and reticulospinal tracts; 4) a central pattern generation model representing pattern generation in the spinal cord; and 5) a muscle reflex model representing the muscle model and its reflex mechanism. The proposed model is able to generate the locomotion of a quadruped robot in flat and natural terrain. The experiment also shows the importance of a postural reflex mechanism when experiencing a sudden obstacle. We show the reflex mechanism when a sudden obstacle is separately detected from both external (retina) and internal (touching afferent) sensory information. We present the biological rationale for supporting the proposed model. Finally, we discuss future contributions, trends, and the importance of the proposed research.

Ipsilateral and Contralateral Interactions in Spinal Locomotor Circuits Mediated by V1 Neurons: Insights from Computational Modeling

We describe and analyze a computational model of neural circuits in the mammalian spinal cord responsible for generating and shaping locomotor-like oscillations. The model represents interacting populations of spinal neurons, including the neurons that were genetically identified and characterized in a series of previous experimental studies. Here, we specifically focus on the ipsilaterally projecting V1 interneurons, their possible role in the spinal locomotor circuitry, and their involvement in the generation of locomotor oscillations. The proposed connections of these neurons and their involvement in different neuronal pathways in the spinal cord allow the model to reproduce the results of optogenetic manipulations of these neurons under different experimental conditions. We suggest the existence of two distinct populations of V1 interneurons mediating different ipsilateral and contralateral interactions within the spinal cord. The model proposes explanations for multiple experimental data concerning the effects of optogenetic silencing and activation of V1 interneurons on the frequency of locomotor oscillations in the intact cord and hemicord under different experimental conditions. Our simulations provide an important insight into the organization of locomotor circuitry in the mammalian spinal cord.

Locomotor Pattern and Force Generation Modulated by Ionic Channels: A Computational Study of Spinal Networks Underlying Locomotion

Locomotion is a fundamental movement in vertebrates produced by spinal networks known as central pattern generators (CPG). During fictive locomotion cat lumbar motoneurons (MNs) exhibit changes in membrane properties, including hyperpolarization of voltage threshold, reduction of afterhyperpolarization and input resistance, and amplification of nonlinear membrane properties. Both modeling and electrophysiological studies suggest that these changes can be produced by upregulating voltage-gated sodium channel (VGSC), persistent sodium (NaP), or L-type calcium channel (LTCC) or downregulating delayed-rectifier potassium (K(DR)) or calcium-dependent potassium channel (KCa) in spinal MNs. Further studies implicate that these channel modulations increase motor output and facilitate MN recruitment. However, it remains unknown how the channel modulation of CPG networks or MN pools affects the rhythmic generation of locomotion and force production of skeletal muscle during locomotion. In order to investigate this issue, we built a two-level CPG model composed of excitatory interneuron pools (Exc-INs), coupled reciprocally with inhibitory interneuron pools (Inh-INs), and projected to the flexor-extensor MN pools innervating skeletal muscles. Each pool consisted of 100 neurons with membrane properties based on cat spinal neurons. VGSC, K(DR), NaP, KCa, LTCC, and H-current channels were included in the model. Simulation results showed that (1) upregulating VGSC, NaP, or LTCC or downregulating KCa in MNs increased discharge rate and recruitment of MNs, thus facilitating locomotor pattern formation, increased amplitude of electroneurogram (ENG) bursting, and enhanced force generation of skeletal muscles. (2) The same channel modulation in Exc-INs increased the firing frequency of the Exc-INs, facilitated rhythmic generation, and increased flexor-extensor durations of step cycles. (3) Contrarily, downregulation of NaP or LTCC in MNs or Exc-INs or both CPG (Exc-INs and Inh-INs) and MNs disrupted locomotor pattern and reduced or even blocked the ENG bursting of MNs and force generation of skeletal muscles. (4) Pharmacological experiments showed that bath application of 25 μM nimodipine or 2 μM riluzole completely blocked fictive locomotion in isolated rat spinal cord, consistent with simulation results. We concluded that upregulation of VGSC, NaP, or LTCC or downregulation of KCa facilitated rhythmic generation and force production during walking, with NaP and LTCC playing an essential role.

Contribution of Afferent Feedback to Adaptive Hindlimb Walking in Cats: A Neuromusculoskeletal Modeling Study

Mammalian locomotion is generated by central pattern generators (CPGs) in the spinal cord, which produce alternating flexor and extensor activities controlling the locomotor movements of each limb. Afferent feedback signals from the limbs are integrated by the CPGs to provide adaptive control of locomotion. Responses of CPG-generated neural activity to afferent feedback stimulation have been previously studied during fictive locomotion in immobilized cats. Yet, locomotion in awake, behaving animals involves dynamic interactions between central neuronal circuits, afferent feedback, musculoskeletal system, and environment. To study these complex interactions, we developed a model simulating interactions between a half-center CPG and the musculoskeletal system of a cat hindlimb. Then, we analyzed the role of afferent feedback in the locomotor adaptation from a dynamic viewpoint using the methods of dynamical systems theory and nullcline analysis. Our model reproduced limb movements during regular cat walking as well as adaptive changes of these movements when the foot steps into a hole. The model generates important insights into the mechanism for adaptive locomotion resulting from dynamic interactions between the CPG-based neural circuits, the musculoskeletal system, and the environment.

Control of Mammalian Locomotion by Somatosensory Feedback

When animals walk overground, mechanical stimuli activate various receptors located in muscles, joints, and skin. Afferents from these mechanoreceptors project to neuronal networks controlling locomotion in the spinal cord and brain. The dynamic interactions between the control systems at different levels of the neuraxis ensure that locomotion adjusts to its environment and meets task demands. In this article, we describe and discuss the essential contribution of somatosensory feedback to locomotion. We start with a discussion of how biomechanical properties of the body affect somatosensory feedback. We follow with the different types of mechanoreceptors and somatosensory afferents and their activity during locomotion. We then describe central projections to locomotor networks and the modulation of somatosensory feedback during locomotion and its mechanisms. We then discuss experimental approaches and animal models used to investigate the control of locomotion by somatosensory feedback before providing an overview of the different functional roles of somatosensory feedback for locomotion. Lastly, we briefly describe the role of somatosensory feedback in the recovery of locomotion after neurological injury. We highlight the fact that somatosensory feedback is an essential component of a highly integrated system for locomotor control. © 2021 American Physiological Society. Compr Physiol 11:1-71, 2021.

Central pattern generator and human locomotion in the context of referent control of motor actions

Unperturbed human locomotion presumably results from feedforward shifts in stable body equilibrium in the environment, thus avoiding falling and subsequent catching considered in alternative theories of locomotion. Such shifts are achieved by relocation of the referent body configuration at which multiple muscle recruitment begins. Rather than being directly specified by a central pattern generator, multiple muscles are activated depending on the extent to which the body is deflected from the referent, threshold body configuration, as confirmed in previous studies. Based on the referent control theory of action and perception, solutions to classical problems in motor control are offered, including the previously unresolved problem of the integration of central and reflex influences on motoneurons and the problem of how posture and movement are related. The speed of locomotion depends on the rate of shifts in the referent body configuration. The transition from walking to running results from increasing the rate of referent shifts. It is emphasised that there is a certain hierarchy between reciprocal and co-activation of agonist and antagonist muscles during locomotion and other motor actions, which is also essential for the understanding of how locomotor speed is regulated. The analysis opens a new avenue in neurophysiological approaches to human locomotion with clinical implications.

Modeling Post-Scratching Locomotion with Two Rhythm Generators and a Shared Pattern Formation

This study aimed to present a model of post-scratching locomotion with two intermixed central pattern generator (CPG) networks, one for scratching and another for locomotion. We hypothesized that the rhythm generator layers for each CPG are different, with the condition that both CPGs share their supraspinal circuits and their motor outputs at the level of their pattern formation networks. We show that the model reproduces the post-scratching locomotion latency of 6.2 ± 3.5 s, and the mean cycle durations for scratching and post-scratching locomotion of 0.3 ± 0.09 s and 1.7 ± 0.6 s, respectively, which were observed in a previous experimental study. Our findings show how the transition of two rhythmic movements could be mediated by information exchanged between their CPG circuits through routes converging in a common pattern formation layer. This integrated organization may provide flexible and effective connectivity despite the rigidity of the anatomical connections in the spinal cord circuitry.

Common and distinct muscle synergies during level and slope locomotion in the cat

Although it is well established that the motor control system is modular, the organization of muscle synergies during locomotion and their change with ground slope are not completely understood. For example, typical reciprocal flexor-extensor muscle synergies of level walking in cats break down in downslope: one-joint hip extensors are silent throughout the stride cycle, whereas hindlimb flexors demonstrate an additional stance phase-related electromyogram (EMG) burst (Smith JL, Carlson-Kuhta P, Trank TV. J Neurophysiol 79: 1702-1716, 1998). Here, we investigated muscle synergies during level, upslope (27°), and downslope (-27°) walking in adult cats to examine common and distinct features of modular organization of locomotor EMG activity. Cluster analysis of EMG burst onset-offset times of 12 hindlimb muscles revealed five flexor and extensor burst groups that were generally shared across slopes. Stance-related bursts of flexor muscles in downslope were placed in a burst group from level and upslope walking formed by the rectus femoris. Walking upslope changed swing/stance phase durations of level walking but not the cycle duration. Five muscle synergies computed using non-negative matrix factorization accounted for at least 95% of variance in EMG patterns in each slope. Five synergies were shared between level and upslope walking, whereas only three of those were shared with downslope synergies; these synergies were active during the swing phase and phase transitions. Two stance-related synergies of downslope walking were distinct; they comprised a mixture of flexors and extensors. We suggest that the modular organization of muscle activity during level and slope walking results from interactions between motion-related sensory feedback, CPG, and supraspinal inputs.NEW & NOTEWORTHY We demonstrated that the atypical EMG activities during cat downslope walking, silent one-joint hip extensors and stance-related EMG bursts in flexors, have many features shared with activities of level and upslope walking. Majority of EMG burst groups and muscle synergies were shared among these slopes, and upslope modulated the swing/stance phase duration but not cycle duration. Thus, synergistic EMG activities in all slopes might result from a shared CPG receiving somatosensory and supraspinal inputs.

Computational Modeling of Spinal Locomotor Circuitry in the Age of Molecular Genetics

Neuronal circuits in the spinal cord are essential for the control of locomotion. They integrate supraspinal commands and afferent feedback signals to produce coordinated rhythmic muscle activations necessary for stable locomotion. For several decades, computational modeling has complemented experimental studies by providing a mechanistic rationale for experimental observations and by deriving experimentally testable predictions. This symbiotic relationship between experimental and computational approaches has resulted in numerous fundamental insights. With recent advances in molecular and genetic methods, it has become possible to manipulate specific constituent elements of the spinal circuitry and relate them to locomotor behavior. This has led to computational modeling studies investigating mechanisms at the level of genetically defined neuronal populations and their interactions. We review literature on the spinal locomotor circuitry from a computational perspective. By reviewing examples leading up to and in the age of molecular genetics, we demonstrate the importance of computational modeling and its interactions with experiments. Moving forward, neuromechanical models with neuronal circuitry modeled at the level of genetically defined neuronal populations will be required to further unravel the mechanisms by which neuronal interactions lead to locomotor behavior.

A Reciprocal Excitatory Reflex Between Extensors Reproduces the Prolongation of Stance Phase in Walking Cats: Analysis on a Robotic Platform

Spinal reflex is essential to the robust locomotion of quadruped animals. To investigate the reflex mechanisms, we developed a quadruped robot platform that emulates the neuromuscular dynamics of animals. The leg is designed to be highly back-drivable, and four Hill-type muscles and neuronal pathways are simulated on each leg using software. By searching for the reflex circuit that contributes to the generation of steady gait in cats through robotic experiments, we found a simple reflex circuit that could produce leg trajectories and a steady gait. In addition, this circuit can reproduce the experimental behavior observed in cats. As a major contribution of this study, we show that the underlying structure of the reflex circuit is the reciprocal coupling between extensor muscles via excitatory neural pathways. In the walking experiments on the robot, a steady gait and experimental behaviors of walking cats emerged from the reflex circuit without any central pattern generators. Furthermore, to take advantage of walking experiments using a neurophysiological robotic platform, we conducted experiments in which a part of the proposed reflex circuit was disconnected for a certain period of time during walking. The results showed that the prolongation of the stance phase caused by the reciprocal excitatory reflex contributed greatly to the generation of a steady gait.

Approaches to revealing the neural basis of muscle synergies: a review and a critique

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.

Neural Interactions in Developing Rhythmogenic Spinal Networks: Insights From Computational Modeling

The mechanisms involved in generation of rhythmic locomotor activity in the mammalian spinal cord remain poorly understood. These mechanisms supposedly rely on both intrinsic properties of constituting neurons and interactions between them. A subset of Shox2 neurons was suggested to contribute to generation of spinal locomotor activity, but the possible cellular basis for rhythmic bursting in these neurons remains unknown. Ha and Dougherty (2018) recently revealed the presence of bidirectional electrical coupling between Shox2 neurons in neonatal spinal cords, which can be critically involved in neuronal synchronization and generation of populational bursting. Gap junctional connections found between functionally-related Shox2 interneurons decrease with age, possibly being replaced by increasing interactions through chemical synapses. Here, we developed a computational model of a heterogeneous population of neurons sparsely connected by electrical or/and chemical synapses and investigated the dependence of frequency of populational bursting on the type and strength of neuronal interconnections. The model proposes a mechanistic explanation that can account for the emergence of a synchronized rhythmic activity in the neuronal population and provides insights into the possible role of gap junctional coupling between Shox2 neurons in the spinal mechanisms for locomotor rhythm generation.

On the Organization of the Locomotor CPG: Insights From Split-Belt Locomotion and Mathematical Modeling

Rhythmic limb movements during locomotion are controlled by central pattern generator (CPG) circuits located in the spinal cord. It is considered that these circuits are composed of individual rhythm generators (RGs) for each limb interacting with each other through multiple commissural and long propriospinal circuits. The organization and operation of each RG are not fully understood, and different competing theories exist about interactions between its flexor and extensor components, as well as about left-right commissural interactions between the RGs. The central idea of circuit organization proposed in this study is that with an increase of excitatory input to each RG (or an increase in locomotor speed) the rhythmogenic mechanism of the RGs changes from "flexor-driven" rhythmicity to a "classical half-center" mechanism. We test this hypothesis using our experimental data on changes in duration of stance and swing phases in the intact and spinal cats walking on the ground or tied-belt treadmills (symmetric conditions) or split-belt treadmills with different left and right belt speeds (asymmetric conditions). We compare these experimental data with the results of mathematical modeling, in which simulated CPG circuits operate in similar symmetric and asymmetric conditions with matching or differing control drives to the left and right RGs. The obtained results support the proposed concept of state-dependent changes in RG operation and specific commissural interactions between the RGs. The performed simulations and mathematical analysis of model operation under different conditions provide new insights into CPG network organization and limb coordination during locomotion.

Candidate Interneurons Mediating the Resetting of the Locomotor Rhythm by Extensor Group I Afferents in the Cat

Sensory information arising from limb movements controls the spinal locomotor circuitry to adapt the motor pattern to demands of the environment. Stimulation of extensor group (gr) I afferents during fictive locomotion in decerebrate cats prolongs the ongoing extension, and terminates ongoing flexion with an initiation of the subsequent extension, i. e. "resetting to extension". Moreover, instead of the classical Ib non-reciprocal inhibition, stimulation of extensor gr I afferents produces a polysynaptic excitation in extensor motoneurons with latencies (∼3.5-4.0 ms) compatible with 3 interposed interneurons. We assume that some interneurons in this pathway actually belong to the rhythm-generating layer of the locomotor Central Pattern Generator (CPG), since their activity was correlated to a resetting of the rhythm. In the present work fictive locomotion was (mostly) induced by i.v. injection of nialamide followed by l-DOPA in paralyzed cats following decerebration and spinalization at C1 level. In some experiments, we extended previous observations during fictive locomotion on the emergence and locomotor state-dependence of polysynaptic excitatory postsynaptic potentials from extensor gr I afferents to ankle extensor motoneurons. However, the main focus was to record location and properties of interneurons (n = 62) that (i) were active during the extensor phase of fictive locomotion and (ii) received short-latency excitation (mono-, di- or polysynaptic) from extensor gr I afferents. We conclude that the interneurons recorded fulfill the characteristics to belong to the neuronal pathway activated by extensor gr I afferents during locomotion, and may contribute to the 'resetting to extension' as part of the locomotor CPG.

Epidural Electrical Stimulation: A Review of Plasticity Mechanisms That Are Hypothesized to Underlie Enhanced Recovery From Spinal Cord Injury With Stimulation

Spinal cord injury (SCI) often results in life-long sensorimotor impairment. Spontaneous recovery from SCI is limited, as supraspinal fibers cannot spontaneously regenerate to form functional networks below the level of injury. Despite this, animal models and humans exhibit many motor behaviors indicative of recovery when electrical stimulation is applied epidurally to the dorsal aspect of the lumbar spinal cord. In 1976, epidural stimulation was introduced to alleviate spasticity in Multiple Sclerosis. Since then, epidural electrical stimulation (EES) has been demonstrated to improve voluntary mobility across the knee and/or ankle in several SCI patients, highlighting its utility in enhancing motor activation. The mechanisms that EES induces to drive these improvements in sensorimotor function remain largely unknown. In this review, we discuss several sensorimotor plasticity mechanisms that we hypothesize may enable epidural stimulation to promote recovery, including changes in local lumbar circuitry, propriospinal interneurons, and the internal model. Finally, we discuss genetic tools for afferent modulation as an emerging method to facilitate the search for the mechanisms of action.

Ambulation in Dogs With Absent Pain Perception After Acute Thoracolumbar Spinal Cord Injury

Acute thoracolumbar spinal cord injury (SCI) is common in dogs frequently secondary to intervertebral disc herniation. Following severe injury, some dogs never regain sensory function to the pelvic limbs or tail and are designated chronically "deep pain negative." Despite this, a subset of these dogs develop spontaneous motor recovery over time including some that recover sufficient function in their pelvic limbs to walk independently without assistance or weight support. This type of ambulation is commonly known as "spinal walking" and can take up to a year or more to develop. This review provides a comparative overview of locomotion and explores the physiology of locomotor recovery after severe SCI in dogs. We discuss the mechanisms by which post-injury plasticity and coordination between circuitry contained within the spinal cord, peripheral sensory feedback, and residual or recovered supraspinal connections might combine to underpin spinal walking. The clinical characteristics of spinal walking are outlined including what is known about the role of patient or injury features such as lesion location, timeframe post-injury, body size, and spasticity. The relationship between the emergence of spinal walking and electrodiagnostic and magnetic resonance imaging findings are also discussed. Finally, we review possible ways to predict or facilitate recovery of walking in chronically deep pain negative dogs. Improved understanding of the mechanisms of gait generation and plasticity of the surviving tissue after injury might pave the way for further treatment options and enhanced outcomes in severely injured dogs.

The Functional Role of Spinal Interneurons Following Traumatic Spinal Cord Injury

Traumatic spinal cord injury (SCI) impedes signal transmission by disrupting both the local neurons and their surrounding synaptic connections. Although the majority of SCI patients retain spared neural tissue at the injury site, they predominantly suffer from complete autonomic and sensorimotor dysfunction. While there have been significant advances in the characterization of the spared neural tissue following SCI, the functional role of injury-induced interneuronal plasticity remains elusive. In healthy individuals, spinal interneurons are responsible for relaying signals to coordinate both sympathetic and parasympathetic functions. However, the spontaneous synaptic loss following injury alters these intricate interneuronal networks in the spinal cord. Here, we propose the synaptopathy hypothesis of SCI based on recent findings regarding the maladaptive role of synaptic changes amongst the interneurons. These maladaptive consequences include circuit inactivation, neuropathic pain, spasticity, and autonomic dysreflexia. Recent preclinical advances have uncovered the therapeutic potential of spinal interneurons in activating the dormant relay circuits to restore sensorimotor function. This review will survey the diverse role of spinal interneurons in SCI pathogenesis as well as treatment strategies to target spinal interneurons.

Fictive Scratching Patterns in Brain Cortex-Ablated, Midcollicular Decerebrate, and Spinal Cats

Background: The spinal cord’s central pattern generators (CPGs) have been explained by the symmetrical half-center hypothesis, the bursts generator, computational models, and more recently by connectome circuits. Asymmetrical models, at odds with the half-center paradigm, are composed of extensor and flexor CPG modules. Other models include not only flexor and extensor motoneurons but also motoneuron pools controlling biarticular muscles. It is unknown whether a preferred model can explain some particularities that fictive scratching (FS) in the cat presents. The first aim of this study was to investigate FS patterns considering the aiming and the rhythmic periods, and second, to examine the effects of serotonin (5HT) on and segmental inputs to FS.

Methods: The experiments were carried out first in brain cortex-ablated cats (BCAC), then spinalized (SC), and for the midcollicular (MCC) preparation. Subjects were immobilized and the peripheral nerves were used to elicit the Monosynaptic reflex (MR), to modify the scratching patterns and for electroneurogram recordings.

Results: In BCAC, FS was produced by pinna stimulation and, in some cases, by serotonin. The scratching aiming phase (AP) initiates with the activation of either flexor or extensor motoneurons. Serotonin application during the AP produced simultaneous extensor and flexor bursts. Furthermore, WAY 100635 (5HT1A antagonist) produced a brief burst in the tibialis anterior (TA) nerve, followed by a reduction in its electroneurogram (ENG), while the soleus ENG remained silent. In SC, rhythmic phase (RP) activity was recorded in the soleus motoneurons. Serotonin or WAY produced FS bouts. The electrical stimulation of Ia afferent fibers produced heteronymous MRes waxing and waning during the scratch cycle. In MCC, FS began with flexor activity. Electrical stimulation of either deep peroneus (DP) or superficial peroneus (SP) nerves increased the duration of the TA electroneurogram. Medial gastrocnemius (MG) stretching or MG nerve electrical stimulation produced a reduction in the TA electroneurogram and an initial MG extensor burst. MRes waxed and waned during the scratch cycle.

Conclusion: Descending pathways and segmental afferent fibers, as well as 5-HT and WAY, can change the FS pattern. To our understanding, the half-center hypothesis is the most suitable for explaining the AP in MCC.

Contribution of Phase Resetting to Adaptive Rhythm Control in Human Walking Based on the Phase Response Curves of a Neuromusculoskeletal Model

Humans walk adaptively in varying environments by manipulating their complicated and redundant musculoskeletal system. Although the central pattern generators in the spinal cord are largely responsible for adaptive walking through sensory-motor coordination, it remains unclear what neural mechanisms determine walking adaptability. It has been reported that locomotor rhythm and phase are regulated by the production of phase shift and rhythm resetting (phase resetting) for periodic motor commands in response to sensory feedback and perturbation. While the phase resetting has been suggested to make a large contribution to adaptive walking, it has only been investigated based on fictive locomotion in decerebrate cats, and thus it remains unclear if human motor control has such a rhythm regulation mechanism during walking. In our previous work, we incorporated a phase resetting mechanism into a motor control model and demonstrated that it improves the stability and robustness of walking through forward dynamic simulations of a human musculoskeletal model. However, this did not necessarily verify that phase resetting plays a role in human motor control. In our other previous work, we used kinematic measurements of human walking to identify the phase response curve (PRC), which explains phase-dependent responses of a limit cycle oscillator to a perturbation. This revealed how human walking rhythm is regulated by perturbations. In this study, we integrated these two approaches using a physical model and identification of the PRC to examine the hypothesis that phase resetting plays a role in the control of walking rhythm in humans. More specifically, we calculated the PRC using our neuromusculoskeletal model in the same way as our previous human experiment. In particular, we compared the PRCs calculated from two different models with and without phase resetting while referring to the PRC for humans. As a result, although the PRC for the model without phase resetting did not show any characteristic shape, the PRC for the model with phase resetting showed a characteristic phase-dependent shape with trends similar to those of the PRC for humans. These results support our hypothesis and will improve our understanding of adaptive rhythm control in human walking.

Phase-Dependent Response to Afferent Stimulation During Fictive Locomotion: A Computational Modeling Study

Central pattern generators (CPGs) in the spinal cord generate rhythmic neural activity and control locomotion in vertebrates. These CPGs operate under the control of sensory feedback that affects the generated locomotor pattern and adapt it to the animal's biomechanics and environment. Studies of the effects of afferent stimulation on fictive locomotion in immobilized cats have shown that brief stimulation of peripheral nerves can reset the ongoing locomotor rhythm. Depending on the phase of stimulation and the stimulated nerve, the applied stimulation can either shorten or prolong the current locomotor phase and the locomotor cycle. Here, we used a mathematical model of a half-center CPG to investigate the phase-dependent effects of brief stimulation applied to CPG on the CPG-generated locomotor oscillations. The CPG in the model consisted of two half-centers mutually inhibiting each other. The rhythmic activity in each half-center was based on a slowly inactivating, persistent sodium current. Brief stimulation was applied to CPG half-centers in different phases of the locomotor cycle to produce phase-dependent changes in CPG activity. The model reproduced several results from experiments on the effect of afferent stimulation of fictive locomotion in cats. The mechanisms of locomotor rhythm resetting under different conditions were analyzed using dynamic systems theory methods.

Spinal V3 Interneurons and Left-Right Coordination in Mammalian Locomotion

Commissural interneurons (CINs) mediate interactions between rhythm-generating locomotor circuits located on each side of the spinal cord and are necessary for left-right limb coordination during locomotion. While glutamatergic V3 CINs have been implicated in left-right coordination, their functional connectivity remains elusive. Here, we addressed this issue by combining experimental and modeling approaches. We employed Sim1^Cre/+; Ai32 mice, in which light-activated Channelrhodopsin-2 was selectively expressed in V3 interneurons. Fictive locomotor activity was evoked by NMDA and 5-HT in the isolated neonatal lumbar spinal cord. Flexor and extensor activities were recorded from left and right L2 and L5 ventral roots, respectively. Bilateral photoactivation of V3 interneurons increased the duration of extensor bursts resulting in a slowed down on-going rhythm. At high light intensities, extensor activity could become sustained. When light stimulation was shifted toward one side of the cord, the duration of extensor bursts still increased on both sides, but these changes were more pronounced on the contralateral side than on the ipsilateral side. Additional bursts appeared on the ipsilateral side not seen on the contralateral side. Further increase of the stimulation could suppress the contralateral oscillations by switching to a sustained extensor activity, while the ipsilateral rhythmic activity remained. To delineate the function of V3 interneurons and their connectivity, we developed a computational model of the spinal circuits consisting of two (left and right) rhythm generators (RGs) interacting via V0_V, V0_D, and V3 CINs. Both types of V0 CINs provided mutual inhibition between the left and right flexor RG centers and promoted left-right alternation. V3 CINs mediated mutual excitation between the left and right extensor RG centers. These interactions allowed the model to reproduce our current experimental data, while being consistent with previous data concerning the role of V0_V and V0_D CINs in securing left–right alternation and the changes in left–right coordination following their selective removal. We suggest that V3 CINs provide mutual excitation between the spinal neurons involved in the control of left and right extensor activity, which may promote left-right synchronization during locomotion.

Variant and Invariant Spatiotemporal Structures in Kinematic Coordination to Regulate Speed During Walking and Running

Humans walk, run, and change their speed in accordance with circumstances. These gaits are rhythmic motions generated by multi-articulated movements, which have specific spatiotemporal patterns. The kinematic characteristics depend on the gait and speed. In this study, we focused on the kinematic coordination of locomotor behavior to clarify the underlying mechanism for the effect of speed on the spatiotemporal kinematic patterns for each gait. In particular, we used seven elevation angles for the whole-body motion and separated the measured data into different phases depending on the foot-contact condition, that is, single-stance phase, double-stance phase, and flight phase, which have different physical constraints during locomotion. We extracted the spatiotemporal kinematic coordination patterns with singular value decomposition and investigated the effect of speed on the coordination patterns. Our results showed that most of the whole-body motion could be explained by only two sets of temporal and spatial coordination patterns in each phase. Furthermore, the temporal coordination patterns were invariant for different speeds, while the spatial coordination patterns varied. These findings will improve our understanding of human adaptation mechanisms to tune locomotor behavior for changing speed.

Flexor and Extensor Ankle Afferents Broadly Innervate Locomotor Spinal Shox2 Neurons and Induce Similar Effects in Neonatal Mice

Central pattern generators (CPGs) in the thoracolumbar spinal cord generate the basic hindlimb locomotor pattern. The locomotor CPG integrates descending commands and sensory information from the periphery to activate, modulate and halt the rhythmic program. General CPG function and response to sensory perturbations are well described in cat and rat models. In mouse, roles for many genetically identified spinal interneurons have been inferred from locomotor alterations following population deletion or modulation. However, the organization of afferent input to specific genetically identified populations of spinal CPG interneurons in mouse remains comparatively less resolved. Here, we focused on a population of CPG neurons marked by the transcription factor Shox2. To directly test integration of afferent signaling by Shox2 neurons, sensory afferents were stimulated during patch clamp recordings of Shox2 neurons in isolated spinal cord preparations from neonatal mice. Shox2 neurons broadly displayed afferent-evoked currents at multiple segmental levels, particularly from caudal dorsal roots innervating distal hindlimb joints. As dorsal root stimulation may activate both flexor- and extensor-related afferents, preparations preserving peripheral nerves were used to provide more specific activation of ankle afferents. We found that both flexor- and extensor-related afferent stimulation were likely to evoke similar currents in a given Shox2 neuron, as assessed by response polarity, latency, duration and amplitude. It has been proposed that Shox2 neurons can be divided into neurons which contribute to rhythm generation and neurons that are premotor by the absence and presence of the V2a marker Chx10, respectively. Response to afferent stimulation did not differ based on Chx10 expression. Although currents evoked in response to flexor and extensor afferent activation did not follow expected functional antagonism, they were consistent with the observation that stimulation of flexor- and extensor-related afferents both reset the phase of ongoing fictive locomotion to flexion in neonatal mice. Together, the data suggest that Shox2 neurons are interposed in multiple sensory pathways and low threshold proprioceptive input reinforces sensory perturbation of ongoing locomotion by similarly activating or inhibiting both the rhythm and patterning layers of the CPG.

Curated Model Development Using NEUROiD: A Web-Based NEUROmotor Integration and Design Platform

Decades of research on neuromotor circuits and systems has provided valuable information on neuronal control of movement. Computational models of several elements of the neuromotor system have been developed at various scales, from sub-cellular to system. While several small models abound, their structured integration is the key to building larger and more biologically realistic models which can predict the behavior of the system in different scenarios. This effort calls for integration of elements across neuroscience and musculoskeletal biomechanics. There is also a need for development of methods and tools for structured integration that yield larger in silico models demonstrating a set of desired system responses. We take a small step in this direction with the NEUROmotor integration and Design (NEUROiD) platform. NEUROiD helps integrate results from motor systems anatomy, physiology, and biomechanics into an integrated neuromotor system model. Simulation and visualization of the model across multiple scales is supported. Standard electrophysiological operations such as slicing, current injection, recording of membrane potential, and local field potential are part of NEUROiD. The platform allows traceability of model parameters to primary literature. We illustrate the power and utility of NEUROiD by building a simple ankle model and its controlling neural circuitry by curating a set of published components. NEUROiD allows researchers to utilize remote high-performance computers for simulation, while controlling the model using a web browser.

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.

Sub-populations of Spinal V3 Interneurons Form Focal Modules of Layered Pre-motor Microcircuits

Layering of neural circuits facilitates the separation of neurons with high spatial sensitivity from those that play integrative temporal roles. Although anatomical layers are readily identifiable in the brain, layering is not structurally obvious in the spinal cord. But computational studies of motor behaviors have led to the concept of layered processing in the spinal cord. It has been postulated that spinal V3 interneurons (INs) play multiple roles in locomotion, leading us to investigate whether they form layered microcircuits. Using patch-clamp recordings in combination with holographic glutamate uncaging, we demonstrate focal, layered modules, in which ventromedial V3 INs form synapses with one another and with ventrolateral V3 INs, which in turn form synapses with ipsilateral motoneurons. Motoneurons, in turn, provide recurrent excitatory, glutamatergic input to V3 INs. Thus, ventral V3 interneurons form layered microcircuits that could function to ensure well-timed, spatially specific movements.

Walking with perturbations: a guide for biped humans and robots

This paper provides an update on the neural control of bipedal walking in relation to bioinspired models and robots. It is argued that most current models or robots are based on the construct of a symmetrical central pattern generator (CPG). However, new evidence suggests that CPG functioning is basically asymmetrical with its flexor half linked more tightly to the rhythm generator. The stability of bipedal gait, which is an important problem for robots and biological systems, is also addressed. While it is not possible to determine how biological biped systems guarantee stability, robot solutions can be useful to propose new hypotheses for biology. In the second part of this review, the focus is on gait perturbations, which is an important topic in robotics in view of the frequent falls of robots when faced with perturbations. From the human physiology it is known that the initial reaction often consists of a brief interruption followed by an adequate response. For instance, the successful recovery from a trip is achieved using some basic reactions (termed elevating and lowering strategies), that depend on the phase of the step cycle of the trip occurrence. Reactions to stepping unexpectedly in a hole depend on comparing expected and real feedback. Implementation of these ideas in models and robotics starts to emerge, with the most advanced robots being able to learn how to fall safely and how to deal with complicated disturbances such as provided by walking on a split-belt.

Functional Recovery of a Locomotor Network after Injury: Plasticity beyond the Central Nervous System

Many animals depend on descending information from the brain for the initiation and proper execution of locomotion. Interestingly, after injury and the loss of such inputs, locomotor function can sometimes be regained without the regrowth of central connections. In the medicinal leech, Hirudo verbana, we have shown that crawling reemerges after removal of descending inputs. Here, we studied the mechanisms underlying this return of locomotion by asking if central pattern generators (CPGs) in crawl-recovered leeches are sufficient to produce crawl-specific intersegmental coordination. From recovered animals, we treated isolated chains of ganglia with dopamine to activate the crawl CPGs (one crawl CPG per ganglion) and observed fictive crawl-like bursting in the dorsal-longitudinal-excitor motoneuron (DE-3), an established crawl-monitor neuron. However, these preparations did not exhibit crawl-specific coordination across the CPGs. Although the crawl CPGs always generated bidirectional activation of adjacent CPGs, we never observed crawl-appropriate intersegmental phase delays. Because central circuits alone were unable to organize crawl-specific coordination, we tested the coordinating role of the peripheral nervous system. In transected leeches normally destined for recovery, we removed afferent information to the anterior-most (lead) ganglion located below the nerve-cord transection site. In these dually treated animals, overt crawling was greatly delayed or prevented. After filling the peripheral nerves with Neurobiotin tracer distal to the nerve-root lesion, we found a perfect correlation between regrowth of peripheral neuronal fibers and crawl recovery. Our study establishes that during recovery after injury, crawl-specific intersegmental coordination switches to a new dependence on afferent information.

State-dependent rhythmogenesis and frequency control in a half-center locomotor CPG

The spinal locomotor central pattern generator (CPG) generates rhythmic activity with alternating flexion and extension phases. This rhythmic pattern is likely to result from inhibitory interactions between neural populations representing flexor and extensor half-centers. However, it is unclear whether the flexor-extensor CPG has a quasi-symmetric organization with both half-centers critically involved in rhythm generation, features an asymmetric organization with flexor-driven rhythmogenesis, or comprises a pair of intrinsically rhythmic half-centers. There are experimental data that support each of the above concepts but appear to be inconsistent with the others. In this theoretical/modeling study, we present and analyze a CPG model architecture that can operate in different regimes consistent with the above three concepts depending on conditions, which are defined by external excitatory drives to CPG half-centers. We show that control of frequency and phase durations within each regime depends on network dynamics, defined by the regime-dependent expression of the half-centers' intrinsic rhythmic capabilities and the operating phase transition mechanisms (escape vs. release). Our study suggests state dependency in locomotor CPG operation and proposes explanations for seemingly contradictory experimental data. NEW & NOTEWORTHY Our theoretical/modeling study focuses on the analysis of locomotor central pattern generators (CPGs) composed of conditionally bursting half-centers coupled with reciprocal inhibition and receiving independent external drives. We show that this CPG framework can operate in several regimes consistent with seemingly contradictory experimental data. In each regime, we study how intrinsic dynamics and phase-switching mechanisms control oscillation frequency and phase durations. Our results provide insights into the organization of spinal circuits controlling locomotion.

Computational modeling of spinal circuits controlling limb coordination and gaits in quadrupeds

Interactions between cervical and lumbar spinal circuits are mediated by long propriospinal neurons (LPNs). Ablation of descending LPNs in mice disturbs left-right coordination at high speeds without affecting fore-hind alternation. We developed a computational model of spinal circuits consisting of four rhythm generators coupled by commissural interneurons (CINs), providing left-right interactions, and LPNs, mediating homolateral and diagonal interactions. The proposed CIN and diagonal LPN connections contribute to speed-dependent gait transition from walk, to trot, and then to gallop and bound; the homolateral LPN connections ensure fore-hind alternation in all gaits. The model reproduces speed-dependent gait expression in intact and genetically transformed mice and the disruption of hindlimb coordination following ablation of descending LPNs. Inputs to CINs and LPNs can affect interlimb coordination and change gait independent of speed. We suggest that these interneurons represent the main targets for supraspinal and sensory afferent signals adjusting gait.

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.

Central control of interlimb coordination and speed-dependent gait expression in quadrupeds

KEY POINTS

Quadrupeds express different gaits depending on speed of locomotion. Central pattern generators (one per limb) within the spinal cord generate locomotor oscillations and control limb movements. Neural interactions between these generators define interlimb coordination and gait. We present a computational model of spinal circuits representing four rhythm generators with left-right excitatory and inhibitory commissural and fore-hind inhibitory interactions within the cord. Increasing brainstem drive to all rhythm generators and excitatory commissural interneurons induces an increasing frequency of locomotor oscillations accompanied by speed-dependent gait changes from walk to trot and to gallop and bound. The model closely reproduces and suggests explanations for multiple experimental data, including speed-dependent gait transitions in intact mice and changes in gait expression in mutants lacking certain types of commissural interneurons. The model suggests the possible circuit organization in the spinal cord and proposes predictions that can be tested experimentally.

ABSTRACT

As speed of locomotion is increasing, most quadrupeds, including mice, demonstrate sequential gait transitions from walk to trot and to gallop and bound. The neural mechanisms underlying these transitions are poorly understood. We propose that the speed-dependent expression of different gaits results from speed-dependent changes in the interactions between spinal circuits controlling different limbs and interlimb coordination. As a result, the expression of each gait depends on (1) left-right interactions within the spinal cord mediated by different commissural interneurons (CINs), (2) fore-hind interactions on each side of the spinal cord and (3) brainstem drives to rhythm-generating circuits and CIN pathways. We developed a computational model of spinal circuits consisting of four rhythm generators (RGs) with bilateral left-right interactions mediated by V0 CINs (V0_D and V0_V sub-types) providing left-right alternation, and conditional V3 CINs promoting left-right synchronization. Fore and hind RGs mutually inhibited each other. We demonstrate that linearly increasing excitatory drives to the RGs and V3 CINs can produce a progressive increase in the locomotor speed accompanied by sequential changes of gaits from walk to trot and to gallop and bound. The model closely reproduces and suggests explanations for the speed-dependent gait expression observed in vivo in intact mice and in mutants lacking V0_V or all V0 CINs. Specifically, trot is not expressed after removal of V0_V CINs, and only bound is expressed after removal of all V0 CINs. The model provides important insights into the organization of spinal circuits and neural control of locomotion.

Intraspinal microstimulation produces over-ground walking in anesthetized cats

OBJECTIVE

Spinal cord injury causes a drastic loss of motor, sensory and autonomic function. The goal of this project was to investigate the use of intraspinal microstimulation (ISMS) for producing long distances of walking over ground. ISMS is an electrical stimulation method developed for restoring motor function by activating spinal networks below the level of an injury. It produces movements of the legs by stimulating the ventral horn of the lumbar enlargement using fine penetrating electrodes (≤50 μm diameter).

APPROACH

In each of five adult cats (4.2-5.5 kg), ISMS was applied through 16 electrodes implanted with tips targeting lamina IX in the ventral horn bilaterally. A desktop system implemented a physiologically-based control strategy that delivered different stimulation patterns through groups of electrodes to evoke walking movements with appropriate limb kinematics and forces corresponding to swing and stance. Each cat walked over an instrumented 2.9 m walkway and limb kinematics and forces were recorded.

MAIN RESULTS

Both propulsive and supportive forces were required for over-ground walking. Cumulative walking distances ranging from 609 to 835 m (longest tested) were achieved in three animals. In these three cats, the mean peak supportive force was 3.5 ± 0.6 N corresponding to full-weight-support of the hind legs, while the angular range of the hip, knee, and ankle joints were 23.1 ± 2.0°, 29.1 ± 0.2°, and 60.3 ± 5.2°, respectively. To further demonstrate the viability of ISMS for future clinical use, a prototype implantable module was successfully implemented in a subset of trials and produced comparable walking performance.

SIGNIFICANCE

By activating inherent locomotor networks within the lumbosacral spinal cord, ISMS was capable of producing bilaterally coordinated and functional over-ground walking with current amplitudes <100 μA. These exciting results suggest that ISMS may be an effective intervention for restoring functional walking after spinal cord injury.

Organization of flexor-extensor interactions in the mammalian spinal cord: insights from computational modelling

KEY POINTS

Alternation of flexor and extensor activity in the mammalian spinal cord is mediated by two classes of genetically identified inhibitory interneurons, V1 and V2b. The V1 interneurons are essential for high-speed locomotor activity. They secure flexor-extensor alternations in the intact cord but lose this function after hemisection, suggesting that they are activated by inputs from the contralateral side of the cord. The V2b interneurons are involved in flexor-extensor alternations in both intact cord and hemicords. We used a computational model of the spinal network, simulating the left and right rhythm-generating circuits interacting via several commissural pathways, and extended this model by incorporating V1 and V2b neuron populations involved in flexor-extensor interactions on each cord side. The model reproduces multiple experimental data on selective silencing and activation of V1 and/or V2b neurons and proposes the organization of their connectivity providing flexor-extensor alternation in the spinal cord.

ABSTRACT

Alternating flexor and extensor activity represents the fundamental property underlying many motor behaviours including locomotion. During locomotion this alternation appears to arise in rhythm-generating circuits and transpires at all levels of the spinal cord including motoneurons. Recent studies in vitro and in vivo have shown that flexor-extensor alternation during locomotion involves two classes of genetically identified, inhibitory interneurons: V1 and V2b. Particularly, in the isolated mouse spinal cord, abrogation of neurotransmission derived by both V1 and V2b interneurons resulted in flexor-extensor synchronization, whereas selective inactivation of only one of these neuron types did not abolish flexor-extensor alternation. After hemisection, inactivation of only V2b interneurons led to the flexor-extensor synchronization, while inactivation of V1 interneurons did not affect flexor-extensor alternation. Moreover, optogenetic activation of V2b interneurons suppressed extensor-related activity, while similar activation of V1 interneurons suppressed both flexor and extensor oscillations. Here, we address these issues using the previously published computational model of spinal circuitry simulating bilateral interactions between left and right rhythm-generating circuits. In the present study, we incorporate V1 and V2b neuron populations on both sides of the cord to make them critically involved in flexor-extensor interactions. The model reproduces multiple experimental data on the effects of hemisection and selective silencing or activation of V1 and V2b neurons and suggests connectivity profiles of these neurons and their specific roles in left-right (V1) and flexor-extensor (both V2b and V1) interactions in the spinal cord that can be tested experimentally.

Quantitative analysis of hindlimbs locomotion kinematics in spinalized rats treated with Tamoxifen plus treadmill exercise

Locomotion recovery after a spinal cord injury (SCI) includes axon regeneration, myelin preservation and increased plasticity in propriospinal and descending spinal circuitries. The combined effects of tamoxifen and exercise after a SCI were analyzed in this study to determine whether the combination of both treatments induces the best outcome in locomotion recovery. In this study, the penetrating injury was provoked by a sharp projectile that penetrates through right dorsal and ventral portions of the T13-L1 spinal segments, affecting propriospinal and descending/ascending tracts. Intraperitoneal application of Tamoxifen and a treadmill exercise protocol, as rehabilitation therapies, separately or combined, were used. To evaluate the functional recovery, angular patterns of the hip, knee and ankle joints as well as the leg pendulum-like movement (PLM) were measured during the unrestricted gait of treated and untreated (UT) animals, previously and after the traumatic injury (15 and 30days post-injury (dpi)). A pattern (curve) comparison analysis was made by using a locally designed Matlab script that determines the Frechet dissimilarity. The SCI magnitude was assessed by qualitative and quantitative histological analysis of the injury site 30days after SCI. Our results showed that all treated groups had an improvement in hindlimbs kinematics compared to the UT group, which showed a poor gait locomotion recovery throughout the rehabilitation period. The group with the combined treatment (tamoxifen+exercise (TE)) presented the best outcome. In conclusion, tamoxifen and treadmill exercise treatments are complementary therapies for the functional recovery of gait locomotion in hemi-spinalized rats.

Interactions between Dorsal and Ventral Root Stimulation on the Generation of Locomotor-Like Activity in the Neonatal Mouse Spinal Cord

We investigated whether dorsal (DR) and ventral root (VR) stimulus trains engage common postsynaptic components to activate the central pattern generator (CPG) for locomotion in the neonatal mouse spinal cord. VR stimulation did not activate the first order interneurons mediating the activation of the locomotor CPG by sacrocaudal afferent stimulation. Simultaneous stimulation of adjacent dorsal or ventral root pairs, subthreshold for evoking locomotor-like activity, did not summate to activate the CPG. This suggests that locomotor-like activity is triggered when a critical class of efferent or afferent axons is stimulated and does not depend on the number of stimulated axons or activated postsynaptic neurons. DR- and VR-evoked episodes exhibited differences in the coupling between VR pairs. In DR-evoked episodes, the coupling between the ipsilateral and contralateral flexor/extensor roots was similar and stronger than the bilateral extensor roots. In VR-evoked episodes, ipsilateral flexor/extensor coupling was stronger than both the contralateral flexor/extensor and the bilateral extensor coupling. For both types of stimulation, the coupling was greatest between the bilateral L1/L2 flexor-dominated roots. This indicates that the recruitment and/or the firing pattern of motoneurons differed in DR and VR-evoked episodes. However, the DR and VR trains do not appear to activate distinct CPGs because trains of DR and VR stimuli at frequencies too low to evoke locomotor-like activity did so when they were interleaved. These results indicate that the excitatory actions of VR stimulation converge onto the CPG through an unknown pathway that is not captured by current models of the locomotor CPG.

Mixed-mode oscillations and population bursting in the pre-Bötzinger complex

This study focuses on computational and theoretical investigations of neuronal activity arising in the pre-Bötzinger complex (pre-BötC), a medullary region generating the inspiratory phase of breathing in mammals. A progressive increase of neuronal excitability in medullary slices containing the pre-BötC produces mixed-mode oscillations (MMOs) characterized by large amplitude population bursts alternating with a series of small amplitude bursts. Using two different computational models, we demonstrate that MMOs emerge within a heterogeneous excitatory neural network because of progressive neuronal recruitment and synchronization. The MMO pattern depends on the distributed neuronal excitability, the density and weights of network interconnections, and the cellular properties underlying endogenous bursting. Critically, the latter should provide a reduction of spiking frequency within neuronal bursts with increasing burst frequency and a dependence of the after-burst recovery period on burst amplitude. Our study highlights a novel mechanism by which heterogeneity naturally leads to complex dynamics in rhythmic neuronal populations.

DOI:http://dx.doi.org/10.7554/eLife.13403.001

Each breath we take removes carbon dioxide from the body and exchanges it for oxygen. A structure called the brainstem, which connects the brain with the spinal cord, generates the breathing rhythm and controls its rate. While this process normally occurs automatically, we can also control our breathing voluntarily, such as when singing or speaking.

Within the brainstem, a group of neurons in the area known as the pre-Bötzinger complex is responsible for ensuring that an animal breathes in at regular intervals. Recordings of the electrical activity from slices of brainstem show that pre-Bötzinger neurons display rhythmic activity with characteristic patterns called “mixed-mode oscillations”. These rhythms consist of bursts of strong activity (“large amplitude bursts”), essential for triggering regular breathing, separated by a series of bursts of weak activity (“small amplitude bursts”). However, it is not clear how mixed-mode oscillations arise.

Bacak, Kim et al. now provide insights into this process by developing two computational models of the pre-Bötzinger complex. The first model consists of a diverse population of 100 neurons joined by a relatively small number of weak connections to form a network. The second model is a simplified version of the first, consisting of just three neurons. By manipulating the properties of the simulated networks, and analysing the data mathematically, Bacak, Kim et al. identify the properties of the neurons that allow them to generate mixed-mode oscillations and thus rhythmic breathing.

The models suggest that mixed-mode oscillations result from the synchronization of many neurons with different levels of activity (excitability). Neurons with low excitability have low bursting frequencies, but generate strong activity and recruit other neurons, ultimately producing large amplitude bursts that cause breathing.

Many parts of the nervous system are also made up of networks of neurons with diverse excitability. A challenge for future studies is thus to investigate whether other networks of neurons similar to the pre-Bötzinger complex generate rhythms that control other repetitive actions, such as walking and chewing.

DOI:http://dx.doi.org/10.7554/eLife.13403.002

Organization of left-right coordination of neuronal activity in the mammalian spinal cord: Insights from computational modelling

KEY POINTS

Coordination of neuronal activity between left and right sides of the mammalian spinal cord is provided by several sets of commissural interneurons (CINs) whose axons cross the midline. Genetically identified inhibitory V0D and excitatory V0V CINs and ipsilaterally projecting excitatory V2a interneurons were shown to secure left-right alternation at different locomotor speeds. We have developed computational models of neuronal circuits in the spinal cord that include left and right rhythm-generating centres interacting bilaterally via three parallel pathways mediated by V0D , V2a-V0V and V3 neuron populations. The models reproduce the experimentally observed speed-dependent left-right coordination in normal mice and the changes in coordination seen in mutants lacking specific neuron classes. The models propose an explanation for several experimental results and provide insights into the organization of the spinal locomotor network and parallel CIN pathways involved in gait control at different locomotor speeds.

ABSTRACT

Different locomotor gaits in mammals, such as walking or galloping, are produced by coordinated activity in neuronal circuits in the spinal cord. Coordination of neuronal activity between left and right sides of the cord is provided by commissural interneurons (CINs), whose axons cross the midline. In this study, we construct and analyse two computational models of spinal locomotor circuits consisting of left and right rhythm generators interacting bilaterally via several neuronal pathways mediated by different CINs. The CIN populations incorporated in the models include the genetically identified inhibitory (V0D ) and excitatory (V0V ) subtypes of V0 CINs and excitatory V3 CINs. The model also includes the ipsilaterally projecting excitatory V2a interneurons mediating excitatory drive to the V0V CINs. The proposed network architectures and CIN connectivity allow the models to closely reproduce and suggest mechanistic explanations for several experimental observations. These phenomena include: different speed-dependent contributions of V0D and V0V CINs and V2a interneurons to left-right alternation of neural activity, switching gaits between the left-right alternating walking-like activity and the left-right synchronous hopping-like pattern in mutants lacking specific neuron classes, and speed-dependent asymmetric changes of flexor and extensor phase durations. The models provide insights into the architecture of spinal network and the organization of parallel inhibitory and excitatory CIN pathways and suggest explanations for how these pathways maintain alternating and synchronous gaits at different locomotor speeds. The models propose testable predictions about the neural organization and operation of mammalian locomotor circuits.

Spinal Cord Stimulation and Augmentative Control Strategies for Leg Movement after Spinal Paralysis in Humans

Severe spinal cord injury is a devastating condition, tearing apart long white matter tracts and causing paralysis and disability of body functions below the lesion. But caudal to most injuries, the majority of neurons forming the distributed propriospinal system, the localized gray matter spinal interneuronal circuitry, and spinal motoneuron populations are spared. Epidural spinal cord stimulation can gain access to this neural circuitry. This review focuses on the capability of the human lumbar spinal cord to generate stereotyped motor output underlying standing and stepping, as well as full weight-bearing standing and rhythmic muscle activation during assisted treadmill stepping in paralyzed individuals in response to spinal cord stimulation. By enhancing the excitability state of the spinal circuitry, the stimulation can have an enabling effect upon otherwise "silent" translesional volitional motor control. Strategies for achieving functional movement in patients with severe injuries based on minimal translesional intentional control, task-specific proprioceptive feedback, and next-generation spinal cord stimulation systems will be reviewed. The role of spinal cord stimulation can go well beyond the immediate generation of motor output. With recently developed training paradigms, it can become a major rehabilitation approach in spinal cord injury for augmenting and steering trans- and sublesional plasticity for lasting therapeutic benefits.

Mechanisms Underlying the Neuromodulation of Spinal Circuits for Correcting Gait and Balance Deficits after Spinal Cord Injury

Epidural electrical stimulation of lumbar segments facilitates standing and walking in animal models and humans with spinal cord injury. However, the mechanisms through which this neuromodulation therapy engages spinal circuits remain enigmatic. Using computer simulations and behavioral experiments, we provide evidence that epidural electrical stimulation interacts with muscle spindle feedback circuits to modulate muscle activity during locomotion. Hypothesis-driven strategies emerging from simulations steered the design of stimulation protocols that adjust bilateral hindlimb kinematics throughout gait execution. These stimulation strategies corrected subject-specific gait and balance deficits in rats with incomplete and complete spinal cord injury. The conservation of muscle spindle feedback circuits across mammals suggests that the same mechanisms may facilitate motor control in humans. These results provide a conceptual framework to improve stimulation protocols for clinical applications.

Neuromusculoskeletal models based on the muscle synergy hypothesis for the investigation of adaptive motor control in locomotion via sensory-motor coordination

Humans and animals walk adaptively in diverse situations by skillfully manipulating their complicated and redundant musculoskeletal systems. From an analysis of measured electromyographic (EMG) data, it appears that despite complicated spatiotemporal properties, muscle activation patterns can be explained by a low dimensional spatiotemporal structure. More specifically, they can be accounted for by the combination of a small number of basic activation patterns. The basic patterns and distribution weights indicate temporal and spatial structures, respectively, and the weights show the muscle sets that are activated synchronously. In addition, various locomotor behaviors have similar low dimensional structures and major differences appear in the basic patterns. These analysis results suggest that neural systems use muscle group combinations to solve motor control redundancy problems (muscle synergy hypothesis) and manipulate those basic patterns to create various locomotor functions. However, it remains unclear how the neural system controls such muscle groups and basic patterns through neuromechanical interactions in order to achieve adaptive locomotor behavior. This paper reviews simulation studies that explored adaptive motor control in locomotion via sensory-motor coordination using neuromusculoskeletal models based on the muscle synergy hypothesis. Herein, the neural mechanism in motor control related to the muscle synergy for adaptive locomotion and a potential muscle synergy analysis method including neuromusculoskeletal modeling for motor impairments and rehabilitation are discussed.