Chemical Senses Advance Access originally published online on November 2, 2005
Chemical Senses 2005 30(9):801-813; doi:10.1093/chemse/bji072
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Computation in the Olfactory System
Department of Neurobiology and Behavior, Cornell University, Ithaca, NY, USA
Correspondence to be sent to: Thomas A. Cleland, Department of Neurobiology and Behavior, Cornell University, Ithaca, NY, 14853 USA. e-mail: tac29{at}cornell.edu
| Abstract |
|---|
|
|
|---|
Computational models are increasingly essential to systems neuroscience. Models serve as proofs of concept, tests of sufficiency, and as quantitative embodiments of working hypotheses and are important tools for understanding and interpreting complex data sets. In the olfactory system, models have played a particularly prominent role in framing contemporary theories and presenting novel hypotheses, a role that will only grow as the complexity and intricacy of experimental data continue to increase. This review will attempt to provide a comprehensive, functional overview of computational ideas in olfaction and outline a computational framework for olfactory processing based on the insights provided by these diverse models and their supporting data.
Key words: coding, contrast enhancement, convergence, dynamics, modeling, normalization, odor representation, synchrony
| Introduction |
|---|
|
|
|---|
In natural environments, airborne chemical stimuli are distributed unpredictably in time and space, and odorants from innumerable sources intermix freely. The olfactory system must be able to detect potential signals of interest within these chemically noisy environments, correctly extract these signals from a complex and changing odor background to form stimulus representations, compare these constructed representations to those of previously experienced odors, differentiate relevant from irrelevant stimuli, and cue an appropriate response. Many of the neural circuit elements comprising the olfactory system have been proposed to contribute to these processes in particular ways; for example, multiple feedback and feedforward interactions among olfactory structures, as well as between olfactory and non-olfactory areas, are thought to contribute to the filtering and construction of olfactory representations. Computational models of olfactory processing have been increasingly utilized to describe and interpret these complex and interrelated phenomena.
Primary olfactory sensory neurons (OSNs) number in the millions in rodents. Their axons are highly convergent, targeting specific, discrete neuropilar synaptic regions within the input layer of the olfactory bulb (OB) called glomeruli. In hamsters, for example, between 1300 and 4700 OSNs expressing the same odorant receptor complement converge upon each glomerulus (Schoenfeld and Knott, 2004
). These large populations of redundant OSNs and their correspondingly high convergence ratios have been proposed to yield advantages such as an improved signal-to-noise ratio, a corresponding increase in effective stimulus sensitivity, and an increased range of tuning to different odorant concentrations (van Drongelen et al., 1978
; Duchamp-Viret et al., 1989
; Meisami, 1989
; Cleland and Linster, 1999
). The molecular receptive ranges or chemical receptive fields of these odorant receptors overlap substantially, such that the identity of odorants is not associated with the activation of a specific receptor but rather is represented by a distributed, combinatorial code (Adrian, 1953
; Moulton, 1967
; Stewart et al., 1979
; Kauer, 1991
), now recognized as a pattern of activation across many receptors. Owing to the specific homotypic convergence of OSN axons, these odor-specific activity patterns can be most clearly observed in imaging studies of OB glomeruli (Friedrich and Korsching, 1997
; Johnson et al., 1998
, 1999
, 2004
; Rubin and Katz, 1999
; Meister and Bonhoeffer, 2001
; Wachowiak et al., 2002
). These overlapping representations underlie two critical properties of the olfactory system that a labeled-line solution would not. First, the number of unique odor representations is not limited to the number of different receptor types (roughly 1000 in mice; Mombaerts, 1996
) but can be estimated as mn, where n denotes the number of receptor types and m denotes the number of recognizable states that each sensor can assume, ultimately limited by the signal-to-noise ratio of the system. Even if only two receptor states, active and inactive, were recognized, this would enable roughly 21000 potential odor stimuli to be discriminated in mice. Second, the fact that structurally and perceptually similar odorant molecules will activate correspondingly overlapping sets of olfactory receptors (ORs) (Linster et al., 2001b
, 2002
; Cleland et al., 2002
) establishes a basis for the recognition of stimulus similarity in the olfactory system. This is a prerequisite for basic postsensory cognitive processes such as generalization (Shepard and Chang, 1963
; Shepard, 1987
; Cleland et al., 2002
) and a tolerance for variance among repeated stimulus samples that a labeled-line system would have no clear means of generating.
Distributed patterns of activity in response to chemical stimuli are transmitted to the OB via OSN axons that terminate in the glomeruli of its input layer. The OB is believed to filter and transform these incoming sensory data, performing normalization, contrast enhancement, and similar operations before conveying the processed olfactory information to several different secondary olfactory structures via mitral cell axon collaterals (Cleland and Linster, 2003
). Notably, the bulb constitutes the last common stage at which olfactory sensory representations can be processed before the signal diverges dramatically into these multiple secondary structures. It is clear from recent investigations that the perceptual qualities of odorants can be predicted, to a limited degree, from the patterns of activation that they evoke at the OB input layer (Linster and Hasselmo, 1999
; Linster et al., 2001b
, 2002
; Cleland et al., 2002
). However, several aspects of odor perception, for example, changes in perception and discrimination capacity due to odor intensity or prior experience, cannot be predicted solely by this first-order representation as reflected in glomerular activation patterns. Nor is the converse true; mitral cell responses in behaving animals cannot be predicted solely by the odor(s) presented but depend substantially on odor contingency (Kay and Laurent, 1999
), as previously suggested in field potential recordings from the OB (Di Prisco and Freeman, 1985
; Gray et al., 1986
; Freeman and Grajski, 1987
; Grajski and Freeman, 1989
). Furthermore, the OB receives substantial centrifugal projections from both cortical and neuromodulatory centers, and its responses to odor presentations are strongly regulated by these centrifugal inputs as well as odor learning and experience (Kay et al., 1996
; Kay, 2003
, 2005
; Ravel et al., 2003
; Wilson and Stevenson, 2003
; Martin et al., 2004
; Wilson et al., 2004
). It is therefore safe to assume that the OB plays an important role in processing incoming sensory information. Accordingly, many models of OB signal processing have been developed, which are grouped here into studies of (1) filtering and contrast enhancement, (2) mechanisms underlying oscillations and spike synchronization, and (3) odor segmentation and associative memory function. In addition, a number of detailed biophysical models of bulbar neurons have been constructed, in many cases to address how their intrinsic properties underlie and interact with network properties.
| Filtering and contrast enhancement |
|---|
|
|
|---|
The high convergence ratio between OSNs expressing a particular odorant receptor and their target glomeruli in the OB (Figure 1A) is believed to improve the signal-to-noise ratio during odorant detection (van Drongelen et al., 1978
|
Contrast enhancement is a common property of sensory systems that narrows (sharpens) sensory representations by specifically inhibiting neurons on the periphery of the representation, thus enhancing the contrast between signal and background (Figure 2A). A number of computational models have investigated the contrast enhancement potential of OB circuitry, most of which, by analogy with the retina, have investigated the potential role of lateral inhibitory projections. Classically, bulb models have emphasized lateral inhibition mediated by mitral cell lateral dendrites (Rall and Shepherd, 1968
|
Contrast enhancement, the effects of which have been directly observed in the OB (Yokoi et al., 1995
| Mechanisms underlying oscillations and spike synchronization |
|---|
|
|
|---|
While recent models have begun to favor glomerular-layer mechanisms for contrast enhancement, granule cell activity also clearly shapes mitral cell response patterns and hence the presumptive odor representations that emerge from the OB. Specifically, several computational models of the OB have suggested that the temporal pattern of spiking among mitral cells may play a role in odor representation (Schild, 1988
|
One of the most widely studied features of OB processing has been the dynamic oscillatory activity patterns observed in the bulb in response to odor stimulation, particularly given the observation that mitral cell spiking is correlated in time with these field oscillations (Freeman and Grajski, 1987
Recently, models of these field oscillatory properties have begun to emphasize their relationship to the regulation of spike timing in mitral cells (Davison et al., 2003
; Margrie and Schaefer, 2003
), suggesting that patterns of spike synchronization among mitral cells responding to the same sensory input are important contributors to the odorant representation at this level. In the example shown in Figure 3C, the average firing rates of the two model mitral cells do not contribute to the differentiation of odors A and B. However, the phase of firing during oscillatory cycles does contribute odor identity information sufficient to discriminate these two odors. This type of coding scheme was first proposed for the insect antennal lobe by Laurent and co-workers (Laurent, 1996
; Laurent et al., 1996
; Stopfer et al., 1997
), supported by behavioral data demonstrating that a reduction in the synchronization of projection neuron action potentials impaired odor discrimination in honeybees (Stopfer et al., 1997
). This phenomenon was subsequently studied in genetically modified mice (Nusser et al., 2001
) and has been modeled by several groups. These models have primarily served to illustrate and emphasize the circuit properties and mechanisms by which the regulation of spike synchronization among secondary neurons could contribute to odor representations. Specifically, some models have explored the biophysical bases of slow temporal patterning and fast oscillogenesis among secondary olfactory neurons (Bazhenov et al., 2001
), while others have demonstrated how the regulation of spike synchronization among secondary neurons can influence the readout of information at the next level of processing (Perez-Orive et al., 2004
; Sivan and Kopell, 2004
), facilitate contrast enhancement, and underlie the salience of olfactory stimuli (Linster and Cleland, 2001
; Cleland and Linster, 2002
).
| Odor segmentation and associative memory function |
|---|
|
|
|---|
Odor segmentation is the general term for the problem of how the olfactory system is able to segregate and identify different odorants that are encountered simultaneously. As most odors comprise multiple separate odorant molecules, it is far from clear how the olfactory system can parse the multitude of odorant stimuli present at any given time and attribute each to appropriately separate sources. One approach has been to hypothesize that odors emitted from different sources can be segregated by OB circuitry based upon their differential fluctuations in time (Fort and Rospars, 1992
Olfactory associative memory functions have been more commonly attributed to the piriform cortex, one of the targets of mitral cell axons projecting from the OB. Specifically, the piriform cortex has been proposed to mediate the associative memory functions necessary for odor-context learning (Haberly and Bower, 1989
; Hasselmo et al., 1990
; Haberly, 2001
) and hierarchical clustering (Ambros-Ingerson et al., 1990
). For example, cortical short-term synaptic depression can be employed to filter out stable background odorants, whereas long-term synaptic plasticity can store associations between neurons responding to the same odor stimuli (Figure 4). Indeed, the extensive intrinsic feedback network in this cortex and its integration with afferent inputs closely resembles the structure of traditional theoretical associative memory networks as first described by Marr (1971)
. Several laboratories have constructed models of piriform cortex, implementing these associative memory functions (Haberly and Bower, 1989
; Ambros-Ingerson et al., 1990
; Hasselmo et al., 1990
; Barkai and Hasselmo, 1994
; Hasselmo et al., 1997
) as well as exploring the role of cortical cholinergic modulation in their regulation (Hasselmo et al., 1990
, 1997
; Patil and Hasselmo, 1999
; Linster and Hasselmo, 2001
; Linster et al., 2003
).
|
| Detailed biophysical models |
|---|
|
|
|---|
Most models of the olfactory system to date have emphasized network-level interactions and the properties of OB and piriform cortical circuitry using simplified cellular models conducive to these larger scale simulations. However, several relatively detailed biophysical models of olfactory neuronsparticularly OSNs and mitral cellshave also been constructed. Biophysical (or compartmental) models focus on membrane and cellular properties rather than network phenomena; to this end, they include substantial morphological, biophysical, and/or physiological detail. For example, membranes are modeled with membrane resistivity and capacitance values drawn from experimental data, with realistic or data-derived lengths, diameters, and morphological branching patterns along with models of passive and active ion channels inserted in appropriate regions of the model neuron. Biophysical models usually can be directly interrelated with electrophysiological data. However, they are often poor choices for large-scale models, due to the large numbers of weakly defined parameters as well as their substantial computational costs. Several detailed OSN models have illustrated how ligandreceptor binding and nonlinear transduction processes can underlie the experimentally observed response properties of these cells (Malaka et al., 1995
| Synthesis |
|---|
|
|
|---|
Computational models of the olfactory system have contributed immensely to the framing of experimental problems and the construction of complex hypotheses regarding its function. Here, we briefly outline a working hypothesis of olfactory system function, integrating the insights derived from the models reviewed above and constrained by their supporting data sets.
Inhaled odorants bind to specific populations of ORs expressed on the apical surface of primary OSNs, which are distributed across the nasal epithelium. ORs are broadly tuned for odorant ligands, such that even simple monomolecular odorants activate a number of receptor types to differing degrees, producing combinatorial patterns of activation among OR classes that reflect odor quality. While the pattern of OSN activation depends primarily upon the ligandreceptor affinities of multiple OR populations for the various molecular moieties (odotopes) of odors' component molecules, it also is likely influenced by other factors such as the fluid dynamics of inhalation, the net molecular sorptiveness of odorant molecules, and the behavioral regulation of odor sampling (reviewed in Schoenfeld and Cleland, 2005
). Other physical factors that affect the pattern of OSN activation include the concentration of odorants as well as interference caused by overlaps among the representations of multiple odotopes that are simultaneously presented to the olfactory epithelium. Different odotopes may compete for receptors for which they have not only different affinities but also different efficacies (Duchamp-Viret et al., 2003
; Araneda et al., 2004
; Oka et al., 2004
; Sanz et al., 2005
), such that a reduction in the activation of an OR genespecific OSN population may as well connote an increased concentration of a relatively antagonistic odotope as a reduced concentration of an agonist odotope. In short, the pattern of OSN activation in any natural scene context is likely to be an unknown composite of multiple, overlapping, and degraded primary odorant representations. It is from this unpromising raw material that the olfactory system must detect and identify relevant stimuli.
OSNs expressing the same OR and hence sharing the same molecular receptive range project their axons to specific glomeruli within the OB input layer (Mombaerts et al., 1996
). Hence, the primary olfactory representation can be conveniently measured by imaging glomeruli: that is, the axonal arborizations of convergent OSNs in the OB glomerular layer. Some glomerular imaging techniques are explicitly presynaptic (Friedrich and Korsching, 1997
; Wachowiak et al., 2002
), even those that are not are likely to reflect predominantly presynaptic activity due to the disproportionate number of OSNs arborizing within each glomerulus compared with the number of bulbar neurons arborizing therein (Shepherd and Greer, 1998
; Schoenfeld and Knott, 2004
); glomeruli, of course, do not contain cell bodies. Under simple experimental conditions, these glomerular response profiles are predictive of odor quality, as measured behaviorally (Johnson and Leon, 2000
; Linster et al., 2001b
; Cleland et al., 2002
; Leon and Johnson, 2003
). However, this concordance is not robust even to changes in odorant concentrations and certainly cannot be expected to persist in the context of a complex olfactory natural scene.
Mitral cells, along with middle and deep tufted cells, are the principal output neurons of the OB. They are directly postsynaptic to OSNs, and as such the pattern of mitral/tufted cell activation across the bulb constitutes the secondary olfactory representation. In mammals, mitral cells sample from only a single glomerulus, hence a glomerulus along with its associated mitral cells and interneurons has been referred to as an "odor column" (Shepherd and Greer, 1998
) that derives its receptive field primarily from a single population of convergent OSNs. However, the mitral/tufted activation pattern also depends on the activity of several classes of bulbar interneurons and hence is substantially transformed with respect to the primary olfactory representation. For example, some form of normalization of stimulus concentrations is clearly evident in the concentrationresponse profiles of mitral cells. Mitral cell activity patterns are relatively stable across concentrations compared with the changes in OSN responses (Chalansonnet and Chaput, 1998
), and when they are affected by odorant concentration changes, these changes are often complex and difficult to predict. Some mitral cells progress from excitation to inhibition with increasing concentrations, while others become more excited, often exhibiting shorter latencies to first spike; other cells display yet other profiles (Harrison and Scott, 1986
; Meredith, 1986
; Wellis et al., 1989
; Chalansonnet and Chaput, 1998
). What is clear, however, is that mitral cells do not monotonically increase their activity in response to increased stimulus intensities, nor, consequently, is there likely to be a substantial broadening of mitral cell activity across the bulb owing to the recruitment of lower affinity receptors. These are two essential ways in which the structure of the secondary (mitral) olfactory representation differs from that of the primary representation as measured among OSNs, and this transformation provides essential constraints to models of bulbar mechanisms.
Normalization, the process by which sets of values are rescaled to a common or tractable range, requires a negative feedback loop that is effectively global in scope; that is, the strength of feedback inhibition should be scaled to the average activity across the bulb rather than to local activation levels if the profiles of relative activation among odor columns are to be preserved. Contrast enhancement, in turn, requires the delivery of inhibition onto odor columns proportional to the activation of columns exhibiting similar receptive fields (molecular receptive ranges). In the OB, the synaptic triad connecting OSN arbors, PG cell dendritic spines, and mitral cell dendrites in close proximity (reviewed by Shepherd and Greer, 1998
) coupled with the lateral excitatory network mediated by external tufted and short axon cells (Aungst et al., 2003
; Hayar et al., 2004
) can effect both normalization and contrast enhancement between the primary and secondary olfactory representations (Cleland and Sethupathy, 2004
). Briefly, OSNs activate mitral cells, external tufted cells, and PG cell dendritic spines in parallel, and these PG spines deliver inhibition onto mitral cell apical dendrites, closely apposing the excitatory OSN inputs. Nontopographical contrast enhancement models (Cleland and Sethupathy, 2004
) predict that PG cellmediated feedforward inhibition of the mitral cell will predominate when the odor column is weakly activated, while direct OSN activation of the mitral cell will predominate when the odor column is strongly activated (Figure 2C), consistent with recordings from rabbit mitral cells (Yokoi et al., 1995
). Normalization is effected via the activation of external tufted cells, which proportionately excite a lateral excitatory network of external tufted and short axon cells, thereby broadly inhibiting mitral cell activity across the bulb via sign-inverting PG cells.
The net activation of mitral cells is then translated into trains of action potentials, the precise timing of which appears to be regulated by coordinated oscillations measurable in field recordings across the OB. These bulbar oscillations are thought to depend on an extensive excitatoryinhibitory network of mitral cell secondary dendrites and granule cell interneurons, as well as reciprocal cortical connections that modulate bulbar activity according to behavioral state (Kay, 2003
; Ravel et al., 2003
; Lagier et al., 2004
; Martin et al., 2004
). Mitral cell spike synchronization patterns can mediate a second level of feature extraction to the odor representation, as follower neurons in diverse central olfactory structures (Cleland and Linster, 2003
) process incoming mitral cell spikes via synaptic learning rules, the best known of which rely upon precise spike timing (Song et al., 2000
; Cleland and Linster, 2002
). Centrifugal modulatory inputs also influence olfactory processing and learning mechanisms within the OB (Sullivan et al., 2000
; Linster and Cleland, 2002
; Yuan et al., 2003
), emphasizing the active role of the OB in shaping and transforming odor signals and the importance of behavioral state.
| Conclusions |
|---|
|
|
|---|
Computational models have an established and growing role within systems neuroscience. As our understanding of neural processing and interactions becomes more sophisticated, computer models of these systems are increasingly necessary in order to understand and interpret experimental results. In the olfactory system in particular, computational modeling will no doubt be essential to understand the integration of the many factors influencing the construction and transformation of odor representations.
| References |
|---|
|
|
|---|
Adrian, E.D. (1953) Sensory messages and sensation; the response of the olfactory organ to different smells. Acta Physiol. Scand., 29, 514.[Medline]
Alkasab, T.K., White, J. and Kauer, J.S. (2002) A computational system for simulating and analyzing arrays of biological and artificial chemical sensors. Chem. Senses, 27, 261275.
Ambros-Ingerson, J., Granger, R. and Lynch, G. (1990) Simulation of paleocortex performs hierarchical clustering. Science, 247, 13441348.
Anton, P.S., Lynch, G. and Granger, R. (1991) Computation of frequency-to-spatial transform by olfactory bulb glomeruli. Biol. Cybern., 65, 407414.[CrossRef][ISI][Medline]
Araneda, R.C., Peterlin, Z., Zhang, X., Chesler, A. and Firestein, S. (2004) A pharmacological profile of the aldehyde receptor repertoire in rat olfactory epithelium. J. Physiol., 555(Pt 3), 743756.
Aungst, J.L., Heyward, P.M., Puche, A.C., Karnup, S.V., Hayar, A., Szabo, G. and Shipley, M.T. (2003) Centre-surround inhibition among olfactory bulb glomeruli. Nature, 426, 623629.[CrossRef][Medline]
Barkai, E. and Hasselmo, M.E. (1994) Modulation of the input/output function of rat piriform cortex pyramidal cells. J. Neurophysiol., 72, 644658.
Bazhenov, M., Stopfer, M., Rabinovich, M., Abarbanel, H.D., Sejnowski, T.J. and Laurent, G. (2001) Model of cellular and network mechanisms for odor-evoked temporal patterning in the locust antennal lobe. Neuron, 30, 569581.[CrossRef][ISI][Medline]
Bhalla, U.S. and Bower, J.M. (1993) Exploring parameter space in detailed single neuron models: simulations of the mitral and granule cells of the olfactory bulb. J. Neurophysiol., 69, 19481965.
Bozza, T., Feinstein, P., Zheng, C. and Mombaerts, P. (2002) Odorant receptor expression defines functional units in the mouse olfactory system. J. Neurosci., 22, 30333043.
Breakspear, M. (2001) Perception of odors by a nonlinear model of the olfactory bulb. Int. J. Neural. Syst., 11, 101124.[Medline]
Castillo, P.E., Carleton, A., Vincent, J.D. and Lledo, P.M. (1999) Multiple and opposing roles of cholinergic transmission in the main olfactory bulb. J. Neurosci., 19, 91809191.
Chalansonnet, M. and Chaput, M.A. (1998) Olfactory bulb output cell temporal response patterns to increasing odor concentrations in freely breathing rats. Chem. Senses, 23, 19.[Abstract]
Chen, W.R., Shen, G.Y., Shepherd, G.M., Hines, M.L. and Midtgaard, J. (2002) Multiple modes of action potential initiation and propagation in mitral cell primary dendrite. J. Neurophysiol., 88, 27552764.
Christie, J.M., Bark, C., Hormuzdi, S.G., Helbig, I., Monyer, H. and Westbrook, G.L. (2005) Connexin36 mediates spike synchrony in olfactory bulb glomeruli. Neuron, 46, 761772.[CrossRef][ISI][Medline]
Claverol, E.T., Brown, A.D. and Chad, J.E. (2002) A large-scale simulation of the piriform cortex by a cell automaton-based network model. IEEE Trans. Biomed. Eng., 49, 921935.[CrossRef][ISI][Medline]
Cleland, T.A. and Linster, C. (1999) Concentration tuning mediated by spare receptor capacity in olfactory sensory neurons: a theoretical study. Neural. Comput., 11, 16731690.
Cleland, T.A. and Linster, C. (2002) How synchronization properties among second-order sensory neurons can mediate stimulus salience. Behav. Neurosci., 116, 212221.[CrossRef][ISI][Medline]
Cleland, T.A. and Linster, C. (2003) Central olfactory processing. In Doty, R.L. (ed.), Handbook of Olfaction and Gustation, 2nd ed. Marcel Dekker, New York, pp. 165180.
Cleland, T.A., Morse, A., Yue, E.L. and Linster, C. (2002) Behavioral models of odor similarity. Behav. Neurosci., 116, 222231.[CrossRef][ISI][Medline]
Cleland, T.A. and Narla, V.A. (2003) Intensity modulation of olfactory acuity. Behav. Neurosci., 117, 14341440.[CrossRef][ISI][Medline]
Cleland, T.A. and Sethupathy, P. (2004) Non-topographical contrast enhancement disambiguates high-dimensional odor representations. Soc. Neurosci. Abstr., 531.2.
Davison, A.P., Feng, J. and Brown, D. (2000) A reduced compartmental model of the mitral cell for use in network models of the olfactory bulb. Brain Res. Bull., 51, 393399.[CrossRef][ISI][Medline]
Davison, A.P., Feng, J. and Brown, D. (2003) Dendrodendritic inhibition and simulated odor responses in a detailed olfactory bulb network model. J. Neurophysiol., 90, 19211935.
Debarbieux, F., Audinat, E. and Charpak, S. (2003) Action potential propagation in dendrites of rat mitral cells in vivo. J. Neurosci., 23, 55535560.
Di Prisco, G.V. and Freeman, W.J. (1985) Odor-related bulbar EEG spatial pattern analysis during appetitive conditioning in rabbits. Behav. Neurosci., 99, 964978.[CrossRef][ISI][Medline]
Djurisic, M., Antic, S., Chen, W.R. and Zecevic, D. (2004) Voltage imaging from dendrites of mitral cells: EPSP attenuation and spike trigger zones. J. Neurosci., 24, 67036714.
Duchamp-Viret, P., Duchamp, A. and Chaput, M.A. (2000) Peripheral odor coding in the rat and frog: quality and intensity specification. J. Neurosci., 20, 23832390.
Duchamp-Viret, P., Duchamp, A. and Chaput, M.A. (2003) Single olfactory sensory neurons simultaneously integrate the components of an odour mixture. Eur. J. Neurosci., 18, 26902696.[CrossRef][ISI][Medline]
Duchamp-Viret, P., Duchamp, A. and Vigouroux, M. (1989) Amplifying role of convergence in olfactory system a comparative study of receptor cell and second-order neuron sensitivities. J. Neurophysiol., 61, 10851094.
Duchamp-Viret, P., Duchamp, A. and Vigouroux, M. (1990) Temporal aspects of information processing in the first two stages of the frog olfactory system: influence of stimulus intensity. Chem. Senses, 15, 349365.
Eeckman, F.H. and Freeman, W.J. (1990) Correlations between unit firing and EEG in the rat olfactory system. Brain. Res., 528, 238244.[CrossRef][ISI][Medline]
Ennis, M., Linster, C., Aroniadou-Anderjaska, V., Ciombor, K. and Shipley, M.T. (1998) Glutamate and synaptic plasticity at mammalian primary olfactory synapses. Ann. N. Y. Acad. Sci., 855, 457466.
Erdi, P., Grobler, T., Barna, G. and Kaski, K. (1993) Dynamics of the olfactory bulb: bifurcations, learning, and memory. Biol. Cybern., 69, 5766.[CrossRef][ISI][Medline]
Ermentrout, G.B. and Kleinfeld, D. (2001) Traveling electrical waves in cortex: insights from phase dynamics and speculation on a computational role. Neuron, 29, 3344.[CrossRef][ISI][Medline]



