Abstract
The central pattern generator (CPG) architecture for rhythm generation remains partly elusive. We compare cat and frog locomotion results, where the component unrelated to pattern formation appears as a temporal grid, and traveling wave respectively. Frog spinal cord microstimulation with N-methyl-D-Aspartate (NMDA), a CPG activator, produced a limited set of force directions, sometimes tonic, but more often alternating between directions similar to the tonic forces. The tonic forces were topographically organized, and sites evoking rhythms with different force subsets were located close to the constituent tonic force regions. Thus CPGs consist of topographically organized modules. Modularity was also identified as a limited set of muscle synergies whose combinations reconstructed the EMGs. The cat CPG was investigated using proprioceptive inputs during fictive locomotion. Critical points identified both as abrupt transitions in the effect of phasic perturbations, and burst shape transitions, had biomechanical correlates in intact locomotion. During tonic proprioceptive perturbations, discrete shifts between these critical points explained the burst durations changes, and amplitude changes occurred at one of these points. Besides confirming CPG modularity, these results suggest a fixed temporal grid of anchoring points, to shift modules onsets and offsets. Frog locomotion, reconstructed with the NMDA synergies, showed a partially overlapping synergy activation sequence. Using the early synergy output evoked by NMDA at different spinal sites, revealed a rostrocaudal topographic organization, where each synergy is preferentially evoked from a few, albeit overlapping, cord regions. Comparing the locomotor synergy sequence with this topography suggests that a rostrocaudal traveling wave would activate the synergies in the proper sequence for locomotion. This output was reproduced in a two-layer model using this topography and a traveling wave. Together our results suggest two CPG components: modules, i.e., synergies; and temporal patterning, seen as a temporal grid in the cat, and a traveling wave in the frog. Animal and limb navigation have similarities. Research relating grid cells to the theta rhythm and on segmentation during navigation may relate to our temporal grid and traveling wave results. Winfree’s mathematical work, combining critical phases and a traveling wave, also appears important. We conclude suggesting tracing, and imaging experiments to investigate our CPG model.
Introduction
This review article will primarily be an account of how, in line with this Frontiers topic, work in a lower vertebrate (frog) spinal cord added insight to previous work investigating the central pattern generator (CPG) for cat locomotion, and how both approaches enrich each other. Although the focus of this review will be influenced by the fact that one of the authors was involved in research in both species, other relevant work will also be reviewed.
We will begin by presenting evidence that the CPG for cat locomotion conceptually consists of two components. We will then present how research in the frog led to insights into each of these two components. We will also review how broader perspectives coming from the hippocampal and medial entorhinal circuitry on one hand, and from mathematical work on phase singularities on the other hand, might further inform insights derived from comparing the results in cat and frog locomotion. We will conclude with a suggestion for future directions.
Results for Cat Locomotion
We begin with some insights into the CPG for cat locomotion, derived from a “black box” approach, using the effects of phasic and tonic proprioceptive inputs on the fictive locomotor output recorded from both forelimbs (Saltiel and Rossignol, 2004a,b). In the fictive preparation, there is no movement as the muscles are paralyzed with curare, and the CPG activity is recorded from muscle nerves (electroneurograms = ENGs). When phasic shoulder perturbations are applied to one forelimb at various times in the locomotor cycle, measured from the onset of the long head of triceps (TriLo, shoulder retractor and elbow extensor), abrupt transitions are seen in their effects, identifying critical points. Retractions monotonically shorten, then lengthen cycle duration (interval between consecutive TriLo onsets) compared to control, but at the flexor to extensor transition (labeled critical point D), they produce a bifurcation where the cycle is either shortened back to control duration, or remains lengthened (Figure 1A). The upper arm of the bifurcation in Figure 1B corresponds to the activation of a second flexor burst, while the contralateral extensor remains active (2:1 rhythm). Protractions during the extensor cycle at first prolong it, but at ~55% of the extensor phase (labeled critical point C), this effect abruptly changes to a shortening (Figure 2A). Conceptually, the plots of Figures 1, 2A are similar to studying phase response curves (Winfree, 1977). The methodological difference is that we plot the perturbed cycle duration against the time of the perturbation, whereas phase response curves plot the phase shift resulting from the change in cycle duration (positive/negative shift, or phase advance/delay, for a shortened/prolonged cycle), against the phase of the perturbation.
Figure 1
Figure 2
Protractions during the flexor phase shorten it, but at ~40% of the contralateral extensor phase (labeled critical point B), an abrupt transition occurs where protraction maximally advances the onset of the ipsilateral extensor to ~55% of the contralateral extensor phase, followed by a linear relationship (Figure 2B). Thus these phasic perturbations identified three critical points in the locomotor cycle. An interplay between these critical points is already seen, since the last result can be reformulated as saying that the effect phasic protractions at point B, is to advance point D to point C.
We found several other lines of evidence strengthening the identification of these critical points.
First, the effects of phasic retractions and protractions on the ipsilateral extensor cycle duration are markedly similar respectively to those produced by a brief hyperpolarizing or depolarizing stimulus on the inter-spike interval of a cyclically firing neuron (Winfree, 1987, his Figures 4.3 and 4.1). Like the effect of retraction on the cycle, a brief hyperpolarizing pulse first shortened, then lengthened the inter-spike interval, and near the end of the interval, produced a phase singularity in the work of Best and Winfree (Best, 1979). And like protraction, a brief depolarization first lengthened, then shortened the inter-spike interval, with a phase singularity in between. Thus it appears that a very basic result found in a single firing neuron is conserved at the level of an entire CPG. A recent study of human locomotion (Funato et al., 2016), using a novel method to obtain phase response curves to treadmill acceleration and deceleration, found the same two key points in the cycle as we did: one at touchdown (like our point D), and one at mid-stance (like our point C). A difference was that these points were not described as phase singularities, but as phases to which the cycle was preferentially reset (no change in phase when perturbation was done at that phase, phase advance and delay for perturbations prior to and after that phase respectively).
Second, some of the critical points correspond to clear-cut events in the morphology of the ENG bursts (Figure 3). Point D essentially corresponds to the transition from the ipsilateral flexor to the ipsilateral extensor phase. Point C corresponds to an abrupt upward inflection in the ipsilateral lateral head of triceps (TriLa, elbow extensor), and the onset of descent of the contralateral cleidobrachialis (ClB, shoulder protractor, and elbow flexor). This correspondence of critical points with ENG burst events further suggests that they are built in the CPG program, and that at least point C assists in interlimb coordination.
Figure 3
Third, the timing of these critical points in the fictive locomotor cycle closely correlates with important biomechanical events in the real locomotor cycle (Manter, 1938). For example, point C corresponds to an upward inflection in the ipsilateral vertical ground reaction force (GRF), a reversal from backwards to forwards in the ipsilateral horizontal GRF, and onset of contralateral elbow extension during swing (transition from F to E1 kinematic phase). Point B corresponds to when, after the initial shoulder retraction and elbow flexion necessary in the early swing phase to clear the ground, the limb movement changes to protraction. Thus not only are the critical points embedded in the CPG program, but they are at times that are biomechanically meaningful.
The most interesting result came from examining the effect of tonic proprioceptive perturbations on the structure of the fictive locomotor cycle (Saltiel and Rossignol, 2004b). First, we found that the interlimb coupling at point C (ipsilateral Trila inflection, contralateral ClB onset of descent) remained preserved during tonic shoulder protractions and retractions, and tonic elbow flexions and extensions. Second, when we synchronized on this preserved coupling to compare the control (limb pendent), and perturbed cycles (limb in a different tonic position), we observed that the reorganization in terms of burst duration was again based on the critical points. For example, tonic shoulder protraction/retraction prolonged/shortened ipsilateral TriLo and contralateral ClB only in their portions before point C. Overall the timing of burst events (e.g., offset, onset) remained unchanged in approximately the half-cycle from point C to the following point B with shoulder perturbations. Tonic elbow flexion instead prolonged the portion of the extensor burst after point C.
More specifically, we observed that shifts between critical points accounted for the changes in burst durations. For example, during tonic protraction, the onset of descent of ClB (point C) is advanced to the timing of point B in the control cycle; and similarly the onset of Trilo (point D) is advanced to the timing of point C in the control cycle. Thus we can summarize this by saying that tonic protraction changes burst durations by producing shifts between critical points: point C is shifted to C’ = B ; point D is shifted to D’ = C (Figure 4A). During tonic extensions and tonic retractions, ipsilateral ClB is prolonged; this happens by delaying the onset of descent of ClB (point C) to the timing of point D in the control cycle, i.e., point C is shifted to C’ = D (Figure 4B). Changes in burst amplitudes seem generally not to be restricted to segments delimited by the critical points. A notable exception is during tonic elbow extension, where the ipsilateral ClB remains of unchanged amplitude initially, and then increases at point B (Figure 4B).
Figure 4
These results suggest that when burst durations are changed by tonic proprioceptive inputs, something important actually remains unchanged, i.e., a temporal structure or grid, which allows the shift of a critical point to the timing of the earlier or later critical point in the sequence.
In other words, the cat locomotor CPG appears to be conceptually made up of two components, ENG burst segments delimited by critical points; and a temporal structure or grid determining where these critical points may be located in time.
Other researchers have suggested two components to the CPG (Lennard, 1985; Burke et al., 2001; Rybak et al., 2006). This is based on a different type of evidence, i.e., observations that spontaneous EMG deletions (Stein and Daniels-McQueen, 2004), or the effects of afferent stimulation, may or may not be associated with permanent resetting of the rhythm. Although formulations of rhythm-generation and pattern-formation in CPG models (Perret and Cabelguen, 1976; Burke et al., 2001; Lafreniere-Roula and McCrea, 2005; McCrea and Rybak, 2008; Zhong et al., 2012) do not involve critical points, they postulate, like our results, a network concerned with temporal organization, that constrains the expression of EMG patterns.
From a broader perspective, it is interesting that Cordo et al. (1993) found, when comparing human reaching and grasping to different heights, or with different hand orientations, that EMG averages only diverged at discrete periodic times in the movement. They concluded that a central timing oscillator seemed involved even in a non-cyclical movement sequence. More recently, Lisman (2014) reviewed evidence for discontinuous control of the motor system, in particular for action selection, by ~10 Hz oscillations.
Results in The Frog
Around the same time that this work (Saltiel and Rossignol, 2004a,b) was being done in the cat, a new approach was being developed at MIT. Inspired by the equilibrium-point hypothesis (Feldman, 1986; Bizzi et al., 1992), the question was whether circuits exist in the spinal cord, that specify limb movement in particular direction, to mediate commands from higher brain centers. Instead of a black-box approach, the method used was a direct electrical micro-stimulation of the intermediate region of the lumbar spinal cord in the spinalized frog. A modular organization was found, where only a limited set of force directions could be produced, with intermediate force directions obtained by co-stimulation of spinal cord sites (Bizzi et al., 1991).
In order to avoid stimulation of fibers of passage, chemical micro-stimulation was next done, using focal intra-spinal iontophoresis (Saltiel et al., 1998). The chemical that was found to work was N-methyl-D-Aspartate (NMDA). Only 24% of tested sites respond to NMDA, and with the parameters used, the spread of NMDA by the average time that forces began was estimated at 150–270 μm radius, which corresponds to less than 1% of the volume of a spinal cord segment.
With NMDA, a minority of the active spinal cord sites responded with a single tonic force, which could be rostral flexion, adduction, caudal extension, or lateral extension. The more common effect was to obtain a rhythm made up of an alternation between these forces, with some preferred combinations of these forces being more frequently observed. The tonic forces were topographically organized, and interestingly, the spinal cord rhythms were obtained from sites near the tonic regions producing their constituent forces, as well as from a common caudal region in the cord. It is not surprising that NMDA often produced a rhythm, since it is a very common agent used to activate CPGs (Wallén and Grillner, 1987; Li, 2011). The above results provide strong evidence that the modules identified in the frog spinal cord could represent building blocks of CPGs. The hypothesis of modularity of CPGs had already been formulated by other investigators (Grillner, 1981; Jordan, 1991).
Thus these modules identified in the frog could be related conceptually to one of the CPG components suggested for the cat, i.e., the production of units of muscle activations, that here would bring the limb in specific directions.
To explore this idea further, the next step was to analyze the EMG recordings produced by NMDA in the frog spinal cord. Looking carefully at EMG patterns, it was realized that these appear composed of smaller units which combine to produce the full pattern (Saltiel et al., 2001). For example, an EMG pattern could sometimes be seen to evolve over time, beginning with a subset of EMGs, to which another subset of EMGs was later added (Figure 5A). In other circumstances, the subsets of EMGs were produced in isolation. These smaller units were termed muscle synergies. EMG patterns evoked in natural behaviors also appeared to represent a combination of synergies (Tresch et al., 1999; d’Avella et al., 2003; Ting and Macpherson, 2004; Cheung et al., 2005).
Figure 5
In order to formally extract these muscle synergies suspected from visual inspection, a computational algorithm was developed (Tresch et al., 1999). Starting from a random set of initial synergies, the algorithm finds the coefficients of activation of synergies whose combinations best reconstruct the data. In the subsequent iterations, the algorithm alternatively updates the synergies, and the coefficients to maximize the EMG data reconstruction; typically a non-negative least-squares method is used. It was found that the combination of seven muscle synergies, labeled as A to G (Figure 5B), explains 91% of the variance of the NMDA-evoked EMGs recorded from 12 muscles (Saltiel et al., 2001).
It is possible to extend the above observations about NMDA evoking rhythms suggestive of CPG activation, by examining the EMG patterns. Thus we found that one of the commonly evoked rhythms, adduction-caudal extension-flexion to body, consists of 2 rhythm subtypes determined on the basis of EMGs (Saltiel et al., 2005). Subtype 1 consists of an adduction based on synergy E, and caudal extension based on synergies A + B; in subtype 2, the adduction is based on synergies E + C, and caudal extension on synergies A + D. Flexion onsets also differ between these 2 subtypes. Topographically, these 2 rhythms appear to originate from interleaved regions in the spinal cord encoding in rostrocaudal order, CE2, ADD1, ADD2 and CE1 (where CE = caudal extension, ADD = adduction and 1 and 2 refer to the two subtypes of linked caudal extensions and adductions). CE1 and CE2 caudal extensions appear to relate respectively to jumping and swimming extensions (Saltiel et al., 2005).
What about locomotion? Frogs are also capable of terrestrial quadrupedal walking (Roh et al., 2011). Here we had the intuition that NMDA-evoked rhythms comprising a lateral extension force, and in particular a commonly observed CE1-lateral force sequence would be the most relevant to locomotion. We reconstructed the locomotor EMG data with the NMDA synergies, and found close similarity between the locomotor synergy sequence, and the synergy sequence of the NMDA CE1-lateral force-flexion cycle (flexion onset without the C synergy). In both cases, an A-B-G-A-F-E-G synergy sequence is seen (Figure 6; Saltiel et al., 2016). Of note, the synergies are not activated in complete isolation in this sequence, but partly overlap; for example, during stance, synergies B and G begin together at the same time as the first synergy A peak, but synergy G outlasts them.
Figure 6
To understand how this sequence might come about, we then turned to the topography of NMDA synergies in the lumbar cord (Saltiel et al., 2016). This was determined by looking at the early EMG responses evoked at 110 spinal cord sites. Although there is considerable overlap, each of the seven synergies is preferentially activated from distinct regions of the spinal cord (Figure 7A), as seen in a pattern of peaks and troughs of activation along the rostrocaudal axis of the cord (Figure 7B). Other researchers using optical stimulation (Hägglund et al., 2013; Levine et al., 2014), or recordings (Tresch and Kiehn, 1999), have also found evidence for a rostrocaudal topography of motor patterns, although generally classified as flexor or extensor, rather than precisely identified synergies.
Figure 7
Interestingly, the temporal sequence order of synergies in locomotion, and the rostrocaudal spatial sequence order of synergies in the spinal cord clearly show similarities (compare Figures 6C, 7B), and this was quantitatively demonstrated as well (Figure 8; Saltiel et al., 2016). Therefore it appears that the synergy topography is such that a rostrocaudal traveling wave of activation would activate them in the proper sequence for locomotion. In our article, this was also investigated as a model, using as in the literature, a pattern-formation layer (PF), here encoding muscle synergies; and a “traveling wave” layer (TW), giving some perpendicular external inputs to the PF layer, in particular to nodes producing synergy A. This model reproduced the locomotor synergy sequence, and allowed to explore possible modifications by afferent inputs. It should be noted that the connectivity between synergy nodes within the PF layer in this model is already sufficient to give a propagating wave of synergy activation. However the external inputs from the TW layer to the PF layer are necessary to obtain the double-peaked synergy A activation during stance, contribute to determining stance onset, and participate in modifications by afferent inputs.
Figure 8
After finding this evidence for a traveling wave in frog locomotion, we noticed that a traveling wave at the interneuronal level had already been reported in the fictive scratch cycle of the cat (Cuellar et al., 2009), as shown by recording of cord dorsum potentials, even after ventral horn removal; and by interneuronal recordings. Prior to this, among limbed vertebrates, traveling waves of motoneuron activation had been reported in the newt (Delvolvé et al., 1997), and the rodent (Bonnot et al., 2002; Cazalets, 2005), and suggested in the cat (Yakovenko et al., 2002), and human (Ivanenko et al., 2006). In vertebrates without limbs, a traveling wave occurs during swimming in the lamprey (Cohen et al., 1992), and zebrafish (Wiggin et al., 2012). Among invertebrates with limbs, a traveling wave occurs in the crustacean swimmeret system (Davis, 1969; Heitler, 1980). More recently, a traveling wave of interneuronal activation has been reported by examining oscillator ensembles in the Aplysia pedal ganglion during locomotion (Bruno et al., 2015). These interneuronal ensembles both show rotational dynamics as seen in a principal component analysis, and physical rotation of activity, demonstrated as an elliptical trajectory when physically mapped on the ganglion. It has also been shown that a two-dimensional spinal cord model could generate traveling waves accounting for the various quadruped gaits (Kaske et al., 2003).
What could the neural substrate be for a rostrocaudal traveling wave? The first idea would be a descending excitatory pathway, since it is generally accepted that an excitatory mechanism is at the core of the locomotor CPG (Grillner and El Manira, 2015). In the last 10–15 years, genetic methods have allowed to identify locomotor functions for several of the V0 to V3 interneurons defined by their transcription factors (Kiehn, 2016), and to suggest what their connectivity might be (Rybak et al., 2015). Among these, the excitatory V2a interneurons, which exist in the mouse, zebrafish, Xenopus tadpole and lamprey, have several features that might make them good candidates to participate in a rostrocaudal wave of excitation within the PF layer. First, V2a interneurons of type II in the mouse, which become silent during non-resetting deletions, are believed to be part of the PF layer and to recruit motoneurons (Zhong et al., 2012; Rybak et al., 2015). Second, V2a interneurons appear from optogenetic experiments to also be connected with other V2a interneurons, and their selective activation is sufficient to produce the rostrocaudal wave of swimming in the zebrafish (Ljunggren et al., 2014). Third, their axons are ipsilaterally descending, at least in the zebrafish (Kimura et al., 2006; Ljunggren et al., 2014). One may thus ask whether V2a type II interneurons might be premotor interneurons projecting both to motoneurons to encode synergies, and to other premotor neurons further caudally, to participate in a rostrocaudal wave of excitation producing the locomotor synergy sequence as in our model (Saltiel et al., 2016).
In the mammalian CPG literature, genetic studies of intralimb coordination have focused on flexor-extensor alternation, probably for two reasons. First, the L2 and L5 ventral roots show alternating activity during in vitro neonatal rat locomotor activity, with L2 in phase with simple flexors, and L5 with extensors, although the EMG pattern is in fact more complex (Kiehn and Kjaerulff, 1996). Second, whereas a longitudinally traveling wave is a well-established concept for limbless vertebrates such as the lamprey or zebrafish, it is only more recently that it has been proposed as a mechanism for locomotion of vertebrates with limbs (see above). Nevertheless, it is interesting that the V1 and V2b inhibitory interneurons, which are key for flexor-extensor alternation in the mouse, are respectively homologous to an ascending inhibitory interneuron (aIN in Xenopus tadpole, and CiA in zebrafish), and to a descending inhibitory interneuron (VeLD; Zhang et al., 2014). Further, V1-derived (ascending) inhibition promoted limb extension, and V2b-derived (descending) inhibition promoted limb flexion in the mouse (Britz et al., 2015). Given the greater rostral representation of flexor synergies, and generally greater caudal representation of extensor synergies in the frog spinal cord (Saltiel et al., 2016), these are the results we would expect if similar longitudinal inhibitory pathways also operate in the adult frog. The advantage of the postulated rostrocaudal traveling wave is that it does not reduce the activation of the seven synergies in the frog locomotor sequence, or the complex EMG patterns in other species, to just flexor-extensor alternation. At the same time, mechanisms for flexor-extensor alternation, and for a traveling wave are not incompatible. In particular, the ascending inhibitory pathway is active during swimming in Xenopus tadpole, and may be a candidate for regulating intersegmental delay (Li et al., 2004; Zhang et al., 2014). The combination of descending excitation, as discussed above (see also Cowley et al., 2010), and ascending inhibition, would be expected to favor a rostrocaudally traveling wave, as pointed out by Kimura et al. (2006). Such combined descending excitation and ascending inhibition between the more rostral flexor, and more caudal extensor synergies, could also relate to the observations of tighter phase-locking at the extension to flexion (E to F) than at the F to E transition (Boothe et al., 2013). These authors suggested a tighter phase-locking would be expected for a transition resulting from an escape or release from inhibition (here postulated from E to F).
In summary, examination of the genetic literature on CPG interneuronal components in different species, suggests hypotheses about how a traveling wave mechanism for locomotion in the limbed vertebrate might be implemented.
There is evidence that traveling waves may relate to the rhythm-generator (RG) component in the two-layer model of the CPG, and thus be involved in time-keeping. In resetting deletions of the motor output of fictive scratch in the cat, where the cycle period is not maintained, the cord dorsum potential rostrocaudally traveling wave simultaneously disappears. By contrast, the wave is maintained during non-resetting deletions. Both observations have been reproduced in a longitudinal two-layer CPG model with a traveling wave (Pérez et al., 2009).
When thinking back about our results in cat locomotion, we had already related the synergy modules in the frog, to the segments of ENG burst activation in the cat. Now with this new result of a traveling wave in frog locomotion, it seems natural to ask whether it could relate in some way to our other result in cat locomotion, i.e., the temporal grid or structure.
It is not straightforward to answer this question. However we would like to suggest that perhaps one way to approach it is to think about parallels with the literature on the hippocampus/medial entorhinal cortex (mEC) circuitry.
Possible Insights from Some Similarities to Hippocampus/Entorhinal Circuitry
One may see elements of similarity between spatial navigation of a locomoting animal, controlled by the hippocampus and mEC, and the navigation of the limbs in space during locomotion, under control of the CPG. One may speak of limb navigation during the swing phase for example because there are constraints such as clearing the ground at take-off, or coordinating different joints during limb protraction, and then during landing. In navigation of the whole animal, place cell activity is determined by both self-motion and environmental inputs. During locomotion, the limb EMGs are both determined by a central program and by afferent input.
This idea of similarity is also supported by comparing what happens to place cell activity when a rat linear trajectory is shortened by putting the start-box closer to the goal location, and what happens in cat fictive locomotion during tonic forelimb protraction which shortens the swing phase. In the former case, one may see compared to the control trajectory, an abrupt shift from place cells active near the start-box to place cells active farther along the trajectory, with disappearance of place cells for intermediate locations, even though these locations are still being traversed (Gothard et al., 1996, their Figure 8). In the latter case, one sees a disappearance of the intermediate segment of ClB (shoulder protractor and elbow flexor) electroneurogram (segment B-C in Figure 4A), corresponding to an abrupt shift directly of critical point C to B, when comparing to the control locomotor output (Saltiel and Rossignol, 2004b).
We already mentioned that we could approximately relate our results in frog and cat locomotion by considering synergy modules in the frog equivalent to the ENG burst segments in the cat. This suggests that the other component for locomotor generation in the frog, i.e., a traveling wave, might be equivalent to the other locomotor component in the cat, i.e., a temporal grid. Perhaps the analogy with the hippocampus and mEC can give glimpses into how a traveling wave could relate to a temporal grid. Grid cells in the mEC are periodically activated during locomotion. There is evidence that the theta rhythm, produced by the medial septum with a frequency modulated by locomotor speed, may be crucial to obtain grid cell spatial periodicity (Brandon et al., 2011; Koenig et al., 2011; Hinman et al., 2016; but see Yartsev et al., 2011; Pastoll et al., 2013). Interestingly this theta oscillation appears organized as a traveling wave in the hippocampus, such that there are many theta inputs from the medial septum which just differ by their phases (Patel et al., 2012; Zhang and Jacobs, 2015). In terms of how a traveling theta wave may result in a grid structure, models suggest that the distributed phases of theta oscillatory input to different grid cells gate which ones respond to other shared periodic inputs, since these can only activate grid cells in a specific phase range. Because these inputs would slightly differ from theta in frequency, the grid cell subset which is in the appropriate phase for activation by the shared input cyclically changes to give successive firing of different grid cells and grid cell periodicity (Hasselmo and Shay, 2014; Shay et al., 2016). This idea of conjunction of inputs could relate to a two-layer model of the locomotor CPG, where periodic recruitment of a synergy sometimes depends on the conjunction of inputs from within a PF layer, with a discrete input from a separate TW layer (Saltiel et al., 2016).
A different suggested role for the theta rhythm in the mEC is to give rise to nested gamma oscillations, which appear regulated tightly enough in time, to be described as clock-like (Pastoll et al., 2013). Thus, although the time scales are different, an analogy could be suggested between the locomotor traveling wave and the temporal grid we respectively found in the frog and the cat, and the theta traveling wave and nested gamma clock-like oscillations in the mEC. In locomotion, the temporal grid of critical points subdivides the cycle in different segments; and critical points allow adjustments of these segments when afferent inputs indicate for example that the actual limb position is ahead or behind that expected from CPG activity. Some of these adjustments can be seen as advancing or delaying the CPG program to a later or earlier segment (Figure 4). In the hippocampus, the gamma cycles nested within the theta cycle may also help to represent the environment in segments (Gupta et al., 2012; Pfeiffer and Foster, 2015). And a change in the gamma frequency (slow vs. fast gamma) may promote a change in focus of the segmentation during the theta cycle, from looking ahead/prediction to encoding of actual location or experience, or to looking behind, when navigating or during episodic memory (Gupta et al., 2012; Bieri et al., 2014; Zheng et al., 2016).
In summary, research on the hippocampus/mEC system appears to give an interesting perspective to our results on the locomotor CPG, and suggests some features of generality. Tentatively, this may include a traveling wave giving rise to a grid structure; and using such a grid to allow flexible, yet discrete, adjustments in the segmentation.
Relationship to Arthur Winfree’s Work on Phase Singularities
In his book, Winfree (1987) first considered a single oscillating cell (squid giant axon, or the sinoatrial node), that was space-clamped, i.e., only time was a consideration. Testing the effect of a brief depolarizing or hyperpolarizing stimulus applied at different times in the inter-spike interval, there is a critical phase and stimulus magnitude where the timing of the post-stimulus spike is no longer predictable. These phase singularities were found both theoretically, e.g., using Hodgkin-Huxley equations, and experimentally. They occurred at a similar time in the squid axon and sinoatrial cycles, at ~38% and ~88% (half-a-cycle apart) of the inter-spike cycle, for the depolarizing and hyperpolarizing stimuli respectively (his Figures 4.3 and 4.1; Best, 1979). The first way that we can relate our work to these results is that we also identified critical points for the effect of phasic protractions and retractions on the Triceps ENG fictive locomotor cycle in the cat (time 0 being the onset of Triceps), and that these occurred at ~38% and ~89% of the cycle, respectively (Figures 1A, 2A). Thus it appears that a very basic result found in a single firing neuron is conserved at the level of an entire CPG.
Winfree (1983) next had the idea that he could easily extend his results to include the consideration of space as well as time. Oscillators were positioned as neighbors inside a rectangular space. By using a traveling wave of activation moving from right to left, a timing gradient was set up along the horizontal dimension of the rectangle. Thus in the space between two consecutive moving vertical wavefronts located some distance apart inside the rectangle, the array of oscillators covered the entire phase spectrum from 0 to 1 along the x-axis of the rectangle. Then applying a stimulus at different positions along the horizontal axis of the rectangle was equivalent to applying stimuli at different phases in time in the previous time-only paradigm. For a depolarizing stimulus, Winfree’s simulations predicted a pair of phase singularities located at the two possible intersections inside the rectangle between a vertical line at the critical phase for that stimulus, and a circle of critical stimulus magnitude from a centrally-applied radially-decaying stimulus applied at that phase (Winfree, 1987, his Figures 6.1 and 6.2). Similar predictions were made for a hyperpolarizing stimulus (his Figures 7.13 and 7.14). These predictions were confirmed in experiments done on the beating heart (Shibata et al., 1988; Winfree, 1989), and also relate to results in oscillating chemical solutions (Zaikin and Zhabotinsky, 1970; Winfree, 1972). Thus the second way that we can relate our work to Winfree results, is that he combined in the same analysis, a phase gradient created by a traveling wave, such as we found in the frog, and critical points at specific phases, such as we found in the cat, to understand where phase singularities will arise in space.
In Winfree’s simulations and related experiments, the phase singularities are the sites of rotors (pivots of a rotating wave) which organize timing in space. For example, Figures 9A,B (adapted from Winfree, 1987, Figures 6.1 and 6.2) are snapshots of the isochrons (curves connecting sites where phase is the same) immediately before and after a depolarizing stimulus at the critical phase. The cycle between the two wavefronts is divided in 12 intervals, and isochrons indicate regions of space that are synchronous in time. Wavefront isochrons are labeled 0 = 12, and isochrons labeled 1, 2, … indicate that the wavefront was there 1, 2, … units of time previously. The rotors have arisen from a stimulus given in Figure 9A, between isochrons 4 and 5 on the x-axis, at ~38% of the cycle between the wavefronts at the left and right, i.e., at the critical phase for a depolarizing stimulus. We notice that between the two rotors, the isochrons that were 4–5 before the stimulus (Figure 9A), are now 10–11 after the stimulus (Figure 9B), while the other regions do not show much change in their isochrons after the stimulus. If we imagine Figure 9B to be the spinal cord rostrocaudally oriented from right to left, with interneurons located in the space between the two rotors encoding a particular synergy, these simulations suggest a way for a stimulus to change the timing of that synergy. This would occur without affecting the timing of other components of the spinal cord output, encoded by other regions where little or no change occurs in isochrons after the stimulus. This kind of outcome is generally similar to the results in the cat where the timing of only some burst components is modified by afferent proprioceptive inputs. And it relates to the concept of two-layer models of CPG, where one layer maintains a temporal structure, while burst components encoded by another layer can be shifted in time, yet with the overall temporal structure remaining intact.
Figure 9
Finally, it is interesting that the critical phase for a hyperpolarizing stimulus at ~88% of the cycle is between isochrons 10–11 (Figure 9C). Perhaps the change as a result of a depolarizing stimulus in Figures 9A,B, from isochrons 4–5 (critical phase for a depolarizing stimulus) to 10–11 (critical phase for a hyperpolarizing stimulus) is an example of shifts between critical points. Winfree alludes to that observation in the legend of Figure 7.14 of his 1987 book. It remains unclear however how to explain the larger set of critical points that we found in the cat.
In summary, while it has long been suggested to view the CPG as coupled oscillators, it seems that useful insights might specifically come from theoretical considerations as in the ground-breaking work of Arthur Winfree, who emphasized phase singularities as a key to understand the organization of timing in space. This approach starts with the same two ingredients as our experimental evidence in locomotion for: (1) critical points, which we found precisely at the same times as reported by Winfree; and (2) a rostrocaudally traveling wave that would create a phase gradient in space. These two elements are then combined to localize phase singularities in the physical space of the spinal cord, with implications for altering the timing of output components such as synergies, while preserving the overall temporal structure of the cycle. Compared to other systems studied in this way (myocardium, chemical oscillations), the spinal cord and locomotion have the added features that the oscillators encode different synergies, with a specific topography in the spinal cord space. We could thus envision that it might be possible to characterize what the critical loci indicated by the vertical dashed lines in Figures 9A,C may be encoding. And expressing the new timing of the critical locus after the perturbation in terms of where the wavefront would have to be located in the unperturbed cord to result in that timing, might allow comparison to the critical point shifts observed when superimposing the perturbed and control locomotor cycles in the cat (Figure 4).
Future Directions
Given that the frog research has shed some light on cat locomotion, it would appear useful to further investigate in the frog the spinal cord circuitry that appears to underlie the locomotor synergy sequence, using our suggested two-layer model (Saltiel et al., 2016) to guide experiments. There are now several available techniques for retrograde tracing that can be limited to monosynaptic spread (Wickersham et al., 2007), or that can be made selective by using the Cre-loxP system. A good point of entry would be the synergy A (key muscle = Gastrocnemius) circuitry, since it is at the intersection of the PF layer and TW layer in our model.
For example, the vesicular stomatitis virus (VSV) pseudotyped with the rabies virus envelope, and encoding a fluorescent protein, can be used for retrograde tracing in a wide variety of species (Mundell et al., 2015). When the VSV-G (glycoprotein) gene is deleted from the virus (VSV-∆G), and the rabies virus G gene (RV-G) is supplied separately as trans, one expects a purely monosynaptic retrograde labeling (Beier et al., 2013). Thus their combined injections in GA muscle should selectively label synergy A premotor interneurons. Repeating an application of the RV-G gene at a spinal cord site evoking synergy A should again complement VSV-∆G, this time in synergy A premotor interneurons, resulting in labeling of their monosynaptic inputs. According to our model, this might include for example synergy G premotor interneurons in the PF layer, and interneurons from the TW layer. Coupling this experiment with a different fluorescent protein to label synergy G premotor interneurons could help in distinguishing between these two inputs, if indeed they correspond to different populations.
Another possible method for retrograde tracing is to use the wheat germ agglutinin (WGA) gene in an adeno-associated virus (AAV) vector to carry the Cre gene from GA muscle to synergy A premotor neurons. WGA labeling of premotor neurons from muscle injections has been shown (Harrison et al., 1986; Minic et al., 2016), and WGA tracing has been used in the frog (Levy et al., 2015). Then applying at a spinal cord site evoking synergy A, a Cre-dependent virus with a double-floxed fused WGA-fluorescent protein (Sugita and Shiba, 2005) or TTC (tetanus toxin fragment C)-fluorescent protein (Maskos et al., 2002) should retrogradely label inputs to synergy A premotor neurons.
Indeed, using anatomical techniques, evidence has been found for connectivity between spinal premotor neurons projecting to different motoneurons (Levine et al., 2014), but this has not yet been related to specific synergies. An example with functional implications is the demonstration of unidirectional connectivity between brainstem premotor neurons responsible for successive phases (pharyngeal and esophageal) of swallowing (Broussard et al., 1998). This would represent connectivity within a PF layer. More recently, evidence for different regional distributions of putative RG and PF interneurons has been reported (Griener et al., 2013).
It would also be interesting to search for more evidence for a traveling wave by using a method based on the CaMPARI genetically encoded calcium indicator that would be delivered widely to the spinal cord, using a virus selectively recovered for widespread transduction ability (Deverman et al., 2016). In this method, active neurons when exposed to a 405 nm light permanently photoconvert from a green to a red color. Shining the light at specific phases of the locomotor cycle would allow to see which neurons are active in that phase; and whether there is a rostrocaudal traveling wave for progressively later phases. Progress is being made in getting the light into greater depth up to 1 mm (Flytzanis et al., 2015).
In conclusion, the results in frog locomotion suggest specific ways in which the CPG circuitry could be anatomically investigated. At the same time, the dual perspective from the results in cat and frog locomotion suggests that theoretical approaches, inspired from other fields and disciplines, might be worthwhile to elucidate how critical points, and a traveling wave may be combined into a CPG architecture. Further experimental data such as the anatomical investigations above, may also provide other clues helpful to the theoretical approach.
Statements
Author contributions
All authors listed have made substantial, direct and intellectual contribution to the work. PS wrote the manuscript, and MCT, KW, AA and EB revised the manuscript.
Acknowledgments
This work was supported by National Institutes of Health (NIH) Grant NS09343, and Office of Naval Research (ONR) Grant: N00014/90/J/1946 to EB; a Canadian Medical Research Council Fellowship to PS; a Howard Hughes predoctoral fellowship to MCT; and the Swiss National Science Foundation to KW. We gratefully acknowledge the guidance of Serge Rossignol, in whose laboratory the work in the cat was performed; Margo Cantor and Sylvester Szczepanowski for technical support; and Charlotte Potak for administrative assistance.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
BeierK. T.SaundersA. B.OldenburgI. A.SabatiniB. L.CepkoC. L. (2013). Vesicular stomatitis virus with the rabies virus glycoprotein directs retrograde transsynaptic transport among neurons in vivo. Front. Neural Circuits7:11. 10.3389/fncir.2013.00011
2
BestE. N. (1979). Null space in the Hodgkin-Huxley equations: a critical test. Biophys. J.27, 87–104. 10.1016/s0006-3495(79)85204-2
3
BieriK. W.BobbittK. N.ColginL. L. (2014). Slow and fast γ rhythms coordinate different spatial coding modes in hippocampal place cells. Neuron82, 670–681. 10.1016/j.neuron.2014.03.013
4
BizziE.HoganN.Mussa-IvaldiF. A.GiszterS. F. (1992). Does the nervous system use equilibrium-point control to guide single and multiple joint movements?Behav. Brain Sci.15, 603–613. 10.1017/S0140525X00072538
5
BizziE.Mussa-IvaldiF. A.GiszterS. F. (1991). Computations underlying the execution of movement: a biological perspective. Science253, 287–291. 10.1126/science.1857964
6
BonnotA.WhelanP. J.MentisG. Z.O’DonovanM. J. (2002). Spatiotemporal pattern of motoneuron activation in the rostral lumbar and the sacral segments during locomotor-like activity in the neonatal mouse spinal cord. J. Neurosci.22:RC203.
7
BootheD. L.CohenA. H.TroyerT. (2013). Phase locking asymmetries at flexor-extensor transitions during fictive locomotion. PLoS One8:e64421. 10.1371/journal.pone.0064421
8
BrandonM. P.BogaardA. R.LibbyC. P.ConnerneyM. A.GuptaK.HasselmoM. E. (2011). Reduction of theta rhythm dissociates grid cell spatial periodicity from directional tuning. Science332, 595–599. 10.1126/science.1201652
9
BritzO.ZhangJ.GrossmannK. S.DyckJ.KimJ. C.DymeckiS.et al. (2015). A genetically defined asymmetry underlies the inhibitory control of flexor-extensor locomotor movements. Elife4:e04718. 10.7554/eLife.13038
10
BroussardD. L.LynnR. B.WiednerE. B.AltschulerS. M. (1998). Solitary premotor neuron projections to the rat esophagus and pharynx: implications for control of swallowing. Gastroenterology114, 1268–1275. 10.1016/s0016-5085(98)70433-0
11
BrunoA. M.FrostW. N.HumphriesM. D. (2015). Modular deconstruction revels the dynamical and physical building blocks of a locomotion motor program. Neuron86, 304–318. 10.1016/j.neuron.2015.03.005
12
BurkeR. E.DegtyarenkoA. M.SimonE. S. (2001). Patterns of locomotor drive to motoneurons and last-order interneurons: clues to the structure of the CPG. J. Neurophysiol.86, 447–462.
13
CazaletsJ. R. (2005). Metachronal propagation of motoneurone burst activation in isolated spinal cord of newborn rat. J. Physiol.568, 583–597. 10.1113/jphysiol.2005.086850
14
CheungV. C.d’AvellaA.TreschM. C.BizziE. (2005). Central and sensory contributions to the activation and organization of muscle synergies during natural motor behaviors. J. Neurosci.25, 6419–6434. 10.1523/JNEUROSCI.4904-04.2005
15
CohenA. H.ErmentroutG. B.KiernelT.KopellN.SigvardtK. A.WilliamsT. L. (1992). Modelling of intersegmental coordination in the lamprey central pattern generator for locomotion. Trends Neurosci.15, 434–438. 10.1016/0166-2236(92)90006-t
16
CordoP.SchieppatiM.BevanL.CarltonL. G.CarltonM. J. (1993). Central and peripheral coordination in movement sequences. Psychol. Res.55, 124–130. 10.1007/bf00419643
17
CowleyK. C.ZaporozhetsE.SchmidtB. J. (2010). Propriospinal transmission of the locomotor command signal in the neonatal rat. Ann. N Y Acad. Sci.1198, 42–53. 10.1111/j.1749-6632.2009.05421.x
18
CuellarC. A.TapiaJ. A.JuárezV.QuevedoJ.LinaresP.MartínezL.et al. (2009). Propagation of sinusoidal electrical waves along the spinal cord during a fictive motor task. J. Neurosci.29, 798–810. 10.1523/JNEUROSCI.3408-08.2009
19
d’AvellaA.SaltielP.BizziE. (2003). Combinations of muscle synergies in the construction of a natural motor behavior. Nat. Neurosci.6, 300–308. 10.1038/nn1010
20
DavisW. J. (1969). The neural control of swimmeret beating in the lobster. J. Exp. Biol.50, 99–117.
21
DelvolvéI.BemT.CabelguenJ. M. (1997). Epaxial and limb muscle activity during swimming and terrestrial stepping in the adult newt, Pleurodeles waltl. J. Neurophysiol.78, 638–650.
22
DevermanB. E.PravdoP. L.SimpsonB. P.KumarS. R.ChanK. Y.BanerjeeA.et al. (2016). Cre-dependent selection yields AAV variants for widespread gene transfer to the adult brain. Nat. Biotechnol.34, 204–209. 10.1038/nbt.3440
23
FeldmanA. G. (1986). Once more on the equilibrium-point hypothesis (lambda model) for motor control. J. Mot. Behav.18, 17–54. 10.1080/00222895.1986.10735369
24
FlytzanisN.ChallisC.LoogerL.SchreiterE.GradinaruV. (2015). Deep brain mapping of functional connectivity in intact circuits via CaMPARI and PACT tissue clearing. Soc. Neurosci. Abstr.736.08.
25
FunatoT.YamamotoY.AoiS.ImaiT.AoyagiT.TomitaN.et al. (2016). Evaluation of the phase-dependent rhythm control of human walking using phase response curves. PLoS Comput. Biol.12:e1004950. 10.1371/journal.pcbi.1004950
26
GothardK. M.SkaggsW. E.NcNaughtonB. L. (1996). Dynamics of mismatch correction in the hippocampal ensemble code for space: interaction between path integration and environmental cues. J. Neurosci.16, 8027–8040.
27
GrienerA.DyckJ.GosgnachS. (2013). Regional distribution of putative rhythm-generating and pattern-forming components of the mammalian locomotor CPG. Neuroscience250, 644–650. 10.1016/j.neuroscience.2013.07.070
28
GrillnerS. (1981). “Control of locomotion in bipods, tetrapods, and fish,” in Handbook of Physiology. The Nervous System. Motor Control, (Vol. 2) ed. BrooksV. B. (Bethesda, MD: Williams & Wilkins), 1179–1236.
29
GrillnerS.El ManiraA. (2015). The intrinsic operation of the networks that make us locomote. Curr. Opin. Neurobiol.31, 244–249. 10.1016/j.conb.2015.01.003
30
GuptaA. S.van der MeerM. A. A.TouretzkyD. S.RedishD. A. (2012). Segmentation of spatial experience by hippocampal theta sequences. Nat. Neurosci.15, 1032–1039. 10.1038/nn.3138
31
HägglundM.DoughertyK. J.BorgiusL.ItoharaS.IwasatoT.KiehnO. (2013). Optogenetic dissection reveals multiple rhythmogenic modules underlying locomotion. Proc. Natl. Acad. Sci. U S A110, 11589–11594. 10.1073/pnas.1304365110
32
HarrisonP. J.JankowskaE.ZytnickiD. (1986). Lamina VIII interneurons interposed in crossed reflex pathways in the cat. J. Physiol.371, 147–166. 10.1113/jphysiol.1986.sp015965
33
HasselmoM. E.ShayC. (2014). Grid cell firing patterns may arise from feedback interaction between intrinsic rebound spiking and transverse traveling waves with multiple heading angles. Front. Syst. Neurosci.8:201. 10.3389/fnsys.2014.00201
34
HeitlerW. J. (1980). “Neural mechanisms of central pattern generation in the crayfish swimmeret system,” in Neurobiology of Invertebrates, (Vol. 23) ed. SalánkiJ. (Budapest: Pergamon Press-Akadémiai Kiadó), 369–383.
35
HinmanJ. R.BrandonM. P.ClimerJ. R.ChapmanG. W.HasselmoM. E. (2016). Multiple running speed signals in medial entorhinal cortex. Neuron91, 666–679. 10.1016/j.neuron.2016.06.027
36
IvanenkoY. P.PoppeleR. E.LacquanitiF. (2006). Spinal cord maps of spatiotemporal α-motoneuron activation in humans walking at different speeds. J. Neurophysiol.95, 602–618. 10.1152/jn.00767.2005
37
JordanL. M. (1991). “Brainstem and spinal cord mechanisms for the initiation of locomotion,” in Neurobiological Basis of Human Locomotion, eds ShimamuraM.GrillnerS.EdgertonV. R. (Tokyo: Japan Scientific Societies Press), 3–20.
38
KaskeA.WinbergG.CösterJ. (2003). Emergence of coherent traveling waves controlling quadruped gaits in a two-dimensional spinal cord model. Biol. Cybern.88, 20–32. 10.1007/s00422-002-0335-0
39
KiehnO. (2016). Decoding the organization of spinal circuits that control locomotion. Nat. Rev. Neurosci.17, 224–238. 10.1038/nrn.2016.9
40
KiehnO.KjaerulffO. (1996). Spatiotemporal characteristics of 5-HT and dopamine-induced rhythmic hindlimb activity in the in vitro neonatal rat. J. Neurophysiol.75, 1472–1482.
41
KimuraY.OkamuraY.HigashijimaS. (2006). alx, a zebrafish homolog of Chx10, marks ipsilateral descending excitatory interneurons that participate in the regulation of spinal locomotor circuits. J. Neurosci.26, 5684–5697. 10.1523/JNEUROSCI.4993-05.2006
42
KoenigJ.LinderA. N.LeutgebJ. K.LeutgebS. (2011). The spatial periodicity of grid cells is not sustained during reduced theta oscillations. Science332, 592–595. 10.1126/science.1201685
43
Lafreniere-RoulaM.McCreaD. A. (2005). Deletions of rhythmic motoneuron activity during fictive locomotion and scratch provide clues to the organization of the mammalian central pattern generator. J. Neurophysiol.94, 1120–1132. 10.1152/jn.00216.2005
44
LennardP. R. (1985). Afferent perturbations during “monopodal” swimming movements in the turtle: phase-dependent cutaneous modulation and proprioceptive resetting of the locomotor rhythm. J. Neurosci.5, 1434–1445.
45
LevineA. J.HinckleyC. A.HildeK. L.DriscollS. P.PoonT. H.MontgomeryJ. M.et al. (2014). Identification of a cellular node for motor control pathways. Nat. Neurosci.17, 586–593. 10.1038/nn.3675
46
LevyS. L.WhiteJ. J.SillitoeR. V. (2015). “Wheat germ agglutinin (WGA) tracing: a classic approach for unraveling neural circuitry,” in Neural Tracing Methods: Tracing Neurons and their Connections, Neuromethods, (Vol. 92) ed. ArenkielB. R. (New York, NY: Springer), 51–66.
47
LiW. C. (2011). Generation of locomotion rhythms without inhibition in vertebrates: the search for pacemaker neurons. Integr. Comp. Biol.51, 879–889. 10.1093/icb/icr021
48
LiW. C.HigashijimaS.ParryD. M.RobertsA.SoffeS. R. (2004). Primitive roles for inhibitory interneurons in developing frog spinal cord. J. Neurosci.24, 5840–5848. 10.1523/JNEUROSCI.1633-04.2004
49
LismanJ. (2014). Two-phase model of the basal ganglia: implications for discontinuous control of the motor system. Philos. Trans. R. Soc. Lond. B Biol. Sci.369:20130489. 10.1098/rstb.2013.0489
50
LjunggrenE. E.HauptS.AusbornJ.AmpatzisK.El ManiraA. (2014). Optogenetic activation of excitatory premotor interneurons is sufficient to generate coordinated locomotor activity in larval zebrafish. J. Neurosci.34, 134–139. 10.1523/JNEUROSCI.4087-13.2014
51
ManterJ. T. (1938). The dynamics of quadrupedal walking. J. Exp. Biol.15, 522–540.
52
MaskosU.KissaK.St ClomentC.BrûletP. (2002). Retrograde trans-synaptic transfer of green fluorescent protein allows the genetic mapping of neuronal circuits in transgenic mice. Proc. Natl. Acad. Sci. U S A99, 10120–10125. 10.1073/pnas.152266799
53
McCreaD. A.RybakI. A. (2008). Organization of mammalian locomotor rhythm and pattern generation. Brain Res. Rev.57, 134–146. 10.1016/j.brainresrev.2007.08.006
54
MinicZ.ZhangY.MaoG.GoshgarianH. G. (2016). Transporter protein-coupled DCPX nanoconjugates induce diaphragmatic recovery after SCI by blocking Adenosine A1 receptors. J. Neurosci.38, 3441–3452. 10.1523/JNEUROSCI.2577-15.2016
55
MundellN. A.BeierK. T.PanY. A.LapanS. W.Göz AytürkD.BerezovskiiV. K.et al. (2015). Vesicular stomatitis virus enables gene transfer and transsynaptic tracing in a wide range of organisms. J. Comp. Neurol.523, 1639–1663. 10.1002/cne.23761
56
PastollH.SolankaL.van RossumM. C. W.NolanM. F. (2013). Feedback inhibition enables theta-nested γ oscillations and grid firing fields. Neuron77, 141–154. 10.1016/j.neuron.2012.11.032
57
PatelJ.FujisawaS.BerényiA.RoyerS.BuszákiG. (2012). Traveling theta waves along the entire septotemporal axis of the hippocampus. Neuron75, 410–417. 10.1016/j.neuron.2012.07.015
58
PérezT.TapiaJ. A.MirassoC. R.García-OjalvoJ.QuevedoJ.CuellarC. A.et al. (2009). An intersegmental neuronal architecture for spinal wave propagation under deletions. J. Neurosci.29, 10254–10263. 10.1523/JNEUROSCI.1737-09.2009
59
PerretC.CabelguenJ. M. (1976). Central and reflex participation in the timing of a locomotor bifunctional muscle, the semi-tendinosus, in the cat. Brain Res.106, 390–395. 10.1016/0006-8993(76)91035-0
60
PfeifferB. E.FosterD. J. (2015). Autoassociative dynamics in the generation of sequences of hippocampal place cells. Science349, 180–183. 10.1126/science.aaa9633
61
RohJ.CheungV. C.BizziE. (2011). Modules in the brainstem and spinal cord underlying motor behaviors. J. Neurophysiol.106, 1363–1378. 10.1152/jn.00842.2010
62
RybakI. A.DoughertyK. J.ShevtsovaN. A. (2015). Organization of the mammalian locomotor CPG: review of computational model and circuit architectures based on genetically identified spinal interneurons. eNeuro2:ENEURO.0069-15.2015. 10.1523/ENEURO.0069-15.2015
63
RybakI. A.ShevtsovaN. A.Lafreniere-RoulaM.McCreaD. A. (2006). Modelling spinal circuitry involved in locomotor pattern generation: insights from deletions during fictive locomotion. J. Physiol.577, 617–639. 10.1113/jphysiol.2006.118703
64
SaltielP.d’AvellaA.Wyler-DudaK.BizziE. (2016). Synergy temporal sequences and topography in the spinal cord: evidence for a traveling wave in frog locomotion. Brain Struct. Funct.221, 3869–3890. 10.1007/s00429-015-1133-5
65
SaltielP.RossignolS. (2004a). Critical points in the forelimb fictive locomotor cycle and motor coordination: effects of phasic retractions and protractions of the shoulder in the cat. J. Neurophysiol.92, 1342–1356. 10.1152/jn.00564.2003
66
SaltielP.RossignolS. (2004b). Critical points in the forelimb fictive locomotor cycle and motor coordination: evidence from the effects of tonic proprioceptive perturbations in the cat. J. Neurophysiol.92, 1329–1341. 10.1152/jn.00563.2003
67
SaltielP.TreschM. C.BizziE. (1998). Spinal cord modular organization and rhythm generation: an NMDA iontophoretic study in the frog. J. Neurophysiol.80, 2323–2339.
68
SaltielP.Wyler-DudaK.d’AvellaA.AjemianR.BizziE. (2005). Localization and connectivity in spinal interneuronal networks: the adduction-caudal extension-flexion rhythm in the frog. J. Neurophysiol.94, 2120–2138. 10.1152/jn.00117.2005
69
SaltielP.Wyler-DudaK.d’AvellaA.TreschM. C.BizziE. (2001). Muscle synergies encoded within the spinal cord: evidence from focal intraspinal NMDA iontophoresis in the frog. J. Neurophysiol.85, 605–619.
70
ShayC. F.FerranteM.ChapmanW. G.IV.HasselmoM. E. (2016). Rebound spiking in layer II medial entorhinal cortex stellate cells: possible mechanism of grid cell function. Neurobiol. Learn. Mem.129, 83–98. 10.1016/j.nlm.2015.09.004
71
ShibataN.ChenP.-S.DixonE. G.WolfP. D.DanieleyN. D.SmithW. M.et al. (1988). Influence of shock strength and timing on induction of ventricular arrhythmias in dogs. Am. J. Physiol.255, H891–H901.
72
SteinP. S.Daniels-McQueenS. (2004). Variations in motor patterns dring fictive rostral scratching in the turtle: knee-related deletions. J. Neurophysiol.91, 2380–2384. 10.1152/jn.01184.2003
73
SugitaM.ShibaY. (2005). Genetic tracing shows segregation of taste neuronal circuitries for bitter and sweet. Science309, 781–785. 10.1126/science.1110787
74
TingL. H.MacphersonJ. M. (2004). A limited set of muscle synergies for force control during a postural task. J. Neurophysiol.93, 609–613. 10.1152/jn.00681.2004
75
TreschM. C.KiehnO. (1999). Coding of locomotor phase in populations of neurons in rostral and caudal segments of the neonatal rat lumbar spinal cord. J. Neurophysiol.82, 3563–3574.
76
TreschM. C.SaltielP.BizziE. (1999). The construction of movement by the spinal cord. Nat. Neurosci.2, 162–167. 10.1038/5721
77
WallénP.GrillnerS. (1987). N-methyl-D-aspartate receptor-induced, inherent oscillatory activity in neurons active during fictive locomotion in the lamprey. J. Neurosci.7, 2745–2755.
78
WickershamI. R.LyonD. C.BarnardR. J.MoriT.FinkeS.ConzelmannK.et al. (2007). Monosynaptic restriction of transsynaptic tracing from single, genetically targeted neurons. Neuron53, 639–647. 10.1016/j.neuron.2007.01.033
79
WigginT. D.AndersonT. M.EianJ.PeckJ. H.MasinoM. A. (2012). Episodic swimming in the larval zebrafish is generated by a spatially distributed spinal network with modular functional organization. J. Neurophysiol.108, 925–934. 10.1152/jn.00233.2012
80
WinfreeA. T. (1972). Spiral waves of chemical activity. Science175, 634–636. 10.1126/science.175.4022.634
81
WinfreeA. T. (1977). Phase control of neural pacemakers. Science197, 761–763. 10.1126/science.887919
82
WinfreeA. T. (1983). Sudden cardiac death: a problem in topology. Sci. Am.248, 144–149, 152–157, 160–161. 10.1038/scientificamerican0583-144
83
WinfreeA. T. (1987). When Time Breaks Down.Princeton, NJ: Princeton University Press.
84
WinfreeA. T. (1989). Electrical instability in cardiac muscle: phase singularities and rotors. J. Theor. Biol.138, 353–405. 10.1016/s0022-5193(89)80200-0
85
YakovenkoS.MushahwarV.VanderHorstV.HolstegeG.ProchazkaA. (2002). Spatiotemporal activation of lumbosacral motoneurons in the locomotor step cycle. J. Neurophysiol.87, 1542–1553. 10.1152/jn.00479.2001
86
YartsevM. M.WitterM. P.UlanovskyN. (2011). Grid cells without theta oscillations in the entorhinal cortex of bats. Nature488, 103–107. 10.1038/nature10583
87
ZaikinA. N.ZhabotinskyA. M. (1970). Concentration wave propagation in two-dimensional liquid-phase self-oscillating system. Nature225, 535–537. 10.1038/225535b0
88
ZhangH.JacobsJ. (2015). Traveling theta waves in the human hippocampus. J. Neurosci.35, 12477–12487. 10.1523/JNEUROSCI.5102-14.2015
89
ZhangJ.LanuzaG. M.BritzO.WangZ.SiembabV. C.ZhangY.et al. (2014). V1 and v2b interneurons secure the alternating flexor-extensor motor activity mice require for limbed locomotion. Neuron82, 138–150. 10.1016/j.neuron.2014.02.013
90
ZhengC.BieriK. W.HsiaoY. T.ColginL. L. (2016). Spatial sequence coding differs during slow and fast γ rhythms in the hippocampus. Neuron89, 398–408. 10.1016/j.neuron.2015.12.005
91
ZhongG.ShevtsovaN. A.RybalI. A.Harris-WarrickR. M. (2012). Neuronal activity in the isolated mouse spinal cord during spontaneous deletions in fictive locomotion: insights in locomotor central pattern generator organization. J. Physiol.590, 4735–4759. 10.1113/jphysiol.2012.240895
Summary
Keywords
central pattern generator, locomotion, critical point shifts, traveling wave, temporal grid, hippocampus, Winfree’s phase singularities, spinal cord
Citation
Saltiel P, d’Avella A, Tresch MC, Wyler K and Bizzi E (2017) Critical Points and Traveling Wave in Locomotion: Experimental Evidence and Some Theoretical Considerations. Front. Neural Circuits 11:98. doi: 10.3389/fncir.2017.00098
Received
18 July 2017
Accepted
20 November 2017
Published
08 December 2017
Volume
11 - 2017
Edited by
Jean-Marie Cabelguen, Institut National de la Santé et de la Recherche Médicale, France
Reviewed by
Elias Manjarrez, Benemérita Universidad Autónoma de Puebla, Mexico; Joseph Ayers, Northeastern University, United States
Updates
Copyright
© 2017 Saltiel, d’Avella, Tresch, Wyler and Bizzi.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Philippe Saltiel saltiel@mit.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.