Insects find food and mates by navigating odorant plumes that can be highly intermittent, with intensities and durations that vary rapidly over orders of magnitude. Much is known about olfactory responses to pulses and steps, but it remains unclear how olfactory receptor neurons (ORNs) detect the intensity and timing of natural stimuli, where the absence of scale in the signal makes detection a formidable olfactory task. By stimulating Drosophila ORNs in vivo with naturalistic and Gaussian stimuli, we show that ORNs adapt to stimulus mean and variance, and that adaptation and saturation contribute to naturalistic sensing. Mean-dependent gain control followed the Weber-Fechner relation and occurred primarily at odor transduction, while variance-dependent gain control occurred at both transduction and spiking. Transduction and spike generation possessed complementary kinetic properties, that together preserved the timing of odorant encounters in ORN spiking, regardless of intensity. Such scale-invariance could be critical during odor plume navigation.https://doi.org/10.7554/eLife.27670.001
Insects follow odor trails carried by the wind to find mates and sources of food. The turbulent motion of the air means that these odors tend to arrive in whiffs with varying intensities and durations, which makes it difficult to distinguish them. Insects use sensory cells called olfactory receptor neurons on their antenna to process odors. Specialized receptor proteins on the surface of these olfactory receptor neurons detect odor molecules and set off a cascade of events in these cells that ends with a signal being sent to the brain.
Much is known about how insects detect and process different kinds of smells, but it remains less clear how their olfactory neurons process the timing and intensity of odor whiffs. Now, Gorur-Shandilya et al. report what happens in the olfactory receptor neurons of fruit flies when they have to compensate for variations in the duration and intensity of odor whiffs.
In the experiments, fruit flies were exposed to two sweet-smelling odors. To do so, Gorur-Shandilya et al. built an apparatus that enabled them to control the airflow with enough precision that they could simulate the variability in the timing and intensity of natural odors in the air. The response of the flies’ olfactory receptor neurons to these smells was recorded. The experiments showed that the neurons could adapt to both the average intensity and the variance in intensity of odor signals.
The ability of these neurons to adapt to the average intensity of the odors followed a specific pattern, which is also seen in sensory cells responsible for vision and touch. Adapting to the average strength of an odor slows down the first of two steps in its processing. However, the second step has a complementary mechanism to speed up signals to the brain, so the timing of an odor whiff is accurately captured regardless of how strong it is. Based on these results, Gorur-Shandilya et al. created a biophysical model that could reproduce the experimental data, including the slowdown in the first step.
The experiments and the model may now help other scientists to investigate how different animals detect and process smells. For example, some insects are pests of agricultural crops, while other insects, such as mosquitos, spread diseases between people. A better understanding of how insects detect odors may help scientists to find ways to interfere with these processes to protect food crops and reduce the spread of tropical diseases.https://doi.org/10.7554/eLife.27670.002
Insects navigate odor landscapes that are often not smooth gradients (Cardé and Willis, 2008; Riffell et al., 2008). Instead, turbulent airflows shape odor plumes into intermittent whiffs separated by stochastic durations of background air (blanks). In the absence of reliable spatial gradients, navigating insects may combine the timing of whiff encounters (Vergassola et al., 2007) with sensation of wind direction (Budick and Dickinson, 2006; Cardé and Willis, 2008; Duistermars et al., 2009) to navigate odor plumes towards mates and food. Insect olfactory systems face dual challenges in detecting natural odor plumes. First, the intensity of whiffs is typically distributed according to a power law (Murlis et al., 1992), with intense whiffs interleaved unpredictably with weak ones (Riffell et al., 2008). Second, whiff durations and blank durations are also distributed as a power law over a wide range of time scales (Celani et al., 2014). The encoding problem is aggravated by shifting local statistics of odor encounters, which change with wind speed (Nagel and Wilson, 2016, 2011), position (Justus et al., 2002), or environment (Murlis et al., 2000). How does the olfactory system manage to encode whiffs of odors whose intensities and timing can vary over such wide ranges?
Several features of the olfactory system contribute to encoding odor stimuli of different intensities. A single odorant can be detected by multiple receptor types, with different sensitivities(Hallem and Carlson, 2006). Static compressive nonlinearities at both olfactory receptor neurons (ORNs) and their post-synaptic targets, the projection neurons (PNs), selectively amplify weak signals and suppress responses to large signals (Bhandawat et al., 2007; de Bruyne et al., 2001). Glomerular mechanisms implement a type of divisive gain control that maintains PN sensitivity within the range of changing ORN responses (Bhandawat et al., 2007; Luo et al., 2010; Olsen et al., 2010; Olsen and Wilson, 2008). Finally, transduction currents in response to odor pulses scale inversely with the intensity of the background signal, consistent with the Weber-Fechner Law (Cafaro, 2016; Cao et al., 2016). However, whether ORN firing follows a similar scaling is unclear (Cafaro, 2016; Martelli et al., 2013). Thus, although it is known that the input-output curve of ORNs to odor stimuli changes with odor background, how ORN gain (from stimulus to firing rate) scales with background signal intensity has not been characterized.
Olfactory responses in insects can be fast. Transduction can be initiated within milliseconds of odor reaching the antenna (Szyszka et al., 2014). The speed of the response is enhanced by ORN spike generation, which emphasizes changes in transduction currents (Nagel and Wilson, 2011), and by PNs (Kim et al., 2015), which maintain fast information transmission from ORNs to PNs (Jeanne and Wilson, 2015; Nagel et al., 2015; Raccuglia et al., 2016). In contrast, adaptation to high intensity stimuli slows down transduction (Cao et al., 2016; Kaissling et al., 1987; Nagel and Wilson, 2011), a property that might make it difficult to reliably encode the timing of odor encounters.
We investigated in vivo how Drosophila ORNs encode encounters with naturalistic odor plumes. To address this question, we first developed an odorant delivery system that reproducibly delivered odorants with naturalistic or Gaussian statistics with controlled means and variances. We simultaneously recorded the odorant stimulus using a fast photo-ionization detector (PID), and recorded extracellularly from identified ORNs.
We found that ORNs encoded broadly distributed naturalistic signals by using two mechansism: front-end nonlinearities that are inherent in receptor binding to ligand, as well as two adaptation mechanisms that are sensitive to the mean and variance of the stimulus. These adaptive mechanisms caused ORNs to rapidly desensitize following encounters with odorant whiffs, dynamically adjusting gain while responding to intermittent odorant stimuli. ORNs adapted to changes in the mean stimulus at the level of transduction by scaling gain inversely with the stimulus intensity, consistent with the Weber-Fechner Law. Variance-dependent gain control took place at both signal transduction and spiking machinery.
While the transduction response time slowed down with increasing stimulus intensity, the spiking machinery sped up to compensate. These complementary kinetic changes caused the firing rate response time to remain invariant with stimulus intensity. This reveals a mechanism that could allow ORNs to preserve information about the precise timing of odor encounters over a wide range of rapidly changing signal intensities. A minimal two-state model of the activity of the Or-Orco complex (olfactory receptor and co-receptor Orco) with an adaptation mechanism that feeds back onto the free energy difference between active and inactive conformations reproduced Weber-Fechner scaling, slowdown of signal transduction kinetics, and responses to intermittent and Gaussian stimuli.
Odorant signals used to study ORN adaptation typically consist of long pulses or constant backgrounds of various intensities (Cafaro, 2016; Cao et al., 2016; Martelli et al., 2013; Nagel and Wilson, 2011). However, airborne stimuli encountered by flying insects can be intermittent with both the intensities of encounters and durations between encounters broadly distributed as power laws (Celani et al., 2014). Since ORN transduction can be adapted by odorant pulses as brief as 35 ms on timescales as fast as 500 ms (Cao et al., 2016), we asked to what extent the gain of ORNs could change dynamically during responses to naturalistic stimuli, amplifying responses to isolated whiffs of odorant, and suppressing responses to whiffs following dense clumps of whiffs.
We measured the responses of ab3A and ab2A ORNs to naturalistic stimuli of ethyl acetate and 2-butanone. These odorants elicit spikes in these neurons (Hallem and Carlson, 2006), and are easy to control and measure (Martelli et al., 2013) (Figure 1—figure supplement 1). We used in vivo extracellular recording to record both the local field potential (LFP) and spikes from a single sensillum, with simultaneous measurement of the stimulus (Figure 1—figure supplement 2). Previous results have shown that: LFP responses are unaffected by the addition of TTX, which eliminates neural spiking, suggesting that LFP signals were generated upstream of the spiking machinery; and that LFP signals are unaffected when the neuron’s partner cell in the sensillum is genetically ablated, when that partner does not sense the odorant, suggesting that LFP signals are generated by the neuron of interest (Nagel and Wilson, 2011). Though the LFP could reflect activity of nearby sensilla, it serves as an imperfect but useful proxy for transduction activity in ORNs (Johnston et al., 1995; Kaissling, 1986; Nagel and Wilson, 2011; Su et al., 2012).
The naturalistic stimulus we used was intermittent and consisted of brief odor whiffs of varied amplitude (Figure 1a–b). Durations of whiffs and blanks were broadly distributed, with a power law of exponent −3/2 to match natural intermittent statistics of odor plumes (Celani et al., 2014) (Figure 1—figure supplement 3). ab2A and ab3A ORNs responded to whiffs with transient decreases in the local field potential (LFP) and corresponding increases in the firing rate (Figure 1a–b).
Even though individual whiff intensities were broadly distributed (First line in Figure 1a) ORN responses to these whiffs were more even, so that responses to faint whiffs were amplified more than those to intense whiffs. To quantify these differences, we defined the gain of the neuron to be the change in the response for a unit change in the stimulus (). Since ORNs do not respond instantaneously to odorant stimuli (Cao et al., 2016; de Bruyne et al., 2001; Martelli et al., 2013; Nagel and Wilson, 2011) we fit linear filters to best predict the LFP and firing rate from the stimulus. We used these filters to make linear predictions of the responses from the stimuli (Figure 1—figure supplement 4). Changes in gain were therefore defined as deviations from the linear prediction of response from the stimulus, similar to (Baccus and Meister, 2002; Kim and Rieke, 2001). We visualized these gain changes by plotting the LFP responses against linear prediction of the LFP (Figure 1c) and the firing rate against the linear prediction of the firing rate (Figure 1d). Each excursion in these plots corresponds to the ORN’s response to a single whiff. Excursions occurred with different slopes, suggesting that ORN gain changed frequently in time. Deviations from linearity persisted even when filters computed from Gaussian inputs were used to project the stimulus, suggesting that the existence of these deviations do not depend on the exact shape of the filter, but rather reflect a property of the ORN response not captured by the linear model (Figure 1—figure supplement 4).
Variations in the gain () clearly do not arise solely from a static output nonlinearity, such as one associated with a linear-nonlinear transformation, since plotting neuron response against the projected stimulus (Figure 1c–d) does not yield a single transformative function. (Dayan and Abbott, 2001). We reasoned that changes in the gain could arise from input nonlinearities due to odor-receptor binding and channel opening. To visualize the nonlinearity between the stimulus and response, we plotted LFP and firing rate responses to each whiff in the naturalistic stimulus as a function of the amplitude of that whiff (Figure 1e–f). A clear sigmoidal dependency is visible in the plot of LFP responses against whiff intensity, consistent with a front-end nonlinearity arising from receptor-odorant binding. However, in both the LFP and firing rates, responses to whiffs with similar intensities varied significantly, deviating from a single sigmoidal dose-response curve (Figure 1e–f).
What causes these deviations from the dose-response curve? One possibility is that these deviations are due to random variability in the responses of the neuron. Another possibility is that these deviations are due to adaptation of the neuron to the stimulus history preceding each whiff, which may vary with every whiff. To distinguish between these possibilities, we collected whiffs that had similar amplitudes, and examined the LFP and firing rate responses they evoked (Figure 2a–b). The amplitude of LFP and firing rate responses elicited by these whiffs varied inversely with the amplitude of the preceding stimulus: whiffs that occurred in isolation (purple) elicited the largest responses, while whiffs that followed earlier, large whiffs (blue, red) elicited the smallest responses, suggesting that ORN responses can be modulated by stimulus history.
To quantify this context-dependent modulation, we estimated deviations of the LFP and firing rate response to each whiff from the median response (see Materials and methods). Deviations in LFP response to each whiff decreased with mean stimulus in the preceding 300 ms (Figure 2c, , Spearman test), and were uncorrelated with the amplitude of the whiff that elicited them (Figure 2c, inset, , Spearman test). Similarly, deviations in the firing rate responses to each whiff decreased with mean stimulus in the preceding 300 ms (Figure 2d. , Spearman test), and were uncorrelated with the amplitude of the whiff that elicited them (Figure 2d, inset, , Spearman test).
To generalize beyond a particular timescale of the stimulus history, we parametrized the stimulus history of each whiff by the amplitude and time since the preceding whiff, and grouped estimated deviations from the median response into positive or negative (Figure 2e–f). When response deviations were negative (smaller than median responses, blue dots), the amplitude of the preceding whiffs tended to be larger, and the time since the previous whiff tended to be shorter. When response deviations were positive (red dots), the amplitude of preceding whiffs tended to be smaller ( is 0.52 V when deviations are positive vs. 1.09 when deviations are negative for firing rate responses and 1.04 V vs. 1.06 V for LFP responses), and the time since the previous whiff tended to be longer ( is 3025 ms vs. 361 ms for firing rate responses, and 2556 ms vs. 566 ms for LFP responses).
What causes this context-dependent suppression of responses following preceding whiffs? One possibility is a bi-lobed filter, with one positive and one negative lobe, followed by a rectifying nonlinearity. Such a filter is partly differentiating, and has been measured in linear models of the firing rate (Kim et al., 2011; Martelli et al., 2013; Nagel and Wilson, 2011) and would lead to attenuated responses to the second of two closely spaced whiffs due to linear superposition. Such a mechanism may partly account for context dependent variation in firing rates. However, stimulus-to-LFP filters, computed for this stimulus and others, are mono-lobed (Figure 2—figure supplement 1), and appear purely integrating (Nagel and Wilson, 2011), ruling out contributions to dynamic modulation of LFP responses by this mechanism. A model with a static front-end nonlinearity and a mono-lobed filter fit to the LFP also cannot reproduce context-dependent adaptation observed in the LFP (Figure 2—figure supplement 1), suggesting that this context-dependent variation in response arises at least in part from ORNs dynamically varying gain in response to naturalistic stimuli.
Since the mean and variance of naturalistic stimuli are correlated over many timescales, (Figure 1—figure supplement 3), it is unclear whether adaptation in this context is sensitive to the mean or the variance (or to some other statistic) of preceding whiffs. To determine how changing one moment of the stimulus distribution changed ORN gain, and to disambiguate the effect of receptor saturation from adaptation, we proceeded to other experiments using Gaussian stimuli with changing means (Figure 3) and variances (Figure 4).
A common strategy used by sensory systems to encode signals over a broad range of background intensities is to scale the response according to the Weber-Fechner law (Fechner, 1860; Stevens, 1957; Weber, 1834), that is, to control gain inversely with stimulus mean. In ORNs, transduction currents elicited by odorant pulses are reduced by preceding pulses (Nagel and Wilson, 2011) and scale inversely with background intensity, consistent with the Weber-Fechner law (Cafaro, 2016; Cao et al., 2016). It remains unclear whether the ORNs’ ultimate output — the firing rate — follows the same Weber-Fechner scaling (Cafaro, 2016; Martelli et al., 2013). We therefore stimulated ab3A ORNs with a set of fluctuating ethyl acetate stimuli with increasing means (Figure 3a) but roughly constant variances (Figure 3b, Figure 2—figure supplement 1a). ORNs responded to the stimulus with the smallest mean by modulating firing rates between 0–60 Hz (Figure 3c). This response range progressively decreased with increasing mean stimulus intensity (Figure 3d), though the mean response remained at ~30 Hz. To estimate ORN input-output curves, we plotted ORN responses against the stimulus projected through the normalized best-fit linear filter for each stimulus, estimated by least-squares fitting (Baccus and Meister, 2002; Chichilnisky, 2001; de Boer and Kuyper, 1968; Rieke et al., 1997) (Figure 3e) Input-output curves grew shallower with increasing mean stimulus. We defined the ORN gain to be the slope of the input-output curve at that mean stimulus, similar to (Baccus and Meister, 2002).
ORN gain in each trial varied with the mean stimulus in that trial as an approximate power law with exponent –1 (Figure 3f). ORN gains could also be estimated by the ratio of standard deviation of the response to the standard deviation of the stimulus. This measure yielded similar values of ORN gain, and also decreased as a power law with exponent −1 (Figure 3—figure supplement 1b). This exponent is consistent with the Weber-Fechner Law, which postulates that the just noticeable difference between two stimuli is inversely proportional to the stimulus magnitude (Stevens, 1957). Rescaling the projected stimulus by the gain predicted by Weber’s Law collapsed all input-output curves onto a single curve (Figure 3g).
Can front-end or back-end nonlinearities reproduce the observed change of input (stimulus)-output (firing rate) curves (Figure 3e)? Clearly, no single output nonlinearity can fit the data shown in (Figure 3e), since a single function cannot fit all the input-output curves. Since a front-end nonlinearity is present (Figure 1), we asked whether a static nonlinear-linear (NL) model could reproduce this data, with the input nonlinearity parameterized by a Hill function where S represents the input, and K the half maximum value. (Figure 3—figure supplement 2a–c). NL model responses increased with mean stimulus (Figure 3—figure supplement 2c), unlike in the data (Figure 3d–e). However, if the half maximum value of the Hill function was allowed to vary with the mean stimulus, the model could qualitatively reproduce the data, suggesting adaptation at the front-end nonlinearity (Figure 3—figure supplement 2d–f).
To determine if similar gain-control relative to mean signal intensity was broadly observed, we tested additional ORNs from the two major olfactory organs of the fly, the antenna and the maxillary palp (ab2A, pb1A), and used ecologically relevant odorants from three different functional groups (ketones: 2-butanone, alcohols: 1-pentanol, esters: isoamyl acetate) in various combinations. In all five cases, the neurons decreased gain with increasing odorant concentration, and obeyed a roughly inverse scaling (Figure 3—figure supplement 1c–f). Rescaling the projected stimulus by the Weber-Fechner relation collapsed all ORN responses onto a single curve, similar to Figure 3g (Figure 3—figure supplement 3). Thus in vivo, for various neurons and odorants, ORN firing rate followed the Weber-Fechner Law.
In other sensory modalities, such as vision, some peripheral neurons adapt not only to the mean but also to the variance of the signal (Baccus and Meister, 2002; Rieke, 2001). We therefore asked whether ORNs adjust their gain in response to changes in the variance of the signal. We stimulated ab3A ORNs with fluctuating ethyl acetate stimuli in which the variance of the signal changed every 5 s (Figure 4a), switching back and forth between high to low values, around a nearly constant mean (Figure 4b; Figure 4—figure supplement 1a), a protocol used to study gain control in visual neurons (Baccus and Meister, 2002; Fairhall et al., 2006; Shapley and Victor, 1978; Rieke, 2001; Smirnakis et al., 1997).
As expected, ORNs responded to input fluctuations by modulating their firing rate. Interestingly, ORN firing rate variance did not vary as much as the stimulus variance between epochs of high and low stimulus variances, suggesting that ORNs actively changed their gain to compensate for such input differences (Figure 4c–d). ORN input-output curves during high variance epochs (red) were shallower than during low variance (blue) epochs (Figure 4e). Trial-wise ORN gain decreased with the variance of the stimulus (Figure 4f). ORN gains estimated by dividing the standard deviation of the response by the standard deviation of the stimulus showed a similar decrease in ORN gain with stimulus variance (Figure 4—figure supplement 1b).
A simple coding strategy maximizes a neuron’s information capacity by matching its input-output curve to the cumulative distribution function (c.d.f) of the stimulus (Laughlin, 1981). Like the c.d.f.s (dashed lines), the input-output curves (solid) are steeper during the low variance epoch. On a trial-by-trial basis, ORN gain was correlated with the c.d.f slope (r2 = 0.7) (Figure 4g). However, as the input variance changed by a factor of 2.5, the gain in the neuron only changed by a factor of 1.7, not as much as would be required for optimal information encoding. In these experiments, the gain changed within ~130 ms following the change in stimulus variance (Figure 4h).
ORN responses arise through two sequential steps: odor transduction followed by spike generation (Nagel and Wilson, 2011). Does each step possess separate gain control mechanisms, or is gain control achieved solely at one step? Previous studies place the mechanism of adaptation to mean stimulus at the level of signal transduction (Cafaro, 2016; Cao et al., 2016; Nagel and Wilson, 2011). How the spiking machinery might influence gain control, and where adaptation to signal variance takes place, remain unknown.
We reanalyzed the responses of ab3A to ethyl acetate signals (Figure 2 and 3) and measured ‘transduction gain’ (stimulus to LFP) and ‘firing gain’ (LFP to firing rate). Changing the stimulus mean alone changed gain in LFP (Figure 5a–b, Figure 5—figure supplement 1). However, gain at the spiking machinery was largely invariant to the ten-fold change in the mean stimulus (p=0.41, Spearman rank correlation), with a 1 mV change in LFP leading to a ~ 10 Hz change in the firing rate, consistent with earlier studies (Nagel and Wilson, 2011) (Figure 5c–d). Transduction gain, like ORN gain, scaled with the Weber-Fechner Law, for a variety of odor-receptor combinations (Figure 5—figure supplement 1a-d) consistent with previous studies (Cafaro, 2016; Cao et al., 2016). In contrast, adaptation to the stimulus variance changed gain both at transduction and at spiking (Figure 5e–h). Both gains changed by a factor of ~1.3 from the high to the low variance epoch (p<0.001, Wilcoxon signed rank test), contributing roughly equally to variance gain control (Figure 5—figure supplement 1f).
Stimulating ORNs with odorants evokes responses from both transduction and spiking machinery, making it hard to determine how independent gain control at the two modules are. To further pin-point the contributions of signal transduction and firing machineries to gain control, we expressed Chrimson channels (Klapoetke et al., 2014) in ab3A ORNs and activated them using red light, either in isolation or in combination with odorants.
First, we used a fluctuating ethyl acetate stimulus to probe transduction and ORN gain while increasing the neuron’s firing rate using increasing backgrounds of red light (Figure 6a–b). While increasing light levels elicited increasing firing rates (Figure 6b inset), ORN and transduction gain did not vary with the intensity of supplemental light (Figure 6a–b). This suggests that constitutive spiking activity does not feed back onto LFP adaptation, or overall ORN gain.
Second, we used a fluctuating light stimulus to probe the spiking gain while stimulating the ORN and its receptors with increasing backgrounds of ethyl acetate odorant (Figure 6c–d). While ethyl acetate backgrounds of increasing intensity increased ORN firing rate (Figure 6d inset), they failed to change gain in the spiking machinery. Increasing odor backgrounds moved input-output curves along the y-axis (Figure 6c), consistent with increasing firing due to background odor, but failed to change the slope of these curves, suggesting that ORN gain to the fluctuating light probe was not changed. This suggests that adaptation at transduction does not affect gain of the spiking machinery, consistent with our result that increasing odor backgrounds decreased gain at transduction, but not spiking (Figure 5a–d). Thus, Weber-Fechner scaling in ORN gain control to stimulus mean is insulated from activity of the spiking machinery.
Variance gain control exists in a wide range of neurons (Baccus and Meister, 2002; Díaz-Quesada and Maravall, 2008; Nagel and Doupe, 2006; Rieke, 2001; Wark et al., 2009; Zaghloul et al., 2005) and in models of spiking neurons (Gaudry and Reinagel, 2007; Hong et al., 2007; Yu and Lee, 2003; Yu et al., 2005), suggesting that variance gain control could be an intrinsic property of spiking neurons. To determine if the spiking machinery alone could give rise to variance gain control, we stimulated ab3A ORNs that express Chrimson with fluctuating light stimuli of different variances at fixed mean. ORN input-output curves were steeper when the variance of the light stimulation was smaller (Figure 6e–f), similar to the curves observed with odor stimulation (cf. Figure 4e). We observed that gain changed by a factor of ~1.5 when the standard deviation of the light stimulus changed by a factor of ~3, consistent with variance gain control occurring partly in the spiking machinery (Figure 5e–h), though Chrimson channels might exhibit their own nonlinear activation properties.
When navigating odor plumes, the precise timing of the encounter with the plume carries important information, which may be lost if adaptation changes the lag between signal and response in a concentration-dependent manner. The kinetics of ORN spiking in response to pulses of odorant are invariant to the pulse intensity and to the background intensity over a range of odorant concentrations (Martelli et al., 2013). Paradoxically, adaptation to background odorants slows transduction current responses to odor pulses (Cao et al., 2016; Kaissling et al., 1987; Nagel and Wilson, 2011). We hypothesized that these seemingly contradictory results might be resolved if the ORN spiking machinery speeds up to compensate for the intensity-dependent slowdown in the LFP.
We characterized responses to odorant stimuli on increasing backgrounds by measuring both ORN spike rates and LFPs. We computed cross correlation functions between the stimulus and the LFP for various stimulus backgrounds. For stimuli on low backgrounds, LFP cross-correlation functions peaked earlier, while for stimuli on larger backgrounds, LFP cross-correlation functions peaked later (Figure 7a), consistent with previous results (Cao et al., 2016; Kaissling et al., 1987; Nagel and Wilson, 2011). Surprisingly, cross-correlation functions from the stimulus to the firing rate were similar between stimuli on low and high backgrounds (Figure 7b), consistent with (Martelli et al., 2013). This selective change in the kinetics of the LFP, but not the firing rate, occurred even though there was no change in the stimulus autocorrelation function from low to high stimulus (Figure 7c). We defined the LFP and firing rate lags relative to the stimulus by the location of the peak of the cross-correlation function, and found that while LFP response lags increased with increasing odorant concentration (p<10−2, Spearman test), firing rate lags remained relatively invariant with odorant concentration (p>0.1, Spearman test) (Figure 7d–g).
For firing rate lags to remain invariant with odorant concentration despite a slowdown in transduction, the kinetics of the spiking machinery need to speed up with increasing input to the cell. To test if this is the case, we stimulated ab3A ORNs expressing Chrimson with Gaussian red light stimuli with increasing means, and measured lags between the applied light stimulus and firing rate. Firing lags decreased with increasing light power concentration (p<10−4, Spearman test), suggesting that the ORN spiking machinery can speed up with increasing mean input currents (Figure 7h), as would accompany increasing odor backgrounds (Cao et al., 2016).
If adaptation to the mean slows down transduction, which is compensated for at spiking, does adaptation to the stimulus variance also lead to similar compensatory kinetics? We found that a three-fold change in the stimulus variance, despite leading to changes in LFP and firing gains (Figures 4–5), did not significantly change kinetics either at LFP or firing rate (Figure 7—figure supplement 1), consistent with our earlier results suggesting that mean and variance gain control have distinct mechanisms.
How do adaptive mechanisms at transduction preserve both the Weber-Fechner Law and lead to response slowdowns? In the following we show that a minimal two-state model of the olfactory receptor-olfactory co-receptor (Or-Orco) complex with an adaptation architecture similar to that of the bacterial chemotaxis system (Asahina et al., 2009; Barkai and Leibler, 1997; Emonet and Cluzel, 2008; Shimizu et al., 2010) can reproduce the LFP responses to naturalistic and Gaussian stimuli, as well as Weber-Fechner Law and its accompanying response slow down.
In our model, Or-Orco complexes can be active or inactive (C and C* in Figure 8a) and the active complex binds odorant S with higher affinity than the inactive complex. We assume that ligand (un)binding is fast compared to (in)activation rates ( and in Figure 8b). The fraction of active complexes therefore obeys the equation
where the rates of activation and inactivation are nonlinear functions of the odor concentration S and of the free energy difference between the unbound active and inactive states (Equations 3-4 in Materials and methods). The LFP is modeled as a linear filter acting on the activity (Figure 8b, Equation 5 in Materials and methods). At steady state, Equation (1) reduces to , where the bar indicates steady state. is a monotonically increasing function of the odor concentration S. Increasing the free energy difference shifts this function towards higher values of S, therefore reducing the sensitivity of the system. We model adaptation by assuming that activity of the Or-Orco complex controls the activity of factors that act on the complex to modify the free energy difference :
where is the rate of adaptation. Importantly, the rate of change of only depends on the activity a but not on the free energy difference . The architecture of this feedback is similar to that of the bacterial chemotaxis system and ensures that for increasing values of S, the changes in compensate for changes in free energy due to ligand binding (Barkai and Leibler, 1997). Thus, adaptation eventually returns a to the adapted value providing Weber-Fechner scaling (Shimizu et al., 2010) (as in Figure 3). We assume that the free energy of the complex can only be changed within a finite range, and that the lower bound is reached for small values of S. Thus, in the absence of ligand, the steady state activity can be smaller than . For non-zero values of S, the steady state activity first increases with background signal intensity (Martelli et al., 2013), before it becomes independent of background intensity once it reaches (Figure 8—figure supplement 1), as seen in Figure 3c.
An important intrinsic property of this model is that adaptation to increasing background of odorant decreases the rates of activation and inactivation, and , of the Or-Orco complex , providing a self-consistent explanation for the slowdown of the response kinetics of the LFP upon adaptation. It is interesting to note that this kinetic property emerges because (1) the switching rates are decreasing functions of the free energy difference and (2) the requirement of Weber-Fechner scaling, which causes the adapted value of to scale with the logarithm of the mean signal intensity (see Materials and methods).
The resulting model (Equations 1–2, and Equations 3–5 in Materials and methods) contains six parameters plus another three for converting the signal from activity to LFP. We fit this model to LFP responses to the Gaussian and naturalistic stimuli. The model decreased gain with the mean stimulus background, consistent with Weber-Fechner Law (Figure 8c), and predicted the observed decrease in the LFP gains well (Figure 8d, ). In addition, response lags of this model with respect to the stimulus increased with the mean stimulus (Figure 8e), similar to the slowdown observed in the LFP responses (cf. Figure 7). Finally, this model can also reproduce LFP responses to naturalistic, intermittent signals, approximating well the time trace (Figure 8f,h) and the dependence on previous whiffs (Figure 8g, compare to Figure 2).
Since the spiking machinery compensates for the slowdown in LFP responses to preserve the timing of odorant encounters, we wondered if a simplification of this model that ignores the slowdown of the LFP kinetics upon adaptation could be used to predict firing rate: where N is a static nonlinearity, is a partially derivative-taking linear filter ( indicates convolution), and the steady state solution of Equation (1) with obeying Equation 2. This simplification reduces this model to a type of adaptive nonlinear-linear-nonlinear (NLN) model, which preserves Weber-Fechner Law and reproduces the firing rates of ORN in response to both naturalistic and Gaussian Stimuli (Figure 8—figure supplement 2). Thus it could be a useful tool in modeling ORN responses received by PNs, or in constructing computational models of the antennal lobe (Assisi et al., 2011; Bazhenov et al., 2001; Berck et al., 2016; Capurro et al., 2012; Chong et al., 2012; Hopfield, 1991; Kee et al., 2015; Koulakov et al., 2007; Luo et al., 2010; Sanda et al., 2016; Satoh et al., 2010; Stevens, 2015). The model reproduced the change in the input-output curves on increasing the mean stimulus (Figure 8—figure supplement 2a, cf. Figure 3e) and decreased gain inversely with the mean stimulus, consistent with the Weber-Fechner Law (Figure 8—figure supplement 2b). The model reproduced the observed decrease in the ORN gains (Figure 8—figure supplement 2c, ), and responses to naturalistic stimuli (Figure 8—figure supplement 2d–f).
We examined how ORNs encode naturalistic odor signals and characterized how ORN gain is dynamically modulated in response to stimuli. By using precisely controlled, repeatable odorant stimuli, and linear modelling, we found that: (1) ORN gain varies dynamically during responses to naturalistic stimuli, suppressing responses to whiffs following earlier whiffs (Figures 1 and 2). (2) Gain varies inversely with the mean stimulus (Weber-Fechner Law) and decreases with increasing stimulus variance (Figures 2 and 3). (3) Variance gain control was distributed across transduction and spiking, while mean gain control occurs only at transduction (Figures 5 and 6). (4) While gain control slows stimulus-to-transduction kinetics with increasing background intensity, this is compensated for by a corresponding speed up of transduction-to-spiking kinetics, which maintains the stimulus-to-firing rate kinetics relatively independent of stimulus intensity (Figure 7). A two-state model of the Or-Orco complex activity supplemented with an adaptive architecture similar to that of the classic bacterial chemotaxis system is sufficient to reproduce several key features of LFP response, firing rate response, and gain control (Figure 8).
The Weber-Fechner Law has been observed in several sensory systems, including vision (Burkhardt, 1994; Laughlin and Hardie, 1978; Nikonov et al., 2006), audition (Riesz, 1928), and somatosensation (Holway and Pratt, 1936). In olfaction, the Weber-Fechner Law was demonstrated at the LFP level (Cafaro, 2016; Cao et al., 2016). Here we directly measured ORN firing rate and stimulus intensity and found that the ORN firing rate exhibited Weber-Fechner gain scaling relative to the mean stimulus intensity for five different odor-receptor combinations (Figure 3, Figure 3—figure supplement 1). These data suggest that olfaction shares the Weber-Fechner Law with other sensory systems.
What is the purpose of front-end Weber-Fechner gain scaling? ORNs are capable of spiking up to ~300 Hz (Hallem and Carlson, 2006); however, we found that with their compressive gain control, ORNs maintained firing rates between 0–50 Hz to fluctuating odor stimuli, even with a ten-fold increase in the mean stimulus. The ORNs’ postsynaptic partners, the projection neurons (PNs)(Olsen et al., 2010), are most sensitive to ORN firing rates of below ~50 Hz (Jeanne and Wilson, 2015). Thus, gain scaling at ORNs could act to maintain ORN firing rates in the range that PNs are most sensitive to for a wide range of concentrations, avoiding saturation of the ORN-PN synapse.
In principle, gain control in sensory systems could be affected by several moments of the stimulus distribution, measured over many timescales. In the visual system, gain control depends on stimulus mean and variance, and some studies have shown little dependence on higher moments like the skew and the kurtosis (Bonin et al., 2006; Tkačik et al., 2014). Cell-intrinsic variance gain control exists in a variety of systems, including the retina (Beaudoin and Manookin, 2008; Zaghloul et al., 2005), lateral geniculate nucleus (Lesica et al., 2007), auditory neurons(Nagel and Doupe, 2006), and cortex (Díaz-Quesada and Maravall, 2008; Ringach and Malone, 2007). Photoreceptors do not exhibit variance gain control, and variance adaptation arises only in the subsequent processing in bipolar cells and ganglion cells (Baccus and Meister, 2002; Kim and Rieke, 2001; Rieke, 2001).
What could be the functional role of variance gain control in olfaction? One possibility is to help ensure that ORN responses occupy a large fraction of their dynamic range (Laughlin, 1981). While we quantified our stimulus in terms of the first and second moments of the stimulus statistics in this study, these moments may not map simply onto the salient features that are most relevant to the fly’s encoding scheme. Variance gain control could therefore be a consequence of an adaptive representation that is important to the coding properties of the ORN, but remains unknown to us. Nonetheless, because variance gain control is distributed between transduction and spiking machinery (Figure 5), and mean gain control slows down transduction (Figure 7) (Cao et al., 2016; Nagel and Wilson, 2011) but variance gain control does not (Figure 7—figure supplement 1), adaptation to the stimulus variance is mechanistically distinct from adaptation to stimulus mean.
The results presented here, and the models that reproduce them, focus on the phenomenology of gain control in ORNs. Do these phenomenological results constrain possible mechanisms that could implement gain control in ORNs? Weber-Fechner gain scaling (Figure 3) can be reproduced by models using feed-forward loops (Clark et al., 2013; Goentoro et al., 2009), integral feedback (Yi et al., 2000), or both (Schulze et al., 2015). A detailed biophysical model of odor-receptor binding and channel opening has been proposed to account for transduction responses to odors (Nagel and Wilson, 2011). While this model can change gain via a negative feedback mechanism, it does not reproduce the Weber-Fechner law, or the slowdown of LFP kinetics upon adaptation. Here we showed that these features emerge if we assume that: (i) the activity a of the Or-Orco complex feeds back onto the free energy difference between the active and inactive state of the unbound Or-Orco complex, which in turns affects both the rates of activation and inactivation of the complex; and (ii) the rate at which is modified only depends on the activity a. In summary, the steady state activity in our model depends nonlinearly on the stimulus, reproducing the effects of saturation in responses to naturalistic stimuli. Adaptation shifts the effective half-maximum of the input nonlinearity to the right, recapitulating Weber-Fechner gain control; and decreases transition rates from active to inactive receptor complexes, reproducing slowing LFP responses with adaptation.
Such an architecture reproduces the Weber-Fechner law and is similar to that of the bacterial chemotaxis system (Barkai and Leibler, 1997). There, adaptation is mediated by two antagonistic factors, one that acts on inactive complexes only, and another one that acts on active complexes (Barkai and Leibler, 1997). While the molecular architecture of the signaling pathway in ORNs has not been fully characterized, several studies have implicated calcium as a slow diffusible factor that could mediate adaptation to the mean stimulus (Deshpande et al., 2000; Störtkuhl et al., 1999). Decreasing extracellular calcium levels, or internal free calcium, breaks Weber-Fechner gain scaling at transduction (Cao et al., 2016). Other mechanisms have also been implicated in adaptation of ORNs, like autoregulation of Orco via cAMP signaling (Getahun et al., 2013). While slower adaptive processes also exist, our data on responses to naturalistic stimuli (Figures 1 and 2), and data from paired-pulse experiments (Cao et al., 2016), suggest that some adaptation mechanisms act on fast timescales of several hundred milliseconds.
Many models that decrease gain with increasing mean stimulus also speed up response kinetics (Clark et al., 2013; De Palo et al., 2012; Nagel and Wilson, 2011; Schulze et al., 2015; Seung, 2003), describing well the phenomenology of other sensory systems where gain and response speed trade off (Baylor and Hodgkin, 1974; Dunn et al., 2007; Nagel and Doupe, 2006; Payne and Howard, 1981). However, in olfactory systems, transduction kinetics slow down with increasing stimulus background, both in insect ORNs (Figure 7, (Cao et al., 2016; Nagel and Wilson, 2011)) and in vertebrate ORNs (Reisert and Matthews, 1999). It is not trivial to devise systems in which kinetics slow down with increasing stimulus background. To exhibit this property, a system must increase its effective timescale of response with stimulus intensity, for example, by decreasing all reaction rates uniformly. Interestingly, our model also exhibits a slowdown in the LFP kinetics upon adaptation. This feature emerges intrinsically from the model architecture because: (i) the feedback of the activity onto the free energy difference affects both the activation and deactivation rates of the complex ( and ); and (ii) the Weber-Fechner gain control causes to scale logarithmically with the stimulus, which in turn causes and to decrease.
Earlier work modelled the transformation from LFP to firing rates using a derivative-taking kernel (Nagel and Wilson, 2011). Here, we show that the temporal structure of these kernels depends on the adaptation state of the ORN, and must take derivatives on shorter timescales at higher stimuli to compensate for slowing transduction kinetics. Consistent with this, we see that the latency of spiking decreases increasing optogenetic drive (Figure 7h). Under these conditions, the change in spiking latency is smaller than the increase in transduction lags we observe (Figure 7d–g). While the mechanism of the speed up in spiking with increasing odor stimulus is not known, the neuron’s ability to spike with shorter latencies relative to transduction could depend on the adapted state of its receptors, the level of intracellular calcium, or the distance of its membrane potential from firing thresholds.
What cellular mechanisms could give rise to gain control that is variance dependent? (Figure 4). We found that both the transduction machinery and the spiking machinery of the ORNs exhibit variance-sensitive gain-control (Figures 5 and 6). Variance gain control after transduction could arise from the spike generating machinery. Hodgkin-Huxley (HH) model neurons exhibit variance gain control (Hong et al., 2008; Lundstrom et al., 2008; Yu and Lee, 2003). Simpler neuron models, like the FitzHugh-Nagamo model (Hong et al., 2007), and the linear integrate-and-fire (LIF) model (Yu and Lee, 2003) also exhibit variance-dependent gain control. In the visual system, non-spiking bipolar neurons show variance gain control (Baccus and Meister, 2002; Rieke, 2001) so mechanisms for variance gain control in the absence of spike generation might be similar between these systems.
Previous studies of olfactory adaptation employed conditioning and probe stimuli (Cafaro, 2016; Martelli et al., 2013; Nagel and Wilson, 2011), which typically adapt neurons over many seconds or minutes before testing response properties with a short probe. Other studies using paired pulse protocols (Cao et al., 2016) found that responses to brief pulses of odorant reduced gain on timescales as brief as 500 ms, which is close to the timescale of the neural response to odors (Kim et al., 2011; Martelli et al., 2013; Nagel and Wilson, 2011). Similar fast timescales of gain control have been observed in the visual system (Burns et al., 2002; Baylor and Hodgkin, 1974). We found that this fast gain control was employed by ORNs to dynamically control gain during responses to naturalistic odorant stimuli (Figures 1–2).
Dynamic gain control allows ORNs to respond to the rapidly changing statistics of natural odor plumes, letting gain decrease quickly in response to a large whiff and then ramp up again to a subsequent small whiff. Dynamic inhibition in the antennal lobe (Nagel et al., 2015; Raccuglia et al., 2016) would permit PNs to remain sensitive to these rapid changes in ORN firing rate, ensuring propagation of information about odor encounters to the brain.
Insects follow odor plumes to their source to find food or reproductive mates (Murlis et al., 1992). For flies, this task is challenging since they fly fast (~30 cm/s) (Tammero and Dickinson, 2002) and odor filaments are narrow (Murlis et al., 1992). Even for relatively broad and static odor plumes, flies are within odor plumes so briefly that they experience plume contact and plume loss in quick succession (10–250 ms) (van Breugel and Dickinson, 2014). Olfactory search behavior in this setting consists of rapid flight surges on encountering odor plumes, and stereotyped crosswind casts on losing odor plumes (van Breugel and Dickinson, 2014). Navigation based on odor intensities alone may not be possible, as odor intensities are not informative about the direction to the odor source at length scales longer than 10 cm (Murlis et al., 1992). Indeed, there is a growing body of evidence underlining the importance of timing in olfaction (Martelli et al., 2013; Rebello et al., 2014; Shusterman et al., 2011; Smear et al., 2013, 2011).
In this context, it may be important for the fly to know precisely when it encountered an odor filament. Previous studies have shown that kinetics of transduction slowed during adaptation (Cao et al., 2016; Nagel and Wilson, 2011), but kinetics of firing rate did not (Martelli et al., 2013). Here we reproduced both findings and resolved this apparent contradiction. We discovered that the spiking machinery speeds the kinetics back up. These complementary kinetic mechanisms mean that the timing of short odorant encounters is preserved in neural encoding, regardless of intensity. Such an encoding scheme could aid insects in navigating odor plumes to their source.
When a system responds identically, in amplitude and in kinetics, to stimuli that are different only in scale, the system is said to show fold change detection (FCD)(Goentoro et al., 2009). FCD thus implies the Weber-Fechner law, but systems can obey the Weber-Fechner Law without showing FCD. Another requirement for FCD is that response kinetics remain invariant with respect to the mean stimulus intensity. Thus, ORN responses are intriguingly similar to the response phenomenology of FCD networks (Goentoro and Kirschner, 2009).
Interestingly, olfactory adaptation is linked to flight in insects. Olfactory receptors (ORs) adapt and have co-evolved with flight (Edwards, 1997; Getahun et al., 2013; Jones et al., 2005), and occur only in flying insects (Missbach et al., 2014). In contrast, the more ancient ionotropic receptors (Missbach et al., 2014), found in all insects, do not appear to adapt to prolonged odor stimuli (Cao et al., 2016). While ORs play an important role in larval olfactory navigation (Hernandez-Nunez et al., 2015; Mathew et al., 2013; Schulze et al., 2015; Gepner et al., 2015), the statistics of odor signals close to surfaces, and in the air, where flying insects encounter them, may be very different (Murlis et al., 1992) (Martelli et al., 2013). Receptors capable of fast adaptation may allow flying insects to detect brief whiffs of airborne odors.
Single sensillum recordings from Drosophila antennae were performed as described previously (de Bruyne et al., 2001; Martelli et al., 2013). The recording electrode was inserted into a sensillum on the antenna of an immobilized Drosophila melanogaster and a reference electrode was placed in the eye. Electrical signals were amplified using an Iso-DAM amplifier (World Precision Instruments). The ab3 sensillum was identified by (1) its size and location on the antenna (2) test pulses of 2-heptanone, to which the B neuron is very sensitive, (3) spike shapes (A spikes are larger than B spikes) and (4) spontaneous firing (ab3B fires at a higher rate than ab3A). Other sensilla were identified using test odors to which either the A or B neuron strongly responded to.
All sensilla recorded from in this study contained two neurons (de Bruyne et al., 1999; Song et al., 2012). Generally, spikes from the ‘A’ neuron are larger than spikes from the ‘B’ neuron. However, spike amplitude and spike shape changed in our experiments with strong odor or light drive, due to a phenomenon called ‘pinching’ (Olsson et al., 2006), and due to small movements of the recording electrode relative to the sensillum. To identify spikes from the A neuron under these challenging conditions, we developed a spike-sorting software package written in MATLAB (Mathworks, Inc.), available at https://github.com/emonetlab/spikesort. A copy is archived at https://github.com/elifesciences-publications/spikesort. This package uses the full spike shape, with various dimensionality reduction and clustering methods to reliably identify spikes from noise, and to sort identified spikes. This spike-sorting package performed with 99.5% accuracy compared to manually sorted data on a test dataset.
The local field potential was recorded by lowering the gain of the amplifier and switching to DC mode, where we recorded the sensillar potential without any low frequency cutoff. Spike detection and sorting was reliable in either mode. Since we were only interested in the deflections of LFP in response to a fluctuating odor stimulus, we band-passed the raw voltage in software to remove spikes and slow fluctuations.
Flies were reared at 25°C on conventional fly medium (Helfand and Carlson, 1989). All experiments were performed on adult female flies 3–5 days post-eclosion. Unless otherwise mentioned, recordings were from ab3A neurons in Canton-S flies. In Figure 6 and Figure 7h, we recorded from ab3A ORNs in w; Or22a-GAL4/+; UAS-Chrimson/+ flies. In these flies, only ab3A ORNs were sensitive to light, while ab3B neurons and nearby ab2 sensilla were not.
We used a Photo-Ionization Detector (PID) (200B, Aurora Scientific) to measure the odor stimulus during every experiment. Stimulus measurements occurred simultaneously with all electrophysiology, and the tip of the PID probe was <1 cm of the odor delivery tube and the fly (Figure 1—figure supplement 1). The PID was calibrated by depleting known volumes of pure odorants, and the response of the PID was found to be approximately linear with odorant flux (Figure 1—figure supplement 1). However, due to gradual changes in the sensitivity of the PID detector, odor intensity measurements are not comparable across experiments.
We measured the intensity of red light that we used to activate Chrimson at the location of the fly using a PM160 light power meter (Thorlabs). We used this to construct a function mapping control signals to our LED to light power in µW, and transformed control signals into light power using this function.
Odorants in gas phase were delivered to the antenna by blowing air over pure monomolecular odorants in liquid phase. The flow rate of air over the liquid odorant determined the gas phase concentration.
Mass Flow Controllers (MFCs) (Aalborg instruments and Controls, Inc. and Alicat Scientific) were used to regulate airflows. Dynamic response parameters of Alicat MFCs were chosen either for high speed, and driven with switching times of up to 20 ms (at the cost of reproducibility) or were chosen for high precision, and driven with switching times of 100 ms (at the cost of very fast stimulus control). An odorized airstream (0–200 mL/min) was fed into a main airstream (2 L/min) that was delivered through a glass tube positioned within 10 mm of the fly’s antenna. The secondary airstream passed through a scintillation vial with a machined plastic screw-top lid containing pure odorant (Figure 1—figure supplement 1). Using pure odorant and gas phase dilution permitted excellent reproducibility of the odor stimulus. All tubing was made of chemical-resistant PTFE tubing (McMaster Carr, stock #5239K24).
By varying the control signals to the MFC bank, steps, pulses, and frozen noise waveforms with arbitrary distributions could be reliably delivered. We wrote a general-purpose acquisition and control system called kontroller (available at https://github.com/emonetlab/kontroller) in MATLAB (Mathworks, Inc.) to control MFCs, valves and LEDs and to collect data from electrophysiology and the stimulus measurement.
To generate naturalistic odor stimuli, we randomly varied flow rates over 0–200 mL/min, and used a small solenoid valve (Lee Co.) to deliver 50 ms whiffs of odorant. We used the same frozen random sequence in subsequent trials.
To generate approximately Gaussian-distributed stimuli that differed in their mean, but with similar variances, we started with the ansatz that air flow rates were proportional to measured gas-phase stimulus. We then defined target Gaussian distributions that differed only in their mean, and a parametric distribution from which we drew control signals to the MFC. Using kontroller to automate the process, we performed a direct search on hardware to find the best distribution of control signals that was closest to the desired Gaussian distribution. Further rounds of off-line numerical optimization using nonparametric models of the delivery system ensured that the resultant stimulus distributions were as close to Gaussians and with variances as similar to one another as possible.
We used two MFCs driven by control signals with different variances to generate two Gaussian-distributed stimuli with differing variances. Solenoid valves (Lee Co.) were used to switch from one airstream to the other every 5 s. Control signals were iteratively optimized using simulations and kontroller till the mean stimulus intensity from the two odor lines was indistinguishable.
Since the intensity distribution in our naturalistic odor stimulus was very broad, we defined whiffs of odor as short excursions of the odorant signal above the noise floor. Blanks were defined as the periods of time between whiffs. Whiff intensities were broadly distributed, and were fit with a functional form proposed in Celani et al. (2014). Both whiff and blank duration distributions were fit with a power law with exponent −3/2, following theoretical calculations for a jet flow (Celani et al., 2014) (Figure 1—figure supplement 3). Mean and variance of naturalistic odor stimuli were computed in non-overlapping windows of length (Figure 1—figure supplement 3e, = 400 ms). Window lengths were varied from = 10 ms to = 10 s (Figure 1—figure supplement 3f).
For every whiff shown in Figure 1e–f, we computed the median response for all whiffs in a bin centered around that whiff’s stimulus amplitude, that encompassed other whiffs if their amplitude was within 10%. The fractional deviation in response to a given whiff i is defined as , and is a dimensionless number that is negative when responses are smaller than the median (). To determine if the time to the previous whiff and the amplitude of the previous whiff of negative deviations and positive deviations were different, we used a 2-sample two-dimensional Kolmogorov-Smirnoff test, based on the method proposed by Peacock (Peacock, 1983).
Linear filters can be reliably estimated even in the presence of a output nonlinearity for white Gaussian inputs (Chichilnisky, 2001). We extracted linear filters from measured Gaussian odor stimuli and ORN responses using least-squares fitting. Given time series of input S and response R, we reshaped S into a matrix where each row of contained the stimulus up to N samples in the past where N is the length of the filter to be calculated. Using this matrix, we computed the stimulus covariance matrix . The linear kernel that is the best linear predictor of R given S is given by . However, since input signals have autocorrelation functions with power approaching 0 at high frequencies, estimated filters were occasionally dominated by high-frequency artifacts. To remove these, we regularized C, using where I is the identity matrix and r is a regularization factor in units of the mean eigenvalue of C. Finally, we normalized the filters as described in (Baccus and Meister, 2002) to preserve the units of the stimulus in the linear projections. The regularization value r was chosen by cross-validation to create filters that best predicted the data.
Linear filters extracted from naturalistic stimuli (Figure 1, Figure 1—figure supplement 4) are not unbiased estimates of the true filter, since the stimulus is not Gaussian. Similarly, if a front-end nonlinearity precedes a linear filtering step in a system, this fitting procedure will not yield an unbiased estimate of the true filter. Nonetheless, linear filters computed by least squares fitting are the filters that best predict R given S, in the least-squares sense.
We defined ORN input-output curves to be the output nonlinearity of a LN model, which were estimated by plotting ORN response vs. the projected stimulus, and then computing a piecewise linear function using 50 bins along the horizontal axis. Computing piecewise linear functions allowed us to visualize the output nonlinearity without making explicit assumptions of the functional form of the nonlinearity. Dashed lines in Figure 4e are the cumulative distributions of the stimulus, computed over all the data.
In general, for any system with a single stimulus and response, we define the gain of the system by measuring the slope of the nonlinearity in the best fit LN model, normalizing the filters to preserve the scale of the stimulus as in Baccus and Meister (2002). With Gaussian stimuli that only gently perturb the system, the nonlinearity is simply a straight line, and the gain is computed by the average slope of a linear fit to the output nonlinearity (Figure 3). When output nonlinearities are strong, we estimated gain by the slope at the midpoint of the nonlinearity (Figures 4 and 5). We measured three different gains: (1) transduction gains, from the stimulus to the LFP, (2) firing gains, from the LFP to the firing rate, and (3) overall ORN gains, from the stimulus to the firing rate. Transduction gain had units of (mV/V) since deflections in LFP are measured in mV and the stimulus is measured in V. Similarly, firing gain had units of Hz/mV. Overall ORN gain had units of Hz/V when stimulating ORNs with fluctuating odor, and had units of Hz/µW when stimulating with light.
In experiments where we changed the variance (Figure 4), low-variance epochs tended to have a mean stimulus ~8% higher than high-variance epochs (Figure 4—figure supplement 1), despite our best efforts to keep the stimulus mean the same. To estimate gain changes solely due to the change in stimulus variance, we divided the projected stimulus by the mean stimulus in each epoch in each trial (Figure 4e). Differences in gain between high- and low-variance epochs remain significant even without this correction (Figure 4—figure supplement 1b–c). A single filter was used to project stimuli in both low- and high-variance epochs; changes in gain from low- to high-variance stimuli thus appear solely in the nonlinearity (Figure 4—figure supplement 1d).
In all our data, we measured the stimulus together with the response of ORNs. This allowed us to estimate transduction and firing lags with respect to the stimulus. In general, we estimated response lags by computing cross-correlation functions from the stimulus to the response. Lag was defined to be the location of the peak of the cross-correlation function. (Figure 7a–b).
To estimate the timescale of variance gain control (Figure 4), we computed input-output curves from the projected stimulus to the firing rate in 50 ms bins, pooling all trials together. This allowed us to estimate gain in 50 ms bins, together with the stimulus contrast (standard deviation/mean). We plotted time series of instantaneous gain and stimulus contrast (Figure 4h), and observed that while the stimulus contrast changed rapidly after the switch from high to low variance (at t = 5 s), the instantaneous gain changed more slowly, but still changed in ~130 ms, suggesting that timescale of variance gain control is relatively rapid.
To determine if ab3A transduction-to-firing gain varied with the mean stimulus (Figure 5d), we used a Spearman rank correlation test. To determine if gains varied significantly from low to high variance epochs (Figures 4 and 5), we first reshaped the raw data into trials 10 s long. Each trial consisted of a high variance epoch followed by a low variance epoch. Each trial was fit with three linear models, and yielded three pairs of gains, for transduction gain, firing gain, and total ORN gain. We discarded all trials where any linear model fit was poorly correlated with data, retaining only trials where all fits had high correlations with data (r2 > 0.8, see Figure 4—figure supplement 1f–g). We used the Wilcoxon signed rank test on these tuples to determine if the difference in gains in the low and high variance epochs was statistically significant.
We assume that Or and Orco form a complex that can exist in two conformations that can bind ligand. The concentration of unbound active complexes is C* and that of unbound inactive one is C. The corresponding concentrations for the bound complexes are C*S and CS (Figure 8a). The fraction of active Or-Orco complexes is therefore (Figure 8a):
(Un)binding of odorant is taken to be much faster than the (in)activation. Thus, the probability to be bound in the active and inactive cases are
Here are the dissociation constants for each state. Let the free energy difference in units of between the active and inactive states be and when unbound and when bound, respectively. For simplicity, we assume detailed balance (this can easily be relaxed (Skoge et al., 2013)), which constrains the free energy difference between C*S and CS to be . The activation kinetics can then be described by
where the rates are
We constrained the energy barrier between the active and inactive conformations by making the simplifying assumption that , where is an intrinsic switching rate (see e.g. Skoge et al., 2013). From these considerations we can then derive Equation (1) from the main text, which describes the dynamics of the activity:
with the rates
Given a steady state signal S, the activity relaxes towards
Where the overbar indicates that Equation (1) is solved at steady state. We model adaptation by supposing that the activity feeds back onto the free energy difference with rates that depend only on the activity (the effective switching rate is constant):
where is the rate of adaptation. Note that the free energy is bounded both from below and from above. In practice, we only need the lower bound . At steady state (for values of S high enough that ) we have which implies that
Thus, at steady state, adaptation causes the free energy difference of the complex to increase with the logarithm of the background signal intensity.
Substituting into the definitions of the (in)activation rates, we get
When , the rates and are decreasing function of S in the range
where the approximation is valid when . In our case this bound on S is large (~40 V in PID measurement units) and in our experiments the rates are decreasing functions of S over the entire range.
The output of the model described above is a time series of the fraction of receptors that are active, . Receptor activation can lead to the opening of other channels, which results in a transduction current that we measure as changes in the LFP. To generate LFP responses from this, we use
where represents a convolution and is a linear time-invariant mono-lobed filter that is given by:
Since firing rates cannot be negative, and since the LFP to spiking transformation has been shown to be partly differentiating (Nagel and Wilson, 2011), we generated firing rate responses from the receptor activity using
where is an output nonlinearity which is a simple threshold linear function ). represents a convolution and is a linear time-invariant filter that is given by the sum of two other kernels:
where each kernel is parameterized by a Gamma function:
Changes in time scale and sensitivity in turtle photoreceptorsThe Journal of Physiology 242:729–758.https://doi.org/10.1113/jphysiol.1974.sp010732
The statistical computation underlying contrast gain controlJournal of Neuroscience 26:6346–6353.https://doi.org/10.1523/JNEUROSCI.0284-06.2006
Free-flight responses of Drosophila Melanogaster to attractive odorsJournal of Experimental Biology 209:3001–3017.https://doi.org/10.1242/jeb.02305
Non-linear blend coding in the moth antennal lobe emerges from random glomerular networksFrontiers in Neuroengineering 5:6.https://doi.org/10.3389/fneng.2012.00006
Navigational strategies used by insects to find distant, wind-borne sources of odorJournal of Chemical Ecology 34:854–866.https://doi.org/10.1007/s10886-008-9484-5
A simple white noise analysis of neuronal light responsesNetwork: Computation in Neural Systems 12:199–213.https://doi.org/10.1080/713663221
Theoretical NeuroscienceThe MIT Press.
A dynamical feedback model for adaptation in the olfactory transduction pathwayBiophysical Journal 102:2677–2686.https://doi.org/10.1016/j.bpj.2012.04.040
Flies require bilateral sensory input to track odor gradients in flightCurrent Biology 19:1301–1307.https://doi.org/10.1016/j.cub.2009.06.022
Intrinsic mechanisms for adaptive gain rescaling in barrel cortexJournal of Neuroscience 28:696–710.https://doi.org/10.1523/JNEUROSCI.4931-07.2008
The evolution of insect flight: implications for the evolution of the nervous systemBrain, Behavior and Evolution 50:8–12.https://doi.org/10.1159/000113317
Selectivity for multiple stimulus features in retinal ganglion cellsJournal of Neurophysiology 96:2724–2738.https://doi.org/10.1152/jn.00995.2005
Elemente Der PsychophysikBreitkopf und Härtel.
Single neuron computation: from dynamical system to feature detectorNeural Computation 19:3133–3172.https://doi.org/10.1162/neco.2007.19.12.3133
Intrinsic gain modulation and adaptive neural codingPLoS Computational Biology 4:e1000119.https://doi.org/10.1371/journal.pcbi.1000119
Foundations of Cellular NeurophysiologyThe MIT Press.
Measurement of Odor-Plume structure in a wind tunnel using a photoionization detector and a tracer gasEnvironmental Fluid Mechanics 2:115–142.https://doi.org/10.1023/A:1016227601019
Adaptation processes in Insect Olfactory receptorsAnnals of the New York Academy of Sciences 510:104–112.https://doi.org/10.1111/j.1749-6632.1987.tb43475.x
Chemo-electrical transduction in insect olfactory receptorsAnnual Review of Neuroscience 9:121–145.https://doi.org/10.1146/annurev.ne.09.030186.001005
Feed-Forward versus feedback inhibition in a basic olfactory circuitPLOS Computational Biology 11:e1004531.https://doi.org/10.1371/journal.pcbi.1004531
System identification of Drosophila olfactory sensory neuronsJournal of Computational Neuroscience 30:143–161.https://doi.org/10.1007/s10827-010-0265-0
Independent optical excitation of distinct neural populationsNature Methods 11:338–346.https://doi.org/10.1038/nmeth.2836
A simple coding procedure enhances a neuron's information capacityZeitschrift Fur Naturforschung. Section C, Biosciences 36:910–912.https://doi.org/10.1515/znc-1981-9-1040
Common strategies for light adaptation in the peripheral visual systems of fly and dragonflyJournal of Comparative Physiology ? A 128:319–340.https://doi.org/10.1007/BF00657606
Fractional differentiation by neocortical pyramidal neuronsNature Neuroscience 11:1335–1342.https://doi.org/10.1038/nn.2212
Intensity invariant dynamics and odor-specific latencies in olfactory receptor neuron responseJournal of Neuroscience 33:6285–6297.https://doi.org/10.1523/JNEUROSCI.0426-12.2013
Odor plumes and how insects use themAnnual Review of Entomology 37:505–532.https://doi.org/10.1146/annurev.en.37.010192.002445
Spatial and temporal structures of pheromone plumes in fields and forestsPhysiological Entomology 25:211–222.https://doi.org/10.1046/j.1365-3032.2000.00176.x
Biophysical mechanisms underlying olfactory receptor neuron dynamicsNature Neuroscience 14:208–216.https://doi.org/10.1038/nn.2725
Physiological features of the S- and M-cone photoreceptors of wild-type mice from single-cell recordingsThe Journal of General Physiology 127:359–374.https://doi.org/10.1085/jgp.200609490
Receptor expression and sympatric speciation: unique olfactory receptor neuron responses in F1 hybrid Rhagoletis populationsJournal of Experimental Biology 209:3729–3741.https://doi.org/10.1242/jeb.02444
Two-dimensional goodness-of-fit testing in astronomyMonthly Notices of the Royal Astronomical Society 202:615–627.https://doi.org/10.1093/mnras/202.3.615
Adaptation of the odour-induced response in frog olfactory receptor cellsThe Journal of Physiology 519 Pt 3:801–813.https://doi.org/10.1111/j.1469-7793.1999.0801n.x
Spikes: Exploring the Neural CodeMIT Press.
Differential intensity sensitivity of the ear for Pure TonesPhysical Review 31:867–875.https://doi.org/10.1103/PhysRev.31.867
Physical processes and real-time chemical measurement of the insect olfactory environmentJournal of Chemical Ecology 34:837–853.https://doi.org/10.1007/s10886-008-9490-7
The operating point of the cortex: neurons as large deviation detectorsJournal of Neuroscience 27:7673–7683.https://doi.org/10.1523/JNEUROSCI.1048-07.2007
Classification of odorants across layers in locust olfactory pathwayJournal of Neurophysiology 115:2303–2316.https://doi.org/10.1152/jn.00921.2015
The Handbook of Brain Theory and Neural NetworksAmplification, attenuation, and integration, The Handbook of Brain Theory and Neural Networks.
The effect of contrast on the transfer properties of cat retinal ganglion cellsThe Journal of Physiology 285:275–298.https://doi.org/10.1113/jphysiol.1978.sp012571
Chemical sensing by Nonequilibrium Cooperative receptorsPhysical Review Letters 110:248102–248105.https://doi.org/10.1103/PhysRevLett.110.248102
Multiple perceptible signals from a single olfactory glomerulusNature Neuroscience 16:1687–1691.https://doi.org/10.1038/nn.3519
Determinants of the Drosophila odorant receptor patternDevelopmental Cell 22:363–376.https://doi.org/10.1016/j.devcel.2011.12.015
De Pulsu, Resorptione, Auditu Et Tactu. Annotationes Anatomicae Et Physiologicae, Auctoreprostat apud C.F. Koehler.
Fred RiekeReviewing Editor; Howard Hughes Medical Institute, University of Washington, United States
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
[Editors’ note: a previous version of this study was rejected after peer review, but the authors submitted for reconsideration. The first decision letter after peer review is shown below.]
Thank you for submitting your work entitled "Olfactory receptor neurons use fast dynamic gain control to encode intermittent odorant stimuli" for consideration by eLife. Your article has been favorably evaluated by a Senior Editor and three reviewers, one of whom, Fred Rieke (Reviewer #1), is a member of our Board of Reviewing Editors. The following individual involved in review of your submission has agreed to reveal their identity: Tim Holy (Reviewer #2).
Our decision has been reached after consultation between the reviewers. Based on these discussions and the individual reviews below, we regret to inform you that your work will not be considered further for publication in eLife, at least without substantial revisions.
Two main points emerged in discussion among the three reviewers about the paper. We highlight these points because we felt they could change how the data is interpreted. All reviewers agreed that these points were important to deal with before we could consider the paper further. Individual reviews, attached below, will both elaborate on these points as well as detail other important issues to consider.
Adaptation vs. saturation.
Adaptation as defined in the paper includes both time-dependent changes in cellular processing as well as static nonlinearities in the response – such as receptor saturation. We felt it was quite important to separate these effects as much as possible. This could, for example, involve fitting NLN models to the responses, providing a "front-end" nonlinearity that could account for receptor saturation. Several related comments in the individual reviews involve how gain is defined and how much the conclusions depend on the definition of gain.
Validity of LN models.
The LN models the form the core analysis tool in the paper are not validated. Several assumptions are central to the analysis. Some of these points may also originate from a lack of clarity in the paper about details of how the LN models were fit. Most importantly, it is not clear in all cases how the linear filter used to project the stimulus was computed. For example, since the standard linear filter calculation assumes symmetry of the stimulus, it is not clear how filters were computed for the naturalistic odor stimuli.
This paper describes adaptation to the mean and variance of odorant stimuli, with a particular focus on adaptation to natural inputs. The question is quite interesting, and for the most part I found the conclusions of the paper well supported by the data presented. Several issues, however, limit my enthusiasm.
How good are LN model fits to the data?
The paper relies extensively on LN models to analyze adaptation. More validation of this approach is needed – e.g. how much of the response variance do LN models describe? Related to this point – are the linear filters used in analysis of responses to naturalistic inputs derived from responses to Gaussian noise inputs (this seems to be the case from the Methods)? If not, how are the filters computed? If the filters are computed for Gaussian noise inputs, are they still valid for naturalistic inputs? Similarly, in Figure 4 – are the low variance and high variance filters assumed to be the same, and is this indeed correct?
LFP as a proxy for transduction.
The LFP is used to estimate properties of odor transduction. It is not clear from the paper how good a proxy the LFP is for transduction, specifically whether it is accurate enough to support the conclusions drawn. More verification of the validity of using the LFP is needed.
The gain in Figure 2 is defined as the slope of the firing rate vs. projected stimulus curve. But doesn't that mean comparing different firing rates for different whiffs – and hence potentially confounding adaptive effects with static nonlinearities in the responses? This issue also comes up in Figure 3B.
Compensation for slowing of transduction.
The conclusion that speeding of post-transduction processes with adaptation compensate for slowing of transduction is not well developed, especially for a major conclusion of the paper. This relates to the point above about using the LFP as a proxy for transduction.
The submission by Gorur-Shandilya and colleagues examines the dynamics of odorants responses in Drosophila ORNs using both naturalistic (plume-based) stimuli and gaussian stimuli of controlled mean and variance. The authors show that adaptation occurs on rapid time scales in a manner consistent with the well-known Weber-Fechner "law." By measuring both field potentials and spiking, the authors attempt to dissociate transduction dynamics from spiking dynamics, and discover that the two have compensatory kinetic characteristics, so that overall the firing preserves the dynamics of the stimulus.
This study brings the analysis of Drosophila ORN dynamics closer to the standard set in, e.g., phototransduction and several other sensory systems, and for that reason has considerable value. Overall, the results suggest a fair degree of correspondence with the early stages of visual encoding, right down to the magnitude of the relevant time scales. As the authors acknowledge, it is primarily a "phenomenology" study, leaving most questions of mechanism for future work. (The approximate dissociation of spiking vs. transduction is something of an exception, however.) However, the history of other systems argues that quantitative characterization of the phenomenology is an essential step along the path to discovering the underlying machinery.
At a big-picture level, the Introduction sets the stage well, the manuscript as a whole appears scholarly, and the claims seem well-supported by the data and analyses. A key component of this work is the impressive control over stimulus delivery, which is much harder for olfactory stimuli than for visual stimuli, and the authors clearly went to great lengths to achieve this control. Consequently, this study has a lot to recommend it, and little to dislike. My concerns are minor, focusing on the clarity of descriptions, particularly in the Materials and methods section (which is not up to the standard set in the rest of the manuscript). In its current state I suspect that readers will have considerable difficulty understanding exactly what was done at several points in the manuscript. There are some places where I question details of the results, but it seems likely that this is due to a misunderstanding of the precise form of the analysis.
Figure 2D: one oddity is that the tangent of the f-s curve even at 0 appears strongly correlated with the ultimate size of the excursion. Naively, this would appear to violate causality, but of course it could reflect correlations in the stimulus. Nevertheless, it begs the question of how much of this analysis is really dominated by the stimulus correlations; there don't appear to be any examples of a large stimulus that were not preceded by some adaptation, and that leads to concern about the thoroughness of this particular test and analysis.
This paper examines the responses of Drosophila olfactory receptor neurons to dynamic fluctuating stimuli and analyzes their responses predominantly using linear filter methods. Based on these analyses the authors argue that ORNs exhibit fast adaptation that allows them to better encode the properties of odor plumes. The experiments appear to be carefully performed and the data appear to be solid. However, I have concerns about the novelty of the approach and findings, as well as some of the analyses and interpretations.
Linear filter methods have been used for many years to approximate the response of a sensory neuron to a complex time-varying stimuli. A common theme of these studies is the finding that the linear approximation of the response varies as a function of stimulus statistics such as the mean and variance (e.g. Fairhall 2001, Baccus and Meister 2002). In particular, the approximate gain of the neuron generally decreases with increased stimulus variance, a phenomenon known as "gain scaling" or "variance scaling". Following on these studies several groups asked whether single neurons could exhibit this property in response to current injection, and found that it could. Subsequent modeling studies have shown that Hodgkin-Huxley models and other spiking models are capable of producing gain scaling (e.g. Lundstrom et al., 2008). Therefore, the finding that Drosophila ORNs exhibit gain scaling, and that some of this scaling arises at the level of spike generation does not seem novel or surprising.
In the Drosophila olfactory neurons studied here, several papers have examined responses to temporally modulated stimuli (e.g. Kim 2010, Nagel 2011, Martelli 2013, Cafaro, 2016, Cao, 2016). Multiple studies have found that transduction response magnitude (measured intracellularly or by LFP) is well described by a Hill equation, (that is, it rises sigmoidally with log odor concentration), as one would predict from the basic biochemistry of receptor activation (Cao, Cafaro, Nagel). Second, they have found that OR-expressing ORNs exhibit calcium-dependent adaptation that dynamically adjusts sensitivity and gain (Cao, 2016). The approach taken here confounds these two phenomena, as it assumes that responses are a linear function of concentration (here shown to be linearly related to PID voltage). Therefore, the "adaptation" reported by the authors likely contains contributions from the nonlinearity inherent in receptor activation, as well as true dynamical adaptation. This might be addressed by computing filters based on the log concentration, rather than linear concentration, and by comparing "adaptation" (estimated from projections onto linear filters) in OR- and IR expressing ORNs, as IR-expressing ORNs do not exhibit calcium-dependent adaptation, but do show activation nonlinearities.
An issue in the paper is what is meant by "adaptation." The authors seem to argue that adaptation is any change in the linear approximation to the response. For example, in the Introduction they state that they are focusing on phenomenology, which they call "a prerequisite to understanding the mechanisms that implement them." In the Discussion, they say that "'adaptation state is a construction to understand nonlinearities in that response." However, many other parts of the paper imply that these "adaptations" are actively produced by the cell in order to optimize coding. For example, in the title "Olfactory receptor neurons use fast dynamic gain control", Results "small isolated whiffs were amplified…suggesting that ORNs adjusted their gain dynamically", "suggesting that ORNs actively changed their gain." The question of whether certain coding phenomena arise from active adaptation processes, or from unavoidable nonlinearities of receptor-ligand interactions seems germane to the framing of the paper, and the question of whether the phenomena they describe are specializations of insect neurons for encoding odor plumes or general properties of receptor-activated neurons.
A second issue is the characterization of phenomenology versus mechanism. Although the authors claim to focus on phenomenology, a significant part of the paper is devoted to the distinguishing phenomena arising at the level of transduction/LFP, and those arising at the level of spikes- an attempt to place phenomenology in a mechanistic context. Given that multiple stages have been experimentally identified in the transduction process (Cao et al. 2016), and that many models of spike generation are available, it does not seem fair to call this phenomenological characterization a "necessary prerequisite" to mechanistic understanding.
The finding that ORN spikes seem to counteract some of the slowing of transduction responses is interesting and novel but not well-developed.
I am concerned about the linear filter method applied to responses to sparse stimuli in Figure 2. As stated in the Materials and methods, linear filters can be reliably estimated in the presence of an output nonlinearity if the stimuli are Gaussian- however, this is generally not true for stimuli that are strongly skewed, such as the one used here.
In general the authors do not report the goodness of fit of their linear model, making it difficult to interpret whether the observed changes in gain (where gain is the slope of the projected versus real response) are significant.
The authors state: "It remains unclear whether the ORN's ultimate output- the firing rate- follows the same Weber-Fechner scaling." However, Cafaro, 2016 does address this question for ORN spikes, as well as PN membrane potential and spikes.
Introduction, last paragraph: It is unclear at this point in the text why the authors choose to focus on a history of 300ms. This should be better motivated.
Subsection “Fast gain control maintains timing information of naturalistic odor signals”, third paragraph. It would be helpful to at least briefly mention how lags were computed in the Results.
Subsection “Variance gain control in olfaction”. No potential mechanisms for gain scaling of transduction responses are given.
Subsection “Fast gain control could aid in naturalistic odor detection”. The authors state that paradigms that employ conditioning and probe stimuli cannot easily quantify the dynamics of gain change. This does not seem accurate. These paradigms explicitly allow the experimenter to distinguish the stimuli used to generate adaptation from those used to test it. Small combinations of pulses are used to examine the dynamics of gain change in transduction in Cao et al. 2016.
Subsection “Stimulus measurement”. It would be helpful to include a diagram of the odor delivery apparatus along with calibration data.
[Editors’ note: what now follows is the decision letter after the authors submitted for further consideration.]
Thank you for submitting your article "Olfactory receptor neurons use gain control and complementary kinetics to encode intermittent odorant stimuli" for consideration by eLife. Your article has been reviewed by two peer reviewers, and the evaluation has been overseen by a Reviewing Editor and a Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Tim Holy (Reviewer #3).
The reviewers have discussed the reviews with one another all agreed that the revisions have improved the paper substantially. There are a few outstanding issues outlined in the reviews below.
The revised manuscript describes encoding of dynamic stimuli by ORNs. The authors have addressed several of the major concerns cited by all reviewers. First, they find that a static front-end nonlinearity (a Hill function) can account for much of the nonlinearity observed in response to naturalistic odor stimuli, but that a dynamic adaptation component remains. Second, they have expanded their analysis of the dynamics of spiking and transduction, to show that speeding of spiking compensates for slowing of transduction with increased stimulus mean, resulting in spike timing that is independent of stimulus mean- an interesting and novel result. Finally, they have developed a receptor-based model to account for the observed forms of adaptation.
Overall the manuscript is much improved. I think it makes a substantial contribution to the literature. I do have some remaining issues with the manuscript that I think should be resolved prior to publication.
1) The Hill function as a source of nonlinearity.
The reviewers previously asked whether many of the nonlinearities described in the data could arise from receptor-ligand interactions, which should give rise to a nonlinear relationship between odor concentration and transduction. The authors have re-analyzed their data and found that responses to natural stimuli are well-described by such a relationship (Figure 1E-F), although context-dependent deviations from this relationship remain (Figure 2). However, it seems that this understanding is not entirely integrated into the text. For example, the title of the first section is "ORN responses to naturalistic odorant stimuli show deviations from linearity that arise from adaptation and saturation" which implies that saturation is the only nonlinearity (not the logarithm implied by the Hill function). Much time is spent on showing that responses deviate from linearity, and then receptor-ligand interactions are suggested as an alternative to output nonlinearities, which seems odd, given that receptor-ligand nonlinearities are an obligate feature of a system composed of odorant receptors that bind and are activated by odorants. A finding that odor concentration is encoded linearly by such a system would be much more surprising. Along these lines, in the Introduction, the authors state that "ORNs employed front-end nonlinearities" which makes it sound like an active choice on the part of ORNs. I think a fairer statement of the findings would be that a front-end nonlinearity, combined with dynamic adaptation that shifts the midpoint of this nonlinearity, can account for ORN responses to a variety of dynamic stimuli.
2) Explanation of Weber-Fechner scaling.
Understanding the origin of Weber-Fechner scaling is a major goal of the paper. Along the lines above I would suggest that Figure 3—figure supplement 2 should be made part of the main text. This shows that the expected front-end nonlinearity, combined with adaptation that shifts the activation curve to the right, is sufficient to account for ORN responses to stimuli with different means. As such a rightward shift with adaptation has been previously reported (Kaissling, 1987, Nagel and Wilson, 2011), this makes for an elegant explanation of the observed responses.
The authors state that the receptor model alone cannot account for Weber-Fechner scaling (response to reviewers 1B). It is true that the receptor model alone does not produce responses that change gain and have similar mean, as in the data. However, it does seem that the receptor model alone (Figure 3—figure supplement 2, panel c) shows a decrease in gain with increased mean. If the output of this simple model were plotted as in Figure 3F what would it look like? Does the change in gain arise from the dynamic adaptation (the shift in half-max) or from the nonlinearity?
3) Receptor-based model.
The authors present a new model, based on a 2-state receptor model, to account for the observed adaptation. As I understand it, their model makes two assumptions: (1) receptor activation/inactivation are the rate limiting steps, and (2) receptor activity feeds back onto these rates, slowing both. The model thus makes specific predictions about what biophysical steps give rise to adaptation.
Given the findings in Figure 3—figure supplement 2, it seems like the main features necessary to explain the results are (1) steady-state output rises as S/S+K, (2) adaptation shifts K to the right, (3) adaptation decreases transition rates, leading to slower responses. The authors should clarify which aspects of their model are essential for reproducing the key features of the data.
4) Use of the LFP as a proxy for transduction.
It is true that it has become somewhat standard in the literature to do this but the LFP should be interpreted with a bit of caution. In Nagel and Wilson 2011, palp LFPs were used because LFPs in this structure were found to more closely reflect activity in a single sensillum (presumably because of the greater spacing between sensilla). A brief caveat noting that LFPs can contain a contribution from nearby sensilla would be welcome.
The revisions thoroughly address my concerns. I am quite satisfied with the revised manuscript.https://doi.org/10.7554/eLife.27670.024
- Srinivas Gorur-Shandilya
- Mahmut Demir
- Thierry Emonet
- Damon A Clark
- Damon A Clark
- Damon A Clark
- Junjiajia Long
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We gratefully acknowledge helpful conversations with John Carlson, Julijana Gjorgjieva, Charles Greer, Joe Howard, Carlotta Martelli, and Steven Zucker. We are thankful to John Carlson for giving us access to some of the instrumentation in his lab. MD, SG-S and TE were partially supported by the Whitehall Foundation. JL was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) Postgraduate Scholarships-Doctoral Program PGSD2-471587-2015. DAC was supported by a Sloan Research Fellowship, a Searle Scholar Award, and the Smith Family Foundation.
- Fred Rieke, Reviewing Editor, Howard Hughes Medical Institute, University of Washington, United States
© 2017, Gorur-Shandilya et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.