The nematode worm C. elegans chooses between bacterial foods as if maximizing economic utility

  1. Abraham Katzen
  2. Hui-Kuan Chung
  3. William T Harbaugh
  4. Christina Della Iacono
  5. Nicholas Jackson
  6. Elizabeth E Glater
  7. Charles J Taylor
  8. Stephanie K Yu
  9. Steven W Flavell
  10. Paul W Glimcher
  11. James Andreoni
  12. Shawn R Lockery  Is a corresponding author
  1. Institute of Neuroscience, University of Oregon, United States
  2. Center for Neural Science, New York University, United States
  3. Neuroscience Institute, New York University School of Medicine, United States
  4. Department of Economics, University of Oregon, United States
  5. Department of Neuroscience, Pomona College, United States
  6. Department of Chemistry, Pomona College, United States
  7. Picower Institute for Learning and Memory, Department of Brain & Cognitive Sciences, Massachusetts Institute of Technology, United States
  8. Department of Economics, University of California, San Diego, United States
10 figures, 1 video, 3 tables and 7 additional files

Figures

Design of a GARP experiment and tests for utility maximization.

(A) Direct violation of utility maximization. Diagonal lines indicate choice sets (n=11). Choice sets are distinguished by having different values of the overall budget and/or different prices for at least one of the goods. Within a choice set, the expenditure implied by each bundle of goods is constant and equals the budget. In A (blue), the budget is $8, oranges are $2 per unit, apples are $1 per unit. In B (red): the budget is $8, oranges are $1 per unit, apples are $2 per unit. Filled circles, chosen amounts; plus signs, available amounts not chosen. Given the choices shown and the more is better rule, adbab and bcaba. Therefore, choices a and b directly violate utility maximization. (B). Indirect violation of utility maximization. Symbols as in A. The choices a and c constitute an indirect violation of utility maximization as described in the text.

Figure 2 with 2 supplements
Edible bacteria act as goods over which worms form preferences through experience.

(A) Food quality training and preference assays. Filled circles represent patches of bacteria as indicated in the key. Stars indicate worm starting locations. (B). Mean preference at t= 60 min. in the open-field accumulation assay for trained and untrained N2 and ceh-36 mutants. Asterisk, see Table 21. Replicates, Trained, strain(N): N2(8), ceh-36(9). Replicates, Untrained, strain(N): N2(8), ceh-36(8). (C). Mean preference vs. time for trained and untrained N2 worms in T-maze accumulation assays, with and without sodium azide in the food patches. (D). Mean preference index vs. time for trained and untrained cat-2(tm2261) mutants and N2 controls in T-maze accumulation assays. N2 data are from C. (B–D). Error bars, 95% CI. For sample size (N), statistical methods used, and significance level, see Table 2, rows 1-15.

Figure 2—source data 1

Edible bacteria act as goods over which worms form preferences through experience.

https://cdn.elifesciences.org/articles/69779/elife-69779-fig2-data1-v2.xlsx
Figure 2—source data 2

Loss of dopamine signaling does not reduce proportion of time on food.

https://cdn.elifesciences.org/articles/69779/elife-69779-fig2-data2-v2.xlsx
Figure 2—figure supplement 1
T-maze diagram.

The teardrop shaped feature acts as worm diode. Worms leave the teardrop more easily than they re-enter it, thus helping to ensure they do not congregate at the starting point. The mask’s ability to confine worms in the test area is improved by floating the maze on the surface of the agarose medium while it is still liquid (M. Brooks, pers. comm.). Dimensions are in mm.

Figure 2—figure supplement 2
Loss of dopamine signaling does not reduce proportion of time on food.

Mutations in the gene cat-2 do not significantly later the ability of worms to find and remain in food patches. Statistics (t-tests): N2 vs. cat-2(n4547), t(38) = –0.95, p=0.35; N2 vs cat-2(e1112), t(40) = –1.38, p=0.18. N2, n = 30; cat-2(n4547), n = 10; cat-2(e1112) n = 12. Error bars, SEM.

Single-worm food choice assays.

(A) Layout of the Y-chip. Asterisks indicate the position of recording electrodes. Ground electrodes (not shown) were inserted into the food and buffer ports to reduce electrical interference. The chip is shown configured for the experiments in Figures 4B and 5B. (B). Area of detail shown in A. The dashed line is the centerline of chip; the white line within the worm is its centerline. The black arrow connects the middle of the neck where it enters the food channel with anterior end of worm’s centerline. Positive values of head angle (θ) indicate displacement toward H food. Colored arrows show direction of flow. (C). Schematic overview of fluidic system. (D). Typical electropharyngeogram. Each pair of excitation (E) and relaxation (R) spikes constitutes one pharyngeal pump.

Validation of the Y-chip for measuring food preferences.

(A) Familiar food effect. Mean pump frequency of worms grown on H or M food and tested on the same or the converse food. Asterisk, see Table 217. Both foods were at OD 1. Pumping was recorded for 12 min. Replicates, Grown on H, tested on (N): Familiar(16), Unfamiliar(6). Replicates, Grown on M, tested on(N): Familiar(23), Unfamiliar(13). Error bars, 95% CI. (B). Food quality learning. Mean fraction of pumps in H food in trained and untrained worms. Asterisk, see Table 220. The dashed line indicates equal preference for H and M food. Both foods were at OD 1. Error bars, 95% CI. (C,D). Time course of pump frequency at four different densities of familiar food. Food enters the chip at t=0 sec. Optical density is indicated next to each trace. The black trace shows pumping in the absence of food. Shading, ± SEM E. Dependence of mean peak pump frequency density of familiar food. (F). Dependence of latency to half-maximal pump frequency on density of familiar food. (E,F). Error bars, 95% CI. For sample size (N), statistical methods used, and significance level, see Table 2, rows 16-000.

Figure 4—source data 1

Validation of the Y-chip for measuring food preferences.

https://cdn.elifesciences.org/articles/69779/elife-69779-fig4-data1-v2.xlsx
Figure 5 with 2 supplements
Economic analysis of food choice in C. elegans.

(A) Demand curves. Mean consumption of familiar food versus its price. Consumption is computed as number of pumps times optical density of bacteria. Price is computed by Equation 2. The data are fit by Equation 8 with ϵ=-2.2 for H food and ϵ=-1.9 for M food. H food price(N): 0.33(13), 1.0 (16), 3.3 (21), 10(15). M food price(N): 0.33(9), 1.0 (23), 3.3 (7), 10(10). (B). Price ratio curves. Food preference, measured as fraction of pumps in H food, versus price ratio for Trained and Untrained worms. Horizontal dashed line: indifference between H and M food; vertical dashed line: H and M food at equal price. Data at log price ratio = 0 are replotted from Figure 4B. Replicates, Trained, point(N): a(9), b(14), c(18), d(27), e(20), f(10), g(6). Replicates, Untrained, point(N): a(12), b(17), c(13), d(12), e(28), f(10), g(7). (C). GARP analysis of C. elegans food preferences. Plotted points show mean consumption of M food versus consumption of H food in Trained and Untrained worms. Lines are choice sets as in Figure 1. The x and y intercepts of each line indicate the amounts of H and M food that would have been consumed if the worm spent all its pumps on one or the other food type. Error bars, 95% CI. (D). Predicted consumption of H and M food in Trained and Untrained animals on a widely-used ensemble containing 11 budget lines (Harbaugh et al., 2001). (A–C). Error bars, 95% CI. For sample size (N), statistical methods used, and significance level, see Table 2, rows 23-26.

Figure 5—source data 1

Economic analysis of food choice in C. elegans.

https://cdn.elifesciences.org/articles/69779/elife-69779-fig5-data1-v2.xlsx
Figure 5—source data 2

Distributions of preference values in trained and untrained animals in Figure 5B.

https://cdn.elifesciences.org/articles/69779/elife-69779-fig5-data2-v2.xlsx
Figure 5—figure supplement 1
Distributions of preference values in trained and untrained animals in Figure 5B.
Figure 5—figure supplement 2
Indirectly revealed preferences inherent in the choices shown in Figure 5C.
Higher order features of utility maximization.

(A) H and M food act as substitutes. Consumption is plotted against price for triplets of cost ratios in which the concentration of one food was constant and the concentration of the other food was variable (see Supplementary file 1). Lower case italic letters: data points in Figure 5B. Capital letters: the food whose density was constant, the consumption of which is plotted on the y-axis. Solid lines: regression slope different from zero (p0.01); dashed lines: slope not different from zero. Error bars, 95% CI. (B). H and M food are not perfect substitutes. Colored contour lines are indifference curves in a perfect substitute model (Equation 10) with β=10/11. Data points from Figure 5C are replotted for comparison, with associated budget lines, according to the conventions of that figure. ‘X’ symbols indicate the point of highest utility on each budget line. (C–D). Best fitting parameterizations of the CES function (Equation 11) for Trained and Untrained animals. Each panel shows the seven the iso-utility lines that are tangent to the budget lines. Goodness of fit can be assessed by observing that the iso-utility lines are tangent to the budget lines at, or near, the data points which indicate mean consumption of H and M food. For sample size (N), statistical methods used, and significance level, see Table 2, rows 27-32.

Behavioral mechanisms of utility maximization.

(A) Pumping-rate model of preference. Equation 13 is evaluated for each worm in Figure 5B and the result is plotted against preference in terms of fraction of pumps in H food for the same animal. (B). Dwell-time model of preference. Same as A but using Equation 14. (C). Regression of Equation 14 against mean head angle as defined in Figure 3B. (A–C). Blue lines: regressions on the data. (D). Diminished ceh-36 function eliminates the effect of food quality training on food preference. Asterisk, see Table 238. (H and M) are at OD = 1. N2 data are from Figure 4B. Replicates, Trained, strain(N): N2(20), ceh-36(7). Replicates, Untrained, strain(N): N2(28), ceh-36(11). Error bars, 95% CI. Correlation coefficients are shown ±95% CI. For sample size (N), statistical methods used, and significance level, see Table 2, rows 33-39.

Figure 8 with 1 supplement
Characterization of AWC’s response to delivery and removal of food.

(A) Ensemble averages of relative fluorescence versus time in response to onset of the indicated food in trained and untrained animals. (B). Summary of data in A, showing mean integrated calcium transients. (C). Ensemble averages of relative fluorescence versus time in response to removal of the indicated food in trained and untrained animals. Asterisk: untrained group, mean peak response, H vs. M food, see Table 242. Cross: trained group, mean peak response, H vs. M food, see Table 243. (D). Summary of data in C, showing mean integrated calcium transients Asterisk, see Table 246. (A–D). OD = 1 for H and M food. Shading, ± SEM. Error bars, 95% CI. For sample size (N), statistical methods used, and significance level, see Table 2, rows 40-47.

Figure 8—source data 1

Characterization of AWC's response to delivery and removal of food.

https://cdn.elifesciences.org/articles/69779/elife-69779-fig8-data1-v2.xlsx
Figure 8—figure supplement 1
Imaging chip.

Fluidic features are shown in black. Dimensions are in mm.

Figure 9 with 1 supplement
Characterization of AWC’s response to bacterial foods in Y-chip assays.

(A) Ensemble average of relative fluorescence versus time in response to transitions from H to M (H→M) or M to H (M→H) food. Asterisk, see Table 248. Both foods were presented at OD = 1. Food was switched at t = 0. H→M, Trained, N = 16; H→M, Untrained N = 14. M→H, Trained, N = 16; M→H, Untrained, N = 15. Shading,±SEM. (B). Typical fluorescence waveform in response to a series of transitions between H and M food. Notation: F0, baseline fluorescence after sustained exposure to H food (≅ 2 min.); FM, maximum H→M fluorescence; FH, minimum M→H fluorescence. The box delineates the time period of presumptive steady-state responses over which mean response amplitudes AF were computed for each recorded worm. (C–E). AF versus utility, preference, and log price ratio. Lines and curves are, respectively, linear and gaussian fits. Italic letters refer to labeled points in Figure 5B. Replicates, Trained, point(N): a(6), b(10), c(18), d(12), e(19), f(10), g(9). Replicates, Untrained. point(N): a(8), b(11), c(10), d(8), e(12), f(8), g(10). Error bars, 95% CI. For sample size (N), statistical methods used, and significance level, see Table 2, rows 48-55.

Figure 9—source data 1

Characterization of AWC’s response to bacterial foods in Y-chip assays.

https://cdn.elifesciences.org/articles/69779/elife-69779-fig9-data1-v2.xlsx
Figure 9—figure supplement 1
Control for mechanical artifacts in the experiment of Figure 9B assays.

This trace shows a typical AWC fluorescence waveform in response to a series of mock transitions between H and M food in which all syringes in the imaging system contained food-free buffer.

Price ratio curves and utility maximization.

(A) A pair of intersecting budget lines wherein choices a and b are governed by a price-ratio curve (not shown) that is monotonic-increasing. To support a direct violation of GARP, a must be to the right of the intersection of lines A and B, and b must be to the left of it. However, by price-ratio monotonicity, b must lie to the right of a, for reasons described in the text. This constraint precludes direct violations of GARP. (B). Choice data from a human participant in a GARP experiment (Camille et al., 2011). The data are consistent with utility maximization. Inset, price ratio data inferred from the experiment. Log price ratios for points a-k are calculated as logXZ/YZ, where XZ and YZ are, respectively, the x and y intercepts of budget line Z. The fraction of total budget spent on good X is computed as xZ/XZ, where xZ is the amount of good X chosen on line Z.

Videos

Video 1
Foraging behavior in the Y-chip.

Simulated Y-chip experiment. The worm is held at its midsection by a vacuum activated clamp, leaving the head (left) and tail free to move. The both fluid streams contain bacteria-free buffer, flowing to the right to left. Food dye was added to the lower stream to visualize the interface between streams. Bubbles originating at the clamp are formed by air that has been pulled through the PDMS walls of the chip by the vacuum. The worm prefers the dyed stream as it contains potassium sorbate, which acts as a chemoattractant.

Tables

Table 1
Volatile organic compounds released by H and M food.
DesignationCAS#CompoundH food (DA1877)M food (DA1885)
a1534-08-3Methylthioacetate*0.01.0
b624-92-0Dimethyl disulfide*19.438.3
c2432-51-1S-Methyl butanethioate0.04.7
d23747-45-7S-Methyl 3-methylbutanethioate0.02.4
e3658-80-8Dimethyltrisulfide*1.56.4
  1. Compounds were identified by gas chromatography-mass spectrometry of headspace of H and M food and confirmed with known standards. Amounts were inferred from the area under elution peaks, averaged across two replicates, and normalized to the amount of methylthioacetate in DA1885.

  2. *

    Identified chemoattractant; other compounds are uncharacterized in chemotaxis assays.

  3. Vapor is nematicidal.

Table 2
Statistics.

Horizontal location of cell entries varies by row.

RowFigureTestEffect or comparison testedUnits of replication or samplingNumber of replicates or samplesStatisticDF 1 or combined DFDF 2pEffect size
12B*t-testTrained N2 vs. Untrained N2Assay platesN=8t2.50142.56E-021.248
22Bt-testN2 Untrained, 60 min mean I>0Assay platesN=8t10.2271.86E-05
32Bt-testceh-36 Untrained, 60 min mean I>0Assay platesN=8t3.5879.01E-03
42Bt-testTrained ceh-36 vs. Untrained N2Assay platesN=8t0.64142.50E-01
52Bt-testUntrained ceh-36 vs. Untrained N2Assay platesN=8t1.87148.22E-02
62Bt-testTrained ceh-36 vs. Untrained ceh-36Assay platesN=8t0.90143.82E-01
72CTwo-factor ANOVA, repeated measures, main effectAzide- Trained vs. Azide- UntrainedAssay platesN≥9 / treatmentF11.281212.98E-030.349
82Ct-testAzide +Trained, 60 min mean I>0Assay platesN=10t6.35108.36E-05
92Ct-testAzide+, Untrained, 60 min mean I>0Assay platesN=10t3.09101.15E-02
102CTwo-factor ANOVA, repeated measures, main effectAzide +Trained vs. Azide +UntrainedAssay platesN=10 / treatmentF10.321204.37E-030.340
112CTwo-factor ANOVA, repeated measures, main effectAzide +Trained, vs. Azide– TrainedAssay platesN≥9 / treatmentF11.171193.43E-030.370
122CTwo-factor ANOVA, repeated measures, main effectAzide +Untrained vs. Azide– UntrainedAssay platesN≥10 / treatmentF38.281223.16E-060.635
132DTwo-factor ANOVA, repeated measures, main effectUntrained N2 vs. Untrained cat-2Assay platesN≥9 / strainF23.251219.14E-050.493
142DTwo-factor ANOVA, repeated measures, main effectTrained N2 vs. Trained cat-2Assay platesN=9 / strainF52.501189.74E-070.207
152DTwo-factor ANOVA, repeated measures, main effectTrained cat-2 vs. Untrained cat-2Assay platesN=9 / treatmentF0.901186.43E-01
164ATwo-factor ANOVA, main effectFamiliar vs. UnfamiliarWormsN≥19 / treatmentF10.461542.10E-030.162
174 A*t-testUnfamiliar, grown in H vs. grown in MWormsN≥6t2.65172.10E-020.876
184Bt-testTrained, f_H>0.5WormsN=19t8.60188.60E-08
194Bt-testUntrained, f_H>0.5WormsN=28t4.35271.70E-04
204B*t-testTrained vs. UntrainedWormsN≥19t2.95445.11E-030.850
214ETwo-factor ANOVAMain effect of optical densityWormsN≥22 / densityF31.5831089.80E-150.467
224FTwo-factor ANOVAMain effect of optical densityWormsN≥22 / densityF3.1031063.00E-020.081
235BTwo-factor ANOVAMain effect of price ratioWormsAvg N=15 / ratioF44.1361957.89E-340.576
245BTwo-factor ANOVAMain effect of trainingWormsAvg N=15 / ratioF36.1611958.82E-090.156
255Bt-testTrained, point a, mean f_H<0.5WormsN=9t6.2282.52E-04
265Bt-testUntrained, point a, mean f_H<0.5WormsN=12t8.29114.66E-06
276ARegression with replication slope testFigure 6A, points abd, Trained slope ≠ 0WormsN≥9 / ratioF118.791471.85E-14
286ARegression with replication slope testFigure 6A, points abd, Untrained slope ≠ 0WormsN≥10 / ratioF28.521394.26E-06
296ARegression with replication slope testFigure 6A, points cef, Trained slope ≠ 0WormsN≥10 / ratioF20.291464.54E-05
306ARegression with replication slope testFigure 6A, points cef, Untrained slope ≠ 0WormsN≥14 / ratioF26.561494.55E-06
316ARegression with replication slope testFigure 6A, points deg, Trained slope ≠ 0WormsN≥6 / ratioF2.341511.32E-01
326ARegression with replication slope testFigure 6A, points deg, Untrained slope ≠ 0WormsN≥7 / ratioF7.821457.56E-03
337ALinear correlationFrequency ratio vs. f_HWormsN=142t6.641416.53E-10
347BLinear correlationDwell time ratio vs. f_HWormsN=203t39.602026.95E-97
357CLinear correlationDwell time ratio vs. mean head angleWormsN=203t35.962022.62E-89
367Dt-testceh-36 Untrained, f_H>0.5WormsN=11t4.20101.84E-03
377DTwo-factor ANOVATreatment ×Strain interactionWormsN≥7 / treatmentF5.031622.85E-020.075
387D*t-testN2, Trained vs. UntrainedWormsN≥20t3.45461.21E-030.186
397Dt-testceh-36, Trained vs. UntrainedWormsN=7 / treatmentt0.58165.83E-01
408BTwo-factor ANOVAMain effect of food typeWormsN≥6 / treatmentF3.561237.20E-02
418DTwo-factor ANOVAMain effect of food typeWormsN≥7 / treatmentF18.421252.00E-040.424
428 C*t-testPeak response, Untrained, H vs. M foodWormsN≥7 / treatmentt2.98121.30E-020.913
438 C×t-testPeak response, Trained, H vs. M foodWormsN≥7 / treatmentt2.29133.96E-021.184
448BTwo-factor ANOVAMain effect of trainingWormsN≥6 / treatmentF0.001239.90E-01
458DTwo-factor ANOVAMain effect of trainingWormsN=7 / treatmentF7.521251.10E-020.883
468D*t-testH food, Trained vs. UntrainedWormsN=7 / treatmentt2.86121.44E-021.528
478Dt-testM food, Trained vs. UntrainedWormsN≥7 / treatmentt0.80134.39E-01
489 A*t-testH → M, peak response, Trained vs. UntrainedWormsN≥14 / treatmentt2.66281.29E-020.972
499At-testM → H, area under the curve, Trained vs. UntrainedWormsN≥15 / treatmentt0.74294.67E-01
519CLinear correlationAWC activation vs. utility, TrainedWormsN≥6 / meant0.1059.27E-01
529CLinear correlationAWC activation vs. utility, UntrainedWormsN≥8 / meant0.1758.75E-01
539DLinear correlationAWC activation vs. preference, TrainedWormsN≥6 / meant1.5751.77E-01
549DLinear correlationAWC activation vs. preference, UntrainedWormsN≥8 / meant0.4456.78E-01
559Ct-testTrained, point e vs. Untrained, point dWormsN≥8 / meant2.14254.23E-020.889
  1. p-values associated with significant results are shown in bold font. Sample size was determined by increasing the number of biological replicates until the coefficient of variation for each means converged. Each experiment was performed once, with the indicated number of biological replicates; the non-stationary nature of the organism precluded technical replicates. No data were censored or excluded.

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Strain, strain background (E. coli)OP50CGC (C. elegans Genetic Center)RRID:WB-STRAIN:WBStrain00041969
Strain, strain background (Comamonas sp.)DA1877CGCRRID:WB-STRAIN:WBStrain00040995
Strain, strain background (Bacillus simplex)DA1885CGCRRID:WB-STRAIN:WBStrain00040997
Genetic reagent (C. elegans)N2, BristolCGCRRID:WB-STRAIN:WBStrain00000001
Genetic reagent (C. elegans)CX5893CGCRRID:WB-Strain00005275
Genetic reagent (C. elegans)JP5651National Resource Project (Japan)Nonecat-2(tm2261)
Genetic reagent (C. elegans)MT15620CGCRRID:WB-Strain00027527cat-2(n4547)
Genetic reagent (C. elegans)CB1112Stern et al., 2017; PMID:29198526RRID:WB-Strain00004246cat-2(e1112)
Genetic reagent (C. elegans)XL322Lockery labNonentIs1703[str-2::GCaMP6s-wcherry; unc-122::dsred2]
Recombinant DNA reagentGCaMP-6s::wCherryZhen labNone
Software, algorithmIgor ProWavemetrics https://www.wavemetrics.com/Version 9.01

Additional files

Supplementary file 1

Design of the analysis in Figure 6A.

Boxes are groups of choice sets in which the price of one food was constant (shading) while the other was variable. Numbers are optical density from which price was computed by Equation 2. Lower case letters refer to data points in Figure 5B.

https://cdn.elifesciences.org/articles/69779/elife-69779-supp1-v2.docx
Supplementary file 2

Linear correlation equations.

These equations show how the economic variables in the left column were computed based on food type, density and preference in trained and untrained animals. The quantities wT and wU are defined in Equation 15 in the main text. The quantities UT and UU are, respectively, utility in trained and untrained animals computed according to the CES function as fitted to the data in Figure 6C and D.

https://cdn.elifesciences.org/articles/69779/elife-69779-supp2-v2.docx
Supplementary file 3

Tests of linear correlations between AWC activation and economic variable it might hypothetically represent.

Significance of correlations was tested using the F distribution. This table shows, for each economic variable tested, the correlation coefficient, value of the F statistic, its two degrees of freedom, and the corresponding p-value. No significant correlations where found. For definitions of variables see Supplementary file 2.

https://cdn.elifesciences.org/articles/69779/elife-69779-supp3-v2.docx
Supplementary file 4

T-maze CAD file.

https://cdn.elifesciences.org/articles/69779/elife-69779-supp4-v2.zip
Supplementary file 5

Y-chip CAD file.

https://cdn.elifesciences.org/articles/69779/elife-69779-supp5-v2.zip
Supplementary file 6

Imaging chip CAD file.

https://cdn.elifesciences.org/articles/69779/elife-69779-supp6-v2.zip
MDAR checklist
https://cdn.elifesciences.org/articles/69779/elife-69779-mdarchecklist1-v2.pdf

Download links

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

Downloads (link to download the article as PDF)

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

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

  1. Abraham Katzen
  2. Hui-Kuan Chung
  3. William T Harbaugh
  4. Christina Della Iacono
  5. Nicholas Jackson
  6. Elizabeth E Glater
  7. Charles J Taylor
  8. Stephanie K Yu
  9. Steven W Flavell
  10. Paul W Glimcher
  11. James Andreoni
  12. Shawn R Lockery
(2023)
The nematode worm C. elegans chooses between bacterial foods as if maximizing economic utility
eLife 12:e69779.
https://doi.org/10.7554/eLife.69779