Abstract
In Drosophila, olfactory information received by olfactory receptor neurons (ORNs) is first processed by an incoherent feed forward neural circuit in the antennal lobe (AL) that consists of ORNs (input), inhibitory local neurons (LNs), and projection neurons (PNs). This “early” olfactory information processing has two important characteristics. First, response of a PN to its cognate ORN is normalized by the overall activity of other ORNs, a phenomenon termed “divisive normalization.” Second, PNs respond strongly to the onset of ORN activities, but they adapt to prolonged or continuously varying inputs. Despite the importance of these characteristics for learning and memory, their underlying mechanisms are not fully understood. Here, we develop a circuit model for describing the ORN-LN-PN dynamics by including key neuron-neuron interactions such as short-term plasticity (STP) and presynaptic inhibition (PI). By fitting our model to experimental data quantitatively, we show that a strong STP balanced between short-term facilitation (STF) and short-term depression (STD) is responsible for the observed nonlinear divisive normalization in Drosophila. Our circuit model suggests that either STP or PI alone can lead to adaptive response. However, by comparing our model results with experimental data, we find that both STP and PI work together to achieve a strong and robust adaptive response. Our model not only helps reveal the mechanisms underlying two main characteristics of the early olfactory process, it can also be used to predict PN responses to arbitrary time-dependent signals and to infer microscopic properties of the circuit (such as the strengths of STF and STD) from the measured input-output relation. Our circuit model may be useful for understanding the role of STP in other sensory systems.
1. Introduction
Sensory systems have evolved different strategies to efficiently represent and process physiologically relevant stimuli in the presence of various biophysical constraints. For example, the olfactory system is confronted with the challenge that there are numerous odors each consisting of multiple volatile molecules with a wide range of concentrations. Yet the olfactory system possesses a remarkable ability to detect and discriminate odors using a relatively small repertoire of odor receptors (ORs) through a combinatorial code, i.e., each odorant is sensed by multiple receptors and each receptor can be activated by many odorants (Hallem and Carlson, ; Saito et al., ; Si et al., ).
The functional organization of the olfactory systems across different species is highly conserved (Su et al., ; Hansson and Stensmyr, ; Uchida et al., ). In both insects and vertebrates, an olfactory receptor neuron (ORN) typically expresses only one type of OR. ORNs that express the same OR converge to the same glomerulus in the olfactory bulb (vertebrates) or the antennal lobe (AL, insects). In Drosophila, peripheral odor information is processed in the AL before transmitted to higher brain areas by projection neurons (PNs) (Masse et al., ; Wilson, ). Each PN typically innervates one glomerulus. The transfer function between ORN and PN is a saturating nonlinear function (Bhandawat et al., ; Kazama and Wilson, ; Olsen and Wilson, ), i.e., a small ORN input is disproportionally amplified while a strong input saturates the response. Lateral inhibition by local interneurons (LNs) in the AL increases the level of ORN input needed to drive PNs to saturation, the strength of inhibition scales with the total forward input to the AL, a phenomenon called “divisive normalization” (Olsen et al., ; Carandini and Heeger, ), which has been widely observed across different sensory modalities and brain regions (Carandini and Heeger, ; Ferguson and Cardin, ). Divisive normalization in the AL was found beneficial for efficient odor coding (Olsen and Wilson, ; Luo et al., ).
Airborne odors are intermittent and have complex spatio-temporal profiles (Murlis et al., ; Vickers et al., ). The ability to detect and respond to temporal variation of odors is crucial for successful odor-guided navigation (David et al., ; Victor et al., ; Demir et al., ). This is partially achieved by the adaptive responses of PNs to time-dependent inputs from ORNs. In particular, PNs respond transiently to the onset of a step-function like ORN input and fall back to low firing rates for the prolonged input, showing highly adaptive response. For more complex time-dependent ORN inputs, the response of PNs depends on both the ORN firing rate and its rate of change (Kim et al., , ). Indeed, one of the hallmark properties of all sensory systems is adaptation, which is crucial for detecting and tracking time-varying signals (Wark et al., ). For adaptive response in Drosophila olfactory system, although a phenomenological linear-nonlinear model was proposed to fit experimental data (Kim et al., ), a mechanistic understanding of how the ORN-PN-LN circuit in the AL leads to the adaptive response is still missing.
The aim of this study is to understand the mechanistic origins of both divisive normalization and adaptive response in Drosophila by modeling dynamics of the AL neural circuit. From previous studies (Kazama and Wilson, ; Martelli and Fiala, ), the synapses between ORNs and PNs in AL exhibit strong short-term plasticity (STP), which is a form of fast activity-dependent modulation of synaptic strength (Stevens and Wang, ; Abbott et al., ; Markram et al., ; Dittman et al., ; Wang et al., ). Another important factor in the AL circuit is that inhibition by LNs is due to presynaptic inhibition (PI) at the axon terminal of ORNs (Olsen and Wilson, ). In this paper, we develop a simple circuit model of the Drosophila AL that includes both STP and PI. By using analytical methods and numerical simulations, we study the effects of STP on divisive normalization and response to time-varying stimuli of PNs. From direct quantitative comparison of our model results to experimental data, we show that STP is essential for the observed highly nonlinear divisive normalization; and both STP and PI determine the adaptive response observed in experiments.
2. Methods
2.1. A Circuit Model of the Antennal Lobe With STP
There are around 50 types of ORNs in Drosophila melanogaster, each of them expresses one OR. ORNs that express the same OR converge to the same glomerulus in the AL. A given odor typically activates several types of ORNs, hence different glomeruli. Each PN innervates one glomeulus and projects to higher brain areas like mushroom body and lateral horn. Since most LNs in AL are GABAergic, we will only consider inhibitory interneurons. Although LNs have distinct morphologies, innervation patterns, and response dynamics to odors (Chou et al., ; Nagel et al., ; Nagel and Wilson, ), for the purpose of this study, we do not differentiate them in our model. Generally speaking, LNs innervate different glomeruli and target the boutons of ORN axons, forming presynatpic inhibition (Figure 1A). We consider the simplified neural circuit of ORN-PN-LN in the AL, as shown in Figure 1B. Since ORNs promote the firing of PNs and LNs, while LNs in turn inhibit the firing of PNs, these neurons forms an incoherent feedforward loop (IFFL), a motif that has been widely observed in biochemical networks (Shen-Orr et al., ; Ma et al., ; Tu and Rappel, ).
Figure 1
For simplicity, we used a mean-field model to describe the firing rates of the neurons (Dayan and Abbott, ; Trappenberg, ; Gerstner et al., ). The dynamics of PN (LN) population firing rate RPN(LN) can be written as: where R is the firing rate of the cognate ORN that responds to a particular (private) odorant. The sum in Equation (2) is over all ORNs connected to the LN, including non-cognate ORNs that respond only to public odorant(s). The timescale τE is the relaxation time of the firing rate. ωEE and ωIE are synaptic weights of the synapses from ORN to PN and LN, respectively, assumed to be homogeneous among different ORNs.
The effect of PI is modeled by Equation (1) with a (dimensionless) variable 0 < p < 1 that represents reduction of the effective ORN firing rate due to presynaptic inhibition by LN. The dynamics of p, simplified from previous studies (Zhang et al., , ), is modeled as: where ρ is a constant and τp is the relaxation time of p. In the limit τp ≪ τE, we can use the quasi-steady state approximation , which suggests that p decreases with RLN.
The effect of STP can be separated into short-term facilitation (STF) and short-term depression (STD), which are modeled by u+ and x in Equation (1), respectively. Following previous work (Tsodyks et al., ), we denote u− (u+) as the neurotransmitter releasing probability just before (after) the arrival of a presynaptic spike; and x as the fraction of available neurotransmitters (Figure 1A). Applying the mean-field model for STP (Tsodyks et al., ), we have the following dynamics for u− and x: where u+ = u− + U(1 − u−) with U as the increment in release probability after each spike. Without any presynaptic firing (R = 0), we have x = 1, u− = 0, and u+ = U at steady state. With presynaptic firing, x decreases and u+ increases before their steady state values are recovered with time constants τD and τF respectively. The strength of STF and STD can be measured by the dimensionless recovery times and normalized by the relaxation time for firing rate τE. The longer the STP recovery time (τF or τD), the stronger the STP effect (STF or STD).
In the rest of the paper, we use the neural circuit model (Equations 1–5) to describe and explain several response properties of PNs including divisive normalization for constant (steady state) inputs and adaptive response to time-varying inputs. In both cases, we compare our model results with existing experiments and focus on understanding the effects of STP and PI on the observed behaviors.
2.2. Response of PNs to Triangular-Shaped Firing Rates of ORNs
In this section, we describe in detail the approximations used when we study PN's response to triangular-shaped inputs. First, we consider a simpler scenario, where ORN's firing rate increases linearly without bound, i.e., R = Kt. Then, Equations (1)–(5) become To analyze the adaptive behavior of PN responses, we derive an approximate solution of the above equations. These equations can be separated into two groups based on the biological mechanisms they describe: the first two equations describe presynaptic inhibition and the other equations are related to STP, which we study separately in the following.
2.2.1. The Presynaptic Inhibition Equations
Solving the first two equations of Equations (6) we obtain: with initial conditions p(0) = 1 and RLN(0) = 0. Here , as defined previously in the steady state solution in the main text. Since the plateau appears at t ≫ τE, the integrand of Equation (7) can be approximated as The validity of this approximation is supported by direct numerical integration of Equation (7) (Supplementary Figure 6). Plugging it back into Equation (7), we obtain: where denotes the exponential integral. By asymptotically expanding Ei(x) as , we have (up to the leading order terms): Multiplying Equation (9) by Kt gives the effective input Reff Equation (14).
2.2.2. The STP Equations
The last four equations in (6) describe the STP mechanism. From Equation (9), in the limit t ≫ τp, τE, we have p ≈ (AR(t))−1, which leaves the effective input Reff = p(t)R(t) approximately constant in time. We use this approximation to investigate the plateau behavior of PN response. Plugging it back into Equation (6), we have where 1/τu ≡ 1/τF + ReffU.
The equation for x now becomes: where . Introducing , we have Since we are interested in the plateau behavior of PNs, where t ≫ τu, τx. We can approximate as . Plugging it back to Equation (10) and neglecting higher order corrections, we have x(t) ≈ τx/τD. Similarly, The magnitude of the plateau response of PN can be approximated by its steady state activity for large t, where Reff = p(t)Kt ≈ 1/A with p(t) given by Equation (9). We see that the plateau magnitude of RPN is related to R(t) only through Reff, which is the final value that pR(t) reaches. It is independent of t and K, therefore the magnitude of PN response plateau is not affected by the input changing rate K, but only depends on intrinsic properties of the system.
2.3. Determination of Model Parameters and Initial Conditions
All model parameters were taken from previous studies or estimated by fitting to experimental data with their values constrained in reasonable physiological ranges. In Figure 2, each line is associated with a public odor at certain concentration. Since the firing rate due to all ORNs are not available in the experiment, and only local field potential was measured. We fitted each line by a independent . The resulting are found proportional to local field potentials measured in experiments (Supplementary Figure 1). The best-fit parameters (Table 1) were obtained by minimizing the residual sum of square for all the 16 (DL5) or 15 (VM7) points for each ORN-PN pair. The qualitative behavior of modeling results in Figure 3 are robust across a wide range of parameters. The parameters used are based on parameters in Figure 2 with slight adjustment within physiological range. Parameters in Figures 4, 5 are the same, which are chosen to fit the plateau height, the peak response and peak time simultaneously. In all the numerical simulations of Equations (1)–(5), all the variables except x and p start from 0. x and p, according to their biological meanings introduced, all start from 1.
Figure 2
Table 1
| Parameter | Meaning | DL5 | VM7 |
|---|---|---|---|
| τE | Time constant for excitatory synapse | 50 ms | 50 ms |
| ωEE | Synaptic weight from ORN to PN | 160 nS | 105 nS |
| ωIE | Synaptic weight from ORN to LN | 10 nS | 10 nS |
| ρ | Intrinsic strength for presynaptic inhibition | 1.9 ms | 2.5 ms |
| U | Increase in release probability for faciliation | 0.31 | 0.24 |
| τD | Time constant for STD | 368 ms | 160 ms |
| τF | Time constant for STF | 339 ms | 150 ms |
Model parameters used in Figure 2.
Figure 3
Figure 4
Figure 5
3. Results
3.1. STP Is Crucial for the Observed Nonlinear Divisive Normalization
A PN's response is suppressed by the firing of LNs, which can be activated by many ORNs, some of which can respond to a public odorant other than the cognate (private) odorant (see Figure 1). This introduces lateral inhibition and reduces the response of PN to a private odorant in the presence of public odorants. In Olsen et al. (), PN responses to a private odorant (which activates only the cognate ORN) with a background of different concentrations of a public odorant (presumably activates many other ORNs) were measured. As shown in Figure 2 (symbols are from experiments), the PN response to the private (cognate) odorant is affected by the level of public odorant (each color represents a different concentration of the public odorant). This divisive normalization effect was fitted phenomenologically by a Hill function in Olsen et al. (): where is defined as the average response over the first 500ms after the stimulus is applied. Note that is different from the adapted (steady state) response due to the finite adaptation time of PNs. In Equation (11), σ is proportional to the sum of the firing rates of non-cognate (public) ORNs; Rmax is the maximum PN firing rate; K1/2 is the firing rate of ORN at which PN has the half maximum response when σ = 0, γa is the Hill coefficient for the average response. Olsen et al. () were able to fit their experimental data with γa ≈ 1.5, which indicates a strong nonlinear effect in divisive normalization (γa > 1). However, the underlying mechanism of this “nonlinear” divisive normalization remains unclear.
Here, we use the circuit model to explain the experimentally observed nonlinear divisive normalization behaviors. By numerically solving our model (Equations 1–5), we computed the 500 ms average response the same way as in the experiments for different public odorant backgrounds. As shown in Figure 2, our model fits well with the measured responses of both PNs that innervate either DL5 or VM7 glomeruli in different backgrounds of public odorant concentrations that are fitted as in the model. The fitted values of are found to be linearly proportional to the measured local field potentials in experiments (Olsen et al., ; see Supplementary Figure 1), which further supports the validity of our model. Other best-fit parameters are listed in Table 1. Most of the parameters for DL5 and VM7 remain approximately the same, however, the two STP timescales (τD and τF) and the lateral inhibition strength ρ, which are intrinsic properties of the specific glomerulus, are different. Quantitatively, both STP timescales (τD and τF) obtained from our model fitting are consistent with the range of these timescales measured in experiments (Abbott et al., ). We also find that the sensitivity to lateral inhibition (ρ) is slightly stronger of VM7 than that of DL5, which is also consistent with the experiments (Olsen et al., ).
The most interesting finding from fitting our model with the experimental data is that both DL5 and VM7 have strong STP effects (SF, SD > 1), and both STF and STD strength are stronger in DL5 than those in VM7. However, the relative strength between STD and STF, r ≡ SD/SF, remains roughly the same for DL5 (r ≈ 1.08) and VM7 (r ≈ 1.07). Thus, our results suggest that strong and balanced STD and STF effects are responsible for the observed nonlinear divisive normalization in both VM7 and DL5. In fact, our model fails to fit experimental data without STP. As shown in Supplementary Figure 3, when we set τD = τF = 0, the PN response curves are roughly linear to ORN input within the range of the experimental data, which confirms the crucial role of STP in the observed “nonlinear” divisive normalization behavior.
To better understand this “nonlinear” effect in divisive normalization, we leveraged the simplicity of our model to derive an analytical expression for the steady state response of PNs. For any constant input (cognate ORN with firing rate R) to the antennal lobe, the output , i.e., the firing rate of the cognate PN, can be determined analytically by solving the steady state in Equations (1–5): where with , and with denoting the total input from public ORNs.
In the absence of STP, i.e., when τD → 0 and τF → 0, x = 1 and u+ = U remain constant, the PN response (Equation 12) reduces to: where both the numerator and the denominator depend linearly on R. We thus refer to Equation (13) as “linear” divisive normalization. In the presence of STP, i.e., when τF ≠ 0 and τD ≠ 0, both the denominator and the numerator in Equation (12) are nonlinear in R, which leads to the “nonlinear” divisive normalization behavior.
Similar to the empirical Hill function (11), the above steady state PN response curve exhibits a sigmoidal shape that can be characterized by three parameters: the maximum response ; the half maximum input R1/2 defined as ; and an effective Hill coefficient . From (13), we have the linear divisive normalization behavior (γ = 1) in the absence of STP. In the presence of STP, we have γ > 1, which can be used to characterize the nonlinearity of the response. In Figures 3A,B, the PN response function (12) with different STP strengths are shown. As expected, for a given input R, STF enhances the response while STD suppresses it, and the latter has a stronger effect. The dependence of γ on the STP strengths SD and SF are shown in Figure 3C. Interestingly, the Hill coefficient γ is enhanced by both STD and STF. Larger values of γ are reached by having roughly the same STD and STF strengths (dotted line in Figure 3C). A careful comparison between the Hill coefficient γa obtained using 500 ms average response of PNs (as in the experiments) and the effective Hill coefficient γ from the steady state response showed that they are highly correlated and both show similar dependence on SF and SD (Supplementary Figure 2). Thus, the steady state analysis confirms the crucial role of STP in generating the experimentally observed nonlinear divisive normalization.
3.2. Both STP and PI Control the Adaptive Responses to Time-Varying Stimuli
Odors in the environment are highly intermittent and dynamic (Murlis et al., ; Vickers et al., ). The ability to detect and respond to temporal variation of odor stimuli is crucial for the survival of many animals. Kim et al. () studied the responses of PNs to different time-dependent ORN signals that follow triangle-shaped temporal patterns with different peak times. The responses of PNs were “asymmetric,” with a faster rising phase, followed by a plateau and a slower decaying phase, depending on the rate of change in ORN's firing rate profile (see Figure 3 in Kim et al., ). Here, we use our model to explain the response patterns to these time-dependent signals.
We start by considering the response of PNs to a signal that increases linearly with time, i.e., R(t) = Kt. In the short time limit t ≪ τp, PI is negligible so the reduction factor p ≈ 1 (see Equation 3). As a result, the response is linearly proportional to the input . In the long time limit t ≫ τp, τE, p ≈ (AR(t))−1 decreases inversely proportional to the input R(t) which reduces the effective input pR(t). In fact, the effective input can be approximated as (see section 2 for detailed derivation): In the long time limit t ≫ τE, τp, the system adapts by adjusting the inhibition factor p so that the effective input reaches a constant as t → ∞, which is independent of the input. The corresponding adapted response can thus be determined analytically: which takes exactly the same form as the response to an effective time-independent signal as in (12). This is supported by numerical simulation as shown in Figure 4. The effect of PI is crucial for canceling out the increasing signal, resulting in a constant effective input . From Equation (15), it is clear that STP controls (modulates) the adapted response in the same way as it affects the response to a constant signal, i.e., STF enhances the adapted response and STD suppresses it.
We now study the PN responses to triangle-shaped input signals similar to those used in experiments (Kim et al., ) with our model. As shown in Figure 5A, the general response dynamics follows closely the experimentally observed behaviors. During the rising phase of the input signal, the PN response activity reaches its peak in a timescale that depends on the rate of change (K) of the input signal. For small K shown as the purple line in Figure 5A, the PN response reaches a plateau before the input reaches its peak due to the adaptive effect of presynaptic inhibition described above. The shaded region in Figure 5A (lower panel) shows the range of the experimentally observed plateau consistent with the model result (purple line). Note that for higher K the plateau activities depend on the rates of change in the input signals because they do not have enough time to reach the adapted value. In the descending phase of the input signal, the PN activity decreases following the input.
Quantitatively, the PN response can be described by two parameters: the peak response, which is defined as the PN activity at the peak time of the input signal, and the time to reach the peak response. In Figures 5B,C, we show that the results for these quantities measured from our model are in excellent agreement with those obtained from experiments (Kim et al., ).
Interestingly, STP alone can also lead to adaptive responses. As shown in Figure 1B, aside from the IFFL in which the negative (inhibitory) arm is formed via presynaptic inhibition (PI) by LN (ORN → LN → PN), there is another IFFL in the circuit with the negative (inhibitory) feedforward arm formed via STD (ORN → x → PN), which can achieve adaptive response even in the absence of PI. Quantitatively, for the experimental system studied by Kim et al. (), we find that the adaptive response is more strongly affected by PI. This can be seen by comparing PN's response to triangle-shaped input in three model variants: the standard model used in Figure 5 with strong PI and moderate STP, a model with PI only and without STP (SD = SF = 0), and a model with only STP without PI (ρ = 0). As shown in Supplementary Figure 4, the model response agrees with experiments qualitatively well even with PI alone, however, STP is required to achieve quantitative agreement with experiments. On the other hand, even though a very strong STD may also cause the system to exhibit adaptive response (Supplementary Figure 4B), the fit of the STP-only model to experimental data remains relatively poor even with fine tuning of the STF and STD strengths (see Supplementary Figures 4B–D for details).
Our results show that both STD and PI control the adaptive response albeit via different routes for the feed-forward inhibition and both of them are needed to explain the experimental results quantitatively. Given their different timescales and the fact that the PI mediated inhibition occurs upstream of the inhibition caused by STD, the more dominant role of PI suggested by our analysis seems reasonable. However, we believe that the double IFFL design in the AL circuit may represent a network architecture that improves the robustness of the system and STP is needed to modulate the adapted activity [see Equation (15)].
4. Discussion
In this study, we developed a simple neural circuit model for the antennal lobe of Drosophila and systematically studied the role of STP in early odor information processing. Combining analytical derivations of a steady-state solution of the model and numerical simulations, we showed that the model can capture key characteristics of PNs' responses to different ORN inputs, in particular divisive normalization and adaptive responses to time-varying signals. Comparison with experimental results revealed that STP is crucial for the observed nonlinear divisive normalization. We found that there are two IFFLs (one via STP and the other via PI) in the AL circuit (Figure 1B) and either of the two IFFLs or together can lead to adaptive response. We speculate that this sequential IFFL design in the AL circuit may be desirable for a robust adaptive behavior. Since both divisive normalization and adaptive response to time-varying inputs are omnipresent in nervous systems, the effects of STP identified here may also apply to other sensory modalities and brain regions.
Previous experiments have suggested that STD largely determines the nonlinear response function of PNs (Kazama and Wilson, ), and PI plays an important role in divisive normalization (Olsen and Wilson, ; Olsen et al., ). Our model extends these studies and showed that STP is crucial for the observed nonlinear divisive normalization. For stable ORN input, STF enhances PN's response while STD suppresses it. Yet, both STF and STD enhance the nonlinearity of the divisive normalization (Figure 3). Although the transient response of PNs to step-like ORN input has been attributed to STD, our model shows that PI predominately determines PN's response to more dynamic ORN input, such as triangle-shaped input patterns (Figure 5). In fact, our model can predict the response properties of PNs to arbitrary time-dependent ORN inputs. For example, when applied to a set of sine-wave ORN inputs with different frequencies, our model predicts that the amplitude of PN's oscillatory response will increase as the frequency gets higher, while the time advance of PN response peak to ORN response peak will decrease (see Supplementary Figure 5). These predictions can be tested by future experiments.
Our model not only reveals the underlying mechanisms for the observed nonlinear divisive normalization behavior and adaptive responses to time-varying signals, it also provides a general framework for relating the microscopic properties of the system such as time scales and strengths of STP and PI to macroscopic behaviors such as the input-output relation. As demonstrated in this work in the cases of VM7 and DL5 glomeruli, we can use our model to infer microscopic properties of the system from the measured input-output relation. More specifically, from our model study, the timescales of both STD (τD) and STF (τF) are predicted to be longer in DL5 than those in VM7, which can be verified by experiments (Abbott et al., ). The model-based analysis of the input-output response can be extended to other glomeruli. Our model can also be used to make predictions for changes in the input-output relation when certain microscopic properties, e.g., the STP strength (τF and τD) or the PI strength (ρ) are perturbed. These predictions can be tested in future experiments.
As we focused on building a minimal model to understand the underlying mechanism for nonlinear divisive normalization and adaptive response, we have made several simplifications in our study. First, our model is a mean-field rate model which neglects the noise. The nonlinear divisive normalization effect is defined based on steady state firing rate of PN (Equation 12), which only depends on the averaged values of variables in Equations (1)–(5). Additive noise (with 0 mean) in the neural dynamics of PN and LN will be averaged out and does not affect the steady state firing of PNs. For the adaptive response of PNs to time-varying stimuli, additive noise does not change the qualitative behavior of peak time, peak firing rate and plateau firing rate of PNs (Supplementary Figure 7).
Second, at the synapse level, our model ignores the STP effect at the ORN-LN synapses (Nagel and Wilson, ), which can affect the LN response especially its temporal dynamics. In our circuit model, we only considered one type of LNs for simplicity. There are several types of LNs in the antennal lobe with diverse innervation patterns and physiological properties (Chou et al., ; Nagel et al., ; Nagel and Wilson, ). For example, a small fraction of LNs are excitatory. Panglomerular LNs innervate all glomeruli and have higher spontaneous firing rates than other LNs. They are inhibited or only weakly excited by odors. Such inhibition of panglomerular LNs tends to dis-inhibit the entire AL in the presence of odors (Chou et al., ). LNs also show distinct response dynamics to odors. This temporal diversity of LNs likely shape the dynamics of PNs to complex ecologically relevant odor stimuli, such as the response time and the synchronization of spikes across PNs (Tanaka et al., ; Nagel and Wilson, ). Flies can sense temporal features of odor plume to help navigation (Álvarez-Salvado et al., ; Demir et al., ). Future study should address how the additional complexity of AL circuit affect PN's response properties and contribute to the olfactory behavior (Kao and Lo, ). Furthermore, we only considered PNs that innervate a single glomerulus (uPNs), some PNs do receive input from multiple glomeruli (mPNs). uPNs and mPNs have different projection patterns and may carry different aspects of odor information to the higher brain regions such as lateral horn and mushroom body (Bates et al., ). Future studies that incorporate these important features (both at the synapse level and the network level) will further our understanding of the rich dynamics in the early olfactory information processing in Drosophila.
Funding
The work was supported by the National Natural Science Foundation of China (Grant Nos. 12090053 and 32088101). The work by YT was supported by a NIH grant (R35GM131734).
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author contributions
YL, QL, SQ, and YT contributed to the study design and prepared the manuscript. YL performed all simulations. CT and YT supervised the study. All authors read and approved the final manuscript.
Acknowledgments
We thank Guangwei Si, Tianyi Wu, and members of Tang lab for helpful discussion.
Conflict of interest
YT was employed by IBM T. J. Watson Research Center. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncom.2021.730431/full#supplementary-material
References
1
AbbottL. F.VarelaJ.SenK.NelsonS. (1997). Synaptic depression and cortical gain control. Science275, 221–224. 10.1126/science.275.5297.221
2
Álvarez-SalvadoE.LicataA. M.ConnorE. G.McHughM. K.KingB. M.StavropoulosN.et al. (2018). Elementary sensory-motor transformations underlying olfactory navigation in walking fruit-flies. Elife7:e37815. 10.7554/eLife.37815
3
BatesA. S.SchlegelP.RobertsR. J.DrummondN.TamimiI. F.TurnbullR.et al. (2020). Complete connectomic reconstruction of olfactory projection neurons in the fly brain. Curr. Biol. 30, 3183–3199. 10.1016/j.cub.2020.06.042
4
BhandawatV.OlsenS. R.GouwensN. W.SchliefM. L.WilsonR. I. (2007). Sensory processing in the drosophila antennal lobe increases reliability and separability of ensemble odor representations. Nat. Neurosci. 10, 1474–1482. 10.1038/nn1976
5
CarandiniM.HeegerD. J. (2012). Normalization as a canonical neural computation. Nat. Rev. Neurosci. 13, 51–62. 10.1038/nrn3136
6
ChouY.-H.SpletterM. L.YaksiE.LeongJ. C.WilsonR. I.LuoL. (2010). Diversity and wiring variability of olfactory local interneurons in the drosophila antennal lobe. Nat. Neurosci. 13:439. 10.1038/nn.2489
7
DavidC.KennedyJ.LudlowA. (1983). Finding of a sex pheromone source by gypsy moths released in the field. Nature303, 804–806. 10.1038/303804a0
8
DayanP.AbbottL. F. (2001). Theoretical Neuroscience: Computational and Mathematical Modeling of Neural Systems. Computational Neuroscience Series. Cambridge, MA: The MIT Press.
9
DemirM.KadakiaN.AndersonH. D.ClarkD. A.EmonetT. (2020). Walking drosophila navigate complex plumes using stochastic decisions biased by the timing of odor encounters. Elife9:e57524. 10.7554/eLife.57524
10
DittmanJ. S.KreitzerA. C.RegehrW. G. (2000). Interplay between facilitation, depression, and residual calcium at three presynaptic terminals. J. Neurosci. 20, 1374–1385. 10.1523/JNEUROSCI.20-04-01374.2000
11
FergusonK. A.CardinJ. A. (2020). Mechanisms underlying gain modulation in the cortex. Nat. Rev. Neurosci. 21, 80–92. 10.1038/s41583-019-0253-y
12
GerstnerW.KistlerW. M.NaudR.PaninskiL. (2014). Neuronal Dynamics: From Single Neurons to Networks and Models of Cognition. Cambridge: Cambridge University Press. 10.1017/CBO9781107447615
13
HallemE. A.CarlsonJ. R. (2006). Coding of odors by a receptor repertoire. Cell125, 143–160. 10.1016/j.cell.2006.01.050
14
HanssonB. S.StensmyrM. C. (2011). Evolution of insect olfaction. Neuron72, 698–711. 10.1016/j.neuron.2011.11.003
15
KaoK.-W.LoC.-C. (2020). Short term depression, presynaptic inhibition and local neuron diversity play key functional roles in the insect antennal lobe. J. Comput. Neurosci. 48, 213–227. 10.1007/s10827-020-00747-4
16
KazamaH.WilsonR. I. (2008). Homeostatic matching and nonlinear amplification at identified central synapses. Neuron58, 401–413. 10.1016/j.neuron.2008.02.030
17
KimA. J.LazarA. A.SlutskiyY. B. (2011). System identification of drosophila olfactory sensory neurons. J. Comput. Neurosci. 30, 143–161. 10.1007/s10827-010-0265-0
18
KimA. J.LazarA. A.SlutskiyY. B. (2015). Projection neurons in drosophila antennal lobes signal the acceleration of odor concentrations. Elife4:e06651. 10.7554/eLife.06651
19
LuoS. X.AxelR.AbbottL. (2010). Generating sparse and selective third-order responses in the olfactory system of the fly. Proc. Natl. Acad. Sci. U.S.A. 107, 10713–10718. 10.1073/pnas.1005635107
20
MaW.TrusinaA.El-SamadH.LimW. A.TangC. (2009). Defining network topologies that can achieve biochemical adaptation. Cell138, 760–773. 10.1016/j.cell.2009.06.013
21
MarkramH.WangY.TsodyksM. (1998). Differential signaling via the same axon of neocortical pyramidal neurons. Proc. Natl. Acad. Sci. U.S.A. 95, 5323–5328. 10.1073/pnas.95.9.5323
22
MartelliC.FialaA. (2019). Slow presynaptic mechanisms that mediate adaptation in the olfactory pathway of drosophila. eLife8:e43735. 10.7554/eLife.43735
23
MasseN. Y.TurnerG. C.JefferisG. S. (2009). Olfactory information processing in drosophila. Curr. Biol. 19, R700–R713. 10.1016/j.cub.2009.06.026
24
MurlisJ.ElkintonJ. S.CardeR. T. (1992). Odor plumes and how insects use them. Annu. Rev. Entomol. 37, 505–532. 10.1146/annurev.en.37.010192.002445
25
NagelK. I.HongE. J.WilsonR. I. (2015). Synaptic and circuit mechanisms promoting broadband transmission of olfactory stimulus dynamics. Nat. Neurosci. 18:56. 10.1038/nn.3895
26
NagelK. I.WilsonR. I. (2016). Mechanisms underlying population response dynamics in inhibitory interneurons of the drosophila antennal lobe. J. Neurosci. 36, 4325–4338. 10.1523/JNEUROSCI.3887-15.2016
27
OlsenS. R.BhandawatV.WilsonR. I. (2010). Divisive normalization in olfactory population codes. Neuron66, 287–299. 10.1016/j.neuron.2010.04.009
28
OlsenS. R.WilsonR. I. (2008). Lateral presynaptic inhibition mediates gain control in an olfactory circuit. Nature452, 956–960. 10.1038/nature06864
29
SaitoH.ChiQ.ZhuangH.MatsunamiH.MainlandJ. D. (2009). Odor coding by a mammalian receptor repertoire. Sci. Signal. 2:ra9. 10.1126/scisignal.2000016
30
Shen-OrrS. S.MiloR.ManganS.AlonU. (2002). Network motifs in the transcriptional regulation network of Escherichia coli. Nat. Genet. 31, 64–68. 10.1038/ng881
31
SiG.KanwalJ. K.HuY.TaboneC. J.BaronJ.BerckM.et al. (2019). Structured odorant response patterns across a complete olfactory receptor neuron population. Neuron101, 950–962. 10.1016/j.neuron.2018.12.030
32
StevensC. F.WangY. (1995). Facilitation and depression at single central synapses. Neuron14, 795–802. 10.1016/0896-6273(95)90223-6
33
SuC.-Y.MenuzK.CarlsonJ. R. (2009). Olfactory perception: receptors, cells, and circuits. Cell139, 45–59. 10.1016/j.cell.2009.09.015
34
TanakaN. K.ItoK.StopferM. (2009). Odor-evoked neural oscillations in drosophila are mediated by widely branching interneurons. J. Neurosci. 29, 8595–8603. 10.1523/JNEUROSCI.1455-09.2009
35
TrappenbergT. (2009). Fundamentals of Computational Neuroscience. Oxford: Oxford University Presss.
36
TsodyksM.PawelzikK.MarkramH. (1998). Neural networks with dynamic synapses. Neural Comput. 10, 821–835. 10.1162/089976698300017502
37
TuY.RappelW.-J. (2018). Adaptation in living systems. Annu. Rev. Condens. Matter Phys. 9, 183–205. 10.1146/annurev-conmatphys-033117-054046
38
UchidaN.PooC.HaddadR. (2014). Coding and transformations in the olfactory system. Annu. Rev. Neurosci. 37, 363–385. 10.1146/annurev-neuro-071013-013941
39
VickersN. J.ChristensenT. A.BakerT. C.HildebrandJ. G. (2001). Odour-plume dynamics influence the brain's olfactory code. Nature410, 466–470. 10.1038/35068559
40
VictorJ. D.BoieS. D.ConnorE. G.CrimaldiJ. P.ErmentroutG. B.NagelK. I. (2019). Olfactory navigation and the receptor nonlinearity. J. Neurosci. 39, 3713–3727. 10.1523/JNEUROSCI.2512-18.2019
41
WangY.MarkramH.GoodmanP. H.BergerT. K.MaJ.Goldman-RakicP. S. (2006). Heterogeneity in the pyramidal network of the medial prefrontal cortex. Nat. Neurosci. 9, 534–542. 10.1038/nn1670
42
WarkB.LundstromB. N.FairhallA. (2007). Sensory adaptation. Curr. Opin. Neurobiol. 17, 423–429. 10.1016/j.conb.2007.07.001
43
WilsonR. I. (2013). Early olfactory processing in drosophila: mechanisms and principles. Annu. Rev. Neurosci. 36, 217–241. 10.1146/annurev-neuro-062111-150533
44
ZhangD.LiY.WuS. (2013). Concentration-invariant odor representation in the olfactory system by presynaptic inhibition. Comput. Math. Methods Med. 2013:507143. 10.1155/2013/507143
45
ZhangD.WuS.RaschM. J. (2015). Circuit motifs for contrast-adaptive differentiation in early sensory systems: the role of presynaptic inhibition and short-term plasticity. PLoS ONE10:e118125. 10.1371/journal.pone.0118125
Summary
Keywords
short-term plasticity, presynaptic inhibition, circuit model, olfactory system, Drosophila, divisive normalization
Citation
Liu Y, Li Q, Tang C, Qin S and Tu Y (2021) Short-Term Plasticity Regulates Both Divisive Normalization and Adaptive Responses in Drosophila Olfactory System. Front. Comput. Neurosci. 15:730431. doi: 10.3389/fncom.2021.730431
Received
24 June 2021
Accepted
23 September 2021
Published
22 October 2021
Volume
15 - 2021
Edited by
Pulin Gong, The University of Sydney, Australia
Reviewed by
Wei Lin, Fudan University, China; Matthias H. Hennig, University of Edinburgh, United Kingdom
Updates
Copyright
© 2021 Liu, Li, Tang, Qin and Tu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shanshan Qin ssqin@seas.harvard.eduYuhai Tu yuhai@us.ibm.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.