- Research article
- Open Access
Navigational choice between reversal and curve during acidic pH avoidance behavior in Caenorhabditis elegans
BMC Neuroscience volume 16, Article number: 79 (2015)
Under experimental conditions, virtually all behaviors of Caenorhabditis elegans are achieved by combinations of simple locomotion, including forward, reversal movement, turning by deep body bending, and gradual shallow turning. To study how worms regulate these locomotion in response to sensory information, acidic pH avoidance behavior was analyzed by using worm tracking system.
In the acidic pH avoidance, we characterized two types of behavioral maneuvers that have similar behavioral sequences in chemotaxis and thermotaxis. A stereotypic reversal-turn-forward sequence of reversal avoidance caused an abrupt random reorientation, and a shallow gradual turn in curve avoidance caused non-random reorientation in a less acidic direction to avoid the acidic pH. Our results suggest that these two maneuvers were each triggered by a distinct threshold pH. A simulation study using the two-distinct-threshold model reproduced the avoidance behavior of the real worm, supporting the presence of the threshold. Threshold pH for both reversal and curve avoidance was altered in mutants with reduced or enhanced glutamatergic signaling from acid-sensing neurons.
C. elegans employ two behavioral maneuvers, reversal (klinokinesis) and curve (klinotaxis) to avoid acidic pH. Unlike the chemotaxis in C. elegans, reversal and curve avoidances were triggered by absolute pH rather than temporal derivative of stimulus concentration in this behavior. The pH threshold is different between reversal and curve avoidance. Mutant studies suggested that the difference results from a differential amount of glutamate released from ASH and ASK chemosensory neurons.
Neuronal circuits transform sensory inputs into behavior in animals. Elucidating the mechanism of the sensory-motor transformation pathway is one of important issues in neurobiology. In Caenorhabditis elegans, practically all behavior under experimental conditions are achieved through sets of simple locomotion, including forward and backward movement, turning by deep body bending, shallow gradual reorientations, and resting [1, 2]. Worms adaptively combine these types of locomotion according to sensory cues. Therefore, the initial step in understanding the mechanism of the sensory-motor transformation pathway in C. elegans is to determine how the worms regulate their locomotion in response to environmental stimuli.
C. elegans shows chemotaxis and thermotaxis toward favorable chemicals and temperature, respectively [3, 4], both of which are achieved through three behavioral strategies: (1) biased random walk (i.e., klinokinesis), a repeated, abrupt, random reorientation until the worm is headed in a favorable direction [5–7], (2) non-random steering (i.e., klinotaxis), a gradual reorientation in a favorable direction [8, 9], and (3) slowing response (i.e., orthotaxis), deceleration allowing the worm to remain in favorable conditions . By the analyses of locomotion of the worms, in which specific neurons or genes were inactivated, combined with defined stimulus application in temporal and spatial manner, function of sensory neurons and the downstream neuronal circuit involved in attraction behaviors were studied in detail [7, 9, 11–16]. However, analysis of locomotion during avoidance behavior to a spatial gradient of repellents is relatively limited to date [17–19].
In this study, locomotion of C. elegans during acidic pH avoidance behavior was analyzed. Similar to the attraction behavior, behavioral maneuvers probably corresponding to biased random walk and non-random steering were observed in the avoidance. Apparently, these behavioral maneuvers were triggered by the threshold pH for each maneuver, rather than by temporal variation in the repellent concentration. The apparent threshold pH was regulated, at least in part, by glutamatergic neurotransmission by ASH and ASK neurons.
Behavioral maneuvers observed during acidic pH avoidance
To explore the locomotory behavior of C. elegans on a spatial gradient of repellent, we used acidic pH as a repulsive sensory stimulus. C. elegans has been shown to avoid acidic pH less than 4.0 . Unlike other water-soluble chemoattractants and repellents, the concentration gradient of protons (i.e., acidic pH gradient) can be visualized by adding pH-sensitive dye to the assay plate (Fig. 1a). Since the acidic yellow border of BPB dye (around pH 4.0) closely corresponds to the acidic pH avoided by worms, we can record worm avoidance behavior along with the direction and position of the worm relative to the color border.
In our observation, C. elegans showed four types of behavior during acidic pH avoidance. In long reversal avoidance (Fig. 1b, Additional file 1: Movie S1), an abrupt long backward locomotion of three or more body undulations followed by reorientation by deep body bending and subsequent forward locomotion in a less acidic direction were observed when the worm came closer to the acidic color border. Gray et al.  described that the initial deep body bending (omega turn) immediately after long reversal cause a reorientation having an average angle of 137° ± 8° toward ventral side. Because of this deep body bending, the direction of the worm locomotion after long reversal avoidance was biased in a less acidic direction. In short reversal avoidance (Fig. 1c), an abrupt short backward locomotion of two or less body undulations, followed by forward locomotion, were observed near the acidic yellow border. During this behavior, directional change after the short backward locomotion was always small (less than 45°), and reorientation bias in a less acidic direction was not observed. When the worm failed to reorient itself in a less acidic direction after the short reversal avoidance, it repeated the avoidance behavior until it could escape from the acidic region, similar to the biased random walk strategy in worm chemotaxis . In the present study, we collectively referred to these two types of behavior as reversal avoidance. The behavioral parameters described below were indistinguishable between long and short reversal. In the gradual curve avoidance (Fig. 1d, Additional file 2: Movie S2), worms avoiding the acidic region kept moving forward with gradual reorientations. The deep curve avoidance (also known as omega turn, Fig. 1e), in which worms bend their body deeply in a less acidic direction, was rarely observed. We collectively referred to these two types of behavior as curve avoidance. In successful curve avoidance, the direction of forward movement after the behavioral maneuver was biased in a less acidic direction, suggesting non-random reorientation during the acidic pH avoidance behavior.
Navigational choice during acidic pH avoidance behavior is dependent on the bearing
In this study, we analyzed two parameters: (1) the angle between the direction of forward locomotion prior to reorientation and the acidic color border (angle of encounter), and (2) the position of the worm relative to the acidic color border when the worm entered furthest into (or most closely approached) the acidic yellow region during the reorientation maneuver (Fig. 1b–e). Figure 2a shows the relationship between the angle of encounter and the two reorientation maneuvers during the avoidance behavior. Behavioral choices between reversal and curve avoidance appeared to have been dependent on the angle of encounter (correlation coefficient: R = 0.92, p < 0.01) (Fig. 2a, b). In other words, the worms that encountered the acidic region at a deep, perpendicular angle approximating 90° chose reversal more often than curve. Conversely the worms that encountered the acidic region at a shallow angle chose curve more often than reversal. The ratio of reversal avoidance gradually increased and that of curve avoidance decreased with an increase in the angle of encounter (Fig. 2b). To see how the worms accomplish angle-dependent choice of reorientation maneuvers, we examined the avoidance behavior under various conditions, as described below.
Effect of steepness of acidic pH gradient on navigational choice
Based on the previous observations on C. elegans chemotaxis behavior, it is widely accepted that the worm can sense temporal changes in attractant concentrations (dC/dt) caused by its own locomotion. Positive dC/dt suppresses reversal, and negative dC/dt promotes reversal [5, 15, 16]. In the acidic pH avoidance behavior, it is possible to think that the dC/dt of repellent is also involved in the angle-dependent navigational choice. When the worms encountered the acidic region at a deep angle, they sensed a larger dC/dt than when they encountered the acidic region at a shallow angle. This hypothesis implies that a large dC/dt causes reversal and a small dC/dt causes curve avoidance. In this scenario, the steepness of the acidic pH gradient should affect navigational choice. In a shallow gradient, dC/dt is consistently lower than in a steep gradient, thus the worm may chose more curve than in a steep gradient. In order to see the effect of the steepness of the gradient, we changed the buffer and made an acidic pH gradient using a lower concentration of HCl.
The gradient became less steep when we used 100 mM, 30 mM HCl on 5 mM sodium acetate buffer (pH 6.0) and 30 mM HCl on 0.2 mM sodium acetate buffer (Fig. 3a). Under these conditions, however, the worms showed essentially the same pattern of angle dependence as in the steep gradient (Figs. 2, 3b, c). The ratio of curve avoidance did not increase even in the shallowest gradient examined. Therefore, dC/dt of the repellent (i.e., steepness of the gradient) may not be a major determinant of the angle-dependent navigational choice in this avoidance behavior.
Two-distinct-threshold model on navigational choice during acidic pH avoidance
In the course of analyzing worm position relative to the color border, we noted that under every condition examined, the mean distance at which the reversal avoidance occurred locates more acidic region than that at which curve avoidance occurred (Fig. 4). Although the difference between the mean distance for reversal and curve avoidance increased in shallower gradients, the pH values at the mean distance estimated from the measured pH gradient (Fig. 3a) were similar, especially for the curve avoidance. This observation led us to hypothesize that there are two distinct threshold pH values for reversal and curve. In this two-threshold model, the worm begins curving when it approaches the less acidic threshold pH for curve. If worm encounters the acidic region at a shallow angle, it can complete the curve avoidance with a single curve or multiple sequential curves without reaching the second, more acidic threshold pH for reversal. However, if the worm encounters the acidic region at a deep angle, it reaches the second threshold for reversal before it finishes its curve, thus executing reversal avoidance. Indeed, small reorientations like the beginning of curve avoidance were frequently observed immediately before reversal avoidance (Additional file 1: Movie S1). In addition, the worms showed small head swings (Additional file 3: Fig. S3), which were different from body undulation during the avoidance maneuvers. The small head swings may reflect the worm’s repeated sampling of the local pH to improve the fidelity of avoidance behavior, supporting the idea that the worm detects the local threshold pH rather than temporal variation of the stimuli caused by its locomotion. Taken together, these results suggest the presence of two distinct threshold pH values for the two avoidance maneuvers. These results also indicate that the apparent angle-dependent choice during acidic pH avoidance behavior is the result of the difference between these two threshold pH values.
The mean distance of the reversal and curve in Fig. 4 roughly corresponds to pH 3.8–4.0 and pH 4.6–4.8, respectively, when we estimated from the gradient measured in Fig. 3a. Each threshold pH may be found near or within the respective ranges. However, we could not estimate the precise threshold pH value because of the low spatial resolution of our pH measurements. Although the steepness of the repellent gradient did not affect the navigational choice, steepness of gradient did have marked effect on the variability of worm position at which avoidance was executed. Both curve and reversal avoidance occurred in a very narrow range of distance when a steep gradient was used, as opposed to that when shallow gradient was used (Fig. 4; Additional file 3: Fig. S4).
Computer simulation of acidic pH avoidance using the two-distinct-threshold model
To evaluate the validity of the two-distinct-threshold model for navigational choice in acidic pH avoidance behavior in C. elegans, we conducted a computer simulation study by constructing a stochastic model for the behavior. As described above, we hypothesized that the decision between curve and reversal was driven by the absolute pH value detected by the worm. Therefore, in this model, we defined the probability of reorientation by two stochastic functions for curve and reversal. These functions were pH-dependent and each had a distinct threshold pH. Behavioral parameters of the model worm, such as speed  and degree of reorientation after reversal  were defined based on observations of the real worm. Parameters used in these functions were explored to reproduce the experimental data of real worms under the condition shown in Fig. 4d. Accordingly, we developed a model that had an apparent angle-dependent choice of behavioral maneuvers. The model was run on different conditions having acidic pH gradients steeper than a gradient used for model construction. The model could reproduce results similar to those observed in the real worm under steeper conditions (Fig. 5b, c; Additional file 3: Fig. S5), demonstrating the wide applicability of this model. Although the mean distances did not exactly replicate those of real worms, the angle distribution and shape of histogram resembled those of the real worm. These results strongly support the idea that worms choose a behavioral maneuver according to the distinct threshold pH value during avoidance behavior.
Roles of glutamatergic neurotransmission on navigational choice
At least four classes of bilaterally symmetric chemosensory neuron pairs, ASE, ADF, ASH, and ASK are involved in acidic pH avoidance behavior in C. elegans . Of these, ASH and ASK neurons are known to be glutamatergic because eat-4, a neuronal vesicular glutamate transporter gene necessary for glutamatergic neurotransmission, expressed in these neurons . To study the role of glutamatergic signaling in the acidic pH avoidance behavior, the eat-4(ky5) mutant was used for behavioral analysis.
Consistent with the notion that the eat-4(ky5) mutant has a defect in spontaneous backward locomotion upon food removal , the frequency of reversal avoidance was obviously reduced in the mutant (Fig. 6a, b, e, f). Moreover, the mean distance at which the reversal and curve avoidance occurred was significantly shifted to the acidic region (Fig. 6f, i). These results indicate that glutamatergic signaling in the C. elegans nervous system affects several aspects of acidic pH avoidance behavior including reversal frequency and worm position relative to the acidic color border, which reflects putative threshold pH for avoidance. Both reversal frequency and mean distance were partially recovered by ASH- and ASK-specific expression of wild-type eat-4 cDNA (Fig. 6b–d, g–i) driven by sra-6 and sra-9 promoter, respectively . Since many other glutamatergic neurons, such as interneurons, might be involved in the avoidance behavior, the sensory neuron specific expression was not sufficient to recover the phenotypic defects completely.
We next examined the acidic pH avoidance behavior of rgs-3(vs19), egl-3(n150) and egl-3(n729), to investigate the role of intensity of glutamatergic signaling from chemosensory neurons. The rgs-3 gene encodes a C. elegans homolog of a mammalian RGS (regulator of G-protein signaling) proteins, and is expressed in a small subset of sensory neurons including ASH and ASK . Loss of the rgs-3 gene products suggested to cause a reduced glutamatergic output from the ASH synapse. In contrast, the egl-3 mutant has enhanced glutamatergic output because the egl-3 gene encodes an enzyme required for processing of the neuropeptide that downregulates glutamatergic output from ASH [26, 27].
The apparent angle dependence of the navigational choice was not affected in rgs-3(vs19), egl-3(n150), and egl-3(n729) mutants (data not shown). However, the positions at which avoidance was executed were significantly different. In the rgs-3(vs19) mutant having reduced glutamatergic signaling, both reversal and curve avoidance occurred at more acidic regions than those in wild-type. Whereas in the egl-3 mutants having enhanced glutamatergic signaling, both reversal and curve avoidance occurred at less acidic regions than wild-type (Fig. 7a). To further confirm the phenotype, the egl-3 mutants were examined using slightly shallower gradient, and showed statistically significant differences from wild-type (Fig. 7b). Similar to the case for eat-4 mutant, ASH- and ASK-specific expression of wild-type rgs-3 recovered the phenotype (Fig. 7c). These results indicate that the glutamatergic signal output from ASH and ASK is involved in the regulation of both reversal and curve avoidance. Although the two threshold pHs are distinct, the neurotransmitter regulating these two behavioral maneuvers is not distinct. Large amounts of glutamate may be released from ASH and ASK in the presence of an intense stimulus, and small amounts of the transmitter may be released upon weak stimulation. Information about the ambient pH may be translated into the amount of glutamate synaptically released. In the rgs-3 mutant having reduced glutamatergic output, ASH and ASK neurons require more intense stimulation to release the threshold amount of glutamate. A weaker stimulus is sufficient for ASH and ASK neurons in the egl-3 mutant having enhanced glutamatergic output to reach the threshold.
C. elegans appears to have two distinct threshold pHs for curve and reversal during acidic pH avoidance behavior
The locomotory sequences observed during acidic pH avoidance were quite similar to those observed in attraction behavior [5–10], behavior upon food removal [21, 23, 28, 29] and spontaneous locomotion [1, 2], suggesting that the avoidance behavior consists of intrinsic locomotory programs similar to other types of behavior in C. elegans.
In the positive chemotaxis and thermotaxis, behavioral strategies such as biased random walk and non-random steering are largely regulated by temporal variations in attractant concentration and ambient temperature caused by the worm locomotion [5, 6, 8, 13, 15, 16, 30]. Here, we propose that in the acidic pH avoidance behavior, both reversal (a form of biased random walk) and curve avoidance (a form of non-random steering) are regulated by a respective threshold pH. Wild-type worms initially executed curve avoidance at the relatively less acidic threshold pH for curve. If reorientation by curve avoidance was not sufficient and the worms reached the more acidic second threshold pH for reversal before they finished curving, they executed reversal (Fig. 8a). Although the behavioral maneuvers that occurred during the avoidance behavior were not completely sequential, and the worm positions at which they avoided the acidic pH were variable possibly because of the stochastic nature of behavioral decision making within the nervous system of the worm, the threshold model is plausible because in a population assay of the acidic pH avoidance behavior on a radial spatial gradient, the worm traces draw a definitive clear zone implying the threshold . Furthermore, our computer simulation study substantiates the effectiveness of the two-threshold method during acidic pH avoidance behavior.
Although temporal variation in acidity seems to have no effect on the navigational choice of behavioral maneuver, temporal variation may play a role during curve avoidance. Forward locomotion of the worms after curve avoidance is biased in a less acidic direction. Worms may accomplish this non-random steering by detecting the variation in acidity with head swings, similar to those observed in isothermal tracking behavior and salt chemotaxis [8, 31]. Temporal variations may also be important when the worm is in the middle of an acidic pH gradient lower than pH 4.0. The threshold method is only effective when close to the threshold range.
The threshold method may be shared in avoidance behavior that has obvious thresholds in population assays, such as heavy metals  and osmotic avoidance . This method, however, may not be common in all avoidance behavior. Iino and Yoshida  described that during sodium chloride avoidance after salt chemotaxis learning [34–36], the pirouette and weathervane responses are reversed. Worms move down the salt gradient as far as possible according to the temporal variation in the sodium chloride concentration.
In an attraction behavior, information about decrements in attractant concentration is as valuable as that about increments in the concentration. Whereas in acidic pH avoidance behavior, information about decreased acidity is not important for worms. At a pH 5.0–6.0, worms show no behavioral responses to ambient pH . This disproportionate importance of information may play a role in the evolution of differences in behavioral regulation (threshold vs. temporal variation).
Decoding of the threshold pH within a neuronal circuit of the worm
Our mutant analysis suggested that both reversal and curve avoidance are regulated, at least in part, by glutamatergic synaptic output from ASH and ASK neurons. As the intensity of the stimulus increases, the amount of glutamate released from these neurons to synaptic clefts may also increase. Neurons postsynaptic to these neurons may decode the information about the intensity of signals by referring an intrinsically encoded threshold amount of glutamate. Backward locomotion is promoted during reversal avoidance, while forward locomotion is promoted during curve avoidance. How does a single neurotransmitter released from a small number of sensory neurons regulate two mutually exclusive behavioral responses? An important insight was suggested by Mellem et al. . Based on their integrative study using genetics and electrophysiology, they proposed a “differential activation model” in which a specific level of glutamate released from ASH causes a distinct response in the postsynaptic neurons. In their model, a small amount of glutamate activates GLRs (and possibly GLC) synaptically-localized and a large amount of glutamate activates perisynaptic NMR receptors in addition to GLR/GLC at the synapse between ASH and backward command interneurons. According to the wiring diagram of the worm neuronal network drawn by White et al. , a bilateral pair of ASH neurons sends their synaptic outputs to the backward command interneurons AVA, AVD, and AVE as well as to the forward command interneuron AVB. In acidic pH avoidance behavior, a small amount of glutamate released by ASH may activate synaptic GLR receptors on both backward and forward command interneurons. In addition, GLC receptors may also be activated in ASH-backward interneuron synapses, thus downregulating neuronal excitation. Overall, forward locomotion is promoted by weaker stimulation during curve avoidance. In this situation, amphid interneurons downstream to ASH and/or ASK may be activated for regulation of non-random steering. In the presence of an intense stimulus, backward command interneurons may be more activated through perisynaptic NMR receptors by glutamate spilled out from the synapse, thereby promoting backward locomotion (Fig. 8b, Additional file 3: Fig. S6). Although the role of ASK is not clear, ASK may contribute to the regulation of reversal avoidance via pathway yet unidentified, because ASK has an important role on promoting backward locomotion .
Cell ablation studies on both sensory and interneurons, analysis of mutants defective in the postsynaptic glutamate receptors, physiological analysis either by electrophysiological or optical methods may help us understanding fully the mechanism of behavioral regulation during the acidic pH avoidance behavior. In this context, it should be noted that the acid-sensing chemosensory neurons ASE and ADF send their synaptic outputs mainly to AIY and AIZ interneurons, respectively, which may be involved in the regulation of non-random steering .
C. elegans avoids acidic pH using two behavioral maneuvers, reversal (klinokinesis) and curve (klinotaxis). It is widely accepted that in positive chemotaxis, C. elegans detects temporal derivative of stimulus concentration (dC/dt) to regulate reversal locomotion. However, in acidic pH avoidance, probability of reversal and curve avoidances were dependent on absolute pH rather than dC/dt. The pH threshold is different between reversal and curve avoidance. Curve avoidance always occurred at less acidic region than reversal, irrespective of steepness of the acidic pH gradient. Mutant studies suggested that the difference in threshold pH results from a differential amount of glutamate released from ASH and ASK chemosensory neurons in response to stimulus intensity.
Strains and culture
The wild-type animal used in this study was Caenorhabditis elegans ver. Bristol, strain N2. Worms were grown under standard conditions at 20 °C . Mutant strains used were CX5 eat-4(ky5) III, LX242 rgs-3(vs19) II, MT1541 egl-3(n729) V, MT150 egl-3(n150) V, obtained from Caenorhabditis Genetics Center (CGC). The transgenic strains used were WKB70 rgs-3(vs19) II; iwtEx50[sra-6::rgs-3cDNA; myo-3::RFP], WKB71 rgs-3(vs19) II; iwtEx51[sra-6::rgs-3cDNA; myo-3::RFP], WKB72 rgs-3(vs19) II; iwtEx52[sra-9::rgs-3cDNA; unc-122::mCherry], WKB73 rgs-3(vs19) II; iwtEx53[sra-9::rgs-3cDNA; unc-122::mCherry],WKB77 eat-4(ky5)III; iwtEx57[sra-9::eat-4cDNA, mec-4::GFP], WKB78 eat-4(ky5)III; iwtEx58[sra-9::eat-4cDNA, mec-4::GFP], WKB79 eat-4(ky5)III; iwtEx59[sra-6::eat-4cDNA, mec-4::GFP], WKB80 eat-4(ky5)III; iwtEx60[sra-6::eat-4cDNA, mec-4::GFP]. We used sra-6 and sra-9 promoters for cell-specific expression of cDNA in ASH and ASK neurons, respectively . Germline transformation of mutants were performed by microinjecting a mixture of rescue construct (20 ng/μl) and marker construct (80 ng/μl) into hermaphrodite gonad according to Mello et al. .
Analysis of acidic pH avoidance behavior
A 9-cm petri dish containing 10 ml of 2 % agar, 5 or 0.2 mM sodium acetate pH 6.0 (pH adjusted with acetic acid), 50 mM NaCl, 1 mM CaCl2, and 1 mM MgSO4 was used as an assay plate. In some experiments, 5 mM potassium phosphate (pH 6.0) was also used instead of sodium acetate. To visualize acidic regions on the plate, 1 ml of water-saturated bromophenol blue (BPB) solution was included in the gel. To form the acidic pH gradient, 2 % acidic agar gels containing 100 mM or 30 mM HCl were prepared by mixing equal volumes of 4 % agar solution and 200 mM or 60 mM HCl solution maintained at 60 °C, then pouring into a rectangular plastic dish (W 8 cm × L 12 cm). Three to four strips of the acidic agar gel (W 5 mm × L 45 mm × H 1 mm) were excised from the dish, arranged in a triangle or rectangle on the assay plate. The assay plate with acidic agar gel strips was left for 30 min at 20 °C until an acidic pH gradient formed (Fig. 1a). Although the acidic yellow region was expanded over time, relationship between yellow color border and worm avoidance behaviors (distance and angle, see below) were not changed with the time.
A single well-fed young adult C. elegans was transferred to the center of the assay plate, and allowed to move freely on the plate for 50 min. The behavior of the worm was tracked using automatic tracking system  and was recorded. Movies were analyzed using custom-made software. Acidic pH avoidance behavior by reversal or curve were identified by watching the movies and were recorded along with the time, angle of encounter to the acidic yellow border, and the distance of the worm position relative to the yellow border when the worm moved furthest into the acidic region during the avoidance behavior (Fig. 1b–e). Reversal locomotion of three or more head swings were referred to long and that of two or less head swings were referred to short reversals, deep body bendings (omega-shaped turns) were visually identified by the head nearly touching the tail as described by Gray et al. . Gradual curves were also visually identified when worms change their direction less than 90° within single bend as described by Kim et al. . Although the wild-type C. elegans consistently avoided the acidic region, it occasionally failed in the process of achieving eventual avoidance. In this situation, the worm showed multiple consecutive avoidance behavior, until it succeeded in avoiding the acidic region. We only analyzed the first occurrence of the serial avoidance to precisely evaluate the relationship between the initial choice of behavioral maneuvers and the two parameters.
Six thin (1 mm in width) strips of agar gel around the acidic yellow border were excised from the assay plates (−3 to +3 mm, as indicated in Fig. 1f) and collected individually at the bottom of separate 1.5 ml sample tubes. The pH of each of the gel strips was measured by a micro pH electrode (9669-10D, HORIBA, Tokyo, Japan).
In this model, the location of the worms were represented as a point (x t , y t ) almost corresponding to their center of mass (Fig. 5a). The model worms moved on a virtual field that had an acidic pH gradient, same as that shown in Fig. 3a. The model worms moved in a straight line in one direction at a constant speed (0.15 mm/s), updating their position and direction every 0.8 s. The model worms had a pH sensor 0.51 mm ahead of its positional point and swung it side by side, as shown in Fig. 5a. The sensor detected the ambient pH every 0.8 s. The initial position of the model worm was randomly determined within the less acidic region of the virtual field (Additional file 3: Fig. S5A). To represent the frequency of the angle of encounter in the real worm, the initial direction θ 0 (−90° < θ 0 ≤ 90°) of the worm movement was selected from random values having 1 − |sinθ| distribution.
The model worm implemented two types of reorientation mechanisms corresponding to curve and reversal. Both reorientation events occurred stochastically, depending on the pH value detected. The probability of curve was defined by the following equation:
Similarly, the probability of reversal was defined by the equation:
where C is a pH value detected by the model worm, α c and β c correspond to deviation and center of probability density function for curve, respectively, and α r and β r correspond to those for reversal, respectively. These parameters conferred threshold pH and variation for each reorientation on the model worm. The model worm reoriented in a new direction when the curve or reversal was determined. If the curve and reversal were determined at the same time point, reversal was selected. The degree of reorientation in reversal was chosen from a normally distributed random value having 137° ± 8° distributions . The direction of reorientation in reversal was determined randomly toward either side. The magnitude of reorientation in curve was chosen from a normally distributed random value having φ c ± Δφ c distribution, in which φ c and Δφ c were determined by parameter optimization (see below). The direction of reorientation in curve, determined at time t n , was selected by comparing the pH detected at t n and t n−1 . The model worm reoriented toward the less acidic side. The behavior of the model worms started at time t 0 and ended when the worms reoriented in a less acidic direction by either of the two reorientation maneuvers. The type of reorientation maneuver selected at the end of the behavior was recorded along with the position of the model worm and direction of movement immediately before the reorientation.
Parameter sets (α c , β c , α r , β r in pH values) in Eqs. (1) and (2) and degrees of curve reorientations (φ c ± Δφ c in angle) were selected from several values. Parameter values used were α c (0.125, 0.25, 0.5), β c (2.8, 3.3, 3.8, 4.3), α r (0.125, 0.25, 0.5), β r (β c − 0.4, β c − 0.8, β c − 01.2, β c − 1.6), φ c (30, 60, 80), Δφ c (10, 30, 60), and initial position of the model worm (<−3.0, <−2.0, <−1.0 in mm). Simulations were performed on 3000 model worms for every possible combination of parameter sets (2916 combinations). A set of parameters that produced results matching the experimental data of the real worm shown in Fig. 4d were selected (Fig. 5d). The same parameter set was used for analysis under different conditions with steeper pH gradients (Fig. 5b, c; Additional file 3: Fig S5A–C).
All programs were written from scratch by using C and Ruby.
Availability of supporting data
The source code for computer simulation study is available from https://github.com/kzmsakata/threshold_turn.
Croll NA. Components and patterns in the behaviour of the nematode Caenorhabditis elegans. J Zool. 1975;176:159–76.
Kim D, Park S, Mahadevan L, Shin JH. The shallow turn of a worm. J Exp Biol. 2011;214:1554–9.
Ward S. Chemotaxis by the nematode Caenorhabditis elegans: identification of attractants and analysis of the response by use of mutants. Proc Natl Acad Sci USA. 1973;70:817–21.
Hedgecock EM, Russell RL. Normal and mutant thermotaxis in the nematode Caenorhabditis elegans. Proc Natl Acad Sci USA. 1975;72:4061–5.
Pierce-Shimomura JT, Morse TM, Lockery SR. The fundamental role of pirouettes in Caenorhabditis elegans chemotaxis. J Neurosci. 1999;19:9557–69.
Ryu WS, Samuel ADT. Thermotaxis in Caenorhabditis elegans analyzed by measuring responses to defined Thermal stimuli. J Neurosci. 2002;22:5727–33.
Chalasani SH, Chronis N, Tsunozaki M, Gray JM, Ramot D, Goodman MB, Bargmann CI. Dissecting a circuit for olfactory behavior in Caenorhabditis elegans. Nature. 2007;450:63–70.
Luo L, Clark DA, Biron D, Mahadevan L, Samuel ADT. Sensorimotor control during isothermal tracking in Caenorhabditis elegans. J Exp Biol. 2006;209:4652–62.
Iino Y, Yoshida K. Parallel use of two behavioral mechanisms for chemotaxis in Caenorhabditis elegans. J Neurosci. 2009;29:5370–80.
Albrecht DR, Bargmann CI. High-content behavioral analysis of Caenorhabditis elegans in precise spatiotemporal chemical environments. Nat Methods. 2011;8:599–605.
Zariwala HA, Miller AC, Faumont S, Lockery SR. Step response analysis of thermotaxis in Caenorhabditis elegans. J Neurosci. 2003;23:4369–77.
Chung SH, Clark DA, Gabel CV, Mazur E, Samuel ADT. The role of the AFD neuron in C. elegans thermotaxis analyzed using femtosecond laser ablation. BMC Neurosci. 2006;11:1–11.
Clark DA, Gabel CV, Gabel H, Samuel ADT. Temporal activity patterns in thermosensory neurons of freely moving Caenorhabditis elegans encode spatial thermal gradients. J Neurosci. 2007;27:6083–90.
Gabel CV, Gabel H, Pavlichin D, Kao A, Clark D, Samuel ADT. Neural circuits mediate electrosensory behavior in Caenorhabditis elegans. J Neurosci. 2007;27:7586–96.
Suzuki H, Thiele TR, Faumont S, Ezcurra M, Lockery SR, Schafer WR. Functional asymmetry in Caenorhabditis elegans taste neurons and its computational role in chemotaxis. Nature. 2008;454:114–7.
Thiele TR, Faumont S, Lockery SR. The neural network for chemotaxis to tastants in Caenorhabditis elegans is specialized for temporal differentiation. J Neurosci. 2009;29:11904–11.
Bretscher AJ, Busch KE, de Bono M. A carbon dioxide avoidance behavior is integrated with responses to ambient oxygen and food in Caenorhabditis elegans. Proc Natl Acad Sci USA. 2008;105:8044–9.
Kimura KD, Fujita K, Katsura I. Enhancement of odor avoidance regulated by dopamine signaling in Caenorhabditis elegans. J Neurosci. 2010;30:16365–75.
Bretscher AJ, Kodama-Namba E, Busch KE, Murphy RJ, Soltesz Z, Laurent P, de Bono M. Temperature, oxygen, and salt-sensing neurons in C. elegans are carbon dioxide sensors that control avoidance behavior. Neuron 2011; 69:1099–1113.
Sambongi Y, Takeda K, Wakabayashi T, Ueda I, Wada Y, Futai M. Caenorhabditis elegans senses protons through amphid chemosensory neurons: proton signals elicit avoidance behavior. NeuroReport. 2000;11:2229–32.
Gray JM, Hill JJ, Bargmann CI. A circuit for navigation in Caenorhabditis elegans. Proc Natl Acad Sci USA. 2005;102:3184–91.
Lee RY, Sawin ER, Chalfie M, Horvitz HR, Avery L. EAT-4, a homolog of a mammalian sodium-dependent inorganic phosphate cotransporter, is necessary for glutamatergic neurotransmission in Caenorhabditis elegans. J Neurosci. 1999;19:159–67.
Hills T. Dopamine and glutamate control area-restricted search behavior in Caenorhabditis elegans. J Neurosci. 2004;24:1217–25.
Troemel ER, Chou JH, Dwyer ND, Colbert Ha, Bargmann CI. Divergent seven transmembrane receptors are candidate chemosensory receptors in C. elegans. Cell 1995; 83:207–218.
Ferkey DM, Hyde R, Haspel G, Dionne HM, Hess HA, Suzuki H, Schafer WR, Koelle MR, Hart AC. C. elegans G protein regulator RGS-3 controls sensitivity to sensory stimuli. Neuron. 2007;53:39–52.
Kass J, Jacob TC, Kim P, Kaplan JM. The EGL-3 proprotein convertase regulates mechanosensory responses of Caenorhabditis elegans. J Neurosci. 2001;21:9265–72.
Mellem JE, Brockie PJ, Zheng Y, Madsen DM, Maricq AV. Decoding of polymodal sensory stimuli by postsynaptic glutamate receptors in C. elegans. Neuron. 2002;36:933–44.
Tsalik EL, Hobert O. Functional mapping of neurons that control locomotory behavior in Caenorhabditis elegans. J Neurobiol. 2003;56:178–97.
Wakabayashi T, Kitagawa I, Shingai R. Neurons regulating the duration of forward locomotion in Caenorhabditis elegans. Neurosci Res. 2004;50:103–11.
Garrity P, Goodman MB, Samuel ADT, Sengupta P. Running hot and cold: behavioral strategies, neural circuits, and the molecular machinery for thermotaxis in C. elegans and Drosophila. Genes Dev. 2010;24:2365–82.
Satoh Y, Sato H, Kunitomo H, Fei X, Hashimoto K, Iino Y. Regulation of experience-dependent bidirectional chemotaxis by a neural circuit switch in Caenorhabditis elegans. J Neurosci. 2014;34:15631–7.
Sambongi Y, Nagae T, Liu Y, Yoshimizu T, Takeda K, Wada Y, Futai M. Sensing of cadmium and copper ions by externally exposed ADL, ASE, and ASH neurons elicits avoidance response in Caenorhabditis elegans. NeuroReport. 1999;10:753–7.
Culotti JG, Russell RL. Osmotic avoidance defective mutants of the nematode Caenorhabditis elegans. Genetics. 1978;90:243–56.
Saeki S, Yamamoto M, Iino Y. Plasticity of chemotaxis revealed by paired presentation of a chemoattractant and starvation in the nematode Caenorhabditis elegans. J Exp Biol. 2001;1764:1757–64.
Hukema RK, Rademakers S, Dekkers MPJ, Burghoorn J, Jansen G. Antagonistic sensory cues generate gustatory plasticity in Caenorhabditis elegans. EMBO J. 2006;25:312–22.
Tomioka M, Adachi T, Suzuki H, Kunitomo H, Schafer WR, Iino Y. The insulin/PI 3-kinase pathway regulates salt chemotaxis learning in Caenorhabditis elegans. Neuron. 2006;51:613–25.
White JG, Southgate E, Thomson JN, Brenner S. The structure of the nervous system of the nematode Caenorhabditis elegans: the mind of a worm. Phil Trans R Soc L. 1986;314:1–340.
Brenner S. The genetics of Caenorhabditis elegans. Genetics. 1974;77:71–94.
Mello CC, Kramer JM, Stinchcomb D, Ambros V. Efficient gene transfer in C. elegans: extrachromosomal maintenance and integration of transforming sequences. EMBO J. 1991;10:3959–70.
Hoshi K, Shingai R. Computer-driven automatic identification of locomotion states in Caenorhabditis elegans. J Neurosci Methods. 2006;157:355–63.
TW participated in the design of the study, carried out the molecular genetic experiments, behavioral assay, data analysis, and wrote the manuscript. KS carried out the simulation study. TT, HI, SS, MW participated in the behavioral assay and data analysis. RS supervised the project. All authors read and approved the final manuscript.
We thank Caenohabditis Genetics Center, which is funded by the National Institutes of Health-National Center for Research Resources, for providing us the worm strains. We also thank Mr. Ryo Yamatsuki for helping with the software for behavioral analysis. This work was supported by MEXT KAKENHI (Innovative Areas: Systems Molecular Ethology ).) to R.S. and by a Research Program for Olfactory Communication in Mammals from Ministry of Education, Culture, Sports, Science, and Technology in Japan to T.W.
The authors declare that they have no competing interests.
Tokumitsu Wakabayashi and Kazumi Sakata contributed equally to this work
A movie for long reversal avoidance of wild-type C. elegans. The bright region on the right shows the acidic yellow region and the dark region shows the less acidic region.
A movie for gradual curve avoidance of wild-type C. elegans. The bright region on the left shows acidic yellow region and the dark region shows the less acidic region.
Supplemental Figures S3–S6.
About this article
Cite this article
Wakabayashi, T., Sakata, K., Togashi, T. et al. Navigational choice between reversal and curve during acidic pH avoidance behavior in Caenorhabditis elegans . BMC Neurosci 16, 79 (2015). https://doi.org/10.1186/s12868-015-0220-0
- Caenorhabditis elegans
- Acidic pH avoidance
- Sensory-motor transformation