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

Abstract
Introduction
Results
Discussion
Materials and methods
References
Article and author information
Metrics

Abstract

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.

https://doi.org/10.7554/eLife.31050.001

Introduction

To effectively move in the environment, limbed animals use different gaits depending on the desired locomotor speed (Grillner, 1981; Hildebrand, 1989). Mice change gait from walk to trot followed by gallop and bound as locomotor speed increases (Clarke and Still, 1999; Herbin et al., 2004; Herbin et al., 2007; Bellardita and Kiehn, 2015; Lemieux et al., 2016). Each gait is characterized by a distinct pattern of limb movements (Hildebrand, 1989), which in turn is defined by neural interactions between spinal circuits controlling each limb.

It is generally accepted that each limb is controlled by a separate central pattern generator (CPGs), because cats are able to walk on split-belt treadmills with limbs stepping at different speeds (Forssberg et al., 1980; Halbertsma, 1983; Frigon et al., 2013, 2015; Thibaudier et al., 2013; Frigon, 2017). Hence, quadrupedal locomotion is controlled by four CPGs (one per limb) located in separate sections of the spinal cord, specifically on the left and right side in the lumbar and cervical enlargements (Kato, 1990; Ballion et al., 2001; Juvin et al., 2005; Juvin et al., 2012). Therefore, the central control of interlimb coordination and gait expression is defined by neurons and neuronal pathways that mediate interactions between circuits controlling each limb. These interlimb pathways are mediated by (a) short projecting commissural interneurons (CINs) that control left-right interactions at the lumbar and, presumably, cervical levels (Talpalar et al., 2013; Bellardita and Kiehn, 2015; Molkov et al., 2015; Shevtsova et al., 2015; Danner et al., 2016), and (b) descending (cervical-to-lumbar) and ascending (lumbar-to-cervical) long propriospinal neurons (LPNs) that control homolateral and diagonal cervical-lumbar interactions (Matsushita et al., 1979; Skinner et al., 1979, 1980; Menétrey et al., 1985; Nathan et al., 1996; Dutton et al., 2006; Brockett et al., 2013; Ruder et al., 2016; Flynn et al., 2017).

The role and speed-dependent involvement of genetically identified CINs, specifically V0_V and V0_D, in alternation and expression of different gaits have been experimentally investigated (Talpalar et al., 2013; Bellardita and Kiehn, 2015). Specifically, it was shown that mutants lacking V0_V CINs selectively lose trot, and mutants lacking both V0_V and V0_D CINs can only bound at any speed (Bellardita and Kiehn, 2015; Kiehn, 2016). Using modeling (Danner et al., 2016), we presented a possible organization of the central pathways mediating interlimb coordination that exhibited speed-dependent gait expression and the correct loss of gaits after selective removal of V0_V and both V0_V and V0_D CINs. These V0 interneurons may contain both short projecting CINs as well as diagonally projecting LPNs (Ruder et al., 2016).

LPNs that directly couple the cervical and lumbar segments have been anatomically identified (Edinger, 1896; Nathan and Smith, 1959; Matsushita et al., 1979; Skinner et al., 1979; Skinner et al., 1980; Menétrey et al., 1985; Nathan et al., 1996; Dutton et al., 2006; Brockett et al., 2013; Flynn et al., 2017). Several studies in vitro found that LPNs can be involved in coordination of left-right and cervical-lumbar activities recorded from four ventral roots (Ballion et al., 2001; Juvin et al., 2005; Juvin et al., 2012; Cowley et al., 2010). However, their role in speed-dependent gait expression in vivo, remains poorly understood. Recently, Ruder et al. (2016) described axonal projections of different LPNs and their synaptic and genetic properties. These data differ from the connections proposed in our previous model (Danner et al., 2016). In addition, Ruder et al. (2016) have shown that deletion of descending (cervical-to-lumbar) LPNs causes a speed-dependent distortion of left-right coordination. These new findings provided us with the information needed to revise our model and investigate the role of different LPN pathways in interlimb coordination and gait control during locomotion.

In locomotion induced by stimulation of the mesencephalic locomotor region (MLR), increasing the intensity of stimulation causes an increase of locomotor speed accompanied by sequential gait transitions (Orlovskiĭ et al., 1966; Shik et al., 1966; Shik and Orlovsky, 1976; Skinner and Garcia-Rill, 1984; Grillner, 1985; Atsuta et al., 1990; Nicolopoulos-Stournaras and Iles, 2009). From these observations, we hypothesized that the same brainstem drive that controls speed by exciting the rhythm generators also affects interlimb coordination by acting on CINs or LPNs (Danner et al., 2016). This hypothesis is in accordance with multiple data confirming that different CINs and LPNs receive supraspinal inputs (Lloyd, 1942; Lloyd and McIntyre, 1948; Bannatyne et al., 2003; Jankowska et al., 2003; Cowley et al., 2010; Zaporozhets et al., 2011). However, the involvement of these neurons in mediating brainstem control of gait has not been demonstrated. Therefore, organization and role of brainstem inputs to these neurons can so far only be based on predictions derived from computational models.

Here, we present a computational model of the spinal locomotor circuits and investigate the potential role of CINs and LPNs in the control of and transition between different locomotor gaits. The model reproduces the speed-dependent gait expression in intact and genetically transformed mice, lacking V2a, V0_V, and all V0 neurons (Crone et al., 2009; Bellardita and Kiehn, 2015), as well as the disruption of hindlimb coordination following ablation of descending LPNs (Ruder et al., 2016). The model proposes the following roles for LPNs: diagonal V0_D LPNs supports walk, diagonal V0_V LPNs together with local V0_V CINs stabilize trot, and homolateral LPNs ensure fore-hind alternation in all gaits. Finally, we show that additional external inputs to CINs and LPNs may affect interlimb coordination and gait expression independent of speed. We suggest that these spinal interneurons are the main targets of supraspinal and somatosensory afferents to adjust interlimb coordination and gait to different environmental and behavioral conditions.

Results

The model

The model consists of four rhythm generators (RGs), each controlling one limb (Figure 1). These RGs interact via several bidirectional pathways mediated by CINs (left-right) and LPNs (homolateral and diagonal). Neural populations were modeled with a simplified, 'activity-dependent’ population model (see Materials and methods). Each RG consisted of a flexor (F) and extensor center (E) that mutually inhibited each other through intermediate inhibitory interneurons (InF and InE, respectively). Similar to our previous models (Rybak et al., 2006a, 2006b, 2013, 2015; Molkov et al., 2015; Shevtsova et al., 2015; Danner et al., 2016), both centers incorporated a slowly inactivating, persistent sodium current ( $I_{NaP}$ ), allowing each center to intrinsically generate rhythmic bursting under certain conditions defined by external tonic drives or level of excitation, which was consistent with experimental data (Hägglund et al., 2013). However, we followed the concept of asymmetric flexion-dominated CPG organization (Pearson and Duysens, 1976; Duysens, 1977; Zhong et al., 2012; Duysens et al., 2013; Machado et al., 2015). Based on this concept and our previous models (Rybak et al., 2015; Molkov et al., 2015; Shevtsova et al., 2015; Danner et al., 2016; Ausborn et al., 2017), we assumed that under normal conditions extensor centers receive relatively high drive that keeps them in the regime of tonic activity. Hence, they exhibited rhythmic bursting only due to rhythmic inhibition from the corresponding intrinsically oscillating flexor centers. In this case, the frequency of oscillation in the model was defined by the brainstem drive to flexor centers and was almost independent of the phasic interactions between RGs mediated by CINs and LPNs (Rybak et al., 2015).

Figure 1

Download asset Open asset

Model concept.

Each limb is controlled by a separate rhythm generator (RG). The local commissural neurons (CINs) as well as homolaterally and diagonally projecting (descending and ascending) long propriospinal neurons (LPNs) couple the four RGs. Brainstem drive acts on the RGs to control locomotor speed and on CINs and diagonal LPNs to control gaits.

https://doi.org/10.7554/eLife.31050.002

Left-right CIN connections at cervical and lumbar levels (Figure 2A and Figure 3) are organized in accordance to our previous models (Rybak et al., 2015; Shevtsova et al., 2015; Danner et al., 2016): the inhibitory V0_D CINs provide direct mutual inhibition between the left and right flexor centers, the excitatory V0_V CINs also provide mutual inhibition between the flexor centers (receiving inputs from excitatory V2a and acting through inhibitory Ini neurons), the excitatory V3 CINs provide mutual excitation between the flexor centers, and the inhibitory CINi CINs provide inhibition from the extensor to the contralateral flexor centers.

Figure 2

Download asset Open asset

Connections within the spinal cord.

(A) Connections between left and right rhythm generators (RG) within each girdle. (B) Connections from the fore to hind RGs via descending (cervical-to-lumbar) long propriospinal neurons (LPNs). (C) Connections from the hind to fore RGs via ascending (lumbar-to-cervical) LPNs. (D) Brainstem drive connections to the extensor and flexor centers, commissural interneurons, and LPNs. Spheres represent neural populations. Excitatory and inhibitory connections are marked by arrowheads and circles, respectively.

https://doi.org/10.7554/eLife.31050.003

Figure 3

Download asset Open asset

Full model schematic.

Spheres represent neural populations and lines synaptic connections. Excitatory and inhibitory connections are marked by arrowheads and circles, respectively. RG, rhythm generator.

https://doi.org/10.7554/eLife.31050.004

Connections between the cervical and lumbar circuits are mediated through mono- and polysynaptic neural pathways. Here, we modeled them as LPN connections. Ruder et al. (2016) reported that there are excitatory and inhibitory descending LPNs, whereas ascending LPNs are only excitatory. Both descending and ascending LPNs project homolaterally and diagonally. There are more diagonal then homolateral excitatory LPNs, while the opposite is the case for inhibitory LPNs. The excitatory homolateral LPNs were identified as Shox2 and the excitatory diagonal LPNs as V0_V neurons. We have hypothesized that diagonally projecting inhibitory LPNs are of V0_D type.

The cervical-to-lumbar connections in our model (Figure 2B) are organized accordingly: the excitatory Shox2 neurons provide excitation from each cervical extensor to its homolateral lumbar flexor center, inhibitory LPNs of unidentified type (LPNi) provide inhibition from each cervical flexor to its homolateral lumbar flexor center, V0_V LPNs provide excitation from each cervical flexor to its diagonal flexor center, and inhibitory V0_D LPNs provide inhibition from each cervical flexor to its diagonal lumbar flexor center.

Lumbar-to-cervical connections (Figure 2C) mirror the excitatory cervical-to-lumbar connections but do not include inhibitory LPNs.

Finally, there were two brainstem drives (Figure 2D): one that excited flexor centers and inhibited the local V0 CINs at cervical and lumbar levels and the descending diagonal V0_D LPNs, promoting left-right and diagonal alternation between the corresponding RGs, and the other that provided constant excitatory drive to the extensor centers. The drive to the flexor centers, CINs, and LPNs was varied to influence locomotor speed and gait. Specifically, the drive to flexor centers defined the locomotor frequency. Based on our suggestion, this drive also inhibited V0 CINs and diagonal V0_D LPNs changing the balance of excitatory and inhibitory interactions between the corresponding RGs. The full model schematic can be seen in Figure 3.

The model was used to analyze drive- and frequency-dependent changes in the locomotor phase durations and the expression of different gaits in an intact system and following the removal of particular neuron types.

Control of locomotor frequency and gait by brainstem drive

We investigated the model performance by changing the parameter $α$ , which defined brainstem drive to flexor centers, CINs, and LPNs (see Materials and methods). The model generated oscillations when parameter $α$ was changed from 0 to 1.05. Within this range, an increase of $α$ led to an increase in locomotor frequency from 2 to 12 Hz (Figure 4A and B, top diagram). This increase in frequency occurred mostly due to shortening of the extension phase, while the flexion phase remained almost constant (Figure 4B, bottom diagram), which is a characteristic property of fictive and real locomotion across species (Halbertsma, 1983; Hildebrand, 1989; Frigon and Gossard, 2009; Danner et al., 2015; Shevtsova et al., 2015). The increased frequency in the model was accompanied by sequential gait changes from walk to trot to gallop and then bound (Figure 4A) consistent with experimental observations (Figure 4C, Bellardita and Kiehn, 2015).

Figure 4

Download asset Open asset

Model performance.

(A) Gait changes during 6 s of simulation with increase of brainstem drive (parameter $α$ ) every second (shown at the bottom). Gait changes were subject to a brief transitional period (shown as transitional colors). (B) Dependency of frequency and phase durations on $α$ . Vertical dashed lines (labeled as $x_{1} - x_{4}$ ) represent selected values of $α$ related to different gaits. Note that the blue line almost completely covers the red one (indicating increase and decrease of $α$ , respectively) in the top diagram. (C) Frequency ranges in which gaits are observed in mice (created from data from Figure 1E of Bellardita and Kiehn, 2015). (**D–G**) Phase differences and step cycle compositions of four characteristic gaits produced by the model that correspond to the $α$ -values ( $x_{1} - x_{4}$ ) in B.

https://doi.org/10.7554/eLife.31050.005

The expressed walking gait was a lateral-sequence walk (Figure 4A,D), in which flexion phases proceeded in the following order: hind-right, fore-right, hind-left and fore-left. This walk occurred at frequencies from 2 to 4 Hz (Figure 4B). At low frequencies, the extension phase was three to four times longer than the flexion phase duration, the RGs were in flexion sequentially, and the homolateral and diagonal phase differences were close to 0.25 and 0.75, respectively (Figure 4D).

Trot (Figure 4A,E) was characterized by diagonal synchronization, left-right alternation and fore-hind alternation. It was expressed over a large range of frequencies (from 4 to 10.5 Hz; Figure 4B). At low frequencies, with longer extension than flexion phase durations, the gait constituted a walking trot and otherwise a running trot.

Gallops are characterized by synchronization or quasi-synchronization of the hind RGs and a non-zero phase difference between the fore RGs. In the model (Figure 4A,B,F), gallop was expressed at frequencies between 9 and 11 Hz (between 9 and 10.5 Hz both gallop and trot occurred). The sequence of flexion phases at cervical or lumbar levels could be either left-right or right-left but was always the same for both fore and hind RGs. Thus, transverse but not rotary gallops were expressed. The left-right phase difference of the hind RGs was always closer to synchronization than of the fore RGs. Some of the gallops exhibited an interval between the end of the hind extension and the beginning of the fore extension, in which all RGs were in flexion at the same time (Figure 4F). This corresponds to an aerial phase, characteristic for gallops.

Bound (Figure 4A,B,G) was characterized by left-right synchronization at cervical and lumbar levels, and homolateral and diagonal alternation. The fore extension phases were directly followed by hind extension phases, which were then followed by all RGs being in flexion (corresponding to an aerial phase). Bound was expressed at frequencies between 11 and 12 Hz. The flexion phase duration was always longer than the extension phase duration (Figure 4B, bottom diagram).

The frequency ranges, phase differences, and step cycle compositions (relative phase durations) of all gaits were consistent with mouse locomotion (Figure 4B–G; Bellardita and Kiehn, 2015; see gait definitions in Materials and methods).

Walk and the transition to trot: diagonal V0_D LPNs ensure stable walk

To investigate control of locomotion and gait transitions by brainstem drive, we used bifurcation diagrams reflecting normalized phase differences between oscillations in hind and fore left-right, homolateral and diagonal RGs with changing parameter $α$ that represented the brainstem drive and was used as bifurcation parameters (see examples in Figure 5). Normalized phase differences of 0.5 correspond to alternation, whereas phase differences of 0 or 1 correspond to synchronization. Gaits were operationally defined based on these phase differences (see Table 2 in Materials and methods) and are marked in the diagrams by colored areas. Blue and red lines indicate the stable phase differences with stepwise increase and decrease of $α$ , respectively. Any discrepancies between the red and blue lines indicate regions of bi- or multistability. For these $α$ values two or more different stable gaits coexisted and could be expressed depending on the initial conditions. Bifurcations can be seen as abrupt changes of the stable phase differences. Figure 5A shows the bifurcation diagrams of the intact model.

Figure 5

Download asset Open asset

Bifurcation diagrams of the intact model (A), after removal of all V0_V neurons (B), of only V0_V LPNs (C), and of all V0_V and V0_D neurons (D).

Normalized phase differences of 0.5 correspond to alternation, whereas phase differences of 0 or 1 correspond to synchronization. Blue and red lines indicate the stable phase differences with stepwise increase and decrease of the bifurcation parameter $α$ , respectively. Colored areas indicate the expressed gait.

https://doi.org/10.7554/eLife.31050.006

At the lowest values of brainstem drive, both left-right (local V0_D and V0_V CINs) and diagonal alternation (diagonal V0_D LPNs) promoting pathways dominated. Together with three to four times longer extension than flexion phase durations (Figure 4B bottom diagram) and homolateral alternation (present over all frequencies), this ensured that only one RG was in flexion at any given time (Figure 4A,D and Figure 5A).

With increasing $α$ , inhibitory brainstem drive caused the activity of the diagonally projecting V0_D LPNs to decrease and the balance to shift from diagonal alternation (promoted by V0_D) to diagonal synchronization (promoted by V0_V). Together with a progressive decrease in the relative duration of the extension phase, this caused the flexion phases of the diagonal RGs to progressively overlap and the diagonal phase difference to move from a quarter-phase lag toward synchronization (Figure 5A). In parallel, the homolateral phase difference moved toward alternation and the left-right phase differences remained at a strict out-of-phase alternation (0.5 normalized phase difference). Together this constitutes a gradual transition from walk to trot.

The sequence of the flexion phases during walk was determined by the hypothetical CINi population that was introduced in our previous model (CINi2 in Danner et al., 2016) to ensure 0.5 normalized left-right phase difference during walk.

Trot and the transition to gallop and bound: diagonal V0_V LPNs stabilize trot

The transition from trot to gallop and bound was governed by brainstem-drive-controlled speed-dependent changes in the left-right and diagonal pathways. The left-right alternation promoting local CINs (V0_D and V0_V) received inhibition from the brainstem that increased with speed, while the left-right synchronization promoting local CINs (V3 and CINi) did not receive brainstem inhibition (Figure 2D, Figure 3). Thus, with increasing brainstem drive (and frequency), the pathway promoting left-right synchronization became stronger than those promoting left-right alternation. This resulted in a transition from trot to gallop and then to bound (Figure 5A). The diagonally projecting V0_V LPNs, promoting diagonal synchronization, acted synergistically with local V0_V CINs.

Two sequential bifurcations occurred: when the gait switched from trot to gallop and then from gallop to bound (Figure 5A). The bifurcation from trot to gallop resembled a subcritical (backwards) pitchfork bifurcation in the hind (and fore) left-right phase differences. The transition exhibited hysteresis: perfect left-right alternation and quarter-phase lags (in either direction) were stable at the same $α$ -values (blue and red lines show different stable left-right phase differences in the orange shaded area in Figure 5A). In other words, when drive was increasing (blue lines), the transition from trot to gallop occurred at a higher $α$ -value than the transition from gallop to trot when drive was decreasing (red lines). Such a hysteresis is a common feature of gait transitions in quadrupeds (Heglund and Taylor, 1988). Stable gallops occurred because the decreasing V0_V CINs and LPNs activities were accompanied by a decrease of the burst duration, while V3 and CINi CINs were active during the whole flexion or extension phase, respectively. Note, that the previous model (Danner et al., 2016) could not reproduce stable gallops; they only existed as a transitional regime. Therefore, the addition of diagonal LPNs in the present model was critical for stable gallop to exist. The bifurcation from trot to gallop occurred when V0_V CINs and LPNs could not support left-right alternation and diagonal synchronization any more, but V3 and CINi CINs were not strong enough to cause perfect left-right synchronization. Therefore, stable states emerged around 0/1 normalized left-right phase differences. The deviation from synchronization was determined by the duration of the bursts of V0_V CINs and LPNs. As the burst durations decreased further, the left-right phase differences transitioned towards synchronization (and the diagonal phase difference towards alternation) via a supercritical pitchfork bifurcation. This bifurcation provided the transition from gallop to bound without hysteresis (blue and red lines of the left-right phase differences in Figure 5A follow the same trajectory).

Gait expression following removal of different spinal interneuron types

Experimental observations showed that mice lacking V0_V neurons lose trot, mice lacking V0_V and V0_D only exhibit bound (Bellardita and Kiehn, 2015), and mice lacking V2a neurons transition to left-right synchronization at lower speeds compared to intact animals (Crone et al., 2009). Thus, to further validate the model, we removed several types of genetically identified neurons and compared the model performance with experimental results.

When all V0_V neurons were deleted (Figure 5B) only walk, gallop and bound were expressed; trot was selectively lost. Left-right synchronization promoting pathways became dominant at a lower $α$ -value (lower frequency), and left-right alternation was only supported by V0_D CINs. Deletion of V2a had the same effect as deletion of V0_V, because V2a relayed all inputs to V0_V CINs and LPNs.

Deletion of only diagonally projecting V0_V LPNs (and not V0_V CINs) did not result in a loss of any gait, but caused the transition from trot to gallop to occur at a lower $α$ -value and frequency (Figure 5C). Furthermore, an additional bifurcation in the left-right and diagonal phase differences occurred during trot. At low trot-frequencies and $α$ -values, a strict out-of-phase left-right alternation (0.5 normalized phase difference) was stable, while at medium trot-frequencies, a pair of stable states near 0.5 normalized phase difference emerged that then transitioned via a region of multistability to a pair of stable states near 0/1 normalized phase difference (expression of gallop). This highlights how local V0_V CINs and LPNs synergistically stabilize trot by providing left-right alternation and diagonal synchronization.

Deletion of both V0_V and V0_D resulted in the removal of all left-right alternation promoting pathways and caused the loss of all alternating gaits (Figure 5D). Bound was stable over all frequencies.

Deletion of descending (cervical-to-lumbar) LPNs affects left-right coordination

Ruder et al. (2016) showed that deletion of descending LPNs influences left-right coordination: at low speeds both the fore and hind limbs were alternating, then at medium speeds the hind limbs, and at high speeds both the fore and hind limbs demonstrated disturbed coordination with both alternating and synchronized activities.

In our model, deletion of the cervical-to-lumbar LPNs resulted in the emergence of new stable states (Figure 6B). At low brainstem drive values walk remained the only stable state and with increasing drive it gradually transitioned towards trot. With further increase of $α$ , trot was first the only stable gait (Figure 6B,C). Then, both trot and gallop were stable at the same $α$ -values (Figure 6B,D1,D2,E1,E2) and finally bound became the single stable state (Figure 6B,F). During the gallop at low $α$ -values, the hind RGs were close to synchronization while the fore RGs remained alternating (Figure 6B,D2). With increasing $α$ , the fore left-right phase differences gradually transitioned towards synchronization, thus gallops at higher $α$ -values exhibited fore and hind left-right phase differences close to synchronization (Figure 6B,E2).

Figure 6

Download asset Open asset

Model performance after deletion of all cervical-to-lumbar long propriospinal neurons (LPNs).

Bifurcation diagrams, indicating the stable phase differences depending on the bifurcation parameter $α$ for the intact model (A) and after deletion of all cervical-to-lumbar LPNs (B). (**C–F**) Extension phases of stable gaits (identified by vertical dashed lines in B labeled as $x_{1} - x_{4}$ ). At $α = {x_{2}, x_{3}}$ two different gaits were stable simultaneously (depicted in D1, D2 and E1, E2, respectively).

https://doi.org/10.7554/eLife.31050.007

The introduced asymmetry by the deletion of pathways in only one direction (from cervical to lumbar but not from lumbar to cervical) caused bifurcations of higher codimension and resulted in several bi- and multistabilities (Figure 6B). Interestingly, the stable states closely approximate the union of the stable states of the intact model Figure 6A and the model after deletion of all diagonally projecting V0_V LPNs (Figure 5C). Thus, the disturbance of the fore and hind left-right phase differences was mainly caused by the deletion of diagonally projecting V0_V LPNs. The remaining homolaterally projecting LPNs maintained fore-hind alternation after deletion. Local V0_V CINs maintained the stable branch of left-right alternation.

Noise causes high step-to-step variability after deletion of cervical-to-lumbar LPNs

In the above section, we described that deletion of the cervical-to-lumbar LPNs resulted in the emergence of new, additional stable states at medium and high brainstem drive values. This means that two gaits were stable at the same $α$ -values (e.g. see $x_{2}$ in Figure 6B and Figure 6D1,D2). To consider how these changes in the number of steady states affect system behavior, such as step-to-step variability, we increased the noisy currents in all populations ( $σ_{Noise}$ from 0.005 pA to 1.75 pA). In the intact model, at medium brainstem drives (related to trot), the increased noise caused phase durations and phase differences to become more variable, but both, the fore and hind RGs remained alternating over all steps (Figure 7A). After deletion of cervical-to-lumbar LPNs, the applied noise caused spontaneous transitions between the steady states: the activity of hind RGs was spontaneously switching between synchronization and alternation while the fore RGs remained alternating (see dashed box in Figure 7B).

Figure 7

Download asset Open asset

Model performance under application of noise before (intact) and after ablation of cervical-to-lumbar long projecting propriospinal neurons (LPNs).

(**A,B**) Exemplary extension phases of the intact model (A) and after removal of cervical-to-lumbar LPNs (B). The dashed box in B indicates synchronization of hind RGs that transiently occurred after removal of cervical-to-lumbar LPNs. (**C–E**,**C1–E1**) Phase differences categorized into three equally sized bins (gray: appropriate for trot, yellow: quarter phase difference off, red: antiphase) before and after ablation of cervical-to-lumbar LPNs for 1000 s simulation in the model (**C–E**) and for experimental data pooled across animals (**C1–E1**) at three different speeds and $α$ -values. C1–E1 were created from data extracted from Figure 5B,C and S5A,D of Ruder et al. (2016).

https://doi.org/10.7554/eLife.31050.008

A similar disruption of interlimb coordination has been described in mice after deletion of cervical-to-lumbar LPNs (Ruder et al., 2016). They reported random occurrences of clusters of misaligned steps and spontaneous transitions between alternation and synchronization, which were similar to those exhibited by our model. Furthermore, the model also reproduced the speed-dependent effect of changes in the proportion of appropriate and misaligned steps: at low speeds hind and fore limb coordination remained appropriate for trot, at medium speeds hind limb coordination was disturbed, and at high speeds both fore and hind limb coordination was affected (Figure 7C–E). Thus, in the regions where both trot and gallop were stable, the applied noise caused spontaneous transitions between the stable states.

Gait changes by brainstem drive and by drive-independent inputs to CINs and LPNs

Mice usually switch gaits and locomotor speeds abruptly with few or no intermediary steps (Bellardita and Kiehn, 2015). To simulate the dynamics of gait transitions, we initialized the model with a constant $α$ -value that defined a certain gait (and frequency) and then abruptly changed $α$ to a value characteristic of another gait. Transitions from walk to trot, from gallop to trot and from walk to gallop and back to walk (Figure 8A–C) mirrored experimental observations closely (Figure 8D–E; Bellardita and Kiehn, 2015): they usually required only one or two cycles to stabilize the new gait and locomotor frequency.

Figure 8

Download asset Open asset

Gait transitions with changes of brainstem drive (**A–C**) and independent of brainstem drive (**G–I**).

(**A–C**) Instantaneous frequency, extension phases before and after parameter $α$ was abruptly changed (bottom trace). Black vertical lines indicate the time when $α$ was changed and the shaded green areas indicate the transitions periods. (A) Transition from walk to trot when $α$ was changed from 0.02 to 0.4. (B) Transition from gallop to trot when $α$ was changed from 0.85 to 0.6. (C) Transition from walk to gallop and back to walk when $α$ was changed from 0.02 to 0.9 and to 0.02. (**D–F**) Experimentally observed instantaneous speed during gait changes in mice corresponding to A–C (created from data extracted from Figure 3A–C of Bellardita and Kiehn 2015). (**G–I**) Instantaneous frequency and extension phases during gait changes caused by additional drives to CINs and LPNs. Bottom trace indicates time course of additional drives. (G) Additional inhibitory drive to all V0_V CINs and LPNs ( $m_{I} = 0.0$ , $b_{I} = 0.2$ ) at $α = 0.5$ and caused a transition from trot to bound. During the transitional period a gallop occurred. (H) An additional excitatory drive to local V0_V CINs ( $m_{E} = 0.0$ , $b_{E} = 0.05$ ) at $α = 0.925$ caused a transition from gallop to trot. (I) Additional excitatory drive to cervical, local V0_V CINs ( $m_{E} = 0.0$ , $b_{E} = 0.1$ ) at $α = 0.975$ caused the left-right phase difference of the fore RGs to change from almost synchronized to a quarter-phase lag during gallop.

https://doi.org/10.7554/eLife.31050.009

Finally, gait changes could also be induced by additional inputs to CINs and LPNs without changing the common brainstem drive controlling locomotor speed (frequency). For example, additional inhibition to V0_V CINs and LPNs during trot induced a transition to bound without significantly affecting the frequency, resulting in a bound at a frequency characteristic of trot (Figure 8G). Excitatory drive to the local V0_V CINs during gallop induced a transition to trot, again without significantly affecting speed (Figure 8H). Excitation of the local cervical V0_V CINs during gallop caused the fore left-right phase differences to move farther from synchronization and the left and right extension phases of the fore RGs to overlap less, while the hind RG coordination remained almost unchanged (Figure 8I). Thus, tonic inputs to CINs and LPNs were able to alter interlimb coordination and gait almost independently of frequency.

Discussion

Control of locomotor speed and speed-dependent gait transitions

Central control of locomotor speed and speed-dependent gait expression involves many neural, biomechanical, metabolic, environmental and behavioral factors. Most of these factors are beyond the scope of the present study. Here, we focused on the potential role of central interactions between rhythm-generating circuits controlling each limb and their regulation by brainstem drives.

In our model, two such drives were considered: constant drive to the extensor centers and variable drive to the flexor centers and various CINs and LPNs. The drive to the extensor centers put them in a tonic activity regime if the flexor centers were not active. This drive may represent a separate descending supraspinal drive to extensor circuits operating through the medial reticulospinal or vestibulospinal tracts and/or inputs from cutaneous afferents and load receptors from extensor muscles (Hiebert and Pearson, 1999; Dietz and Duysens, 2000; Bouyer and Rossignol, 2003; Rossignol et al., 2006).

The brainstem drive to the flexor centers put them into a rhythmic bursting regime and strong inhibition between flexor and extensor centers mediated by inhibitory interneurons caused rhythmic flexor-extensor alternation (Zhang et al., 2014; Britz et al., 2015; Rybak et al., 2015; Shevtsova et al., 2015; Danner et al., 2016; Shevtsova and Rybak, 2016). Progressive increase of this drive increased locomotor frequency, in accordance with studies in decerebrate cats and rats, where stimulation of the MLR induced locomotor activity that increased in frequency (speed) with increasing stimulation (Orlovskiĭ et al., 1966; Shik et al., 1966; Shik and Orlovsky, 1976; Skinner and Garcia-Rill, 1984; Grillner, 1985; Atsuta et al., 1990; Nicolopoulos-Stournaras and Iles, 2009). In many of these studies, the increase in speed was accompanied by sequential changes in gait from walk to trot and to gallop and bound. The exact mechanisms for these speed-dependent gait transitions are currently unknown, but it is reasonable to suggest that these transitions are initiated and supported by inputs and drives from different excitatory and inhibitory neurons in the medullary and pontomedullary reticular formation that project to the spinal cord and can directly affect different CINs and LPNs (Alstermark et al., 1987; Jankowska et al., 2003, 2005; Matsuyama et al., 2004; Mitchell et al., 2016; Ruder et al., 2016).

Following our previous model (Danner et al., 2016), we suggest that gait transitions are caused by brainstem drives acting on neurons coupling the RGs and affecting the balance between the activities of pathways promoting synchronization and alternation. In the previous model, only homolateral LPNs and local CINs were considered. The transition from alternating (trot) to synchronized gaits (gallop and bound) resulted from excitatory drive to local V3 CINs, promoting left-right synchronization in both cervical and lumbar levels, so that with increasing drive the left-right interactions through these pathways overcame interactions via pathways supporting left-right alternation.

In the current model, the same conceptual idea was implemented by increasing inhibitory influence of the brainstem drive to local V0 CINs (V0_D and V0_V) and diagonal V0_D LPNs (Figure 2D and Figure 3). This increasing inhibition influenced both the walk to trot and the trot to gallop and bound transitions: the walk to trot transition occurred when diagonal V0_V became stronger than diagonal V0_D LPNs, and the trot to gallop and bound transition occurred when local V3 became stronger than V0_V CINs and LPNs. Additionally, this inhibitory drive allowed for stable solutions for gallop, which in the previous model only occurred as transient solutions (Danner et al., 2016).

Unfortunately, there is no experimental evidence for or against brainstem excitation of V3 or inhibition of V0 CINs and LPNs. Either mechanisms could exist and operate in reality. Moreover, they could even coexist and cooperate in multiple ways. This issue can be resolved in future experiments by recording from identified CINs and LPNs and analyzing their response to MLR stimulation.

Long propriospinal neurons and the organization of interactions between rhythm-generating circuits

The left-right coordination in the present model is mediated by both the local CINs, coupling left and right RGs at cervical and lumbar levels, and the diagonal LPNs, that couple cervical and lumbar RGs. The local CIN interactions (Figure 2A and Figure 3) were implemented based on earlier experimental studies (Bellardita and Kiehn, 2015) and followed our previous computational model (Danner et al., 2016). The homolateral and diagonal LPN-mediated pathways (Figure 2B,C and Figure 3) were incorporated to be consistent with (Ruder et al., 2016). The bidirectional homolateral excitatory (from each extensor to the corresponding flexor center, mediated by Shox2 interneurons) and the descending inhibitory (between the flexor centers, mediated by LPNi interneurons) LPN pathways acted synergistically to ensure fore-hind alternation over all frequencies and gaits. The descending diagonal inhibitory LPNs (V0_D) together with the local inhibitory V0_D CINs ensured stable quarter-phase lags between homolateral and diagonal RGs during walk. The diagonal ascending and descending excitatory LPNs (V0_V) promoted diagonal synchronization and (together with the inhibitory homolateral interactions) supported left-right alternation during trot.

The presence of only cervical-to-lumbar but not lumbar-to-cervical inhibitory LPNs (Ruder et al., 2016) and stronger lumbar-to-cervical excitatory influence than cervical-to-lumbar (Juvin et al., 2005; Brockett et al., 2013) introduced a fore-hind asymmetry. Thus, the dynamics of our model were consistent with left-right but not fore-hind interchange (Schöner et al., 1990; Golubitsky et al., 1999). The left-right symmetry (because of symmetric LPN and CIN interactions) resulted in stable normalized left-right phase differences of 0.5 (in walk and trot), or 0/1 (in bound), or pairs of steady states symmetric around 0.5 (in gallops). These symmetry properties resulted in the stability of both left and right leading gallops. The fore-hind asymmetry caused the lateral-sequence walk to be the only stable walk, which is the preferred walk of several quadrupeds, including mice, dogs, cats, and humans crawling (Patrick et al., 2009; Bellardita and Kiehn, 2015; Righetti et al., 2015; Frigon, 2017), as well as the correct extension phase sequences during gallop and bound. Afferent inputs to LPNs could potentially affect the symmetry properties and result in different stable solutions. For example, cats exhibit diagonal-sequence walk when walking on a split belt treadmill with the fore and hind limbs at a different speed (Thibaudier et al., 2013).

Role of genetically identified spinal interneurons in interlimb coordination

We used our model to simulate the effects of ablating some genetically identified spinal interneurons on the speed-dependent expression and transition of different gaits. The removal of V0_V neurons abolished trot (Figure 5B) and only bound could be expressed after the removal of both V0_V and V0_D (Figure 5D), which is fully consistent with the experimental data of Bellardita and Kiehn (2015). Also, because in our model V2a neurons mediate input to V0_V neurons (Figure 2 and Figure 3; see Crone et al., 2008; Talpalar et al., 2013; Rybak et al., 2015; Shevtsova et al., 2015), removal of V2a had the same effect as removing V0_V: transition to gallop and bound at a lower frequency, consistent with Crone et al. (2009).

It is necessary to note that the role of V0_V neurons in the expression and support of trot in the present model significantly differs from our previous model (Danner et al., 2016). The previous model had been developed before the new data on LPN organization (Ruder et al., 2016) became available and only included homolateral LPNs and not diagonal ones. Therefore, in that model, the role of V0_V neurons in supporting trot was entirely based on their involvement in support of left-right alternation. In the present model, both V0_V CINs and LPNs support trot, by providing left-right alternation (CINs) and diagonal synchronization (LPNs). The model suggests that removal of one type of these pathways (e.g. the diagonal LPNs, Figure 5C) does not fully eliminate trot, but shifts transition to gallop to a lower locomotor frequency. Transition to gallop results from the increasing inhibitory influence of the brainstem drive to V0_V neurons with increasing frequency. The predicted inhibitory influence of brainstem drive (e.g. from MLR) to identified V0_V neurons (both local CINs and LPNs) awaits experimental testing.

Effect of ablation of the descending cervical-to-lumbar LPN pathways

Simulating deletion of cervical-to-lumbar connections was performed for additional model validation in relation to the recent experimental data: Ruder et al. (2016) showed that deletion of these connections, did not affect interlimb coordination at low speeds, disturbed left-right coordination between hind limbs at medium speeds, exhibiting spontaneous switching to left-right synchronization, and caused the same disturbances in left-right coordination between both the fore and hind limbs at high speed. In the model, this deletion resulted in additional steady states with almost synchronous left-right activities of hind RGs at the medium drive values (frequencies) and slow migration of such additional steady states to left-right synchronization of fore RGs (Figure 6B).

Introduction of noise led to random transitions between these coexisting stable states, and resulted in spontaneous changes in phase relationships between hind RGs (including occasional synchronization) at medium frequencies and both hind and fore RGs at higher frequencies (Figure 7). This qualitatively corresponded to the experimental results on ablation of cervical-to-lumber connections (Ruder et al., 2016) and hence provided additional validation to our model.

Commissural and long propriospinal neurons can serve as main targets for supraspinal and sensory afferent signals controlling limb coordination

Although limbed animals, including quadrupeds express gaits according to their speed, there are many other factors affecting selection of the appropriate gaits, including energetic/metabolic (Hoyt and Taylor, 1981), biomechanical (Alexander and Jayes, 1983; Biewener, 1990; Farley and Taylor, 1991; Fukuoka et al., 2015), environmental and behavioral factors. For example, animals need to maintain balance when maneuvering in dynamic environments, or when chasing prey or escaping from a predator. As locomotor activities are controlled by networks within the spinal cord, signals affecting interlimb coordination to accommodate these factors can be (a) common or specific supraspinal inputs or drives, reflecting additional visual, auditory, vestibular and different goal-directed factors, and/or (b) proprioceptive and cutaneous afferent feedbacks from the periphery that informs the spinal controller about various biomechanical constraints. The supraspinal control of gait transitions has been confirmed by multiple experiments during fictive or real locomotion controlled by brainstem stimulation (Orlovskiĭ et al., 1966; Shik et al., 1966; Shik and Orlovsky, 1976; Skinner and Garcia-Rill, 1984; Grillner, 1985; Atsuta et al., 1990; Nicolopoulos-Stournaras and Iles, 2009). At the same time, there are experimental data confirming the important role that proprioceptive and exteroceptive feedback plays in the control of interlimb coordination and gait (reviewed by Pearson, 1995; Pearson, 2004; Frigon, 2017).

The locomotor gait of quadrupeds is defined by the phase relationships of the movements of the different limbs. Assuming that each limb is controlled by a separate RG, each gait can be represented as a set of phase differences between the rhythmic activities generated by the RGs. In quadrupeds, all left-right interactions within the cervical and lumbar spinal cord, including interactions between left and right RG circuits, are mediated by local CINs, and all diagonal and homolateral interactions between the cervical and lumbar circuits, including interactions between the corresponding RGs, are mediated by LPNs with descending and ascending, diagonal and homolateral projections, respectively. Therefore, the activities of CINs and LPNs located within the four sections of the cord (left and right hemicords at cervical and lumbar levels) explicitly define interactions between the four rhythm-generating circuits (or RGs; Figure 9). An integrated pattern of activity of these neurons defines interlimb coordination and gait. Multiple studies confirm that spinal CINs and LPNs receive excitatory and inhibitory inputs from supraspinal structures (Skinner et al., 1979; Skinner et al., 1980; Alstermark et al., 1987; Bannatyne et al., 2003; Jankowska et al., 2003; Jankowska et al., 2005; Jankowska et al., 2006; Matsuyama et al., 2004; Cowley et al., 2008; Cowley et al., 2010; Jankowska, 2008; Szokol et al., 2011; Zaporozhets et al., 2011; Juvin et al., 2012; Mitchell et al., 2016; Ruder et al., 2016), which could explicitly provide speed-dependent gait control as in our model. At the same time, as we showed in our simulations (Figure 8G–I), inputs to CINs and LPNs other than supraspinal drives controlling locomotor frequency can also induce gait changes independent of speed. In particular, these inputs can come from sensory afferents. This is consistent with multiple experimental data showing that CINs and LPNs receive direct and indirect input from sensory afferents (Lloyd and McIntyre, 1948; Edgley et al., 2003; Bannatyne et al., 2006; Jankowska et al., 2009; Flynn et al., 2017). In accord with this is the notion that afferent feedback from each limb projects mostly to the spinal cord section controlling the same limb (Frigon, 2017) and hence its influence on spinal circuits and RGs controlling other limbs should be mostly mediated by CINs and LPNs located in this section.

Figure 9

Download asset Open asset

Schematic representation of control of limb coordination and gait by different local commissural interneurons (CINs) and long propriospinal neurons (LPNs).

(A) Circuits controlling four limbs. Each limb is controlled by its own rhythm generator (RG). Local CINs and homolateral and diagonal LPNs couple the RGs and coordinate their activities. Supraspinal signals and sensory feedbacks (directly or also indirectly through dorsal horn interneurons; not shown) provide inputs to the RGs, CINs and LPNs. Inputs to the RGs affect the locomotor frequency (speed), while inputs to CINs and LPNs affect interlimb coordination and gait expression. (B) A more detailed representation of local CINs and LPNs that are located in one cord section (here left cervical) and integrate the corresponding intraspinal and supraspinal inputs and sensory information (from the fore left limb) to mediate the effect of this limb’s activity on spinal circuits controlling the other three limbs.

https://doi.org/10.7554/eLife.31050.010

Therefore, we suggest that CINs and LPNs represent the main neural targets for different local/intraspinal, supraspinal, and sensory inputs to control interlimb coordination and adjust gait to various internal and external conditions. Inputs to other neural structures, such as the RG centers, could also alter interlimb coordination by transsynaptically influencing the balance of the CIN and LPN pathways, but would also affect other parameters, such as locomotor frequency. Separate (independent of brainstem drives) signals to CINs and LPNs only (e.g. from sensory afferents) can affect limb coordination and change gait independent of speed (as shown in Figure 8G–I). Figure 9 summarizes the neural circuitry within the spinal cord responsible for interlimb coordination and possible inputs from supraspinal and peripheral afferents. For simplicity, in this schematic, the same CIN and LPN populations receive all intraspinal, supraspinal and afferent inputs. In real cord, the effects of these inputs are rather mediated by separate or overlapping populations of CINs and LPNs, which may interact with each other (e.g. via presynaptic inhibition) competing for their influence on gait expression.

Model limitations

In this study, we focused only on central interactions, without considering biomechanical constraints and afferent feedbacks from the limbs. Mechanical coupling between the limbs and body can also influence interlimb coordination and interact with the neural control (Nishikawa et al., 2007). Moreover, gaits are not only characterized by the phase differences between the limbs but also by a change in kinematics and muscle activation patterns (Bellardita and Kiehn, 2015). Investigating these issues requires simulation of limb and body biomechanics in addition to the spinal circuits, which is outside the scope of the current study. We also did not consider spinal circuits operating below the RGs, such as pattern formation networks, circuits related to muscle afferent inputs, reflex circuits including Ia, and Ib interneurons, Renshaw cells and motoneurons (see Rybak et al., 2006a, 2006b; McCrea and Rybak, 2007, 2008; Zhong et al., 2012), which play an important role in controlling the limbs. Moreover, how inputs from various supraspinal structures interact with different spinal mechanisms and afferent feedback at different spinal levels in the context of controlling interlimb coordination remains largely unknown. All the above will be the focus of our future investigations.

Materials and methods

Neuron model

Request a detailed protocol

The model consists of a network of interconnected populations of neurons. Each population was represented by a non-spiking, ‘activity-based’ model (Ermentrout, 1994) and described by

C \cdot d V / d t = - I_{NaP} - I_{L} - I_{SynE} - I_{SynI} - I_{Noise}

for flexor and extensor centers, and by

C \cdot d V / d t = - I_{L} - I_{SynE} - I_{SynI} - I_{Noise}

for all other populations. In Equations 1 and 2, $V$ represents the average membrane potential, $C$ the membrane capacitance, $I_{NaP}$ the persistent sodium current, $I_{L}$ the leakage current, $I_{SynE}$ and $I_{SynI}$ excitatory and inhibitory synaptic currents, respectively, and $I_{Noise}$ a noisy current.

The leakage current was described by

I_{L} = g_{L} \cdot (V - E_{L}),

where $g_{L}$ is the leakage conductance and $E_{L}$ the leakage reversal potential.

Excitatory and inhibitory synaptic currents ( $I_{SynE}$ and $I_{SynI}$ , respectively) of population $i$ were described by:

\begin{aligned} (4) & I_{S y n E, i} & = g_{S y n E} \cdot {\sum_{j} [S (w_{j i}) \cdot f (V_{j})] + D_{E, i}} \cdot (V_{i} - E_{S y n E}); \\ (5) & I_{S y n I, i} & = g_{S y n I} \cdot {\sum_{j} [S (- w_{j i}) \cdot f (V_{j})] + D_{I, i}} \cdot (V_{i} - E_{S y n I}); \end{aligned}

where $g_{SynE}$ and $g_{SynI}$ are synaptic conductances, and $E_{SynE}$ and $E_{SynI}$ the reversal potentials of the excitatory and inhibitory synapses; $w_{j i}$ is the synaptic weight from population $j$ to $i$ ( $w_{j i} > 0$ for excitatory connections and $w_{j i} < 0$ for inhibitory connections); and function $S$ is defined as

S (x) = {\begin{matrix} x, & i f x \geq 0 \\ 0, & o t h e r w i s e \end{matrix} .

The output function $f (V)$ translates $V$ into the integrated population activity or neural output and was defined by the linear piecewise function:

f (V) = {\begin{aligned} 0, & i f V < V_{t h r} \\ (V - V_{t h r}) / (V_{m a x} - V_{t h r}), & i f V_{t h r} \leq V < V_{m a x} \\ 1, & i f V \geq V_{m a x} \end{aligned} .

Excitatory $D_{E, i}$ and inhibitory $D_{I, i}$ drives to population $i$ were modeled as a linear function of the free parameter $α$ :

D_{{E, I}, i} (α) = m_{{E, I}, i} \cdot α + b_{{E, I}, i},

where $m_{{E, I}, i}$ is the slope, and $b_{{E, I}, i}$ the intercept. $α$ represents the value of the variable brainstem drive. If not otherwise specified, $m_{{E, I}, i}$ and $b_{{E, I}, i}$ were set to 0 (lack of drive input to the corresponding neuron population).

The persistent sodium current in the flexor and extensor centers of the RG was described by

I_{NaP} = {\bar{g}}_{NaP} \cdot m \cdot h \cdot (V - E_{Na}),

where ${\bar{g}}_{NaP}$ is the maximal conductance, $m$ and $h$ activation and inactivation variables of $I_{NaP}$ , respectively, and $E_{Na}$ the sodium reversal potential. Activation of the persistent sodium current was considered instant and was modeled by

m (V) = {1 + \exp [(V - V_{1 / 2, m}) / k_{m}]}^{- 1}

and slow inactivation obeyed the following differential equation

τ_{h} (V) \cdot d h / d t = h_{\infty} (V) - h

with

\begin{aligned} (12) & h_{\infty} (V) & = {1 + \exp [(V - V_{1 / 2, h}) / k_{h}]}^{- 1}; \\ (13) & τ_{h} (V) & = τ_{0} + (τ_{m a x} - τ_{o}) / \cosh [(V - V_{1 / 2, τ}) / k_{τ}] . \end{aligned}

In Equations 10–13, $V_{1 / 2}$ and $k$ represent half-voltage and slope of the corresponding variable, respectively, and $τ_{0}$ and $τ_{\max}$ are the baseline and maximum of inactivation time constant $τ_{h}$ , respectively.

The noisy current $I_{Noise}$ was described as an Ornstein-Uhlenbeck process

d I_{Noise} / d t = - I_{Noise} / τ_{Noise} + σ_{Noise} \sqrt{2 / τ_{Noise}} \cdot ξ_{i} (t)

where $τ_{Noise}$ is the time constant, $σ_{Noise}$ the standard deviation or strength of the noise, and $ξ_{i} (t)$ (in $1 / \sqrt{s}$ ) was the population-specific normalized Gaussian noise.

Model parameters

Request a detailed protocol

The model parameters were adapted from our previous model (Danner et al., 2016): $C = 10$ pF; $g_{L} = 4.5$ nS for RG centers and $g_{L} = 2.8$ nS for all other neurons; ${\bar{g}}_{NaP} = 4.5$ nS; $g_{SynE} = g_{SynI} = 10.0$ nS; $E_{L} = - 62.5$ mV for RG centers and $E_{L} = - 60$ mV for all other neurons; $E_{Na} = 50$ mV; $E_{SynE} = - 10$ mV; $E_{SynI} = - 75$ mV; $V_{thr} = - 50$ mV; $V_{\max} = 0$ mV; $V_{1 / 2, m} = - 40$ mV; $V_{1 / 2, h} = - 45$ mV; $k_{m} = - 6$ mV; $k_{h} = 4$ mV; $τ_{0} = 80$ ms; $τ_{\max} = 160$ ms; $V_{1 / 2, τ} = - 35$ mV; $k_{τ} = 15$ mV; and $τ_{Noise} = 10$ ms.

Following drive parameters were used: $m_{E} = 0.0$ and $b_{E} = 0.1$ for the extensor centers; $m_{E} = 0.1$ and $b_{E} = 0.0$ for the flexor centers; $m_{I} = 0.75$ and $b_{I} = 0.0$ for all V0_D CINs; and $m_{I} = 0.15$ and $b_{I} = 0.0$ for the local homologous V0_V CINs. Thus, the extensor centers received constant drive (independent of $α$ ) and the flexor centers and V0 CINs received variable drive whose changes were proportional to $α$ . The connection weights are listed in Table 1. Weights of connections within the RGs were adapted from our previous model (Danner et al., 2016) and other weights were selected within their operating ranges and tuned to allow gait transitions to occur at experimentally observed locomotor frequencies (Bellardita and Kiehn, 2015). To test the robustness of the model, we simultaneously varied all connection weights; each weight was multiplied by a normally distributed random number with a mean of 1 and standard deviation $σ_{p}$ of 0.01, 0.02, 0.05, and 0.10. For each $σ_{p}$ , 100 random models were created and bifurcation diagrams were calculated. With $σ_{p} \leq 0.02$ all randomized models retained all stable regimes and their sequential transitions with changes of $α$ , at $σ_{p} = 0.05$ 17% and at $σ_{p} = 0.10$ 40% of the models lost some stable solutions (gaits such as bound or trot). Thus, the final model represents a coarse system allowing parameter variations without dramatic (qualitative) changes in behavior. To simulate deletion of particular populations, their outputs were set to 0.

Table 1

Connection weights.

https://doi.org/10.7554/eLife.31050.011

Source	Target ( $w_{i j}$ )
	Within girdle and side of the cord
RG-F	i-Ini-F (0.40), i-V0D (0.70), i-V2a (1.00), i-V3 (0.35), i-V2a-diag (0.50)
f-RG-F	i-Ini-Hom (0.70), i-V0D-diag (0.50)
RG-E	i-Ini-E (0.40), i-CINi (0.40), i-Sh2-Hom (0.50)
Ini-F	i-RG-E (–1.00)
Ini-E	i-RG-F (–0.08)
V2a	i-V0V (1.00)
V2a-diag	i-V0V-diag (0.90)
IniV0V	i-RG-F (–0.07)
	Between left and right homologous circuits
V0D	c-RG-F (–0.07)
V0V	c-IniV0V (0.60)
V3	c-RG-F (0.03)
CINi	c-RG-F (–0.03)
	Between fore and hind homolateral circuits
f-Ini-Hom	h-RG-F (–0.01)
f-Sh2-Hom	h-RG-F (0.01)
h-Sh2-Hom	f-RG-F (0.125)
	Between diagonal circuits
f-V0D-diag	dh-RG-F (–0.075)
f-V0V-diag	dh-RG-F (0.02)
h-V0V-diag	df-RG-F (0.065)

i-, ipsilateral; c-, contralateral; f-, fore; h-, hind; d, diagonal. CINi, inhibitory commissural interneurons. Ini, regular inhibitory interneurons. RG-F, flexor center; RG-E, extensor center.

Computer simulations

Request a detailed protocol

The set of differential equations were solved using the odeint (Ahnert et al., 2011) implementation of the Runge-Kutta-Fehlberg 7–8 variable step-size solver of the boost C++ library (version 1.63.0). $I_{N o i s e}$ was calculated before the simulation with the Forward Euler method and for a fixed step size of 1 ms. The custom C++ code was compiled with clang 800.0.42.1 (Apple LLVM 8.0.0) for macOS 10.12.3. Simulation results were analyzed using Matlab 2016b. Source code and Matlab scripts for all simulations are available in ModelDB (McDougal et al., 2017) at http://modeldb.yale.edu/234101.

Data analysis

Request a detailed protocol

Activities of the flexor and extensor center were used to assess the model behavior. A RG was considered as being in flexion when $f (V)$ of its flexor center was greater or equal to 0.1, otherwise it was considered as being in extension. The locomotor period was defined as the duration between two consecutive onsets of the left hind flexor center; the frequency as the reciprocal of the period; extension and flexion phase durations as the duration between onset and offset of extension and flexion, respectively. Normalized phase differences were calculated as the delay between the onsets of the extension phases of a pair of rhythm generators divided by the locomotor period. Four normalized phase differences were calculated: hind left-right (right hind – left hind), fore left-right (right fore – left fore), homolateral (left fore – left hind) and diagonal (right fore – left hind). Gaits were operationally defined based on the phase differences (see Table 2).

Table 2

Operational definition of gaits.

https://doi.org/10.7554/eLife.31050.012

	Normalized phase differences
Gait	Left-right hind	Homolateral	Diagonal
Walk*	[0.25,0.75]	[0.1,0.4) and (0.6,0.9]	(0.1,0.4] and [0.6,0.9)
Trot	[0.25,0.75]	[0.25,0.75]	[0.0,0.1] and [0.9,1.0)
Gallop	(0.025,0.25] and [0.75,0.975)	[0.25,0.75]	[0.25,0.75]
Bound	[0.0,0.025] and [0.975,1.0)	[0.25,0.75]	[0.25,0.75]

*Classification of walk additionally required longer extension than flexion phase durations.

Analysis of model performance

Request a detailed protocol

To identify stable solutions of the model, bifurcation diagrams were built for the four normalized phase differences. To this end, $α$ was increased from 0.0 to 1.05 and then decreased back to 0.0 in 1000 equally spaced steps (step size of 0.00105). At each step, the simulation was run for 10 s and the final state was used as the initial condition for the next step. Circular statistics were evaluated for the last five locomotor cycles. At each step, the circular standard deviation of the normalized phase differences of the last five steps was smaller than 0.001. Thus, the mean phase differences can be regarded as the stable solutions. The sequential increase and decrease of $α$ was performed to uncover regions where multiple states are stable. An additional simulation run with stepwise change of $α$ in the opposite direction was initialized when a discrete change of phase differences between two $α$ -values occurred, to uncover potentially missed stable trajectories. Initial conditions were randomized and multiple runs were evaluated. $σ_{Noise}$ was set to 0.005 pA, making noisy currents several orders of magnitude smaller than any other current. The noise did not affect the model behavior other than to ensure that the system did not remain on an unstable trajectory.

To model step-to-step variability as presented in decabper, simulations with increased noisy currents ( $σ_{Noise} = 1.75$ pA) were performed. The free parameter $α$ was set to 0.3, 0.6, and 0.75. At each value, the simulation was run for 1000 s. Phase differences were calculated and partitioned into three bins (cf. Figure 7), depending on their appropriateness for trot and to ensure comparability with experimental data (Ruder et al., 2016). The counts per bin were then divided by the total number of locomotor cycles and multiplied by 100.

To model dynamics of gait transitions, extension phases and the instantaneous frequency are illustrate before and after an abrupt change in $α$ -value or additionally introduced drives. To this end, the model was initialized with a predefined $α$ -value and ten seconds of settling period were allowed before the abrupt parameter change. Parameters used are specified in the legend of Figure 8.

References

(2011) Odeint – solving ordinary differentialequations in C++
AIP Conference Proceedings 1389:1586–1589.

https://doi.org/10.1063/1.3637934
- Google Scholar
1. Alexander RM
2. Jayes AS
(1983) A dynamic similarity hypothesis for the gaits of quadrupedal mammals
Journal of Zoology 201:135–152.

https://doi.org/10.1111/j.1469-7998.1983.tb04266.x
- Google Scholar
(1987) Subpopulations and functions of long C3-C5 propriospinal neurones
Brain Research 404:395–400.

https://doi.org/10.1016/0006-8993(87)91402-8
- PubMed
- Google Scholar
(1990)
Characteristics of electrically induced locomotion in rat in vitro brain stem-spinal cord preparation

Journal of Neurophysiology 64:727–735.
- PubMed
- Google Scholar
(2017) State-Dependent Rhythmogenesis and Frequency Control in a Half-Center Locomotor CPG
Journal of Neurophysiology jn.00550.2017.

https://doi.org/10.1152/jn.00550.2017
- PubMed
- Google Scholar
(2001) Forelimb locomotor generators and quadrupedal locomotion in the neonatal rat
European Journal of Neuroscience 14:1727–1738.

https://doi.org/10.1046/j.0953-816x.2001.01794.x
- PubMed
- Google Scholar
(2003) Networks of inhibitory and excitatory commissural interneurons mediating crossed reticulospinal actions
European Journal of Neuroscience 18:2273–2284.

https://doi.org/10.1046/j.1460-9568.2003.02973.x
- PubMed
- Google Scholar
(2006) Differential projections of excitatory and inhibitory dorsal horn interneurons relaying information from group II muscle afferents in the cat spinal cord
Journal of Neuroscience 26:2871–2880.

https://doi.org/10.1523/JNEUROSCI.5172-05.2006
- PubMed
- Google Scholar
1. Bellardita C
2. Kiehn O
(2015) Phenotypic characterization of speed-associated gait changes in mice reveals modular organization of locomotor networks
Current Biology 25:1426–1436.

https://doi.org/10.1016/j.cub.2015.04.005
- PubMed
- Google Scholar
1. Biewener AA
(1990) Biomechanics of mammalian terrestrial locomotion
Science 250:1097–1103.

https://doi.org/10.1126/science.2251499
- PubMed
- Google Scholar
1. Bouyer LJ
2. Rossignol S
(2003) Contribution of cutaneous inputs from the hindpaw to the control of locomotion. II. Spinal cats
Journal of Neurophysiology 90:3640–3653.

https://doi.org/10.1152/jn.00497.2003
- PubMed
- Google Scholar
1. Britz O
2. Zhang J
3. Grossmann KS
4. Dyck J
5. Kim JC
6. Dymecki S
7. Gosgnach S
8. Goulding M
(2015) A genetically defined asymmetry underlies the inhibitory control of flexor–extensor locomotor movements
eLife 4:e04718.

https://doi.org/10.7554/eLife.04718
- Google Scholar
(2013) Ascending and descending propriospinal pathways between lumbar and cervical segments in the rat: evidence for a substantial ascending excitatory pathwaq
Neuroscience 240:83–97.

https://doi.org/10.1016/j.neuroscience.2013.02.039
- PubMed
- Google Scholar
1. Clarke KA
2. Still J
(1999) Gait analysis in the mouse
Physiology & Behavior 66:723–729.

https://doi.org/10.1016/S0031-9384(98)00343-6
- PubMed
- Google Scholar
(2008) Propriospinal neurons are sufficient for bulbospinal transmission of the locomotor command signal in the neonatal rat spinal cord
The Journal of Physiology 586:1623–1635.

https://doi.org/10.1113/jphysiol.2007.148361
- PubMed
- Google Scholar
(2010) Propriospinal transmission of the locomotor command signal in the neonatal rat
Annals of the New York Academy of Sciences 1198:42–53.

https://doi.org/10.1111/j.1749-6632.2009.05421.x
- PubMed
- Google Scholar
1. Crone SA
2. Quinlan KA
3. Zagoraiou L
4. Droho S
5. Restrepo CE
6. Lundfald L
7. Endo T
8. Setlak J
9. Jessell TM
10. Kiehn O
11. Sharma K
(2008) Genetic ablation of V2a ipsilateral interneurons disrupts left-right locomotor coordination in mammalian spinal cord
Neuron 60:70–83.

https://doi.org/10.1016/j.neuron.2008.08.009
- PubMed
- Google Scholar
(2009) In mice lacking V2a interneurons, gait depends on speed of locomotion
Journal of Neuroscience 29:7098–7109.

https://doi.org/10.1523/JNEUROSCI.1206-09.2009
- PubMed
- Google Scholar
1. Danner SM
2. Hofstoetter US
3. Freundl B
4. Binder H
5. Mayr W
6. Rattay F
7. Minassian K
(2015) Human spinal locomotor control is based on flexibly organized burst generators
Brain 138:577–588.

https://doi.org/10.1093/brain/awu372
- PubMed
- Google Scholar
(2016) Central control of interlimb coordination and speed-dependent gait expression in quadrupeds
The Journal of Physiology 594:6947–6967.

https://doi.org/10.1113/JP272787
- PubMed
- Google Scholar
1. Dietz V
2. Duysens J
(2000) Significance of load receptor input during locomotion: a review
Gait & Posture 11:102–110.

https://doi.org/10.1016/S0966-6362(99)00052-1
- PubMed
- Google Scholar
(2006) Long ascending propriospinal projections from lumbosacral to upper cervical spinal cord in the rat
Brain Research 1119:76–85.

https://doi.org/10.1016/j.brainres.2006.08.063
- PubMed
- Google Scholar
1. Duysens J
(1977)
Reflex control of locomotion as revealed by stimulation of cutaneous afferents in spontaneously walking premammillary cats

Journal of Neurophysiology 40:737–751.
- PubMed
- Google Scholar
(2013) The flexion synergy, mother of all synergies and father of new models of gait
Frontiers in Computational Neuroscience 7:14.

https://doi.org/10.3389/fncom.2013.00014
- PubMed
- Google Scholar
(2003) Both dorsal horn and lamina VIII interneurones contribute to crossed reflexes from feline group II muscle afferents
The Journal of Physiology 552:961–974.

https://doi.org/10.1113/jphysiol.2003.048009
- PubMed
- Google Scholar
1. Edinger L
(1896)
Verlag Von F.C.W. Vogel

Vorlesungen über den Bau der nervösen Centralorgane des Menschen und der Thiere, Verlag Von F.C.W. Vogel, 5 edn, Leipzig, Deutschland.
- Google Scholar
1. Ermentrout B
(1994) Reduction of conductance-based models with slow synapses to neural nets
Neural Computation 6:679–695.

https://doi.org/10.1162/neco.1994.6.4.679
- Google Scholar
1. Farley CT
2. Taylor CR
(1991) A mechanical trigger for the trot-gallop transition in horses
Science 253:306–308.

https://doi.org/10.1126/science.1857965
- PubMed
- Google Scholar
1. Flynn JR
2. Conn VL
3. Boyle KA
4. Hughes DI
5. Watanabe M
6. Velasquez T
7. Goulding MD
8. Callister RJ
9. Graham BA
(2017) Anatomical and molecular properties of long descending propriospinal neurons in mice
Frontiers in Neuroanatomy 11:5.

https://doi.org/10.3389/fnana.2017.00005
- PubMed
- Google Scholar
(1980) The locomotion of the low spinal cat. II. Interlimb coordination
Acta Physiologica Scandinavica 108:283–295.

https://doi.org/10.1111/j.1748-1716.1980.tb06534.x
- PubMed
- Google Scholar
1. Frigon A
2. Gossard JP
(2009) Asymmetric control of cycle period by the spinal locomotor rhythm generator in the adult cat
The Journal of Physiology 587:4617–4628.

https://doi.org/10.1113/jphysiol.2009.176669
- PubMed
- Google Scholar
(2013) Split-belt walking alters the relationship between locomotor phases and cycle duration across speeds in intact and chronic spinalized adult cats
Journal of Neuroscience 33:8559–8566.

https://doi.org/10.1523/JNEUROSCI.3931-12.2013
- PubMed
- Google Scholar
(2015) Modulation of forelimb and hindlimb muscle activity during quadrupedal tied-belt and split-belt locomotion in intact cats
Neuroscience 290:266–278.

https://doi.org/10.1016/j.neuroscience.2014.12.084
- PubMed
- Google Scholar
1. Frigon A
(2017) The neural control of interlimb coordination during mammalian locomotion
Journal of Neurophysiology 117:2224–2241.

https://doi.org/10.1152/jn.00978.2016
- PubMed
- Google Scholar
1. Fukuoka Y
2. Habu Y
3. Fukui T
(2015) A simple rule for quadrupedal gait generation determined by leg loading feedback: a modeling study
Scientific Reports 5:8169.

https://doi.org/10.1038/srep08169
- PubMed
- Google Scholar
(1999) Symmetry in locomotor central pattern generators and animal gaits
Nature 401:693–695.

https://doi.org/10.1038/44416
- PubMed
- Google Scholar
Book
1. Grillner S
(1981)
Control of locomotion in bipeds, tetrapods, and fish

In: Brookhart J. M, Mountcastle V. B, editors. Handbook of Physiology—The Nervous System II. American Physiological Society. John Wiley and Sons, Inc. pp. 1179–1236.
- Google Scholar
1. Grillner S
(1985) Neurobiological bases of rhythmic motor acts in vertebrates
Science 228:143–149.

https://doi.org/10.1126/science.3975635
- PubMed
- Google Scholar
1. Halbertsma JM
(1983)
The stride cycle of the cat: the modelling of locomotion by computerized analysis of automatic recordings

Acta Physiologica Scandinavica. Supplementum 521:1–75.
- PubMed
- Google Scholar
1. Heglund NC
2. Taylor CR
(1988)
Speed, stride frequency and energy cost per stride: how do they change with body size and gait?

The Journal of Experimental Biology 138:301–318.
- PubMed
- Google Scholar
(2004) Symmetrical and asymmetrical gaits in the mouse: patterns to increase velocity
Journal of Comparative Physiology A 190:895–906.

https://doi.org/10.1007/s00359-004-0545-0
- PubMed
- Google Scholar
1. Herbin M
2. Hackert R
3. Gasc JP
4. Renous S
(2007) Gait parameters of treadmill versus overground locomotion in mouse
Behavioural Brain Research 181:173–179.

https://doi.org/10.1016/j.bbr.2007.04.001
- PubMed
- Google Scholar
1. Hiebert GW
2. Pearson KG
(1999)
Contribution of sensory feedback to the generation of extensor activity during walking in the decerebrate Cat

Journal of Neurophysiology 81:758–770.
- PubMed
- Google Scholar
1. Hildebrand M
(1989) The quadrupedal gaits of vertebrates
BioScience 39:766–775.

https://doi.org/10.2307/1311182
- Google Scholar
1. Hoyt DF
2. Taylor CR
(1981) Gait and the energetics of locomotion in horses
Nature 292:239–240.

https://doi.org/10.1038/292239a0
- Google Scholar
(2013) Optogenetic dissection reveals multiple rhythmogenic modules underlying locomotion
PNAS 110:11589–11594.

https://doi.org/10.1073/pnas.1304365110
- PubMed
- Google Scholar
(2003)
Neuronal basis of crossed actions from the reticular formation on feline hindlimb motoneurons

Journal of Neuroscience 23:1867–1878.
- PubMed
- Google Scholar
(2005) Functional differentiation and organization of feline midlumbar commissural interneurones
The Journal of Physiology 565:645–658.

https://doi.org/10.1113/jphysiol.2005.083014
- PubMed
- Google Scholar
(2006) Neuronal relays in double crossed pathways between feline motor cortex and ipsilateral hindlimb motoneurones
The Journal of Physiology 575:527–541.

https://doi.org/10.1113/jphysiol.2006.112425
- PubMed
- Google Scholar
1. Jankowska E
(2008) Spinal interneuronal networks in the cat: elementary components
Brain Research Reviews 57:46–55.

https://doi.org/10.1016/j.brainresrev.2007.06.022
- PubMed
- Google Scholar
(2009) Commissural interneurons with input from group I and II muscle afferents in feline lumbar segments: neurotransmitters, projections and darget cells
The Journal of Physiology 587:401–418.

https://doi.org/10.1113/jphysiol.2008.159236
- PubMed
- Google Scholar
(2005) Propriospinal circuitry underlying interlimb coordination in mammalian quadrupedal locomotion
Journal of Neuroscience 25:6025–6035.

https://doi.org/10.1523/JNEUROSCI.0696-05.2005
- PubMed
- Google Scholar
1. Juvin L
2. Le Gal JP
3. Simmers J
4. Morin D
(2012) Cervicolumbar coordination in mammalian quadrupedal locomotion: role of spinal thoracic circuitry and limb sensory inputs
Journal of Neuroscience 32:953–965.

https://doi.org/10.1523/JNEUROSCI.4640-11.2012
- PubMed
- Google Scholar
1. Kato M
(1990) Chronically isolated lumbar half spinal cord generates locomotor activities in the ipsilateral hindlimb of the cat
Neuroscience Research 9:22–34.

https://doi.org/10.1016/0168-0102(90)90042-D
- PubMed
- Google Scholar
1. Kiehn O
(2016) Decoding the organization of spinal circuits that control locomotion
Nature Reviews Neuroscience 17:224–238.

https://doi.org/10.1038/nrn.2016.9
- PubMed
- Google Scholar
1. Lloyd DP
2. McIntyre AK
(1948)
Analysis of forelimb-hindlimb reflex activity in acutely decapitate cats

Journal of Neurophysiology 11:455–470.
- PubMed
- Google Scholar
1. Lemieux M
2. Josset N
3. Roussel M
4. Couraud S
5. Bretzner F
(2016) Speed-dependent modulation of the locomotor behavior in adult mice reveals attractor and transitional gaits
Frontiers in Neuroscience 10:42.

https://doi.org/10.3389/fnins.2016.00042
- PubMed
- Google Scholar
1. Lloyd DPC
(1942)
Mediation of descending long spinal reflex activity

Journal of Neurophysiology 5:435–458.
- Google Scholar
(2015) Primacy of flexor locomotor pattern revealed by ancestral reversion of motor neuron identity
Cell 162:338–350.

https://doi.org/10.1016/j.cell.2015.06.036
- PubMed
- Google Scholar
(1979) The location of spinal neurons with long descending axons (long descending propriospinal tract neurons) in the cat: a study with the horseradish peroxidase technique
The Journal of Comparative Neurology 184:63–79.

https://doi.org/10.1002/cne.901840105
- PubMed
- Google Scholar
1. Matsuyama K
2. Nakajima K
3. Mori F
4. Aoki M
5. Mori S
(2004) Lumbar commissural interneurons with reticulospinal inputs in the cat: morphology and discharge patterns during fictive locomotion
The Journal of Comparative Neurology 474:546–561.

https://doi.org/10.1002/cne.20131
- PubMed
- Google Scholar
1. McCrea DA
2. Rybak IA
(2007) Modeling the mammalian locomotor CPG: insights from mistakes and perturbations
Progress in Brain Research 165:235–253.

https://doi.org/10.1016/S0079-6123(06)65015-2
- PubMed
- Google Scholar
1. McCrea DA
2. Rybak IA
(2008) Deletions of rhythmic motoneuron activity during fictive locomotion and scratch provide clues to the organization of the mammalian central pattern generator
Brain Research Reviews 57:134–146.

https://doi.org/10.1016/j.brainresrev.2007.08.006
- Google Scholar
1. McDougal RA
2. Morse TM
3. Carnevale T
4. Marenco L
5. Wang R
6. Migliore M
7. Miller PL
8. Shepherd GM
9. Hines ML
(2017) Twenty years of ModelDB and beyond: building essential modeling tools for the future of neuroscience
Journal of Computational Neuroscience 42:1–10.

https://doi.org/10.1007/s10827-016-0623-7
- PubMed
- Google Scholar
(1985) Propriospinal fibers reaching the lumbar enlargement in the rat
Neuroscience Letters 58:257–261.

https://doi.org/10.1016/0304-3940(85)90174-0
- PubMed
- Google Scholar
(2016) Corticospinal and reticulospinal contacts on cervical commissural and long descending propriospinal neurons in the adult rat spinal cord; evidence for powerful reticulospinal connections
PLoS One 11:e0152094.

https://doi.org/10.1371/journal.pone.0152094
- PubMed
- Google Scholar
(2015) Mechanisms of left-right coordination in mammalian locomotor pattern generation circuits: a mathematical modeling view
PLoS Computational Biology 11:e1004270.

https://doi.org/10.1371/journal.pcbi.1004270
- PubMed
- Google Scholar
1. Nathan PW
2. Smith MC
(1959) Fasciculi proprii of the spinal cord in man
Brain 82:610–668.

https://doi.org/10.1093/brain/82.4.610
- PubMed
- Google Scholar
(1996) Vestibulospinal, reticulospinal and descending propriospinal nerve fibres in man
Brain 119:1809–1833.

https://doi.org/10.1093/brain/119.6.1809
- PubMed
- Google Scholar
1. Nicolopoulos-Stournaras S
2. Iles JF
(2009) Hindlimb muscle activity during locomotion in the rat (Rattus norvegicus) (Rodentia: Muridae)
Journal of Zoology 203:427–440.

https://doi.org/10.1111/j.1469-7998.1984.tb02342.x
- Google Scholar
1. Nishikawa K
2. Biewener AA
3. Aerts P
4. Ahn AN
5. Chiel HJ
6. Daley MA
7. Daniel TL
8. Full RJ
9. Hale ME
10. Hedrick TL
11. Lappin AK
12. Nichols TR
13. Quinn RD
14. Satterlie RA
15. Szymik B
(2007) Neuromechanics: an integrative approach for understanding motor control
Integrative and Comparative Biology 47:16–54.

https://doi.org/10.1093/icb/icm024
- PubMed
- Google Scholar
(1966)
[Locomotion induced by stimulation of the mesencephalon]

Doklady Akademii nauk SSSR 169:1223–1226.
- PubMed
- Google Scholar
(2009) Interlimb coordination in human crawling reveals similarities in development and neural control with quadrupeds
Journal of Neurophysiology 101:603–613.

https://doi.org/10.1152/jn.91125.2008
- PubMed
- Google Scholar
Book
1. Pearson KG
2. Duysens J
(1976)
Neural control of locomotion

In: Herman R. M, Grillner S, Stein P. S. G, Stuart D. G, editors. Function of Segmental Reflexes in the Control of Stepping in Cockroaches and Cats. Boston: Springer. pp. 519–537.
- Google Scholar
1. Pearson KG
(1995) Proprioceptive regulation of locomotion
Current Opinion in Neurobiology 5:786–791.

https://doi.org/10.1016/0959-4388(95)80107-3
- PubMed
- Google Scholar
1. Pearson KG
(2004) Generating the walking gait: role of sensory feedback
Progress in Brain Research 143:123–129.

https://doi.org/10.1016/S0079-6123(03)43012-4
- PubMed
- Google Scholar
(2015) Kinematic and gait similarities between crawling human infants and other quadruped mammals
Frontiers in Neurology 6:17.

https://doi.org/10.3389/fneur.2015.00017
- PubMed
- Google Scholar
(2006) Dynamic sensorimotor interactions in locomotion
Physiological Reviews 86:89–154.

https://doi.org/10.1152/physrev.00028.2005
- PubMed
- Google Scholar
(2016) Long-distance descending spinal neurons ensure quadrupedal locomotor stability
Neuron 92:1063–1078.

https://doi.org/10.1016/j.neuron.2016.10.032
- PubMed
- Google Scholar
(2006a) Modelling spinal circuitry involved in locomotor pattern generation: insights from deletions during fictive locomotion
The Journal of Physiology 577:617–639.

https://doi.org/10.1113/jphysiol.2006.118703
- PubMed
- Google Scholar
(2006b) Modelling spinal circuitry involved in locomotor pattern generation: insights from the effects of afferent stimulation
The Journal of Physiology 577:641–658.

https://doi.org/10.1113/jphysiol.2006.118711
- PubMed
- Google Scholar
(2013) Modelling genetic reorganization in the mouse spinal cord affecting left-right coordination during locomotion
The Journal of Physiology 591:5491–5508.

https://doi.org/10.1113/jphysiol.2013.261115
- PubMed
- Google Scholar
(2015) Organization of the mammalian locomotor CPG: Review of computational model and circuit architectures based on genetically identified spinal interneurons(1,2,3)
eNeuro 2:ENEURO.0069-15.2015.

https://doi.org/10.1523/ENEURO.0069-15.2015
- PubMed
- Google Scholar
(1990) A synergetic theory of quadrupedal gaits and gait transitions
Journal of Theoretical Biology 142:359–391.

https://doi.org/10.1016/S0022-5193(05)80558-2
- PubMed
- Google Scholar
(2015) Organization of left-right coordination of neuronal activity in the mammalian spinal cord: Insights from computational modelling
The Journal of Physiology 593:2403–2426.

https://doi.org/10.1113/JP270121
- PubMed
- Google Scholar
1. Shevtsova NA
2. Rybak IA
(2016) Organization of flexor-extensor interactions in the mammalian spinal cord: insights from computational modelling
The Journal of Physiology 594:6117–6131.

https://doi.org/10.1113/JP272437
- PubMed
- Google Scholar
(1966)
[Control of walking and running by means of electric stimulation of the midbrain]

Biofizika 11:659–666.
- PubMed
- Google Scholar
1. Shik ML
2. Orlovsky GN
(1976)
Neurophysiology of locomotor automatism

Physiological reviews 56:465–501.
- PubMed
- Google Scholar
(1979) Cells of origin of long descending propriospinal fibers connecting the spinal enlargements in cat and monkey determined by horseradish peroxidase and electrophysiological techniques
The Journal of Comparative Neurology 188:443–454.

https://doi.org/10.1002/cne.901880307
- PubMed
- Google Scholar
(1980) Responses of long descending propriospinal neurons to natural and electrical types of stimuli in cat
Brain Research 196:387–403.

https://doi.org/10.1016/0006-8993(80)90403-5
- PubMed
- Google Scholar
1. Skinner RD
2. Garcia-Rill E
(1984) The mesencephalic locomotor region (MLR) in the rat
Brain Research 323:385–389.

https://doi.org/10.1016/0006-8993(84)90319-6
- PubMed
- Google Scholar
(2011) Organization of functional synaptic connections between medullary reticulospinal neurons and lumbar descending commissural interneurons in the neonatal mouse
Journal of Neuroscience 31:4731–4742.

https://doi.org/10.1523/JNEUROSCI.5486-10.2011
- PubMed
- Google Scholar
1. Talpalar AE
2. Bouvier J
3. Borgius L
4. Fortin G
5. Pierani A
6. Kiehn O
(2013) Dual-mode operation of neuronal networks involved in left-right alternation
Nature 500:85–88.

https://doi.org/10.1038/nature12286
- PubMed
- Google Scholar
(2013) Coordination between the fore- and hindlimbs is bidirectional, asymmetrically organized, and flexible during quadrupedal locomotion in the intact adult cat
Neuroscience 240:13–26.

https://doi.org/10.1016/j.neuroscience.2013.02.028
- PubMed
- Google Scholar
(2011) Neurochemical excitation of propriospinal neurons facilitates locomotor command signal transmission in the lesioned spinal cord
Journal of Neurophysiology 105:2818–2829.

https://doi.org/10.1152/jn.00917.2010
- PubMed
- Google Scholar
1. Zhang J
2. Lanuza GM
3. Britz O
4. Wang Z
5. Siembab VC
6. Zhang Y
7. Velasquez T
8. Alvarez FJ
9. Frank E
10. Goulding M
(2014) V1 and v2b interneurons secure the alternating flexor-extensor motor activity mice require for limbed locomotion
Neuron 82:138–150.

https://doi.org/10.1016/j.neuron.2014.02.013
- PubMed
- Google Scholar
(2012) Neuronal activity in the isolated mouse spinal cord during spontaneous deletions in fictive locomotion: insights into locomotor central pattern generator organization
The Journal of Physiology 590:4735–4759.

https://doi.org/10.1113/jphysiol.2012.240895
- PubMed
- Google Scholar

Article and author information

Author details

Simon M Danner

Department of Neurobiology and Anatomy, Drexel University College of Medicine, Philadelphia, United States

Contribution
Conceptualization, Software, Formal analysis, Investigation, Visualization, Methodology, Writing—original draft, Writing—review and editing

For correspondence
simon.danner@gmail.com

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-4642-7064
Natalia A Shevtsova

Department of Neurobiology and Anatomy, Drexel University College of Medicine, Philadelphia, United States

Contribution
Investigation, Methodology, Writing—review and editing

Competing interests
No competing interests declared
Alain Frigon

Department of Pharmacology-Physiology, Université de Sherbrooke, Sherbrooke, Canada

Contribution
Conceptualization, Writing—review and editing

Competing interests
No competing interests declared
Ilya A Rybak

Department of Neurobiology and Anatomy, Drexel University College of Medicine, Philadelphia, United States

Contribution
Conceptualization, Supervision, Funding acquisition, Investigation, Visualization, Methodology, Project administration, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-3461-349X

Funding

National Institutes of Health (R01 NS081713)

Ilya A Rybak

National Institutes of Health (R01 NS090919)

Ilya A Rybak

National Institutes of Health (R01 NS095366)

Natalia A Shevtsova

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

2,971

views
476

downloads
103

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Article PDF

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Simon M Danner
Natalia A Shevtsova
Alain Frigon
Ilya A Rybak

(2017)

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

eLife 6:e31050.

https://doi.org/10.7554/eLife.31050

Share this article

Cite this article

Model concept.

Connections within the spinal cord.

Full model schematic.

Model performance.

Bifurcation diagrams of the intact model (A), after removal of all V0V neurons (B), of only V0V LPNs (C), and of all V0V and V0D neurons (D).

Model performance after deletion of all cervical-to-lumbar long propriospinal neurons (LPNs).

Model performance under application of noise before (intact) and after ablation of cervical-to-lumbar long projecting propriospinal neurons (LPNs).

Gait transitions with changes of brainstem drive (A–C) and independent of brainstem drive (G–I).

Schematic representation of control of limb coordination and gait by different local commissural interneurons (CINs) and long propriospinal neurons (LPNs).

Connection weights.

Operational definition of gaits.

Author details

Simon M Danner

Contribution

For correspondence

Competing interests

Natalia A Shevtsova

Contribution

Competing interests

Alain Frigon

Contribution

Competing interests

Ilya A Rybak

Contribution

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Further reading

Bifurcation diagrams of the intact model (A), after removal of all V0_V neurons (B), of only V0_V LPNs (C), and of all V0_V and V0_D neurons (D).