Response spectra of the larval odor receptor repertoire
A consensus set of 25 larval Or receptors was identified on the basis of three studies of
Or gene expression (
Couto et al., 2005;
Kreher et al., 2005). The expression of these 25 genes in larval ORNs has been documented either by in situ hybridization (
Fishilevich et al., 2005) or, in the case of
Or49a, by reporter gene expression driven by an
Or promoter via the
GAL4/UAS system (
Fishilevich et al., 2005;
Kreher et al., 2005). We have expressed all 25 genes in the empty neuron system and tested them with a diverse panel of odorants. A previous work found that the response spectrum, response mode, i.e. excitation or inhibition, and temporal dynamics of antennal receptors were the same in the empty neuron as in the ORN in which they were normally expressed (
Hallem et al., 2004b). The present analysis includes 14 receptors that were not tested in our earlier study of larval odor receptors (
Kreher et al., 2005), and it extends the analysis of the other 11 receptors.
Of the 25
Or genes, 21 conferred odor responses. Alleles of the remaining four genes
(Or1a, Or33a, Or63a, and
Or83a) in our laboratory Canton-S strain appeared to be non-functional. They did not confer a regular pattern of spontaneous action potentials, and they did not confer any odor-evoked responses when expressed from genomic constructs in the empty neuron system. Or1a and Or33a each contain an amino acid polymorphism, in comparison to another wild-type strain, Oregon-R. The Oregon-R alleles conferred a regular pattern of spontaneous action potentials; the Oregon-R allele of Or33a did not confer strong responses (≥100 spikes/s) to any odorant of our panel or to mixtures representing 100 odorants (
Figure S1A and data not shown). The Oregon-R allele of Or1a yielded a response >50 spikes/s only to a mixture containing methyl benzoate, which activates an Or1a-expressing ORN in that strain (
Fishilevich et al., 2005;
Louis et al., 2008). Not only did Or63a not function when expressed from Canton-S genomic DNA, but it did not function when expressed from Canton-S cDNA, or from Oregon-R genomic DNA. We recovered a second allele of
Or83a from Canton-S that contains a premature stop codon. In summary, the simplest interpretation of our results is that in our Canton-S laboratory strain, 21
Or genes encode functional odor receptors.
The 21 receptors were systematically tested with 26 diverse odorants, each examined at two concentrations, and CO
2; thus 1113 receptor-odorant combinations were tested (). The odorants were selected to include ketones, aromatics, alcohols, esters, aldehydes, a terpene, and an organic acid, and they include compounds of varying chain lengths. We initially tested the odorants as “10
-2” dilutions (see Experimental Procedures); this dosage has been used in many previous studies (
De Bruyne et al., 2001;
Hallem et al., 2004b;
Kreher et al., 2005;
Larsson et al., 2004;
Wilson et al., 2004) and allows a convenient assessment of the response spectra of receptors. Lower concentrations were administered as 10
-4 dilutions (
Hallem and Carlson, 2006). At the higher concentration tested, 19% of the receptor-odorant combinations yielded an excitatory response of ≥50 spikes/s, 12% yielded a strong response, defined here as ≥100 spikes/s, 8% yielded a response of ≥150 spikes/s, and 4% yielded a response of ≥200 spikes/s ( and
Figure S2A).
Receptors varied a great deal in the number of odorants that excited them strongly. At one extreme, Or67b was strongly excited by eight odorants (30% of the panel) and Or42a and Or85c were strongly excited by seven at the higher concentration. At the other extreme, four of the 21 receptors were not strongly excited by any odorants. The receptors showed a smooth distribution between these two extremes (
Figure S2A).
Odorants varied according to the number of receptors that they strongly excited. 1-hexanol elicited a strong excitatory response from six receptors at the higher concentration, whereas four compounds did not strongly excite any receptor ().
Inhibitory responses, defined as a reduction in the firing rate to ≤50% of the spontaneous level, were also observed (,
Figure S2A). At the higher doses, most odorants inhibited at least one receptor, and most receptors were inhibited by at least one odorant. Our estimate of the extent of inhibition should be regarded as an approximation. For example, it is difficult to observe inhibition in the receptors with low spontaneous firing rates. We also note that we have adopted a conservative definition of inhibition.
A small fraction of the odorant-receptor combinations tested here were also evaluated in a previous study of antennal receptors (
Hallem and Carlson, 2006). The results are generally in good agreement (see Experimental Procedures).
At lower odorant concentrations, strong excitatory responses are much more sparse. When the concentration was reduced by two orders of magnitude, the percentage of combinations producing ≥50 spikes/s fell from 19% to 4%, the percentage producing ≥100 spikes/s fell from 12% to 2%, and the percentage producing ≥200 spikes/s fell from 4% to 1% (). Of the nine receptors that responded strongly to at least one odorant at low concentrations, six responded strongly to one odorant only. No odorant elicited a strong response from more than one receptor.
We have established tuning curves for each receptor of the larval repertoire, at each concentration tested (). When tested with 10-2 dilutions, the tuning curves showed a range of tuning breadths, with some receptors appearing very broadly tuned. By contrast, when tested with 10-4 dilutions, none of the receptors appeared broadly tuned, and many appeared more narrowly tuned.
Behavioral responses
Having thus examined the physiological responses generated by each odor of the panel, we then measured the behavioral responses of the larva to each odor. We used a simple behavioral paradigm in which larvae migrate across a Petri plate toward a source of odorant (). A response index (RI) was calculated at the end of a 5 min test period, as in earlier studies (
Monte et al., 1989;
Rodrigues and Siddiqi, 1978), by counting the number of animals, S, on the half of the plate containing the odorant, subtracting the number, C, on the control half, and dividing by the total; thus RI=(S-C)/(S+C). The RI ranges from 1 (complete attraction) to 0 (no response) to -1 (complete repulsion).
The odors of the panel elicited a broad range of responses (). Most odors were attractive at the test concentration, with E2-hexenal eliciting the strongest response: RI= 0.74±0.07 (n=10 trials). Among repellents, the strongest was geranyl acetate: RI= -0.40±0.07 (n=10 trials; for statistical analysis of repellency, see Experimental Procedures). Odorants of diverse chemical classes elicited strong responses of the same response mode, i.e. the strongest attractants (E2-hexenal, propyl acetate, 2-heptanone) are an aldehyde, ester, and ketone, respectively. Likewise, the repellents included a terpene, a ketone, an alcohol, and aromatics.
Predicting perceptual quality from odor space: odor masking
We used the responses of the receptor repertoire to examine perceptual relationships among odorants. We constructed a 21-dimensional odor space in which each dimension represents the response of one odor receptor in spikes/s. Odorants were plotted in this space based on the responses they elicited at a 10-2 dilution. We then calculated distances pairwise between all odorants. Two kinds of distances were determined. First, Euclidean distances were calculated pairwise between all odorants. Second, angular distances were calculated, based on the angle between the two vectors that extend from the origin of the 21-dimensional space to each of the two odorants (see Experimental Procedures).
From the matrix of distances we constructed two separate three-dimensional projections of odor space, one based on the Euclidean distances and one based on the angular distances, using Multi-Dimensional Scaling (MDS; ). The closest and most distant pairs of odorants in terms of Euclidean distance are shown in . Members of the closest pairs are structurally related. (We exclude from consideration three odorants that fail to elicit a response of >50 spikes/s from any receptor and that map near each other for this reason.) The most distant relationships, by contrast, are in each case between an aliphatic and an aromatic odorant. In the case of angular distance, the most closely related pairs are again structurally related, and the distant pairs are again structurally dissimilar (). The two closest pairs in angular distance are also two of the closest in Euclidean distance. Overall, the Euclidean and angular distances were generally concordant for pairs of odors that were close in odor space, but showed less agreement for pairs of odors that were more distant (
Figure S3).
We note that in general the Euclidean distance is more sensitive to the magnitudes of receptor responses, whereas the angular distance is more sensitive to the “pattern” of receptor firing, i.e. to the identity of the responding receptors. (Consider, for example, two vectors in a three-dimensional space whose directions are very similar; increasing the magnitude of one vector will have an effect on the Euclidean distance but not on the angular distance.) The Euclidean distances among odors were used as the basis of a hierarchical cluster analysis (
Figure S4), which shows that in a number of cases structurally similar odors clustered together, although in no case did all odorants of a particular structural class cluster together.
We then asked whether two odorants that are close in odor space are close in perceptual qualities. There are several paradigms for examining perceptual relationships between two odorants. We have examined odor masking, a phenomenon that is convenient to measure and may be of direct importance to animals in their natural environment by affecting what odors they perceive.
The odor masking paradigm is based on the behavioral paradigm shown in . Responses to a point source of odorant A are measured in the presence of a background of odorant B; odorant B is uniformly distributed across a filter that covers nearly the entire lid of the Petri dish in which the assay is conducted. We then determine whether the background of odorant B decreases the response of the animals to odorant A. This paradigm was used in a previous study of
Drosophila larvae (
Rodrigues, 1980), and a very similar paradigm has been used in
C. elegans (
L'Etoile and Bargmann, 2000;
Wes and Bargmann, 2001): if the worms responded to odor A in a background of B, then they were inferred to be able to discriminate A from B; if they cannot discriminate A from B, then a background of A was expected to block the response to B. We will interpret the results not in terms of discrimination but in terms of masking. The ability to identify odors in a background of other odors has also been described as odor segmentation (
Wilson and Mainen, 2006).
As an initial test we asked whether the response to a point source of odorant A was decreased by the presence of a background of odorant A. We found that the response to a 10-4 dilution of ethyl acetate was in fact decreased by a background of ethyl acetate (). The mean response to the ethyl acetate source decreased progressively as the dose of the background odor was increased, over six orders of magnitude: a severe decrease was observed when the dilution of the masking odorant reached 10-4, and the response was abolished when the masking odorant reached a 10-3 dilution.
We then systematically tested the odorants of our panel at a 10-2 dilution for their ability to block the response to ethyl acetate. The background odorants varied across a broad range in their effects, from complete masking to no masking (). In addition to ethyl acetate, a background of propyl acetate, ethyl butyrate, or 1-hexanol reduced the RI to below 0.2.
The analysis was then extended to examine masking of responses to five additional odorants that elicited strong behavioral responses: 2-heptanone, 3-octanol, E2-hexenal, ethyl butyrate, and 2,3-butanedione. We tested the ability of ethyl acetate and other odorants to mask the responses to these five odorants ().
We found that at least one odor masked the response to a point source of each of the five tested odorants. When ethyl butyrate was used as the point source, ethyl acetate masked the response as potently as did ethyl butyrate itself. Reciprocally, among the five odorants used as a point source, ethyl butyrate most potently masked response to ethyl acetate (). When all the masking data () were pooled, yielding 48 odorant pairs, there were 13 cases (27%) in which the RI to point source odorants was reduced by ≥50% compared to the control value observed when paraffin oil was used as the masking odor.
We then asked whether there was a correlation between the extent to which odor A masked odor B and the distance between them in odor space. We defined a masking index, MI, as a ratio of the RI in the absence and presence of the masking odorant (see Experimental Procedures). Thus, a high MI value indicates a high degree of masking, and a low MI value indicates a low degree of masking.
We found that odor masking correlated with distance in odor space. When angular distances were used to measure distance in odor space, we found an R2 value of 44.0% by regression analysis (p<0.0009; , black line). Moreover, for most of the point source odors in this analysis, the masking odors do not provide a broad sampling of odor space; when we limit the analysis to the point source odor ethyl acetate, in which all odors of the panel were sampled as masking odors, R2=55.9% (p<0.0009; , orange line). For 2-heptanone, which was also tested with a broader sampling of masking odors, R2 was even higher: 58.4% (p<0.0009; , red line).
When Euclidean distances were used to analyze the entire data set, a predictive relationship was still clearly observed, but the variation explained was lower than when angular distances were used: R2=35.1%, as opposed to R2=44.0% for angular distances.
Behavioral response to one odorant across a concentration range depends on two receptors
Having observed a relationship between perception and the responses of the entire receptor repertoire, we next addressed how the responses of individual receptors are integrated to engender behavioral responses. We began by examining the response to ethyl acetate, which elicits behavioral responses over doses spanning more than six orders of magnitude (), but which elicits strong responses from only two receptors, Or42a and Or42b ().
Or42a and Or42b show a striking difference in their sensitivity to ethyl acetate, as measured physiologically (). The half-maximal response of Or42b occurs at a dose that is below the threshold of Or42a, suggesting that Or42b is a high-affinity receptor and that Or42a is a low-affinity receptor. Thus Or42b is most informative at low ethyl acetate concentrations and Or42a is most informative at high concentrations.
We next examined the roles of these receptors in driving behavioral response to ethyl acetate at concentrations ranging over three orders of magnitude, a range in which the behavioral response is relatively uniform (). To determine how the behavioral response depends on the activities of the two receptors, we examined mutants of each. Mutants of Or42b, i.e. Or42a+Or42b-, which lack the receptor that is sensitive to low concentrations, show a reduced response at low concentrations, but an approximately normal response at high concentrations (). Mutants of Or42a, i.e. Or42a-Or42b+, give a response index that is approximately normal at low concentrations, but reduced at high concentrations. These results suggest that the behavioral response to low concentrations is driven primarily by Or42b, and the response at high concentrations is dependent on Or42a.
Interestingly, the response of
Or42a-Or42b+ mutants declines and becomes repellent as the concentration increases. One possible interpretation, among others, is that hyperactivation of Or42b, which is very sensitive to ethyl acetate, triggers a repulsion circuit, but that in wild type this circuit is overridden or suppressed by activation of Or42a. In any case, this result is consistent with the existence of non-linearity in the olfactory circuitry and the possibility of interactions among ORNs in the larval antennal lobe via lateral connections, which have been documented in the more complex olfactory system of the adult fly (
Bhandawat et al., 2007;
Kazama and Wilson, 2008;
Olsen et al., 2007;
Olsen and Wilson, 2008;
Shang et al., 2007).
As a control for these experiments, we have shown that mutations of Or42a and Or42b do not have a general effect on odor response: both respond normally to 2-heptanone and propyl acetate ().
In summary, the receptors that have low and high thresholds, respectively, for physiological response to ethyl acetate are the receptors that are required at low and high concentrations, respectively, for behavioral response. Thus two receptors are required in order for the animal to respond strongly to a broad range of ethyl acetate concentrations. We note that the consistency between the results of the physiological analysis and the behavioral genetic analysis provides additional validation of the empty neuron system.
Predicting the behavioral response index from receptor responses
Having shown how the response of two receptors can be integrated to engender the behavioral response to a single odor, we next sought to expand our focus and consider how the responses of the entire receptor repertoire are integrated. As an initial step, we summed the total number of action potentials elicited by each odorant from all 21 receptors of the receptor repertoire, and plotted each sum against the response index obtained for that odorant. Although we had not expected to find a simple relationship, we found a modest but clear correlation between the total spike input and the behavioral response: greater spike input correlated with greater RI
L (
Figure S5A,B). (RI
L is used for statistical rigor; the RI was expressed as a log-odds ratio, i.e. logit transformed to yield RI
L, as described in Experimental Procedures, to best satisfy the normality assumptions of the analysis. A plot of RI v. total spike input is nearly superimposable.) The correlation was observed when we plotted the total number of spikes elicited by a 10
-2 dilution of each odorant (R
2=0.33, p=0.002) or by a 10
-4 dilution (R
2=0.28, p=0.006). We found comparable correlation coefficients for attractants and repellents when considered separately, in the case of both concentrations (not shown).
This relationship encouraged us to ask whether a more refined function of receptor response might provide a more powerful prediction of behavior. We asked whether there exists a set of coefficients {a, b, c….,} which when multiplied by the physiological responses of receptors in spikes/s {R1, R2, R3…,} would yield products that sum to equal the behavioral response, i.e. RIL = aR1 + bR2 + cR3 +…. In other words, we asked whether the behavioral response can be predicted as a weighted integration of input from the respective receptors.
We tested all possible subsets of receptors to identify subsets whose activities together predicted behavior. A stepwise linear regression model was used, along with several criteria for model selection. It should be noted that the modeling approach used here attempts to predict the behavioral outcome dictated by the olfactory circuitry, and does not presume to necessarily reveal the form of the circuitry. In particular, the number of receptors in the model arises from statistical procedures to minimize overfitting the data (see Experimental Procedures).
We were surprised to find that the responses of small subsets of receptors were powerful predictors of behavioral response. Remarkably, 81% of the behavioral variation could be explained by the activity of only five receptors (Or42a, Or45a, Or74a, Or82a, and Or85c), as follows:
RIL= 0.0509 + 0.0061 Or42a + 0.0022 Or45a - 0.0039 Or74a - 0.0113 Or82a + 0.0050 Or85c
We note that while the model above was ranked most highly in our analysis, we identified other models that explained similar levels of variation. However, of the ten models that explained the greatest degree of variation, while controlling for overfitting, all contain Or42a, Or82a and Or85c, and each model is based on six or fewer receptors. The models were selected using the criteria described in Experimental Procedures.
How well does the five-receptor model defined above predict the behavioral response to the odors of the panel? We assessed the predictive power of the model in two ways: first with a drop-out analysis and second with a novel set of odorants. For the drop-out analysis we systematically withdrew the data for each of the 26 odorants, fitted the regression using the responses of the five receptors to the remaining 25 odorants, and then predicted the behavioral response for the withdrawn odorant. A comparison of the predicted and observed behavior is shown for each odorant in . Overall, the model predicted 74% of the behavioral variation when each odorant was withdrawn in turn.
We then tested the model with an entirely different set of odorants. We chose a set of 14 odorants that were not used in the original regression model. They were selected to have a similar representation of functional groups and chain lengths as the original set. We measured the activities they elicited from the five receptors of the model and the behavioral responses they elicited from the larva (
Fig. S6). The model again predicted the behavioral responses well. Specifically, the model predicted 55% of the total behavioral variation, and a single odorant, 1-pentanol, accounted for 46% of the prediction error. If this odorant had been excluded, then the model would have predicted 71% of the behavioral variation among the remaining 13 odorants.
Predicting repellency from receptor responses
Perhaps the most fundamental classification of odor responses is the division into attractive and repellent behavioral responses. Little is known about the mechanistic basis by which an individual odor stimulus elicits attraction v. repulsion. In addition to the intrinsic scientific interest of this problem, the basis of repellency is of great practical interest: the ability to predict repellency could be of value, for example, in the development of new insect repellents.
Among the initial set of 26 odorants examined, six had negative RI values and can be considered as repellents (for statistical analysis, see Experimental Procedures). We wondered whether there were any common features among the physiological responses that these odorants evoked.
We noticed a striking relationship between repellency and inhibitory responses across the repertoire: of the three odors that elicit the greatest number of inhibitory responses from the receptors of the repertoire, all are repellents. The 26 odors of our panel elicited between 0 and 8 inhibitory responses from the 21 receptors when tested at a 10
-2 dilution ( and
S2A). The two odors that elicited 8 inhibitory responses were 2-methyl phenol and fenchone, which are among the strongest repellents; the one odor that elicited 7 inhibitory responses was 1-nonanol, also a repellent ( and Cobb et al., 1992). The one odor that elicited 6 inhibitory responses, moreover, was 4-methyl phenol, one of the weakest attractants. Overall, the mean number of receptors inhibited by the repellents was 5.8, whereas only 2.9 inhibitory responses were elicited by non-repellents (p<0.01, Mann-Whitney test).
Heretofore we have used a stringent definition of inhibition: a reduction in response to <50% of the spontaneous rate. As another test of the relationship between repellency and inhibition, we adopted a less stringent criterion: that the mean response rate is less than the mean spontaneous rate. Using this criterion we found that repellents inhibited 10.0 receptors v. 5.3 receptors for non-repellents (p<0.01, Mann-Whitney test).
While there is a clear difference in the mean number of receptors inhibited by repellents v. non-repellents, there is also overlap in the distributions. We therefore asked whether repellents could be predicted from the responses of the entire receptor repertoire using a multivariate classification procedure. We first used a linear discriminant function that provided a clear separation of repellents from non-repellents (
Figure S7), but with some cross-validation error. We then found that perfect discrimination of the six repellents from the 20 non-repellents could be achieved using a non-parametric kernel density discrimination function (SAS, Proc DISCRIM), with zero cross-validation error. The number of repellents in the data set is low, and thus the results must be interpreted with caution. Moreover, the form of the discriminant function does not necessarily reflect the form of the circuitry driving the repellent response. However, this analysis suggests that a useful statistical prediction of repellent behavior can be developed from the physiological responses of the 21-receptor repertoire.
Odor space is largely conserved between two olfactory systems
Our analysis of the larval odor receptor repertoire allowed us to address the conservation of odor space among olfactory systems. Do two odorants that elicit similar patterns of activity among the receptors of one receptor repertoire elicit similar patterns among the receptors of another?
We examined the 16 odorants that have been tested against both the larval receptors and against adult antennal receptors {
Hallem, 2006 #1491}. We computed both Euclidean distances and angular distances between each pair of odorants within each olfactory system and then compared the distances between the two olfactory systems by matrix correlations. The matrix correlation using angular distances was 0.728 and for Euclidean distances it was 0.559; both correlations were significant at p < 0.0001 by a matrix permutation test (Mantel test) (
Manly, 1991); this test shows that the elements of two matrices containing corresponding distances are significantly correlated). Thus, even though only 6 receptors overlap between the two sets of tested receptors (n=21 larval receptors; n=24 adult receptors), the overall distance relationships between the 16 tested odors, as measured by electrophysiological activity patterns, remain largely consistent between the two receptor repertoires.
To investigate further the similarity between larval and adult odor space, we performed the analysis again after excluding the six receptors that are common to both sets of receptors (i.e. we considered 15 larval receptors and 18 adult receptors). For the angular distance the matrix correlation was 0.675 (p < 0.0001) and for the Euclidean distance the matrix correlation was 0.606 (p < 0.0001). Thus the overall distance relationships between the 16 tested odors, based on these larval receptors and these adult receptors, are also largely consistent.