Motor systems must adapt to perturbations and changing conditions both within and outside the body. We refer to the ability of a system to maintain performance despite perturbations as “robustness,” and the ability of a system to deploy alternative strategies that improve fitness as “flexibility.” Different classes of pattern-generating circuits yield dynamics with differential sensitivities to perturbations and parameter variation. Depending on the task and the type of perturbation, high sensitivity can either facilitate or hinder robustness and flexibility. Here we explore the role of multiple coexisting oscillatory modes and sensory feedback in allowing multiphasic motor pattern generation to be both robust and flexible. As a concrete example, we focus on a nominal neuromechanical model of triphasic motor patterns in the feeding apparatus of the marine mollusk Aplysia californica. We find that the model can operate within two distinct oscillatory modes and that the system exhibits bistability between the two. In the “heteroclinic mode,” higher sensitivity makes the system more robust to changing mechanical loads, but less robust to internal parameter variations. In the “limit cycle mode,” lower sensitivity makes the system more robust to changes in internal parameter values, but less robust to changes in mechanical load. Finally, we show that overall performance on a variable feeding task is improved when the system can flexibly transition between oscillatory modes in response to the changing demands of the task. Thus, our results suggest that the interplay of sensory feedback and multiple oscillatory modes can allow motor systems to be both robust and flexible in a variable environment.

Many behaviors require reliably generating sequences of motor activity while adapting the activity to incoming sensory information. This process has often been conceptually explained as either fully dependent on sensory input (a chain reflex) or fully independent of sensory input (an idealized central pattern generator, or CPG), although the consensus of the field is that most neural pattern generators lie somewhere between these two extremes. Many mathematical models of neural pattern generators use limit cycles to generate the sequence of behaviors, but other models, such as a heteroclinic channel (an attracting chain of saddle points), have been suggested. To explore the range of intermediate behaviors between CPGs and chain reflexes, in this paper we describe a nominal model of swallowing in Aplysia californica. Depending upon the value of a single parameter, the model can transition from a generic limit cycle regime to a heteroclinic regime (where the trajectory slows as it passes near saddle points). We then study the behavior of the system in these two regimes and compare the behavior of the models with behavior recorded in the animal in vivo and in vitro. We show that while both pattern generators can generate similar behavior, the stable heteroclinic channel can better respond to changes in sensory input induced by load, and that the response matches the changes seen when a load is added in vivo. We then show that the underlying stable heteroclinic channel architecture exhibits dramatic slowing of activity when sensory and endogenous input is reduced, and show that similar slowing with removal of proprioception is seen in vitro. Finally, we show that the distributions of burst lengths seen in vivo are better matched by the distribution expected from a system operating in the heteroclinic regime than that expected from a generic limit cycle. These observations suggest that generic limit cycle models may fail to capture key aspects of Aplysia feeding behavior, and that alternative architectures such as heteroclinic channels may provide better descriptions.

The rodent hippocampus and entorhinal cortex contain spatially-modulated cells that serve as the basis for spatial coding. Both medial entorhinal grid cells and hippocampal place cells have been shown to encode spatial information across multiple spatial scales that increase along the dorsoventral axis of these structures. Place cells near the dorsal pole possess small, stable, and spatially selective firing fields, while ventral cells have larger, less stable and less spatially selective firing fields. One possible explanation for these dorsoventral changes in place field properties is that they arise as a result of similar dorsoventral differences in the properties of the grid cell inputs to place cells. Here we test the alternative hypothesis that dorsoventral place field differences are due to higher amounts of non-spatial inputs to ventral hippocampal cells. We use a computational model of the entorhinal-hippocampal network to assess the relative contributions of grid scale and non-spatial inputs in determining place field size and stability. In addition, we assess the consequences of grid node firing rate heterogeneity on place field stability. Our results suggest that dorsoventral differences in place cell properties can be better explained by changes in the amount of non-spatial inputs, rather than by changes in the scale of grid cell inputs, and that grid node heterogeneity may have important functional consequences. The observed gradient in field size may reflect a shift from processing primarily spatial information in the dorsal hippocampus to processing more non-spatial, contextual and emotional information near the ventral hippocampus.

An important problem in neuroscience is that of constructing quantitative measures of the similarity between neural spike trains. These measures can be used, for example, to assess the reliability of the response of a single neuron to repeated stimulus presentations, or to uncover relationships in the firing patterns of multiple neurons in a population. While several similarity measures have been proposed, the extent to which they take into account various biologically important spike train features such as bursts of spikes, or periods of inactivity remains poorly understood. Here we compare these measures using tests specifically designed to assess the sensitivity to bursts and silent periods. In addition, we propose two new measures. The first is designed to detect periods of shared silence between spike trains, while the second is designed to emphasize the presence of common bursts. To assist researchers in determining which measure is best suited to their particular data analysis needs, we also show how these measures can be combined and how their parameters can be determined on the basis of physiologically relevant quantities.