Shifts in a single muscle's control potential of body dynamics are determined by mechanical feedback

A hierarchical model for external electrical control of an insect, accounting for inter-individual variation of muscle force properties

Cyborg control of insect movement is promising for developing miniature, high-mobility, and efficient biohybrid robots. However, considering the inter-individual variation of the insect neuromuscular apparatus and its neural control is challenging. We propose a hierarchical model including inter-individual variation of muscle properties of three leg muscles involved in propulsion (retractor coxae), joint stiffness (pro- and retractor coxae), and stance-swing transition (protractor coxae and levator trochanteris) in the stick insect Carausius morosus. To estimate mechanical effects induced by external muscle stimulation, the model is based on the systematic evaluation of joint torques as functions of electrical stimulation parameters. A nearly linear relationship between the stimulus burst duration and generated torque was observed. This stimulus-torque characteristic holds for burst durations of up to 500ms, corresponding to the stance and swing phase durations of medium to fast walking stick insects. Hierarchical Bayesian modeling revealed that linearity of the stimulus-torque characteristic was invariant, with individually varying slopes. Individual prediction of joint torques provides significant benefits for precise cyborg control.

The force response of muscles to activation and length perturbations depends on length history

Recent studies have demonstrated that muscle force is not determined solely by activation under dynamic conditions, and that length history has an important role in determining dynamic muscle force. Yet, the mechanisms for how muscle force is produced under dynamic conditions remain unclear. To explore this, we investigated the effects of muscle stiffness, activation and length perturbations on muscle force. First, submaximal isometric contraction was established for whole soleus muscles. Next, the muscles were actively shortened at three velocities. During active shortening, we measured muscle stiffness at optimal muscle length (L0) and the force response to time-varying activation and length perturbations. We found that muscle stiffness increased with activation but decreased as shortening velocity increased. The slope of the relationship between maximum force and activation amplitude differed significantly among shortening velocities. Also, the intercept and slope of the relationship between length perturbation amplitude and maximum force decreased with shortening velocity. As shortening velocities were related to muscle stiffness, the results suggest that length history determines muscle stiffness and the history-dependent muscle stiffness influences the contribution of activation and length perturbations to muscle force. A two-parameter viscoelastic model including a linear spring and a linear damper in parallel with measured stiffness predicted history-dependent muscle force with high accuracy. The results and simulations support the hypothesis that muscle force under dynamic conditions can be accurately predicted as the force response of a history-dependent viscoelastic material to length perturbations.

Cockroaches adjust body and appendages to traverse cluttered large obstacles

To traverse complex terrain, animals often transition between locomotor modes. It is well-known that locomotor transitions can be induced by switching in neural control circuits or driven by a need to minimize metabolic energetic cost. Recent work discovered that locomotor transitions in complex 3-D terrain cluttered with large obstacles can emerge from physical interaction with the environment controlled by the nervous system. For example, to traverse cluttered, stiff grass-like beams, the discoid cockroach often transitions from using a strenuous pitch mode pushing across to using a less strenuous roll mode rolling into and through the gaps, and this transition requires overcoming a potential energy barrier. Previous robotic physical modeling demonstrated that kinetic energy fluctuation of body oscillation from self-propulsion can help overcome the barrier and facilitate this transition. However, the animal was observed to transition even when the barrier still exceeded kinetic energy fluctuation. Here, we further studied whether and how the cockroach makes active adjustments to facilitate this transition to traverse cluttered beams. The animal repeatedly flexed its head and abdomen, reduced hind leg sprawl, and depressed one hind leg and elevated the other during the pitch-to-roll transition, which were absent when running on a flat ground. Using a refined potential energy landscape with additional degrees of freedom to model these adjustments, we found that head flexion did not substantially reduce the transition barrier, whereas leg sprawl reduction did so dramatically. We speculate that head flexion is for sensing the terrain to guide the transition via sensory feedback control.

Flight muscle power increases with strain amplitude and decreases with cycle frequency in zebra finches (Taeniopygia guttata)

Birds that use high flapping frequencies can modulate aerodynamic force by varying wing velocity, which is primarily a function of stroke amplitude and wingbeat frequency. Previous measurements from zebra finches (Taeniopygia guttata) flying across a range of speeds in a wind tunnel demonstrate that although the birds modulated both wingbeat kinematic parameters, they exhibited greater changes in stroke amplitude. These two kinematic parameters contribute equally to aerodynamic force, so the preference for modulating amplitude over frequency may instead derive from limitations of muscle physiology at high frequency. We tested this hypothesis by developing a novel in situ work loop approach to measure muscle force and power output from the whole pectoralis major of zebra finches. This method allowed for multiple measurements over several hours without significant degradation in muscle power. We explored the parameter space of stimulus, strain amplitude and cycle frequencies measured previously from zebra finches, which revealed overall high net power output of the muscle, despite substantial levels of counter-productive power during muscle lengthening. We directly compared how changes to muscle shortening velocity via strain amplitude and cycle frequency affected muscle power. Increases in strain amplitude led to increased power output during shortening with little to no change in power output during lengthening. In contrast, increases in cycle frequency did not lead to increased power during shortening but instead increased counter-productive power during lengthening. These results demonstrate why at high wingbeat frequency, increasing wing stroke amplitude could be a more effective mechanism to cope with increased aerodynamic demands.

In vitro virtual reality: an anatomically explicit musculoskeletal simulation powered by in vitro muscle using closed-loop tissue-software interaction

Muscle force-length dynamics are governed by intrinsic contractile properties, motor stimulation and mechanical load. Although intrinsic properties are well characterised, physiologists lack in vitro instrumentation to account for combined effects of limb inertia, musculoskeletal architecture and contractile dynamics. We introduce in vitro virtual reality (in vitro-VR) which enables in vitro muscle tissue to drive a musculoskeletal jumping simulation. In hardware, muscle force from a frog plantaris was transmitted to a software model where joint torques, inertia and ground reaction forces were computed to advance the simulation at 1 kHz. To close the loop, simulated muscle strain was returned to update in vitro length. We manipulated (1) stimulation timing and (2) the virtual muscle's anatomical origin. This influenced interactions among muscular, inertial, gravitational and contact forces dictating limb kinematics and jump performance. We propose that in vitro-VR can be used to illustrate how neuromuscular control and musculoskeletal anatomy influence muscle dynamics and biomechanical performance.

Nanometer-scale structure differences in the myofilament lattice spacing of two cockroach leg muscles correspond to their different functions

Muscle is highly organized across multiple length scales. Consequently, small changes in the arrangement of myofilaments can influence macroscopic mechanical function. Two leg muscles of a cockroach have identical innervation, mass, twitch responses, length-tension curves and force-velocity relationships. However, during running, one muscle is dissipative (a 'brake'), while the other dissipates and produces significant positive mechanical work (bifunctional). Using time-resolved X-ray diffraction in intact, contracting muscle, we simultaneously measured the myofilament lattice spacing, packing structure and macroscopic force production of these muscles to test whether structural differences in the myofilament lattice might correspond to the muscles' different mechanical functions. While the packing patterns are the same, one muscle has 1 nm smaller lattice spacing at rest. Under isometric stimulation, the difference in lattice spacing disappeared, consistent with the two muscles' identical steady-state behavior. During periodic contractions, one muscle undergoes a 1 nm greater change in lattice spacing, which correlates with force. This is the first identified structural feature in the myofilament lattice of these two muscles that shares their whole-muscle dynamic differences and quasi-static similarities.

History-dependent perturbation response in limb muscle

Muscle mediates movement but movement is typically unsteady and perturbed. Muscle is known to behave non-linearly and with history-dependent properties during steady locomotion, but the importance of history dependence in mediating muscle function during perturbations remains less clear. To explore the capacity of muscles to mitigate perturbations during locomotion, we constructed a series of perturbations that varied only in kinematic history, keeping instantaneous position, velocity and time from stimulation constant. We found that the response of muscle to a perturbation is profoundly history dependent, varying 4-fold as baseline frequency changes, and dissipating energy equivalent to ∼6 times the kinetic energy of all the limbs in 5 ms (nearly 2400 W kg^-1). Muscle energy dissipation during a perturbation is predicted primarily by the force at the onset of the perturbation. This relationship holds across different frequencies and timings of stimulation. This history dependence behaves like a viscoelastic memory producing perturbation responses that vary with the frequency of the underlying movement.

Information-based centralization of locomotion in animals and robots

The centralization of locomotor control from weak and local coupling to strong and global is hard to assess outside of particular modeling frameworks. We developed an empirical, model-free measure of centralization that compares information between control signals and both global and local states. A second measure, co-information, quantifies the net redundancy in global and local control. We first validate that our measures predict centralization in simulations of phase-coupled oscillators. We then test how centralization changes with speed in freely running cockroaches. Surprisingly, across all speeds centralization is constant and muscle activity is more informative of the global kinematic state (the averages of all legs) than the local state of that muscle's leg. Finally we use a legged robot to show that mechanical coupling alone can change the centralization of legged locomotion. The results of these systems span a design space of centralization and co-information for biological and robotic systems.

Body stiffness and damping depend sensitively on the timing of muscle activation in lampreys

Unlike most manmade machines, animals move through their world using flexible bodies and appendages, which bend due to internal muscle and body forces, and also due to forces from the environment. Fishes in particular must cope with fluid dynamic forces that not only resist their overall swimming movements but also may have unsteady flow patterns, vortices, and turbulence, many of which occur more rapidly than what the nervous system can process. Has natural selection led to mechanical properties of fish bodies and their component tissues that can respond very quickly to environmental perturbations? Here, we focus on the mechanical properties of isolated muscle tissue and of the entire intact body in the silver lamprey, Ichthyomyzon unicuspis. We developed two modified work loop protocols to determine the effect of small perturbations on the whole body and on isolated segments of muscle as a function of muscle activation and phase within the swimming cycle. First, we examined how the mechanical properties of the whole lamprey body change depending on the timing of muscle activity. Relative to passive muscle, muscle activation can modulate the effective stiffness by about two-fold and modulate the effective damping by >10-fold depending on the activation phase. Next, we performed a standard work loop test on small sections of axial musculature while adding low-amplitude sinusoidal perturbations at specific frequencies. We modeled the data using a new system identification technique based on time-periodic system analysis and harmonic transfer functions (HTFs) and used the resulting models to predict muscle function under novel conditions. We found that the effective stiffness and damping of muscle varies during the swimming cycle, and that the timing of activation can alter both the magnitude and timing of peak stiffness and damping. Moreover, the response of the isolated muscle was highly nonlinear and length dependent, but the body's response was much more linear. We applied the resulting HTFs from our experiments to explore the effect of pairs of antagonistic muscles. The results suggest that when muscles work against each other as antagonists, the combined system has weaker nonlinearities than either muscle segment alone. Together, these results begin to provide an integrative understanding of how activation timing can tune the mechanical response properties of muscles, enabling fish to swim effectively in their complex and unpredictable environment.

Millisecond Spike Timing Codes for Motor Control

Millisecond variations in spiking patterns can radically alter motor behavior, suggesting that traditional rate-based theories of motor control require revision. The importance of spike timing in sensorimotor control arises from dynamic interactions between the nervous system, muscles, and the body. New mechanisms, model systems, and theories are revealing how these interactions shape behavior.

Endogenous rhythm and pattern-generating circuit interactions in cockroach motor centres

Cockroaches are rapid and stable runners whose gaits emerge from the intricate, and not fully resolved, interplay between endogenous oscillatory pattern-generating networks and sensory feedback that shapes their rhythmic output. Here we studied the endogenous motor output of a brainless, deafferented preparation. We monitored the pilocarpine-induced rhythmic activity of levator and depressor motor neurons in the mesothoracic and metathoracic segments in order to reveal the oscillatory networks’ architecture and interactions. Data analyses included phase relations, latencies between and overlaps of rhythmic bursts, spike frequencies, and the dependence of these parameters on cycle frequency. We found that, overall, ipsilateral connections are stronger than contralateral ones. Our findings revealed asymmetries in connectivity among the different ganglia, in which meta-to-mesothoracic ascending coupling is stronger than meso-to-metathoracic descending coupling. Within-ganglion coupling between the metathoracic hemiganglia is stronger than that in the mesothoracic ganglion. We also report differences in the role and mode of operation of homologue network units (manifested by levator and depressor nerve activity). Many observed characteristics are similar to those exhibited by intact animals, suggesting a dominant role for feedforward control in cockroach locomotion. Based on these data we posit a connectivity scheme among components of the locomotion pattern generating system.

Summary: Detailed analysis of fictive motor patterns unveils endogenous characteristics of the cockroach thoracic locomotion control networks and their interrelations and enables an explanatory connectivity model.

How do treadmill speed and terrain visibility influence neuromuscular control of guinea fowl locomotion?

Locomotor control mechanisms must flexibly adapt to both anticipated and unexpected terrain changes to maintain movement and avoid a fall. Recent studies revealed that ground birds alter movement in advance of overground obstacles, but not treadmill obstacles, suggesting context-dependent shifts in the use of anticipatory control. We hypothesized that differences between overground and treadmill obstacle negotiation relate to differences in visual sensory information, which influence the ability to execute anticipatory manoeuvres. We explored two possible explanations: (1) previous treadmill obstacles may have been visually imperceptible, as they were low contrast to the tread, and (2) treadmill obstacles are visible for a shorter time compared with runway obstacles, limiting time available for visuomotor adjustments. To investigate these factors, we measured electromyographic activity in eight hindlimb muscles of the guinea fowl (Numida meleagris, N=6) during treadmill locomotion at two speeds (0.7 and 1.3 m s⁻¹) and three terrain conditions at each speed: (i) level, (ii) repeated 5 cm low-contrast obstacles (<10% contrast, black/black), and (iii) repeated 5 cm high-contrast obstacles (>90% contrast, black/white). We hypothesized that anticipatory changes in muscle activity would be higher for (1) high-contrast obstacles and (2) the slower treadmill speed, when obstacle viewing time is longer. We found that treadmill speed significantly influenced obstacle negotiation strategy, but obstacle contrast did not. At the slower speed, we observed earlier and larger anticipatory increases in muscle activity and shifts in kinematic timing. We discuss possible visuomotor explanations for the observed context-dependent use of anticipatory strategies.

Summary: Guinea fowl (Numida meleagris) show speed-dependent shifts in neuromuscular control during obstacle negotiation, characterized by a greater reliance on anticipatory modulation and stride-to-stride neural adjustments at slow speed, shifting towards feedforward activation and intrinsic mechanical stability at high speed.

Intersegmental coupling and recovery from perturbations in freely running cockroaches

Cockroaches are remarkably stable runners, exhibiting rapid recovery from external perturbations. To uncover the mechanisms behind this important behavioral trait, we recorded leg kinematics of freely running animals in both undisturbed and perturbed trials. Functional coupling underlying inter-leg coordination was monitored before and during localized perturbations, which were applied to single legs via magnetic impulses. The resulting transient effects on all legs and the recovery times to normal pre-perturbation kinematics were studied. We estimated coupling architecture and strength by fitting experimental data to a six-leg-unit phase oscillator model. Using maximum-likelihood techniques, we found that a network with nearest-neighbor inter-leg coupling best fitted the data and that, although coupling strengths vary among preparations, the overall inputs entering each leg are approximately balanced and consistent. Simulations of models with different coupling strengths encountering perturbations suggest that the coupling schemes estimated from our experiments allow animals relatively fast and uniform recoveries from perturbations.

Feedback control as a framework for understanding tradeoffs in biology

Control theory arose from a need to control synthetic systems. From regulating steam engines to tuning radios to devices capable of autonomous movement, it provided a formal mathematical basis for understanding the role of feedback in the stability (or change) of dynamical systems. It provides a framework for understanding any system with regulation via feedback, including biological ones such as regulatory gene networks, cellular metabolic systems, sensorimotor dynamics of moving animals, and even ecological or evolutionary dynamics of organisms and populations. Here, we focus on four case studies of the sensorimotor dynamics of animals, each of which involves the application of principles from control theory to probe stability and feedback in an organism's response to perturbations. We use examples from aquatic (two behaviors performed by electric fish), terrestrial (following of walls by cockroaches), and aerial environments (flight control by moths) to highlight how one can use control theory to understand the way feedback mechanisms interact with the physical dynamics of animals to determine their stability and response to sensory inputs and perturbations. Each case study is cast as a control problem with sensory input, neural processing, and motor dynamics, the output of which feeds back to the sensory inputs. Collectively, the interaction of these systems in a closed loop determines the behavior of the entire system.

Instantaneous kinematic phase reflects neuromechanical response to lateral perturbations of running cockroaches

Instantaneous kinematic phase calculation allows the development of reduced-order oscillator models useful in generating hypotheses of neuromechanical control. When perturbed, changes in instantaneous kinematic phase and frequency of rhythmic movements can provide details of movement and evidence for neural feedback to a system-level neural oscillator with a time resolution not possible with traditional approaches. We elicited an escape response in cockroaches (Blaberus discoidalis) that ran onto a movable cart accelerated laterally with respect to the animals' motion causing a perturbation. The specific impulse imposed on animals (0.50 [Formula: see text] 0.04 m s[Formula: see text]; mean, SD) was nearly twice their forward speed (0.25 [Formula: see text] 0.06 m s[Formula: see text]. Instantaneous residual phase computed from kinematic phase remained constant for 110 ms after the onset of perturbation, but then decreased representing a decrease in stride frequency. Results from direct muscle action potential recordings supported kinematic phase results in showing that recovery begins with self-stabilizing mechanical feedback followed by neural feedback to an abstracted neural oscillator or central pattern generator. Trials fell into two classes of forward velocity changes, while exhibiting statistically indistinguishable frequency changes. Animals pulled away from the side with front and hind legs of the tripod in stance recovered heading within 300 ms, whereas animals that only had a middle leg of the tripod resisting the pull did not recover within this period. Animals with eight or more legs might be more robust to lateral perturbations than hexapods.

Abdicating power for control: a precision timing strategy to modulate function of flight power muscles

Muscles driving rhythmic locomotion typically show strong dependence of power on the timing or phase of activation. This is particularly true in insects' main flight muscles, canonical examples of muscles thought to have a dedicated power function. However, in the moth (Manduca sexta), these muscles normally activate at a phase where the instantaneous slope of the power-phase curve is steep and well below maximum power. We provide four lines of evidence demonstrating that, contrary to the current paradigm, the moth's nervous system establishes significant control authority in these muscles through precise timing modulation: (i) left-right pairs of flight muscles normally fire precisely, within 0.5-0.6 ms of each other; (ii) during a yawing optomotor response, left-right muscle timing differences shift throughout a wider 8 ms timing window, enabling at least a 50 per cent left-right power differential; (iii) timing differences correlate with turning torque; and (iv) the downstroke power muscles alone causally account for 47 per cent of turning torque. To establish (iv), we altered muscle activation during intact behaviour by stimulating individual muscle potentials to impose left-right timing differences. Because many organisms also have muscles operating with high power-phase gains (Δ(power)/Δ(phase)), this motor control strategy may be ubiquitous in locomotor systems.

Abdicating power for control: a precision timing strategy to modulate function of flight power muscles

Muscles driving rhythmic locomotion typically show strong dependence of power on the timing or phase of activation. This is particularly true in insects' main flight muscles, canonical examples of muscles thought to have a dedicated power function. However, in the moth (Manduca sexta), these muscles normally activate at a phase where the instantaneous slope of the power-phase curve is steep and well below maximum power. We provide four lines of evidence demonstrating that, contrary to the current paradigm, the moth's nervous system establishes significant control authority in these muscles through precise timing modulation: (i) left-right pairs of flight muscles normally fire precisely, within 0.5-0.6 ms of each other; (ii) during a yawing optomotor response, left-right muscle timing differences shift throughout a wider 8 ms timing window, enabling at least a 50 per cent left-right power differential; (iii) timing differences correlate with turning torque; and (iv) the downstroke power muscles alone causally account for 47 per cent of turning torque. To establish (iv), we altered muscle activation during intact behaviour by stimulating individual muscle potentials to impose left-right timing differences. Because many organisms also have muscles operating with high power-phase gains (Δ(power)/Δ(phase)), this motor control strategy may be ubiquitous in locomotor systems.

Temperature gradients drive mechanical energy gradients in the flight muscle of Manduca sexta

A temperature gradient throughout the dominant flight muscle (dorsolongitudinal muscle, DLM(1)) of the hawkmoth Manduca sexta, together with temperature-dependent muscle contractile rates, demonstrates that significant spatial variation in power production is possible within a single muscle. Using in situ work-loop analyses under varying muscle temperatures and phases of activation, we show that regional differences in muscle temperature will induce a spatial gradient in the mechanical power output throughout the DLM(1). Indeed, we note that this power gradient spans from positive to negative values across the predicted temperature range. Warm ventral subunits produce positive power at their in vivo operating temperatures, and therefore act as motors. Concurrently, as muscle temperature decreases dorsally, the subunits produce approximately zero mechanical power output, acting as an elastic energy storage source, and negative power output, behaving as a damper. Adjusting the phase of activation further influences the temperature sensitivity of power output, significantly affecting the mechanical power output gradient that is expressed. Additionally, the separate subregions of the DLM(1) did not appear to employ significant physiological compensation for the temperature-induced differences in power output. Thus, although the components of a muscle are commonly thought to operate uniformly, a significant within-muscle temperature gradient has the potential to induce a mechanical power gradient, whereby subunits within a muscle operate with separate and distinct functional roles.

Lessons for circuit function from large insects: towards understanding the neural basis of motor flexibility

Motor behaviors result from information processing that occurs in multiple neural networks acting at all levels from the initial selection of the behavior to its final generation. A long-standing research interest is how single neural networks can help generate different motor behaviors, that is, the origin of motor flexibility. Modern experimental techniques allow studying neural network activity during the production of multiple motor behaviors. Recent data provide strong evidence that the neural networks controlling insect legs are individually modified in task-dependent and finely tuned fashions. Understanding the mechanistic basis of these neural network modifications will be of particular interest in the upcoming years.

A bio-robotic platform for integrating internal and external mechanics during muscle-powered swimming

To explore the interplay between muscle function and propulsor shape in swimming animals, we built a robotic foot to mimic the morphology and hind limb kinematics of Xenopus laevis frogs. Four foot shapes ranging from low aspect ratio (AR = 0.74) to high (AR = 5) were compared to test whether low-AR feet produce higher propulsive drag force resulting in faster swimming. Using feedback loops, two complementary control modes were used to rotate the foot: force was transmitted to the foot either from (1) a living plantaris longus (PL) muscle stimulated in vitro or (2) an in silico mathematical model of the PL. To mimic forward swimming, foot translation was calculated in real time from fluid force measured at the foot. Therefore, bio-robot swimming emerged from muscle-fluid interactions via the feedback loop. Among in vitro-robotic trials, muscle impulse ranged from 0.12 ± 0.002 to 0.18 ± 0.007 N s and swimming velocities from 0.41 ± 0.01 to 0.43 ± 0.00 m s(-1), similar to in vivo values from prior studies. Trends in in silico-robotic data mirrored in vitro-robotic observations. Increasing AR caused a small (∼10%) increase in peak bio-robot swimming velocity. In contrast, muscle force-velocity effects were strongly dependent on foot shape. Between low- and high-AR feet, muscle impulse increased ∼50%, while peak shortening velocity decreased ∼50% resulting in a ∼20% increase in net work. However, muscle-propulsion efficiency (body center of mass work/muscle work) remained independent of AR. Thus, we demonstrate how our experimental technique is useful for quantifying the complex interplay among limb morphology, muscle mechanics and hydrodynamics.

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.

Mechanics, modulation and modelling: how muscles actuate and control movement

Animal movement is often complex, unsteady and variable. The critical role of muscles in animal movement has captivated scientists for over 300 years. Despite this, emerging techniques and ideas are still shaping and advancing the field. For example, sonomicrometry and ultrasound techniques have enhanced our ability to quantify muscle length changes under in vivo conditions. Robotics and musculoskeletal models have benefited from improved computational tools and have enhanced our ability to understand muscle function in relation to movement by allowing one to simulate muscle-tendon dynamics under realistic conditions. The past decade, in particular, has seen a rapid advancement in technology and shifts in paradigms related to muscle function. In addition, there has been an increased focus on muscle function in relation to the complex locomotor behaviours, rather than relatively simple (and steady) behaviours. Thus, this Theme Issue will explore integrative aspects of muscle function in relation to diverse locomotor behaviours such as swimming, jumping, hopping, running, flying, moving over obstacles and transitioning between environments. Studies of walking and running have particular relevance to clinical aspects of human movement and sport. This Theme Issue includes contributions from scientists working on diverse taxa, ranging from humans to insects. In addition to contributions addressing locomotion in various taxa, several manuscripts will focus on recent advances in neuromuscular control and modulation during complex behaviours. Finally, some of the contributions address recent advances in biomechanical modelling and powered prostheses. We hope that our comprehensive and integrative Theme Issue will form the foundation for future work in the fields of neuromuscular mechanics and locomotion.

Shifts in a single muscle's control potential of body dynamics are determined by mechanical feedback

Muscles are multi-functional structures that interface neural and mechanical systems. Muscle work depends on a large multi-dimensional space of stimulus (neural) and strain (mechanical) parameters. In our companion paper, we rewrote activation to individual muscles in intact, behaving cockroaches (Blaberus discoidalis L.), revealing a specific muscle's potential to control body dynamics in different behaviours. Here, we use those results to provide the biologically relevant parameters for in situ work measurements. We test four hypotheses about how muscle function changes to provide mechanisms for the observed control responses. Under isometric conditions, a graded increase in muscle stress underlies its linear actuation during standing behaviours. Despite typically absorbing energy, this muscle can recruit two separate periods of positive work when controlling running. This functional change arises from mechanical feedback filtering a linear increase in neural activation into nonlinear work output. Changing activation phase again led to positive work recruitment, but at different times, consistent with the muscle's ability to also produce a turn. Changes in muscle work required considering the natural sequence of strides and separating swing and stance contributions of work. Both in vivo control potentials and in situ work loops were necessary to discover the neuromechanical coupling enabling control.

A single muscle's multifunctional control potential of body dynamics for postural control and running

A neuromechanical approach to control requires understanding how mechanics alters the potential of neural feedback to control body dynamics. Here, we rewrite activation of individual motor units of a behaving animal to mimic the effects of neural feedback without concomitant changes in other muscles. We target a putative control muscle in the cockroach, Blaberus discoidalis (L.), and simultaneously capture limb and body dynamics through high-speed videography and a micro-accelerometer backpack. We test four neuromechanical control hypotheses. We supported the hypothesis that mechanics linearly translates neural feedback into accelerations and rotations during static postural control. However, during running, the same neural feedback produced a nonlinear acceleration control potential restricted to the vertical plane. Using this, we reject the hypothesis from previous work that this muscle acts primarily to absorb energy from the body. The conversion of the control potential is paralleled by nonlinear changes in limb kinematics, supporting the hypothesis that significant mechanical feedback filters the graded neural feedback for running control. Finally, we insert the same neural feedback signal but at different phases in the dynamics. In this context, mechanical feedback enables turning by changing the timing and direction of the accelerations produced by the graded neural feedback.