Nonlinear muscles, passive viscoelasticity and body taper conspire to create neuromechanical phase lags in anguilliform swimmers

Mechanosensory Feedback in Lamprey Swimming Models and Applications in the Field of Spinal Cord Regeneration

The central pattern generator (CPG) in anguilliform swimming has served as a model for examining the neural basis of locomotion. This system has been particularly valuable for the development of mathematical models. As our biological understanding of the neural basis of locomotion has expanded, so too have these models. Recently, there have been significant advancements in our understanding of the critical role that mechanosensory feedback plays in robust locomotion. This work has led to a push in the field of mathematical modeling to incorporate mechanosensory feedback into CPG models. In this perspective piece, we review advances in the development of these models and discuss how newer complex models can support biological investigation. We highlight lamprey spinal cord regeneration as an area that can both inform these models and benefit from them.

Computational Study of Stiffness-Tuning Strategies in Anguilliform Fish

Biological evidence demonstrates that fish can tune their body stiffness to improve thrust and efficiency during swimming locomotion. However, the stiffness-tuning strategies that maximize swimming speed or efficiency are still unclear. In the present study, a musculo-skeletal model of anguilliform fish is developed to study the properties of variable stiffness, in which the planar serial-parallel mechanism is used to model the body structure. The calcium ion model is adopted to simulate muscular activities and generate muscle force. Further, the relations among the forward speed, the swimming efficiency, and Young's modulus of the fish body are investigated. The results show that for certain body stiffness, the swimming speed and efficiency are increased with the tail-beat frequency until reaching the maximum value and then decreased. The peak speed and efficiency are also increased with the amplitude of muscle actuation. Anguilliform fish tend to vary their body stiffness to improve the swimming speed and efficiency at a high tail-beat frequency or small amplitude of muscle actuation. Furthermore, the midline motions of anguilliform fish are analyzed by the complex orthogonal decomposition (COD) method, and the discussions of fish motions associated with the variable body stiffness and the tail-beat frequency are also presented. Overall, the optimal swimming performance of anguilliform fish benefits from the matching relationships among the muscle actuation, the body stiffness, and the tail-beat frequency.

Effect of internal damping on locomotion in frictional environments

The gaits of undulating animals arise from a complex interaction of their central nervous system, muscle, connective tissue, bone, and environment. As a simplifying assumption, many previous studies have often assumed that sufficient internal force is available to produce observed kinematics, thus not focusing on quantifying the interconnection between muscle effort, body shape, and external reaction forces. This interplay, however, is critical to locomotion performance in crawling animals, especially when accompanied by body viscoelasticity. Moreover, in bioinspired robotic applications, the body's internal damping is indeed a parameter that the designer can tune. Still, the effect of internal damping is not well understood. This study explores how internal damping affects the locomotion performance of a crawler with a continuous, viscoelastic, nonlinear beam model. Crawler muscle actuation is modeled as a traveling wave of bending moment propagating posteriorly along the body. Consistent with the friction properties of the scales of snakes and limbless lizards, environmental forces are modeled using anisotropic Coulomb friction. It is found that by varying the crawler body's internal damping, the crawler's performance can be altered, and distinct gaits could be achieved, including changing the net locomotion direction from forward to back. We will discuss this forward and backward control and identify the optimal internal damping for peak crawling speed.

Flexible tensegrity wing design and insights in principles of swimming kinematics of batoid rays

A novel tensegrity wing design is first proposed which can emulate the kinematic waves of the pectoral fin of batoid rays and has a simple structure for manufacture. The attitude control and the regulation of wing natural frequency are realized by wing morphing. Then analytical insights in batoid ray swimming are gained by analyzing the analytical wing (cable)-fluid interaction model, whose parameters are determined based on the biological data. The stride length (traveled distance per cycle normalized by the body length (BL)) is shown to be almost invariant among different-sized rays if the phase and amplitude of wing flexion angles remain unchanged. This result is supported by biological data, 1.5 and 1.47 respectively for the manta ray and cownose ray, though their flapping frequencies (0.15-0.45 Hz and 0.64-1.25 Hz respectively) and body sizes (1.25 m and 0.15 m respectively) are very different, and similar to the expression for the carangiform fish swimming. In other words, the swimming kinematics of two different swimming forms are described by a similar analytical equation when the body resonance is exploited. The fluid force and cable tension are both found to be proportional to the fourth power of the body size and the square of the wing flapping frequency, which may tell that the flapping frequency of the manta ray (BL = 1.25 m) is much smaller than that of the cownose ray (BL = 0.15 m) is to avoid both the large actuation tension and fluid force density due to the size increase.

A New Muscle Activation Dynamics Model, That Simulates the Calcium Kinetics and Incorporates the Role of Store-Operated Calcium Entry Channels, to Enhance the EMG-Driven Hill-type Models

Hill-type models are frequently used in biomechanical simulations. They are attractive for their low computational cost and close relation to commonly measured musculotendon parameters. Still, more attention is needed to improve the activation dynamics of the model specifically because of the nonlinearity observed in the EMG-Force relation. Moreover, one of the important and practical questions regarding the assessment of the model's performance is how adequately can the model simulate any fundamental type of human movement without modifying model parameters for different tasks? This paper tries to answer this question by proposing a simple physiologically based activation dynamics model. The model describes the ?kinetics of the calcium dynamics while activating and deactivating the muscle contraction process. Hence, it allowed simulating the recently discovered role of store-operated calcium entry (SOCE) channels as immediate counter-flux to calcium loss across the tubular system during excitation-contraction coupling. By comparing the ability to fit experimental data without readjusting the parameters, the proposed model has proven to have more steady performance than phenomenologically based models through different submaximal isometric contraction levels. This model indicates that more physiological insights is key for improving Hill-type model performance.

Energetics of optimal undulatory swimming organisms

Energy consumption is one of the primary considerations in animal locomotion. In swimming locomotion, a number of questions related to swimming energetics of an organism and how the energetic quantities scale with body size remain open, largely due to the difficulties with modeling and measuring the power production and consumption. Based on a comprehensive theoretical framework that incorporates cyclic muscle behavior, structural dynamics and swimming hydrodynamics, we perform extensive computational simulations and show that many of the outstanding problems in swimming energetics can be explained by considering the coupling between hydrodynamics and muscle contraction characteristics, as well as the trade-offs between the conflicting performance goals of sustained swimming speed U and cost of transport COT. Our results lead to three main conclusions: (1) in contrast to previous hypotheses, achieving optimal values of U and COT is independent of producing maximal power or efficiency; (2) muscle efficiency in swimming, in contrast to that in flying or running, decreases with increasing body size, consistent with muscle contraction characteristics; (3) the long-standing problem of two disparate patterns of longitudinal power output distributions in swimming fish can be reconciled by relating the two patterns to U-optimal or COT-optimal swimmers, respectively. We also provide further evidence that the use of tendons in caudal regions is beneficial from an energetic perspective. Our conclusions explain and unify many existing observations and are supported by computational data covering nine orders of magnitude in body size.

3D computational models explain muscle activation patterns and energetic functions of internal structures in fish swimming

How muscles are used is a key to understanding the internal driving of fish swimming. However, the underlying mechanisms of some features of the muscle activation patterns and their differential appearance in different species are still obscure. In this study, we explain the muscle activation patterns by using 3D computational fluid dynamics models coupled to the motion of fish with prescribed deformation and examining the torque and power required along the fish body with two primary swimming modes. We find that the torque required by the hydrodynamic forces and body inertia exhibits a wave pattern that travels faster than the curvature wave in both anguilliform and carangiform swimmers, which can explain the traveling wave speeds of the muscle activations. Notably, intermittent negative power (i.e., power delivered by the fluid to the body) on the posterior part, along with a timely transfer of torque and energy by tendons, explains the decrease in the duration of muscle activation towards the tail. The torque contribution from the body elasticity further clarifies the wave speed increase or the reverse of the wave direction of the muscle activation on the posterior part of a carangiform swimmer. For anguilliform swimmers, the absence of the aforementioned changes in the muscle activation on the posterior part is consistent with our torque prediction and the absence of long tendons from experimental observations. These results provide novel insights into the functions of muscles and tendons as an integral part of the internal driving system, especially from an energy perspective, and they highlight the differences in the internal driving systems between the two primary swimming modes.

For undulatory swimming, fish form posteriorly traveling waves of body bending by activating their muscles sequentially along the body. However, experimental observations have shown that the muscle activation wave does not simply match the bending wave. Researchers have previously computed the torque required for muscles along the body based on classic hydrodynamic theories and explained the higher wave speed of the muscle activation compared to the curvature wave. However, the origins of other features of the muscle activation pattern and their variation among different species are still obscure after decades of research. In this study, we use 3D computational fluid dynamics models to compute the spatiotemporal distributions of both the torque and power required for eel-like and mackerel-like swimming. By examining both the torque and power patterns and considering the energy transfer, storage, and release by tendons and body viscoelasticity, we can explain not only the features and variations in the muscle activation patterns as observed from fish experiments but also how tendons and body elasticity save energy. We provide a mechanical picture in which the body shape, body movement, muscles, tendons, and body elasticity of a mackerel (or similar) orchestrate to make swimming efficient.

Emergence of behavior through morphology: a case study on an octopus inspired manipulator

The complex motion abilities of the Octopus vulgaris have been an intriguing research topic for biologists and roboticists alike. Various studies have been conducted on the underlying control architectures employed by these high dimensional biological organisms. Researchers have attempted to replicate these architectures on robotic systems. Contrary to previous approaches, this study focuses on a robotic system, which is only morphologically similar to the Octopus vulgaris, and how it would behave under different control policies. This sheds light on the underlying optimality principles that these biological systems employ. Open loop control policies are obtained through a trajectory optimization method on a learned forward dynamic model. The motion patterns emerging from variations in morphology and environment were then derived to study the role of the body and environment. Results show that for the specific case of dynamic reaching with a soft appendage, the invariance in motion profile is a fundamental constraint imposed by the morphology and environment, independent from the controller. This suggests how morphological design can simplify stable control even for highly dimensional nonlinear dynamical systems and can provide insights into design of new soft robotic mechanisms.

Adaptive hindlimb split-belt treadmill walking in rats by controlling basic muscle activation patterns via phase resetting

To investigate the adaptive locomotion mechanism in animals, a split-belt treadmill has been used, which has two parallel belts to produce left–right symmetric and asymmetric environments for walking. Spinal cats walking on the treadmill have suggested the contribution of the spinal cord and associated peripheral nervous system to the adaptive locomotion. Physiological studies have shown that phase resetting of locomotor commands involving a phase shift occurs depending on the types of sensory nerves and stimulation timing, and that muscle activation patterns during walking are represented by a linear combination of a few numbers of basic temporal patterns despite the complexity of the activation patterns. Our working hypothesis was that resetting the onset timings of basic temporal patterns based on the sensory information from the leg, especially extension of hip flexors, contributes to adaptive locomotion on the split-belt treadmill. Our hypothesis was examined by conducting forward dynamic simulations using a neuromusculoskeletal model of a rat walking on a split-belt treadmill with its hindlimbs and by comparing the simulated motions with the measured motions of rats.

Signatures of proprioceptive control in Caenorhabditis elegans locomotion

Animal neuromechanics describes the coordinated self-propelled movement of a body, subject to the combined effects of internal neural control and mechanical forces. Here we use a computational model to identify effects of neural and mechanical modulation on undulatory forward locomotion of Caenorhabditis elegans, with a focus on proprioceptively driven neural control. We reveal a fundamental relationship between body elasticity and environmental drag in determining the dynamics of the body and demonstrate the manifestation of this relationship in the context of proprioceptively driven control. By considering characteristics unique to proprioceptive neurons, we predict the signatures of internal gait modulation that contrast with the known signatures of externally or biomechanically modulated gait. We further show that proprioceptive feedback can suppress neuromechanical phase lags during undulatory locomotion, contrasting with well studied advancing phase lags that have long been a signature of centrally generated, feed-forward control.This article is part of a discussion meeting issue 'Connectome to behaviour: modelling C. elegans at cellular resolution'.

The role of curvature feedback in the energetics and dynamics of lamprey swimming: A closed-loop model

Like other animals, lampreys have a central pattern generator (CPG) circuit that activates muscles for locomotion and also adjusts the activity to respond to sensory inputs from the environment. Such a feedback system is crucial for responding appropriately to unexpected perturbations, but it is also active during normal unperturbed steady swimming and influences the baseline swimming pattern. In this study, we investigate different functional forms of body curvature-based sensory feedback and evaluate their effects on steady swimming energetics and kinematics, since little is known experimentally about the functional form of curvature feedback. The distributed CPG is modeled as chains of coupled oscillators. Pairs of phase oscillators represent the left and right sides of segments along the lamprey body. These activate muscles that flex the body and move the lamprey through a fluid environment, which is simulated using a full Navier-Stokes model. The emergent curvature of the body then serves as an input to the CPG oscillators, closing the loop. We consider two forms of feedback, each consistent with experimental results on lamprey proprioceptive sensory receptors. The first, referred to as directional feedback, excites or inhibits the oscillators on the same side, depending on the sign of a chosen gain parameter, and has the opposite effect on oscillators on the opposite side. We find that directional feedback does not affect beat frequency, but does change the duration of muscle activity. The second feedback model, referred to as magnitude feedback, provides a symmetric excitatory or inhibitory effect to oscillators on both sides. This model tends to increase beat frequency and reduces the energetic cost to the lamprey when the gain is high and positive. With both types of feedback, the body curvature has a similar magnitude. Thus, these results indicate that the same magnitude of curvature-based feedback on the CPG with different functional forms can cause distinct differences in swimming performance.

When animals move, they receive sensory inputs, which in turn are used to modulate the movement. Relatively little is known about how these inputs affect performance during steady locomotion. Using a computational model of a swimming lamprey, we investigated two different types of feedback, both consistent with experimental data. Both have strong, but different, effects on swimming speed and energy consumption, suggesting that sensory feedback is crucial not just for responding to perturbations, but also for high performance steady locomotion.

Transition and formation of the torque pattern of undulatory locomotion in resistive force dominated media

In undulatory locomotion, torques along the body are required to overcome external forces from the environment and bend the body. These torques are usually generated by muscles in animals and closely related to muscle activations. In previous studies, researchers observed a single traveling wave pattern of the torque or muscle activation, but the formation of the torque pattern is still not well understood. To elucidate the formation of the torque pattern required by external resistive forces and the transition as kinematic parameters vary, we use simplistic resistive force theory models of self-propelled, steady undulatory locomotors and examine the spatio-temporal variation of the internal torque. We find that the internal torque has a traveling wave pattern with a decreasing speed normalized by the curvature speed as the wave number (the number of wavelengths on the locomotor's body) increases from 0.5 to 1.8. As the wave number increases to 2 and greater values, the torque transitions into a two-wave-like pattern and complex patterns. Using phasor diagram analysis, we reveal that the formation and transitions of the pattern are consequences of the integration and cancellation of force phasors.

IB2d: a Python and MATLAB implementation of the immersed boundary method

The development of fluid-structure interaction (FSI) software involves trade-offs between ease of use, generality, performance, and cost. Typically there are large learning curves when using low-level software to model the interaction of an elastic structure immersed in a uniform density fluid. Many existing codes are not publicly available, and the commercial software that exists usually requires expensive licenses and may not be as robust or allow the necessary flexibility that in house codes can provide. We present an open source immersed boundary software package, IB2d, with full implementations in both MATLAB and Python, that is capable of running a vast range of biomechanics models and is accessible to scientists who have experience in high-level programming environments. IB2d contains multiple options for constructing material properties of the fiber structure, as well as the advection-diffusion of a chemical gradient, muscle mechanics models, and artificial forcing to drive boundaries with a preferred motion.

A review on locomotion robophysics: the study of movement at the intersection of robotics, soft matter and dynamical systems

Discovery of fundamental principles which govern and limit effective locomotion (self-propulsion) is of intellectual interest and practical importance. Human technology has created robotic moving systems that excel in movement on and within environments of societal interest: paved roads, open air and water. However, such devices cannot yet robustly and efficiently navigate (as animals do) the enormous diversity of natural environments which might be of future interest for autonomous robots; examples include vertical surfaces like trees and cliffs, heterogeneous ground like desert rubble and brush, turbulent flows found near seashores, and deformable/flowable substrates like sand, mud and soil. In this review we argue for the creation of a physics of moving systems-a 'locomotion robophysics'-which we define as the pursuit of principles of self-generated motion. Robophysics can provide an important intellectual complement to the discipline of robotics, largely the domain of researchers from engineering and computer science. The essential idea is that we must complement the study of complex robots in complex situations with systematic study of simplified robotic devices in controlled laboratory settings and in simplified theoretical models. We must thus use the methods of physics to examine both locomotor successes and failures using parameter space exploration, systematic control, and techniques from dynamical systems. Using examples from our and others' research, we will discuss how such robophysical studies have begun to aid engineers in the creation of devices that have begun to achieve life-like locomotor abilities on and within complex environments, have inspired interesting physics questions in low dimensional dynamical systems, geometric mechanics and soft matter physics, and have been useful to develop models for biological locomotion in complex terrain. The rapidly decreasing cost of constructing robot models with easy access to significant computational power bodes well for scientists and engineers to engage in a discipline which can readily integrate experiment, theory and computation.

Similarities and Differences for Swimming in Larval and Adult Lampreys

The spinal locomotor networks controlling swimming behavior in larval and adult lampreys may have some important differences. As an initial step in comparing the locomotor systems in lampreys, in larval animals the relative timing of locomotor movements and muscle burst activity were determined and compared to those previously published for adults. In addition, the kinematics for free swimming in larval and adult lampreys was compared in detail for the first time. First, for swimming in larval animals, the neuromechanical phase lag between the onsets or terminations of muscle burst activity and maximum concave curvature of the body increased with increasing distance along the body, similar to that previously shown in adults. Second, in larval lampreys, but not adults, absolute swimming speed (U; mm s(-1)) increased with animal length (L). In contrast, normalized swimming speed (U'; body lengths [bl] s(-1)) did not increase with L in larval or adult animals. In both larval and adult lampreys, U' and normalized wave speed (V') increased with increasing tail-beat frequency. Wavelength and mechanical phase lag did not vary significantly with tail-beat frequency but were significantly different in larval and adult animals. Swimming in larval animals was characterized by a smaller U/V ratio, Froude efficiency, and Strouhal number than in adults, suggesting less efficient swimming for larval animals. In addition, during swimming in larval lampreys, normalized lateral head movements were larger and normalized lateral tail movements were smaller than for adults. Finally, larval animals had proportionally smaller lateral surface areas of the caudal body and fin areas than adults. These differences are well suited for larval sea lampreys that spend most of the time buried in mud/sand, in which swimming efficiency is not critical, compared to adults that would experience significant selection pressure to evolve higher-efficiency swimming to catch up to and attach to fish for feeding as well as engage in long-distance migration during spawning. Finally, the differences in swim efficiency for larval and adult lampreys are compared to other animals employing the anguilliform mode of swimming.

Gait and speed selection in slender inertial swimmers

Inertial swimmers use flexural movements to push water and generate thrust. We quantify this dynamical process for a slender body in a fluid by accounting for passive elasticity and hydrodynamics and active muscular force generation and proprioception. Our coupled elastohydrodynamic model takes the form of a nonlinear eigenvalue problem for the swimming speed and locomotion gait. The solution of this problem shows that swimmers use quantized resonant interactions with the fluid environment to enhance speed and efficiency. Thus, a fish is like an optimized diode that converts a prescribed alternating transverse motion to forward motion. Our results also allow for a broad comparative view of swimming locomotion and provide a mechanistic basis for the empirical relation linking the swimmer's speed U, length L, and tail beat frequency f, given by U/L ~ f [Bainbridge R (1958) J Exp Biol 35:109-133]. Furthermore, we show that a simple form of proprioceptive sensory feedback, wherein local muscle activation is function of body curvature, suffices to drive elastic instabilities associated with thrust production and leads to a spontaneous swimming gait without the need for a central pattern generator. Taken together, our results provide a simple mechanistic view of swimming consistent with natural observations and suggest ways to engineer artificial swimmers for optimal performance.

Strategies for swimming: explorations of the behaviour of a neuro-musculo-mechanical model of the lamprey

Experiments were performed on a neuro-musculo-mechanical model of a lamprey, to explore the strategies for controlling swimming speed. The muscle component of the model was based on previous experiments on isolated lamprey muscle. The patterns of muscle activation were those found in EMG studies on swimming lampreys. The fluid mechanics were modelled with G.I. Taylor's simplification. Tail beat frequencies of 2-6 sec(-1) were combined with muscle activation strengths of 0.1% to 20% of maximum tetanic isometric strength. The resulting forward swimming speed and changing body shape were recorded. From the changing body shape the speed of the backward-travelling wave of curvature was calculated, as well as the ratio between the speeds of the waves of activation and curvature. For any given activation strength there was a tail beat frequency that gave maximal forward speed. Furthermore, for all the combinations of activation strength and tail beat frequency that gave such maximum swimming speeds, the ratio of the speed of the wave of curvature to the wave of muscle activation was approximately 0.75. This is similar to the ratio found in swimming lampreys.

The role of mechanical resonance in the neural control of swimming in fishes

The bodies of many fishes are flexible, elastic structures; if you bend them, they spring back. Therefore, they should have a resonant frequency: a bending frequency at which the output amplitude is maximized for a particular input. Previous groups have hypothesized that swimming at this resonant frequency could maximize efficiency, and that a neural circuit called the central pattern generator might be able to entrain to a mechanical resonance. However, fishes swim in water, which may potentially damp out many resonant effects. Additionally, their bodies are elongated, which means that bending can occur in complicated ways along the length of the body. We review previous studies of the mechanical properties of fish bodies, and then present new data that demonstrate complex bending properties of elongated fish bodies. Resonant peaks in amplitude exist, but there may be many of them depending on the body wavelength. Additionally, they may not correspond to the maximum swimming speed. Next, we describe experiments using a closed-loop preparation of the lamprey, in which a preparation of the spinal cord is linked to a real-time simulation of the muscle and body properties, allowing us to examine resonance entrainment as we vary the simulated resonant frequency. We find that resonance entrainment does occur, but is rare. Gain had a significant, though weak, effect, and a nonlinear muscle model produced resonance entrainment more often than a linear filter. We speculate that resonance may not be a critical effect for efficient swimming in elongate, anguilliform swimmers, though it may be more important for stiffer carangiform and thunniform fishes.

A forced damped oscillation framework for undulatory swimming provides new insights into how propulsion arises in active and passive swimming

A fundamental issue in locomotion is to understand how muscle forcing produces apparently complex deformation kinematics leading to movement of animals like undulatory swimmers. The question of whether complicated muscle forcing is required to create the observed deformation kinematics is central to the understanding of how animals control movement. In this work, a forced damped oscillation framework is applied to a chain-link model for undulatory swimming to understand how forcing leads to deformation and movement. A unified understanding of swimming, caused by muscle contractions (“active” swimming) or by forces imparted by the surrounding fluid (“passive” swimming), is obtained. We show that the forcing triggers the first few deformation modes of the body, which in turn cause the translational motion. We show that relatively simple forcing patterns can trigger seemingly complex deformation kinematics that lead to movement. For given muscle activation, the forcing frequency relative to the natural frequency of the damped oscillator is important for the emergent deformation characteristics of the body. The proposed approach also leads to a qualitative understanding of optimal deformation kinematics for fast swimming. These results, based on a chain-link model of swimming, are confirmed by fully resolved computational fluid dynamics (CFD) simulations. Prior results from the literature on the optimal value of stiffness for maximum speed are explained.

The damped harmonic oscillator framework has been applied to interrogate numerous engineering systems like the tuned mass damper used in power transmission, automobiles, and buildings to reduce vibrations. We apply the same framework to undulatory swimming to understand the emergence of movement due to muscular and/or environmental forcing. It helps elucidate why flexible bodies can indeed be propelled forward by not only the internal muscular forcing but also by external fluid forces as reported earlier in which dead trout were found to swim in the wake of cylinder. We show how forcing triggers the first few deformation modes of the swimmer similar to how the appropriate forcing triggers the fundamental deformation modes on a guitar string. We show how the lower deformation modes produce the best forward propulsion of the body. This insight reveals that swimming is viable for small enough frequencies of neuromuscular activation relative to the natural frequencies of the body and for sufficiently stiff elastic properties. Thus, these results identify the key mechanistic parameters that would have been crucial to the evolutionary emergence of swimming animals. The proposed framework can help understand neural control of movement in swimming as well as control of underwater vehicles.

Emergence of the advancing neuromechanical phase in a resistive force dominated medium

Undulatory locomotion, a gait in which thrust is produced in the opposite direction of a traveling wave of body bending, is a common mode of propulsion used by animals in fluids, on land, and even within sand. As such, it has been an excellent system for discovery of neuromechanical principles of movement. In nearly all animals studied, the wave of muscle activation progresses faster than the wave of body bending, leading to an advancing phase of activation relative to the curvature toward the tail. This is referred to as "neuromechanical phase lags" (NPL). Several multiparameter neuromechanical models have reproduced this phenomenon, but due to model complexity, the origin of the NPL has proved difficult to identify. Here, we use perhaps the simplest model of undulatory swimming to predict the NPL accurately during sand-swimming by the sandfish lizard, with no fitting parameters. The sinusoidal wave used in sandfish locomotion, the friction-dominated and noninertial granular resistive force environment, and the simplicity of the model allow detailed analysis, and reveal the fundamental mechanism responsible for the phenomenon: the combination of synchronized torques from distant points on the body and local traveling torques. This general mechanism should help explain the NPL in organisms in other environments; we therefore propose that sand-swimming could be an excellent system with which to generate and test other neuromechanical models of movement quantitatively. Such a system can also provide guidance for the design and control of robotic undulatory locomotors in complex environments.

Environmental interaction influences muscle activation strategy during sand-swimming in the sandfish lizard Scincus scincus

Animals like the sandfish lizard (Scincus scincus) that live in desert sand locomote on and within a granular medium whose resistance to intrusion is dominated by frictional forces. Recent kinematic studies revealed that the sandfish utilizes a wave of body undulation during swimming. Models predict that a particular combination of wave amplitude and wavelength yields maximum speed for a given frequency, and experiments have suggested that the sandfish targets this kinematic waveform. To investigate the neuromechanical strategy of the sandfish during walking, burial and swimming, here we use high-speed X-ray and visible light imaging with synchronized electromyogram (EMG) recordings of epaxial muscle activity. While moving on the surface, body undulation was not observed and EMG showed no muscle activation. During subsurface sand-swimming, EMG revealed an anterior-to-posterior traveling wave of muscle activation which traveled faster than the kinematic wave. Muscle activation intensity increased as the animal swam deeper into the material but was insensitive to undulation frequency. These findings were in accord with empirical force measurements, which showed that resistance force increased with depth but was independent of speed. The change in EMG intensity with depth indicates that the sandfish targets a kinematic waveform (a template) that models predict maximizes swimming speed and minimizes the mechanical cost of transport as the animal descends into granular media. The differences in the EMG pattern compared with EMG of undulatory swimmers in fluids can be attributed to the friction-dominated intrusion forces of granular media.

Whole-organ cell shape analysis reveals the developmental basis of ascidian notochord taper

Here we use in toto imaging together with computational segmentation and analysis methods to quantify the shape of every cell at multiple stages in the development of a simple organ: the notochord of the ascidian Ciona savignyi. We find that cell shape in the intercalated notochord depends strongly on anterior-posterior (AP) position, with cells in the middle of the notochord consistently wider than cells at the anterior or posterior. This morphological feature of having a tapered notochord is present in many chordates. We find that ascidian notochord taper involves three main mechanisms: Planar Cell Polarity (PCP) pathway-independent sibling cell volume asymmetries that precede notochord cell intercalation; the developmental timing of intercalation, which proceeds from the anterior and posterior towards the middle; and the differential rates of notochord cell narrowing after intercalation. A quantitative model shows how the morphology of an entire developing organ can be controlled by this small set of cellular mechanisms.

Optimal shape and motion of undulatory swimming organisms

Undulatory swimming animals exhibit diverse ranges of body shapes and motion patterns and are often considered as having superior locomotory performance. The extent to which morphological traits of swimming animals have evolved owing to primarily locomotion considerations is, however, not clear. To shed some light on that question, we present here the optimal shape and motion of undulatory swimming organisms obtained by optimizing locomotive performance measures within the framework of a combined hydrodynamical, structural and novel muscular model. We develop a muscular model for periodic muscle contraction which provides relevant kinematic and energetic quantities required to describe swimming. Using an evolutionary algorithm, we performed a multi-objective optimization for achieving maximum sustained swimming speed U and minimum cost of transport (COT)--two conflicting locomotive performance measures that have been conjectured as likely to increase fitness for survival. Starting from an initial population of random characteristics, our results show that, for a range of size scales, fish-like body shapes and motion indeed emerge when U and COT are optimized. Inherent boundary-layer-dependent allometric scaling between body mass and kinematic and energetic quantities of the optimal populations is observed. The trade-off between U and COT affects the geometry, kinematics and energetics of swimming organisms. Our results are corroborated by empirical data from swimming animals over nine orders of magnitude in size, supporting the notion that optimizing U and COT could be the driving force of evolution in many species.

Proprioceptive feedback reinforces centrally generated stepping patterns in the cockroach

The relative importance of sensory input for the production of centrally generated motor patterns is crucial to our understanding of how animals coordinate their body segments to locomote. In legged locomotion, where terrain heterogeneity may require stride-by-stride changes in leg placement, evidence suggests that sensory information is essential for the timing of leg movement. In a previous study we showed that in cockroaches, renowned for rapid and stable running, a coordinated pattern can be elicited from the motor centres driving the different legs in the absence of sensory feedback. In the present paper, we assess the role of movement-related sensory inputs in modifying this central pattern. We studied the effect of spontaneous steps as well as imposed transient and periodic movements of a single intact leg, and demonstrate that, depending on the movement properties, the resulting proprioceptive feedback can significantly modify phase relationships among segmental oscillators of other legs. Our analysis suggests that feedback from front legs is weaker but more phasically precise than from hind legs, selectively transferring movement-related information in a manner that strengthens the inherent rhythmic pattern and modulates local perturbations.

Mechanisms underlying rhythmic locomotion: interactions between activation, tension and body curvature waves

Undulatory animal locomotion arises from three closely related propagating waves that sweep rostrocaudally along the body: activation of segmental muscles by motoneurons (MNs), strain of the body wall, and muscle tension induced by activation and strain. Neuromechanical models that predict the relative propagation speeds of neural/muscle activation, muscle tension and body curvature can reveal crucial underlying control features of the central nervous system and the power-generating mechanisms of the muscle. We provide an analytical explanation of the relative speeds of these three waves based on a model of neuromuscular activation and a model of the body-fluid interactions for leech anguilliform-like swimming. First, we deduced the motoneuron spike frequencies that activate the muscle and the resulting muscle tension during swimming in intact leeches from muscle bending moments. Muscle bending moments were derived from our video-recorded kinematic motion data by our body-fluid interaction model. The phase relationships of neural activation and muscle tension in the strain cycle were then calculated. Our study predicts that the MN activation and body curvature waves have roughly the same speed (the ratio of curvature to MN activation speed ≈0.84), whereas the tension wave travels about twice as fast. The high speed of the tension wave resulting from slow MN activation is explained by the multiplicative effects of MN activation and muscle strain on tension development. That is, the product of two slower waves (activation and strain) with appropriate amplitude, bias and phase can generate a tension wave with twice the propagation speed of the factors. Our study predicts that (1) the bending moment required for swimming is achieved by minimal MN spike frequency, rather than by minimal muscle tension; (2) MN activity is greater in the mid-body than in the head and tail regions; (3) inhibitory MNs not only accelerate the muscle relaxation but also reduce the intrinsic tonus tension during one sector of the swim cycle; and (4) movements of the caudal end are passive during swimming. These predictions await verification or rejection through further experiments on swimming animals.

Spikes alone do not behavior make: why neuroscience needs biomechanics

Neural circuits do not function in isolation; they interact with the physical world, accepting sensory inputs and producing outputs via muscles. Since both these pathways are constrained by physics, the activity of neural circuits can only be understood by considering biomechanics of muscles, bodies, and the exterior world. We discuss how animal bodies have natural stable motions that require relatively little activation or control from the nervous system. The nervous system can substantially alter these motions, by subtly changing mechanical properties such as body or leg stiffness. Mechanics can also provide robustness to perturbations without sensory reflexes. By considering a complete neuromechanical system, neuroscientists and biomechanicians together can provide a more integrated view of neural circuitry and behavior.

Mechanisms underlying rhythmic locomotion: dynamics of muscle activation

We have studied the dynamical properties of tension development in leech longitudinal muscle during swimming. A new method is proposed for modeling muscle properties under functionally relevant conditions where the muscle is subjected to both periodic activation and rhythmic length changes. The 'dual-sinusoid' experiments were conducted on preparations of leech nerve cord and body wall. The longitudinal muscle was activated periodically by injection of sinusoidal currents into an identified motoneuron. Simultaneously, sinusoidal length changes were imposed on the body wall with prescribed phase differences (12 values equally spaced over 2π radians) with respect to the current injection. Through the singular value decomposition of appropriately constructed tension data matrices, the leech muscle was found to have a multiplicative structure in which the tension was expressed as the product of activation and length factors. The time courses of activation and length factors were determined from the tension data and were used to develop component models. The proposed modeling method is a general one and is applicable to contractile elements for which the effects of series elasticity are negligible.

Mechanisms underlying rhythmic locomotion: body-fluid interaction in undulatory swimming

Swimming of fish and other animals results from interactions of rhythmic body movements with the surrounding fluid. This paper develops a model for the body-fluid interaction in undulatory swimming of leeches, where the body is represented by a chain of rigid links and the hydrodynamic force model is based on resistive and reactive force theories. The drag and added-mass coefficients for the fluid force model were determined from experimental data of kinematic variables during intact swimming, measured through video recording and image processing. Parameter optimizations to minimize errors in simulated model behaviors revealed that the resistive force is dominant, and a simple static function of relative velocity captures the essence of hydrodynamic forces acting on the body. The model thus developed, together with the experimental kinematic data, allows us to investigate temporal and spatial (along the body) distributions of muscle actuation, body curvature, hydrodynamic thrust and drag, muscle power supply and energy dissipation into the fluid. We have found that: (1) thrust is generated continuously along the body with increasing magnitude toward the tail, (2) drag is nearly constant along the body, (3) muscle actuation waves travel two or three times faster than the body curvature waves and (4) energy for swimming is supplied primarily by the mid-body muscles, transmitted through the body in the form of elastic energy, and dissipated into the water near the tail.

Intersegmental coordination of cockroach locomotion: adaptive control of centrally coupled pattern generator circuits

Animals’ ability to demonstrate both stereotyped and adaptive locomotor behavior is largely dependent on the interplay between centrally generated motor patterns and the sensory inputs that shape them. We utilized a combined experimental and theoretical approach to investigate the relative importance of CPG interconnections vs. intersegmental afferents in the cockroach: an animal that is renowned for rapid and stable locomotion. We simultaneously recorded coxal levator and depressor motor neurons (MN) in the thoracic ganglia of Periplaneta americana, while sensory feedback was completely blocked or allowed only from one intact stepping leg. In the absence of sensory feedback, we observed a coordination pattern with consistent phase relationship that shares similarities with a double-tripod gait, suggesting central, feedforward control. This intersegmental coordination pattern was then reinforced in the presence of sensory feedback from a single stepping leg. Specifically, we report on transient stabilization of phase differences between activity recorded in the middle and hind thoracic MN following individual front-leg steps, suggesting a role for afferent phasic information in the coordination of motor circuits at the different hemiganglia. Data were further analyzed using stochastic models of coupled oscillators and maximum likelihood techniques to estimate underlying physiological parameters, such as uncoupled endogenous frequencies of hemisegmental oscillators and coupling strengths and directions. We found that descending ipsilateral coupling is stronger than ascending coupling, while left–right coupling in both the meso- and meta-thoracic ganglia appear to be symmetrical. We discuss these results in comparison with recent findings in stick insects that share similar neural and body architectures, and argue that the two species may exemplify opposite extremes of a fast–slow locomotion continuum, mediated through different intersegment coordination strategies.

Interactions between internal forces, body stiffness, and fluid environment in a neuromechanical model of lamprey swimming

Animal movements result from a complex balance of many different forces. Muscles produce force to move the body; the body has inertial, elastic, and damping properties that may aid or oppose the muscle force; and the environment produces reaction forces back on the body. The actual motion is an emergent property of these interactions. To examine the roles of body stiffness, muscle activation, and fluid environment for swimming animals, a computational model of a lamprey was developed. The model uses an immersed boundary framework that fully couples the Navier-Stokes equations of fluid dynamics with an actuated, elastic body model. This is the first model at a Reynolds number appropriate for a swimming fish that captures the complete fluid-structure interaction, in which the body deforms according to both internal muscular forces and external fluid forces. Results indicate that identical muscle activation patterns can produce different kinematics depending on body stiffness, and the optimal value of stiffness for maximum acceleration is different from that for maximum steady swimming speed. Additionally, negative muscle work, observed in many fishes, emerges at higher tail beat frequencies without sensory input and may contribute to energy efficiency. Swimming fishes that can tune their body stiffness by appropriately timed muscle contractions may therefore be able to optimize the passive dynamics of their bodies to maximize peak acceleration or swimming speed.

Reflexes and preflexes: on the role of sensory feedback on rhythmic patterns in insect locomotion

Neuromuscular systems are stabilized and controlled by both feedforward and feedback signals. Feedforward pathways driven by central pattern generators (CPGs), in conjunction with preflexive mechanical reaction forces and nonlinear muscle properties, can produce stable stereotypical gaits. Feedback is nonetheless present in both slow and rapid running, and preflexive mechanisms can join with neural reflexes originating in proprioceptive sensors to yield robust behavior in uncertain environments. Here, we develop a single degree-of-freedom neuromechanical model representing a joint actuated by an agonist/antagonist muscle pair driven by motoneurons and a CPG in a periodic rhythm characteristic of locomotion. We consider two characteristic feedback modes: phasic and tonic. The former encodes states such as position in the timing of individual spikes, while the latter can transmit graded measures of force and other continuous variables as spike rates. We use results from phase reduction and averaging theory to predict phase relationships between CPG and motoneurons in the presence of feedback and compare them with simulations of the neuromechanical model, showing that both phasic and tonic feedback can shift motoneuronal timing and thereby affect joint motions. We find that phase changes in neural activation can cooperate with preflexive displacement and velocity effects on muscle force to compensate for externally applied forces, and that these effects qualitatively match experimental observations in the cockroach.

A new model for force generation by skeletal muscle, incorporating work-dependent deactivation

A model is developed to predict the force generated by active skeletal muscle when subjected to imposed patterns of lengthening and shortening, such as those that occur during normal movements. The model is based on data from isolated lamprey muscle and can predict the forces developed during swimming. The model consists of a set of ordinary differential equations, which are solved numerically. The model's first part is a simplified description of the kinetics of Ca(2+) release from sarcoplasmic reticulum and binding to muscle protein filaments, in response to neural activation. The second part is based on A. V. Hill's mechanical model of muscle, consisting of elastic and contractile elements in series, the latter obeying known physiological properties. The parameters of the model are determined by fitting the appropriate mathematical solutions to data recorded from isolated lamprey muscle activated under conditions of constant length or rate of change of length. The model is then used to predict the forces developed under conditions of applied sinusoidal length changes, and the results compared with corresponding data. The most significant advance of this model is the incorporation of work-dependent deactivation, whereby a muscle that has been shortening under load generates less force after the shortening ceases than otherwise expected. In addition, the stiffness in this model is not constant but increases with increasing activation. The model yields a closer prediction to data than has been obtained before, and can thus prove an important component of investigations of the neural-mechanical-environmental interactions that occur during natural movements.