Introduction

The ability to adapt to sensory cues is key to the survival of many living organisms18. The goal of this sensory adaptation is to condense neural responses to a persisting or recurring stimulus in order to focus on novel cues. Such condensed neural responses are thought to increase system efficiency by reducing the metabolic requirements of information encoding912. However, is sensory information lost when neural responses are condensed during adaptation? While this computational task appears relatively elementary, any solution should satisfy two important requirements. First, attenuation of stimulus-evoked responses due to adaptation should not alter the identity of the stimulus. Second, information regarding the intensity of the stimulus may have to be preserved when odor intensity conveys information of behavioral importance. The latter requirement is paramount when sensory stimuli have different behavioral significance at varying intensities. For example, some odorants are attractive at lower intensities but switch to being aversive beyond certain threshold values1315. Alternately, sensitivity to intensity gradients has to be maintained for sustained periods to allow for navigation towards or away from the stimulus source1618. Such biological contexts necessitate the preservation of odor intensity information across repeated encounters with the odorant. However, the mechanism by which information about stimulus intensity is robustly maintained is not fully understood.

Sensory adaptation occurs over both long- and short-time scales12,1921. Long-term adaptation requires gene transcription and protein translation that occurs on the time scale of hours to days6,2224, whereas short-term adaptation occurs on a millisecond to minute time scale2,11,25. Short-term adaptation results from a wide range of mechanisms with prominent models in the olfactory system being vesicle depletion and facilitation of lateral inhibition2532. While both vesicle depletion and lateral inhibition are likely to occur in the olfactory system, these complementary mechanisms will manifest in two different,but testable, outcomes. Vesicle depletion will impact the most prominent neurons with elevated odorevoked responses, whereas facilitation of lateral inhibition would progressively suppress activity in the relatively weak responders. Irrespective of which mechanisms underlie olfactory adaptation, the neural responses evoked by a stimulus will change over time and across repeated encounters. This raises the question, does adaptation then corrupt or alter information regarding the odorant?

The ambiguity that neural adaptation introduces has been documented in visual33 and auditory systems34,35. However, prior studies have shown that odor identity is robustly maintained by the spatiotemporally patterned neural responses and can allow precise recognition across repeated encounters with the odorant30,3638. But is intensity information corrupted by the changes in neural activity due to sensory adaptation? Both reductions in stimulus intensity and sensory adaptation lead to a reduction in the overall spiking responses in most neural circuits. Do these neural response variations due to neural adaptation and intensity decrements be disentangled? We sought to examine this in the relatively simple insect olfactory system.

Our results show that odor intensity changes result in subtle variations in the overall combination of activated neurons. While sensory adaptation results in an overall reduction in the population neural responses, the activated ensemble is robustly maintained to allow simultaneous and precise decoding of odor identity and intensity. Using an innate behavioral assay, we show that physiological adaptation coincided with enhanced behavioral differentiation between odorant intensities. Thus, our work reveals a simple but elegant approach to achieve adaptation-invariant intensity coding and recognition in the insect olfactory system.

Results

Behavioral responses to an odorant vary with repetition and intensity

We began by examining how behavioral responses change with repetitions of an odorant. We used an innate behavioral assay where the presentation of certain odorants triggered an appetitive palp-opening response (POR) in starved locusts (Fig. 1a; Supp. Fig. 1). The probability of an odor eliciting a POR (P(POR)) was regarded as an indicator for stimulus perception (Fig. 1b). Note that the P(POR) was computed by averaging the behavioral responses across locusts for each trial or repetition number (Fig. 1c, d). Intriguingly, we found that the P(POR) systematically increased with repetitions of the same odorant and stabilized after about five stimulus presentations.

Behavioral response facilitation upon odor repetition.

(a)Schematic of the innate palp-opening response (POR) experimental paradigm. Opening of the locust’s maxillary palps within 15 s of the odor onset was considered a POR.

(b)Schematic of the block odor stimulation protocol. Each block consisted of ten trials and a 4s odor pulse was presented in each trial. A 15-minute no-stimulation period separated the blocks.

(c)Response matrices are shown summarizing individual locust PORs (rows) for blocks of ten trials (columns). White boxes indicate PORs in a specific trial, while black boxes indicate the absence of PORs in that trial. Locusts were sorted such that the least responsive locusts are shown at the top and the most responsive ones are near the bottom. PORs varied between locusts, across trials.

(d)The probability of PORs across locusts (P(POR)) is shown as a function of trial number for two odorants: iaa and bza at high intensities (1% v/v). P(POR)s as a function of trial number is shown for four odorants (oct, hex, iaa, bza) at two different intensities (0.1% v/v – low and 1% v/v – high). The trials were split into early (trials 1-5) and late (trials 6-10) trials. While the P(POR) of both odorant intensities increased over trials, the P(POR) values for high-intensity odor exposures were notably higher than P(PORs) for low-intensity odor presentations.

(f) Bar plots compare the average P(POR) for early and late trials for the low- and the high-intensity presentations of the same odor. P(POR) values observed at high and low intensities of an odorant were similar for early trials (p>0.05, one-tailed t-test). For three out of four odors tested the P(PORs) for high-intensity odor exposures were significantly higher (p<0.05, one-tailed t-test). This suggests that locusts can better perceptually distinguish between high and low odor intensities during later trials.

Previous results have shown that the POR is a function of both stimulus identity and intensity3941. Therefore, we presented each locust with randomized blocks of the same odorant at either high (H, 1% v/v) or low (L, 0.1% v/v) intensities (Fig. 1e). Our results indicate that the P(POR) evoked by high and low stimulus intensities of an odorant were not significantly different in the earlier trials (Fig. 1f; p > 0.05 for hex, iaa, and bza; one-tailed t-test). However, the P(POR) became significantly different for the high and low-intensity exposures of the same odorant during the later trials (p ≤ 0.05 for hex, oct, and iaa; one-tailed t-test). In other words, high and low stimulus intensities were more behaviorally discriminable in the later trials. Hence, these results counterintuitively suggest that behavioral responses to an odorant at two different intensities were different after multiple encounters with the same stimulus.

Odor-evoked response adaptation leads to complex changes in individual neurons

How do neural responses change as a function of stimulus repetition? To understand this, we examined the odor-evoked responses of individual projection neurons (PNs) in the locust antennal lobe as different odorants at two different intensities were repeatedly presented (Fig. 2a). Odor-evoked responses in individual PNs varied subtly across trials (Fig. 2b). In general, most PN responses to the odorant were stronger in the first trial and reduced during later trials of the same stimulus. Notably, the changes in PN responses in two different blocks of trials, where the same odorant was presented, showed repeatable changes in odor-evoked spiking patterns and the overall response strength (Fig. 2c).

Ensemble neural activity systematically changes across repeated encounters with an odorant.

(a)Schematic of the olfactory stimulation protocol. Each block consisted of twenty-five trials with a four seconds odor pulse delivered in each trial. The inter-trial interval was 56 s. Two datasets were collected. Each dataset consisted of five randomized blocks of four odorants (dataset 1: hex (H), hex (L), oct (H), oct (L), hex (H)-repeat; dataset 2: iaa (H), iaa (L), bza (H), bza (L), iaa (H)-repeat). A fifteen-minute no-odor stimulation period separated blocks of trials.

(b)Raster plots of eight representative projection neurons (PNs) in the locust antennal lobe are shown. Spiking activities are shown for twenty-five trials (rows) with earlier trials shown at the top and later trials at the bottom. The shaded region indicates the four seconds odor stimulation period, and the identity of the stimulus is indicated in each plot.

(c)Raster plots are shown for four representative PNs during two blocks of trials. The same odorant was presented in both blocks. Note that spiking activity changes are repeatable across the block of trials.

(d) Schematic showing how vesicular depletion and lateral inhibition facilitation models would change spiking activity in individual neurons. Y-axis represents the change in response over trials ((1st trial response – 25th trial response)/ (1st trial response + 25th trial response)). Positive numbers indicate that the first trial had a stronger response, hence response reduction. Negative numbers indicate that the last trial had a stronger response, hence response facilitation. Along the x-axis, the response in the 1st trial is shown. Vesicular depletion should impact the strongly activated neurons more, whereas lateral inhibition facilitation should progressively suppress weak responders.

(e)Following the schematic in panel d, response changes observed in PNs are shown. Only PNs with significant odor-evoked responses are included in this plot. Responsiveness of PN was determined based on whether their peak odor-evoked response in the first trial exceeded the mean + 6.5 SDs of spontaneous activity preceding odor stimulation in that first trial. In total hex (L/H), n = 26, 39 PNs; oct (L/H), n = 35, 28 PNs; iaa (L/H), n = 39, 35 PNs; bza (L/H), n = 33, 40 PNs were included in the plot. Different colors and symbols are used to denote the identity of the odorant and its intensity (low-intensity trials – triangles and high-intensity trials – circles). The dotted line shows the regression fit revealing a linear relationship, if any, between the activation strength in the first trial (x-axis) and the response reduction observed (y-axis). As can be noted the R2 values are low indicating that both vesicle depletion and lateral inhibition facilitation models fail to capture the adaptation trends in our datasets.

(f)Changes in PN spiking activity over trials during high-intensity odor exposures (y-axis) are plotted against response changes observed for the same PN during low-intensity exposures of the same odorant. The dotted line indicates the regression fit revealing a linear relationship between the response changes observed during high and low-intensity exposures of the same odorant (hex, R2=0.034; oct, R2=0.032; iaa R2=0.067; bza R2=0.0016). The poor linear regression fit here indicates that reductions in neural response amplitude for one odor intensity do not model reductions in response amplitude for another odor intensity.

Next, we examined how individual PNs adapted during stimulus repetitions. Based on existing literature, two widely considered models of adaptation in the olfactory system are vesicular depletion and facilitation of lateral inhibition. While the vesicular depletion model predicts that the strongest firing neurons would adapt the most, in the lateral inhibition facilitation model the weakly activated neurons would progressively be suppressed more (Fig. 2d). To test these two opposing but complementary theories, we calculated the percentage change in the total spike counts observed during first and the 25th odor presentations (Fig. 2e, y-axis; positive values indicate weaker response in the later trial and hence response reduction, whereas negative values indicate facilitation of odor-evoked responses). We plotted these normalized response changes against the total spike counts observed in the very first trial (Fig. 2e, x-axis). The dotted lines represent the best linear model that fit the two datasets we recorded. If vesicular depletion was the predominant form of adaptation, the slope of the regression line would be positive. Whereas a negative slope would indicate that lateral inhibition was the predominant form of adaptation. However, as can be noted, across both datasets with two concentrations of the four odorants tested, neither model captured the adaptation dynamics in our datasets (R2 = 0.05, 0.01, respectively).

Are trial-to-trial changes in stimulus-evoked neural responses at two different concentrations of the same odorant correlated? If this is the case, then it would be reasonable to expect neurons that adapt more at the higher concentration to also change in a similar fashion at a lower concentration and vice versa. To assess this, we calculated and plotted the normalized change in spike counts for each PN at both high and low concentrations of the same odorant (Fig. 2f). Notably, the PNs with the greatest response change during repetitions of the higher-intensity odor exposures were not the ones that adapted heavily when encountering the same stimulus at a lower intensity (R2 = 0.034 (hex), 0.032 (oct), 0.067 (iaa), 0.0016 (bza)). In sum, these results indicate that change due to adaptation in individual neurons is not a simple function of its response strength. Further, how individual neurons adapted is a function of both odor identity and intensity.

Adaptation and intensity decrements both reduce ensemble neural response strength

What is the overall change in the spiking response across the neural ensemble? To examine this, we calculated the average spiking response across all PNs recorded and plotted the spike counts as a function of time for each trial (Fig. 3a). Expectedly, the overall response strength was stronger in the first trial and reduced when the same stimulus was repeatedly encountered. A similar reduction in overall response strengths was also observed when the odorant intensity was reduced (Fig. 3b)

Stimulus repetition and intensity decrements reduce spiking responses.

(a) Peristimulus time histograms (PSTHs) across all PNs (hex and oct, n = 80 PNs; iaa and bza, n = 81 PNs) for each odorant. PSTHs of trials 1 and 25 are shown.

(b) PSTHs across all PNs are compared between high and low-intensity odor exposures. Response during the first trial is shown.

(c) Summed spike counts (across PNs and four-second odor presentation) are calculated and shown as a function of the trial number. The dotted line indicates the spike count of the 25th trial of the high-intensity odor exposure.

Are the response reductions during stimulus repetition comparable to those observed with stimulus intensity decrements? To determine this, we plotted the total spike counts summed across all PNs over the entire stimulus presentation duration and plotted it as a function of the trial number (Fig. 3c). We compared the changes over trials with the total spike counts observed during the high (H) and low (L) intensities of the same odorant. Notably, the total PN ensemble spike counts decreased due to both stimulus intensity reduction and repetition of the stimulus. Do these results indicate that adaptation could potentially confound information regarding stimulus intensity?

Adaptation-invariant encoding of odor intensity

Could the information regarding stimulus intensity be robustly encoded in the population neural responses? As noted earlier, both spiking activities in individual neurons and the volume of spikes at the ensemble level vary across trials. If adaptation alters the overall spiking activity amplitude (i.e. vector length; Fig. 4a), while changes in stimulus intensity alter the combination of activated neurons (i.e. vector direction; Fig. 4b), then information about odor identity and intensity can be preserved in an adaptationinvariant fashion. To test whether the recorded ensemble of neurons in our datasets encoded information in such an adaption-invariant fashion, we first visualized the high-dimensional neural activities in each trial using a dimensionality reduction approach (see Methods). We plotted the ensemble responses in each 50 ms time bin during the odor presentation time window (4 s), and linked them based on the order of their occurrence, to generate trial-by-trial odor response trajectories (Fig. 4c, d). Note that each trial generated a single loop response trajectory after dimensionality reduction. Further, neural responses in different trials evoked closed-loop trajectories that systematically reduced in length with increasing repetition but maintained their direction indicating that the combination of neurons activated were mostly conserved. Different odorants evoked responses that were markedly different in the combination of activated neurons and therefore their response trajectories were oriented in different directions. Notably, different concentrations of the same odorant also generated response trajectories that subtly varied in direction. As a direct consequence, even though the overall spiking activities reduced over trials, the odor response trajectories reliably maintained their direction in all the trials. Hence, information about odor identity and intensity was robustly maintained across all trials.

Odor identity and intensity information is maintained across trials.

(a)A schematic showing how the ensemble neural activity might change between two odorants and across multiple trials or repetition. The combination of neural activated should differ between different odorants, and therefore the odor identity should be represented by population response vectors that differ in their direction. Repetitions should reduce response strength without altering the combination of neurons activated. If this were the case, the later repetitions of the same odorant would evoke a response that can be represented by vectors that maintain directions while progressively becoming lesser in length or magnitude.

(b)A schematic coding scheme for achieving adaptation invariant intensity coding. The combination of neurons activated changes markedly with odor identity and subtly with odor intensity. These responses become less intense without altering the combination of neurons activated. If this were the case, then the ensemble vector direction would change with both odor identity and intensity, and the vector direction robustly maintained even though the vector length continues to change with repetition.

(c)Trial-by-trial odor-evoked ensemble PN response trajectories after dimensionality-reduction are shown for hex and oct. Only high-concentration exposures of both odorants are included in this plot. For each odor, the temporal response trajectories for the 1st, 5th, 10th, 15th, 20th, and 25th trials are shown (color gradient from light to dark as trial number increases). Note that the odor-response trajectories become increasingly smaller during later trials, but the direction of the trajectories is robustly maintained.

(d)Trial-by-trial response trajectories for low concentrations exposures of the hex and oct are shown. Similar convention as panel c. (f)Similar plot as in panel c, but now comparing the odor-evoked response trajectories elicited by hex at high and low intensities.

(g)Similar plot as in panel d but comparing the response trajectories elicited by oct at high and low intensities.

To quantify these observations, we performed a correlation analysis using high-dimensional PN response vectors. If a similar combination of PNs (high-dimensional activity vectors) were activated in two different trials, then their correlation will be high. Correlations across different trials and concentrations of the same odorant are shown in Fig. 5a. Confirming results from the dimensionality reduction analysis (Fig. 4), we found odor-evoked responses across different trials of an odorant at a particular intensity were highly correlated (diagonal blocks on each panel). The similarity between odorevoked responses evoked by different concentrations of the same odorant was considerably less correlated (off-diagonal blocks). A hierarchical clustering analysis revealed that the odor-evoked responses nicely clustered based on odor identity and then by intensity (Fig. 5b)

High and low stimulus intensities activate distinct ensembles.

(a)Correlations between neural responses observed in different trials are shown. Each pixel/matrix-element represents a similarity between mean neural responses in one trial versus those in another trial. Diagonal blocks reveal the correlation between trials when the same odorant at a specific intensity was repeatedly presented.

(b)A dendrogram was generated using a correlation distance metric comparing trial-by-trial ensemble spiking activities evoked by two different stimuli at two different intensities (see Methods). Two major response clusters that correspond to stimulus identity and intensity were identified. The number at the leaf node represents the trial number.

In sum, our results reveal that a combinatorial code could encode information regarding odor identity and intensity in an adaptation-invariant fashion. Since the same information about a stimulus is represented with fewer spikes in the later trials, we conclude that adaptation refines the odor codes by making them more efficient.

Neural response suppression vs. behavioral output facilitation

Neural responses evoked by an odorant reduced with repetition (Fig. 3). However, as we noted earlier, the behavioral responses increased upon stimulus repetition (Fig. 1). This negative correlation is evident when the neural responses (total spiking activities evoked across neurons) was plotted against the behavioral POR response probability during a particular trial (Fig. 6). This negative correlation between neural suppression and behavioral enhancement was observed for all odorants and both concentrations examined. It is worth noting that the last trial of the high-intensity exposure evoked neural responses that were comparable in strength to the first trial of the low-intensity exposure of the same stimulus. Yet, the behavioral responses observed during these two trials were markedly different. Hence, our results also indicate that the probability of the POR response is not simply a function of the total number of spikes elicited.

Neural adaptation inversely correlates with behavioral facilitation.

(a) The probability of odor-evoked POR (P(POR)) for a given trial is plotted against the total spike counts elicited during the trial or repetition number. Therefore each point represents a single trial. Symbols with light and dark colors are used for differential trials of high and low-concentration odor exposures. The line represents the regression fit between the behavioral and neural responses (hex (L), R2 = 0.65; hex (H), R2 = 0.79). The negative slope of the regression line indicates that while the neural responses diminish over trials, the behavioral responses increase over trials.

(b-d) Similar plots as in panel a but for oct (L/H), iaa (L/H), and bza (L/H).

Neural and behavioral changes due to adaptation are odor specific

Finally, we wondered if the neural response reduction and behavioral enhancements were a global, non-specific state change in the olfactory system brought about by the repetition of any odorant, or are the observed neural and behavioral response changes odor specific. To examine this, we used a ‘catch-trial paradigm.’ (Fig. 7a) In these sets of experiments, we repeatedly presented an odorant to induce adaptation. After a substantial reduction in the neural response or stabilization of the behavioral response, a deviant or catch stimulus was presented. Persistence of the observed reduction in neural responses and enhancement in behavioral responses when the catch stimulus is introduced would reveal that the adaptation-induced changes are not odor specific.

Adaptation results in odor-specific neural response reductions and behavioral output facilitation.

(a)Two blocks of trials were used. First, a block of 30 trials where one odorant (hex) was presented in all trials except the 26th trial (the catch trial). During the catch trial, a deviant stimulus (iaa or app) was presented. After a 15-minute no-odor reset period, a second block of ten trials of the deviant stimulus was presented (trials 31 – 40). This was done to determine the unadapted (first trial) and adapted (later trials) responses of the same set of PNs to the stimulus used in the catch trial.

(b)Summed spike counts (during 4 s odor presentation period) across all PNs were calculated and shown as a function of trials (hex repetitions). Repeated presentations of odor A resulted in a substantial reduction in odor-evoked spike counts. However, the presentation of the deviant odor (iaa) during the catch trial (trial # 26) resulted in a marked increase in spike counts. The response strength was similar to the non-adapted responses (Trial # 31) for iaa (H). (right) Similar results are shown when hex was repeatedly presented and app was used as the deviant stimulus. Note that the response to app during the catch trial did not recover as was noted for iaa (left).

(c)Correlations between odor-evoked responses with hex response in the first trial are shown as a function of trial number.

(d)A similar catch-trial paradigm was followed for behavioral experiments. Each locust was presented with a repeating odorant for eight trials. A deviant stimulus was presented in the ninth trial. This was followed by six more repetitions of the recurring stimulus.

(e)Response matrix summarizing POR responses of individual locust PORs are shown. Same convention as in Fig. 1. (f)The P(POR) across locusts was calculated and plotted as a function of the trial number. Notably, when the deviant stimulus was presented there is a marked decrease in P(POR). The dotted line indicates the unadapted P(POR) value for the deviant stimulus.

We used the catch-trial experiments with two different catch stimuli. Our neural recordings reveal that when iaa (ester) was used as the catch stimulus after repeated presentations of hex (alcohol), the neural responses strength recovered and increased to an overall higher-spike count (Fig. 7b). The observed catchiaa response strength matched the unadapted response level to iaa when it was separately presented following a 15-minute reset window (i.e. Fig. 7b, trial# 31). In contrast, when apple (a complex mixture) was used as the catch stimulus, the neural responses did not recover to higher levels indicating significant cross-adaptation (Fig. 7b). Nevertheless, in both cases, the ensemble responses evoked by the catch stimulus were highly similar to the unadapted responses it usually generates (Fig. 7c). Therefore, these results indicate that there is stimulus specificity to the adaptation-induced changes in neural responses and that cross-adaptation does not corrupt the identity of the deviant stimulus.

Behaviorally, we chose odor pairs where the component odorants had markedly different levels of POR responses (Fig. 7d). Repeated presentations of an odorant (bza or iaa) increased the P(POR) across the locusts as observed before (Fig. 7e). Notably when the catch stimulus was presented, the behavioral 1responses dropped close to the first-trial P(POR) elicited by that odorant (Fig. 7f; dotted line). These findings were also observed when the recurring and catch odorants were switched (Fig 7f).

In sum, these neural and behavioral results indicate that the observed neural and behavioral response changes brought about by the repeated encounters with an odorant are at least partially odorspecific and are not global non-specific changes in the olfactory system.

Discussion

Sensory systems often have to deal with constraints that are of an opposing nature. For example, the ability to maintain stable representations is key to robust recognition of sensory stimuli4244. Often variations in the external environment, such as humidity and temperature, may introduce changes in neural responses that have to be filtered out to preserve the identity of the stimuli. Additionally, the strength or intensity of the stimulus can also change widely. Compressing or even removing these variations might be necessary to recognize an odorant independent of its intensity4247. Alternately, in certain cases, it might also become important to maintain some of these variations so they can be exploited for guiding or altering the behavioral response. As an example, it might be worth noting that many chemicals that are pleasant at lower intensities can become repelling or even harmful to the organism at extremely high intensities1315. Hence, the behavioral responses should accordingly vary with intensity. Further, many organisms have shown that they can follow an odorant to their source (source localization). Such behaviors are of ecological importance as they allow the organisms to find the food source or potential mates. In all these cases, it becomes important to maintain variations with stimulus intensity to ensure the organisms are able to approach the odor source. In our earlier studies, we have focused on examining how robust odor recognition is achieved37,40,48. Here we examined the latter problem of how to maintain information about the stimulus intensity while the responses in the neural circuits are becoming weaker due to adaptation.

Consistent with earlier studies30,49,50, our data revealed that the response of individual projection neurons in the antennal lobe varied with the repetition of the same stimulus. While a reduction in odorevoked responses was observed in many neurons, an increase in spiking activity and changes in temporal firing patterns was also observed. Notably, these changes were not random but highly repeatable (Fig. 2c). Hence, these results indicate that the effect of adaptation in a neural circuit is highly consistent and reliable.

In addition to spike rate reduction, oscillatory synchronization, and inter-neuronal coherence build up in the locust antennal lobe over repeated encounters with an odorant30,51,52. This short-term memory was not observed in the antenna but endured in the antennal lobe for several minutes after the termination of the stimulus. Similar results were also reported in the mouse olfactory bulb43. The inputs from sensory neurons were shown to be consistent across repetition whereas the outputs of the olfactory bulb were systematically reduced with stimulus repetition. These results suggest that loci of short-term adaptation, at least when the stimulus is discontinuous and recurring, is in the antennal lobe/olfactory bulb circuitry. Although the precise mechanism for achieving the same is not fully understood49,53.

What is the computational significance of this neural adaptation? Sensory adaptation has been implicated in several important computations such as high-pass filtering, matching neural response to stimulus statistics, generating sparser codes, and optimal representations33,35,5456. In the olfactory system, adaptation has also been suggested to stabilize neural representation such that recognition of odor identity from ensemble neural responses improved after repeated encounters with the stimulus30. However, adaptation also introduces a potential confound whereby the absolute neural response strength becomes an unreliable indicator of stimulus strength. Our results extend these earlier studies by revealing how this potential confound could be resolved in the olfactory system.

Increases in stimulus intensity have also been shown to result in higher response amplitudes in firing neurons and in the recruitment of additional activated neurons35,42,44,57. While our results show that this is indeed true, we also found that several neurons responded preferentially during the lower concentration exposures but did not respond when the higher intensity of the same stimulus was presented (Supp. Fig. 2). Hence, the activated combination of neurons differed both with stimulus identity and intensity. Furthermore, even though adaptation reduced the overall spiking activities across neurons, the combination of neurons activated was robustly maintained. In other words, the population response vector direction (determined by the combination of neurons activated) was consistent. Only the vector length (determined by the total spikes across neurons) is reduced with the repetition of the stimulus. This simple approach was sufficient to encode the information about the odor intensity in an adaptation-invariant fashion. Additionally, since equivalent information about stimulus identity and intensity was preserved with fewer spikes, neural adaptation can be considered as a method of optimizing the neural representation2,3,9,10,33,34.

Does the observed persistence of concentration information in the neural activities after adaptation translate to observable differences in behavioral outcomes? Using an innate palp-opening response (POR) we examined how behavioral responses varied with stimulus repetition. Unlike the neural responses, we found that the PORs increased with repetition and after about five trials settled into odor-specific levels of responses. Notably, the PORs in the later trials were significantly higher for higher stimulus intensities for most of the odorants used. In sum, our results showed that locusts are better able to behaviorally differentiate stimulus intensity after adaptation.

Alternately, adaptation due to the repetition of an odorant could alter the perceived intensity of the recurring odorant or other related test odorants. The compromise of stimulus intensity information due to adaptation has been observed in behavioral discrimination tasks in rats and even in psychophysical studies in humans58,59. Such a result can arise when the combinatorial response patterns evoked by the adapting and the test stimulus have significant overlap. However, our results indicate that when changes in odor intensities alter the combination of neurons activated, it is feasible to maintain information about odor intensity levels in a robust fashion.

Finally, we note that the anticorrelated link between neural suppression and behavioral enhancements (called ‘repetition priming’) is well established in other animal models60. How can weaker neural responses lead to enhanced behavioral outcomes? Our results indicate that the antennal lobe neural responses in later trials were not sharpened by pruning out weaker responses (Supp. Fig. 2). The weakening of the response strength could potentially lead to a sparser response in the downstream mushroom body. In addition, as we noted earlier, neural response synchrony increases with stimulus repetition. Whether weaker but synchronized activity in neurons drives sparser and more selective responses which then produce enhanced behavioral outcomes remains to be carefully examined.

Acknowledgements

We thank members of the Raman Lab (Washington University in St. Louis) for feedback on the manuscript. We thank Pearl Olsen for insect care. This research was supported by NSF (1707221, 1724218, 2021795) and ONR (N00014-19-1-2049, N00014-21-1-2343) grants to B.R. and an Imaging Sciences Fellowship to D.L.

Author contributions

BR conceived the study and designed the experiments/analyses. DL performed the behavioral experiments and associated analysis. DS performed all the electrophysiology experiments and collected the data. AC collected the electrophysiology data for the catch-trial paradigm. DL and LZ performed all the neural data analysis. DL and BR wrote the paper taking inputs from all the authors. BR supervised all aspects of the work.

Methods

Odor stimulation

For delivering odorants, we followed a protocol described in our earlier work. Briefly, odorants were diluted in mineral oil to either 1% or 0.1% concentration (v/v) and sealed in glass bottles (60 ml) with an air inlet and outlet. A pneumatic picopump (WPI Inc., PV-820) was used to displace a constant volume (0.1 L/min) of the static headspace above the diluted odor-mineral oil mixture into a desiccated carrier air stream (0.75 L/min) directed toward one of the locust’s antennae. A vacuum funnel placed behind the locust preparation continuously removed the delivered odors.

The first set of experiments included multiple blocks, with twenty-five trials each when one odorant at one intensity was repeatedly presented. Each trial in the block included a four-second stimulus presentation window. The inter-stimulus interval (between trials) was 60 seconds. Different odorants at different intensities were repeatedly presented in different blocks. A 15-minute window, when no stimulus was presented, separated two consecutive blocks of trials. This window was included to reset any short-term memory that may have formed due to repeated presentation of the same stimulus30.

The second set of experiments involved two blocks of trials. The first block of trials included 30 trials and the second block consisted of ten trials. In the first 25 trials of block 1, hexanol at 1% v/v was repeatedly presented. In the 26th trial, a puff of either isoamyl acetate 1% v/v or apple 1% v/v (a ‘deviant’ odorant) was presented. In trials 27-30, hexanol at 1% v/v was again presented. This was followed by a 15-minute reset window when no stimulus was presented. In the second block of ten trials (trials 31-40), the deviant stimulus was repeatedly presented. All inter-stimulus intervals (between trials) were 60 seconds and all odor presentations were four seconds long.

Electrophysiology

Post-fifth instar adult locusts (Schistocerca americana) were reared in a crowded colony with a 12-hour - 12 hour light-dark cycle. both males and females were used for electrophysiological experiments. First, the locusts were immobilized with both antennae intact. Then the primary olfactory region of their brain, the antennal lobes, were exposed, desheathed, and perfused with room-temperature locust saline. Extracellular multiunit recordings of projection neurons (PNs) were performed with a 16-channel, 4x4 silicon probe (NeuroNexus) that was superficially inserted in the antennal lobe (AL). Prior to each experiment, all probes were electroplated with gold to achieve impedances in the range of 200 to 300 kΩ. The recordings were acquired with a custom 16-channel amplifier (Biology Electronics Ship; Caltech, Pasadena, CA). The signals were amplified with a 10k gain, bandpass filtered (0.3 to 6 kHz), and sampled at 15 kHz using a LabView data acquisition system. A visual demonstration of this protocol is available online61.

PN spike sorting

To obtain single-unit PN responses, spike sorting was performed offline using the four best recording channels and conservative statistical principles62. Spikes belonging to single PNs were identified as described in earlier work37,63. The following criteria were used to identify single units: cluster separation >5x noise standard deviations, number of cluster spikes within 20ms < 6.5% of total spikes, and spike waveform variance < 6.5x noise standard deviations. In total, 161 PNs from 40 locusts were identified. Two datasets were collected. In the first dataset responses of 80 PNs to hex and iaa were recorded, and in the second dataset responses of 81 PNs to iaa and bza were monitored.

Correlation analysis

The PN spikes were binned in 50 ms non-overlapping time bins, and spike counts of different PNs were concatenated to obtain a population spike count vector. Pearson correlation coefficients between two PN ensemble spike count vectors were calculated using Equation 1.

Here, Xtrial i and Xtrial j are time-averaged high-dimensional activity in trials i and j, respectively. σtrial i and σtrial j are the standard deviations of X rial i and Xrial j, respectively.

Each pixel/matrix-element in the correlation plot shown in Fig. 5a indicates the similarity between PN spike count vectors observed in the ith and jth trials.

Tensor-based data decomposition

We first organized neural response data as a three-way array (Neuron × Time × Trials; the stimulus information was also blended into trial dimension) and then employed a direct 3-way tensor decomposition approach64,65. Here, the 3-d data cube was approximated using three loading matrices, A, B, and C with elements ai f(neuron dimension), bjf (time dimension), and ckf (trial dimension). ei j k was the residual element (see the equation below). The tri-linear model was found using alternating least squares.

where i, j, k denote the three different dimensions, and F indicates the total number of factors used for the analysis that was determined by the core consistency diagnostics64,65. In our case, when F = 3, the core consistency was above 50 %, while it dropped to below 40% when F = 4. Therefore, we used three factors for our data decomposition.

Trial-to-trial odor trajectory

For this analysis, we first reconstructed the dataset by computing the outer product of the loading matrices that were obtained by the tensor decomposition. The reconstructed 3-d tensor was then unfolded into a concatenated matrix (i.e. along the trial dimension). After unfolding, the ensemble projection neuron responses were arranged as time series data of n dimensions (where n is the number of neurons) and m steps (the number of 50 ms time bins × the number of trials). Note that only the projection neuron activities during the four-second stimulus presentation window in each trial were used for this analysis. The ensemble projection neuron response vectors (in a given 50 ms time bin) were projected onto the three eigenvectors of the response covariance matrix that accounted for the most variance in the dataset, using principal component analysis. Finally, the low-dimensional points were connected in a temporal order to visualize neural response trajectories to different stimuli on a trial-to-trial basis. All trajectory plots shown in Fig. 4 were generated after smoothing with a 3-point running average low-pass filter.

Neural response similarity and dendrogram analysis

First, we calculated the summed spike counts during the 4 s odor presentation window for each projection neuron. Then the correlation similarity between two spike count profiles across projection neurons was calculated using Equation 1. Similarly, in this analysis, xj and x j represent a n × 1 vector (n = 80 for hex-2oct; n = 81 for iaa-bza) for trial i and j, respectively. The dendrogram was generated by hierarchical clustering of all stimulus identities, intensities, and individual trials based on the correlation distance. The dendrogram was created in such a way that the furthest pairwise distance between any two samples assigned to an individual cluster was minimized.

Behavior experiments

Experiments were performed on post-fifth instar locusts of either sex that were starved for approximately 24 hours. All behavioral experiments occurred between 10 am and 2 pm. The protocols for the innate behavioral preference experiments and palp tracking used in this study were published in previous studies37,41,48,63. The odor delivery setup and stimulus sequences for the single odorants were similar to the ones described for the electrophysiological experiments (Figs. 1a, 2a). Any palp movement that occurred within 15 seconds of the odor stimulus onset was considered a palp-opening response (POR).

POR responses of locusts across trials and odor intensities.

(a) POR responses of each locust (rows) across trials (columns) are shown for both high and low-intensity hexanol exposures. Black boxes indicate no POR responses, and white boxes indicate POR responses in that trial. The probability of palp-opening responses across locusts is shown as a function of trial number. The average POR in the first five and last five trials is compared between high and low-intensity hexanol exposures (p>0.05, one-tailed t-test).

(b-d) Similar plots characterizing locust POR responses for the other three odorants.

Combinatorial response patterns across two intensities of the same odorant.

Comparison of odor-evoked responses between high and low intensities exposures of all four odorants used in the study. Spiking activities of individual PNs were summed over the entire odor presentation window, sorted based on their response to the high-intensity presentations, and plotted to reveal the combinatorial response. Mean ± S.D. across the twenty-five trials is shown.