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

Excitatory Spinal Lhx9-Derived Interneurons Modulate Locomotor Frequency in Mice

Locomotion allows us to move and interact with our surroundings. Spinal networks that control locomotion produce rhythm and left-right and flexor-extensor coordination. Several glutamatergic populations, Shox2 non-V2a, Hb9-derived interneurons, and, recently, spinocerebellar neurons have been proposed to be involved in the mouse rhythm generating networks. These cells make up only a smaller fraction of the excitatory cells in the ventral spinal cord. Here, we set out to identify additional populations of excitatory spinal neurons that may be involved in rhythm generation or other functions in the locomotor network. We use RNA sequencing from glutamatergic, non-glutamatergic, and Shox2 cells in the neonatal mice from both sexes followed by differential gene expression analyses. These analyses identified transcription factors that are highly expressed by glutamatergic spinal neurons and differentially expressed between Shox2 neurons and glutamatergic neurons. From this latter category, we identified the Lhx9-derived neurons as having a restricted spinal expression pattern with no Shox2 neuron overlap. They are purely glutamatergic and ipsilaterally projecting. Ablation of the glutamatergic transmission or acute inactivation of the neuronal activity of Lhx9-derived neurons leads to a decrease in the frequency of locomotor-like activity without change in coordination pattern. Optogenetic activation of Lhx9-derived neurons promotes locomotor-like activity and modulates the frequency of the locomotor activity. Calcium activities of Lhx9-derived neurons show strong left-right out-of-phase rhythmicity during locomotor-like activity. Our study identifies a distinct population of spinal excitatory neurons that regulates the frequency of locomotor output with a suggested role in rhythm-generation in the mouse alongside other spinal populations.

Distinguishing subtypes of spinal locomotor neurons to inform circuit function and dysfunction

Locomotion is a complex motor task executed by spinal neurons. Given the diversity of spinal cord neurons, linking neuronal cell type to function is a challenge. Molecular identification of broad spinal interneuronal classes provided a great advance. Recent studies have used other classifiers, including location, electrophysiological properties, and connectivity, in addition to gene profiling, to narrow the acuity with which groups of neurons can be related to specific functions. However, there are also functional populations without a clear identifier, as exemplified by rhythm generating neurons. Other considerations, including experience or plasticity, add a layer of complexity to the definition of functional subpopulations of spinal neurons, but spinal cord injury may provide insight.

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.

Synaptic connectivity amongst components of the locomotor central pattern generator

In the past two decades we have learned an enormous amount of information regarding the identity of functional components of the neural circuitry responsible for generating locomotor activity in mammals. Molecular techniques, combined with classic electrophysiological and anatomical approaches, have resulted in the identification of a handful of classes of genetically defined interneuronal populations, and a delineation of the specific function of many of these during stepping. What lags behind at this point is a clear picture of the synaptic connectivity of each population, this information is key if we are to understand how the interneuronal components that are responsible for locomotor activity work together to form a functional circuit. In this mini review I will summarize what is, and what is not, known regarding the synaptic connectivity of each genetically defined interneuronal population that is involved in locomotion.

Identification of adult spinal Shox2 neuronal subpopulations based on unbiased computational clustering of electrophysiological properties

Spinal cord neurons integrate sensory and descending information to produce motor output. The expression of transcription factors has been used to dissect out the neuronal components of circuits underlying behaviors. However, most of the canonical populations of interneurons are heterogeneous and require additional criteria to determine functional subpopulations. Neurons expressing the transcription factor Shox2 can be subclassified based on the co-expression of the transcription factor Chx10 and each subpopulation is proposed to have a distinct connectivity and different role in locomotion. Adult Shox2 neurons have recently been shown to be diverse based on their firing properties. Here, in order to subclassify adult mouse Shox2 neurons, we performed multiple analyses of data collected from whole-cell patch clamp recordings of visually-identified Shox2 neurons from lumbar spinal slices. A smaller set of Chx10 neurons was included in the analyses for validation. We performed k-means and hierarchical unbiased clustering approaches, considering electrophysiological variables. Unlike the categorizations by firing type, the clusters displayed electrophysiological properties that could differentiate between clusters of Shox2 neurons. The presence of clusters consisting exclusively of Shox2 neurons in both clustering techniques suggests that it is possible to distinguish Shox2+Chx10- neurons from Shox2+Chx10+ neurons by electrophysiological properties alone. Computational clusters were further validated by immunohistochemistry with accuracy in a small subset of neurons. Thus, unbiased cluster analysis using electrophysiological properties is a tool that can enhance current interneuronal subclassifications and can complement groupings based on transcription factor and molecular expression.

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.

The role of V3 neurons in speed-dependent interlimb coordination during locomotion in mice

Speed-dependent interlimb coordination allows animals to maintain stable locomotion under different circumstances. The V3 neurons are known to be involved in interlimb coordination. We previously modeled the locomotor spinal circuitry controlling interlimb coordination (Danner et al., 2017). This model included the local V3 neurons that mediate mutual excitation between left and right rhythm generators (RGs). Here, our focus was on V3 neurons involved in ascending long propriospinal interactions (aLPNs). Using retrograde tracing, we revealed a subpopulation of lumbar V3 aLPNs with contralateral cervical projections. V3OFF mice, in which all V3 neurons were silenced, had a significantly reduced maximal locomotor speed, were unable to move using stable trot, gallop, or bound, and predominantly used a lateral-sequence walk. To reproduce this data and understand the functional roles of V3 aLPNs, we extended our previous model by incorporating diagonal V3 aLPNs mediating inputs from each lumbar RG to the contralateral cervical RG. The extended model reproduces our experimental results and suggests that locally projecting V3 neurons, mediating left-right interactions within lumbar and cervical cords, promote left-right synchronization necessary for gallop and bound, whereas the V3 aLPNs promote synchronization between diagonal fore and hind RGs necessary for trot. The model proposes the organization of spinal circuits available for future experimental testing.

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.

Asymmetric and transient properties of reciprocal activity of antagonists during the paw-shake response in the cat

Mutually inhibitory populations of neurons, half-center oscillators (HCOs), are commonly involved in the dynamics of the central pattern generators (CPGs) driving various rhythmic movements. Previously, we developed a multifunctional, multistable symmetric HCO model which produced slow locomotor-like and fast paw-shake-like activity patterns. Here, we describe asymmetric features of paw-shake responses in a symmetric HCO model and test these predictions experimentally. We considered bursting properties of the two model half-centers during transient paw-shake-like responses to short perturbations during locomotor-like activity. We found that when a current pulse was applied during the spiking phase of one half-center, let’s call it #1, the consecutive burst durations (BDs) of that half-center increased throughout the paw-shake response, while BDs of the other half-center, let’s call it #2, only changed slightly. In contrast, the consecutive interburst intervals (IBIs) of half-center #1 changed little, while IBIs of half-center #2 increased. We demonstrated that this asymmetry between the half-centers depends on the phase of the locomotor-like rhythm at which the perturbation was applied. We suggest that the fast transient response reflects functional asymmetries of slow processes that underly the locomotor-like pattern; e.g., asymmetric levels of inactivation across the two half-centers for a slowly inactivating inward current. We compared model results with those of in-vivo paw-shake responses evoked in locomoting cats and found similar asymmetries. Electromyographic (EMG) BDs of anterior hindlimb muscles with flexor-related activity increased in consecutive paw-shake cycles, while BD of posterior muscles with extensor-related activity did not change, and vice versa for IBIs of anterior flexors and posterior extensors. We conclude that EMG activity patterns during paw-shaking are consistent with the proposed mechanism producing transient paw-shake-like bursting patterns found in our multistable HCO model. We suggest that the described asymmetry of paw-shaking responses could implicate a multifunctional CPG controlling both locomotion and paw-shaking.

The existence of multifunctional central pattern generators (CPGs), circuits which control more than one rhythmic motor behavior, is an intriguing hypothesis. We suggest that the cat paw-shaking response could be a transient response of the locomotor CPG. Our general prediction is that this CPG is multifunctional, and in addition to the locomotor rhythm, it can generate a transient, ten-times faster, paw-shake-like response to a stimulus. In our multistable half-center oscillator (HCO) CPG model, we applied perturbations to the locomotor pattern which resulted in a transient paw-shake-like pattern that eventually returned back to the locomotor pattern. We showed that the inactivation of the slow inward current that drives the locomotor rhythm produced asymmetry of the transient flexor and extensor activity in a symmetric HCO model. To test predictions from our model about the transient nature of the paw-shake response, we compared burst durations (BDs) and interburst intervals (IBIs) of the model half-centers in consecutive cycles of paw-shake-like responses with the BD and IBI of electromyographic (EMG) activity bursts of cat hindlimb flexors and extensors recorded during a paw-shake response. In both cases, we found similar asymmetric trends in the BD and IBI throughout a paw-shake response.

Optogenetic Activation of V1 Interneurons Reveals the Multimodality of Spinal Locomotor Networks in the Neonatal Mouse

In the spinal cord, classes of interneurons have been studied in vitro to determine their role in producing or regulating locomotion. It is unclear whether all locomotor behaviors are produced by the same circuitry or engage different subsets of neurons. Here, in neonatal mice of either sex, we test this idea by comparing the actions of a class of spinal, inhibitory interneuron (V1) expressing channelrhodopsin driven by the engrailed-1 transcription factor on the rhythms elicited by different methods. We find that, although the overall locomotor activities in vitro are similar, V1 interneuron depolarization produces opposite effects depending of the mode of activation of the locomotor circuitry. The differential behavior of V1 neurons suggests that their function depends on how the locomotor rhythm is activated and is consistent with the idea that the functional organization of the corresponding locomotor networks also differs.SIGNIFICANCE STATEMENT The neural networks dictating the execution of fictive locomotion are located in the spinal cord. It is generally assumed that the mode of activation of these spinal networks should not change the recruitment or function of neurons. Here, we manipulated the activity of a class of interneuron (V1), which targets these networks and found that their activation induces opposite effects depending on the mode of activation. This suggests that the mode of activation of the spinal networks differentially recruits either V1 interneurons or other interneurons, or both.

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.

Unusual Quadrupedal Locomotion in Rat during Recovery from Lumbar Spinal Blockade of 5-HT7 Receptors

Coordination of four-limb movements during quadrupedal locomotion is controlled by supraspinal monoaminergic descending pathways, among which serotoninergic ones play a crucial role. Here we investigated the locomotor pattern during recovery from blockade of 5-HT₇ or 5-HT_2A receptors after intrathecal application of SB269970 or cyproheptadine in adult rats with chronic intrathecal cannula implanted in the lumbar spinal cord. The interlimb coordination was investigated based on electromyographic activity recorded from selected fore- and hindlimb muscles during rat locomotion on a treadmill. In the time of recovery after hindlimb transient paralysis, we noticed a presence of an unusual pattern of quadrupedal locomotion characterized by a doubling of forelimb stepping in relation to unaffected hindlimb stepping (2FL-1HL) after blockade of 5-HT₇ receptors but not after blockade of 5-HT_2A receptors. The 2FL-1HL pattern, although transient, was observed as a stable form of fore-hindlimb coupling during quadrupedal locomotion. We suggest that modulation of the 5-HT₇ receptors on interneurons located in lamina VII with ascending projections to the forelimb spinal network can be responsible for the 2FL-1HL locomotor pattern. In support, our immunohistochemical analysis of the lumbar spinal cord demonstrated the presence of the 5-HT₇ immunoreactive cells in the lamina VII, which were rarely 5-HT_2A immunoreactive.

Spinal Inhibitory Interneurons: Gatekeepers of Sensorimotor Pathways

The ability to sense and move within an environment are complex functions necessary for the survival of nearly all species. The spinal cord is both the initial entry site for peripheral information and the final output site for motor response, placing spinal circuits as paramount in mediating sensory responses and coordinating movement. This is partly accomplished through the activation of complex spinal microcircuits that gate afferent signals to filter extraneous stimuli from various sensory modalities and determine which signals are transmitted to higher order structures in the CNS and to spinal motor pathways. A mechanistic understanding of how inhibitory interneurons are organized and employed within the spinal cord will provide potential access points for therapeutics targeting inhibitory deficits underlying various pathologies including sensory and movement disorders. Recent studies using transgenic manipulations, neurochemical profiling, and single-cell transcriptomics have identified distinct populations of inhibitory interneurons which express an array of genetic and/or neurochemical markers that constitute functional microcircuits. In this review, we provide an overview of identified neural components that make up inhibitory microcircuits within the dorsal and ventral spinal cord and highlight the importance of inhibitory control of sensorimotor pathways at the spinal level.

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.

Inhibitory control of active expiration by the Bötzinger complex in rats

KEY POINTS

Contraction of abdominal muscles at the end of expiration during metabolic challenges (such as hypercapnia and hypoxia) improves pulmonary ventilation. The emergence of this active expiratory pattern requires the recruitment of the expiratory oscillator located on the ventral surface of the medulla oblongata. Here we show that an inhibitory circuitry located in the Bötzinger complex is an important source of inhibitory drive to the expiratory oscillator. This circuitry, mediated by GABAergic and glycinergic synapses, provides expiratory inhibition that restrains the expiratory oscillator under resting condition and regulates the formation of abdominal expiratory activity during active expiration. By combining experimental and modelling approaches, we propose the organization and connections within the respiratory network that control the changes in the breathing pattern associated with elevated metabolic demand.

ABSTRACT

The expiratory neurons of the Bötzinger complex (BötC) provide inhibitory inputs to the respiratory network, which, during eupnoea, are critically important for respiratory phase transition and duration control. Here, we investigated how the BötC neurons interact with the expiratory oscillator located in the parafacial respiratory group (pFRG) and control the abdominal activity during active expiration. Using the decerebrated, arterially perfused in situ preparations of juvenile rats, we recorded the activity of expiratory neurons and performed pharmacological manipulations of the BötC and pFRG during hypercapnia or after the exposure to short-term sustained hypoxia - conditions that generate active expiration. The experimental data were integrated in a mathematical model to gain new insights into the inhibitory connectome within the respiratory central pattern generator. Our results indicate that the BötC neurons may establish mutual connections with the pFRG, providing expiratory inhibition during the first stage of expiration and receiving excitatory inputs during late expiration. Moreover, we found that application of GABAergic and glycinergic antagonists in the BötC caused opposing effects on abdominal expiratory activity, suggesting complex inhibitory circuitry within the BötC. Using mathematical modelling, we propose that the BötC network organization and its interactions with the pFRG restrain abdominal activity under resting conditions and contribute to abdominal expiratory pattern formation during active expiration observed during hypercapnia or after the exposure to short-term sustained hypoxia.

Emergence of Different Gaits in Infancy: Relationship Between Developing Neural Circuitries and Changing Biomechanics

How does gait-specific pattern generation evolve in early infancy? The idea that neural and biomechanical mechanisms underlying mature walking and running differ to some extent and involve distinct spinal and supraspinal neural circuits is supported by various studies. Here we consider the issue of human gaits from the developmental point of view, from neonate stepping to adult mature gaits. While differentiating features of the walk and run are clearly distinct in adults, the gradual and progressive developmental bifurcation between the different gaits suggests considerable sharing of circuitry. Gaits development and their biomechanical determinants also depend on maturation of the musculoskeletal system. This review outlines the possible overlap in the neural and biomechanical control of walking and running in infancy, supporting the idea that gaits may be built starting from common, likely phylogenetically conserved elements. Bridging connections between movement mechanics and neural control of locomotion could have profound clinical implications for technological solutions to understand better locomotor development and to diagnose early motor deficits. We also consider the neuromuscular maturation time frame of gaits resulting from active practice of locomotion, underlying plasticity of development.

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.

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.

Expanding our horizons: central pattern generation in the context of complex activity sequences

Central pattern generators (CPGs) are central nervous system (CNS) networks that can generate coordinated output in the absence of patterned sensory input. For decades, this concept was applied almost exclusively to simple, innate, rhythmic movements with essentially identical cycles that repeat continually (e.g. respiration) or episodically (e.g. locomotion). But many natural movement sequences are not simple rhythms, as they include different elements in a complex order, and some involve learning. The concepts and experimental approaches of CPG research have also been applied to the neural control of complex movement sequences, such as birdsong, though this is not widely appreciated. Experimental approaches to the investigation of CPG networks, both for simple rhythms and for complex activity sequences, have shown that: (1) brief activation of the CPG elicits a long-lasting naturalistic activity sequence; (2) electrical stimulation of CPG elements alters the timing of subsequent cycles or sequence elements; and (3) warming or cooling CPG elements respectively speeds up or slows down the rhythm or sequence rate. The CPG concept has also been applied to the activity rhythms of populations of mammalian cortical neurons. CPG concepts and methods might further be applied to a variety of fixed action patterns typically used in courtship, rivalry, nest building and prey capture. These complex movements could be generated by CPGs within CPGs (‘nested’ CPGs). Stereotypical, non-motor, non-rhythmic neuronal activity sequences may also be generated by CPGs. My goal here is to highlight previous applications of the CPG concept to complex but stereotypical activity sequences and to suggest additional possible applications, which might provoke new hypotheses and experiments.

Mapping Connectivity Amongst Interneuronal Components of the Locomotor CPG

The basic rhythmic activity characteristic of locomotion in mammals is generated by a neural network, located in the spinal cord, known as the locomotor central pattern generator (CPG). Although a great deal of effort has gone into the study of this neural circuit over the past century, identification and characterization of its component interneurons has proven to be challenging, largely due to their location and distribution. Recent work incorporating a molecular approach has provided a great deal of insight into the genetic identity of interneurons that make up this neural circuit, as well as the specific roles that they play during stepping. Despite this progress we still know relatively little regarding the manner in which these neuronal populations are interconnected. In this article we review the interneuronal populations shown to be involved in locomotor activity, briefly summarize their specific function, and focus on experimental work that provides insight into their synaptic connectivity. Finally, we discuss how recently developed viral approaches can potentially be incorporated to provide further insight into the network structure of this neural circuit.

V1 interneurons regulate the pattern and frequency of locomotor-like activity in the neonatal mouse spinal cord

In the mouse spinal cord, V1 interneurons are a heterogeneous population of inhibitory spinal interneurons that have been implicated in regulating the frequency of the locomotor rhythm and in organizing flexor and extensor alternation. By introducing archaerhodopsin into engrailed-1-positive neurons, we demonstrate that the function of V1 neurons in locomotor-like activity is more complex than previously thought. In the whole cord, V1 hyperpolarization increased the rhythmic synaptic drive to flexor and extensor motoneurons, increased the spiking in each cycle, and slowed the locomotor-like rhythm. In the hemicord, V1 hyperpolarization accelerated the rhythm after an initial period of tonic activity, implying that a subset of V1 neurons are active in the hemicord, which was confirmed by calcium imaging. Hyperpolarizing V1 neurons resulted in an equalization of the duty cycle in flexor and extensors from an asymmetrical pattern in control recordings in which the extensor bursts were longer than the flexor bursts. Our results suggest that V1 interneurons are composed of several subsets with different functional roles. Furthermore, during V1 hyperpolarization, the default state of the locomotor central pattern generator (CPG) is symmetrical, with antagonist motoneurons each firing with an approximately 50% duty cycle. We hypothesize that one function of the V1 population is to set the burst durations of muscles to be appropriate to their biomechanical function and to adapt to the environmental demands, such as changes in locomotor speed.

An optogenetic study in mice shows that inhibitory neurons that express engrailed-1 regulate the pattern and frequency of locomotor-like activity in the developing mouse spinal cord.

Current Principles of Motor Control, with Special Reference to Vertebrate Locomotion

The vertebrate control of locomotion involves all levels of the nervous system from cortex to the spinal cord. Here, we aim to cover all main aspects of this complex behavior, from the operation of the microcircuits in the spinal cord to the systems and behavioral levels and extend from mammalian locomotion to the basic undulatory movements of lamprey and fish. The cellular basis of propulsion represents the core of the control system, and it involves the spinal central pattern generator networks (CPGs) controlling the timing of different muscles, the sensory compensation for perturbations, and the brain stem command systems controlling the level of activity of the CPGs and the speed of locomotion. The forebrain and in particular the basal ganglia are involved in determining which motor programs should be recruited at a given point of time and can both initiate and stop locomotor activity. The propulsive control system needs to be integrated with the postural control system to maintain body orientation. Moreover, the locomotor movements need to be steered so that the subject approaches the goal of the locomotor episode, or avoids colliding with elements in the environment or simply escapes at high speed. These different aspects will all be covered in the review.

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.

Computational modeling of brainstem circuits controlling locomotor frequency and gait

A series of recent studies identified key structures in the mesencephalic locomotor region and the caudal brainstem of mice involved in the initiation and control of slow (exploratory) and fast (escape-type) locomotion and gait. However, the interactions of these brainstem centers with each other and with the spinal locomotor circuits are poorly understood. Previously we suggested that commissural and long propriospinal interneurons are the main targets for brainstem inputs adjusting gait (Danner et al., 2017). Here, by extending our previous model, we propose a connectome of the brainstem-spinal circuitry and suggest a mechanistic explanation of the operation of brainstem structures and their roles in controlling speed and gait. We suggest that brainstem control of locomotion is mediated by two pathways, one controlling locomotor speed via connections to rhythm generating circuits in the spinal cord and the other providing gait control by targeting commissural and long propriospinal interneurons.

Spinal Shox2 interneuron interconnectivity related to function and development

Neuronal networks generating hindlimb locomotion are located in the spinal cord. The mechanisms underlying spinal rhythmogenesis are unknown but network activity and interconnectivity of excitatory interneurons likely play prominent roles. Here, we investigate interconnectivity within the Shox2 interneuron population, a subset of which has been suggested to be involved in locomotor rhythm generation, using paired recordings in isolated spinal cords or slices from transgenic mice. Sparse unidirectional connections consistent with chemical synaptic transmission and prominent bidirectional connections mediated by electrical synapses were present within distinct subsets of Shox2 interneurons. Moreover, bidirectional electrical connections were preferentially found between functionally-related Shox2 interneurons. Though prevalent in neonatal mice, electrical coupling began to decline in incidence and strength in mice ~ 3 weeks of age. Overall, our data suggest that gap junctional coupling promotes synchronization of Shox2 interneurons, and may be implicated in locomotor rhythmicity in developing mice.

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.

Control of transitions between locomotor-like and paw shake-like rhythms in a model of a multistable central pattern generator

The ability of the same neuronal circuit to control different motor functions is an actively debated concept. Previously, we showed in a model that a single multistable central pattern generator (CPG) could produce two different rhythmic motor patterns, slow and fast, corresponding to cat locomotion and paw shaking. A locomotor-like rhythm (~1 Hz) and a paw shake-like rhythm (~10 Hz) did coexist in our model, and, by applying a single pulse of current, we could switch the CPG from one regime to another (Bondy B, Klishko AN, Edwards DH, Prilutsky BI, Cymbalyuk G. In: Neuromechanical Modeling of Posture and Locomotion, 2016). Here we investigated the roles of slow intrinsic ionic currents in this multistability. The CPG is modeled as a half-center oscillator circuit comprising two reciprocally inhibitory neurons. Each neuron is equipped with two slow inward currents, a Na⁺ current ( I_NaS) and a Ca²⁺ current ( I_CaS). I_CaS inactivates much more slowly and at more hyperpolarized voltages than I_NaS. We demonstrate that I_NaS is the primary current driving the paw shake-like bursting. I_CaS is crucial for the locomotor-like bursting, and it is inactivated during the paw shake-like activity. We investigate the sensitivity of the bursting regimes to perturbations, using a pulse of current to induce a switch from one regime to the other, and we demonstrate that the transition duration is dependent on pulse amplitude and application phase. We also investigate the modulatory roles of the strength of various currents on characteristics of these rhythms and show that their effects are regime specific. We conclude that a multistable CPG is physiologically plausible and derive testable predictions of the model. NEW & NOTEWORTHY Little is known about how a single central pattern generator could produce multiple rhythms. We describe a novel mechanism for multistability of bursting regimes with strongly distinct periods. The proposed mechanism emphasizes the role of intrinsic cellular dynamics over synaptic dynamics in the production of multistability. We describe how the temporal characteristics of multiple rhythms could be controlled by neuromodulation and how single pulses of current could produce a switch between regimes in a functional fashion.

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.

Modular organization of muscle activity patterns in the leading and trailing limbs during obstacle clearance in healthy adults

Human locomotor patterns require precise adjustments to successfully navigate complex environments. Studies suggest that the central nervous system may control such adjustments through supraspinal signals modifying a basic locomotor pattern at the spinal level. To explore this proposed control mechanism in the leading and trailing limbs during obstructed walking, healthy young adults stepped over obstacles measuring 0.1 and 0.2 m in height. Unobstructed walking with no obstacle present was also performed as a baseline. Full body three-dimensional kinematic data were recorded and electromyography (EMG) was collected from 14 lower limb muscles on each side of the body. EMG data were analyzed using two techniques: by mapping the EMG data to the approximate location of the motor neuron pools on the lumbosacral enlargement of the spinal cord and by applying a nonnegative matrix factorization algorithm to unilateral and bilateral muscle activations separately. Results showed that obstacle clearance may be achieved not only with the addition of a new activation pattern in the leading limb, but with a temporal shift of a pattern present during unobstructed walking in both the leading and trailing limbs. An investigation of the inter-limb coordination of these patterns suggested a strong bilateral linkage between lower limbs. These results highlight the modular organization of muscle activation in the leading and trailing limbs, as well as provide a mechanism of control when implementing a locomotor adjustment when stepping over an obstacle.

The neural control of interlimb coordination during mammalian locomotion

Neuronal networks within the spinal cord directly control rhythmic movements of the arms/forelimbs and legs/hindlimbs during locomotion in mammals. For an effective locomotion, these networks must be flexibly coordinated to allow for various gait patterns and independent use of the arms/forelimbs. This coordination can be accomplished by mechanisms intrinsic to the spinal cord, somatosensory feedback from the limbs, and various supraspinal pathways. Incomplete spinal cord injury disrupts some of the pathways and structures involved in interlimb coordination, often leading to a disruption in the coordination between the arms/forelimbs and legs/hindlimbs in animal models and in humans. However, experimental spinal lesions in animal models to uncover the mechanisms coordinating the limbs have limitations due to compensatory mechanisms and strategies, redundant systems of control, and plasticity within remaining circuits. The purpose of this review is to provide a general overview and critical discussion of experimental studies that have investigated the neural mechanisms involved in coordinating the arms/forelimbs and legs/hindlimbs during mammalian 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.