1. Computational and Systems Biology
  2. Neuroscience

A stochastic neuronal model predicts random search behaviors at multiple spatial scales in C. elegans

  1. William M Roberts
  2. Steven B Augustine
  3. Kristy J Lawton
  4. Theodore H Lindsay
  5. Tod R Thiele
  6. Eduardo J Izquierdo
  7. Serge Faumont
  8. Rebecca A Lindsay
  9. Matthew Cale Britton
  10. Navin Pokala
  11. Cornelia I Bargmann
  12. Shawn R Lockery  Is a corresponding author
  1. University of Oregon, United States
  2. University of Pennsylvania, United States
  3. Reed College, United States
  4. California Institute of Technology, United States
  5. University of Toronto, Canada
  6. Indiana University, United States
  7. Children's Hospital Los Angeles, United States
  8. University of Minnesota, United States
  9. New York Institute of Technology, United States
  10. Howard Hughes Medical Institute, Rockefeller University, United States
https://doi.org/10.7554/eLife.12572
Share this article
Cite this article
  1. William M Roberts
  2. Steven B Augustine
  3. Kristy J Lawton
  4. Theodore H Lindsay
  5. Tod R Thiele
  6. Eduardo J Izquierdo
  7. Serge Faumont
  8. Rebecca A Lindsay
  9. Matthew Cale Britton
  10. Navin Pokala
  11. Cornelia I Bargmann
  12. Shawn R Lockery
(2016)
A stochastic neuronal model predicts random search behaviors at multiple spatial scales in C. elegans
eLife 5:e12572.
https://doi.org/10.7554/eLife.12572
  2. Download BibTeX
  3. Download .RIS
9 figures, 4 videos and 7 tables

Figures

Figure 1 with 2 supplements
Download asset Open asset
Descriptive statistics of wild type worm tracks.

(A) (x,y)-coordinates of a worm during 10 min of foraging. Inset: Image of a worm showing the black spot (arrow) used for optical tracking (scale bar = 200 µm). (B) The speed distribution computed from the distance moved between successive video frames had a peak at 180 µm/s, which includes both forward and reverse locomotion. A second peak at 14 µm/s corresponds to pauses. The decreased probability of observing speeds <14 µm/s (<0.47 μm/frame) is due to noise in the position measurement. (C) At least three time constants were required to fit (red) the speed autocovariance function (black; grey shading shows ± 1 sem). (D) The worm’s heading remained nearly constant for ~10 s except for a transient peak at 1.4 s (▼), which corresponds to the period of one half cycle of undulation during sinusoidal locomotion. The dashed line shows random reorientation; shading shows ± 1 sem. (E) Example of v(t) showing periods of forward locomotion, reverse locomotion and pauses of various durations. Upward triangles (▲) mark forward-pause-forward (FPF) events; the downward triangle (▼) marks a reverse-pause-reverse (RPR) event. (F) Velocity distributions for the 5 wild type cohorts (5 colors) analyzed in this study. (G) Ensemble-averaged velocity during FPR transitions. All FPR transitions in all wild type cohorts were aligned at the end of forward movement, grouped according to the duration of the pause (2–9 frames), and averaged. Such transitions were defined using a threshold criterion of |v| < 50 μm/s to identify state P (Rakowski et al., 2013). Pauses lasting ≤ 1 frame are not shown because of ambiguity in state identification; pauses lasting ≥ 10 frames are omitted for clarity. (H) Identical to G except RPF transitions are shown. (I) Cumulative probability distributions for dwell time in the pause state defined as in G and H for all FPR and RPF transitions of duration >1 frame in wild type worms.

https://doi.org/10.7554/eLife.12572.003
Figure 1—figure supplement 1
Download asset Open asset
Optical tracking error.

Position data recorded for 1 min each from 4 dead worms spotted by the usual procedure. Under these conditions the stage is stationary. Data from each worm are shown in a different color. The circle encloses 1 standard deviation of the combined 2-D distribution. Optical tracking is more precise than the resolution of the stage position encoder, and thus does not limit the overall resolution of the position measurement.

https://doi.org/10.7554/eLife.12572.004
Figure 1—figure supplement 2
Download asset Open asset
The worm’s search behavior closely resembles a Brownian random walk on time scales longer than 10 s, but not on shorter time scales.

(A) the velocity autocorrelation function (AV) averaged across the 5 wild-type cohorts shows that movements become uncorrelated after ~10 s, primarily as a result of random reorientation during transitions from reverse to forward motion. The period of the damped oscillations in AV corresponds to the period of sinusoidal undulations during locomotion. (B) The observed linear increase in mean-squared distance travelled with time (black; mean ± sem, averaged data from all wild-type worms) shows that on this time scale search behavior approximates a Brownian random walk. (C) At shorter times the observed (black) relation curves upward because worms travel in relatively straight lines during runs. Assuming that the behavior is stationary over the 10 min observation period, the dependence of r2¯ on t can be calculated from the velocity autocorrelation function (red curves in B and C):

r2¯t≅r2t=2∫0t(t−τ)Av(τ)dτ=2∫0t(t−τ)Vτ·Vt0dτ

where ⟨.⟩ denotes statistical expectation (eq. 2.5.12 of Boon and Yip, 1980). Thus, the worm’s movements approximate Brownian motion on a time scale that is longer than the persistence of the velocity autocorrelation, but not at shorter times.

https://doi.org/10.7554/eLife.12572.005
Figure 2 with 3 supplements
Download asset Open asset
Assumptions with supporting data for the Stochastic Switch Model.

(A) Connectivity of forward and reverse command neurons. Arrows with single heads are monosynaptic connections inferred from the C. elegans connectome (White et al., 1986; Varshney et al., 2011) line thickness is proportional to the number of presynaptic specializations seen in the reconstruction of each pairwise connection. Open, double-headed arrows indicate synaptic pathways from or to the indicated pool of neurons outside the network. (B) Voltage recording from the command neuron AVA in the absence of injected current. In this neuron, quasi-stable membrane potentials are seen at -17 and -32 mV. These results differ from previously published AVA recordings, which were made in the presence of hyperpolarizing current (5–10 pA) that kept the membrane potential near -55 mV (Lindsay et al., 2011). (C) Neuronal representation of the Stochastic Switch Model. Forward and reverse command neurons are represented as single binary neuron-like units ℱ and ℛ, respectively. Arrows depicting cross connections (Wℱℛ, Wℛℱ) represent functional (net mono- and polysynaptic) connections between forward and reverse units. Self-connections (Wℱℱ, Wℛℛ) represent synaptic connections between neurons comprising a given unit, voltage dependent currents in these neurons, and polysynaptic recurrent pathways involving non-command neurons. Downward arrows (hℱ, hℛ) represent the combined effects of input from presynaptic neurons, including sensory neurons, and neuromodulation. (D) Markov model representation of the command neuron network. The color of a unit indicates its state of activation (red on, white off). In addition to the forward state F and the reverse state R, there are two pause states, X and Y. Arrows, with their associated rate constants, indicate transitions in which a single unit changes state. Transitions in which two units change state simultaneously have probability zero because single-unit transitions are assumed to be statistically independent. (E) The most likely sequence of states in the hidden Markov model (computed using the Viterbi algorithm) for a representative data segment.

https://doi.org/10.7554/eLife.12572.010
Figure 2—figure supplement 1
Download asset Open asset
Effects of synaptic input on transition rates in stochastic units ℱ and ℛ.

Rate constants are exponential functions of the unit’s net synaptic input S. Excitatory input increases the on rate and decreases the off rate, whereas inhibitory input has the opposite effect. When synaptic input equals zero, the unit switches stochastically between on and off states with rate constant A.

https://doi.org/10.7554/eLife.12572.011
Figure 2—figure supplement 2
Download asset Open asset
Velocity distributions in forward, reverse and pause states.

For each cohort, the velocity distribution g(v) (black; binwidth 2 µm/s) was smoothed by 10 passes of a 1-2-1 binomial smoothing algorithm, then separated into three overlapping velocity distributions: gF(v) (dashed green), gR(v) (dashed blue), gP(v) (dashed red).

For gP(v) we used a Cauchy distribution (half width 18 µm/s) scaled to fit g(0). We estimated gF(v) and gR(v) by subtracting gP(v) from g(v) and restricting the F and R distributions to positive and negative velocities, respectively. The sum gF(v)+gR(v)+gP(v) is shown in solid orange.

https://doi.org/10.7554/eLife.12572.012
Figure 2—figure supplement 3
Download asset Open asset
Cumulative dwell time distributions in states F, R, X and Y.

Comparisons between cumulative dwell time distributions (solid lines) and the exponential distributions predicted by the hidden Markov model (1−exp⁡(−t/dS)gR(v), where dS is the mean dwell time in state S; Table 1). The observed dwell times were tabulated from the most likely sequence of states obtained using the Viterbi algorithm. The origin on the time axis corresponds to one frame.

https://doi.org/10.7554/eLife.12572.013
Figure 3
Download asset Open asset
Fate diagram of the model. 

The system typically cycles clockwise through states F, X, R, Y, with state F frequently interrupted by FXF transitions, leading to state sequences of the form …(FX)nRY(FX)nRY…. Nearly unidirectional transitions out of a given state are shown by red arrows; blue arrows indicate nearly equiprobable transitions. The width of the arrows and the numbers beside them show the probability that the transition out of the state at the tail of the arrow is into the state at the head. The area of each circle is proportional to the probability of the corresponding state (Table 1, column A).

https://doi.org/10.7554/eLife.12572.015
Figure 4
Download asset Open asset
Relationship between pauses and posture.

(A) Average track curvature upon entry in to the pause state in wild type worms. Prior to computing curvature, tracks of individual worms were mirror-imaged as needed such that positive curvature corresponds to a ventral bend. Tracks in the vicinity of pause events were aligned according to the location of the tracking spot in the pause state, converted to curvature, then averaged over all FX transitions (solid blue line; n = 1907), and all RY transitions (red; n = 295) for which the track length was >1.5 mm; shading shows ± 1 S.D. The trace depicts the curvature of the worm posterior to the tracking spot at the end of forward movement (FX transitions) and anterior to the tracking spot at the end of reverse movement (RY transitions). The dashed blue line shows the average curvature at FXR transitions (i.e., excluding FXF stutters). (B) Locomotory phases at which FX transitions occurred, plotted as blue dots on the unit circle. The phase at each FX transition was computed as φ=2πz1/(z2−z1), where z1 and z2 are the positions of the two downward zero crossings of curvature preceding the pause as indicated in panel A, right. The uniform distribution of points around the circle, and therefore the small magnitude of the vector strength (r=0.14; arrow), shows that there was only a small (but statistically significant) phase preference at the end of forward motion (p<10−16; Rayleigh test). (C) Same as B, but for RY transitions. Vector strength is large (r=0.71), indicating a strong tendency to end reverse runs at a particular phase (p<10−63), with a ventral bend in the middle of the body. (D) Average posture at FXR transitions, calculated by integrating the average curvature, computed over all tracks that persisted for >1.5 mm in state F before the pause and >1 mm in state R after the pause. Arrows indicate direction of motion along the track (blue, forward; red, reverse). FXR transitions were typically a simple reversal along the same track. (E) Same as D but for RYF transitions that persisted for >1.5 mm in state R before the pause and >1 mm in state F after the pause. RYF transitions at the end of reverse runs that persisted for >1.5 mm were usually associated with a ventral bend that resulted in a ~180° change of direction as previously described (Gray et al., 2005).

https://doi.org/10.7554/eLife.12572.016
Figure 5
Download asset Open asset
Ablation of command neurons.

(A) Velocity distribution of ablated cohorts (red) compared to sham operated controls (grey) when the indicated command neuron was killed. Stars indicate significant reduction in velocity for the indicated peak (p<0.05 without (⋆) or with (⋆⋆) correction for multiple comparisons; Table 3). (B) Dwell times in F, R, and P in ablated (red) and sham operated animals (grey). Stars indicate significant differences from sham (as defined in Table 4). Horizontal lines indicate the estimated range of d0, the dwell time in the uncoupled state. Each group of ablated animals was tested in parallel with a distinct set of sham operated controls to minimize the effects of variation between populations. Error bars for dwell times are not shown because statistical significance was calculated using the likelihood ratio test (see Table 4 legend), which does not generate sem estimates, and calculation of confidence intervals would have required an excessive amount of computation time. Stars indicate p<0.05 without (⋆) or with (⋆⋆) correction for multiple comparisons (Table 4).

https://doi.org/10.7554/eLife.12572.017
Figure 6
Download asset Open asset
The Stochastic Switch Model correctly predicts the sign and strength of synaptic connections.

(A) Synaptic weights (mean ± sem, n = 5 cohorts) from maximum likelihood fits to velocity data from wild type worms. (B, C) Left, synaptic current in AVB or AVA when the indicated presynaptic neuron was photoactivated (blue line). Right, mean synaptic current during the first 100 ms of the stimulus plotted against holding potential in the postsynaptic neuron (I-V curve). Lines show linear fits to the data at negative holding potentials which were used to estimate vRev. (D), Zero-current holding potential and reversal potential of synaptic currents (mean ± sem) in the indicated postsynaptic neuron (paired t-tests: AVA to AVB, p= 0.043, n = 9; AVB to AVA, p= 0.019, n = 17). (E), Scatter plot of synaptic currents recorded at a holding potential of -15 mV (unpaired t-test: p= 0.010, n ≥ 25).

https://doi.org/10.7554/eLife.12572.021
Figure 7
Download asset Open asset
Predicted and observed effects of HYP and DEP mutations on dwell times.

(A) Predicted effects of changes in membrane potential. (B ) Predicted effects of changes in input resistance. (C) Dwell times in F, R, and P for cohorts of HYP mutants, DEP mutants, and wild type animals. Stars indicate significant change in dwell time (p<0.05 without (⋆) or with (⋆⋆) correction for multiple comparisons; Table 6). In A-C wild type dwell times are indicated by gray bars. Horizontal lines indicate the estimated range of d0 , the dwell time in the uncoupled state. In the ΔV model, h terms were made more negative to model HYP mutants and more positive to model DEP mutants by subtracting or adding a constant Δh = 0.6; qualitatively similar results were obtained for 0 < ∆h ≤ 0.8. In the Δr model, h and w terms were scaled by (1 + f) to model HYP mutants and by (1 - f) to model DEP mutants, with f = 0.6; qualitatively similar results were obtained for 0 < f ≤ 1. Strains, HYP A: DA572 eat-4(ad572); HYP B: MT6308 eat-4(ky5); HYP C: KP4 glr-1(n2461); DEP A: VM1136 lin-15(n765); akIs9 [lin-15(+), Pglr-1::GLR-1(A/T)]; DEP B: VM188 lin-15(n765); akEx52[lin-15(+), Pnmr-1::GLR-1(A/T)].

https://doi.org/10.7554/eLife.12572.023
Figure 8 with 2 supplements
Download asset Open asset
The Stochastic Switch Model accounts for the three main modes of random search in C. elegans.

(A) Plot of mean forward run length versus the weights hℱ and hℛ, illustrating a minimal model of search-scale regulation. (B-H) Calculated effects on search mode of the weights indicated in parentheses. The frequency of reversals (fFPR) is plotted against mF while these three weights are scanned from -6 to 6 weight units in steps of 0.4. Each point was categorized as cropping (magenta), local search (green), ranging (blue), or indeterminate (grey) according to value of fFPR and mF, and whether or not the associated value of mR (not shown) indicated a short or long reversal; see Materials and methods for definitions of search modes. Yellow diamonds mark the scanned points modeled in Figure 8—figure supplement 1. A = 1 Hz; similar results were obtained for A=Amax and A=Amin (Table 7).

https://doi.org/10.7554/eLife.12572.025
Figure 8—figure supplement 1
Download asset Open asset
Simulated worm tracks illustrating cropping, local search, and ranging as defined in the model.

( A-C ) Simulated time is 600 sec with four replicated per panel, each in a different color. A=Amax; similar results were obtained for A=Amin.

https://doi.org/10.7554/eLife.12572.026
Figure 8—figure supplement 2
Download asset Open asset
Extension of the Stochastic Switch Model to chemotaxis.

(A) Circuitry. Behavioral state (F, R, or P) was determined by a modified version of the Stochastic Switch Model in which on and off chemosensory neurons regulated the values of the inputs to the network. During movement up the gradient, the activation states of the on and off cells were set to 1 and 0, respectively, such that hℱ(t)=hℱ+Δhℱ and hℛ(t)=hℛ+Δhℛ, where hℱ and hℛ are the values of the inputs to the network during local search (Table 2). Conversely, during movement down the gradient, the on and off activation states, and the signs of Δh, were reversed. In the tracks shown, ΔhF and Δhℛ were ± 2.6, the value that optimized chemotaxis performance given the speed of the model worm and standard deviation of the gradient. (B) Simulated chemotaxis. The concentration gradient of chemical attractant was modeled as a two dimensional Gaussian (std. dev. = 1.6 cm) originating at the center of a circular arena. Similar tracks were obtained across the full range of values of A, the fundamental time scale of the model.

https://doi.org/10.7554/eLife.12572.027
Figure 9
Download asset Open asset
Extension of the Stochastic Switch Model to deterministic behaviors.

(A) Three functional circuit motifs for deterministic escape behavior initiated by the nociceptive neuron ASH. (B) Predicted steady-state probability of reversal behavior in the resting state and the activated state of the three motifs shown in A. Plotted values are means across the five wild type cohorts shown in Figure 1F. Error bars are ± sem. Numbers in parenthesis are predicted mean first latency to a reversal response. (C) Left, synaptic current in AVB when ASH was photoactivated (blue line). Right, mean synaptic current during the first 100 ms of the stimulus plotted against holding potential in AVB. The line is fit to the data at negative holding potentials. (D) Mean zero-current holding potential and mean reversal potential of synaptic currents (± sem) in AVB (paired t-test: p= 0.013, n = 4).

https://doi.org/10.7554/eLife.12572.029

Videos

Video 1
Download asset
Forward-Pause-Forward transition.

The worm is crawling on a foodless agar plate. The microscope stage moves continuously to keep the tracking spot near the center of the frame. Stage movement can be assessed by monitoring the white streaks in the background, which are segments of the worm's track at earlier times. Behavioral state is indicated in the upper left corner of the frame. The indicated behavioral transition is shown at normal speed, and slowed down by a factor of 5. The worm is paused when the tracking spot is stationary relative to the streaks.

https://doi.org/10.7554/eLife.12572.006
Video 2
Download asset
Reverse-Pause-Reverse transition.
https://doi.org/10.7554/eLife.12572.007
Video 3
Download asset
Forward-Pause-Reverse transition.
https://doi.org/10.7554/eLife.12572.008
Video 4
Download asset
Reverse-Pause-Forward transition.
https://doi.org/10.7554/eLife.12572.009

Tables

Table 1
Maximum likelihood fits of transition rates in wild type C. elegans.

Each cohort was fitted separately; values are expressed as mean ± sem (n = 5 cohorts). Data from wild type cohorts were obtained on the same days as the experimental cohorts for which they served as controls (Tables 3 and 4), but experimental cohorts in this study were separated by weeks to months. All transition rates were constrained to be ≥0. Transition rates that were calculated using the synaptic constraints (Equation 35) are shaded orange; other constrained values are shaded grey. Mean dwell times and state probabilities were calculated from the transition rates. Column A shows fits using the standard model, which has 8 rate constants with two synaptic constraints, resulting in 6 free parameters that determine the 6 synaptic weights (Figure 2C,D; Materials and methods Equations 31–35). Column B shows fits to a model that has only one pause state (X); this model was derived from the standard model by imposing two more constraints: aFY=aRY≅0, yielding 4 free parameters. To allow comparison of models A and B by the likelihood ratio test, which requires tht model B be a special case of model A, aRY and aFY were set slightly >0 (10-10 s-1), thereby avoiding infinite values for aYF and aYR when applying the synaptic constraints, while maintaining a vanishingly small probability of being in state Y (pY< 10-18). The loge likelihood (summed over the 5 cohorts) for model B was 1854 less than for model A, with 30 degrees of freedom for model A (6 per cohort × 5 cohorts) and 20 degrees of freedom for model B (4 per cohort × 5 cohorts). Applying the likelihood ratio test, the difference was highly significant (p<10-100; p=Chi-squared(2L, df), where L  = 1854 and df  = 30–20 =10. Model C is the most general 3-state (F, R, P) model, which allows all six transitions between the three states. The fitted transition rates for model C were nearly identical to model B. Likelihood values are relative to model A.

https://doi.org/10.7554/eLife.12572.014
ABC
2 pause states6 free parameters1 pause state4 free parameters1 pause state6 free parameters
Δ loge likelihood0-1854-1836
Degrees of freedom302030
mean ± sem (n  = 5)mean ± sem (n  = 5)mean ± sem (n  = 5)
aXR (s-1)1.201 ± 0.0991.019 ± 0.0851.008 ± 0.090
aXF (s-1)1.115 ± 0.0871.915 ± 0.1521.914 ± 0.152
aRX (s-1)0.025 ± 0.0080.507 ± 0.0130.507 ± 0.013
aRY (s-1)0.490 ± 0.01510-10
aFX (s-1)0.182 ± 0.0070.198 ± 0.0090.196 ± 0.008
aFY (s-1)0.007 ± 0.00210-10
aYR (s-1)0.411 ± 0.019>109
aYF (s-1)4.575 ± 0.533>109
aFR (s-1)0.001 ± 0.001
aRF (s-1)0.000 ± 0.000
dF (s)5.329 ± 0.2455.096 ± 0.2355.135 ± 0.227
dR (s)1.945 ± 0.0431.975 ± 0.0491.976 ± 0.049
dX (s)0.441 ± 0.0320.349 ± 0.0260.351 ± 0.027
dY (s)0.208 ± 0.019<10-9
pF0.762 ± 0.0150.7641 ± 0.0150.764 ± 0.014
pR0.158 ± 0.0070.158 ± 0.0070.155 ± 0.007
pX0.063 ± 0.0060.081 ± 0.0080.080 ± 0.008
pY0.017 ± 0.002<10-18
Table 2
Synaptic weights derived from the transition rate constants.

The rate constants were taken from Table 1, column A. Two values of the fundamental switching time, A, corresponding to the minimum (0.40 Hz) and maximum (0.86 Hz) values consistent with the ablation results were used in Materials and methods, Equations 36–38 to calculate the corresponding synaptic weights.

https://doi.org/10.7554/eLife.12572.018
 A= 0.4 Hzmean ± sem (n = 5)A= 0.86 Hzmean ± sem (n = 5)
hℱ1.01 ± 0.080.25 ± 0.08
hℛ1.09 ± 0.080.32 ± 0.08
wℱℛ-5.40 ± 0.43-5.40 ± 0.43
Wℛℱ-0.81 ± 0.06-0.81 ± 0.06
Wℱℱ-0.22 ± 0.061.31 ± 0.06
Wℱℛ1.90 ± 0.333.43 ± 0.33
Table 3
Effects of command neuron ablations on undulation frequency, forward velocity and reverse velocity.

Values were computed separately for each worm and are shown as mean ± sem (n = 19–29 ). Undulation frequency was estimated as one-half of the reciprocal of the time of the first local minimum in the heading autocorrelation function. All p-values are from two-tailed t-tests and are shown without correction for multiple comparisons. Blue denotes significance at p<0.05. Red denotes significance at p<0.05 after Bonferroni correction for 15 comparisons.

https://doi.org/10.7554/eLife.12572.019
Undulation frequency (Hz)Forward velocity (μm/s)Reverse velocity (μm/s)
NeuronClassShamAblatedp <ShamAblatedp <ShamAblatedp <
AVBForward0.355 ± 0.0090.230 ± 0.0077x10-11236 ± 6109 ± 45x10-20-327 ± 7-302 ± 80.04
PVCForward0.283 ± 0.0110.290 ± 0.0100.5187 ± 7192 ± 70.7-253 ± 8-248 ± 60.7
AVDReverse0.270 ± 0.0080.236 ± 0.0080.009173 ± 6141 ± 50.0002-243 ± 4-229 ± 50.06
AVAReverse0.302 ± 0.0050.254 ± 0.0094x10-5195 ± 5155 ± 74x10-5-293 ± 7-69 ± 33x10-22
AVEReverse0.264 ± 0.0070.256 ± 0.0080.6165 ± 4160 ± 50.5-235 ± 4-211 ± 60.003
Table 4
Effects of command neuron ablations on model parameters.

The sign of the change (Δ) caused by the ablation is shown as “+” if the value moved away from 0, “–” if the value moved towards 0. Significance was determined using the likelihood ratio test (Weisstein, Eric W. "Likelihood Ratio." From MathWorld--A Wolfram Web Resource. http://mathworld.wolfram. com/LikelihoodRatio.html), which is based on the reduction in likelihood caused by constraining one of the parameters to have the same value in both the ablated cohort and the corresponding sham cohort. The unconstrained fit thus had 12 free parameters (6 for each of the 2 cohorts being compared), while the constrained fit had 11 free parameters. For example, to test the significance of the change in the mean dwell time in the pause state (dP=(pXdX+pYdY)/(pX+pY)dp=(pxdx+pYdY)/(px+py) caused by ablation of the AVA neuron pair, two cohorts (ablated and sham) were grown and tested under identical conditions. The ln likelihood with 12 free parameters was found to be 894794.075. When dp was constrained to be the same for both cohorts, the ln likelihood for the 11 parameter fit was found to be 894784.676. The test statistic D=2×(894794.075−894784.676)=18.798 was assumed to come from a chi-squared distribution with one degree of freedom, which yielded p=1.45×10−5 (shown in the table as p<10−4). The constrained fitting process was repeated in turn for each ablation/sham pair for each of the 9 rows shown in the table. All p-values are shown without correction for multiple comparisons. Blue denotes significance at p<0.05. Red denotes significance at p<0.05 after Bonferroni correction for 27 comparisons

https://doi.org/10.7554/eLife.12572.020
 REVERSEFORWARD
 AVEAVDAVAAVBPVC
 ShamAblate∆p<ShamAblate∆p<ShamAblate∆p<ShamAblate∆p<ShamAblate∆p<
dF (s)5.4555.221–0.25.1584.081–10-156.7303.143–10-997.2892.642–10-996.0586.558+0.02
dR (s)3.0192.436–10-62.5402.367–0.052.3591.243–10-412.1271.681–10-62.8422.396–0.0005
dX (s)0.5480.548–10.5140.520+0.60.4800.582+10-70.3700.437+10-60.4570.508+0.004
dY (s)0.2290.263+0.0020.2290.241+0.070.2140.331+10-70.1970.144–10-70.2200.226+0.5
dp (s)0.4950.496+10.4600.468+0.50.4370.510+10-40.3310.416+10-110.4100.457+0.003
pF0.7200.747+0.0040.7230.689–0.00050.8090.704–10-240.8180.745–10-150.7490.787+0.0002
pR0.1920.158–10-40.1880.203+0.050.1220.137+0.040.1290.120–0.20.1810.139–10-9
pX0.0730.078+0.070.0720.088+10-70.0580.113+10-90.0410.125+10-990.0560.061+0.2
pY0.0140.017+0.020.0170.020+0.0050.0110.046+10-230.0120.010–0.010.0140.013–0.5
Table 5
Effects of mutations on mean undulation frequency, mean forward velocity and mean reverse velocity.

Values were computed separately for each worm and are shown as mean ± sem (n = 25–31). Undulation frequency was estimated as one-half of the reciprocal of the time of the first local minimum in the heading autocorrelation function. All p-values are from two-tailed t-tests and are shown without correction for multiple comparisons. Blue denotes significance at p<0.05. Red denotes significance at p<0.05 after Bonferroni correction for 15 comparisons.

https://doi.org/10.7554/eLife.12572.022
Undulation frequency (Hz)Forward velocity (μm/s)Reverse velocity (μm/s)
GenotypeClassWild typeMutantp <Wild typeMutantp <Wild typeMutantp <
eat-4(ad572)HYP A0.272 ± 0.0110.222 ± 0.0074x10-4156 ± 5122 ± 41x10-5-228 ± 5-236 ± 90.5
eat-4(ky5)HYP B0.317 ± 0.0110.256 ± 0.0092x10-4184 ± 7143 ± 65x10-5-262 ± 10-271 ± 70.5
glr-1(n2461)HYP C0.294 ± 0.0080.291 ± 0.0100.9158 ± 5166 ± 60.3-236 ± 6-236 ± 51
glr-1::glr-1(A/T)DEP A0.272 ± 0.0110.642 ± 0.0296x10-13156 ± 5112 ± 53x10-7-228 ± 5-143 ± 52x10-15
nmr-1::glr-1(A/T)DEP B0.294 ± 0.0080.695 ± 0.0372x10-12158 ± 5138 ± 50.011-236 ± 6-144 ± 57x10-15
Table 6
Effects of mutations on model parameters.

Significance was determined using the likelihood ratio test as described in Table 4. The sign of the change (Δ) caused by the mutation is shown as “+” if the value moved away from 0, “–” if the value moved towards 0. All p-values are shown without correction for multiple comparisons. Blue denotes significance at p<0.05. Red denotes significance at p<0.05 after Bonferroni correction for 27 comparisons.

https://doi.org/10.7554/eLife.12572.024
HYPDEP
HYP A: eat-4(ad572)HYP B: eat-4(ky5)HYP C: glr-1(n2461)DEP A: glr-1::glr-1(A/T)DEP B: nmr-1::glr-1(A/T)
ControlMutant∆p <ControlMutant∆p<ControlMutant∆p<ControlMutant∆p<ControlMutant∆p<
dF (s)4.7719.564+10-874.9568.643+10-645.1817.871+10-334.7710.940–10-995.1810.742–10-99
dR (s)2.0432.821+10-71.9102.769+10-122.0183.004+10-162.0450.875–10-992.0180.709–10-9 9
dX (s)0.4811.040+0.0050.5290.844+10-430.4590.727+10-390.4820.328–10-490.4600.235–10-99
dY (s)0.2470.382+10-50.2380.290+0.0050.2210.164–10-50.2470.097–10-930.2210.079–10-99
dp (s)0.4280.982+10-990.4660.793+10-520.4090.677+10-430.4280.286–10-610.4090.204–10-99
pF0.7290.839+10-270.7340.832+10-260.7550.785+0.0030.7280.410–10-990.7550.407–10-99
pR0.1770.062–10-370.1670.079–10-290.1610.135–10-40.1770.389+10-990.1610.404+10-99
pX0.0730.090+10-50.0770.081+0.30.0670.073+0.030.0730.164+10-990.0670.151+10-99
pY0.0220.009–10-130.0210.008–10-250.0180.007–10-270.0220.037+10-210.0180.037+10-41
Table 7
Regulation of search mode.

The weights in each subspace were scanned from -6 to 6 weight units in steps of 0.4 with A=Amin or A=Amax. The letter x means that the indicated search mode was present for at least one point in the subspace when A=Amin and when A=Amax; the letters y and z mean that the mode was present only when A=Amin or A=Amax, respectively. See Materials and methods for definitions of search modes. Control-point weights as defined by the theoretical relationship between weights and search scale (Equation 5) are shown in bold. Only the three-weight subspaces are sufficient for producing all three search modes and full coverage of the plane defined by reversal frequency and mF plane as shown in Figure 8.

https://doi.org/10.7554/eLife.12572.028
SubspaceCroppingDwellingRangingCoverage
hℱhℱ x ⋅ ⋅ ⋅
hℛ xx⋅ ⋅ ⋅
wFF xx⋅ ⋅ ⋅
wℛℛ x . . .
wℛℱ x . . .
wℱℛ x . . .
hℱ, wℛℛyxx. . .
hℱ, wℱℛ xx. . .
hℱ, wℛℱ xx. . .
hℛ, wℛℛ xx. . .
hℛ, wℱℛxxx. . .
hℛ, wℛℱ xx. . .
wFF, wℛℱ xx. . .
wFF, wℱℛ xx. . .
wFF, wℛℛ xx. . .
hℛ, wℱℱ xx. . .
hℱ, hℛ xx. . .
hℱ, wℱℱ xx. . .
wℛℛ, wℱℛxx .......
wℛℛ, wℛℱ y . . .
wℱℛ, wℛℱ x .......
hℱ, wℛℱ, wℱℛxxx.......
hℛ, wℛℱ, wℱℛxxx.......
wFF, wℛℱ, wℱℛxxx.......
hℱ, wℛℛ, wℱℛxxx.......
hℛ, wℛℛ, wℱℛxxx.......
wFF, wℛℛ, wℱℛxxx.......
hℱ, wℛℛ, wℛℱxxx.......
hℛ, wℛℛ, wℛℱzxx.......
wFF, wℛℛ, wℛℱxxx.......
hℱ, hℛ, wℱℛxxx.......
hℱ, hℛ, wℛℛxxx.......
hℱ, hℛ, wℛℱxxx.......
hℱ, wFF, wℛℛxxx.......
hℱ, wFF, wℱℛ xx.......
hℱ, wFF, wℛℱ xx.......
hℛ, wFF, wℱℛxxx.......
hℛ, wFF, wℛℱxxx.......
hℛ, wFF, wℛℛyxx.......
hℱ, hℛ, wFFxxx.......
wℛℛ, wℛℱ, wℱℛxx .......

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. William M Roberts
  2. Steven B Augustine
  3. Kristy J Lawton
  4. Theodore H Lindsay
  5. Tod R Thiele
  6. Eduardo J Izquierdo
  7. Serge Faumont
  8. Rebecca A Lindsay
  9. Matthew Cale Britton
  10. Navin Pokala
  11. Cornelia I Bargmann
  12. Shawn R Lockery
(2016)
A stochastic neuronal model predicts random search behaviors at multiple spatial scales in C. elegans
eLife 5:e12572.
https://doi.org/10.7554/eLife.12572