Exploratory search during directed navigation in C. elegans and Drosophila larva
Figures

Diffusive searching in C. elegans and Drosophila larva.
(A) Sample trajectories (left) from 18 worms under isotropic conditions, for 60 min. Tracks have been shifted to all begin at the same location for clarity. A single track (red) is magnified (right) to show the abrupt changes in crawling orientation, flagged as ‘turns’ (gray circles), with comparatively straight-crawling ‘runs’ in between. Runs are characterized by their duration , and turns by their change in heading . (B) Diffusion in the - and -directions and turning rate reduction over time. Worms diffuse in both directions (left), while the rate of turning events steadily decreases over 3000 s (right); it fits an exponential decay (red) with a 765 s time constant. (C) Diffusion in the - and -directions, splitting the time window into s (left) and s (right). (A’) Sample trajectories from 25 Drosophila larvae navigating under isotropic conditions for 15 min. (B’) Diffusion and turning rate over time for crawling Drosophila. and reach similar values, substantially higher than C. elegans (left), and the turning rate does not show a dramatic drop, instead increasing for the first ~2 min. before stabilizing. (C’) Split and graphs covering the first ( s) and second ( s) halves of the trajectory time, converging to similar values in each. Results for (B,C) are based on seven experiments, with 56 tracks and a total of 1608 turns. Results for (B’,C’) are based on 30 experiments, with 434 tracks and 11,294 turns. The directions of overall population drift for (A,A’) are indicated by green arrows, with the numbers indicating the dimensionless drift velocity, in this case extremely small (see Materials and methods). Error bars for (B,B’) are s.e.m.
-
Figure 1—source data 1
Values and s.e.m. for diffusion coefficient vs. time plots
- https://doi.org/10.7554/eLife.30503.004

Consecutive run durations for individual C. elegans tracks under isotropic conditions.
56 sequences are shown (red), with a blue vertical line indicating the change point, found by maximizing the difference in average below and above a particular time in the experiment. Change point times have a wide range, in many cases occurring near the end of trajectories, indicating a steady rise with high variance, rather than an abrupt transition in mean run duration. We note that in more than half the tracks, the difference in mean before and after the change point is less than twice the standard deviation of the data prior to the change point.

Diffusion and navigation in C. elegans and Drosophila thermotaxis.
(A) Schematic of the apparatus (left), where animals crawl atop an agar substrate while exposed to a 1D linear temperature gradient. Sample thermotaxis trajectories (right) from 18 C. elegans crawling for 60 min and 25 Drosophila crawling for 15 min, where the red dots indicate the starting position for all trajectories. Worms start at 20°C and exhibit negative thermotaxis, moving toward their 15°C cultivation temperature, while larvae start at 17.5°C and crawl away from aversive cold temperatures. The wheel indicates labels for crawling direction ranges, with octant 1 parallel to the gradient heading toward warmer temperatures and octant 5 heading antiparallel towards cooler temperatures. (B) Diffusion over time (left) in the x- and y-directions for C. elegans. throughout the experiment, indicating diminished, but highly significant, diffusion along the navigation direction. The average turning rate (upper right) diminishes over time, as in the Figure 1B,B’ isotropic case; the turning rate ratio (lower right) remains below 1 and nearly constant throughout the experiment. (C) Equilibrium probability distributions in the x (top) and y (bottom) directions (both use the same scale), extrapolated from empirical trajectories using the Markov state model (MSM). The lag time used is 750 s. Red lines are smoothed traces to guide the eye. The free-energy picture of equilibrium conditions (right), to place the analysis in context with the protein folding analysis tools employed here. Lower-free energy corresponds to higher population as , where here . (B’) Diffusion, and turning rates for Drosophila larvae, also showing . The average turning rate stabilizes early, and the turning rate ratio , the primary behavioral modulation underlying thermotaxis, remains nearly constant. (C’) Equilibrium probability distributions for Drosophila larvae (left), and the corresponding free-energy landscape (right), determined from the MSM. Analysis is based on 30 experiments, 131 tracks, and 3061 turns for worms, and 20 experiments, 303 tracks, and 7771 turns for larvae. The directions of overall population drift for (A) are indicated by green arrows, with the numbers indicating the dimensionless drift velocity, approximately 10 times greater than for the isotropic navigation cases (see Materials and methods). Error bars are s.e.m. where shown; in the other cases the error bar size is smaller than the line thickness, and therefore not seen.
-
Figure 2—source data 1
Values and s.e.m. for diffusion coefficient vs. time plots
- https://doi.org/10.7554/eLife.30503.007

Reversability, detailed balance, and sampling intervals in the Markov state model approach.
(A) Analysis of the reversibility of the dynamics. Fluxes in positive (red) and negative (blue dashed) directions are shown. The topmost lines show and ; the lines are almost identical. The lowermost lines show and for ; the lines agree, with residual difference due to statistical noise. The intermediate lines show that increasing the statistics decreases the noise and improves the agreement. (B) Equilibrium probability distribution as a function of lag time used to construct the MSM. With increasing lag time the dynamics becomes more Markovian and the determined equilibrium probabilities converge to the limiting one.

Diffusion and navigation in C. elegans salt chemotaxis.
(A) Sample trajectories (25 worms each) of crawling on an agar substrate with a salt concentration gradient increasing toward the right, under a low-salt concentration baseline (25 mM, left) and high-salt concentration baseline (75 mM, right). The wheel indicates labels for crawling direction ranges, with octant 1 pointing directly to higher salt concentrations and octant 5 directly toward low concentrations. The directions of overall population drift are indicated by green arrows, with the numbers indicating the dimensionless drift velocity, approximately 10 times greater than for the isotropic navigation cases (see Materials and methods). (B) Diffusion over time in the and directions for a low-salt baseline (left) and high-salt baseline (right). Error bars are s.e.m. (C) Average turning rate across all crawling directions (top) decays over time for both low and high baseline salt concentration gradients. The turning rate ratio between octant 5 (toward lower concentration) and octant 1 (toward higher concentration) does not stabilize as clearly as for thermotaxis, likely indicative of a reduced level of movement toward preferred salt conditions and a smaller data set. Error bars are s.e.m. (D) Equilibrium probability distributions for worms in low- (left) and high (right)-salt concentration environments, determined by using the Markov state model (MSM). Analysis is based on 14 experiments, 126 tracks, and 4422 turns for low-salt concentration, and 16 experiments, 166 tracks, and 6159 turns for high-salt concentration.
-
Figure 3—source data 1
Values and s.e.m. for diffusion coefficient vs. time plots
- https://doi.org/10.7554/eLife.30503.009

Isothermal tracking in C. elegans includes diffusion in the -direction.
(A) Experimental trajectories under temperature gradient conditions, with worms initially placed at their cultivation temperature of 15°C. Colors are used for the reader to distinguish between individual tracks. (B) The diffusion constants along the and directions, with dominating, but highly significant. Error bars are s.e.m. (C) MSM-generated equilibrium probabilities along the temperature gradient ( axis) showing the long time scale distribution of the population of worms, and perpendicular to the temperature gradient along isotherms ( axis). In , worms avoid regions with low temperatures, but freely explore regions with higher temperature; in , worms explore the axis with approximately equal probability. Analysis based on 25 experiments, with 688 tracks.
-
Figure 4—source data 1
Values and s.e.m. for diffusion coefficient vs. time plots
- https://doi.org/10.7554/eLife.30503.011
Additional files
-
Source code 1
Scripts used in conjunction with MAGAT Analyzer software
- https://doi.org/10.7554/eLife.30503.012
-
Transparent reporting form
- https://doi.org/10.7554/eLife.30503.013