Evaluating evasion strategies in zebrafish larvae.

An effective evasion strategy allows prey to survive encounters with predators. Prey are generally thought to escape in a direction that is either random or serves to maximize the minimum distance from the predator. Here, we introduce a comprehensive approach to determine the most likely evasion strategy among multiple hypotheses and the role of biomechanical constraints on the escape response of prey fish. Through a consideration of six strategies with sensorimotor noise and previous kinematic measurements, our analysis shows that zebrafish larvae generally escape in a direction orthogonal to the predator's heading. By sensing only the predator's heading, this orthogonal strategy maximizes the distance from fast-moving predators, and, when operating within the biomechanical constraints of the escape response, it provides the best predictions of prey behavior among all alternatives. This work demonstrates a framework for resolving the strategic basis of evasion in predator-prey interactions, which could be applied to a broad diversity of animals.

Spontaneous phase coordination and fluid pumping in model ciliary carpets.

Ciliated tissues, such as in the mammalian lungs, brains, and reproductive tracts, are specialized to pump fluid. They generate flows by the collective activity of hundreds of thousands of individual cilia that beat in a striking metachronal wave pattern. Despite progress in analyzing cilia coordination, a general theory that links coordination and fluid pumping in the limit of large arrays of cilia remains lacking. Here, we conduct in silico experiments with thousands of hydrodynamically interacting cilia, and we develop a continuum theory in the limit of infinitely many independently beating cilia by combining tools from active matter and classical Stokes flow. We find, in both simulations and theory, that isotropic and synchronized ciliary states are unstable. Traveling waves emerge regardless of initial conditions, but the characteristics of the wave and net flows depend on cilia and tissue properties. That is, metachronal phase coordination is a stable global attractor in large ciliary carpets, even under finite perturbations to cilia and tissue properties. These results support the notion that functional specificity of ciliated tissues is interlaced with the tissue architecture and cilia beat kinematics and open up the prospect of establishing structure to function maps from cilium-level beat to tissue-level coordination and fluid pumping.

Soft skeletons transmit force with variable gearing.

A hydrostatic skeleton allows a soft body to transmit muscular force via internal pressure. A human's tongue, an octopus' arm and a nematode's body illustrate the pervasive presence of hydrostatic skeletons among animals, which has inspired the design of soft engineered actuators. However, there is a need for a theoretical basis for understanding how hydrostatic skeletons apply mechanical work. We therefore modeled the shape change and mechanics of natural and engineered hydrostatic skeletons to determine their mechanical advantage (MA) and displacement advantage (DA). These models apply to a variety of biological structures, but we explicitly consider the tube feet of a sea star and the body segments of an earthworm, and contrast them with a hydraulic press and a McKibben actuator. A helical winding of stiff, elastic fibers around these soft actuators plays a critical role in their mechanics by maintaining a cylindrical shape, distributing forces throughout the structure and storing elastic energy. In contrast to a single-joint lever system, soft hydrostats exhibit variable gearing with changes in MA generated by deformation in the skeleton. We found that this gearing is affected by the transmission efficiency of mechanical work (MA×DA) or, equivalently, the ratio of output to input work. The transmission efficiency changes with the capacity to store elastic energy within helically wrapped fibers or associated musculature. This modeling offers a conceptual basis for understanding the relationship between the morphology of hydrostatic skeletons and their mechanical performance.

Teamwork in the viscous oceanic microscale.

Nutrient acquisition is crucial for oceanic microbes, and competitive solutions to solve this challenge have evolved among a range of unicellular protists. However, solitary solutions are not the only approach found in natural populations. A diverse array of oceanic protists form temporary or even long-lasting attachments to other protists and marine aggregates. Do these planktonic consortia provide benefits to their members? Here, we use empirical and modeling approaches to evaluate whether the relationship between a large centric diatom, Coscinodiscus wailesii, and a ciliate epibiont, Pseudovorticella coscinodisci, provides nutrient flux benefits to the host diatom. We find that fluid flows generated by ciliary beating can increase nutrient flux to a diatom cell surface four to 10 times that of a still cell without ciliate epibionts. This cosmopolitan species of diatom does not form consortia in all environments but frequently joins such consortia in nutrient-depleted waters. Our results demonstrate that symbiotic consortia provide a cooperative alternative of comparable or greater magnitude to sinking for enhancement of nutrient acquisition in challenging environments.

Elastic Snap-Through Instabilities Are Governed by Geometric Symmetries.

Many elastic structures exhibit rapid shape transitions between two possible equilibrium states: umbrellas become inverted in strong wind and hopper popper toys jump when turned inside out. This snap through is a general motif for the storage and rapid release of elastic energy, and it is exploited by many biological and engineered systems from the Venus flytrap to mechanical metamaterials. Shape transitions are known to be related to the type of bifurcation the system undergoes, however, to date, there is no general understanding of the mechanisms that select these bifurcations. Here we analyze numerically and analytically two systems proposed in recent literature in which an elastic strip, initially in a buckled state, is driven through shape transitions by either rotating or translating its boundaries. We show that the two systems are mathematically equivalent, and identify three cases that illustrate the entire range of transitions described by previous authors. Importantly, using reduction order methods, we establish the nature of the underlying bifurcations and explain how these bifurcations can be predicted from geometric symmetries and symmetry-breaking mechanisms, thus providing universal design rules for elastic shape transitions.

A fundamental propulsive mechanism employed by swimmers and flyers throughout the animal kingdom.

Even casual observations of a crow in flight or a shark swimming demonstrate that animal propulsive structures bend in patterned sequences during movement. Detailed engineering studies using controlled models in combination with analysis of flows left in the wakes of moving animals or objects have largely confirmed that flexibility can confer speed and efficiency advantages. These studies have generally focused on the material properties of propulsive structures (propulsors). However, recent developments provide a different perspective on the operation of nature's flexible propulsors, which we consider in this Commentary. First, we discuss how comparative animal mechanics have demonstrated that natural propulsors constructed with very different material properties bend with remarkably similar kinematic patterns. This suggests that ordering principles beyond basic material properties govern natural propulsor bending. Second, we consider advances in hydrodynamic measurements demonstrating suction forces that dramatically enhance overall thrust produced by natural bending patterns. This is a previously unrecognized source of thrust production at bending surfaces that may dominate total thrust production. Together, these advances provide a new mechanistic perspective on bending by animal propulsors operating in fluids - either water or air. This shift in perspective offers new opportunities for understanding animal motion as well as new avenues for investigation into engineered designs of vehicles operating in fluids.

Multisynchrony in Active Microfilaments.

Biological microfilaments exhibit a variety of synchronization modes. Recent experiments observed that a pair of isolated eukaryotic flagella, coupled solely via the fluid medium, display synchrony at nontrivial phase lags in addition to in-phase and antiphase synchrony. Using an elastohydrodynamic filament model in conjunction with numerical simulations and a Floquet-type theoretical analysis, we demonstrate that it is possible to reach multiple synchronization states by varying the intrinsic activity of the filament and the strength of hydrodynamic coupling between the two filaments. Then, we derive an evolution equation for the phase difference between the two filaments at weak coupling, and use a Kuramoto-style phase sensitivity analysis to reveal the nature of the bifurcations underlying the transitions between these different synchronized states.

Motile cilia create fluid-mechanical microhabitats for the active recruitment of the host microbiome.

We show that mucociliary membranes of animal epithelia can create fluid-mechanical microenvironments for the active recruitment of the specific microbiome of the host. In terrestrial vertebrates, these tissues are typically colonized by complex consortia and are inaccessible to observation. Such tissues can be directly examined in aquatic animals, providing valuable opportunities for the analysis of mucociliary activity in relation to bacteria recruitment. Using the squid-vibrio model system, we provide a characterization of the initial engagement of microbial symbionts along ciliated tissues. Specifically, we developed an empirical and theoretical framework to conduct a census of ciliated cell types, create structural maps, and resolve the spatiotemporal flow dynamics. Our multiscale analyses revealed two distinct, highly organized populations of cilia on the host tissues. An array of long cilia ([Formula: see text]25 [Formula: see text]m) with metachronal beat creates a flow that focuses bacteria-sized particles, at the exclusion of larger particles, into sheltered zones; there, a field of randomly beating short cilia ([Formula: see text]10 [Formula: see text]m) mixes the local fluid environment, which contains host biochemical signals known to prime symbionts for colonization. This cilia-mediated process represents a previously unrecognized mechanism for symbiont recruitment. Each mucociliary surface that recruits a microbiome such as the case described here is likely to have system-specific features. However, all mucociliary surfaces are subject to the same physical and biological constraints that are imposed by the fluid environment and the evolutionary conserved structure of cilia. As such, our study promises to provide insight into universal mechanisms that drive the recruitment of symbiotic partners.

Morphological transitions of axially-driven microfilaments.

The interactions of microtubules with motor proteins are ubiquitous in cellular and sub-cellular processes that involve motility and cargo transport. In vitro motility assays have demonstrated that motor-driven microtubules exhibit rich dynamical behaviors from straight to curved configurations. Here, we theoretically investigate the dynamic instabilities of elastic filaments, with free-ends, driven by single follower forces that emulate the action of molecular motors. Using the resistive force theory at low Reynolds number, and a combination of numerical techniques with linear stability analysis, we show the existence of four distinct regimes of filament behavior, including a novel buckled state with locked curvature. These successive instabilities recapitulate the full range of experimentally-observed microtubule behavior, implying that neither structural nor actuation asymmetry are needed to elicit this rich repertoire of motion.

Phoretic and hydrodynamic interactions of weakly confined autophoretic particles.

Phoretic particles self-propel using self-generated physico-chemical gradients at their surface. Within a suspension, they interact hydrodynamically by setting the fluid around them into motion and chemically by modifying the chemical background seen by their neighbours. While most phoretic systems evolve in confined environments due to buoyancy effects, most models focus on their interactions in unbounded flows. Here, we propose a first model for the interaction of phoretic particles in Hele-Shaw confinement and show that in this limit, hydrodynamic and phoretic interactions share not only the same scaling but also the same form, albeit in opposite directions. In essence, we show that phoretic interactions effectively reverse the sign of the interactions that would be obtained for swimmers interacting purely hydrodynamically. Yet, hydrodynamic interactions cannot be neglected as they significantly impact the magnitude of the interactions. This model is then used to analyse the behavior of a suspension. The suspension exhibits swirling and clustering collective modes dictated by the orientational interactions between particles, similar to hydrodynamic swimmers, but here governed by the surface properties of the phoretic particle; the reversal in the sign of the interaction tends to slow down the swimming motion of the particles.

Model of Collective Fish Behavior with Hydrodynamic Interactions.

Fish schooling is often modeled with self-propelled particles subject to phenomenological behavioral rules. Although fish are known to sense and exploit flow features, these models usually neglect hydrodynamics. Here, we propose a novel model that couples behavioral rules with far-field hydrodynamic interactions. We show that (1) a new "collective turning" phase emerges, (2) on average, individuals swim faster thanks to the fluid, and (3) the flow enhances behavioral noise. The results of this model suggest that hydrodynamic effects should be considered to fully understand the collective dynamics of fish.

Activity-induced instability of phonons in 1D microfluidic crystals.

One-dimensional crystals of passively-driven particles in microfluidic channels exhibit collective vibrational modes reminiscent of acoustic 'phonons'. These phonons are induced by the long-range hydrodynamic interactions among the particles and are neutrally stable at the linear level. Here, we analyze the effect of particle activity - self-propulsion - on the emergence and stability of these phonons. We show that the direction of wave propagation in active crystals is sensitive to the intensity of the background flow. We also show that activity couples, at the linear level, transverse waves to the particles' rotational motion, inducing a new mode of instability that persists in the limit of large background flow, or, equivalently, vanishingly small activity. We then report a new phenomenon of phonons switching back and forth between two adjacent crystals in both passively-driven and active systems, similar in nature to the wave switching observed in quantum mechanics, optical communication, and density stratified fluids. These findings could have implications for the design of commercial microfluidic systems and the self-assembly of passive and active micro-particles into one-dimensional structures.

Cortical activity predicts good variation in human motor output.

Human movement patterns have been shown to be particularly variable if many combinations of activity in different muscles all achieve the same task goal (i.e., are goal-equivalent). The nervous system appears to automatically vary its output among goal-equivalent combinations of muscle activity to minimize muscle fatigue or distribute tissue loading, but the neural mechanism of this "good" variation is unknown. Here we use a bimanual finger task, electroencephalography (EEG), and machine learning to determine if cortical signals can predict goal-equivalent variation in finger force output. 18 healthy participants applied left and right index finger forces to repeatedly perform a task that involved matching a total (sum of right and left) finger force. As in previous studies, we observed significantly more variability in goal-equivalent muscle activity across task repetitions compared to variability in muscle activity that would not achieve the goal: participants achieved the task in some repetitions with more right finger force and less left finger force (right > left) and in other repetitions with less right finger force and more left finger force (left > right). We found that EEG signals from the 500 milliseconds (ms) prior to each task repetition could make a significant prediction of which repetitions would have right > left and which would have left > right. We also found that cortical maps of sites contributing to the prediction contain both motor and pre-motor representation in the appropriate hemisphere. Thus, goal-equivalent variation in motor output may be implemented at a cortical level.

Density Shock Waves in Confined Microswimmers.

Motile and driven particles confined in microfluidic channels exhibit interesting emergent behavior, from propagating density bands to density shock waves. A deeper understanding of the physical mechanisms responsible for these emergent structures is relevant to a number of physical and biomedical applications. Here, we study the formation of density shock waves in the context of an idealized model of microswimmers confined in a narrow channel and subject to a uniform external flow. Interestingly, these density shock waves exhibit a transition from "subsonic" with compression at the back to "supersonic" with compression at the front of the population as the intensity of the external flow increases. This behavior is the result of a nontrivial interplay between hydrodynamic interactions and geometric confinement, and it is confirmed by a novel quasilinear wave model that properly captures the dependence of the shock formation on the external flow. These findings can be used to guide the development of novel mechanisms for controlling the emergent density distribution and the average population speed, with potentially profound implications on various processes in industry and biotechnology, such as the transport and sorting of cells in flow channels.

Cooperative transport in sea star locomotion.

It is unclear how animals with radial symmetry control locomotion without a brain. Using a combination of experiments, mathematical modeling, and robotics, we tested the extent to which this control emerges in sea stars (Protoreaster nodosus) from the local control of their hundreds of feet and their mechanical interactions with the body. We discovered that these animals compensate for an experimental increase in their submerged weight by recruiting more feet that synchronize in the power stroke of the locomotor cycle during their bouncing gait. Mathematical modeling of the mechanics of a sea star replicated this response to loading without a central controller. A robotic sea star was found to similarly recruit more actuators under higher loads through purely decentralized control. These results suggest that an array of biological or engineered actuators are capable of cooperative transport where the actuators are dynamically recruited by the mechanics of the body. In particular, the body's vertical oscillations serve to recruit feet in greater numbers to overcome the weight to propel the body forward. This form of distributed control contrasts the conventional view of animal locomotion as governed by the central nervous system and offers inspiration for the design of engineered devices with arrays of actuators.

Flow Physics Explains Morphological Diversity of Ciliated Organs.

Organs that pump fluids by the coordinated beat of motile cilia through the lumen are integral to animal physiology. Such organs include the human airways, brain ventricles, and reproductive tracts. Although cilia organization and duct morphology vary drastically in the animal kingdom, ducts are typically classified as either carpet or flame designs. The reason behind this dichotomy and how duct design relates to fluid pumping remain unclear. Here, we demonstrate that two structural parameters -- lumen diameter and cilia-to-lumen ratio -- organize the observed duct diversity into a continuous spectrum that connects carpets to flames across all animal phyla. Using a unified fluid model, we show that carpet and flame designs maximize flow rate and pressure generation, respectively. We propose that convergence of ciliated organ designs follows functional constraints rather than phylogenetic distance, along with universal design rules for ciliary pumps.

STRUCTURE-FUNCTION RELATIONSHIPS OF MUCOCILIARY CLEARANCE IN HUMAN AIRWAYS.

Our study focuses on the intricate connection between tissue-level organization and ciliated organ function in humans, particularly in understanding the morphological organization of airways and their role in mucociliary clearance. Mucociliary clearance is a key mechanical defense mechanism of human airways, and clearance failure is associated with many respiratory diseases, including chronic obstructive pulmonary disease (COPD) and asthma. While single-cell transcriptomics have unveiled the cellular complexity of the human airway epithelium, our understanding of the mechanics that link epithelial structure to clearance function mainly stem from animal models. This reliance on animal data limits crucial insights into human airway barrier function and hampers the human-relevant in vitro modeling of airway diseases. This study, for the first time, maps the distribution of ciliated and secretory cell types along the airway tree in both rats and humans, noting species-specific differences in ciliary function and elucidates structural parameters of airway epithelia that predict clearance function in both native and in vitro tissues alike. By uncovering how tissue organization influences ciliary function, we can better understand disruptions in mucociliary clearance, which could have implications for various ciliated organs beyond the airways.

STRUCTURE-FUNCTION RELATIONSHIPS OF MUCOCILIARY CLEARANCE IN HUMAN AIRWAYS.

Mucociliary clearance is a key mechanical defense mechanism of human airways, and clearance failure is linked to major respiratory diseases, such as chronic obstructive pulmonary disease (COPD) and asthma. While single-cell transcriptomics have unveiled the cellular complexity of the human airway epithelium, our understanding of the mechanics that link epithelial structure to clearance function mainly stem from animal models. This reliance on animal data limits crucial insights into human airway barrier function and hampers the human-relevant in vitro modeling of airway diseases. Our study fills this crucial knowledge gap and for the first time (1) maps the distribution of ciliated and secretory cell types on the mucosal surface along the proximo-distal axis of the rat and human airway tree, (2) identifies species-specific differences in ciliary beat and clearance function, and (3) elucidates structural parameters of airway epithelia that predict clearance function in both native and in vitro tissues alike. Our broad range of experimental approaches and physics-based modeling translate into generalizable parameters to quantitatively benchmark the human-relevancy of mucociliary clearance in experimental models, and to characterize distinct disease states.

Dynamic behavior of elastic strips near shape transitions.

Elastic strips provide a general motif for studying shape transitions. When actuated through rotation of its boundaries, a buckled strip exhibits, depending on the direction of rotation, three types of shape transitions: buckling, algebraic snap-through, or exponential snap-through. The transition dynamics is linked to the character of the bifurcation, which, in turn, is disclosed by the normal form of the system, but deriving normal forms is challenging. Recent work has used asymptotic methods to obtain this form for algebraic snap-through, but, to date, there is no methodology for extending this analysis to other transitions. Here we introduce a method to analyze the dynamic characteristics of an elastic strip near a transition and extend, in a straightforward manner, the previously proposed asymptotic analysis to exponential snap-through and buckling transitions. Importantly, we show that these normal forms dictate all the dynamic characteristics of the elastic strip near a shape transition. Our analysis provides reliable tools to diagnose and anticipate elastic shape transitions.

Evaluating Evasion Strategies in Zebrafish Larvae.

An effective evasion strategy allows prey to survive encounters with predators. Prey are generally thought to escape in a direction that is either random or serves to maximize the minimum distance from the predator. Here we introduce a comprehensive approach to determine the most likely evasion strategy among multiple hypotheses and the role of biomechanical constraints on the escape response of prey fish. Through a consideration of six strategies with sensorimotor noise and previous kinematic measurements, our analysis shows that zebrafish larvae generally escape in a direction orthogonal to the predator's heading. By sensing only the predator's heading, this orthogonal strategy maximizes the distance from fast-moving predators, and, when operating within the biomechanical constraints of the escape response, it provides the best predictions of prey behavior among all alternatives. This work demonstrates a framework for resolving the strategic basis of evastion in predator-prey interactions, which could be applied to a broad diversity of animals.