Figures and data

Odors alone evoke PER.
(A) Schematic for recording odor-evoked PER. (B) Six points on the head tracked by DeepLabCut. The movement of the proboscis was characterized by the positions of three segments of the proboscis and the angles between the segments, namely the rostrum angle (R), haustellum angle (H), and labellum angle (L). (C) Example odor-evoked PER. Top, images of a proboscis in a retracted (0.2 s), fully extended (3.1 s) and partially extended (5.5 s) state. Bottom, the three angles over time. Green bar indicates an odor application period. Odor is ethyl butyrate. (D) PER to nine odors in an example fly. Each row corresponds to a trial and each tick mark indicates the timing of PER. The odors used are as follows: ethyl butyrate (EBR), 2-pentanone (2PT), isopentyl acetate (IPA), 4-methylcyclohexanol (MCH), benzaldehyde (BNZ), 3-octanol (OCT) and mineral oil (MO). (E) PER probability in wild type flies in response to different concentrations of odors. The same stimulus was used for air, mineral oil, and water controls. n = 21, 26, 27, and 28 flies for 10-4, 10-2, 10-1, and 0.5 concentration groups, respectively. Error bar, standard error of the mean. PER probability was different between odor concentrations and identities (p = 3.0e-05 and 2.0e-7 for odor concentration and identity factors, Scheirer–Ray–Hare test). (F) PER duration in wild type flies. The black lines indicate the mean and the shaded areas indicate standard error of the mean. n = 28 flies. (G) Integrated PER duration in wild type flies in response to different concentrations of odors. The data is from the same flies as in E. Integrated PER duration was different between odor concentrations and identities (p = 3.6e-13 and 3.5e-11 for odor concentration and odor identity factors, Scheirer–Ray–Hare test). Box plots indicate the median (gray line), mean (black dot), quartiles (box), and 5-95% range (bar).

Odors evoke PER through GRNs.
(A) Integrated PER duration in intact (n = 28) and olfactory organs removed (n = 32) wild type flies. Flies with olfactory organs removed show less PER compared to control flies (p = 0.0079 and 0.10 for group and odor factors, Scheirer–Ray–Hare test). (B-F) Integrated PER duration in UAS control (UAS-Kir2.1/+, n = 20, B), GRN silenced (Gr5a-Gal4>UAS-Kir2.1, n = 25 for C; Gr66a-Gal4>UAS-Kir2.1, n = 24 for D; Ppk28-Gal4>UAS-Kir2.1, n = 28 for E; Ir94e- Gal4>UAS-Kir2.1, n = 23 for F) and Gal4 control flies (Gr5a-Gal4/+, n = 31 for C; Gr66a- Gal4/+, n = 21 for D; Ppk28-Gal4/+, n = 28 for E; Ir94e-Gal4/+, n = 30 for F). Each silenced line was compared to the UAS control in B (p = 4.9e-5, 0.49, 0.17, 0.11 and 0.12, 0.0014, 0.55, 1.2e-5 for genotype and odor factors for C-F, Scheirer–Ray–Hare test). Data were also compared between each silenced line and its Gal4 control (p = 1.0e-15, 0.28, 0.00062, 0.14 and 0.71, 0.0012, 0.56, 5.6e-9 for genotype and odor factors for C-F, Scheirer–Ray–Hare test). (G) Integrated PER duration in control (UAS-Kir2.1/+, n = 21) and Gr66a GRN silenced (Gr66a- Gal4>UAS-Kir2.1, n = 21) flies in a fed state. Silencing of Gr66a GRNs enhances odor-evoked PER in a fed state (p = 1.8e-4 and 0.041 for genotype and odor factors, Scheirer–Ray–Hare test). Odor concentration was 10-1. (H) Integrated PER duration in control (n = 29) and Gr5a mutant (n = 29) flies. PER is severely reduced in mutant flies (p = 1.1e-4 and 0.32 for genotype and odor factors, Scheirer–Ray–Hare t est). Control: EP(x)496 flies. Mutant: ΔEP(x)-5 flies. (I-K) Integrated PER duration in control and Obp RNAi flies. Obp57d/e RNAi flies show reduced (I, n = 33 and 25 for control and RNAi, p = 0.039 and 0.48 for genotype and odor factors, Scheirer– Ray–Hare test), Obp49a RNAi flies show enhanced (J, n = 24 and 20 for control and RNAi, p = 0.0066 and 0.16 for genotype and odor factors, Scheirer–Ray–Hare test), and Obp19b RNAi flies show similar level of PER (K, n = 24 and 21 for control and RNAi, p = 0.24 and 0.26 for genotype and odor factors, Scheirer–Ray–Hare test) as compared to their respective controls, which are tubulin-Gal4/y,w1118;P{attP,y[+],w[3’]} (for I), tubulin- Gal4/w1118;P{VDRCsh60200}attP40 (for J), and tubulin-Gal4/w1118 (for K). Box plots indicate the median (gray and red lines), mean (black and red dots), quartiles (box), and 5-95% range (bar).

Sweet- and bitter-sensing GRNs directly respond to odors,
(A) Schematic for two photon calcium imaging of GRN axon termini in the subesophageal zone. The top image shows an anterior view of a brain of Gr5a-Gal4>UAS-GCaMP6s fly. White rectangle indicates the target region of calcium imaging. Scale bar: 50 μm. (B) Example response (ΔF/F) to ethyl butyrate. Green bar indicates an odor application period. Odor concentration was 10-1 for all the odors in all the experiments described in this figure. (C) Summary responses of Gr5a GRNs to nine odors. Green bar indicates an odor application period. ΔF/F of GCaMP6s fluorescence is color coded according to the scale bar (n = 5). (D) Trial averaged peak responses of Gr5a GRNs to individual odors. Covering the labella with glue reduced odor responses as compared to control (n = 5 and 4 for control and labella covered, p = 2.6e-15 and 5.1e-3 for organ condition and odor factors, mixed two-way ANOVA) but covering the maxillary palps had no effects (n = 5 and 4 for control and palps covered, p = 0.69 and 0.055 for organ condition and odor factors, mixed two-way ANOVA). (E, F) Same as in C, D, but for Gr66a GRNs. As in Gr5a GRNs, covering the labella with glue reduced odor responses as compared to control (n = 6 and 5 for control and labella covered, p = 3.0e-16 and 1.1e-7 for organ condition and odor factors, mixed two-way ANOVA) but covering the maxillary palps had no effects (n = 6 and 5 for control and palps covered, p = 0.48 and 4.2e-11 for organ condition and odor factors, mixed two-way ANOVA). (G) A linear combination of Gr5a and Gr66a GRN responses well predicted the magnitude of odor-evoked PER. PER was measured in Gr5a-Gal4>UAS-GCaMP flies. Each dot represents an odor. Coefficient of determination = 0.81. Weights for Gr5a and Gr66a GRNs are 5.0e-4 and −2.3e-5, indicating that the former and the latter contributes to enhancement and suppression of PER, respectively. (H) Odor responses are severely reduced in Gr5a GRNs in a Gr5a receptor mutant background (n = 5 and 8 for control and mutant, p = 8.0e-06 and 0.037 for genotype and odor factors, mixed two-way ANOVA). Box plots indicate the median (gray and red lines), mean (black and red dots), quartiles (box), and 5-95% range (bar).

Individual taste sensilla sense odors.
(A) Example responses of individual taste sensilla recorded using tip recording. Left, example responses to TCC (30 mM), MCH (1%), banana odor (1%), and sucrose (100 mM) recorded in the same L-type sensillum; Right, example responses to the same set of stimuli in the same I- type sensillum. Traces show the activity between 200 and 700 ms after the stimulus onset. Each gray dot indicates a spike. Bar indicates 50 ms. (B) Average responses of L-type sensilla. (n = 20, 6, 17, 20 flies for TCC, MCH, banana odor, and sucrose, respectively). Responses to stimuli are significantly larger than those to TCC (p = 3.0e-18, one-way ANOVA followed by Tukey’s HSD, p = 0.0051, 0.0021, and 2.7e-7 for MCH, banana, and sucrose compared to TCC). (C) Average responses of I-type sensilla. (n = 30, 9, 22, 30 flies for TCC, MCH, banana odor, and sucrose, respectively). Responses to stimuli are significantly larger than those to TCC (p = 1.4e- 12, one-way ANOVA followed by Tukey’s HSD, p = 3.9e-6, 2.0e-7, and 5.9e-12 for MCH, banana, and sucrose compared to TCC). Box plots indicate the median (gray line), mean (black dot), quartiles (box), and 5-95% range (bar).

Odor-taste integration in single taste sensillum.
(A) Example responses of an L-type taste sensillum to sucrose (0.25%, top), banana odor (0.1%, middle) and the mixture of the two (bottom). Traces show the activity between 200 and 700 ms after the stimulus onset. Each gray dot indicates a spike. Bar indicates 50 ms. (B) Average responses of L-type sensilla. Responses to the mixture are significantly larger than those to components (n = 14 flies, p = 0.001, repeated measures ANOVA followed by paired t-tests with Bonferroni correction, p = 0.0030 and 0.0032 for sucrose and banana compared to mixture). Box plots indicate the median (gray line), mean (black dot), quartiles (box), and 5-95% range (bar).

GRNs integrate multisensory input to enhance PER and food consumption.
(A) Schematic for examining PER in response to odor-taste multimodal stimuli. Stimuli are presented to the labellum with a wick (arrowhead) immersed in a sucrose solution with or without odors. A constant air stream was applied from behind the wick. (B) PER probability in response to sucrose at different concentrations with (green) or without (black) banana odor in wild type flies. Odor concentration was 10-4. Banana odor enhances sucrose-evoked PER not only in intact flies (n = 13 and 11 for control and banana odor, p = 4.6e-4 and 2.1e-7 for odor identity and sucrose concentration factors, Scheirer–Ray–Hare test) but also in flies without olfactory organs (n = 11 and 11 for control and banana odor, p = 0.035 and 0.0070 for odor identity and sucrose concentration factors, Scheirer–Ray–Hare test). Error bar, standard error of the mean. (C) PER probability in response to sucrose is increased with EBR in both intact and olfactory organs-removed flies (n = 25 and 25 for intact and olfactory organs-removed flies, p = 0.0052 and 0.0068 for intact and olfactory organ-removed flies, Wilcoxon signed-rank test). (D) Same as in C, but for MCH (n = 21 and 36, p = 0.0072 and 0.00074 for intact and olfactory organs-removed flies, Wilcoxon signed-rank test). The concentration of odors and sucrose was 10-4 and 0.25% for C and D. Error bar, standard error of the mean. (E) Schematic for quantifying sucrose consumption in the presence or absence of odors. A filter paper on the top was impregnated with sucrose solution containing blue dye and a small piece of filter paper placed in the mesh tube near the top was impregnated with odor or control solution (see Methods). (F) The amount of sucrose consumed by flies with or without the presence of banana odor was quantified with a feeding score (see Methods). The amount of feeding was significantly larger in the presence of odor for both intact (n = 20 and 22 for control and banana odors, p = 0.005, Wilcoxon signed-rank test) and olfactory organs-removed flies (n = 13 and 15 for control and banana odors, p = 0.007, Wilcoxon signed-rank test). (G) The same as in F but with EBR and MCH odors and only for olfactory organs-removed flies (EBR experiment, n = 16 and 17 for control and EBR, p = 0.03, Wilcoxon signed-rank test; MCH experiment, n = 11 and 14 for control and MCH, p = 0.24, Wilcoxon signed-rank test). (H) Summary of mechanisms underlying PER and food intake in a multisensory environment. Error bars are standard error of the mean in bar plots. Box plots indicate the median (gray line), mean (black dot), quartiles (box), and 5-95% range (bar).


Genotypes used to generate data presented in each figure.

Summary of statistical tests.

Characteristics of PER in wild type flies/
(A) Latency to the first odor-evoked PER in wild type flies. Odor concentration was 0.5 (n = 28 flies). Green bar indicates the odor application period. p = 1.0e-4, one-way ANOVA. Box plots indicate the median (gray line), mean (black dot), quartiles (box), and 5-95% range (bar). (B) K- medoid clustering of odor- and sucrose-evoked proboscis extensions based on the rostrum and haustellum angles over time in wild type flies. All the trajectories during PER and partial proboscis extensions were clustered into five clusters (color coded), and the representative trajectory closest to each cluster medoid was plotted for each stimulus. The green trajectories represent full extensions that qualify the definition of PER (see Methods). Note that these green trajectories are observed in response to sucrose and all the odors (except for the controls), suggesting that the movement of proboscis during odor-evoked PER is similar to that during sucrose-evoked PER. (C) Integrated PER duration in response to different odors in wild type flies analyzed by DeepLabCut or manual scoring. The values measured by the two methods were similar (p = 0.94 and 5.5e-3 for scoring method and odor identity factors, Scheirer-Ray-Hare test). Odor concentration is 0.1. Box plots indicate the median (gray line), mean (black dot), quartiles (box), and 5-95% range (bar).

PER evoked by additional odors.
Integrated PER duration in wild type flies in response to 1-hexanol, ethyl acetate, yeast odor, fenchone, isoamyl alcohol, phenethyl alcohol, air, mineral oil and water. In this experiment, flies were fixed in a vertical position following a protocol in a previous study (Oh et al., 2021). Box plots indicate the median (gray line), mean (black dot), quartiles (box), and 5-95% range (bar).

Odor-evoked PER in flies with different genetic backgrounds and in another species.
Integrated PER duration in different flies (n = 20, 28 and 25 in A, B and C). Box plots indicate the median (gray and red lines), mean (black and red dots), quartiles (box), and 5-95% range (bar).

Dependence of odor-evoked PER on a feeding state.
(A) Integrated PER duration in starved (n = 28) and fed (n = 19) wild type flies (p = 0.88 and 5.6e-3 for state and odor factors, Scheirer–Ray–Hare test). Odor concentration is 0.5. (B) Same as in a but for odor concentration of 10-1. PER is enhanced in a starved state (p = 1e-6 and 0.74 for state and odor factors, Scheirer–Ray–Hare test). (C) PER probability in wild type fed flies in response to different concentrations of odors. n = 19 and 28 for 10-1 and 0.5 concentration groups, respectively. Error bar, standard error of the mean. Box plots indicate the median (gray and red lines), mean (black and red dots), quartiles (box), and 5-95% range (bar).

Wings and legs are not required for odor-evoked PER.
(A) Integrated PER duration in wing-removed Orco-Gal4>UAS-Kir2.1 flies (n = 26). GRNs in wings were dispensable for PER to certain odors (p = 6.0e-9, Kruskal-Wallis test, p = 0.0020, 0.0011, 0.75, 0.089, 0.00071, 0.58, 0.84, and 0.94 for comparison between mineral oil (MO) control vs EBR, 2PT, IPA, MCH, BNZ, OCT, Air, and Water, Dunn’s post-hoc tests with Benjamini–Hochberg false discovery rate correction). (B) Integrated PER duration in wing- and leg-removed, Orco-Gal4>UAS-Kir2.1 flies (n = 27). GRNs in wings were dispensable for PER to certain odors (p = 0.004, Kruskal-Wallis test, p = 0.39, 0.029, 0.18, 0.18, 0.029, 0.44, 0.88, and 0.71 for comparison between mineral oil (MO) control vs EBR, 2PT, IPA, MCH, BNZ, OCT, Air, and Water, Dunn’s post-hoc tests with Benjamini–Hochberg false discovery rate correction). Box plots indicate the median (gray line), mean (black dot), quartiles (box), and 5-95% range (bar).

Contributions of OBPs to odor-evoked PER.
(A, B) Integrated PER duration in control (tubulin-Gal4/w1118) and Obp56g RNAi flies. These flies show similar level of PER (n = 24 and 23 for control and RNAi, p = 0.58 and 0.0017 for group and odor factors, Scheirer–Ray–Hare test). (C-H) Integrated PER duration in control (tubulin-Gal4/y,w1118;P{attP,y[+],w[3’]}) and OBP RNAi flies. Obp18a RNAi and Obp57e RNAi flies show reduced level of PER as compared to control (n = 33, 25, 22, 25, 27, and 28 flies for C-H, p = 0.050, 0.36, 0.18, 1.9e-4, 0.92 and 0.10, 0.48, 0.063, 0.045, 0.022 for group and odor factors for D-H, Scheirer–Ray–Hare test). Box plots indicate the median (gray and red lines), mean (black and red dots), quartiles (box), and 5-95% range (bar).

Responses of GRNs to odors and sucrose.
(A) Image of a brain of Gr5a-Gal4>UAS-GCaMP6s fly, stained with anti-GFP to reveal Gr5a neurons (green) and nc82 to reveal synapses (magenta). Scale bar, 50 μm. (B) Baseline florescence intensity of GCaMP6s in Gr5a GRNs. Fluorescence (arbitrary unit) is color coded according to the scale bar. (C) Peak florescence intensity of GCaMP6s in the same fly as shown in B after stimulation with 4-methylcyclohexanol. Fluorescence is color coded according to the scale bar. The ROI used for ΔF/F calculation is indicated by a white rectangle. Scale bar, 20 μm. (D-F) Same as in A-C, but for Gr66a-Gal4>UAS-GCaMP6s fly and benzaldehyde stimulation. (G) Trial averaged peak responses of Gr5a GRNs to different concentrations of ethyl butyrate (EBR) and benzaldehyde (BNZ). n = 5 flies. (H) Same as in G, but for Gr66a GRNs. (I) Responses to 4-methylcyclohexanol are positively correlated with that to sucrose (0.25%) in Gr5a GRNs (R2 = 0.73, p = 0.031). Each dot corresponds to a fly. (J) Trial averaged peak responses of Ppk28 GRNs to individual odors. These neurons do not respond to odors. n = 5 flies. (K) Same as in J, but for Ir94e GRNs. n = 9 flies. Box plots indicate the median (gray line), mean (black dot), quartiles (box), and 5-95% range (bar).

Individual responses of GRNs to odors.
(A) Time course of calcium responses (ΔF/F) to odors in Gr5a GRNs. Responses to nine odors in each fly are shown in one panel (n = 5 flies). Odor application period is 5 to 7 s, and each odor was applied three times (different trials are shown in different colors). (B) Same as in A, but for Gr66a GRNs.

PER and GRN responses to various concentrations of banana odor.
(A) Integrated PER duration in response to a banana odor presented at various concentrations in wild type flies. n = 21 flies. Box plots indicate the median (gray line), mean (black dot), quartiles (box), and 5-95% range (bar). (B) Trial averaged peak responses of Gr5a GRNs in response to various concentrations of banana odor in Gr5a-Gal4>UAS-GCaMP6s flies. n = 5 flies. Error bar, standard error of the mean. Box plots indicate the median (gray line), mean (black dot), quartiles (box), and 5-95% range (bar).

Correlation between odor and taste responses in individual sensilla.
Spike responses to banana odor (0.01%) are positively correlated with those to sucrose (0.25%) in L-type sensilla. Each dot corresponds to a fly (R2 = 0.62, p = 0.00044).