Direct measurements of drag forces in C. elegans crawling locomotion

A neuromechanical model for Drosophila larval crawling based on physical measurements

BACKGROUND

Animal locomotion requires dynamic interactions between neural circuits, the body (typically muscles), and surrounding environments. While the neural circuitry of movement has been intensively studied, how these outputs are integrated with body mechanics (neuromechanics) is less clear, in part due to the lack of understanding of the biomechanical properties of animal bodies. Here, we propose an integrated neuromechanical model of movement based on physical measurements by taking Drosophila larvae as a model of soft-bodied animals.

RESULTS

We first characterized the kinematics of forward crawling in Drosophila larvae at a segmental and whole-body level. We then characterized the biomechanical parameters of fly larvae, namely the contraction forces generated by neural activity, and passive elastic and viscosity of the larval body using a stress-relaxation test. We established a mathematical neuromechanical model based on the physical measurements described above, obtaining seven kinematic values characterizing crawling locomotion. By optimizing the parameters in the neural circuit, our neuromechanical model succeeded in quantitatively reproducing the kinematics of larval locomotion that were obtained experimentally. This model could reproduce the observation of optogenetic studies reported previously. The model predicted that peristaltic locomotion could be exhibited in a low-friction condition. Analysis of floating larvae provided results consistent with this prediction. Furthermore, the model predicted a significant contribution of intersegmental connections in the central nervous system, which contrasts with a previous study. This hypothesis allowed us to make a testable prediction for the variability in intersegmental connection in sister species of the genus Drosophila.

CONCLUSIONS

We generated a neurochemical model based on physical measurement to provide a new foundation to study locomotion in soft-bodied animals and soft robot engineering.

Efficient sliding locomotion of three-link bodies

We study the efficiency of sliding locomotion for three-link bodies with prescribed joint angle motions. The bodies move with no inertia, under dry (Coulomb) friction that is anisotropic (different in the directions normal and tangent to the links) and directional (different in the forward and backward tangent directions). Friction coefficient space can be partitioned into several regions, each with distinct types of efficient kinematics. These include kinematics resembling lateral undulation with very anisotropic friction, small-amplitude reciprocal kinematics, very large-amplitude kinematics near isotropic friction, and kinematics that are very asymmetric about the flat state. In the two-parameter shape space, zero net rotation for elliptical trajectories occurs mainly with bilateral or antipodal symmetry. These symmetric subspaces have about the same peak efficiency as the full space but with much smaller dimension. Adding the second or third harmonics greatly increases the numbers of local optimal for efficiency, but only modestly increases the peak efficiency. Random ensembles with higher harmonics have efficiency distributions that peak near a certain nonzero value and decay rapidly up to the maximum efficiency. A stochastic optimization algorithm is developed to compute optima with higher harmonics. These are simple closed curves, sharpened versions of the elliptical optima in most cases, and achieve much higher efficiencies mainly for small normal friction. With a linear (viscous) resistance law, the optimal trajectories are similar in much of the friction coefficient space, and relative efficiencies are much lower except with very large normal friction.

A three-dimensional habitat for C. elegans environmental enrichment

As we learn more about the importance of gene-environment interactions and the effects of environmental enrichment, it becomes evident that minimalistic laboratory conditions can affect gene expression patterns and behaviors of model organisms. In the laboratory, Caenorhabditis elegans is generally cultured on two-dimensional, homogeneous agar plates abundantly covered with axenic bacteria culture as a food source. However, in the wild, this nematode thrives in rotting fruits and plant stems feeding on bacteria and small eukaryotes. This contrast in habitat complexity suggests that studying C. elegans in enriched laboratory conditions can deepen our understanding of its fundamental traits and behaviors. Here, we developed a protocol to create three-dimensional habitable scaffolds for trans-generational culture of C. elegans in the laboratory. Using decellularization and sterilization of fruit tissue, we created an axenic environment that can be navigated throughout and where the microbial environment can be strictly controlled. C. elegans were maintained over generations on this habitat, and showed a clear behavioral bias for the enriched environment. As an initial assessment of behavioral variations, we found that dauer populations in scaffolds exhibit high-frequency, complex nictation behavior including group towering and jumping behavior.

Intermittent sliding locomotion of a two-link body

We study the possibility of efficient intermittent locomotion for two-link bodies that slide by changing their interlink angle periodically in time. We find that the anisotropy ratio of the sliding friction coefficients is a key parameter, while solutions have a simple scaling dependence on the friction coefficients' magnitudes. With very anisotropic friction, efficient motions involve coasting in low-drag states, with rapid and asymmetric power and recovery strokes. As the anisotropy decreases, burst-and-coast motions change to motions with long power strokes and short recovery strokes, and roughly constant interlink angle velocity on each. These motions are seen in the spaces of sinusoidal and power-law motions described by two and five parameters, respectively. Allowing the duty cycle to vary greatly increases the motions' efficiency compared to the case of symmetric power and recovery strokes. Allowing further variations in the concavity of the power and recovery strokes improves the efficiency further only when friction is very anisotropic. Near isotropic friction, a variety of optimally efficient motions are found with more complex waveforms. Many of the optimal sinusoidal and power-law motions are similar to those that we find with an optimization search in the space of more general periodic functions (truncated Fourier series). When we increase the resistive force's power-law dependence on velocity, the optimal motions become smoother, slower, and less efficient, particularly near isotropic friction.

Dynamics of pattern formation and emergence of swarming in Caenorhabditis elegans

Many animals collectively form complex patterns to tackle environmental difficulties. Several biological and physical factors, such as animal motility, population densities, and chemical cues, play significant roles in this process. However, very little is known about how sensory information interplays with these factors and controls the dynamics of pattern formation. Here, we study the direct relation between oxygen sensing, pattern formation, and emergence of swarming in active Caenorhabditis elegans aggregates. We find that when thousands of animals gather on food, bacteria-mediated decrease in oxygen level slows down the animals and triggers motility-induced phase separation. Three coupled factors—bacterial accumulation, aerotaxis, and population density—act together and control the entire dynamics. Furthermore, we find that biofilm-forming bacterial lawns including Bacillus subtilis and Pseudomonas aeruginosa strongly alter the collective dynamics due to the limited diffusibility of bacteria. Additionally, our theoretical model captures behavioral differences resulting from genetic variations and oxygen sensitivity.

Modelling the mechanics of exploration in larval Drosophila

The Drosophila larva executes a stereotypical exploratory routine that appears to consist of stochastic alternation between straight peristaltic crawling and reorientation events through lateral bending. We present a model of larval mechanics for axial and transverse motion over a planar substrate, and use it to develop a simple, reflexive neuromuscular model from physical principles. The mechanical model represents the midline of the larva as a set of point masses which interact with each other via damped translational and torsional springs, and with the environment via sliding friction forces. The neuromuscular model consists of: 1. segmentally localised reflexes that amplify axial compression in order to counteract frictive energy losses, and 2. long-range mutual inhibition between reflexes in distant segments, enabling overall motion of the model larva relative to its substrate. In the absence of damping and driving, the mechanical model produces axial travelling waves, lateral oscillations, and unpredictable, chaotic deformations. The neuromuscular model counteracts friction to recover these motion patterns, giving rise to forward and backward peristalsis in addition to turning. Our model produces spontaneous exploration, even though the nervous system has no intrinsic pattern generating or decision making ability, and neither senses nor drives bending motions. Ultimately, our model suggests a novel view of larval exploration as a deterministic superdiffusion process which is mechanistically grounded in the chaotic mechanics of the body. We discuss how this may provide new interpretations for existing observations at the level of tissue-scale activity patterns and neural circuitry, and provide some experimental predictions that would test the extent to which the mechanisms we present translate to the real larva.

We investigate the relationship between brain, body and environment in the exploratory behaviour of fruitfly larva. A larva crawls forward by propagating a wave of compression through its segmented body, and changes its crawling direction by bending to one side or the other. We show first that a purely mechanical model of the larva’s body can produce travelling compression waves, sideways bending, and unpredictable, chaotic motions. For this body to locomote through its environment, it is necessary to add a neuromuscular system to counteract the loss of energy due to friction, and to limit the simultaneous compression of segments. These simple additions allow our model larva to generate life-like forward and backward crawling as well as spontaneous turns, which occur without any direct sensing or control of reorientation. The unpredictability inherent in the larva’s physics causes the model to explore its environment, despite the lack of any neural mechanism for rhythm generation or for deciding when to switch from crawling to turning. Our model thus demonstrates how understanding body mechanics can generate and simplify neurobiological hypotheses as to how behaviour arises.

Micropipette force sensors for in vivo force measurements on single cells and multicellular microorganisms

Measuring forces from the piconewton to millinewton range is of great importance for the study of living systems from a biophysical perspective. The use of flexible micropipettes as highly sensitive force probes has become established in the biophysical community, advancing our understanding of cellular processes and microbial behavior. The micropipette force sensor (MFS) technique relies on measurement of the forces acting on a force-calibrated, hollow glass micropipette by optically detecting its deflections. The MFS technique covers a wide micro- and mesoscopic regime of detectable forces (tens of piconewtons to millinewtons) and sample sizes (micrometers to millimeters), does not require gluing of the sample to the cantilever, and allows simultaneous optical imaging of the sample throughout the experiment. Here, we provide a detailed protocol describing how to manufacture and calibrate the micropipettes, as well as how to successfully design, perform, and troubleshoot MFS experiments. We exemplify our approach using the model nematode Caenorhabditis elegans, but by following this protocol, a wide variety of living samples, ranging from single cells to multicellular aggregates and millimeter-sized organisms, can be studied in vivo, with a force resolution as low as 10 pN. A skilled (under)graduate student can master the technique in ~1-2 months. The whole protocol takes ~1-2 d to finish.

Physical exertion exacerbates decline in the musculature of an animal model of Duchenne muscular dystrophy

Duchenne muscular dystrophy (DMD) is a genetic disorder caused by loss of the protein dystrophin. In humans, DMD has early onset, causes developmental delays, muscle necrosis, loss of ambulation, and death. Current animal models have been challenged by their inability to model the early onset and severity of the disease. It remains unresolved whether increased sarcoplasmic calcium observed in dystrophic muscles follows or leads the mechanical insults caused by the muscle's disrupted contractile machinery. This knowledge has important implications for patients, as potential physiotherapeutic treatments may either help or exacerbate symptoms, depending on how dystrophic muscles differ from healthy ones. Recently we showed how burrowing dystrophic (dys-1) C. elegans recapitulate many salient phenotypes of DMD, including loss of mobility and muscle necrosis. Here, we report that dys-1 worms display early pathogenesis, including dysregulated sarcoplasmic calcium and increased lethality. Sarcoplasmic calcium dysregulation in dys-1 worms precedes overt structural phenotypes (e.g., mitochondrial, and contractile machinery damage) and can be mitigated by reducing calmodulin expression. To learn how dystrophic musculature responds to altered physical activity, we cultivated dys-1 animals in environments requiring high intensity or high frequency of muscle exertion during locomotion. We find that several muscular parameters (e.g., size) improve with increased activity. However, longevity in dystrophic animals was negatively associated with muscular exertion, regardless of effort duration. The high degree of phenotypic conservation between dystrophic worms and humans provides a unique opportunity to gain insight into the pathology of the disease as well as the initial assessment of potential treatment strategies.

From head to tail: a neuromechanical model of forward locomotion in Caenorhabditis elegans

With 302 neurons and a near-complete reconstruction of the neural and muscle anatomy at the cellular level, Caenorhabditis elegans is an ideal candidate organism to study the neuromechanical basis of behaviour. Yet despite the breadth of knowledge about the neurobiology, anatomy and physics of C. elegans, there are still a number of unanswered questions about one of its most basic and fundamental behaviours: forward locomotion. How the rhythmic pattern is generated and propagated along the body is not yet well understood. We report on the development and analysis of a model of forward locomotion that integrates the neuroanatomy, neurophysiology and body mechanics of the worm. Our model is motivated by experimental analysis of the structure of the ventral cord circuitry and the effect of local body curvature on nearby motoneurons. We developed a neuroanatomically grounded model of the head motoneuron circuit and the ventral nerve cord circuit. We integrated the neural model with an existing biomechanical model of the worm's body, with updated musculature and stretch receptors. Unknown parameters were evolved using an evolutionary algorithm to match the speed of the worm on agar. We performed 100 evolutionary runs and consistently found electrophysiological configurations that reproduced realistic control of forward movement. The ensemble of successful solutions reproduced key experimental observations that they were not designed to fit, including the wavelength and frequency of the propagating wave. Analysis of the ensemble revealed that head motoneurons SMD and RMD are sufficient to drive dorsoventral undulations in the head and neck and that short-range posteriorly directed proprioceptive feedback is sufficient to propagate the wave along the rest of the body.This article is part of a discussion meeting issue 'Connectome to behaviour: modelling C. elegans at cellular resolution'.

Three-dimensional simulation of the Caenorhabditis elegans body and muscle cells in liquid and gel environments for behavioural analysis

To better understand how a nervous system controls the movements of an organism, we have created a three-dimensional computational biomechanical model of the Caenorhabditis elegans body based on real anatomical structure. The body model is created with a particle system-based simulation engine known as Sibernetic, which implements the smoothed particle-hydrodynamics algorithm. The model includes an elastic body-wall cuticle subject to hydrostatic pressure. This cuticle is then driven by body-wall muscle cells that contract and relax, whose positions and shape are mapped from C. elegans anatomy, and determined from light microscopy and electron micrograph data. We show that by using different muscle activation patterns, this model is capable of producing C. elegans-like behaviours, including crawling and swimming locomotion in environments with different viscosities, while fitting multiple additional known biomechanical properties of the animal. This article is part of a discussion meeting issue 'Connectome to behaviour: modelling C. elegans at cellular resolution'.

NemaFlex: a microfluidics-based technology for standardized measurement of muscular strength of C. elegans

Muscle strength is a functional measure of quality of life in humans. Declines in muscle strength are manifested in diseases as well as during inactivity, aging, and space travel. With conserved muscle biology, the simple genetic model C. elegans is a high throughput platform in which to identify molecular mechanisms causing muscle strength loss and to develop interventions based on diet, exercise, and drugs. In the clinic, standardized strength measures are essential to quantitate changes in patients; however, analogous standards have not been recapitulated in the C. elegans model since force generation fluctuates based on animal behavior and locomotion. Here, we report a microfluidics-based system for strength measurement that we call 'NemaFlex', based on pillar deflection as the nematode crawls through a forest of pillars. We have optimized the micropillar forest design and identified robust measurement conditions that yield a measure of strength that is independent of behavior and gait. Validation studies using a muscle contracting agent and mutants confirm that NemaFlex can reliably score muscular strength in C. elegans. Additionally, we report a scaling factor to account for animal size that is consistent with a biomechanics model and enables comparative strength studies of mutants. Taken together, our findings anchor NemaFlex for applications in genetic and drug screens, for defining molecular and cellular circuits of neuromuscular function, and for dissection of degenerative processes in disuse, aging, and disease.

The control structure of the nematode Caenorhabditis elegans: Neuro-sensory integration and proprioceptive feedback

We develop a biophysically realistic model of the nematode C. elegans that includes: (i) its muscle structure and activation, (ii) key connectomic activation circuitry, and (iii) a weighted and time-dynamic proprioception. In combination, we show that these model components can reproduce the complex waveforms exhibited in C. elegans locomotive behaviors, chiefly omega turns. This is achieved via weighted, time-dependent suppression of the proprioceptive signal. Though speculative, such dynamics are biologically plausible due to the presence of neuromodulators which have recently been experimentally implicated in the escape response, which includes an omega turn. This is the first integrated neuromechanical model to reveal a mechanism capable of generating the complex waveforms observed in the behavior of C. elegans, thus contributing to a mathematical framework for understanding how control decisions can be executed at the connectome level in order to produce the full repertoire of observed behaviors.

Durotaxis in Nematode Caenorhabditis elegans

Durotaxis is a process where cells are able to sense the stiffness of substrates and preferentially migrate toward stiffer regions. Here, we show that the 1-mm-long nematode, Caenorhabditis elegans are also able to detect the rigidity of underlying substrates and always migrate to regions of higher stiffness. Our results indicate that C. elegans are able to judiciously make a decision to stay on stiffer regions. We found that the, undulation frequency, and wavelength of worms, crawling on surfaces show nonmonotonic behavior with increasing stiffness. A number of control experiments were also conducted to verify whether C. elegans are really able to detect the rigidity of substrates or whether the migration to stiffer regions is due to other factors already reported in the literature. As it is known that bacteria and other single-celled organisms exhibit durotaxis toward stiffer surfaces, we conjecture that durotaxis in C. elegans may be one of the strategies developed to improve their chances of locating food.

Predicting path from undulations for C. elegans using linear and nonlinear resistive force theory

A basic issue in the physics of behaviour is the mechanical relationship between an animal and its surroundings. The model nematode C. elegans provides an excellent platform to explore this relationship due to its anatomical simplicity. Nonetheless, the physics of nematode crawling, in which the worm undulates its body to move on a wet surface, is not completely understood and the mathematical models often used to describe this phenomenon are empirical. We confirm that linear resistive force theory, one such empirical model, is effective at predicting a worm's path from its sequence of body postures for forward crawling, reversing, and turning and for a broad range of different behavioural phenotypes observed in mutant worms. Worms recently isolated from the wild have a higher effective drag anisotropy than the laboratory-adapted strain N2 and most mutant strains. This means the wild isolates crawl with less surface slip, perhaps reflecting more efficient gaits. The drag anisotropies required to fit the observed locomotion data (70  ±  28 for the wild isolates) are significantly larger than the values measured by directly dragging worms along agar surfaces (3-10 in Rabets et al (2014 Biophys. J. 107 1980-7)). A proposed nonlinear extension of the resistive force theory model also provides accurate predictions, but does not resolve the discrepancy between the parameters required to achieve good path prediction and the experimentally measured parameters. We confirm that linear resistive force theory provides a good effective model of worm crawling that can be used in applications such as whole-animal simulations and advanced tracking algorithms, but that the nature of the physical interaction between worms and their most commonly studied laboratory substrate remains unresolved.

The effects of groove height and substrate stiffness on C. elegans locomotion

The physical environment surrounding an animal has a significant impact on its behavior. The nematode Caenorhabditis elegans (C. elegans) has proved to be an excellent choice for understanding the adaptability of organisms crawling on soft surfaces. In this work, we investigate the modulation of C. elegans' behavioral kinematics in response to changes in the stiffness of the substrate and study the effect of grooves incised by the worms on their locomotion speed and efficiency. We measure the height of the grooves created by the animals on surfaces with different rigidity using confocal microscopy. Our results indicate that the kinematic properties of C. elegans, including amplitude (A), wavelength (λ) and frequency (f) of head turns depend strongly on surface properties and the height of the grooves created by them. During crawling, we observe that the animal assumes two distinct shapes depending on the stiffness of substrates. As the stiffness increases, the worm's body shape changes gradually from a 'W' shape, which is characterized by low amplitude curvature to the more common 'S' shape, which is characterized by high amplitude curvature, at intermediate values and back to 'W' on stiffer substrates. Although the efficiency is found to vary monotonically with surface stiffness, the forward velocity shows a non-monotonic behavior with the maximum on a surface, where the animal makes the 'S' shape.

The whole worm: brain-body-environment models of C. elegans

Brain, body and environment are in continuous dynamical interaction, and it is becoming increasingly clear that an animal's behavior must be understood as a product not only of its nervous system, but also of the ongoing feedback of this neural activity through the biomechanics of its body and the ecology of its environment. Modeling has an essential integrative role to play in such an understanding. But successful whole-animal modeling requires an animal for which detailed behavioral, biomechanical and neural information is available and a modeling methodology which can gracefully cope with the constantly changing balance of known and unknown biological constraints. Here we review recent progress on both optogenetic techniques for imaging and manipulating neural activity and neuromechanical modeling in the nematode worm Caenorhabditis elegans. This work demonstrates both the feasibility and challenges of whole-animal modeling.

The nematode C. elegans as a complex viscoelastic fluid

The viscoelastic material properties of the model organism C. elegans were probed with a micropipette deflection technique and modelled with the standard linear solid model. Dynamic relaxation measurements were performed on the millimetric nematode to investigate its viscous characteristics in detail. We show that the internal properties of C. elegans can not be fully described by a simple Newtonian fluid. Instead, a power-law fluid model was implemented and shown to be in excellent agreement with experimental results. The nematode exhibits shear thinning properties and its complex fluid characteristics were quantified. The bending-rate dependence of the internal damping coefficient of C. elegans could affect its gait modulation in different external environments.