The significance of dynamical architecture for adaptive responses to mechanical loads during rhythmic behavior

Tradeoffs in the energetic value of neuromodulation in a closed-loop neuromechanical system

Rhythmic motor behaviors controlled by neuromechanical systems, consisting of central neural circuitry, biomechanics, and sensory feedback, show efficiency in energy expenditure. The biomechanical elements (e.g., muscles) are modulated by peripheral neuromodulation which may improve their strength and speed properties. However, there are relatively few studies on neuromodulatory control of muscle function and metabolic mechanical efficiency in neuromechanical systems. To investigate the role of neuromodulation on the system's mechanical efficiency, we consider a neuromuscular model of motor patterns for feeding in the marine mollusk Aplysia californica. By incorporating muscle energetics and neuromodulatory effects into the model, we demonstrate tradeoffs in the energy efficiency of Aplysia's rhythmic swallowing behavior as a function of the level of neuromodulation. A robust efficiency optimum arises from an intermediate level of neuromodulation, and excessive neuromodulation may be inefficient and disadvantageous to an animal's metabolism. This optimum emerges from physiological constraints imposed upon serotonergic modulation trajectories on the energy efficiency landscape. Our results may lead to experimentally testable hypotheses of the role of neuromodulation in rhythmic motor control.

Stable heteroclinic channels as a decision-making model: overcoming low signal-to-noise ratio with mutual inhibition

Bio-inspired robot controllers are becoming more complex as we strive to make them more robust to, and flexible in, noisy, real-world environments. A stable heteroclinic network (SHN) is a dynamical system that produces cyclical state transitions using noisy input. SHN-based robot controllers enable sensory input to be integrated at the phase-space level of the controller, thus simplifying sensor-integrated, robot control methods. In this work, we investigate the mechanism that drives branching state trajectories in SHNs. We liken the branching state trajectories to decision-splits imposed into the system, which opens the door for more sophisticated controls-all driven by sensory input. This work provides guidelines to systematically define an SHN topology, and increase the rate at which desired decision states in the topology are chosen. Ultimately, we are able to control the rate at which desired decision states activate for input signal-to-noise ratios across six orders of magnitude.

Variational analysis of sensory feedback mechanisms in powerstroke-recovery systems

Although the raison d'etre of the brain is the survival of the body, there are relatively few theoretical studies of closed-loop rhythmic motor control systems. In this paper we provide a unified framework, based on variational analysis, for investigating the dual goals of performance and robustness in powerstroke-recovery systems. To demonstrate our variational method, we augment two previously published closed-loop motor control models by equipping each model with a performance measure based on the rate of progress of the system relative to a spatially extended external substrate-such as a long strip of seaweed for a feeding task, or progress relative to the ground for a locomotor task. The sensitivity measure quantifies the ability of the system to maintain performance in response to external perturbations, such as an applied load. Motivated by a search for optimal design principles for feedback control achieving the complementary requirements of efficiency and robustness, we discuss the performance-sensitivity patterns of the systems featuring different sensory feedback architectures. In a paradigmatic half-center oscillator-motor system, we observe that the excitation-inhibition property of feedback mechanisms determines the sensitivity pattern while the activation-inactivation property determines the performance pattern. Moreover, we show that the nonlinearity of the sigmoid activation of feedback signals allows the existence of optimal combinations of performance and sensitivity. In a detailed hindlimb locomotor system, we find that a force-dependent feedback can simultaneously optimize both performance and robustness, while length-dependent feedback variations result in significant performance-versus-sensitivity tradeoffs. Thus, this work provides an analytical framework for studying feedback control of oscillations in nonlinear dynamical systems, leading to several insights that have the potential to inform the design of control or rehabilitation systems.

Full Hill-type muscle model of the I1/I3 retractor muscle complex in Aplysia californica

The coordination of complex behavior requires knowledge of both neural dynamics and the mechanics of the periphery. The feeding system of Aplysia californica is an excellent model for investigating questions in soft body systems' neuromechanics because of its experimental tractability. Prior work has attempted to elucidate the mechanical properties of the periphery by using a Hill-type muscle model to characterize the force generation capabilities of the key protractor muscle responsible for moving Aplysia's grasper anteriorly, the I2 muscle. However, the I1/I3 muscle, which is the main driver of retractions of Aplysia's grasper, has not been characterized. Because of the importance of the musculature's properties in generating functional behavior, understanding the properties of muscles like the I1/I3 complex may help to create more realistic simulations of the feeding behavior of Aplysia, which can aid in greater understanding of the neuromechanics of soft-bodied systems. To bridge this gap, in this work, the I1/I3 muscle complex was characterized using force-frequency, length-tension, and force-velocity experiments and showed that a Hill-type model can accurately predict its force-generation properties. Furthermore, the muscle's peak isometric force and stiffness were found to exceed those of the I2 muscle, and these results were analyzed in the context of prior studies on the I1/I3 complex's kinematics in vivo.

On relaxation times of heteroclinic dynamics

Heteroclinic dynamics provide a suitable framework for describing transient dynamics such as cognitive processes in the brain. It is appreciated for being well reproducible and at the same time highly sensitive to external input. It is supposed to capture features of switching statistics between metastable states in the brain. Beyond the high sensitivity, a further desirable feature of these dynamics is to enable a fast adaptation to new external input. In view of this, we analyze relaxation times of heteroclinic motion toward a new resting state, when oscillations in heteroclinic networks are arrested by a quench of a bifurcation parameter from a parameter regime of oscillations to a regime of equilibrium states. As it turns out, the relaxation is underdamped and depends on the nesting of the attractor space, the size of the attractor's basin of attraction, the depth of the quench, and the level of noise. In the case of coupled heteroclinic units, it depends on the coupling strength, the coupling type, and synchronization between different units. Depending on how these factors are combined, finite relaxation times may support or impede a fast switching to new external input. Our results also shed some light on the discussion of how the stability of a system changes with its complexity.

Tutorial: using NEURON for neuromechanical simulations

The dynamical properties of the brain and the dynamics of the body strongly influence one another. Their interaction generates complex adaptive behavior. While a wide variety of simulation tools exist for neural dynamics or biomechanics separately, there are few options for integrated brain-body modeling. Here, we provide a tutorial to demonstrate how the widely-used NEURON simulation platform can support integrated neuromechanical modeling. As a first step toward incorporating biomechanics into a NEURON simulation, we provide a framework for integrating inputs from a "periphery" and outputs to that periphery. In other words, "body" dynamics are driven in part by "brain" variables, such as voltages or firing rates, and "brain" dynamics are influenced by "body" variables through sensory feedback. To couple the "brain" and "body" components, we use NEURON's pointer construct to share information between "brain" and "body" modules. This approach allows separate specification of brain and body dynamics and code reuse. Though simple in concept, the use of pointers can be challenging due to a complicated syntax and several different programming options. In this paper, we present five different computational models, with increasing levels of complexity, to demonstrate the concepts of code modularity using pointers and the integration of neural and biomechanical modeling within NEURON. The models include: (1) a neuromuscular model of calcium dynamics and muscle force, (2) a neuromechanical, closed-loop model of a half-center oscillator coupled to a rudimentary motor system, (3) a closed-loop model of neural control for respiration, (4) a pedagogical model of a non-smooth "brain/body" system, and (5) a closed-loop model of feeding behavior in the sea hare Aplysia californica that incorporates biologically-motivated non-smooth dynamics. This tutorial illustrates how NEURON can be integrated with a broad range of neuromechanical models.

Code available at

https://github.com/fietkiewicz/PointerBuilder.

Heteroclinic cycling and extinction in May-Leonard models with demographic stochasticity

May and Leonard (SIAM J Appl Math 29:243-253, 1975) introduced a three-species Lotka-Volterra type population model that exhibits heteroclinic cycling. Rather than producing a periodic limit cycle, the trajectory takes longer and longer to complete each "cycle", passing closer and closer to unstable fixed points in which one population dominates and the others approach zero. Aperiodic heteroclinic dynamics have subsequently been studied in ecological systems (side-blotched lizards; colicinogenic Escherichia coli), in the immune system, in neural information processing models ("winnerless competition"), and in models of neural central pattern generators. Yet as May and Leonard observed "Biologically, the behavior (produced by the model) is nonsense. Once it is conceded that the variables represent animals, and therefore cannot fall below unity, it is clear that the system will, after a few cycles, converge on some single population, extinguishing the other two." Here, we explore different ways of introducing discrete stochastic dynamics based on May and Leonard's ODE model, with application to ecological population dynamics, and to a neuromotor central pattern generator system. We study examples of several quantitatively distinct asymptotic behaviors, including total extinction of all species, extinction to a single species, and persistent cyclic dominance with finite mean cycle length.

Variational and phase response analysis for limit cycles with hard boundaries, with applications to neuromechanical control problems

Motor systems show an overall robustness, but because they are highly nonlinear, understanding how they achieve robustness is difficult. In many rhythmic systems, robustness against perturbations involves response of both the shape and the timing of the trajectory. This makes the study of robustness even more challenging. To understand how a motor system produces robust behaviors in a variable environment, we consider a neuromechanical model of motor patterns in the feeding apparatus of the marine mollusk Aplysia californica (Shaw et al. in J Comput Neurosci 38(1):25-51, 2015; Lyttle et al. in Biol Cybern 111(1):25-47, 2017). We established in (Wang et al. in SIAM J Appl Dyn Syst 20(2):701-744, 2021. https://doi.org/10.1137/20M1344974 ) the tools for studying combined shape and timing responses of limit cycle systems under sustained perturbations and here apply them to study robustness of the neuromechanical model against increased mechanical load during swallowing. Interestingly, we discover that nonlinear biomechanical properties confer resilience by immediately increasing resistance to applied loads. In contrast, the effect of changed sensory feedback signal is significantly delayed by the firing rates' hard boundary properties. Our analysis suggests that sensory feedback contributes to robustness in swallowing primarily by shifting the timing of neural activation involved in the power stroke of the motor cycle (retraction). This effect enables the system to generate stronger retractor muscle forces to compensate for the increased load, and hence achieve strong robustness. The approaches that we are applying to understanding a neuromechanical model in Aplysia, and the results that we have obtained, are likely to provide insights into the function of other motor systems that encounter changing mechanical loads and hard boundaries, both due to mechanical and neuronal firing properties.

Sensory feedback expands dynamic complexity and aids in robustness against noise

It has been hypothesized that sensory feedback is a critical component in determining the functionality of a central pattern generator. To test this, Yu and Thomas's recent work Yu and Thomas (Biol Cybern 115(2):135-160, 2021) built a model of a half-center oscillator coupled to a simple muscular model with sensory feedback. They showed that sensory feedback increases robustness against external noise, while simultaneously expanding the potential repertoire of functions the half-center oscillator can perform. However, they show that this comes at the cost of robustness against internal noise.

Quantitative comparison of the mean-return-time phase and the stochastic asymptotic phase for noisy oscillators

Seminal work by A. Winfree and J. Guckenheimer showed that a deterministic phase variable can be defined either in terms of Poincaré sections or in terms of the asymptotic (long-time) behaviour of trajectories approaching a stable limit cycle. However, this equivalence between the deterministic notions of phase is broken in the presence of noise. Different notions of phase reduction for a stochastic oscillator can be defined either in terms of mean-return-time sections or as the argument of the slowest decaying complex eigenfunction of the Kolmogorov backwards operator. Although both notions of phase enjoy a solid theoretical foundation, their relationship remains unexplored. Here, we quantitatively compare both notions of stochastic phase. We derive an expression relating both notions of phase and use it to discuss differences (and similarities) between both definitions of stochastic phase for (i) a spiral sink motivated by stochastic models for electroencephalograms, (ii) noisy limit-cycle systems-neuroscience models, and (iii) a stochastic heteroclinic oscillator inspired by a simple motor-control system.

A homeostasis criterion for limit cycle systems based on infinitesimal shape response curves

Homeostasis occurs in a control system when a quantity remains approximately constant as a parameter, representing an external perturbation, varies over some range. Golubitsky and Stewart (J Math Biol 74(1-2):387-407, 2017) developed a notion of infinitesimal homeostasis for equilibrium systems using singularity theory. Rhythmic physiological systems (breathing, locomotion, feeding) maintain homeostasis through control of large-amplitude limit cycles rather than equilibrium points. Here we take an initial step to study (infinitesimal) homeostasis for limit-cycle systems in terms of the average of a quantity taken around the limit cycle. We apply the "infinitesimal shape response curve" (iSRC) introduced by Wang et al. (SIAM J Appl Dyn Syst 82(7):1-43, 2021) to study infinitesimal homeostasis for limit-cycle systems in terms of the mean value of a quantity of interest, averaged around the limit cycle. Using the iSRC, which captures the linearized shape displacement of an oscillator upon a static perturbation, we provide a formula for the derivative of the averaged quantity with respect to the control parameter. Our expression allows one to identify homeostasis points for limit cycle systems in the averaging sense. We demonstrate in the Hodgkin-Huxley model and in a metabolic regulatory network model that the iSRC-based method provides an accurate representation of the sensitivity of averaged quantities.

Generation of Direct-, Retrograde-, and Source-Wave Gaits in Multi-Legged Locomotion in a Decentralized Manner via Embodied Sensorimotor Interaction

Multi-legged animals show several types of ipsilateral interlimb coordination. Millipedes use a direct-wave gait, in which the swing leg movements propagate from posterior to anterior. In contrast, centipedes use a retrograde-wave gait, in which the swing leg movements propagate from anterior to posterior. Interestingly, when millipedes walk in a specific way, both direct and retrograde waves of the swing leg movements appear with the waves' source, which we call the source-wave gait. However, the gait generation mechanism is still unclear because of the complex nature of the interaction between neural control and dynamic body systems. The present study used a simple model to understand the mechanism better, primarily how local sensory feedback affects multi-legged locomotion. The model comprises a multi-legged body and its locomotion control system using biologically inspired oscillators with local sensory feedback, phase resetting. Each oscillator controls each leg independently. Our simulation produced the above three types of animal gaits. These gaits are not predesigned but emerge through the interaction between the neural control and dynamic body systems through sensory feedback (embodied sensorimotor interaction) in a decentralized manner. The analytical description of these gaits' solution and stability clarifies the embodied sensorimotor interaction's functional roles in the interlimb coordination.

Dynamical consequences of sensory feedback in a half-center oscillator coupled to a simple motor system

We investigate a simple model for motor pattern generation that combines central pattern generator (CPG) dynamics with a sensory feedback (FB) mechanism. Our CPG comprises a half-center oscillator with conductance-based Morris-Lecar model neurons. Output from the CPG drives a push-pull motor system with biomechanics based on experimental data. A sensory feedback conductance from the muscles allows modulation of the CPG activity. We consider parameters under which the isolated CPG system has either "escape" or "release" dynamics, and we study both inhibitory and excitatory feedback conductances. We find that increasing the FB conductance relative to the CPG conductance makes the system more robust against external perturbations, but more susceptible to internal noise. Conversely, increasing the CPG conductance relative to the FB conductance has the opposite effects. We find that the "closed-loop" system, with sensory feedback in place, exhibits a richer repertoire of behaviors than the "open-loop" system, with motion determined entirely by the CPG dynamics. Moreover, we find that purely feedback-driven motor patterns, analogous to a chain reflex, occur only in the inhibition-mediated system. Finally, for pattern generation systems with inhibition-mediated sensory feedback, we find that the distinction between escape- and release-mediated CPG mechanisms is diminished in the presence of internal noise. Our observations support an anti-reductionist view of neuromotor physiology: Understanding mechanisms of robust motor control requires studying not only the central pattern generator circuit in isolation, but the intact closed-loop system as a whole.

Shape versus timing: linear responses of a limit cycle with hard boundaries under instantaneous and static perturbation

When dynamical systems that produce rhythmic behaviors operate within hard limits, they may exhibit limit cycles with sliding components, that is, closed isolated periodic orbits that make and break contact with a constraint surface. Examples include heel-ground interaction in locomotion, firing rate rectification in neural networks, and stick-slip oscillators. In many rhythmic systems, robustness against external perturbations involves response of both the shape and the timing of the limit cycle trajectory. The existing methods of infinitesimal phase response curve (iPRC) and variational analysis are well established for quantifying changes in timing and shape, respectively, for smooth systems. These tools have recently been extended to nonsmooth dynamics with transversal crossing boundaries. In this work, we further extend the iPRC method to nonsmooth systems with sliding components, which enables us to make predictions about the synchronization properties of weakly coupled stick-slip oscillators. We observe a new feature of the isochrons in a planar limit cycle with hard sliding boundaries: a nonsmooth kink in the asymptotic phase function, originating from the point at which the limit cycle smoothly departs the constraint surface, and propagating away from the hard boundary into the interior of the domain. Moreover, the classical variational analysis neglects timing information and is restricted to instantaneous perturbations. By defining the "infinitesimal shape response curve" (iSRC), we incorporate timing sensitivity of an oscillator to describe the shape response of this oscillator to parametric perturbations. In order to extract timing information, we also develop a "local timing response curve" (lTRC) that measures the timing sensitivity of a limit cycle within any given region. We demonstrate in a specific example that taking into account local timing sensitivity in a nonsmooth system greatly improves the accuracy of the iSRC over global timing analysis given by the iPRC.

Control for multifunctionality: bioinspired control based on feeding in Aplysia californica

Animals exhibit remarkable feats of behavioral flexibility and multifunctional control that remain challenging for robotic systems. The neural and morphological basis of multifunctionality in animals can provide a source of bioinspiration for robotic controllers. However, many existing approaches to modeling biological neural networks rely on computationally expensive models and tend to focus solely on the nervous system, often neglecting the biomechanics of the periphery. As a consequence, while these models are excellent tools for neuroscience, they fail to predict functional behavior in real time, which is a critical capability for robotic control. To meet the need for real-time multifunctional control, we have developed a hybrid Boolean model framework capable of modeling neural bursting activity and simple biomechanics at speeds faster than real time. Using this approach, we present a multifunctional model of Aplysia californica feeding that qualitatively reproduces three key feeding behaviors (biting, swallowing, and rejection), demonstrates behavioral switching in response to external sensory cues, and incorporates both known neural connectivity and a simple bioinspired mechanical model of the feeding apparatus. We demonstrate that the model can be used for formulating testable hypotheses and discuss the implications of this approach for robotic control and neuroscience.

Rapid Adaptation to Changing Mechanical Load by Ordered Recruitment of Identified Motor Neurons

As they interact with their environment and encounter challenges, animals adjust their behavior on a moment-to-moment basis to maintain task fitness. This dynamic process of adaptive motor control occurs in the nervous system, but an understanding of the biomechanics of the body is essential to properly interpret the behavioral outcomes. To study how animals respond to changing task conditions, we used a model system in which the functional roles of identified neurons and the relevant biomechanics are well understood and can be studied in intact behaving animals: feeding in the marine mollusc Aplysia We monitored the motor neuronal output of the feeding circuitry as intact animals fed on uniform food stimuli under unloaded and loaded conditions, and we measured the force of retraction during loaded swallows. We observed a previously undescribed pattern of force generation, which can be explained within the appropriate biomechanical context by the activity of just a few key, identified motor neurons. We show that, when encountering load, animals recruit identified retractor muscle motor neurons for longer and at higher frequency to increase retraction force duration. Our results identify a mode by which animals robustly adjust behavior to their environment, which is experimentally tractable to further mechanistic investigation.

A neuromechanical model exploring the role of the common inhibitor motor neuron in insect locomotion

In this work, we analyze a simplified, dynamical, closed-loop, neuromechanical simulation of insect joint control. We are specifically interested in two elements: (1) how slow muscle fibers may serve as temporal integrators of sensory feedback and (2) the role of common inhibitory (CI) motor neurons in resetting this integration when the commanded position changes, particularly during steady-state walking. Despite the simplicity of the model, we show that slow muscle fibers increase the accuracy of limb positioning, even for motions much shorter than the relaxation time of the fiber; this increase in accuracy is due to the slow dynamics of the fibers; the CI motor neuron plays a critical role in accelerating muscle relaxation when the limb moves to a new position; as in the animal, this architecture enables the control of the stance phase speed, independent of swing phase amplitude or duration, by changing the gain of sensory feedback to the stance phase muscles. We discuss how this relates to other models, and how it could be applied to robotic control.

Sensitive Finite-State Computations Using a Distributed Network With a Noisy Network Attractor

We exhibit a class of smooth continuous-state neural-inspired networks composed of simple nonlinear elements that can be made to function as a finite-state computational machine. We give an explicit construction of arbitrary finite-state virtual machines in the spatiotemporal dynamics of the network. The dynamics of the functional network can be completely characterized as a "noisy network attractor" in phase space operating in either an "excitable" or a "free-running" regime, respectively, corresponding to excitable or heteroclinic connections between states. The regime depends on the sign of an "excitability parameter." Viewing the network as a nonlinear stochastic differential equation where a deterministic (signal) and/or a stochastic (noise) input is applied to any element, we explore the influence of the signal-to-noise ratio on the error rate of the computations. The free-running regime is extremely sensitive to inputs: arbitrarily small amplitude perturbations can be used to perform computations with the system as long as the input dominates the noise. We find a counter-intuitive regime where increasing noise amplitude can lead to more, rather than less, accurate computation. We suggest that noisy network attractors will be useful for understanding neural networks that reliably and sensitively perform finite-state computations in a noisy environment.

Design process and tools for dynamic neuromechanical models and robot controllers

We present a serial design process with associated tools to select parameter values for a posture and locomotion controller for simulation of a robot. The controller is constructed from dynamic neuron and synapse models and simulated with the open-source neuromechanical simulator AnimatLab 2. Each joint has a central pattern generator (CPG), whose neurons possess persistent sodium channels. The CPG rhythmically inhibits motor neurons that control the servomotor's velocity. Sensory information coordinates the joints in the leg into a cohesive stepping motion. The parameter value design process is intended to run on a desktop computer, and has three steps. First, our tool FEEDBACKDESIGN uses classical control methods to find neural and synaptic parameter values that stably and robustly control servomotor output. This method is fast, testing over 100 parameter value variations per minute. Next, our tool CPGDESIGN generates bifurcation diagrams and phase response curves for the CPG model. This reveals neural and synaptic parameter values that produce robust oscillation cycles, whose phase can be rapidly entrained to sensory feedback. It also designs the synaptic conductance of inter-joint pathways. Finally, to understand sensitivity to parameters and how descending commands affect a leg's stepping motion, our tool SIMSCAN runs batches of neuromechanical simulations with specified parameter values, which is useful for searching the parameter space of a complicated simulation. These design tools are demonstrated on a simulation of a robot, but may be applied to neuromechanical animal models or physical robots as well.

Robustness, flexibility, and sensitivity in a multifunctional motor control model

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.

Motor neuronal activity varies least among individuals when it matters most for behavior

How does motor neuronal variability affect behavior? To explore this question, we quantified activity of multiple individual identified motor neurons mediating biting and swallowing in intact, behaving Aplysia californica by recording from the protractor muscle and the three nerves containing the majority of motor neurons controlling the feeding musculature. We measured multiple motor components: duration of the activity of identified motor neurons as well as their relative timing. At the same time, we measured behavioral efficacy: amplitude of grasping movement during biting and amplitude of net inward food movement during swallowing. We observed that the total duration of the behaviors varied: Within animals, biting duration shortened from the first to the second and third bites; between animals, biting and swallowing durations varied. To study other sources of variation, motor components were divided by behavior duration (i.e., normalized). Even after normalization, distributions of motor component durations could distinguish animals as unique individuals. However, the degree to which a motor component varied among individuals depended on the role of that motor component in a behavior. Motor neuronal activity that was essential for the expression of biting or swallowing was similar among animals, whereas motor neuronal activity that was not essential for that behavior varied more from individual to individual. These results suggest that motor neuronal activity that matters most for the expression of a particular behavior may vary least from individual to individual. Shaping individual variability to ensure behavioral efficacy may be a general principle for the operation of motor systems.