Olfactory receptor neurons use gain control and complementary kinetics to encode intermittent odorant stimuli
Abstract
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.001eLife digest
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.002Introduction
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.
Results
ORN responses to naturalistic odorant stimuli show deviations from linearity that arise from adaptation and front end nonlinearities
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).
ORNs adapt to stimulus background by decreasing gain according to the Weber-Fechner Law
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.
Fast variance-dependent gain control in ORNs
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).
Mean and variance gain control occur at different stages of odor encoding, and are mechanistically distinct
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).
Modularity of gain control at transduction and spiking
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.
Despite slowdown in transduction, ORN firing rate preserves timing of odor encounters
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.
An adaptive two-state receptor-complex model reproduces Weber-Fechner scaling, slow down of LFP kinetics, and responses to intermittent and Gaussian stimuli
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).
Discussion
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 in olfaction
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.
Variance gain control in olfaction
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.
Models and mechanisms of ORN response and gain control
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.
Dynamic gain control could aid in naturalistic odor detection
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.
Invariant firing rate kinetics could improve odor-guided flight behavior
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.
Materials and methods
Electrophysiology
Single sensillum recordings
Request a detailed protocolSingle 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.
Spike sorting
Request a detailed protocolAll 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.
Local field potentials
Request a detailed protocolThe 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.
Fly stocks and genetic strategies
Request a detailed protocolFlies 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.
Stimulus measurement
Request a detailed protocolWe 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.
Odor stimulus generation
General principle
Request a detailed protocolOdorants 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.
Controlling air flows
Request a detailed protocolMass 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.
Naturalistic stimulus (Figure 1)
Request a detailed protocolTo 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.
Stimulus with changing mean (Figures 2 and 4)
Request a detailed protocolTo 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.
Stimulus with changing variance (Figures 3–4)
Request a detailed protocolWe 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.
Numerical methods
Statistics of naturalistic odor
Request a detailed protocolSince 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).
Estimating deviations in response to naturalistic stimulus
Request a detailed protocolFor 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).
Filter extraction
Request a detailed protocolLinear 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.
ORN input output curves (solid lines in Figures 3e and 4e etc.)
Request a detailed protocolWe 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.
Estimation of gain
Request a detailed protocolIn 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).
Measuring lags
Request a detailed protocolIn 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).
Estimation of variance gain control timescale
Request a detailed protocolTo 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.
Statistical tests
Request a detailed protocolTo 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.
Modelling
Stimulus binding and the activity of the Or-Orco complexes
Request a detailed protocolWe 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.
Kinetic slowdown upon adaptation
Request a detailed protocolSubstituting 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.
Receptor activity to LFP
Request a detailed protocolThe 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:
Receptor activation to firing rates
Request a detailed protocolSince 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:
Code availability
Request a detailed protocolspikesort. MATLAB toolbox to sort spikes from extracellular recordings. Available at https://github.com/emonetlab/spikesort.
kontroller. MATLAB toolbox to acquire data and run experiments on National Instruments hardware. Available at https://github.com/emonetlab/kontroller
References
-
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
-
Light adaptation and photopigment bleaching in cone photoreceptors in situ in the retina of the turtleJournal of Neuroscience 14:1091–1105.
-
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
-
Odor Landscapes in turbulent environmentsPhysical Review X 4:041015–041017.https://doi.org/10.1103/PhysRevX.4.041015
-
A simple white noise analysis of neuronal light responsesNetwork: Computation in Neural Systems 12:199–213.https://doi.org/10.1080/713663221
-
Dynamical adaptation in photoreceptorsPLoS Computational Biology 9:e1003289.https://doi.org/10.1371/journal.pcbi.1003289
-
Triggered correlationIEEE Transactions on Biomedical Engineering 15:169–179.https://doi.org/10.1109/TBME.1968.4502561
-
Odor coding in a model olfactory organ: the Drosophila maxillary palpJournal of Neuroscience 19:4520–4532.
-
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
-
Contrast adaptation in a nonadapting LGN modelJournal of Neurophysiology 98:1287–1296.https://doi.org/10.1152/jn.00618.2006
-
The Weber ratio for intensive discriminationPsychological Review 43:322–340.https://doi.org/10.1037/h0059748
-
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
-
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
-
Temporal contrast adaptation in the input and output signals of salamander retinal ganglion cellsJournal of Neuroscience 21:287–299.
-
Independent optical excitation of distinct neural populationsNature Methods 11:338–346.https://doi.org/10.1038/nmeth.2836
-
Olfactory coding with all-or-nothing glomeruliJournal of Neurophysiology 98:3134–3142.https://doi.org/10.1152/jn.00560.2007
-
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
-
Temporal contrast adaptation in salamander bipolar cellsJournal of Neuroscience 21:9445-54.
-
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
-
Precise olfactory responses tile the sniff cycleNature Neuroscience 14:1039–1044.https://doi.org/10.1038/nn.2877
-
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
-
Olfactory adaptation depends on the trp Ca2+ channel in DrosophilaJournal of Neuroscience 19:4839–4846.
-
The influence of visual landscape on the free flight behavior of the fruit fly Drosophila MelanogasterThe Journal of Experimental Biology 205:327–343.
-
BookDe Pulsu, Resorptione, Auditu Et Tactu. Annotationes Anatomicae Et Physiologicae, Auctoreprostat apud C.F. Koehler.
-
Robust perfect adaptation in bacterial chemotaxis through integral feedback controlPNAS 97:4649–4653.
Article and author information
Author details
Funding
Whitehall Foundation
- Srinivas Gorur-Shandilya
- Mahmut Demir
- Thierry Emonet
Sloan Research Fellowship
- Damon A Clark
Searle Scholar Award
- Damon A Clark
Smith Family Foundation
- Damon A Clark
Natural Sciences and Engineering Research Council of Canada (PGSD2-471587-2015)
- Junjiajia Long
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
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.
Copyright
© 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.
Metrics
-
- 3,245
- views
-
- 567
- downloads
-
- 89
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Neuroscience
Hypothalamic kisspeptin (Kiss1) neurons are vital for pubertal development and reproduction. Arcuate nucleus Kiss1 (Kiss1ARH) neurons are responsible for the pulsatile release of gonadotropin-releasing hormone (GnRH). In females, the behavior of Kiss1ARH neurons, expressing Kiss1, neurokinin B (NKB), and dynorphin (Dyn), varies throughout the ovarian cycle. Studies indicate that 17β-estradiol (E2) reduces peptide expression but increases Slc17a6 (Vglut2) mRNA and glutamate neurotransmission in these neurons, suggesting a shift from peptidergic to glutamatergic signaling. To investigate this shift, we combined transcriptomics, electrophysiology, and mathematical modeling. Our results demonstrate that E2 treatment upregulates the mRNA expression of voltage-activated calcium channels, elevating the whole-cell calcium current that contributes to high-frequency burst firing. Additionally, E2 treatment decreased the mRNA levels of canonical transient receptor potential (TPRC) 5 and G protein-coupled K+ (GIRK) channels. When Trpc5 channels in Kiss1ARH neurons were deleted using CRISPR/SaCas9, the slow excitatory postsynaptic potential was eliminated. Our data enabled us to formulate a biophysically realistic mathematical model of Kiss1ARH neurons, suggesting that E2 modifies ionic conductances in these neurons, enabling the transition from high-frequency synchronous firing through NKB-driven activation of TRPC5 channels to a short bursting mode facilitating glutamate release. In a low E2 milieu, synchronous firing of Kiss1ARH neurons drives pulsatile release of GnRH, while the transition to burst firing with high, preovulatory levels of E2 would facilitate the GnRH surge through its glutamatergic synaptic connection to preoptic Kiss1 neurons.
-
- Computational and Systems Biology
Degree distributions in protein-protein interaction (PPI) networks are believed to follow a power law (PL). However, technical and study bias affect the experimental procedures for detecting PPIs. For instance, cancer-associated proteins have received disproportional attention. Moreover, bait proteins in large-scale experiments tend to have many false-positive interaction partners. Studying the degree distributions of thousands of PPI networks of controlled provenance, we address the question if PL distributions in observed PPI networks could be explained by these biases alone. Our findings are supported by mathematical models and extensive simulations and indicate that study bias and technical bias suffice to produce the observed PL distribution. It is, hence, problematic to derive hypotheses about the topology of the true biological interactome from the PL distributions in observed PPI networks. Our study casts doubt on the use of the PL property of biological networks as a modeling assumption or quality criterion in network biology.