|
|
||||||||
1Department of Neurobiology and Behavior, Cornell University, Ithaca, New York; 2The Rockefeller University, New York, New York; and 3Department of Entomology, University of California, Riverside, California
Submitted 15 December 2004; accepted in final form 19 January 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
| METHODS |
|---|
|
|
|---|
Similar to our previous models of the honeybee AL (Linster and Masson 1996
; Linster and Smith 1997
), the model incorporates olfactory sensory neurons (OSNs), AL interneurons (LNs), and AL projection neurons (PNs). The evolution of each neuron's membrane potential v(t) around its resting potential is described by a first-order differential equation
![]() | (1) |
is the charging time constant of the neuron and Iext(t) is the total input current at time t. For OSNs, the net total input is directly proportional to the olfactory stimulation that they receive. For LNs and PNs, the input from a particular synapse at time t is computed as a function of the connection strength wij between the presynaptic neurons j and the postsynaptic neuron i, the conductance change g(t) due to a presynaptic event xj(t0) at time t0, and the difference between the Nernst potential EN,ij of the associated channel type and the current membrane potential vi(t) of the postsynaptic neuron
![]() | (2) |
![]() | (3) |
![]() | (4) |
min is the spiking threshold and
max is the saturation threshold (at which the maximum spike rate occurs). After any occurrence of a spike, the membrane potential was reset to the resting potential. The differential equations were iterated using exponential Euler with a time step of 1 ms. Model architecture
The computational model of the honeybee AL presented here was designed to incorporate and further investigate results from calcium imaging obtained in the AL. The basic architecture of the model, designed to replicate some of the main known features of AL anatomy and physiology, has been adapted from our previous models (Cleland and Linster 2002
; Linster et al. 1994
; Linster and Cleland 2001
; Linster and Masson 1996
; Linster and Smith 1997
). One of the major adjustments made in the present version of the model is the simulation of "real" input patterns for each odor stimulus as directly extracted from calcium-imaging data. Consequently, we could simulate real odorant response patterns instead of constructing fictional odorant representations in a hypothetical multi-dimensional space as is usually done. Existing data-derived features of the model, e.g., its field-oscillatory properties and their dependence on inhibitory interactions, have been carefully preserved in this version (Cleland and Linster 2002
; Stopfer et al. 1997
). We here focus on the question of how inhibitory connections should be organized to best replicate the measured input-output function of the AL (Sachse and Galizia 2003
). The model serves to test how the transformation from input (glomerular input pattern) to output (glomerular output pattern), as measured in calcium-imaging experiments, can best be achieved by the AL neural network.
In the honeybee AL, all synaptic interactions are located in areas of high synaptic density called glomeruli (Gascuel and Masson 1991
). In the present model, the number of these glomeruli is reduced from 160 to 20, corresponding to the number of glomeruli that could be visualized and for which calcium-imaging data had been obtained from both the input and output (PNs) in response to a common set of odorants (Sachse and Galizia 2002
, 2003
; Sachse et al. 1999
) (see Fig. 1A for detailed description of the model architecture). Each model glomerulus received input from one sensory neuron (OSN, representing the summed activity of many convergent olfactory sensory neurons) and contained one local interneuron which receives excitatory input in that glomerulus and inhibits all other glomeruli [called heteroLN by Fonta et al. (1993)
, Fig. 1B], also representing the summed activity of many functionally attuned neurons and one PN. A second class of local interneurons [called homoLN by Fonta et al. (1993)
] is included in the model; these homoLNs receive input in all glomeruli, project heavily onto themselves, and inhibit all other LNs and all PNs in a homogeneous manner (not shown in Fig. 1A). The feedback connections of these homoLNs onto themselves are likely to provide the oscillatory dynamics that have been observed experimentally (Stopfer et al. 1997
). When stimulated with simulated odorants, the AL network exhibited oscillatory dynamics similar to those described experimentally and PN spiking was phase-locked with the population oscillation. The oscillatory dynamics and phase-locking in the model are dependent on homoLN feedback inhibitory connections (Cleland and Linster 2002) and therefore disappear when these inhibitory interactions are blocked as was shown experimentally when the GABAergic antagonist picrotoxin was introduced into the AL (MacLeod and Laurent 1996
; Stopfer et al. 1997
).
|
|
|
| RESULTS |
|---|
|
|
|---|
|
|
The best results were obtained when we created an interglomerular projection pattern based on functional correlations between glomerular input activity patterns to odorants (functional network, Figs. 1Ciii and 3C). Experimentally, glomerular input patterns were obtained in response to 60 odorants (Galizia et al. 1999
; Sachse et al. 1999
). We here considered the possibility that glomeruli would inhibit each other proportionally to the similarity of their olfactory response profiles, i.e., that a functional inhibitory network would exist. As a consequence, glomeruli activated by common groups of odorants would be more likely to inhibit each other than those activated by noncommon groups of odorants. To determine the possible inhibitory projection patterns, we used glomerular input activity data from the calcium imaging experiments in response to all 60 odors. The glomerular input patterns of the 20 glomeruli considered here were correlated pairwise, resulting in a single variable c (1 < c < 1) describing the olfactory response profile similarity between every pair of glomeruli. The positive values from this matrix were used directly to scale the inhibitory synaptic weights between pairs of glomeruli, whereas negative values were set to zero (Table 1); consequently, the inhibitory connections between any two glomeruli had synaptic strengths of c*w (with c*w
0). Using this established pattern of interglomerular inhibition, we varied w between w = 0.0 and w = 1.0 with a stepsize of 0.025 and ran 50 simulations for each value of w. As described in the preceding text, the average resulting PN activation patterns were then correlated with the average glomerular output patterns from the experimental data. The maximal average correlation of r = 0.90 between the simulated and experimentally measured PN responses to the three odor stimulations was obtained when the synaptic weight w was set to w = 0.15. The average values of r obtained over the range of synaptic weight values chosen ranged from 0.75 to 0.90. Figure 3C shows the comparisons between the experimentally measured input and output patterns and the simulated PN activation patterns in response to 1-hexanol at this best value of w. In comparison (see Fig. 4A), the correlations between simulated and experimentally measured output patterns resulting from the functional inhibition network were significantly higher than those obtained with the morphological or stochastic inhibitory network [morphological: F(1,4) = 30.11, P < 0.005; stochastic: F(1,4) = 18.4, P < 0.05], whereas the latter two were not significantly different from each other [F(1,4) = 0.07, P > 0.5].
|
Comparisons of odor activity patterns at the glomerular input and output of the AL have suggested that the interglomerular network serves to decrease the correlations among odor activation patterns (Sachse and Galizia 2002
) (Fig. 4B, compare experimental inputoutput), presumably facilitating behavioral discrimination among odorants. Indeed, pairs of odorants for which PN responses were more strongly decorrelated (Sachse and Galizia 2002
) were more readily discriminated by bees in a behavioral experiment than those for which PN responses were highly correlated (Laska et al. 1999
). In our simulations, only the interglomerular inhibitory connections based on glomerular response profiles reduced the correlations between pairs of odorants at the glomerular output level with respect to their correlations at the glomerular input level (Fig. 4B, functional inhibitionoutput). Interestingly, AL networks with stochastic inhibitory connections reproduced the correlations between input activation patterns at the level of PN output (Fig. 4B, stochastic inhibitioninput), whereas the network in which physically close glomeruli inhibited each other resulted in somewhat lower correlations between output activation patterns as compared with input activation patterns (Fig. 4B, morphological inhibition-output). These results indicate, as suggested previously (Linster and Cleland 2004
), that to ensure contrast enhancement among odorants with similar and overlapping glomerular activation patterns, interglomerular inhibition should be constructed as a function of glomerular response profiles rather than on their physical proximity.
We then tested how our three "best" networks dealt with data from a different set of experiments, in which hexanol, octanol, and nonanol were each presented at six different concentrations (Sachse and Galizia 2003
). These data were obtained under dissimilar conditions than those used in the previous simulations and are therefore ideally suited to test the generality of our results. To test how well the three "best" networks reproduced this new data set, we used the glomerular activation patterns in response to each of the 18 odorant stimulations (3 odorants at 6 different concentrations each) as input patterns to the AL network and compared the resulting PN activation patterns to those recorded experimentally. The network with functionally organized inhibition reproduced the experimentally measured input-output function significantly better than the two networks with stochastically and morphologically organized inhibition (functional compared with morphological P < 0.05; functional compared with stochastic P < 0.02; Fig. 4C). Interestingly, the average correlations with the experimentally recorded PN responses where best at intermediate concentrations (101 to 103), which corresponds to the working range of odor responses in PNs, i.e., to the steep part of their dose-response function. Correlations decreased considerably at very low concentrations (105), where odor-evoked patterns often consist of just a single activated glomerulus (Sachse and Galizia 2003
). Whether the good performance of the morphologically designed network at these low concentrations is functionally relevant remains to be elucidated.
| DISCUSSION |
|---|
|
|
|---|
The data that we used in the present study result from estimations of input activity derived from bath-applied dye (Sachse et al. 1999
), whereas output activity estimates are measured as postsynaptic calcium increases (Sachse and Galizia 2003
). Because the two techniques monitor different calcium responses, they cannot be compared directly. Rather, each response profile requires calibration across glomeruli as has been done here. We have avoided comparing input-output activity patterns directly; rather, we have determined the interglomerular network most likely to reproduce the experimentally observed transfer function by comparing simulated output with experimentally measured output.
We have previously argued that a positive functional correlation between inhibitory projection patterns and the underlying olfactory similarity space is necessary to render interglomerular inhibition an effective mechanism of contrast enhancement (Linster and Cleland 2004
). This is in contrast to lateral inhibition in the retina, for example, in which neurons with similar (spatial) receptive fields are also physically adjacent, enabling effective contrast enhancement in these two physical dimensions. In the AL/OB network, odorants activate a wide, dispersed array of olfactory glomeruli (Galizia and Menzel 2001
; Johnson et al. 1999
), hence, lateral inhibition (in the classical sense in which neurons with similar response profiles are morphologically adjacent and inhibit one another) cannot exist; as a consequence, inhibition between glomeruli would necessarily be widely dispersed and be functionally rather than morphologically organized. This raises the question of whether efficient processing within the AL can depend on the spatial position of olfactory glomeruli. This view has been strongly advocated for the mammalian OB, in which local interneurons may mediate a center-surround inhibitory network (Aungst et al. 2003
).
In many brain areas, specific spatial arrangements optimize processing and/or total wiring length, a phenomenon best explored within the visual system. In this context, the organization of the inhibitory network correlates with response profiles of the neurons in that neurons responding to similar stimuli tend to inhibit each other more than those not responding to similar stimuli. It is the organization rather than the overall strength of the inhibitory connections that determines the functionality of the inhibitory network. In many systems, such as the visual system, response profiles can be predicted from physical location so that the organization of inhibitory interactions correlates with both physical proximity and response profiles ("lateral inhibition"). In a multidimensional olfactory world, however, two dimensions are not sufficient; hence, neighborhood relationships may be relevant for subgroups of glomeruli but not as a general rule. The simulation studies presented in this work strongly suggest that inhibitory connections within the AL are not spatially constrained and are not local, uniform, or stochastic but rather that inhibitory connections within the AL are dictated by the functional response properties of glomeruli. That is, with increasing overlap in their molecular response profiles, glomeruli increase the strength of inhibitory connections between them. Given the high-dimensional nature of olfactory representations, we speculate that similar principles will hold in the vertebrate OB. Additional mechanisms, some of which have been described in the OB, may also be involved in differentiating olfactory stimuli. For example, intraglomerular excitatory mechanisms (Schoppa and Westbrook 2002
; Urban and Sakmann 2002
) may contribute to enhancing the response of the most strongly activated glomeruli. To date, no excitatory synaptic interactions mediated by local interneurons have been described in the honeybee AL.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: C. Linster, Cornell University, Mudd Hall, Ithaca, NY 14853 (E-mail: CL243{at}cornell.edu)
| REFERENCES |
|---|
|
|
|---|
Christensen TA, Pawlowski VM, Lei H, and Hildebrand JG. Multi-unit recordings reveal context-dependent modulation of synchrony in odor-specific neural ensembles. Nat Neurosci 3: 927931, 2000.[CrossRef][ISI][Medline]
Cleland TA and Linster C. How synchronization properties among second-order sensory neurons can mediate stimulus salience. Behav Neurosci 116: 212221, 2002.[CrossRef][ISI][Medline]
Duchamp-Viret P, and Duchamp A. GABAergic control of odour-induced activity in the frog olfactory bulb: possible GABAergic modulation of granule cell inhibitory action. Neuroscience 56: 905914, 1993.[CrossRef][ISI][Medline]
Faber T, Joerges J, and Menzel R. Associative learning modifies neural representations of odors in the insect brain. Nat Neurosci 2: 7478, 1999.[CrossRef][ISI][Medline]
Fonta C, Sun XJ, and Masson C. Morphology and spatial distribution of bee antennal lobe interneurons responsive to odours. Chem Senses 18: 101119, 1993.
Galizia CG, McIlwrath SL, and Menzel R. A digital three-dimensional atlas of the honeybee antennal lobe based on optical sections acquired by confocal microscopy. Cell Tissue Res 295: 383394, 1999.[CrossRef][ISI][Medline]
Galizia CG and Menzel R. Odor perception in honeybees: coding information in glomerular patterns. Curr Opin Neurobiol 10: 504510, 2000.[CrossRef][ISI][Medline]
Galizia CG and Menzel R. The role of glomeruli in the neural representation of odors: results from optical recording studies. J Insect Physiol 47: 115130, 2001.[CrossRef][ISI][Medline]
Gascuel J and Masson C. A quantitative ultrastructural study of the honeybee antennal lobe. Tissue Cell 23: 341355, 1991.[Medline]
Hansson BS, Carlsson MA, and Kalinova B. Olfactory activation patterns in the antennal lobe of the sphinx moth, Manduca sexta. J Comp Physiol [A] 189: 301308, 2003.
Joerges J, Kuettner A, Galizia G, and Menzel R. Representations of odors and odor mixtures visualized in the honeybee brain. Nature 387: 285288, 1997.
Johnson BA, Woo CC, Hingco EE, Pham KL, and Leon M. Multidimensional chemotopic responses to n-aliphatic acid odorants in the rat olfactory bulb. J Comp Neurol 409: 529548, 1999.[CrossRef][ISI][Medline]
Laska M, Galizia CG, Giurfa M, and Menzel R. Olfactory discrimination ability and odor structure-activity relationships in honeybees. Chem Senses 24: 42938, 1999.
Linster C and Cleland TA. How spike synchronization among olfactory neurons can contribute to sensory discrimination. J Comput Neurosci 10: 187193, 2001.[CrossRef][ISI][Medline]
Linster C and Cleland TA. Configurational and elemental odor mixture perception can arise from local inhibition. J Comput Neurosci 16: 3947, 2004.[CrossRef][ISI][Medline]
Linster C and Gervais R. Investigation of the role of interneurons and their modulation by centrifugal fibers in a neural model of the olfactory bulb. J Comput Neurosci 3: 225246, 1996.[CrossRef][ISI][Medline]
Linster C and Hasselmo M. Modulation of inhibition in a model of olfactory bulb reduces overlap in the neural representation of olfactory stimuli. Behav Brain Res 84: 117127, 1997.[CrossRef][ISI][Medline]
Linster C, Johnson B, Yue E, Morse A, Xu Z, Hingco E, Choi Y, Choi M, Mehissa A, and Leon M. Perceptual correlates of neural representations evoked by odorant enantiomers. J Neurosci 21: 98379843, 2001.
Linster C and Masson C. A neural model of olfactory sensory memory in the honey bee's antennal lobe. Neural Comput 8: 94114, 1996.
Linster C, Marsan D, Masson C, and Kerszberg M. Odor processing in the bee: a preliminary study of the role of central input to the antennal lobe. In: Advances in Neural Information Processing Systems 6, edited by Cowan JD, Tesauro G, and Alspector J. Morgan Kaufmann Publishers 1994, p. 527534.
Linster C and Smith BH. A computational model of the response of honey bee antennal lobe circuitry to odor mixtures: overshadowing, blocking and unblocking can arise from lateral inhibition. Behav Brain Res 87: 114, 1997.[CrossRef][ISI][Medline]
MacLeod K and Laurent G. Distinct mechanisms for synchronization and temporal patterning of odor-encoding neural assemblies. Science 274: 976979, 1996.
Malun D. Inventory and distribution of synapses of identified uniglomerular projection neurons in the antennal lobe of Periplaneta americana. J Comp Neurol 305: 348360, 1991a.[CrossRef][ISI][Medline]
Malun D. Synaptic relationships between GABA-immunoreactive neurons and an identified uniglomerular projection neuron in the antennal lobe of Periplaneta americana: a double-labeling electron microscopic study. Histochemistry 96: 197207, 1991b.[CrossRef][ISI][Medline]
Meister M and Bonhoeffer T. Tuning and topography in an odor map in the rat olfactory bulb. J Neurosci 15: 13511360, 2001.
Meredith M. Patterned response to odor in mammalian olfactory bulb: the influence of intensity. J Neurophysiol 56: 572597, 1986.
Ng M, Roorda RD, Lima SQ, Zemelman BV, Morcillo P, and Miesenbock G. Transmission of olfactory information between three populations of neurons in the antennal lobe of the fly. Neuron 36: 463474, 2002.[CrossRef][ISI][Medline]
Rubin BD and Katz LC. Optical imaging of odorant representations in the mammalian olfactory bulb. Neuron 23: 499511, 1999.[CrossRef][ISI][Medline]
Sachse S and Galizia CG. Role of inhibition for temporal and spatial odor representation in olfactory output neurons: a calcium imaging study. J Neurophysiol 87: 11061117, 2002.
Sachse S and Galizia CG. The coding of odour-intensity in the honeybee antennal lobe: local computation optimizes odor representation. Eur J Neurosci 18: 21192132, 2003.[CrossRef][ISI][Medline]
Sachse S, Rappert A, and Galizia CG. The spatial representation of chemical structures in the antennal lobe of honeybees: steps towards the olfactory code. Eur J Neurosci 11: 39703982, 1999.[CrossRef][ISI][Medline]
Schoppa NE and Westbrook GL. AMPA autoreceptors drive correlated spiking in olfactory bulb glomeruli. Nat Neurosci 5: 11941202, 2002.[CrossRef][ISI][Medline]
Stopfer M, Bhagavan S, Smith BH, and Laurent G. Impaired odor discrimination on desynchronization of odour-encoding neural assemblies. Nature 390: 7074, 1997.[CrossRef][Medline]
Urban NN. Lateral inhibition in the olfactory bulb and in olfaction. Physiol Behav 77: 607612, 2002.[CrossRef][Medline]
Urban NN and Sakmann B. Reciprocal intraglomerular excitation and intra- and interglomerular lateral inhibition between mouse olfactory bulb mitral cells. J Physiol 542: 355367, 2002.
Wachowiak M, Cohen LB, and Zochowski MR. Distributed and concentration-invariant spatial representations of odorants by receptor neuron input to the turtle olfactory bulb. J Neurophysiol 87: 10351045, 2002.
Wachowiak M, Zochowski M, Cohen LB, and Falk CX. The spatial representation of odors by olfactory receptor neuron input to the olfactory bulb is concentration invariant. Biol Bull 199: 162163, 2000.[ISI][Medline]
Wang JW, Wong AM, Flores J, Vosshall LB, and Axel R. Two-photon calcium imaging reveals an odor-evoked map of activity in the fly brain. Cell 112: 271282, 2003.[CrossRef][ISI][Medline]
Wilson RI, Turner GC, and Laurent G. Transformation of olfactory representations in the Drosophila antennal lobe. Science 303: 366370, 2004.
Xu F, Greer CA, and Shepherd GM. Odor maps in the olfactory bulb. J Comp Neurol 422: 489495, 2000.[CrossRef][ISI][Medline]
Yokoi M, Mori K, and Nakanishi S. Refinement of odor molecule tuning by dendrodendritic synaptic inhibition in the olfactory bulb. Proc Natl Acad Sci USA 92: 33713375, 1995.
This article has been cited by other articles:
![]() |
C. E. Reisenman, T. Heinbockel, and J. G. Hildebrand Inhibitory Interactions Among Olfactory Glomeruli Do Not Necessarily Reflect Spatial Proximity J Neurophysiol, August 1, 2008; 100(2): 554 - 564. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Schmuker and G. Schneider Processing and classification of chemical data inspired by insect olfaction PNAS, December 18, 2007; 104(51): 20285 - 20289. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. F. Silbering and C. G. Galizia Processing of Odor Mixtures in the Drosophila Antennal Lobe Reveals both Global Inhibition and Glomerulus-Specific Interactions J. Neurosci., October 31, 2007; 27(44): 11966 - 11977. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. C. Willhite, K. T. Nguyen, A. V. Masurkar, C. A. Greer, G. M. Shepherd, and W. R. Chen Viral tracing identifies distributed columnar organization in the olfactory bulb PNAS, August 15, 2006; 103(33): 12592 - 12597. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. A. Cleland and C. Linster Computation in the Olfactory System Chem Senses, November 1, 2005; 30(9): 801 - 813. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. I. Wilson and G. Laurent Role of GABAergic Inhibition in Shaping Odor-Evoked Spatiotemporal Patterns in the Drosophila Antennal Lobe J. Neurosci., October 5, 2005; 25(40): 9069 - 9079. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. E. Reisenman, T. A. Christensen, and J. G. Hildebrand Chemosensory Selectivity of Output Neurons Innervating an Identified, Sexually Isomorphic Olfactory Glomerulus J. Neurosci., August 31, 2005; 25(35): 8017 - 8026. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||