Calcium dynamics regulating the timing of decision-making in C. elegans

Abstract
eLife digest
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

Brains regulate behavioral responses with distinct timings. Here we investigate the cellular and molecular mechanisms underlying the timing of decision-making during olfactory navigation in Caenorhabditis elegans. We find that, based on subtle changes in odor concentrations, the animals appear to choose the appropriate migratory direction from multiple trials as a form of behavioral decision-making. Through optophysiological, mathematical and genetic analyses of neural activity under virtual odor gradients, we further find that odor concentration information is temporally integrated for a decision by a gradual increase in intracellular calcium concentration ([Ca²⁺]_i), which occurs via L-type voltage-gated calcium channels in a pair of olfactory neurons. In contrast, for a reflex-like behavioral response, [Ca²⁺]_i rapidly increases via multiple types of calcium channels in a pair of nociceptive neurons. Thus, the timing of neuronal responses is determined by cell type-dependent involvement of calcium channels, which may serve as a cellular basis for decision-making.

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

eLife digest

Animals use information from their environment to make decisions, like where to go, what to eat, and with whom to mate. This information may be changing or confusing, and decisions may be quick when the sensory information is clear, or slower when the sensory clues are muddled. Scientists often study this kind of decision-making in monkeys and rodents, but it can be hard to pinpoint the exact decision-making mechanisms because these animals have hundreds of millions neurons in their brains.

Studying the mechanisms that underlie decision-making can be easier in a simpler organism with fewer neurons. A tiny roundworm called Caenorhabditis elegans is one such creature, with only 302 neurons. These worms avoid noxious odors, by first wandering around when they detect the odor, and then fleeing. About 80% of the time the worms flee in the correct direction to escape the foul smell. However, it was not clear how the worms decided which direction to flee.

Now, Tanimoto, Yamazoe-Umemoto et al. show that the worms choose which direction to move by mathematically calculating information about odor concentrations. In the experiments, a robotic microscope simultaneously measured nerve activity and the worm’s behavior while an odor was presented. Specifically, the amount of calcium in the neurons was measured. The experiments showed that when the worms were wandering to determine which direction to flee the amount of calcium in the neurons changed in proportion to how much the concentration of the odor changed overtime.

The experiments suggest that the animals use a mathematical process called integration to add up the changes in the concentration of the odor over time, and when the total reaches a certain threshold the animal successfully moves away from the source. Tanimoto, Yamazoe-Umemoto et al. also identified the gene that enables these calculations. More complicated animals make similar calculations that take into account environmental changes over time when making a decision. Future experiments are needed to determine if more complex animals also use the same mechanism as C. elegans, and whether the same gene is responsible.

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

Introduction

Brains process sensory information to generate various kinds of physiological responses with different timings (i.e., with different latencies to respond): For example, motor control, foraging, decision-making and the sleep-wake cycle range on timescales from milliseconds to days (Buhusi and Meck, 2005; Richelle and Lejeune, 1980). In decision-making, animals choose one from multiple behavioral options based on environmental sensory information, where a temporal delay is associated with the certainty of sensory information. In primates and rodents, increases in neural activity during the delay period according to the sensory information (‘evidence accumulation’) has been described as a key physiological basis for the timing of decision-making (Carandini and Churchland, 2013; Gold and Shadlen, 2007; Schall, 2001; Shadlen and Newsome, 2001). For example, clear sensory information causes a faster rise in neuronal firing rate to a threshold and faster behavioral choice, whereas uncertain information causes a slower rise in the firing rate and slower behavioral choice. Despite their essential roles in the timing of decision-making, the neural mechanisms that generate evidence accumulation still need to be clarified. Theoretical studies suggest that evidence accumulation is mediated by recurrent neural circuits (Gold and Shadlen, 2007; Wang, 2008) while intracellular mechanisms, such as calcium signaling via N-methyl-D-aspartate (NMDA) receptors, calcium-activated nonspecific cation (CAN) channels, and/or voltage-gated calcium channels (VGCCs), have also been proposed (Curtis and Lee, 2010; Major and Tank, 2004).

The cellular and molecular bases of simple decision-making have been studied in invertebrate animals because of the simplicity and accessibility of their nervous systems (Kristan, 2008). For example, a neuron that biases decisions was identified based on optical monitoring of the neuronal activities in the medical leech (Briggman et al., 2005), and genetic analyses identified a neuropeptide and a catecholamine receptor that underpin simple decision-making tasks in Drosophila melanogaster and in Caenorhabditis elegans, respectively (Bendesky et al., 2011; Yang et al., 2008). In addition, a recent behavioral study showed that Drosophila temporally accumulates sensory evidence for decision-making (DasGupta et al., 2014). However, the physiological mechanisms of these decision-making tasks have not been elucidated, and thus not discussed in terms of their possible commonalities with the mechanisms in mammals.

In this study, we reveal the cellular and molecular mechanisms of the timing of decision-making in C. elegans. We first show that the animals migrate in an appropriate direction with unexpectedly high efficiency, likely based on the gradient of a repulsive odorant. From simultaneous monitoring of behavior and neural activity in virtual odor gradients, we find that two pairs of sensory neurons regulate this behavioral response in an opposing manner with different temporal dynamics. A pair of ASH nociceptive neurons exhibits a time-differential-like response to an increase in the odor concentration, which leads to a bout of turns in random directions similar to a ‘reflex’. In contrast, a pair of AWB olfactory neurons exhibits a time-integral response to a decrease in the odor concentration, which leads to turn suppression with a temporal delay resembling ‘deliberation’. The AWB response is mediated by a gradual calcium influx mainly via L-type VGCCs whereas the ASH response is mediated by a rapid calcium influx via multiple types of calcium channels. Thus, our results indicate that the timing of a sensory response, such as deliberate decision-making or rapid reflex, is determined by cell type-dependent involvement of calcium channels.

Results

Efficient directional choice in odor avoidance behavior

C. elegans avoids the odorant 2-nonanone, and this odor avoidance behavior is regulated by stimulus-dependent transitions between the two behavioral states, ‘pirouette’ (a period of short migrations divided by turns) and ‘run’ (a period of long straight migration) (Figure 1A and Figure 1—figure supplement 1A) (Bargmann et al., 1993; Kimura et al., 2010). The pirouette strategy (a form of ‘biased random walk’) is the major behavioral strategy used by these animals for chemotaxis and thermotaxis, and the choice of migratory direction at run initiation is considered to be random in this strategy (Lockery, 2011; Pierce-Shimomura et al., 1999). We found, however, that C. elegans appropriately chose the migratory direction in 2-nonanone avoidance: 78.4% of the migratory directions at run initiation and 83.5% during runs were away from the odor source (Figure 1B and Figure 1—figure supplement 1B). Thus, during odor avoidance, animals chose the appropriate direction about four times more frequently than the inappropriate direction (~80% versus~20%). This probability of run initiation in the appropriate direction for odor avoidance was far higher than that in salt-taxis, the best-studied chemosensory behavior of the animal (59.3% in Figure 1—figure supplement 1C and D, and ~56% in Pierce-Shimomura et al., 1999), and the probability exceeded, or was at least comparable to, that of odor-taxis in Drosophila larvae (73.9%) (Gomez-Marin et al., 2011; Louis et al., 2008). Appropriate directions were chosen from multiple exploratory short migrations in random directions during pirouettes (Figure 1C and D). These results indicate that C. elegans can efficiently choose the appropriate migratory direction on a repulsive odor gradient (Figure 1E).

Figure 1 with 1 supplement see all

Download asset Open asset

*C.elegans* selectively initiates runs away from the odor source.

(A) Examples of the tracks of 2 animals during 12 min of 2-nonanone avoidance assay, overlaid on a schematic drawing of a 9 cm plate. One of the tracks is magnified below. In the magnified view, pirouettes are red and runs are blue. Arrow heads and arrows indicate the directions of run initiation and those during runs, respectively. (B) Histogram indicates the bearings at run initiation during 2-nonanone avoidance (*i.e.*, the bearing of the arrow heads in panel A, and the initial bearing of the blue arrows in panel E). The bearing was determined as B = 0° when animals migrated directly away from the odor source (= down the gradient) and ±180° when they migrated directly toward the source (= up the gradient). Migration away from the odor source (*i.e.,* within ±90° bearings; red arrow) comprised 78.4% of all data. (C) Histogram indicates all the bearings after a turn during pirouettes, including those that later switched to runs (*i.e.*, the initial bearings of the red plus blue arrows in panel E). (D) Percentages of turn numbers per pirouette during 2-nonanone avoidance (solid line and filled circles) and in salt-taxis (dashed lines and open circles). The average number of turns in a pirouette was significantly larger during 2-nonanone avoidance than in salt-taxis (5.2 vs. 2.2 indicated by arrows, respectively; p<0.001 by Mann-Whitney test). (E) Schematic drawings of salt-taxis and 2-nonanone avoidance. Typical chemotaxis such as salt-taxis in *C. elegans* is regulated by the pirouette (*i.e.*, biased random walk) strategy, where the animal initiates runs mostly in a random direction after a pirouette (Pierce-Shimomura et al., 1999). In contrast, in 2-nonanone avoidance, an appropriate migratory direction is chosen from multiple trials in a pirouette. We refer to this as the ‘pivot-and-go’ strategy. All the statistical details are shown in Supplementary file 1. For panels B, C and D (2-nonanone avoidance), the data are from 100 wild-type animals. Panel D (salt-taxis) is from 64 wild-type animals. The following figure supplement is available for Figure 1.

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

C. elegans regulates its behavioral state based on subtle increases and decreases in odor concentration

Next, to understand the correlations between sensory information and appropriate directional choice, we developed a method for measuring and determining the dynamic spatio-temporal pattern of the naturally evaporating and diffusing odor gradient (Figure 2A and Figure 2—figure supplement 1A–D, and Video 1). The measured odor gradient was then used to calculate the odor concentrations that each animal experienced at each position at every second (C_worm in the middle panels of Figure 2B and Figure 2—figure supplement 1E). We found that the temporal changes in C_worm (dC_worm/dt) were strongly correlated with the two behavioral states of the animals (the bottom panels of Figure 2B and Figure 2—figure supplement 1E, and Figure 2C). During pirouettes, the dC_worm/dt values were mostly positive because the animals did not migrate much while the odor concentration on the plate was increasing due to sustained evaporation of the odor from the source (Figure 2—figure supplement 1D and Video 1). During runs, in contrast, dC_worm/dt values were mostly negative because the animals migrated down the gradient. These correlations may suggest that pirouettes and runs are caused by positive and negative dC_worm/dt, respectively. In addition, we further found that the relationships between animals' responsiveness and instantaneous dC_worm/dt largely differed between pirouettes and runs (Figure 2D), suggesting that pirouettes and runs are physiologically distinct behavioral states in terms of sensory response. Taken together, these results suggest that the efficient transitions between discrete behavioral states based on odor concentration information may lead to the appropriate choice of migratory direction as a simple form of decision-making.

Figure 2 with 1 supplement see all

Download asset Open asset

Pirouettes and runs are distinct behavioral states, which are associated with positive and negative *dC_worm/dt*, respectively.

(A) Fitted odor gradient over the assay plate at 12 min, based on the actual measurements shown in Figure 2—figure supplement 1D. (B) (Top) Same with the magnified view of an animal's trajectory in Figure 1A. (Middle, bottom) Graphs showing the 2-nonanone concentration (*C_worm*: middle) or temporal changes in it (*dC_worm/dt*: bottom) at this animal's position at each second during the odor avoidance behavior. As in Figure 1A, pirouettes and runs are red and blue, respectively. Most of the animals did not migrate much during the first 2 min and were excluded from the analysis (green) because of the rapid increases in the odor concentration during this period. See also Figure 2—figure supplement 1E for another example. (C) Correlation between *dC_worm/dt* and pirouettes or runs. *dC_worm/dt* of 2-nonanone was plotted against the animal's x position for each second during pirouettes (top) or runs (bottom). The bars represent the median ± quartiles for each 5 mm fraction. (D) The responsiveness to the instantaneous *dC_worm/dt* differed between pirouettes and runs. The turning rate was determined as the relationship between *dC_worm/dt* during one second of migration and the probability of turning in the next second. Average turning rates ± SEM for every 0.005 μM/s bin during pirouettes (red line) or runs (blue line) are shown. The data in panels C and D are from the same 100 wild-type animals as in Figure 1. The following figure supplement is available for Figure 2.

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

Video 1

Download asset

posterframe for video — Time-course changes in the fitted 2-nonanone concentration.

Although the odor sources were two circles of ~5 mm diameter in the real experiment, they were treated as points in the simulation.

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

Two sensory pathways mediate odor increases or decreases on a virtual odor gradient

To determine whether sensory stimuli are the causal reason for the behavioral response during odor avoidance, and the neural mechanisms linking stimulus and behavior, we developed a novel integrated microscope system for the quantitative optophysiological analyses of freely moving C. elegans on virtual odor gradients. We integrated the OSB system, an auto-tracking microscope system for calcium imaging and optogenetic analyses of neuronal activity in C. elegans (Video 2), with an odor-delivery sub-system (Busch et al., 2012; Tanimoto et al., 2016). Using this new OSB2 system, a moving C. elegans was continuously exposed to an odor flowing from syringe pumps (Figure 3A and Figure 3—figure supplement 1A, and Video 3). The odor concentration in the flow changed according to a predefined program, based on estimated values experienced in the plate assay paradigm (Figure 2). Using the OSB2 system, we first tested whether the dC/dt of the odor itself could regulate an animal's behavior. When animals experienced temporal increases in the odor concentration, they exhibited more frequent turns like in pirouettes (Figure 3B). Conversely, when they experienced temporal decreases, they suppressed turns and exhibited long migrations similar to runs (Figure 3B). During the odor-down phase, no statistical bias was detected in migratory direction (Figure 3—figure supplement 1B). These results indicate that C. elegans responds to temporal changes in the odor concentration and regulates turning rates.

Figure 3 with 1 supplement see all

Download asset Open asset

AWB and ASH sensory neuron pairs regulate turning rate in response to *dC/dt* of 2-nonanone.

(A) Schematic drawing of the OSB2 system. (B) Behavioral response to temporal changes in the 2-nonanone concentration. (Top) Track of a wild-type animal. The first 60 s (gray) is a period of no odor (‘odor-zero phase’ in the bottom left panel), 60–90 s (black) is a period with a constant increase in odor concentration from 0 to 1 μM (‘odor-up phase’), and 90–180 s (gray) is a period with a constant decrease in odor concentration from 1 to 0 μM (‘odor-down phase’). Red dots indicate turns. (Bottom left) Rastergram of turns. The upper portion of the panel shows the measured 2-nonanone concentration in the flow. In the lower portion of the panel, each turn is denoted by a dot, and each row represents the behavioral record of a single animal during the 180 s of analysis. The results of 20 animals are shown. (Bottom right) Ensemble averages ± SEM for the turning rate (turns per second) during each phase in the left panel (n = 20). The turning rates in the three conditions differed significantly from each other (***p<0.001, Kruskal-Wallis test with *post hoc* Steel-Dwass test). (C) Response of AWB neurons. The averages ± SEM of *∆F/F₀* for GCaMP3 (middle) and rastergram of turns of the animals (bottom) are shown (n = 9). (D) Effect of optogenetic activation of AWB neurons on the turning rate in the absence of odor. Transgenic animals expressing ChR2(C128S), a bi-stable variant of ChR2, in AWB were cultivated in the absence or presence of all-*trans*-retinal (ATR) (dashed or solid bars, respectively); exogenous ATR is required for functional ChR in *C. elegans* (Nagel et al., 2005). Average turning rates ± SEM of the 30 s periods before (gray bars) or after blue light illumination for 3 s (blue bars) are shown (n = 20 each). **p<0.01 (Mann-Whitney test). (E) Calcium imaging of ASH neurons using GCaMP3 (n = 7). (F) Effect of optogenetic silencing of ASH neurons on the turning rate. The transgenic animals expressing Arch in ASH neurons were cultivated in the absence or presence of ATR (dashed or solid bars, respectively) and illuminated with green light during the up-phase (n = 16 each). The odor pattern was the same as that in panels B, C, and E. *p<0.05 (Mann-Whitney test). (G) Model of the regulation of odor avoidance by ASH and AWB neurons. All the statistical details are shown in Supplementary file 1. The following figure supplement is available for Figure 3.

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

Video 2

Download asset

Video 3

Download asset

We then performed calcium imaging of neural activity and revealed that two pairs of sensory neurons regulate the avoidance behavior by responding to either increases or decreases in the odor concentration. In C. elegans' sensory neurons, dynamic changes in [Ca²⁺]_i in the cell bodies are generally similar to those in axons, which cause neurotransmitter release from the neurons (Kato et al., 2014; Zahratka et al., 2015). First, we investigated a pair of AWB olfactory neurons. AWB neurons are known to be primarily responsible for 2-nonanone avoidance and to exhibit an odor-OFF response when animals were stimulated with 2-nonanone-saturated buffer in a stepwise manner (Ha et al., 2010; Troemel et al., 1997). However, it was not clear whether and how AWB neurons respond to gradual and/or subtle changes in odor concentration in the air phase. We found that AWB neurons were gradually and continuously activated during the odor-down phase for 90 s (Figures 3C, 90–180 s). Optogenetic activation of AWB neurons by the bi-stable variant of the light-gated cation channel ChR2(C128S) (Berndt et al., 2009) significantly suppressed turns in the absence of the odor (Figure 3D and Video 4). In addition, we also found that ASH neurons, a pair of nociceptive neurons (Bargmann, 2006; Kaplan, 1996), were activated by an odor increase (Figures 3D, 60–90 s). Consistently, optogenetic inactivation of ASH neurons during the odor-up phase by the light-driven H⁺ pump Arch (Chow et al., 2010) significantly suppressed the increase in the turning rate (Figure 3F). Taken together, these results indicate that two distinct sensory pathways respond to the physiological range of odor concentration changes and opposingly regulate the avoidance behavior (Figure 3G).

Video 4

Download asset

Temporal differentiation and integration of odor information regulate avoidance behavior

Through mathematical analyses of ASH and AWB responses, we found, unexpectedly, that the two neuron pairs decode the temporal dynamics of odor concentration information using different computations. When C. elegans was stimulated with different rates of positive dC/dt, the ASH neurons were always activated rapidly (Figure 4A, black lines in middle panels). These responses peaked soon after the onset of the odor-up phase and were mostly maintained during this phase, which can be approximated by time-differentials of the odor concentrations (i.e., dC/dt: Figure 4A, red lines in middle panels; see also Table 1). The time-differential response of sensory neurons is consistent with recent studies (Lockery, 2011; Schulze et al., 2015) and the general idea that sensory neurons respond in a phasic manner by detecting stimulus changes (Delcomyn, 1998). A portion of ASH response (e.g., decay kinetics during the odor-plateau phase) was not well fitted by the time-differential equation, which reflected the fact that ASH response is regulated by multiple mechanisms (see below).

Figure 4 with 3 supplements see all

Download asset Open asset

ASH neurons are activated according to *dC/dt* for initiating turns, and AWB neurons are activated according to the leaky integration of the negative *dC/dt* for suppressing turns with a *dC/dt*-dependent delay.

(A) ASH responses (middle panels) and turns (lower panels) in response to odor concentration increases from 0 to 1 μM in 45 s (left most; n = 32), 90 s (middle left; n = 39), 180 s (middle right; n = 35) or no odor increase (‘no odor control’; right most; n = 41) are shown. (Middle panels) In the three conditions with odor increases, the response patterns of ASH neurons (the average ± SEM: black lines and gray shadows, respectively) were fitted by time-differentials of the average of measured odor concentration (red lines), calculated by the rightmost equations. (Lower panels) Ensemble averages of the turning rates ± SEM in each 10 s bin were calculated. The original data is shown in the raster plot. Black horizontal dashed lines in lower panels indicate the upper limit of 99% prediction interval of all the turning rates during the odor-zero phase (t = −60 ~ 0; gray area). Red vertical dotted lines indicate the time when each turning rate first exceeded the upper limit of prediction interval. In the first bin of odor-up phase (indicated by a black horizontal bar in the lower panels), the turning rate in the 45 s and 90 s conditions increased significantly compared to the no-odor control (***p<0.001, Kruskal-Wallis test with *post hoc* Steel-Dwass test). (B) AWB responses and turns in response to odor concentration decreases from 1 μM (odor-plateau phase; gray area) to 0 μM (odor-zero phase) in 45 s (left-most; n = 42), in 90 s (middle left; n = 43), in 180 s (middle right; n = 48) and no odor decrease (‘odor-plateau control’; right-most; n = 38). AWB responses to different *dC/dt* rates were fitted by a leaky integrator equation of negative *dC/dt* (red lines). Black horizontal dashed lines in the lower panels indicate the lower limit of the 99% prediction interval of the turning rates during the odor-plateau phase. The times when the turning rates became lower than the limit (red vertical dotted lines) were delayed when the *dC/dt* rate was smaller (*p<0.005 and ***p<0.001, Kruskal-Wallis test with *post hoc* Steel-Dwass test). The statistical test was performed in the first 3 bins of the odor-down phase (a black horizontal bar in lower panels). (C), Effect of optogenetic silencing of the AWB response on negative *dC/dt*. Transgenic animals expressing Arch in AWB neurons were cultivated in the absence (n = 18) or presence (n = 19) of ATR and illuminated with green light during the odor-down phase. *p<0.05, **p<0.01, and ***p<0.001 (Mann-Whitney test). In the panels, the gray shading means the period with no odor change, in which the prediction intervals were calculated. All the statistical details are shown in Supplementary file 1. The following figure supplements are available for Figure 4.

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

Table 1

Models and parameters used in the fitting of ASH responses with time-differential equations.

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

Conditions	ASH, wild-type			ASH, odr-3(n2150)	ASH, odr-3(n1605)	ASH, wild-type slow component
Durations of up/down phase	Up 45 s	Up 90 s	Up 180 s	Up 90 s	Up 90 s	Up 90 s
Model	$X (t) = k I$	$X (t) = k I$	$X (t) = k I$	$X (t) = k I$	$X (t) = k I$	$\frac{d X (t)}{d t} = k I - \frac{1}{τ} X (t)$
	$I = \frac{d C (t)}{d t}$	$I = \frac{d C (t)}{d t}$	$I = \frac{d C (t)}{d t}$	$I = \frac{d C (t)}{d t}$	$I = \frac{d C (t)}{d t}$	$I = \frac{d C (t)}{d t}$
Parameters used for fitting to ΔF/F₀
k	56.1 [μM⁻¹·s]	80.9 [μM⁻¹·s]	121.9 [μM⁻¹·s]	49.1 [μM⁻¹·s]	46.7 [μM⁻¹·s]	2.96 [μM⁻¹]
τ	-	-	-	-	-	10.4 [s]
Parameters used for fitting to estimated calcium concentration
k	3.9 [s]	6.2 [s]	7.9 [s]	2.8 [s]	2.9 [s]	-
τ	-	-	-	-	-	-
f_max	9.4	9.5	9.2	9.1	8.7	-
f_min	0.79	0.79	0.77	0.76	0.73	-
X_base	103.8 [nM]	92.2 [nM]	107.0 [nM]	108.2 [nM]	109.7 [nM]	-

The behavioral response also changed rapidly: At the onset of the odor-up phase, the turning rates exceeded the upper limit of the 99% prediction interval for the rate during the odor-zero phase (Figure 4A, red vertical dotted lines in lower panels): The 99% prediction interval is a value in which future data (odor-up phase in this case) will fall with 99% probability based on the observed data (odor-zero phase) (Montgomery and Runger, 2002). In the 45 s and 90 s odor-up conditions (left-most and middle left panels), the turning rates were significantly different from the no odor control (see Supplementary file 1 for statistical details). It should be noted that, even with the slowest odor concentration increase (middle right panel), at which the overall turning rate was around the threshold, the turning rate increased and surpassed the threshold at the onset of the odor-up phase although no statistical difference was detected. These results suggest that the nociceptive ASH neurons compute the time-differential of odor concentration to rapidly cause an aversive response based on a small change in odor concentration.

In contrast, the AWB responses to decreases in odor concentration were time-integral. The responses gradually increased and peaked with a considerable delay, which depended on the dC/dt rates after the onset of the odor-down phase (Figure 4B, black lines in middle panels). These responses could not be fitted by the dC/dt itself but by a leaky integrator equation, in which -dC/dt acted as the input (Figure 4B, red lines in middle panels and Figure 4—figure supplement 1; see also Table 2 for parameters and Table 3 for goodness of fit). In the leaky integrator model, the neuronal response is given by the sum of the past inputs, which decreases exponentially due to leakage. Furthermore, we estimated [Ca²⁺]_i by including a sigmoidal relationship between calcium concentration and fluorescence intensity into the model (Akerboom et al., 2012; Kato et al., 2014). The estimated calcium concentration exhibited a similar pattern to the measured fluorescence change with an estimated basal concentration of about 60–90 nM (Figure 4—figure supplement 2 and Table 2), similar to that generally reported for resting neurons in general (~100 nM) (Clapham, 2007). Taken together, these results suggest that AWB neurons compute the leaky integration of the negative dC/dt to temporally integrate the sensory information as the accumulation of [Ca²⁺]_i.

Table 2

Models and parameters used in the fitting of AWB responses with leaky integrator equations.

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

Conditions	AWB, wild-type			AWB, unc-13(e51)	AWB, unc-31(e928)	AWB, odr-3(n2150)	AWB, odr-3(n1605)
Durations of up/down phase	Down 45 s	Down 90 s	Down 180 s	Down 90 s	Down 90 s	Down 90 s	Down 90 s
Model	$\frac{d X (t)}{d t} = k I - \frac{1}{τ} X (t)$	$\frac{d X (t)}{d t} = k I - \frac{1}{τ} X (t)$	$\frac{d X (t)}{d t} = k I - \frac{1}{τ} X (t)$	$\frac{d X (t)}{d t} = k I - \frac{1}{τ} X (t)$	$\frac{d X (t)}{d t} = k I - \frac{1}{τ} X (t)$	$\frac{d X (t)}{d t} = k I - \frac{1}{τ} X (t)$	$\frac{d X (t)}{d t} = k I - \frac{1}{τ} X (t)$
	$I = - \frac{d C (t)}{d t}$	$I = - \frac{d C (t)}{d t}$	$I = - \frac{d C (t)}{d t}$	$I = - \frac{d C (t)}{d t}$	$I = - \frac{d C (t)}{d t}$	$I = - \frac{C (t) - C (t - Δ t)}{Δ t}$	$I = - \frac{C (t) - C (t - Δ t)}{Δ t}$
						(∆t = 67 s)	(∆t = 66 s)
Parameters used for fitting to ΔF/F₀
k	4.6 [μM⁻¹]	3.6 [μM⁻¹]	3.0 [μM⁻¹]	2.6 [μM⁻¹]	4.9 [μM⁻¹]	8.7 [μM⁻¹]	7.7 [μM⁻¹]
τ	19.7 [s]	28.1 [s]	25.0 [s]	34.5 [s]	28.1 [s]	25.2 [s]	23.5 [s]
Parameters used for fitting to estimated calcium concentration
k	0.40	0.37	0.33	0.26	0.44	1.00	0.87
τ	21.1 [s]	28.5 [s]	25.0 [s]	37.3 [s]	30.7 [s]	17.7 [s]	17.7 [s]
f_max	9.7	10.1	8.5	9.6	8.9	10.9	10.4
f_min	0.81	0.84	0.71	0.80	0.74	0.91	0.87
X_base	66.1 [nM]	56.7 [nM]	89.6 [nM]	66.7 [nM]	81.2 [nM]	41.8 [nM]	51.5 [nM]

Table 3

Parameters and goodness of fit results for mathematical models of ASH and AWB responses.

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

Conditions	ASH, wild-type			AWB, wild-type
Durations of up/down phase	Up 45 s	Up 90 s	Up 180 s	Down 45 s	Down 90 s	Down 180 s
Number of samples (frames) used for calculation of BIC	N = 135 (t = −60 ~ 75 s)	N = 180 (t = −60 ~ 120 s)	N = 270 (t = −60 ~ 210 s)	N = 135 (t = −60 ~ 75 s)	N = 180 (t = −60 ~ 120 s)	N = 270 (t = −60 ~ 210 s)
$X (t) = k I$ $I = \frac{d C (t)}{d t}$	k = 56.1 [μM⁻¹·s] BIC = −222.2	k = 80.9 [μM⁻¹·s] BIC = −446.5	k = 121.9 [μM⁻¹·s] BIC = −637.8	k = −58.4 [μM⁻¹·s] BIC = −170.6	k = −73.4 [μM⁻¹·s] BIC = −372.5	k = −67.3 [μM⁻¹·s] BIC = −1157
$\frac{d X (t)}{d t} = k I - \frac{1}{τ} X (t)$ $I = \frac{d C (t)}{d t}$	k = 5.8 [μM⁻¹] τ = 12.0 [s] BIC = −331.1	k = 18.9 [μM⁻¹] τ = 4.4 [s] BIC = −472.8	k = 11.5 [μM⁻¹] τ = 11.7 [s] BIC = −804.9	k = −4.56 [μM⁻¹] τ = 19.7 [s] BIC = −591.9	k = −3.58 [μM⁻¹] τ = 28.1 [s] BIC = −806.2	k = −3.01 [μM⁻¹] τ = 25.0 [s] BIC = −1458

The dC/dt rate-dependent delay in AWB activation was correlated with the behavioral responses. The time when the turning rate became lower than the 99% prediction interval of the odor-plateau phase was delayed according to the dC/dt rate, which was associated with statistical differences from the odor-plateau control (Figure 4B, red vertical dotted lines in lower panels). This result suggests that the turning rate is suppressed when the AWB activity exceeds a certain value (Figure 4B, red horizontal arrows in middle panels). Interestingly, even when animals were stimulated with the most shallow odor gradient in which the average turning rate appeared gradually decreased (middle right panels in Figure 4B), turn intervals of each animal could be classified into two groups like pirouettes and runs (compare Figure 1—figure supplement 1E with 1A). This result suggests that, even with subtle changes in odor concentration, the behavioral response of C. elegans does not change gradually, but instead transits between high- and low-turning states. Optogenetic suppression of the AWB neurons with Arch during the odor-down phase significantly affected the transition (Figure 4C). Taken together, these results indicate that AWB neurons regulate the transition from a pirouette to a run based on temporal integration of the negative dC/dt.

Computer simulations reveal the behavioral effects of the temporal integration of sensory information

To further understand whether the neural computations of odor concentration contribute to the choice of migratory direction during decision-making, we performed computer simulations of the odor avoidance behavior (Figure 5A). During pirouettes, the model animal frequently repeated turns and short migrations in random directions, and a pirouette transited to a run when the time-integral of -dC_worm/dt over time reached a threshold value. The model produced similar migration patterns to the odor avoidance behavior seen in real animals, in which most of the runs were down the odor gradient (Figure 5B). Interestingly, when we changed the model animal's computation for transition of a pirouette to a run from the temporal integration of dC/dt (the integration model) to the simple dC/dt itself (the differentiation model), the directional choice at run initiation significantly worsened (Figure 5C; p<0.001, Mardia-Watson-Wheeler test). This is because, in the differentiation model, the animals initiated runs even when they transiently sensed dC/dt <0, due to the periodical head swing for example (Iino and Yoshida, 2009; Yamazoe-Umemoto et al., 2015). In contrast, animals in the leaky integration model initiated runs only when they sensed dC/dt <0 for a certain period of time, resulting in an appropriate directional choice. In reality, the movement of the animal's anterior end (where the sensory endings of the ASH and AWB neurons are exposed to the environment) is more random than that in our model (Videos 2 and 4), suggesting that the temporally integrating property of AWB neurons may play a significant role in robust directional choice based on a noisy sensory input.

Figure 5

Download asset Open asset

A computer model reproduced the directional choice in the odor avoidance task in a temporal integration-dependent manner.

(A) Model of the behavioral transition in 2-nonanone avoidance. During a pirouette, a model animal frequently repeated turns and short migrations. When a model animal initiated a migration away from the odor source and sensed *dC/dt* <0, the high-turning state transited to a low-turning state (*i.e.*, a run) when the leaky integrator equation exceeded a threshold. During a run, the animals turned with much lower frequency than in pirouette and at a probability related to *dC/dt* >0 due to straying away from the original direction. (B) Three typical tracks of real (left) and model (right) animals. As shown in Figure 1A, the odor source is on the left side. (C) Histograms of the initial bearings of runs in the model animals. The high-turning-to-low-turning transition was dependent on temporal integration (left) or on differentiation (right). p<0.001 (Mardia-Watson-Wheeler test). All the statistical details are shown in Supplementary file 1.

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

Temporal integration of sensory information occurs in AWB neurons

Temporally integrating neural activity for decision-making and working memory is considered to be based on recurrent synaptic circuits in vertebrates (Aksay et al., 2007; Wang, 2008). To determine whether a synaptic circuit input is required for the temporally integrating property of AWB neurons, we abolished the transmission of synaptic and/or dense core vesicles in most, if not all, neurons via a mutation in unc-13 or unc-31, the orthologs of mammalian Munc-13 and CAPS (calcium activated protein for secretion), respectively (Richmond, 2005). In the unc-13 and unc-31 mutants, however, the AWB responses were essentially similar to that in wild-type animals and fitted with the same leaky integrator equations (Figure 6A and Table 2), indicating that the temporally integrating neural activity occurs within the AWB neuron itself and does not require a recurrent circuit.

Figure 6 with 1 supplement see all

Download asset Open asset

Cell-autonomous computations in AWB neurons.

(A) The AWB responses of *unc-13* (left) and *unc-31* (right) mutants to the odor decreases, which are the same as those shown in the middle left panel of Figure 4B, were essentially similar to those of wild-type animals and fitted by the leaky integrator equations (n = 31 and 17, respectively). (B) The AWB responses of *odr-3(n2150)* mutants (n = 24) were fitted by the leaky integrator equation (red solid line) with the time interval (∆t) being much longer than that of wild-type AWB (red dotted line). The parameters are described in Table 2. In the experiments for panels A and B, the behavioral responses were not analyzed because the *unc-13*, *unc-31*, and *odr-3* mutations affect the activities of multiple neurons, including ASH and AWB.

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

Molecular mechanisms of neural computation in the sensory neurons

To reveal the molecular mechanisms of the temporal integration of sensory information in AWB neurons, we performed genetic manipulations. First, we identified the gene required for the calculation of dC/dt. Previous genetic studies suggest that the signaling pathway mediating odor response in AWB neurons is similar to that in mammalian olfactory/photoreceptor neurons. Upon binding to an odorant, a G-protein coupled receptor activates Gα protein, which eventually leads to the opening of cyclic nucleotide-gated cation channels for depolarization (Bargmann, 2006). ODR-3, one of the 20 Gα homologs in C. elegans, is expressed in AWB, ASH and a few other pairs of sensory neurons. It is localized at the sensory ending, and its dysfunction causes severe chemotaxis defects, suggesting that ODR-3 transduces the sensory signal from odorant receptor (Roayaie et al., 1998). Two alleles of loss-of-function mutations of the odr-3 gene, however, caused gradual activations of AWB neurons, whose peaks were even larger than those of wild-type animals (Figure 6B and Figure 6—figure supplement 1A, black lines), suggesting that ODR-3 may play an inhibitory role in sensory signaling. Interestingly, while the AWB response patterns in wild-type animals were fitted by a leaky integrator equation with a time interval for the input (∆t)=1 s (solid red lines in Figure 4B and dotted red line in 6B), AWB response patterns in odr-3 mutants were better fitted by an equation with ∆t longer than 1 min than by ∆t = 1 s (solid red line in Figure 6B, Figure 6—figure supplement 1A; see also Table 2 for parameters and Table 4 for goodness of fit). Thus, the time-differential property for the sensory input was greatly degraded with the odr-3 mutations while the time-integral property was not affected. This result suggests that the time-differential computation of odor concentration depends mostly on the ODR-3 Gα protein (see Discussion for details). The result also supports the idea that AWB neurons possess a temporally integrating property for sensory inputs. The overall response of ASH neurons, which also express ODR-3, was partially affected by the odr-3 mutations (Figure 6—figure supplement 1B), suggesting that ODR-3 may play a primary role in sensory signaling in ASH neurons.

Table 4

Parameters and goodness of fit for mathematical models of AWB responses in odr-3 mutants.

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

Conditions	AWB, odr-3(n2150)	AWB, odr-3(n1605)
Durations of up/down phase	down 90 s	down 90 s
Number of samples (frames) used for calculation of BIC	N = 180 (t = −60 ~ 120 s)	N = 180 (t = −60 ~ 120 s)
$X (t) = k I$ $I = - \frac{d C (t)}{d t}$	k = 89.1 [μM⁻¹·s] BIC = 25.1	k = 77.3 [μM⁻¹·s] BIC = −39.2
$\frac{d X (t)}{d t} = k I - \frac{1}{τ} X (t)$ $I = - \frac{d C (t)}{d t}$	k = 2.08 [μM⁻¹] τ → ∞ BIC = −477.7	k = 1.77 [μM⁻¹] τ → ∞ BIC = −578.5
$\frac{d X (t)}{d t} = k I - \frac{1}{τ} X (t)$ $I = - \frac{C (t) - C (t - Δ t)}{Δ t}$	k = 8.71 [μM⁻¹] τ = 25.2 [s] ∆t = 67 [s] BIC = −645.6	k = 7.73 [μM⁻¹] τ = 23.5 [s] ∆t = 66 [s] BIC = −779.9

Furthermore, we genetically identified the genes responsible for the time-integral calcium accumulation. In general, [Ca²⁺]_i increases in neurons depend on its influx through the plasma membrane via VGCCs, which can then trigger rapid calcium release from the endoplasmic reticulum (ER) via IP₃ receptors (IP₃Rs) and/or ryanodine receptors (RyRs) as calcium-induced calcium release (CICR) (Catterall, 2011; Clapham, 2007). The pore-forming α₁ subunits of VGCCs are classified into L-type, N/P/Q-type, and T-type subgroups. Previous studies have demonstrated the requirement for EGL-19 and/or UNC-2, the sole orthologs of L- and N-type channels, respectively, in the responses of sensory neurons of C. elegans (Busch et al., 2012; Frøkjaer-Jensen et al., 2006; Hilliard et al., 2005; Kato et al., 2014; Larsch et al., 2015; Suzuki et al., 2003; Zahratka et al., 2015) although their precise roles in intracellular calcium dynamics have remained unclear. We found that an egl-19 reduction-of-function mutation (rf) (Lee et al., 1997) significantly affected the AWB response (Figure 7A). Consistently, treatment with the EGL-19 antagonist Nemadipine-A (NemA) (Kwok et al., 2006) exhibited a similar, or possibly even stronger, effect. In contrast, mutations in unc-2, itr-1 (the IP₃R ortholog) or unc-68 (the RyR ortholog) (Dal Santo et al., 1999; Sakube et al., 1997; Schafer and Kenyon, 1995) did not significantly affect the response (Figure 7A and C left panel). These results suggest that the calcium accumulation in AWB neurons mostly depends on influx through the EGL-19 L-type VGCCs, but not other calcium channels.

Figure 7 with 1 supplement see all

Download asset Open asset

Calcium channels are involved in the dynamic regulation of [Ca²⁺]_i in a cell type-dependent manner.

(A and B) Responses of AWB (panel A) or ASH (panel B) neurons in strains with genetic and/or pharmacological suppression of N/P/Q-type VGCC UNC-2, T-type VGCC CCA-1, L-type VGCC EGL-19, IP₃R ITR-1, and RyR UNC-68. lf is a loss-of-function mutation, and rf is a reduction-of-function mutation. Note that loss-of-function mutations of *egl-19* and *itr-1* are not available due to possible lethality. (C) Suppression of the AWB response during the odor-down phase (left) and of the initial ASH response (10–15 s of the odor-up phase) (right) with the genetic mutation (M) and/or the drug treatment (D). ***p<0.001 (Kruskal-Wallis test with *post hoc* Steel-Dwass test). (D) *egl-19* is also responsible for the slow time-integral component in ASH neurons. Addition of a putative time-integral response using the leaky integrator equation to the transient response in the ‘*egl-19(drug)*’ reproduced the ASH response of the naive wild-type animals. Black line: wild-type response shown in panel B; black dashed line: *egl-19(drug)* response in panel B; red dashed line: the time-integral model of positive *dC/dt*; red line: sum of the black and red dashed lines. The parameters are described in Table 1. All the statistical details are shown in Supplementary file 1.

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

In ASH neurons, whose response pattern was also not affected by synaptic connections (Figure 7—figure supplement 1), the suppression of EGL-19 with NemA mainly affected the later phase of the response, but did not significantly affect the initial response (Figure 7B and C right panel). The mathematical addition of the AWB-like time-integral response to the NemA-treated ASH response resembled the wild-type ASH response (Figure 7D), suggesting that the ASH response consists of both fast and slow components. Because increases in the turning rate occurred at the onset of the odor-up phase (Figure 4A), the fast component may have a major influence on the behavioral response. This fast component was not affected by a loss-of-function mutation in unc-2 or in a T-type VGCC homolog, cca-1 (Steger et al., 2005), suggesting that this calcium response is not mediated by the typical VGCCs. Suppressing CICR by mutations of itr-1 or unc-68, or using dantorolene, a RyR antagonist (Krause et al., 2004), partially affected the magnitude, but not the temporal pattern, of the ASH response regardless of NemA treatment (Figure 7B and C right panel). Taken together, these results suggest that the ASH response may be mediated by rapid calcium influx via an unidentified type of calcium channels and by slow influx via EGL-19 L-type VGCCs, both of which are amplified by CICR. In previous studies, ASH response was significantly suppressed by NemA-treatment as well as by the egl-19(n582) mutation upon stimulation with another repulsive odorant (1-octanol) (Zahratka et al., 2015) although ASH response was essentially unaffected by the same egl-19 mutation upon stimulation with glycerol (Pokala et al., 2014), suggesting that different sensory stimuli may be processed differently in the polymodal ASH neurons.

Discussion

In this study, we quantitatively analyzed olfactory behavior in C. elegans and found that the animals appear to choose the appropriate migratory direction efficiently on the gradient of a repulsive odor based on subtle changes in odor concentration. This result suggests a novel behavioral strategy for navigating sensory gradients (Figure 1E). One of the problems in analyzing animals' responses in a traditional behavioral arena is that the behavior of animals occurs in a closed loop configuration, i.e., the change in movement directly feeds back onto sensory input. Therefore, when the changes in sensory stimulus and behavior are simultaneously observed, it is difficult to distinguish whether sensory input or behavior is the cause of the change. To solve this problem, we established a rigorous quantitative platform to understand the sensory behavior of C. elegans: We measured the temporal changes in sensory stimulus in a traditional behavioral paradigm and reproduced it in an integrated microscope system (the OSB2 system) in an open loop configuration. This made it possible to quantitatively measure behavioral responses caused by sensory stimuli, and the neural activities mediating them.

Based on the results obtained in our system, we consider that the olfactory behavior of the animals comprises two essential features of decision-making at the behavioral level. (1) Discrete behavioral choice based on sensory information: Based on dC/dt, C. elegans choose either a high-turning state (pirouette) or low-turning state (run), which were shown to be physiologically distinct (Figure 2D and Figure 1—figure supplement 1A and E). (2) Timing of decision: The time required for the transition from a high-turning state to a low-turning state is shorter when the odor concentration decreases rapidly, and is longer when the concentration decreases slowly (Figure 4B). Moreover, we found that [Ca²⁺]_i is accumulated in AWB neurons according to a leaky integrator equation of odor concentration change, which is correlated with the behavioral transition (Figure 4B and Figure 4—figure supplement 2). Because odor concentration is a critical piece of information for the choice of behavioral states, the temporal integration in AWB neurons can be regarded as the ‘evidence accumulation’ for decision-making.

Transition between two behavioral states is regulated by different computations of odor concentration

The behavioral transitions in 2-nonanone avoidance are regulated by ASH neurons for odor increase and AWB neurons for odor decrease (Figure 8A). ASH neurons respond to unfavorable changes in the odor concentration in a time-differential manner and initiate turns rapidly. Such a ‘reflex’-like response is consistent with the nociceptive features of ASH neurons and their synaptic connectivity, whereby the animals respond to various noxious stimuli by inducing turns and reversals (Bargmann, 2006; Kaplan, 1996). Time-differential properties have been suggested for the transient responses of ASH and other sensory neurons with all-or-none stimuli in C. elegans (Lockery, 2011). This property has also been reported in a recent study of Drosophila olfactory neurons (Schulze et al., 2015) and is in agreement with the general idea that dynamic signal changes are mediated by phasic sensory receptors (Delcomyn, 1998). Mathematically, the ASH activity can also be approximated by the same leaky integrator equation as AWB, whose fitness is higher than the time-differential equation according to the Bayesian Information Criterion (BIC) (Figure 4—figure supplement 3 and Table 3). This is consistent with the recent report that ASH neurons temporally integrate sensory signals over several seconds (Kato et al., 2014). However, the onset of responses was not sufficiently reproduced with the model (Figure 4—figure supplement 3, red arrows). Our result in Figure 7D suggests that the ASH response consists of the time-integral component mediated by EGL-19 L-type VGCCs and the fast and transient component, possibly mediated by as-yet-unidentified channels. Nociceptive ASH neurons may have developed specialized mechanisms to quickly cause aversive responses based on a slight change in the undesirable signal.

Figure 8

Download asset Open asset

Physiological and molecular models of decision-making by *C. elegans* during odor avoidance.

(A) Computations of ASH and AWB neurons during odor avoidance behavior. (B) Model of the molecular mechanisms for temporal computation of odor information in AWB and ASH neurons. (Left) In AWB neurons, odor decreases likely cause the activation of Gα proteins as an odor-OFF response (Bargmann, 2006; Usuyama et al., 2012), where an unidentified Gα positively transduces the signal and ODR-3 inhibits the signal for the time-differential computation. The Gα signaling is transmitted to the cGMP-gated cation channel TAX-2/TAX-4 (Bargmann, 2006) to cause depolarization. Depolarization then triggers calcium influx via EGL-19 at the cell body, which causes the gradual accumulation of [Ca²⁺]_i. (Right) In ASH neurons, the depolarization at the sensory ending triggers an unidentified rapid and transient calcium channels, as well as EGL-19. The calcium influx through these channels is amplified by CICR via RyR (UNC-68) and IP₃R (ITR-1).

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

In contrast, AWB neurons extract long-term information about favorable changes in odor concentration using the time-integral property, which leads to ‘deliberate’ transitions from pirouettes to runs. This time-integral property has not been characterized previously in the sensory response of C. elegans. The sensory responses of the animals to gradual signal changes have been reported in a few cases, although they were time-differential, stochastic or tonic, rather than time-integral (Biron et al., 2008; Kimura et al., 2004; Luo et al., 2014). In reality, sensory signals likely change with noise rather than at a constant rate. Moreover, an animal's movement itself causes fluctuations in the sensory input. Thus, it is reasonable to conclude that animals extract sensory information over a longer time window by using a time-integral property, such as that described here for AWB neurons. To trigger behavioral transitions according to neuronal activity, a threshold for temporal integration is required. AIZ interneurons, the major postsynaptic targets of AWB neurons (White et al., 1986), have been shown to elicit digital-like excitability and trigger turns upon activation (Li et al., 2014), suggesting that AWB neurons inactivate AIZ neurons in an all-or-none manner at a certain threshold.

Molecular mechanisms of the temporal integration of sensory information

Through genetic analyses, we investigated the molecular mechanisms underlying the computations in AWB and ASH neurons. First, we found that the temporal integration of sensory information occurs within the AWB neuron itself and is not dependent on synaptic connections. To the best of our knowledge, this result is the first experimental demonstration of evidence accumulation for decision-making within a neuron (see below). We also found that the ODR-3 Gα protein possibly plays a time-differential role in the sensory signaling of AWB neurons. The result was surprising because ODR-3 has been considered to play a major stimulatory role in sensory signal transduction, not in its inhibitory modification (Bargmann, 2006; Roayaie et al., 1998). Because ODR-3 is homologous to the inhibitory Gαi proteins (Bargmann, 2006), it is possible that the sensory signal via another unidentified Gα might be downregulated by the activity of ODR-3—such an inhibitory feed-forward loop could compute the time-differential of the input (Figure 8B, left) (Alon, 2007). In ASH neurons, ODR-3 appears to transduce the sensory signal itself (likely with other Gα); the different roles of ODR-3 in AWB and ASH neurons may be related to the fact that AWB and ASH neurons use different molecular pathways for signal sensation (Bargmann, 2006).

Together with the finding that the time-integral calcium increase in AWB neurons is mainly mediated by EGL-19 L-type VGCCs, these results suggest the following model (Figure 8B, left). A constant change in odor concentration leads to a constant and persistent depolarization with a time-differential activity of ODR-3 Gα at the sensory ending. This depolarization is conducted to the cell body and causes the constant activation of EGL-19 L-type VGCCs and the constant influx of extracellular calcium, which leads to the gradual accumulation of [Ca²⁺]_i. In agreement with this model, L-type, but not N/P/Q-type or T-type, VGCCs are known to be activated continuously during depolarization (Hille, 2001). This gradual [Ca²⁺]_i accumulation causes a temporal delay and deliberate decision, acting as a low-pass filter to smooth out noisy sensory inputs. By contrast, in ASH neurons, the depolarization causes the activation of EGL-19 L-type VGCCs and other type(s) of rapid and transient calcium channels, where a small influx of calcium leads to its amplification with CICR and to a rapid behavioral response (Figure 8B, right).

In principle, temporal integration of sensory information for behavioral choice in C. elegans during odor avoidance is analogous to evidence accumulation in perceptual decision-making in mammals. In both cases, animals extract a long-term trend in sensory information by temporally integrating the information for decision-making (Gold and Shadlen, 2007; Schall, 2001). While AWBs are sensory neurons, their integration property may reflect a central nervous system-like function; C. elegans' sensory neurons are known to possess higher-order functions, such as learning and memory, because of the quite small number of neurons in their nervous system (de Bono and Maricq, 2005; Sasakura and Mori, 2013). In higher animals, temporally integrating neural activity is required for decision-making as well as for other brain functions, such as working memory, and is generally considered to be mediated by recurrent neural circuits (Gold and Shadlen, 2007; Schall, 2001; Wang, 2008). This temporally integrating activity has also been suggested to be mediated by intracellular calcium signaling within a cell (Curtis and Lee, 2010; Loewenstein and Sompolinsky, 2003; Major and Tank, 2004); however, this had not been demonstrated experimentally. We speculate that single-cell temporal integrators with L-type VGCCs, such as the AWB neurons, may also be involved in decision-making and working memory in higher animals. Interestingly, the odor avoidance behavior of C. elegans is regulated by dopamine (Kimura et al., 2010), a neuromodulator involved in mammalian decision-making (Schultz, 2007). Thus, further genetic analysis of the behavioral choice during odor avoidance behavior in C. elegans may identify previously undescribed evolutionarily conserved molecular mechanisms that are responsible for decision-making.

Materials and methods

Strains

The techniques used for culturing and handling C. elegans were essentially as described previously (Brenner, 1974). The C. elegans wild-type Bristol strain N2, RRID:WB-STRAIN:JD21 cca-1(ad1650), RRID:WB-STRAIN:MT1212 egl-19(n582), RRID:WB-STRAIN:JT73 itr-1(sa73), RRID:WB-STRAIN:CX3222 odr-3(n1605), RRID:WB-STRAIN:CX2205 odr-3(n2150), RRID:WB-STRAIN:CB55 unc-2(e55), RRID:WB-STRAIN:MT7929 unc-13(e51), RRID:WB-STRAIN:DA509 unc-31(e928), and RRID:WB-STRAIN:CB540 unc-68(e540) were obtained from the Caenorhabditis Genetics Center (University of Minnesota, USA). In all the behavioral and physiological experiments, young adult hermaphrodites were used.

Multi-worm tracking and analysis of 2-nonanone avoidance on a 9 cm plate

Request a detailed protocol

Quantitative analysis of the 2-nonanone avoidance of wild-type animals in the 9 cm agar plate was carried out as previously described (Kimura et al., 2010; Yamazoe-Umemoto et al., 2015). In brief, several adult animals per assay were transferred to the center of a 9 cm NGM agar plate either directly from a standard 6 cm nematode growth medium (NGM) plate with the food bacteria OP-50 (‘fed’) or after a 1 hr starvation on the NGM plate without OP-50 (‘starved’). In the following analysis, we used a data set that is a mixture of 50 fed and 50 starved animals. Although we did not find a significant difference between fed and starved animals in 2-nonanone avoidance (Kimura et al., 2010), the feeding state (fed or starved) could affect some aspects of C. elegans' behavior (Bargmann, 2006) and we wanted to focus on feeding state-independent behavioral aspects of the animals. Two μL of 30% 2-nonanone (diluted in EtOH) were put in two spots on the surface of the agar plate (Figures 1A and 2A), and images of the animals during the avoidance behavior were captured at 1 Hz for 12 min by our multi-worm tracking system with a high-resolution camera in a fixed position (Kawazoe et al., 2013; Yamazoe-Umemoto et al., 2015). In this study, we used a CMOS camera CSB4000F-10 (Toshiba Teli Corp., Japan) equipped with a C mount adaptor and a Nikkor 50 mm f/1.2 lens (Nikon Corp., Japan). Because the camera captures the entire area of the 9 cm plate with a resolution of 2008 × 2044 pixels, an animal of length ~1 mm and width ~0.05 mm is depicted in ~25 pixels. x−y coordinates of the centroids of the animals in each image were measured by Move-tr/2D software (Library Inc., Japan), and were further analysed by Excel2010 (Microsoft) or R (The R Project). Because the animals did not initiate avoidance during the first 2 min on average (Kimura et al., 2010) (Figures 1A and 2B and Figure 2—figure supplement 1E), data between 121–720 s were used for the analysis (Tanimoto et al., 2017).

Definition of pirouettes and runs

Request a detailed protocol

A pirouette is a period of frequent turns and migrations whose duration is shorter than a threshold value (Pierce-Shimomura et al., 1999). The animal's behavioral state in one second was classified as a turn if the absolute value of angle change in migratory vector of the animal's centroid from the previous second (i.e., during 1 s) was larger than 90° or if the migratory velocity was smaller than 0.1 mm/s in the following frames after the large angle change. According to this definition, the reverse and the omega turn (Gray et al., 2005) were recognized as turns. A distribution of turn intervals (i.e., migratory durations) during 2-nonanone avoidance was well-fitted by a sum of two exponentials for shorter and longer intervals (Figure 1—figure supplement 1A). A period at which the numbers of the short and long intervals were equal was 13.1 s and determined as t_crit according to the original definition (Pierce-Shimomura et al., 1999). Migrations whose turn interval was longer than t_crit were classified as runs, and migrations shorter than t_crit as well as turns were classified in pirouettes.

Bearings at the initiation of and during runs

Request a detailed protocol

The directions of animal migrations for 1 s were defined in terms of the bearing, B, with respect to the 2-nonanone gradient, where B = 0° indicates migration directly away from the odor source (i.e., down the gradient) and B = ±180° indicates migration directly toward the odor source (i.e., up the gradient). Bearing at run initiation in salt-taxis by the previous study was calculated from the results of taxis toward NH₄Cl (56.0%) and biotin (55.2%) in Figure 9 of the report (Pierce-Shimomura et al., 1999).

Calibration curve for 2-nonanone measurement

View detailed protocol

Calibration curve for 2-nonanone measurement is described in more detail at Bio-protocol (Yamazoe-Umemoto et al., 2018). To measure local concentrations of gaseous 2-nonanone in the assay plate, we used a gas chromatograph (GC) with a sensitive semiconductor detector, SGVA-N2, which was optimized for 2-nonanone detection (FIS Inc., Japan). To make a calibration curve for the measurement, 0.36, 1.07, 3.56, 35.6, 59.4, 97.2, and 200 μL of liquid 2-nonanone (Wako Pure Chemical, Japan) were vaporized in a 50 L tank DT-T1 (FIS Inc.), each corresponding to 0.04, 0.12, 0.4, 4.0, 6.8, 11.1, and 22.9 μM in the gas phase, respectively. After the volatilization period, 0.2 mL of the gas was sampled with a 2 mL plastic syringe with a needle from an outlet of the tank and was immediately injected into the GC. The volatilization periods were determined for each amount of the liquid to maximize the 2-nonanone signal. Synthetic air Alphagaz 1 (Air Liquide, Japan) was used as a carrier gas. With 260 s retention time, a single large peak of signal intensity (mV) was detected as 2-nonanone signal (Figure 2—figure supplement 1B). The experiments were repeated 3–4 times for each concentration. The correlations between the peak height of the signal and the gaseous 2-nonanone concentration in a log-log plot were well-fitted by two simple regression lines for lower and higher concentrations (Figure 2—figure supplement 1C; R² >0.999 for both). In general, for semiconductor detectors, the correlation between the peak height of the signal and signal concentration in a log-log plot are well-fitted by two simple regression lines for lower and higher concentrations.

Measuring odor gradient by gas chromatograph

View detailed protocol

Measuring odor gradient by gas chromatograph is described in more detail at Bio-protocol (Yamazoe-Umemoto et al., 2018). For the odor sampling, a hole of 1 mm in diameter was made through the bottom of the plastic plate and the agar. Because the molecular weight of 2-nonanone (FW 142.2) is larger as a volatile compound, it did not leak easily from such a small hole. 1, 3, 6, 9 and 12 min after placing the odor at the two spots, a 2 mL plastic syringe, which is the same type as the one used in the calibration, was inserted into the plate through the hole from the bottom, and 0.2 mL of the gas phase was sampled (Figure 2—figure supplement 1A). Each plate was used only once to avoid disturbance of the gradient by the sampling. The sampled gas was immediately injected into the GC for measurement. The concentration of 2-nonanone was calculated from the height of the signal peak according to the regression line for the calibration. For each data point, the measurements were repeated 7–9 times and median and quartile was calculated for the fitting.

Fitting the odor gradient and calculation of C_worm

View detailed protocol

Fitting the odor gradient and calculation of C_worm are described in more detail at Bio-protocol (Yamazoe-Umemoto et al., 2018). The least squares method was used to fit the measured concentration. In the closed plate, the odor concentration asymptotically approaches a constant value. Therefore the measured concentrations were fitted to a phenomenological curve with two exponential saturation functions: C(x, y, t) = a(r₁)(1-exp(-b(r₁)t)) + a(r₂)(1-exp(-b(r₂)t)). r₁ and r₂ are the distances from the position (x, y) on the agar to the two odor sources. The asymptotic concentration a(r) and the increasing rate b(r) are functions of the distance r such as a(r)=a₀ exp(-a₁r - a₂r²) and b(r)=b₀ exp(-b₁r - b₂r²). The assumption that C(x, y, t) is given by the sum of the two independent functions is valid for the low concentration regions x > 0. The fitting parameters a₀ = 20.68 μM, a₁ = 0.7355 cm⁻¹, a₂ = −0.05408 cm⁻², b₀ = 0.8384 min⁻¹, b₁ = 0.7835 cm⁻¹ and b₂ = −0.05761 cm⁻² were determined by the Levenberg-Marquardt method (Press et al., 1992). We consider the measured and the fitted odor gradient as reliable because it is consistent with the fact that the amount of 2-nonanone at the source was apparently reduced by 20–30% after 12 min and with a theoretically calculated simulation (Yamazoe-Umemoto et al., 2015). The 2-nonanone concentration at a given temporal and spatial point of an animal’s centroid was calculated from the fitting curve and was designated as C_worm. Turning rate shown in Figure 2D was determined as the relationship between the dC_worm/dt during one second of migration and the probability of turning in the next second.

Molecular biology and germline transformation

Request a detailed protocol

For the cell-specific expression of mCherry (Shaner et al., 2004), GCaMP3 (Tian et al., 2009), ChR2(C128S) (Berndt et al., 2009) and Arch (Chow et al., 2010), str-1 (Troemel et al., 1997) or srd-23 promoter (Colosimo et al., 2004) was used for AWB-expression, and sra-6 promoter (Troemel et al., 1995) was used for ASH-expression. Germline transformation was performed using microinjection (Mello et al., 1991). The plasmids and strains used in this study are listed in Supplementary file 2 and 3. In Figures 3, 4 and 6B and Figure 6—figure supplement 1, multiple transgenic lines were used for each type of experiment, and the different lines produced similar results. The representative transgenes (i.e., extra chromosomal arrays) were used for genetic analyses in Figures 6A and 7.

Behavioral tracking with the integrated microscope system

Request a detailed protocol

For the OOSaCaBeN (Olfactory and Optogenetic Stimulation associated with Calcium imaging on Behaving Nematode, or OSB2) system, we integrated an auto-tracking microscope system for calcium imaging and optogenetic manipulation with an odor-delivery subsystem (Busch et al., 2012; Tanimoto et al., 2016). Briefly, a wild-type C. elegans (N2) on a NGM plate was placed on a motorized stage HV-STU02 (HawkVision, Japan) combined with an upright microscope and illuminated with infrared light. Bright field images of the animal were acquired by a charge-coupled device (CCD) camera at 200 Hz to regulate the motorized stage for maintaining the region-of-interest (ROI) of a freely moving animal in the center of the view field of the microscope (‘ROI-tracking’, Video 2) (Maru et al., 2010). A ROI was set around the head neuropil. The system also allowed us to maintain the centroid of a whole animal in the center of the view field (‘centroid-tracking’, Video 4). ROI-tracking was used for calcium imaging with a 20× objective lens, and centroid-tracking was used for wild-type behavioral analyses and optogenetic behavioral analyses with a 10× objective lens. The behaving animal was continuously exposed to an odor flowing from two syringe pumps (Video 3), which changed the odor concentration according to a predefined program.

Odor delivery

View detailed protocol

For odor delivery, 4 μM of 2-nonanone was sampled from the 50 L vaporizing tank with a 25 ml Gastight Syringe (Hamilton, USA). Two such syringes were set on a syringe pump HV-SSP01 (HawkVision, Japan) that was controlled by the same program for the auto-tracking. Adapting the gas delivery strategy described previously (Busch et al., 2012), one pump was used for 2-nonanone and the other one was for air. The pump speeds were programmed to deliver a constant gas flow of 8 mL/min from the end of the tube, but with varying combinations from each pump to make the temporal gradient of 2-nonanone concentration. For example, when the pump speed of 2-nonanone syringes was changed from 2 ml/min to 5 ml/min, the air syringes went from 6 ml/min to 3 ml/min during the same period. The programs of the pump speeds were designed so that the magnitude of dC/dt was similar to that which animals experienced during the odor avoidance assay in the plate (Figure 2). The actual concentration of 2-nonanone was monitored at the end of the tube by the same type of semiconductor sensor as the one in the gas chromatograph (GC), and the values were recorded with a PC via a digital multimeter MAS345 (Mastech, Hong Kong) before and after the behavioral assays for each day. The sensor was calibrated every day with a similar method as the GC, with calibration concentrations of 0.5, 1, 2, and 4 μM.

We consider that the measured odor concentrations (Figures 3, 4, 6 and 7) closely matched to the actual odor concentration that the animals experienced during the odor avoidance behavior (Figures 1 and 2) because of the following reasons. (1) The tube end was always maintained at ~1 mm from the freely-moving animal during tracking, and the entire body of an animal was exposed to essentially uniform odor flow without significant diffusion and/or turbulence (Figure 3—figure supplement 1A and Video 3). (2) With the flow (8 ml/min), the animals exhibit robust behavioral response reproducibly through multiple trials (Figures 3 and 4) while the flow itself did not affect the animal's behavior. (3) The animals likely sense the odor concentration in air phase but not in water phase (i.e., agar surface) because of the high hydrophobicity of 2-nonanone (a nine carbon ketone).

Quantitative behavioral experiments on the OSB2 system

Request a detailed protocol

The behavior of animals was calculated from records of displacement of the auto-tracking stage and from the position of the ROI or centroid of the view field. For ROI-tracking, the trajectory of an animal’s behavior was wavier than for centroid tracking because the ROI was usually set around the animal's head, which moves in a sinusoidal pattern (Video 2). To compensate for the wavy pattern, the $x - y$ coordinates for ROI-tracking were calculated as a moving-average for ±10 frames at 10 Hz (i.e., ±1 s). This gave similar results to centroid-tracking on the quantitative behavioral analysis. The migratory trajectory from either of the tracking methods was sampled at each second (1 Hz), and a change in the migratory vector for 1 s larger than 90° was recognized as a turn. In the Figures, ensemble averages in each 10 s bin are shown.

In Figure 4A and B, we investigated the time when the turning rate changed based on the rate of increase or decrease in odor concentration. In order to investigate the timing, it was necessary to finely set the time window. However, since a turn is an uncommon occurrence (a turning rate of 0.1 is once in 10 s), narrowing the time window increased the variation. In order to obtain the same number of turns as the 60 s using a time window of 10 s, six times as much sampling had to be performed. Furthermore, even more samples were required for performing multiple tests. Therefore, we used the prediction interval, a criterion in the field of statistical inference. The 99% prediction interval is an interval in which future data will fall with 99% probability, if it obeys the same probability distribution as the previously observed data (in this case, odor-zero or odor-plateau phase). A 100(1-α)% prediction interval on a single future observation $(X_{n + 1})$ from a normal distribution is given by the following formula:

\bar{x} - t_{\frac{α}{2}, n - 1} s \sqrt{1 + \frac{1}{n}} \leq X_{n + 1} \leq \bar{x} + t_{\frac{α}{2}, n - 1} s \sqrt{1 + \frac{1}{n}}

where $\bar{x}$ is the sample mean, $n$ is the number of previously observed data, $t_{α / 2, n - 1}$ is the 100(1-α/2) percentage point of a t-distribution with $n - 1$ degrees of freedom, and s is the sample standard deviation (Montgomery and Runger, 2002). Using this criterion, we analyzed ‘timing when the unexpected value appears for the first time’ in Figure 4A and B.

Calcium imaging

View detailed protocol

The details of calcium imaging with the OSB system were previously described (Tanimoto et al., 2016). In brief, the sample was exposed to excitation light from a MiLSS (Multi-independent Light Stimulation System, Aska Company, Japan) (Sakai et al., 2013). The images for GCaMP and mCherry were split and simultaneously captured side-by-side on an EM-CCD camera ImagEM with W-View system (Hamamatsu, Japan). Images were taken at a 32.6 ms exposure time and 100 ms sampling interval with 2 × 2 binning. The cell body was tracked off-line with another custom-made program for the centering (Video 2), and signal intensities of particular regions were measured by ImageJ (NIH). The data from frames where the cell body was not centered were omitted. The signal intensity of the background was subtracted from that of the cell body, and the value was moving-averaged for ±1 frames and further analysed. The average of fluorescence intensity of GCaMP during 1 min before the odor increment or decrement was defined as the baseline F₀. Because ∆F/F₀ of GCaMP and the ratio between fluorescence intensities of GCaMP and mCherry (GCaMP/mCherry: R) exhibited similar tendencies, and because and the noise level was smaller in ∆F/F₀ than in R, the data of ∆F/F₀ were used in the figures. In Figure 7, ∆R was used because the mutations in itr-1 or unc-68 could affect the baseline as well as the response calcium levels of the neurons. Also in Figure 7, the animals were immobilized with the acetylcholine receptor agonist levamisole for high-throughput analysis, in which multiple animals were stimulated and imaged simultaneously. Even with the levamisole treatment, the responses of AWB and ASH neurons in the naive wild-type animals were essentially similar to those in the freely moving animals (Figure 4A and B).

Optogenetic analysis

Request a detailed protocol

Animals were raised in the presence or absence of ATR according to the previous report (Kawazoe et al., 2013), and transferred to an NGM plate on the OSB2 system and maintained under the objective lens by auto-tracking. For ChR2(C128S) experiments in the absence of a 2-nonanone stimulus (Figure 3D), after 1 min without light stimulation, the animal was transiently illuminated with blue light (3 s) for activation through BP460-495 and DM505 with ND25 (~0.8 mW/mm²). Turning rates of 30–60 s and 65–95 s were calculated as before or after the blue light illumination, respectively. The turns of 60–65 s were not included in the calculation because blue light illumination (60–63 s) appeared to somewhat affect the animals’ locomotion for a few seconds (Ward et al., 2008). For Arch experiments in the presence of a 2-nonanone stimulus (Figures 3F and 4C), green light was delivered through BP530-550 and DM570 at ~1.0 mW/mm², and turning rates were calculated. The optical filters were from Olympus.

Mathematical modeling of neuronal responses

Request a detailed protocol

For the time-differential models of neuronal responses, the following time-differential equation was used:

X (t) = k \frac{d C (t)}{d t}

where $X (t)$ is neuronal response, $k$ is the conversion factor, and $C (t)$ is the measured odor concentration. The $d C (t) / d t$ was calculated as the central difference of $C (t)$ . This equation indicates that the neuronal response $X (t)$ responds to the odor gradient $d C (t) / d t$ at each time. The value of $k$ was determined by the least squares method to fit $X (t)$ to the measured $Δ F / F_{0}$ in response to the odor gradients.

For the time-integral models of neuronal responses, the following leaky integrator equation was used:

\frac{d X (t)}{d t} = k I (t) - \frac{1}{τ} X (t)

where external input was given by temporal odor change; $I (t) = d C (t) / d t$ . $τ$ is the time constant of leaky integration. $k$ and $τ$ were determined by the least squares method to fit $X (t)$ to the measured $Δ F / F_{0}$ responded to the odor gradients. This differential equation was numerically integrated by the Euler method with a time-step of 1 s. The initial value was $X (t)$ = 0 which corresponds to $Δ F / F_{0}$ = 0 in the basal state. For odr-3 mutants, on the other hand, external input was $I (t) = - (C (t) - C (t - Δ t)) / Δ t$ in the leaky integrator equation, and the values $k, τ$ , and $Δ t$ were determined to fit X(t) to the measured $Δ F / F_{0}$ of odr-3.

Estimation of intracellular calcium concentration in Figure 4—figure supplement 2 was conducted as follows: Since the relationship between fluorescence signals and calcium concentration is non-linear, a change in the neuronal activity to stimulation is properly evaluated, not by the fluorescence intensity of the calcium indicator, but by the calcium concentration itself. Taking the non-linear relationship into account, intracellular calcium concentration [Ca²⁺] was estimated by the Hill equation; $(F - F_{m i n}) / (F_{m a x} - F_{m i n}) = [C a^{2 +}]^{h} / ([C a^{2 +}]^{h} + {K_{d}}^{h})$ . The $F$ is the measured fluorescence intensity, $F_{m i n}$ and $F_{m a x}$ are the fluorescence intensities under Ca²⁺-free and Ca²⁺-saturated conditions, respectively. The $h$ is the Hill coefficient and $K_{d}$ is the dissociation constant. For GCaMP3, the values of $h$ and $K_{d}$ were reported previously (Akerboom et al., 2012). In each experiment, Ca²⁺ response to stimulation is expressed as the ratio of the fluorescence response to the basal fluorescence intensity $F_{0}, Δ F / F_{0} = (F - F_{0}) / F_{0}$ . By solving the Hill equation for [Ca²⁺] in terms of the fluorescence intensities, the following equation to calculate the intracellular calcium concentration from the measured ratio $Δ F / F_{0}$ was obtained:

[C a^{2 +}] = K_{d} {(\frac{1 + Δ F / F_{0} - f_{m i n}}{f_{m a x} - 1 - Δ F / F_{0}})}^{1 / h}

where $f_{m i n} = F_{m i n} / F_{0}$ and $f_{m a x} = F_{m a x} / F_{0}$ are the minimum and maximum fluorescence intensities relative to $F_{0}$ , respectively. For GCaMP3, $f_{m a x} = 12 f_{m i n}$ since the dynamic range $F_{m a x} / F_{m i n}$ (i.e., $f_{m a x} / f_{m i n}$ ) is reported to be ~12 fold (Tian et al., 2009). The time delay of fluorescence response to a calcium concentration change was not taken into account since the temporal resolution of the odor concentration measurement was of the second order, while the association and dissociation time constants of GCaMP3 are of the sub-second order (Tian et al., 2009). When $X (t)$ corresponds to the calcium concentration, the basal value of $X (t)$ in the steady state is not zero since the intracellular calcium concentration is not reduced to zero even in the basal state. Therefore, the leaky integrator equation was generalized as follows:

\frac{d X (t)}{d t} = k I (t) - \frac{1}{τ} (X (t) - X_{b a s e})

where $X_{b a s e}$ corresponds to the basal calcium concentration in the steady state and takes a positive value. For AWB and ASH neurons, unknown model parameters $k, τ, f_{m i n}$ and $X_{b a s e}$ were determined by the least squares method to fit $X (t)$ calculated by the generalized leaky integrator equation to the calcium concentrations estimated from $Δ F / F_{0}$ . Similar estimation of non-linear property of GCaMP3 has been reported previously (Kato et al., 2014). The time-differential and time-integral models reasonably approximated the neural responses under the conditions used in this study. However, with stronger odor concentration changes, input saturation may need to be considered, in which case the input could be put through a logistic sigmoid function for example. The values of the fitting parameters are shown in Tables 1, 2, 3 and 4.

Computer simulation of 2-nonanone avoidance behavior

Request a detailed protocol

The previous algorithms (Iino and Yoshida, 2009; Yamazoe-Umemoto et al., 2015) were modified as follows to simulate 2-nonanone avoidance behavior (Figure 5). The parameters for simulation were based on the migratory statistics of real wild-type animals and contained no free parameters unless otherwise indicated. The model animal moved at a speed of 0.14 mm/s. In the low-turning state, the model animal moved forward with fluctuations in migratory direction, which was randomly chosen from the Gaussian distribution of −0.065 ± 5.14° (mean ± SD). The odor signal periodically fluctuated because of the sinusoidal movement of the animal. The position of animal’s anterior end, where the sensory endings of ASH and AWB neurons are located, was calculated as a sine curve along the animal’s track. The amplitude and frequency of the sine curve was 0.1 mm and 0.5 Hz, respectively (Kimura et al., 2004; Shen et al., 2012). The track of the anterior end was used for the calculation of C_worm. A turn occurred based on the pirouette initiation rate of 0.0326/ (0.200 + exp(−231 ×dC/dt))+0.0260, which is relatively constant (~0.03 s⁻¹) when dC/dt <0 and increases when dC/dt >0; The dC/dt-dependency in the pirouette initiation rate was determined from the probability of pirouette initiation after 2 s of the step for real animals. The turning duration was 3 s. After a turn, the model animal was in the high-turning state and initiated a migration, whose deviation in direction from the direction just before the turn was randomly chosen from a pool of the measured values in real animals. In the high-turning state, the model animals turned at a constant rate of 0.2 s⁻¹, which results in ~95% of migratory duration shorter than the threshold value t_crit (13.1 s). Therefore, in the high-turning state, most of the migrations were classified as pirouettes. When the model animals happened to migrate down the gradient and experienced dC/dt <0, their state transited from high- to low-turning according to the leaky integration of dC/dt described in the previous section. When the leaky integration of dC/dt became higher than 0.18, the high turning state was switched to the low turning state. The threshold value 0.18 was chosen as a value similar to the one associated with the turn suppression in Figure 4B. For the ‘differentiation model,’ dC/dt itself was used for the initiation of a low turning state instead of leaky integration, and the high turning state was switched to the low turning state when dC/dt was negative. The simulation was repeated 100 times for each model animal condition. Time was discretized into intervals with ∆t = 1 s.

Data analysis and statistics

Request a detailed protocol

For the experiments with the OSB2 system, the data were obtained on multiple days from approximately 20–50 animals for each condition. We chose this sample number because a large scale behavioral analysis of C. elegans concluded that 20 animals would discriminate single SD in a behavioral phenotype at over 80% power, and 24 ± 14 (average ± SD) animals per condition were used in the study (Yemini et al., 2013). For Figures 1 and 2 and Figure 1—figure supplements 1, 100 animals were used because we investigated various aspects of behavior in detail. After the sample acquisition, the data of some animals for Figures 3, 4, 6 and 7 were excluded when any of the following problems were found: (1) trials interrupted by errors in auto-tracking, (2) animals with too weak intensity of basal mCherry or GCaMP3 fluorescence for off-line tracking, (3) animals with abnormal sudden transient activation of AWB during odor-zero or odor-up phases (8 out of 279 animals tested for AWB), or (4) animals with basal locomotion speeds slower than 0.02 mm/s (average speed ± SD of normal animals was 0.15 ± 0.04 mm/s). Experimental conditions, such as the presence/absence of ATR, light stimulation, odor gradient, or different strains were randomized on a daily basis.

A Kruskal-Wallis test with a post hoc Steel-Dwass test was used for multiple comparisons in Figures 3B, 4A, B and 7C, while a Mann-Whitney test was used for single comparison in Figures 1D, 3D, F and 4C, and Figure 2—figure supplement 1D (right panel) using R (The R Project) or Prism ver. 5.0 for Mac OSX (GraphPad Software, San Diego, CA). The Mardia-Watson-Wheeler test was used for Figure 5C and Watson's U² test was used for Figure 3—figure supplement 1B by using Oriana ver. 3 (Kovach Computing Service, Wales, UK). All the statistical details are shown in Supplementary file 1.

The Bayesian information criterion (BIC) was used to assess mathematical model fitting in Figures 4 and 6. In BIC, the goodness of fit for the model including a penalty term to discourage overfitting is given by the following equation:

B I C = N l n (\frac{R S S}{N}) + M l n (N)

where N is the number of samples (frames) used for the fitting, RSS is the residual sum of squares obtained from fitting by the least squares method, and M is the number of free parameters in a given model, respectively. M = 1 for the differentiation model (k), M = 2 for the leaky integration model (k, τ), and M = 3 for the leaky integration model for odr-3 (k, τ, ∆t). The lower BIC value means a better fitting with a model.

Data availability

The following data sets were generated

(2017) Data from: Calcium dynamics regulating the timing of decision-making in C. elegans
Available at Dryad Digital Repository under a CC0 Public Domain Dedication.

http://dx.doi.org/10.5061/dryad.h3680

References

1. Akerboom J
2. Chen TW
3. Wardill TJ
4. Tian L
5. Marvin JS
6. Mutlu S
7. Calderón NC
8. Esposti F
9. Borghuis BG
10. Sun XR
11. Gordus A
12. Orger MB
13. Portugues R
14. Engert F
15. Macklin JJ
16. Filosa A
17. Aggarwal A
18. Kerr RA
19. Takagi R
20. Kracun S
21. Shigetomi E
22. Khakh BS
23. Baier H
24. Lagnado L
25. Wang SS
26. Bargmann CI
27. Kimmel BE
28. Jayaraman V
29. Svoboda K
30. Kim DS
31. Schreiter ER
32. Looger LL
(2012) Optimization of a GCaMP calcium Indicator for neural activity imaging
Journal of Neuroscience 32:13819–13840.

https://doi.org/10.1523/JNEUROSCI.2601-12.2012
- PubMed
- Google Scholar
1. Aksay E
2. Olasagasti I
3. Mensh BD
4. Baker R
5. Goldman MS
6. Tank DW
(2007) Functional dissection of circuitry in a neural integrator
Nature Neuroscience 10:494–504.

https://doi.org/10.1038/nn1877
- PubMed
- Google Scholar
1. Alon U
(2007) Network motifs: theory and experimental approaches
Nature Reviews Genetics 8:450–461.

https://doi.org/10.1038/nrg2102
- PubMed
- Google Scholar
1. Bargmann C
(2006) Chemosensation in C. elegans
WormBook pp. 1–29.

https://doi.org/10.1895/wormbook.1.123.1
- Google Scholar
(1993) Odorant-selective genes and neurons mediate olfaction in C. elegans
Cell 74:515–527.

https://doi.org/10.1016/0092-8674(93)80053-H
- PubMed
- Google Scholar
(2011) Catecholamine receptor polymorphisms affect decision-making in C. elegans
Nature 472:313–318.

https://doi.org/10.1038/nature09821
- PubMed
- Google Scholar
(2009) Bi-stable neural state switches
Nature Neuroscience 12:9–234.

https://doi.org/10.1038/nn.2247
- PubMed
- Google Scholar
(2008) An olfactory neuron responds stochastically to temperature and modulates Caenorhabditis elegans thermotactic behavior
PNAS 105:11002–11007.

https://doi.org/10.1073/pnas.0805004105
- PubMed
- Google Scholar
1. Brenner S
(1974)
The genetics of Caenorhabditis elegans

Genetics 77:4.
- Google Scholar
(2005) Optical imaging of neuronal populations during decision-making
Science 307:896–901.

https://doi.org/10.1126/science.1103736
- PubMed
- Google Scholar
1. Buhusi CV
2. Meck WH
(2005) What makes Us tick? functional and neural mechanisms of interval timing
Nature Reviews Neuroscience 6:755–765.

https://doi.org/10.1038/nrn1764
- PubMed
- Google Scholar
1. Busch KE
2. Laurent P
3. Soltesz Z
4. Murphy RJ
5. Faivre O
6. Hedwig B
7. Thomas M
8. Smith HL
9. de Bono M
(2012) Tonic signaling from O₂ sensors sets neural circuit activity and behavioral state
Nature Neuroscience 15:581–591.

https://doi.org/10.1038/nn.3061
- PubMed
- Google Scholar
1. Carandini M
2. Churchland AK
(2013) Probing perceptual decisions in rodents
Nature Neuroscience 16:824–831.

https://doi.org/10.1038/nn.3410
- PubMed
- Google Scholar
1. Catterall WA
(2011) Voltage-gated calcium channels
Cold Spring Harbor Perspectives in Biology 3:a003947–23.

https://doi.org/10.1101/cshperspect.a003947
- PubMed
- Google Scholar
1. Chow BY
2. Han X
3. Dobry AS
4. Qian X
5. Chuong AS
6. Li M
7. Henninger MA
8. Belfort GM
9. Lin Y
10. Monahan PE
11. Boyden ES
(2010) High-performance genetically targetable optical neural silencing by light-driven proton pumps
Nature 463:98–102.

https://doi.org/10.1038/nature08652
- PubMed
- Google Scholar
1. Clapham DE
(2007) Calcium signaling
Cell 131:1047–1058.

https://doi.org/10.1016/j.cell.2007.11.028
- PubMed
- Google Scholar
1. Colosimo ME
2. Brown A
3. Mukhopadhyay S
4. Gabel C
5. Lanjuin AE
6. Samuel AD
7. Sengupta P
(2004) Identification of thermosensory and olfactory neuron-specific genes via expression profiling of single neuron types
Current Biology 14:2245–2251.

https://doi.org/10.1016/j.cub.2004.12.030
- PubMed
- Google Scholar
1. Curtis CE
2. Lee D
(2010) Beyond working memory: the role of persistent activity in decision making
Trends in Cognitive Sciences 14:216–222.

https://doi.org/10.1016/j.tics.2010.03.006
- PubMed
- Google Scholar
(1999) The inositol trisphosphate receptor regulates a 50-second behavioral rhythm in C. elegans
Cell 98:757–767.

https://doi.org/10.1016/S0092-8674(00)81510-X
- PubMed
- Google Scholar
(2014) FoxP influences the speed and accuracy of a perceptual decision in Drosophila
Science 344:901–904.

https://doi.org/10.1126/science.1252114
- PubMed
- Google Scholar
1. de Bono M
2. Maricq AV
(2005) Neuronal substrates of complex behaviors in C. elegans
Annual Review of Neuroscience 28:1–501.

https://doi.org/10.1146/annurev.neuro.27.070203.144259
- PubMed
- Google Scholar
1. Delcomyn F
(1998)
Foundations of Neurobiology

209–231, Properties of Sensory Systems, Foundations of Neurobiology, New York, W H Freeman & Company.
- Google Scholar
(2006) Effects of voltage-gated calcium channel subunit genes on calcium influx in cultured C. elegans mechanosensory neurons
Journal of Neurobiology 66:1125–1139.

https://doi.org/10.1002/neu.20261
- PubMed
- Google Scholar
1. Gold JI
2. Shadlen MN
(2007) The neural basis of decision making
Annual Review of Neuroscience 30:535–574.

https://doi.org/10.1146/annurev.neuro.29.051605.113038
- PubMed
- Google Scholar
(2011) Active sampling and decision making in Drosophila chemotaxis
Nature Communications 2:441.

https://doi.org/10.1038/ncomms1455
- PubMed
- Google Scholar
(2005) A circuit for navigation in Caenorhabditis elegans
PNAS 102:3184–3191.

https://doi.org/10.1073/pnas.0409009101
- PubMed
- Google Scholar
1. Ha HI
2. Hendricks M
3. Shen Y
4. Gabel CV
5. Fang-Yen C
6. Qin Y
7. Colón-Ramos D
8. Shen K
9. Samuel AD
10. Zhang Y
(2010) Functional organization of a neural network for aversive olfactory learning in Caenorhabditis elegans
Neuron 68:1173–1186.

https://doi.org/10.1016/j.neuron.2010.11.025
- PubMed
- Google Scholar
1. Hille B
(2001)
Ion Channels of Excitable Membranes

Voltage-Gated Calcium Channels, Ion Channels of Excitable Membranes, Sunderland MA, Sinauer.
- Google Scholar
(2005) In vivo imaging of C. elegans ASH neurons: cellular response and adaptation to chemical repellents
The EMBO Journal 24:63–72.

https://doi.org/10.1038/sj.emboj.7600493
- PubMed
- Google Scholar
1. Iino Y
2. Yoshida K
(2009) Parallel use of two behavioral mechanisms for chemotaxis in Caenorhabditis elegans
Journal of Neuroscience 29:5370–5380.

https://doi.org/10.1523/JNEUROSCI.3633-08.2009
- PubMed
- Google Scholar
1. Kaplan JM
(1996) Sensory signaling in Caenorhabditis elegans
Current Opinion in Neurobiology 6:494–499.

https://doi.org/10.1016/S0959-4388(96)80055-9
- PubMed
- Google Scholar
1. Kato S
2. Xu Y
3. Cho CE
4. Abbott LF
5. Bargmann CI
(2014) Temporal responses of C. elegans chemosensory neurons are preserved in behavioral dynamics
Neuron 81:616–628.

https://doi.org/10.1016/j.neuron.2013.11.020
- PubMed
- Google Scholar
(2013) A simple optogenetic system for behavioral analysis of freely moving small animals
Neuroscience Research 75:8.

https://doi.org/10.1016/j.neures.2012.04.011
- PubMed
- Google Scholar
(2010) Enhancement of odor avoidance regulated by dopamine signaling in Caenorhabditis elegans
Journal of Neuroscience 30:16365–16375.

https://doi.org/10.1523/JNEUROSCI.6023-09.2010
- PubMed
- Google Scholar
(2004) The C. elegans thermosensory neuron AFD responds to warming
Current Biology 14:1291–1295.

https://doi.org/10.1016/j.cub.2004.06.060
- PubMed
- Google Scholar
(2004) Dantrolene--a review of its pharmacology, therapeutic use and new developments
Anaesthesia 59:364–373.

https://doi.org/10.1111/j.1365-2044.2004.03658.x
- PubMed
- Google Scholar
1. Kristan WB
(2008) Neuronal decision-making circuits
Current Biology 18:R928–R932.

https://doi.org/10.1016/j.cub.2008.07.081
- PubMed
- Google Scholar
1. Kwok TC
2. Ricker N
3. Fraser R
4. Chan AW
5. Burns A
6. Stanley EF
7. McCourt P
8. Cutler SR
9. Roy PJ
(2006) A small-molecule screen in C. elegans yields a new calcium channel antagonist
Nature 441:91–95.

https://doi.org/10.1038/nature04657
- PubMed
- Google Scholar
1. Larsch J
2. Flavell SW
3. Liu Q
4. Gordus A
5. Albrecht DR
6. Bargmann CI
(2015) A circuit for Gradient Climbing in C. elegans Chemotaxis
Cell Reports 12:1748–1760.

https://doi.org/10.1016/j.celrep.2015.08.032
- PubMed
- Google Scholar
1. Lee RY
2. Lobel L
3. Hengartner M
4. Horvitz HR
5. Avery L
(1997) Mutations in the alpha1 subunit of an L-type voltage-activated ^Ca2+ channel cause myotonia in Caenorhabditis elegans
The EMBO Journal 16:6066–6076.

https://doi.org/10.1093/emboj/16.20.6066
- PubMed
- Google Scholar
1. Li Z
2. Liu J
3. Zheng M
4. Xu XZ
(2014) Encoding of both analog- and digital-like behavioral outputs by one C. elegans interneuron
Cell 159:751–765.

https://doi.org/10.1016/j.cell.2014.09.056
- PubMed
- Google Scholar
1. Lockery SR
(2011) The computational worm: spatial orientation and its neuronal basis in C. elegans
Current Opinion in Neurobiology 21:782–790.

https://doi.org/10.1016/j.conb.2011.06.009
- PubMed
- Google Scholar
1. Loewenstein Y
2. Sompolinsky H
(2003) Temporal integration by calcium dynamics in a model neuron
Nature Neuroscience 6:961–967.

https://doi.org/10.1038/nn1109
- PubMed
- Google Scholar
1. Louis M
2. Huber T
3. Benton R
4. Sakmar TP
5. Vosshall LB
(2008) Bilateral olfactory sensory input enhances chemotaxis behavior
Nature Neuroscience 11:187–199.

https://doi.org/10.1038/nn2031
- PubMed
- Google Scholar
1. Luo L
2. Wen Q
3. Ren J
4. Hendricks M
5. Gershow M
6. Qin Y
7. Greenwood J
8. Soucy ER
9. Klein M
10. Smith-Parker HK
11. Calvo AC
12. Colón-Ramos DA
13. Samuel AD
14. Zhang Y
(2014) Dynamic encoding of perception, memory, and movement in a C. elegans chemotaxis circuit
Neuron 82:1115–1128.

https://doi.org/10.1016/j.neuron.2014.05.010
- PubMed
- Google Scholar
1. Major G
2. Tank D
(2004) Persistent neural activity: prevalence and mechanisms
Current Opinion in Neurobiology 14:675–684.

https://doi.org/10.1016/j.conb.2004.10.017
- PubMed
- Google Scholar
Conference
1. Maru M
2. Igarashi Y
3. Arai S
4. Hashimoto K
(2010)
Fluorescent microscope system to track a particular region of C. elegans

IEEE/SICE International Symposium on System Integration.
- Google Scholar
(1991)
Efficient gene transfer in C.elegans: extrachromosomal maintenance and integration of transforming sequences

The EMBO Journal 10:3959–3970.
- PubMed
- Google Scholar
Book
1. Montgomery DC
2. Runger GC
(2002)
Applied Statistics and Probability for Engineers

New York: Wiley.
- Google Scholar
(2005) Light activation of channelrhodopsin-2 in excitable cells of Caenorhabditis elegans triggers rapid behavioral responses
Current Biology 15:2279–2284.

https://doi.org/10.1016/j.cub.2005.11.032
- PubMed
- Google Scholar
(1999)
The fundamental role of pirouettes in Caenorhabditis elegans chemotaxis

Journal of Neuroscience 19:7–9569.
- PubMed
- Google Scholar
1. Pokala N
2. Liu Q
3. Gordus A
4. Bargmann CI
(2014) Inducible and titratable silencing of Caenorhabditis elegans neurons in vivo with histamine-gated chloride channels
PNAS 111:2770–2775.

https://doi.org/10.1073/pnas.1400615111
- PubMed
- Google Scholar
Book
(1992)
Numerical Recipes in C

Cambridge: Cambridge University Press.
- Google Scholar
Book
1. Richelle M
2. Lejeune H
(1980)
Time in Animal Behaviour

Oxford: Pergamon Press.
- Google Scholar
1. Richmond J
(2005) Synaptic function
WormBook pp. 1–14.

https://doi.org/10.1895/wormbook.1.69.1
- PubMed
- Google Scholar
(1998) The Gα protein ODR-3 mediates olfactory and nociceptive function and Controls Cilium Morphogenesis in C. elegans olfactory neurons
Neuron 20:55–67.

https://doi.org/10.1016/S0896-6273(00)80434-1
- Google Scholar
1. Sakai S
2. Ueno K
3. Ishizuka T
4. Yawo H
(2013) Parallel and patterned optogenetic manipulation of neurons in the brain slice using a DMD-based projector
Neuroscience Research 75:59–64.

https://doi.org/10.1016/j.neures.2012.03.009
- PubMed
- Google Scholar
1. Sakube Y
2. Ando H
3. Kagawa H
(1997) An abnormal ketamine response in mutants defective in the ryanodine receptor gene ryr-1 (unc-68) of Caenorhabditis elegans
Journal of Molecular Biology 267:849–864.

https://doi.org/10.1006/jmbi.1997.0910
- PubMed
- Google Scholar
1. Sasakura H
2. Mori I
(2013) Behavioral plasticity, learning, and memory in C. elegans
Current Opinion in Neurobiology 23:92–99.

https://doi.org/10.1016/j.conb.2012.09.005
- PubMed
- Google Scholar
1. Schafer WR
2. Kenyon CJ
(1995) A calcium-channel homologue required for adaptation to dopamine and serotonin in Caenorhabditis elegans
Nature 375:73–78.

https://doi.org/10.1038/375073a0
- PubMed
- Google Scholar
1. Schall JD
(2001) Neural basis of deciding, choosing and acting
Nature Reviews Neuroscience 2:33–42.

https://doi.org/10.1038/35049054
- PubMed
- Google Scholar
1. Schultz W
(2007) Behavioral dopamine signals
Trends in Neurosciences 30:203–210.

https://doi.org/10.1016/j.tins.2007.03.007
- PubMed
- Google Scholar
1. Schulze A
2. Gomez-Marin A
3. Rajendran VG
4. Lott G
5. Musy M
6. Ahammad P
7. Deogade A
8. Sharpe J
9. Riedl J
10. Jarriault D
11. Trautman ET
12. Werner C
13. Venkadesan M
14. Druckmann S
15. Jayaraman V
16. Louis M
(2015) Dynamical feature extraction at the sensory periphery guides chemotaxis
eLife 4:1129.

https://doi.org/10.7554/eLife.06694
- PubMed
- Google Scholar
1. Shadlen MN
2. Newsome WT
(2001)
Neural basis of a perceptual decision in the parietal cortex (area LIP) of the rhesus monkey

Journal of Neurophysiology 86:6–1936.
- PubMed
- Google Scholar
(2004) Improved monomeric red, orange and yellow fluorescent proteins derived from Discosoma sp. red fluorescent protein
Nature Biotechnology 22:1567–1572.

https://doi.org/10.1038/nbt1037
- PubMed
- Google Scholar
(2012) Undulatory locomotion of Caenorhabditis elegans on wet surfaces
Biophysical Journal 102:2772–2781.

https://doi.org/10.1016/j.bpj.2012.05.012
- PubMed
- Google Scholar
1. Steger KA
2. Shtonda BB
3. Thacker C
4. Snutch TP
5. Avery L
(2005) The C. elegans T-type calcium channel CCA-1 boosts neuromuscular transmission
Journal of Experimental Biology 208:2191–2203.

https://doi.org/10.1242/jeb.01616
- PubMed
- Google Scholar
1. Suzuki H
2. Kerr R
3. Bianchi L
4. Frøkjaer-Jensen C
5. Slone D
6. Xue J
7. Gerstbrein B
8. Driscoll M
9. Schafer WR
(2003) In vivo imaging of C. elegans mechanosensory neurons demonstrates a specific role for the MEC-4 channel in the process of gentle touch sensation
Neuron 39:1005–1017.

https://doi.org/10.1016/j.neuron.2003.08.015
- PubMed
- Google Scholar
Data
1. Tanimoto Y
2. Yamazoe-Umemoto A
3. Fujita K
4. Kawazoe Y
5. Miyanishi Y
6. Yamazaki SJ
7. Fei X
8. Busch KE
9. Gengyo-Ando K
10. Nakai J
11. Iino Y
12. Iwasaki Y
13. Hashimoto K
14. Kimura KD
(authors) (2017) Data from: Calcium dynamics regulating the timing of decision-making in C. elegans
Dryad Digital Repository.

https://doi.org/10.5061/dryad.h3680
1. Tanimoto Y
2. Zheng YG
3. Fei X
4. Fujie Y
5. Hashimoto K
6. Kimura KD
(2016) In actio optophysiological analyses reveal functional diversification of dopaminergic neurons in the nematode C. elegans
Scientific Reports 6:7.

https://doi.org/10.1038/srep26297
- PubMed
- Google Scholar
1. Tian L
2. Hires SA
3. Mao T
4. Huber D
5. Chiappe ME
6. Chalasani SH
7. Petreanu L
8. Akerboom J
9. McKinney SA
10. Schreiter ER
11. Bargmann CI
12. Jayaraman V
13. Svoboda K
14. Looger LL
(2009) Imaging neural activity in worms, flies and mice with improved GCaMP calcium indicators
Nature Methods 6:875–881.

https://doi.org/10.1038/nmeth.1398
- PubMed
- Google Scholar
(1995) Divergent seven transmembrane receptors are candidate chemosensory receptors in C. elegans
Cell 83:207–218.

https://doi.org/10.1016/0092-8674(95)90162-0
- PubMed
- Google Scholar
(1997) Reprogramming chemotaxis responses: sensory neurons define olfactory preferences in C. elegans
Cell 91:161–169.

https://doi.org/10.1016/S0092-8674(00)80399-2
- PubMed
- Google Scholar
(2012) A model of the intracellular response of an olfactory neuron in Caenorhabditis elegans to odor stimulation
PLoS One 7:e42907.

https://doi.org/10.1371/journal.pone.0042907
- PubMed
- Google Scholar
1. Wang XJ
(2008) Decision making in recurrent neuronal circuits
Neuron 60:215–234.

https://doi.org/10.1016/j.neuron.2008.09.034
- PubMed
- Google Scholar
1. Ward A
2. Liu J
3. Feng Z
4. Xu XZ
(2008) Light-sensitive neurons and channels mediate phototaxis in C. elegans
Nature Neuroscience 11:916–922.

https://doi.org/10.1038/nn.2155
- PubMed
- Google Scholar
(1986) The structure of the nervous system of the nematode Caenorhabditis elegans
Philosophical Transactions of the Royal Society B: Biological Sciences 314:1–340.

https://doi.org/10.1098/rstb.1986.0056
- PubMed
- Google Scholar
(2015) Modulation of different behavioral components by neuropeptide and dopamine signalings in non-associative odor learning of Caenorhabditis elegans
Neuroscience Research 99:22–33.

https://doi.org/10.1016/j.neures.2015.05.009
- PubMed
- Google Scholar
(2018) Measuring Spatiotemporal Dynamics of Odor Gradient for Small Animals by Gas Chromatography
BIO-PROTOCOL 8:e2797.

https://doi.org/10.21769/BioProtoc.2797
- Google Scholar
1. Yang CH
2. Belawat P
3. Hafen E
4. Jan LY
5. Jan YN
(2008) Drosophila egg-laying site selection as a system to study simple decision-making processes
Science 319:1679–1683.

https://doi.org/10.1126/science.1151842
- PubMed
- Google Scholar
1. Yemini E
2. Jucikas T
3. Grundy LJ
4. Brown AE
5. Schafer WR
(2013) A database of Caenorhabditis elegans behavioral phenotypes
Nature Methods 10:877–879.

https://doi.org/10.1038/nmeth.2560
- PubMed
- Google Scholar
(2015) Serotonin differentially modulates Ca²⁺ transients and depolarization in a C. elegans nociceptor
Journal of Neurophysiology 113:1041–1050.

https://doi.org/10.1152/jn.00665.2014
- PubMed
- Google Scholar

Article and author information

Author details

Yuki Tanimoto

Department of Biological Sciences, Graduate School of Science, Osaka University, Toyonaka, Japan

Contribution
YT, Conceptualization, Investigation, Methodology, Writing—original draft, Writing—review and editing

Contributed equally with
Akiko Yamazoe-Umemoto

Competing interests
The authors declare that no competing interests exist.

"This ORCID iD identifies the author of this article:" 0000-0002-7434-678X
Akiko Yamazoe-Umemoto

Department of Biological Sciences, Graduate School of Science, Osaka University, Toyonaka, Japan

Contribution
AY-U, Conceptualization, Investigation, Methodology, Writing—original draft

Contributed equally with
Yuki Tanimoto

Competing interests
The authors declare that no competing interests exist.
Kosuke Fujita

Department of Biological Sciences, Graduate School of Science, Osaka University, Toyonaka, Japan

Present address
Department of Ophthalmology, Graduate School of Medicine, Tohoku University, Sendai, Japan

Contribution
KF, Investigation, Writing—review and editing

Competing interests
The authors declare that no competing interests exist.
Yuya Kawazoe

Department of Biological Sciences, Graduate School of Science, Osaka University, Toyonaka, Japan

Contribution
YK, Investigation, Writing—review and editing

Competing interests
The authors declare that no competing interests exist.
Yosuke Miyanishi

Department of Biological Sciences, Graduate School of Science, Osaka University, Toyonaka, Japan

Contribution
YM, Investigation, Writing—review and editing

Competing interests
The authors declare that no competing interests exist.
Shuhei J Yamazaki

Department of Biological Sciences, Graduate School of Science, Osaka University, Toyonaka, Japan

Contribution
SJY, Validation, Investigation, Writing—review and editing

Competing interests
The authors declare that no competing interests exist.
Xianfeng Fei

Faculty of Science and Technology, Tohoku Bunka Gakuen University, Sendai, Japan

Contribution
XF, Software, Investigation, Methodology, Writing—review and editing

Competing interests
The authors declare that no competing interests exist.
Karl Emanuel Busch

Centre for Integrative Physiology, The University of Edinburgh, Edinburgh, United Kingdom

Contribution
KEB, Investigation, Methodology, Writing—review and editing

Competing interests
The authors declare that no competing interests exist.
Keiko Gengyo-Ando

Graduate Shool of Science and Engineering, Brain and Body System Science Institute, Saitama University, Saitama, Japan

Contribution
KG-A, Funding acquisition, Investigation, Methodology, Writing—review and editing

Competing interests
The authors declare that no competing interests exist.
Junichi Nakai

Graduate Shool of Science and Engineering, Brain and Body System Science Institute, Saitama University, Saitama, Japan

Contribution
JN, Funding acquisition, Investigation, Methodology, Writing—review and editing

Competing interests
The authors declare that no competing interests exist.
Yuichi Iino

Department of Biological Sciences, Graduate School of Science, The University of Tokyo, Tokyo, Japan

Contribution
YIi, Software, Funding acquisition, Investigation, Methodology, Writing—review and editing

Competing interests
The authors declare that no competing interests exist.
Yuishi Iwasaki

Department of Intelligent Systems Engineering, Ibaraki University, Hitachi, Japan

Contribution
YIw, Software, Validation, Investigation, Visualization, Methodology, Writing—review and editing

Competing interests
The authors declare that no competing interests exist.
Koichi Hashimoto

Graduate School of Information Sciences, Tohoku University, Sendai, Japan

Contribution
KH, Software, Funding acquisition, Investigation, Methodology, Writing—review and editing

Competing interests
The authors declare that no competing interests exist.
Koutarou D Kimura

Department of Biological Sciences, Graduate School of Science, Osaka University, Toyonaka, Japan

Contribution
KDK, Conceptualization, Supervision, Funding acquisition, Investigation, Methodology, Writing—original draft, Writing—review and editing

For correspondence
kokimura-lab@umin.ac.jp

Competing interests
The authors declare that no competing interests exist.

"This ORCID iD identifies the author of this article:" 0000-0002-3359-1578

Funding

Japan Society for the Promotion of Science (Grant-in-Aid for JSPS fellows)

Yuki Tanimoto
Akiko Yamazoe-Umemoto

Japan Society for the Promotion of Science (Interdisciplinary graduate school program for systematic understanding of health and disease)

Yuki Tanimoto
Shuhei J Yamazaki

Japan Agency for Medical Research and Development (Brain Mapping by Integrated Neurotechnologies for Disease Studies (Brain/MINDS))

Keiko Gengyo-Ando
Junichi Nakai

Ministry of Education, Culture, Sports, Science and Technology (Regional Innovation Cluster Program (City Area Type, Central Saitama Area))

Keiko Gengyo-Ando
Junichi Nakai

Ministry of Education, Culture, Sports, Science and Technology (KAKENHI JP23115704)

Keiko Gengyo-Ando

Ministry of Education, Culture, Sports, Science and Technology (KAKENHI JP21115504)

Junichi Nakai

Ministry of Education, Culture, Sports, Science and Technology (KAKENHI JP16H06536)

Junichi Nakai
Koichi Hashimoto

Ministry of Education, Culture, Sports, Science and Technology (KAKENHI JP23115703)

Koichi Hashimoto

Ministry of Education, Culture, Sports, Science and Technology (KAKENHI JP21115502)

Koichi Hashimoto

Ministry of Education, Culture, Sports, Science and Technology (KAKENHI JP23115711)

Koutarou D Kimura

Ministry of Education, Culture, Sports, Science and Technology (KAKENHI JP16H06545)

Koutarou D Kimura

Ministry of Education, Culture, Sports, Science, and Technology (The Osaka University Life Science Young Independent Researcher Support Program)

Koutarou D Kimura

Japan Society for the Promotion of Science (PRESTO)

Koutarou D. Kimura

Mitsubishi Foundation

Koutarou D. Kimura

Shimadzu Science Foundation

Koutarou D. Kimura

Takeda Science Foundation

Koutarou D. Kimura

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank K Tanaka (FIS Inc., Japan), Y Mori (HawkVision Inc., Japan) and K Ishida (Aska Company, Japan) for generous technical support for the hardwares. We also thank Drs. Y Fujie, I Takeuchi, M Usuyama, M Hendricks, Y Zhang, A Gottschalk, P Sengupta, H Kagoshima, M Koga, S Oda, T Obara, M Chen for materials and technical support, C Bargmann, S Nakanishi, Y Komura, T Kikuchi, M S Kitazawa, O Hobert, K Fujimoto and the Kimura laboratory members for suggestions and comments. Nematode strains were provided by the Caenorhabditis Genetics Center (funded by the NIH Office of Research Infrastructure Programs P40 OD010440), and neuronal information was provided by WormBase (funded by National Human Genome Research Institute grant #41 HG002223) and by WormAtlas (http://www.wormatlas.org).

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.