- Research article
- Open Access
Non-topographical contrast enhancement in the olfactory bulb
BMC Neurosciencevolume 7, Article number: 7 (2006)
Contrast enhancement within primary stimulus representations is a common feature of sensory systems that regulates the discrimination of similar stimuli. Whereas most sensory stimulus features can be mapped onto one or two dimensions of quality or location (e.g., frequency or retinotopy), the analogous similarities among odor stimuli are distributed high-dimensionally, necessarily yielding a chemotopically fragmented map upon the surface of the olfactory bulb. While olfactory contrast enhancement has been attributed to decremental lateral inhibitory processes among olfactory bulb projection neurons modeled after those in the retina, the two-dimensional topology of this mechanism is intrinsically incapable of mediating effective contrast enhancement on such fragmented maps. Consequently, current theories are unable to explain the existence of olfactory contrast enhancement.
We describe a novel neural circuit mechanism, non-topographical contrast enhancement (NTCE), which enables contrast enhancement among high-dimensional odor representations exhibiting unpredictable patterns of similarity. The NTCE algorithm relies solely on local intraglomerular computations and broad feedback inhibition, and is consistent with known properties of the olfactory bulb input layer. Unlike mechanisms based upon lateral projections, NTCE does not require a built-in foreknowledge of the similarities in molecular receptive ranges expressed by different olfactory bulb glomeruli, and is independent of the physical location of glomeruli within the olfactory bulb.
Non-topographical contrast enhancement demonstrates how intrinsically high-dimensional sensory data can be represented and processed within a physically two-dimensional neural cortex while retaining the capacity to represent stimulus similarity. In a biophysically constrained computational model of the olfactory bulb, NTCE successfully mediates contrast enhancement among odorant representations in the natural, high-dimensional similarity space defined by the olfactory receptor complement and underlies the concentration-independence of odor quality representations.
Primary olfactory sensory neurons (OSNs) line the nasal epithelium and respond to the presence of odors that diffuse through the nasal mucus layer and bind to olfactory receptors expressed on OSN cilia. Each OSN expresses only one or a few species of olfactory receptor, which define the molecular receptive range , or chemical receptive field, of that OSN. The axons of several thousand OSNs expressing the same receptors converge onto a few discrete glomeruli in the olfactory bulb, within which they form glutamatergic synapses with the dendrites of mitral cells, periglomerular cells, and external tufted cells. Both mitral and external tufted cells also form excitatory synapses onto periglomerular cells within the same glomerulus, whereas the periglomerular cells inhibit local mitral cells as well as presynaptically inhibiting OSN output synapses [2, 3] (Figure 1). The mitral cells (along with middle and deep tufted cells) are the primary output neurons of the olfactory bulb, projecting via axon collaterals to multiple central structures within the brain . A detailed description of these and other synaptic relationships within the olfactory bulb has been provided by Shipley et al. .
Due to the convergence of OSNs expressing the same olfactory receptor proteins, the pattern of activated glomeruli on the surface of the olfactory bulb reflects the pattern of activated olfactory receptors, identifying the constellation of chemical qualities that together constitute the presented odor. Many studies have attempted to elucidate the organization of this chemosensory map on the surface of the olfactory bulb. However, while great efforts have been made to infer meaningful patterns in these data, essentially all such studies have revealed fragmented, patchy maps without a clear, predictive ordering of the characteristics of chemical stimuli. This result stands in sharp contrast to analogous sensory maps in, for example, somatosensory, visual, and auditory cortices, in which neurons with overlapping receptive fields are located in predictable locations adjacent to one another . Fragmented stimulus quality maps such as those observed in the olfactory bulb pose unique problems for stimulus processing.
Contrast enhancement is a common property of sensory systems that serves to narrow, or sharpen, sensory representations by specifically inhibiting neurons on the periphery of the representation – e.g., the edges of a retinal image – in order to enhance its figure-background contrast. In many sensory systems, contrast enhancement is mediated by lateral inhibitory projections. Crucially, the effectiveness of lateral inhibition in this context depends upon the topographical mapping of stimulus similarities within the relevant brain region, such that the physical proximity of neurons reliably reflects the similarity of the information that they mediate. In the retina, for example, the spatial contrast of visual images is enhanced by lateral inhibitory projections within the two dimensions of the retinal field . Overlapping regions of the visual field are sampled by adjacent photoreceptors; hence, physically neighboring neurons mediate similar sensory information. Consequently, the projection of inhibition by neurons onto their physical neighbors is an effective means of projecting inhibition onto those neurons that mediate similar sensory information. In the auditory system, frequency tuning in the inferior colliculus and medial geniculate body is similarly sharpened along the single dimension of frequency [8, 9]. One-dimensional frequency tuning maps in auditory structures ensure that physically neighboring neurons will encode similar sound frequencies; hence, nearest-neighbor lateral inhibition again can effect mutual inhibition among similarly-tuned neurons. Critically, this mechanism of contrast enhancement is effective in these two modalities only because they are both low-dimensional. As neural cortices are layered structures, and thus functionally two-dimensional, contrast enhancement via nearest-neighbor lateral inhibition is effective only for modalities in which the relevant similarities among stimulus elements can be mapped continuously onto two or fewer dimensions.
While odor quality representations within the olfactory bulb are clearly sharpened by contrast enhancement , the fragmentation of odor maps indicates that the physical proximity of glomeruli cannot reliably connote similarity in their molecular receptive ranges, rendering nearest-neighbor lateral inhibition ineffective at sharpening odor stimulus representations (see below). Indeed, quantitative analyses of calcium imaging data from the analogous insect antennal lobe have demonstrated that glomerular tuning similarities do not correlate with their physical proximity, and that a pattern of functional inhibition reflecting the former, not the latter, better reproduces the observed input-output function of antennal lobe glomeruli . While such a proximity-independent functional topology is theoretically achievable using targeted lateral projections (i.e., selective inhibitory projections to similarly-tuned glomeruli irrespective of physical proximity), the combinatorial requirements for the formation of such a projection pattern are intractable [12, 13]. Furthermore, the probability of receptive field overlap between odor receptors is not only a function of these receptors' structures, but also of the statistics of the chemical environment. As the distribution of odor stimuli in natural environments is constantly in flux, quantitative measures of tuning similarity among glomeruli (and hence the appropriate target glomeruli for lateral inhibition) will not remain stable over time. Hence, any contrast enhancement mechanisms dependent upon specific morphological projections will be poorly optimized for most odor environments. An alternative neural circuit algorithm is required, one that can both perform similarity-dependent computations on fragmented representations and continually remain optimized with respect to changing chemical environments.
We here propose and demonstrate a novel neural mechanism, non-topographical contrast enhancement (NTCE), which permits the regulation of olfactory tuning curves in an arbitrarily high-dimensional sensory space, using established olfactory bulb circuitry, without reliance on specific connections among glomeruli. This lack of dependence on the anatomical proximity of glomeruli enables the bulb to remain optimized with respect to the statistics of changing odor environments, and additionally renders the system robust to the generation of novel glomeruli over the course of development, evolution, or experimental manipulation [14–18]. The NTCE algorithm is consistent with, and in some cases explains, diverse physiological datasets from the olfactory bulb, and in particular replicates the canonical contrast enhancement results of Yokoi and colleagues . We first present the core principles of NTCE in an illustrative model, and then demonstrate that it can be implemented by established glomerular layer circuitry using a biophysically constrained compartmental model of the olfactory bulb glomerular layer.
Results and discussion
Fragmented maps and the dimensionality of stimulus representations
All auditory stimuli can be mapped as spectra along a single axis of frequency. Retinotopic images of visual stimuli can be unambiguously represented as two-dimensional patterns of ganglion cell activation. However, substantially higher-dimensional spaces are required to map the analogous sensory space in which odor stimuli are distributed [11, 19–21]. Heuristically, this reflects the large number of different ways that one can gradually alter the molecular structure and charge distributions of odorous molecules or the component ratios of multicomponent odors; more concretely, it reflects the large number of independent sensors (i.e., types of odorant receptor) that provide the raw material for olfactory representations. Each glomerulus represents a convergent population of olfactory sensory neurons expressing the same odorant receptor protein, which determines its molecular receptive range; due to this convergence, glomeruli can be treated as selective, low-noise chemosensory units . In mice, there are roughly 1000 such receptor types and hence roughly 1000 chemoreceptively distinct sets of glomeruli [23–25]. Odorant stimulation activates characteristic, odor-specific groups of glomeruli, and perceptual similarities in odor quality correlate with the similarities in these odor-evoked patterns of glomerular activation [26, 27]. As each glomerulus is capable of being independently activated to differing degrees, whether by natural or artificial stimuli , there are roughly 1000 potential dimensions in mice along which odor quality representations may be gradually varied. Consequently, in order for nearest-neighbor relationships to reliably reflect odorant similarity in mice, a topographical map of odor representations would need to occupy an approximately 1000-dimensional similarity space. In layered neural structures, however, only two physical dimensions are available for topographical mapping.
As the olfactory sensory neuron complement converges onto the glomerular layer of the olfactory bulb, these high-dimensional representations of odorant similarity are projected onto the two physical dimensions of the olfactory bulb input layer, a process facilitated by the glomerular architecture of that layer. Mathematically, whenever high-dimensional representational patterns are projected onto lower-dimensional spaces, they form discrete maps, that is, patchy maps with embedded discontinuities (Figure 2). Neural network algorithms such as the self-organizing map [29, 30], which perform such transformations adaptively, are able to retain local similarity relationships up to certain limits imposed by the topology of projection, at which point discontinuities are inevitable. One familiar example of this phenomenon in the nervous system is the Penfield-Rasmussen somatosensory homunculus [31, 32], essentially a one-dimensional representation of the two-dimensional external body surface mapped coronally along the posterior parietal cortex. While body surface contiguity is maintained within discrete segments of this map, there are also sharp discontinuities mandated by its projection onto a one-dimensional surface as well as the presence of anisotropic elaborations (such as limbs) that disrupt contiguity relationships. For example, thumb and neck are represented in adjacent cortical regions, as are toes and genitalia, along with many more subtle examples, while representations of the arms and hands are interposed between those of the face and the rest of the head. The existence of such discontinuities in the neural representation is inevitable, as illustrated in Figure 2; furthermore, the frequency of the discontinuities will increase along with the difference between the intrinsic dimensionality of the representation and that of the neural substrate. Consequently, as the high-dimensional odor space defined by an animal's olfactory receptor complement and current chemical environment projects onto the two physical dimensions of the olfactory bulb glomerular layer, the chemotopic similarity map expressed thereon will be unavoidably fragmented, incorporating many embedded discontinuities. This is exactly what has been observed in numerous studies of odor representations across the glomerular layer of the olfactory bulb, although these commonalities have not always been emphasized by the different laboratories authoring these studies.
Substantial efforts have been made by several laboratories to elucidate an underlying order to the distribution of glomeruli across the surface of the olfactory bulb. While the multiple glomeruli activated by a given odorant stimulus are typically scattered, imaging studies of odorant-specific glomerular activation patterns have observed a broad tendency for the centroids of these activated groups to migrate to progressively more ventral regions as the odorant molecules become larger – e.g., straight-chain aldehydes of increasing aliphatic chain lengths – and more hydrophobic [33–35]. This phenomenon has been attributed to the differential sorption of odorants along the nasal inspiratory path (reviewed in ). Similarly, within the dorsal region of the bulb (the aspect most accessible to in vivo imaging studies), longer chain lengths have been associated with centroids of activation positioned more rostrally and laterally [34, 37–39], though concentration artifacts may have influenced this finding . Further studies have established ad hoc modules or clusters within the olfactory bulb that exhibit greater densities of glomeruli sensitive to particular odorant moieties than would be expected by chance [40–44]. However, these broad chemotopic profiles across the bulbar surface must not be confused with a reliable topography of feature similarity across the olfactory bulb. The chemotopic maps that have been described do not extend to the finer scale of individual glomeruli upon which contrast enhancement mechanisms operate [37, 45]. In each of the imaging studies cited above, many glomeruli responding to similar odorant features are located at considerable distances from one another, and there are innumerable examples of strongly activated glomeruli located adjacent to unresponsive glomeruli in response to any given odorant – a pattern indicative of a highly fragmented map of odor quality. While both of these phenomena have been consistently observed in bulbar imaging studies, they have rarely attracted substantial discussion. The important consequence of these data is that, while some neighboring glomeruli are similarly tuned for odorant stimuli, others are not; hence, the physical proximity of glomeruli is not a reliable indicator of similarity in their molecular receptive ranges. For center-surround lateral inhibition to reliably mediate contrast enhancement, each of the odorant stimuli that substantially activate any given glomerulus must also activate all of its physical neighbors, albeit to a greater or lesser degree. This is not the case in the olfactory bulb.
Glomerular circuitry implements non-topographical contrast enhancement (NTCE)
Contrast enhancement can be functionally defined as a process of competition between neurons proportional to the similarity of the information that they mediate; this definition is, of course, agnostic to mechanism. The olfactory bulb input layer provides precisely the architecture necessary to implement a form of non-topographical contrast enhancement that is not limited by its anatomically two-dimensional architecture. The NTCE algorithm described here depends upon two interacting computational mechanisms derived from established bulbar circuitry. The first mechanism operates independently within each glomerulus at the synaptic triad between OSN axons, periglomerular cell spines, and mitral cell apical dendrites [46–48]. Direct sensory input from OSN axonal arborizations activates both periglomerular neurons and mitral cell primary dendrites, while the periglomerular cells in turn inhibit local mitral cell primary dendrites via dendrodendritic synapses (Figure 1). This architecture creates a contrast enhancement generator element within each glomerulus that transforms mitral cell activity along a half-hat function (Equation 2, Figure 3; see below) such that mitral cells connected to moderately activated glomeruli – i.e., the "edges" of the high-dimensional odor representation – are specifically inhibited out of the active ensemble. When all glomeruli are combined, their collective output activity reflects the input pattern of glomerular activation filtered through a Mexican hat function with a maximum dimensionality equal to the number of glomeruli and a topology of similarity inherited from the odor environment. In other words, NTCE generates a pattern of inhibition equivalent to that which would be generated by a competitive lateral inhibitory network in which the strength of inhibition between any two glomeruli was proportional to the similarity in their molecular receptive ranges. This mechanism is illustrated below in an abstract model and in a uniglomerular compartmental model.
The second mechanism of NTCE is a global negative feedback loop utilizing the lateral excitatory network of external tufted and short-axon cells (the ET/SA network; Figure 1). This network feeds the average level of bulbar activity back upon all mitral cells as inhibition, which is necessary in order to moderate the influence of odor concentration and maintain the first NTCE mechanism within its effective dynamic range. This mechanism is illustrated below in a multiglomerular compartmental model.
Computational principles: illustrative model
Irrespective of modality or dimensionality, successful contrast enhancement implies that the single or few most strongly activated units will yield robust output in response to a given stimulus, while output activity in more modestly activated units will be specifically inhibited and minimally-activated units will remain inactive. That is, output activity as a function of input activity yields a half-hat function (Figure 3, Mi out ; see Methods), the generator of the familiar Mexican hat function (Figure 3, inset) and a signature function of contrast enhancement. As shown in a simple illustrative model (see Methods), this half-hat function can be generated independently within each glomerulus by driving parallel sigmoidal excitatory and inhibitory processes with the same sensory input, given that the inhibitory process both is more sensitive to that input than the parallel excitatory process and saturates at a lower activity level. Both of these conditions are favored by the much smaller volume and higher input resistance of periglomerular cell spines compared with mitral cell dendrites. Under these conditions, each glomerulus inhibits its own mitral cell-mediated output in scaled proportion to its sensory input, such that its net mitral output level exhibits a half-hat function along the axis of ligand-receptor affinity (Figure 3, Mi out ). Hence, mitral cells innervating glomeruli that are the most strongly tuned to a given odorant feature will be activated, while mitral cells innervating moderately well-tuned glomeruli will be inhibited. Because of this natural attunement to relative degrees of glomerular activation irrespective of their location or interconnectivity, NTCE naturally inherits the intrinsic topology of the external chemosensory environment (as filtered by the animal's complement of primary odorant receptors). This is a critical property for a system likely to encounter environments in which the distribution of relevant odorants is unpredictable.
Uniglomerular compartmental model
In order to determine whether the established electrophysiological and cytoarchitectonic properties of olfactory glomerular circuitry could indeed mediate these computations, we constructed a compartmental model of olfactory glomerular circuitry using the simulation language NEURON [49–51]. We first implemented a uniglomerular model in order to demonstrate the first of the two mechanisms underlying NTCE. The model included representations of mitral cells [52, 53], periglomerular cells constrained to available anatomical and biophysical data and exhibiting appropriate single-cell physiological responses to inputs [2, 47, 54–58], and olfactory sensory neuron terminals designed to convey the average activation level of a coactivated and highly convergent olfactory sensory neuron population. To better illustrate the NTCE phenomenon, we limited the number of unconstrained parameters in our model and entirely omitted post-glomerular circuitry such as granule cells and mitral-granule reciprocal connectivity. Accordingly, we made no effort to model mitral cell bistability or complex spike patterning; the temporal window of interest in this model is the period preceding the first mitral cell spike, determining whether and when a mitral cell will first fire and hence preceding the recruitment of granule cells into the olfactory response. Mitral cell spike rates are used herein solely as a convenient measure of the degree of neural activation.
First, identical levels of olfactory sensory input were delivered to periglomerular dendritic spines and mitral distal dendritic tufts via glutamatergic synapses in a uniglomerular version of the model. We did not alter the relative sensitivities of the two neuron types to sensory input other than that afforded by the periglomerular spine's greater input resistance and smaller volume, though such a sensitivity difference could have been readily generated by manipulating postsynaptic receptor densities. The periglomerular cell spine also inhibited the mitral cell dendrite via a thresholded, graded mechanism. Stimulation of the OSN terminals of this network with low-affinity odorant stimuli generated no activity in mitral cells (Figure 4Aviii), while the mitral cell activation evoked by medium-affinity odorants was dissipated by periglomerular inhibition before spike initiation (Figure 4Avi–vii). Stimulation with the highest-affinity odorants, however, evoked activity in mitral cells that overcame PG-mediated inhibition and produced action potentials (Figure 4Av). Plotting mitral cell activation as a function of odor ligand-receptor affinity generated a half-hat function (Figure 4B) corresponding to that in the illustrative model (Figure 3, Mi out ), demonstrating a successful, classical contrast enhancement generator element within a single glomerulus. In contrast, mitral cells exhibited monotonically increasing activation in response to presentation of higher-affinity odorants when PG-mediated inhibition was absent (Figure 4Aiv–i).
In our analysis, the established biophysical characteristics of glomerular neurons and synapses materially contribute to the effectiveness of NTCE. The synaptic triad between OSNs, periglomerular cell spines, and mitral cells [46–48] is optimal for enabling periglomerular inhibition to effectively shunt the synaptic excitation of mitral cells, owing to the close proximity of excitatory and inhibitory inputs [59–61]. The high input resistance, thin dendritic processes, and tiny spines of PG neurons effect a a a rapid local depolarization of membrane upon excitation by sensory input (Figure 4Aix–x). This enables the fast release of GABA onto mitral cell primary dendrites before the concomitant excitation of these larger, leakier dendrites can evoke an action potential in the electrotonically distant soma or somewhat nearer primary dendrite (Figure 4Aiii–iv; [52, 53, 62]. Indeed, recordings from mitral cells in vivo have shown fast inhibitory postsynaptic potentials (ipsps) preceding mitral cell spike initiation even in the most rapidly-responding mitral cells [58, 63, 64]. While this chloride-mediated shunt inhibition may often be insufficiently hyperpolarizing to evoke clearly visible ipsps, it is highly effective at blocking spike initiation by multiplicatively decrementing the cumulative effect of excitatory sensory input (Figure 4Avi–vii;[61, 65]; see also ). Furthermore, the activation of low-threshold T-type calcium currents in PG cells , while not essential to contrast enhancement, enhances its efficacy by potentiating the effect of small inputs on local PG inhibitory output (Figure 4Aix). Finally, preintegrative mechanisms such as NTCE, in which inhibition is directed to mitral cell inputs prior to their integration into a spiking response, enjoy theoretical computational advantages that postintegrative lateral inhibitory mechanisms (such as those mediated by mitral-granule interactions) cannot provide, including an improved capacity to represent arbitrary combinations of overlapping inputs and to recognize partial input patterns [67, 68].
Multiglomerular compartmental model
The second mechanism required by the NTCE algorithm is essentially a global inhibitory feedback loop. This mechanism, based on the ET/SA lateral excitatory network and inhibitory PG cells (Figure 1), mitigates the effects of absolute odor concentration, thus enabling mitral cell activity patterns to better reflect relative odor ligand-receptor affinities and efficacies and hence preserve the representation of odor quality across a reasonable concentration range. We constructed a multiglomerular version of this model in order to directly illustrate non-topographical contrast enhancement among multiple glomeruli differentially activated by the same odorant, and to assess the ability of the olfactory bulb ET/SA network to normalize bulbar responses to different odor concentrations via global feedback inhibition (see Methods). In the present model, ET/SA network activity was simulated by projecting slow synaptic inhibition from each of the other glomeruli onto mitral cell primary dendrites (Figure 4, inset), effectively translating the NTCE half-hat function of each glomerulus on its abscissa in proportion to average bulbar input intensity such that only a small and consistent population of the most strongly activated mitral cells overcame this activity-dependent inhibition. Ten glomeruli were interconnected via the ET/SA network such that all mitral cells received identically-weighted inhibitory inputs from each of the other nine glomeruli. The ten glomeruli were assigned different affinities for each test odorant such that each odorant representation constituted a unique ten-dimensional vector. Addition of intraglomerular inhibition mediated by periglomerular cells and global normalizing inhibition mediated by the ET/SA network strongly disambiguated the mitral cell response profiles evoked by two similar odorants (Figure 5A,B).
Mechanisms of concentration compensation
The persistence of odor quality across concentrations is a difficult and continuing problem in olfaction. The most plausible explanation for this property is that common elements are retained among the representations of different concentrations of the same odorant, suggesting that some sort of normalization for concentration is likely to occur within the olfactory bulb. Indeed, such a normalization process has been proposed for the olfactory bulb based on imaging data showing that normalized odor-specific glomerular activity maps are relatively similar across stimulus concentrations [69, 70], and that the degree of similarity in these normalized maps predicts similarities in perceived odor quality across concentrations more reliably than do non-normalized maps . Consistent with these proposals, the NTCE circuit mechanism inherently compensates for concentration as the primary glomerular representation, measured predominantly from OSN presynaptic arbors, is transformed into the secondary representation based on patterns of activated mitral cells. Specifically, global feedback inhibition mediated via the ET/SA network serves to improve stability in odor quality representations among mitral cells across reasonable concentration ranges. When a given model odorant was applied at a series of concentrations (with synaptic weights unchanged), the higher concentrations recruited increasingly broad distributions of (presynaptically) activated glomeruli, as has been observed in imaging studies [38, 69, 71]. With intraglomerular inhibition intact, but in the absence of interglomerular inhibition, mitral cell activation patterns also broadened with increased odorant concentrations, as these higher concentrations recruited more weakly-tuned glomeruli into the ensemble of excited mitral cells (Figure 5C). However, with an intact interglomerular ET/SA network, each odorant concentration in the series could evoke comparable levels of mitral cell activity in a consistent, odor-specific subpopulation of mitral cells, conferring a degree of concentration-independence upon odor-specific mitral cell response profiles (Figure 5D). While the synaptic strengths underlying interglomerular inhibition could be tuned so as to produce either gradual recruitment or gradual elimination of mitral cells from the ensemble as odor concentration increased, a relative constancy in the ensemble of activated mitral cells is conservative in that odor-evoked mitral cell spike counts do not vary consistently with changing odorant concentrations [72–74], and certainly do not scale with the monotonically increasing presynaptic activation of olfactory glomeruli as shown in imaging studies. More specifically, given realistic odors that present multiple, diverse ligands for olfactory receptors, this constancy is only an average; individual mitral cells' responses would be expected to change qualitatively as their relative degrees of input activation change with respect to the population average. Indeed, mitral cell response profiles measured across concentration series have presented appropriately complex results; for example, nonresponsive mitral cells often become inhibited as odor concentration increases [72–74], while other mitral cells change from a net inhibitory response to a net excitatory response .
The normalization of olfactory representations for concentration does not imply that concentration information is lost. Model mitral cells excited by an odor stimulus exhibited shorter spike latencies as concentration increased (Figure 4), consistent with experimental data [63, 75]. Concentration may also be represented in the degree of stimulus-evoked spike synchronization among mitral cells, as has been suggested for the insect antennal lobe . Finally, middle and deep tufted cells, as well as the trigeminal system, may contribute to the representation and perception of odor intensity . NTCE does not eliminate concentration information, but simply mitigates its impact on odor quality representations.
Contrast enhancement tested with systematically varying odor stimuli
The clearest demonstration of contrast enhancement in the olfactory bulb has been provided by Yokoi and colleagues . Briefly, these authors showed that an individual mitral cell could respond to systematically varied odorant stimuli with a classical Mexican hat function: e.g., one featured mitral cell exhibited no response to a 3-carbon n-aldehyde, inhibition to a 4-carbon n-aldehyde, excitation to 5-, 6-, and 7-carbon n-aldehydes, inhibition to an 8-carbon n-aldehyde, and no response to 9-, 10-, or 11-carbon n-aldehydes. Furthermore, administration of the GABAA antagonist bicuculline appeared to broaden this Mexican hat response, transforming an inhibitory mitral cell odor response into an excitatory response. While the authors attributed this effect to the blockade of granule cell-mediated inhibition, it is also consistent with a blockade of GABAergic periglomerular cell-mediated inhibition and with NTCE. We modeled these authors' data by presenting a systematically varying odorant series to our ten-glomerulus model (see Methods). As NTCE is independent of specific lateral projections, the model did not have to be adjusted to perform contrast enhancement along the particular axes of variation used to model the homologous series of odorants. NTCE was able to replicate the contrast enhancement function within model mitral cells based on periglomerular inhibition alone (Figure 6A; compare with Figure 2A in reference ). Furthermore, blockade of this periglomerular inhibition transformed an inhibitory mitral cell odor response into an excitatory response (Figure 6B; compare with Figure 5D in reference ).
Interpretations of existing data
Contrast enhancement among the representations of structurally similar odorants has been widely attributed to a process of reciprocal inhibition between mitral cells mediated via mutual dendrodendritic synapses with the dendrites of inhibitory granule cells [10, 78–81], though lateral interactions via periglomerular neurons have also been proposed to mediate this function [2, 82, 83]. Lateral projections and their putative formation of an inhibitory surround have been both explicitly and implicitly presumed to reflect the two physical dimensions afforded by the columnar structure of the olfactory bulb such that physically neighboring mitral cells mutually inhibit one another more strongly than do more distantly neighboring cells [2, 10, 38, 79, 80, 84, 85], perhaps by analogy with the retina . However, as explained above, physical proximity-based solutions such as decremental lateral inhibition cannot effectively mediate similarity-dependent computations such as contrast enhancement in a high-dimensional modality such as olfaction.
This problem has been specifically investigated in the analogous honeybee antennal lobe, using calcium imaging of odor-evoked glomerular activity among convergent OSNs and projection neurons (PNs; analogous to vertebrate mitral cells). Increasing odorant concentrations increased glomerular activation levels and recruited additional glomeruli, as in the vertebrate OB, while PNs tended to become more inhibited as odorant concentrations increased, particularly within reasonable concentration ranges . Furthermore, in a systematic study of the transformation between OSN and PN odor-evoked activity patterns, PN activity patterns were significantly better fit by a model using a functionally-based pattern of inhibition (comparable to that generated by NTCE) than by a model employing nearest-neighbor lateral inhibition . Indeed, the best fits based on nearest-neighbor inhibitory projections were similar to control fits based on random interglomerular projections.
Even within the vertebrate olfactory bulb, the hypothesis of proximity-based representations of odor similarity, while popular, is contraindicated in the literature. Mitral cells innervating neighboring glomeruli exhibit unrelated odorant response profiles , in contrast to mitral cells innervating the same glomerulus, which have very similar profiles . Furthermore, as odor quality is gradually changed, using common, behaviorally validated models for the systematic variation of odor quality [10, 27, 90, 91], the glomeruli recruited and lost from the gradually shifting odor representation are broadly distributed across the bulb rather than being consistently adjacent [42, 45, 71]. Quantitative measurements of response-correlated maps have revealed broad regions of odor-evoked inhibition among mitral cells that do not directly sample from strongly activated glomerular domains , a result consistent with numerous studies suggesting that a net inhibition is a common response of mitral cells to odor or olfactory nerve stimulation [54, 63, 72, 74, 75, 89, 92]. While this result alone does not exclude the possibility of (broad) center-surround inhibition, it is more consistent with a widespread, locally-induced, nontopographical pattern of inhibition such as that produced by NTCE. Indeed, in support of the latter hypothesis, strong odor-evoked inhibition of mitral cells has been clearly observed at locations far removed from the region(s) of peak activation evoked by that odor , and in vivo recordings coupled with calcium imaging data have suggested that action potential amplitudes and evoked calcium transients in mitral cell lateral dendrites do not decline in amplitude with increasing distance from the soma, but instead are propagated nondecrementally, suggesting an axon-like mechanism . Both of these results argue against a decremental, topographically localized projection of lateral inhibition. Finally, while granule cell-mediated lateral inhibition can affect spike timing in mitral cells , these inhibitory synaptic inputs, widely distributed along mitral cell lateral dendrites, may not be capable of preventing spike initiation in mitral cells, particularly given that mitral cell spikes may be initiated in the primary dendrite and even potentially the glomerular tuft [52, 62, 95]. This implicates glomerular circuitry, rather than external plexiform layer circuitry, as the primary determinant of whether a given mitral cell will respond to a given odor input with action potentials or with inhibition. However, the long projections of mitral cell lateral dendrites and their dense interconnectivity with granule cell dendrites are likely to underlie a network of mutual inhibition among mitral cells approaching all-to-all connectivity, and hence are likely to be capable of performing secondary computations upon arbitrarily high-dimensional maps. The mitral-granule network is therefore potentially well suited to transform glomerular output in accordance with odor learning, presumably through the learned modification of mitral-granule reciprocal synapses and learning-correlated granule cell neurogenesis [96–99], and perhaps involving the generation of field oscillations and the synchronization of mitral cell spikes [94, 100–103]. Indeed, a contrast enhancement transformation based on local field oscillations and spike synchronization has been modeled in antennal lobe PNs, analogous to mitral cells [76, 104].
In the interests of simplicity, some known glomerular circuit elements were omitted from the present model. First, periglomerular cells deliver GABAB-ergic and D2 dopaminergic presynaptic inhibition onto the axon terminals of convergent OSNs [105–109]. This feedback inhibition regulates OSN glutamate release and has been proposed to help normalize OSN synaptic output with respect to concentration; such an effect would be expected to facilitate NTCE. Second, mutual excitation among mitral cells innervating the same glomerulus [110–113] would be expected to underlie cooperativity in the glomerular output response, sharpening the distinction between highly activated glomeruli in which mitral cells generate spikes and more modestly-activated glomeruli in which incoming OSN activity is dissipated without evoking mitral spikes. Third, self- (and potentially mutual) inhibition among PG cells, while hindering PG cell spike initiation, actually depolarizes PG cells and hence may have a potentiating effect on their graded inhibition of mitral cell apical dendrites . This putative positive feedback loop may further contribute to the efficacy of PG inhibition of mitral cell spike initiation. Finally, some PG cells have axons that project to nearby glomeruli; while sparse, these projections have served as a basis for some models of lateral inhibition [82, 83]. The utility of these projections is not clear; they potentially may serve to improve the functional differentiation of neighboring glomeruli with very similar molecular receptive ranges, such as might be generated via activity-dependent segregation .
The compartmentalized architecture of olfactory glomeruli may reflect a specific adaptation for the representation and processing of high-dimensional computational spaces within the physical constraints of two-dimensional cortical layers. If true, this would explain the observation that a functionally similar glomerular architecture has evolved independently several times in diverse phyla , and indeed is a nearly universal characteristic of complex olfactory systems. The limited diversity observed in these analogous olfactory structures (olfactory bulb, antennal lobe) may even suggest fundamental constraints on the effective processing of olfactory information above a certain level of complexity. This principle may also extend to other brain regions possessing glomerular architectures, notably the cerebellar cortex, implying that these circuits are also involved in the processing of high-dimensional representations.
The NTCE algorithm presented here is the basis of a novel theory of bulbar function, integrating diverse data sets gathered by several laboratories. Fundamentally, it is a winner-take-most algorithm utilizing local feed-forward inhibition and global feedback inhibition to generate the competitive interactions among glomeruli typically associated with lateral inhibitory projections. Most importantly, it solves the problem of how to represent and process intrinsically high-dimensional sensory data within a physically two-dimensional neural cortex, while offering an explanation of the patchy, discontinuous odor quality maps observed in the olfactory bulb. As a model of contrast enhancement, it replicates the results of Yokoi et al.  illustrating the canonical Mexican hat function within single mitral cells and the dependence of this function on bulbar GABAA receptors. NTCE innately distributes inhibition among mitral cells according to the similarities in their molecular receptive ranges, an intractable problem for mechanisms based on lateral inhibitory projections embedded in discontinuous feature maps. Unlike mechanisms based on decremental lateral inhibitory projections, NTCE is independent of the physical proximity of similarly-tuned glomeruli, an important advantage given the broadly distributed representations of odor stimuli observed in the olfactory bulb. It offers a mechanism for the normalization of odorant representations with respect to concentration, reflecting the observation that mitral cell spike counts do not monotonically increase in proportion to primary sensory neuron spike counts as odorant concentrations are increased. Rather, it predicts complex effects on mitral cell activity patterns, including fast initial inhibition in response to odors and a variety of responses to odor concentration series owing in part to the unpredictable pattern of associations of multiple odor ligands with different odorant receptors. Finally, NTCE depends only on known properties of the olfactory bulb input layer, including its cellular morphologies and connectivity, passive membrane properties, active electrophysiological and pharmacological responses, and specialized glomerular architecture, to enable robust contrast enhancement in an objective, externally-defined, and high-dimensional similarity space.
Illustrative model of NTCE
NTCE depends upon two interacting computational mechanisms, each of which is consistent with available anatomical and physiological data. The first of these mechanisms, based on the rapid shunt inhibition of activated mitral cells by coactivated periglomerular spines, narrows the population of activated mitral cells with respect to the population of activated glomeruli such that only the most strongly activated subset of activated glomeruli generates action potentials in its associated mitral cells. Crucially, this contrast enhancement mechanism is independent of interglomerular projections, depending only upon relative degrees of glomerular activation and hence naturally inheriting the topology of odor similarity that is generated by the olfactory receptor complement. The second mechanism, essentially a global negative feedback loop, mitigates the effects of absolute odor concentration, thus enabling mitral cell activity patterns to reflect relative odor ligand-receptor affinities and efficacies and hence preserve the representation of odor quality across a reasonable concentration range.
The model depicted in Figure 3 was composed in Matlab (Mathworks, Natick, MA) and illustrates the first of these two mechanisms. For simplicity of illustration, odor concentrations are held constant. Mitral cells (Miin) and periglomerular cells (PG) within a given glomerulus independently receive identical levels of direct sensory input (depending on the odor ligand-receptor affinity of their OSNs; abscissa), although they respond to this input differently based on their passive membrane properties (see Results and Discussion). As neuronal activation levels in this illustration (ordinate) depend upon ligand-receptor binding between the odorant and the olfactory receptor species projecting to the glomerulus in question, the direct activations of both Mi and PG cells vary sigmoidally as a function of odor ligand-receptor affinity:
where K is a scaling constant determining relative excitability (due to input resistance), Y1/2 is the odor ligand-receptor affinity at which half-maximal output activity would be evoked in the mitral or PG cell via direct sensory input, A denotes the odor ligand-receptor affinity variable (abscissa), and m is the cooperativity (Hill equivalent; ) of each process. Net mitral cell output (Miout) – i.e., the generation of spiking activity – depends on both direct mitral cell activation by sensory input (Miin) and on the secondary inhibitory influence of periglomerular cell activity (PG), resulting in a net activation function for mitral cell output that depends on the difference between these two sigmoids:
in which subscript e denotes the parameters of mitral cell direct excitation (Miin) and subscript i denotes the inhibitory effect of periglomerular cell activity (PG). In the example shown in Figure 3, K e = 1, K i = 0.6, Y1/2(e)= 10-5 M, Y1/2(i)= 10-4 M, and both m terms are unity.
Contrast enhancement is typically described with an on-center/inhibitory surround, or Mexican hat, function ; Figure 3, inset). Traditional Mexican hat functions in one dimension are typically modeled as a difference of Gaussians, or as the difference between two differences of sigmoids. Differences between sigmoids of equal maxima (Ke, Ki) and cooperativities (me, mi) yield bell-like curves; the difference between two bell-like curves (whether gaussians or differences of sigmoids) generates a Mexican hat function given that the excitatory difference is of greater maximum magnitude and the inhibitory difference is of greater breadth. The corresponding function in NTCE is one half of a Mexican hat, modeled as a difference between sigmoids in which Ke > Ki and Y1/2(e) < Y1/2(i) and referred to as a half-hat function (Figure 3, Mi out ; equation 2). Mexican hat functions of arbitrary dimensionality can be generated by rotating this half-hat function around the coordinate of its maximum activation level in a given context. When Ke > Ki and Y1/2(e) < Y1/2(i) – both conditions which are favored by the smaller volume and higher input resistance of periglomerular cell spines compared with mitral cell dendrites  – the output process Miout (equation 2) exhibits contrast enhancement along the axis of odor ligand-receptor affinity A (see Results and Discussion).
Compartmental model of NTCE
The compartmental model depicted in Figures 4 and 5 was composed in the neural simulation language NEURON [49–51]. Three cell types were explicitly modeled: mitral cells, periglomerular cells, and olfactory sensory neuron (OSN) axonal arbors. To minimize unconstrained parameters, the ET/SA network  was modeled as a set of delayed inhibitory synapses from each OSN to all mitral cells sampling from glomeruli other than the one to which that OSN projected, hence abstracting the projection of interglomerular inhibition mediated through the chain composed of local external tufted cells, short-axon cells, and periglomerular cells in other glomeruli. Granule cells and other post-glomerular circuitry were omitted. For simplicity, AMPA and NMDA glutamatergic synapses were combined into a single mechanism. All synapses were thresholded and graded so as to reflect either genuinely graded synaptic communication in vivo (e.g., PG to Mi synapses) or the spike density functions of large neuronal populations (e.g., convergent OSN inputs to PG and Mi cells), and approached steady-state exponentially.
The mitral cell model was modified from that of Chen and colleagues [52, 53], which is based on a simple complement of Hodgkin-Huxley currents (INa, IK) with heterogeneous expression levels and optimized for its spike initiation properties. Model morphology and kinetic parameters were as in reference . Specifically, two apical dendritic tufts fed into a single primary dendrite, which connected with a soma, an axon hillock, and two secondary dendrites, all of which expressed INa at relatively low levels. The axon hillock in turn was connected to an initial segment and an axon with five node/internode pairs; the initial segment and nodes contained high levels of INa, facilitating spike initiation. For the present simulations, glutamatergic excitatory and GABAA-ergic inhibitory synaptic inputs were inserted at the distal ends of apical dendrites.
The periglomerular soma was connected to two narrow dendrites (20 um length, 1 um diam) from which protruded spine shafts and bodies (1 um diam). All periglomerular synaptic inputs and outputs were localized on spine bodies. While the axon contained only Hodgkin-Huxley currents, the dendrites and spines additionally contained a complement of currents capable of replicating low-threshold calcium bursts as observed in periglomerular cell recordings ; see also . Each current mechanism was constructed in the NMODL language , and in some cases adapted from existing mechanisms as follows: mammalian hippocampal Hodgkin-Huxley sodium and potassium currents , low-threshold T-type calcium current , periglomerular cell hyperpolarization-activated (H-type) cation current , A-type inactivating potassium current , and a calcium diffusion model . The kinetic parameters of each of these model mechanisms were retained, and have been made available on ModelDB  by their respective authors.
OSN axonal arbor
The OSN arbor was a passive, single-compartment model, serving as a source for graded glutamatergic excitation of mitral and periglomerular neurons; this graded excitation represented the spike density function of a large and highly convergent population of similarly tuned spiking neurons.
The uniglomerular compartmental model included only mitral and periglomerular cells and OSN axonal arbors, and is directly comparable to the illustrative model described above and in Figure 3. OSN inputs excited both periglomerular cell spines and mitral cell distal dendrites via graded glutamatergic synapses, while periglomerular cell spines inhibited mitral cell distal dendrites via graded GABAA-ergic synapses (EGlu = 0 mV; EGABA = -80 mV). All synaptic time constants were set to 1 ms to reflect multiple, imprecisely timed inputs. The uniglomerular network was stimulated with four different levels of olfactory sensory neuronal input (12.6, 20, 90, or 140 pA injected into the OSN terminal); for each level of sensory input, the evoked activity was measured in (a) mitral cell somata in the absence of a PG-to-Mi inhibitory synapse, (b) mitral cells in the presence of a PG-to-Mi inhibitory synapse, and (c) periglomerular cells (Figure 4).
Glomerular activation levels contribute to odor representations in that they reflect the pattern of relative affinities that the olfactory receptor complement has for odor ligands. However, ligands of lower affinity presented at higher concentrations can evoke levels of glomerular activation indiscriminable from those evoked by high-affinity ligands at lower concentrations . This enables changing odor concentrations to degrade the effectiveness of contrast enhancement and the integrity of the olfactory representation. This problem can be resolved by normalizing glomerular output with respect to the sum of activity across all glomeruli , thereby dynamically regulating the stringency of the bulb's selectivity for odorant features . In the olfactory bulb, this function is attributable to the ET/SA lateral excitatory network formed by external tufted and short-axon cells [2, 3]. Short-axon cells (the term is a misnomer) are indirectly activated by sensory input (Figure 1), ramify broadly and isotropically within the glomerular layer for distances of up to 850 um (interconnecting areas separated by up to 30 glomerular radii), and form excitatory synapses upon distant periglomerular neurons, which in turn inhibit local mitral cells . While the original study describing this ET/SA network proposed a center-surround topology , the observation that short-axon cells also form excitatory synapses upon one another and upon (excitatory) external tufted cells in multiple glomeruli, which in turn excite other SA cells, indicates the existence of a broad lateral excitatory network across the bulb. As short-axon cells are few in number and diverge broadly, sampling from 2–4 glomeruli and exhibiting extensive axonal branching throughout the glomerular layer, the ET/SA network is a wholly appropriate mediator of a broad and relatively indiscriminate projection of inhibition upon the mitral cell population. In the present model, ET/SA network activity fed back inhibition upon each mitral cell from all other glomeruli, hence translating the NTCE half-hat function on its abscissa in proportion to net bulbar input intensity such that only a small and consistent population of the most strongly activated mitral cells overcame this inhibition. In this way, NTCE could also establish a degree of concentration-independence in the secondary olfactory representation (i.e., mitral cell ensemble activity), in a manner consistent with experimental data and likely to facilitate the recognition of odor representations across concentrations. Signals reflecting absolute odor concentration could still be discerned from other response properties such as mitral cell spike latencies or tufted cell activity ; however, this question was not addressed in the present report.
Contrast enhancement reflected in individual mitral cells
In order to replicate Figure 2A from Yokoi et al. , the canonical demonstration of olfactory bulb contrast enhancement, the ten-glomerulus model was adapted in two ways. First, a 0.5 Hz sinusoidal input current was injected into the mitral cell distal glomerular tuft to replicate the artificial respiration-associated background activity depicted in that paper. Odors were then presented normally. Second, the sequential presentation of a homologous series of odorants was modeled as follows. All molecular receptive ranges were modeled as normal distributions of odor ligand-receptor affinities. As each subsequent odorant in the homologous series was presented, the affinity of the glomerulus of primary interest for that odorant first increased and then decreased along the trajectory of a normal distribution. The affinities of the other nine glomeruli for each odorant were also drawn from normal distributions, but randomly rather than sequentially, simulating a population of glomeruli with unknown tuning curves, the activity of which would vary unpredictably with changes in odor stimulation. Similar results were obtained when the affinities of the other nine glomeruli were simply held constant (data not shown). Data were high-pass filtered after generation to emphasize spiking activity and hence resemble extracellular recordings. Yokoi et al.  also illustrated that blockade of bulbar GABAA receptors with bicuculline could transform an mitral cell inhibitory odor response into a neutral or excitatory response. Using the ten-glomerulus model, we applied an odorant stimulus with moderate affinity for our glomerulus of interest, evoking an inhibitory response. We then blocked periglomerular synapses onto mitral cells and applied this stimulus again in order to replicate Figure 5D from Yokoi et al. .
Mori K, Shepherd GM: Emerging principles of molecular signal processing by mitral/tufted cells in the olfactory bulb. Semin Cell Biol. 1994, 5 (1): 65-74. 10.1006/scel.1994.1009.
Aungst JL, Heyward PM, Puche AC, Karnup SV, Hayar A, Szabo G, Shipley MT: Centre-surround inhibition among olfactory bulb glomeruli. Nature. 2003, 426 (6967): 623-629. 10.1038/nature02185.
Hayar A, Karnup SV, Shipley MT, Ennis M: Olfactory bulb glomeruli: external tufted cells intrinsically burst at theta frequency and are entrained by patterned olfactory input. J Neurosci. 2004, 24 (5): 1190-1199. 10.1523/JNEUROSCI.4714-03.2004.
Cleland TA, Linster C: Central olfactory processing. Handbook of olfaction and gustation, 2nd ed. Edited by: Doty RL. 2003, New York , Marcel Dekker, 165-180.
Shipley MT, McLean JH, Ennis M: Olfactory system. The rat nervous system. Edited by: Paxinos G. 1995, San Diego , Academic Press, 899-928. 2nd
Mountcastle VB: The columnar organization of the neocortex. Brain. 1997, 120 ( Pt 4): 701-722. 10.1093/brain/120.4.701.
Cook PB, McReynolds JS: Lateral inhibition in the inner retina is important for spatial tuning of ganglion cells. Nat Neurosci. 1998, 1 (8): 714-719. 10.1038/3714.
Suga N, Zhang Y, Yan J: Sharpening of frequency tuning by inhibition in the thalamic auditory nucleus of the mustached bat. J Neurophysiol. 1997, 77 (4): 2098-2114.
Yang L, Pollak GD, Resler C: GABAergic circuits sharpen tuning curves and modify response properties in the mustache bat inferior colliculus. J Neurophysiol. 1992, 68 (5): 1760-1774.
Yokoi M, Mori K, Nakanishi S: Refinement of odor molecule tuning by dendrodendritic synaptic inhibition in the olfactory bulb. Proc Natl Acad Sci U S A. 1995, 92 (8): 3371-3375.
Linster C, Sachse S, Galizia G: Computational modeling suggests that response properties rather than spatial position determine connectivity between olfactory glomeruli. J Neurophysiol. 2005, 93 (6): 3410-3417. 10.1152/jn.01285.2004.
Poirazi P, Mel BW: Impact of active dendrites and structural plasticity on the memory capacity of neural tissue. Neuron. 2001, 29 (3): 779-796. 10.1016/S0896-6273(01)00252-5.
Chklovskii DB, Mel BW, Svoboda K: Cortical rewiring and information storage. Nature. 2004, 431 (7010): 782-788. 10.1038/nature03012.
Tsuboi A, Yoshihara S, Yamazaki N, Kasai H, Asai-Tsuboi H, Komatsu M, Serizawa S, Ishii T, Matsuda Y, Nagawa F, Sakano H: Olfactory neurons expressing closely linked and homologous odorant receptor genes tend to project their axons to neighboring glomeruli on the olfactory bulb. J Neurosci. 1999, 19 (19): 8409-8418.
Zheng C, Feinstein P, Bozza T, Rodriguez I, Mombaerts P: Peripheral olfactory projections are differentially affected in mice deficient in a cyclic nucleotide-gated channel subunit. Neuron. 2000, 26 (1): 81-91. 10.1016/S0896-6273(00)81140-X.
Tozaki H, Tanaka S, Hirata T: Theoretical consideration of olfactory axon projection with an activity-dependent neural network model. Mol Cell Neurosci. 2004, 26 (4): 503-517. 10.1016/j.mcn.2004.04.006.
Serizawa S, Ishii T, Nakatani H, Tsuboi A, Nagawa F, Asano M, Sudo K, Sakagami J, Sakano H, Ijiri T, Matsuda Y, Suzuki M, Yamamori T, Iwakura Y: Mutually exclusive expression of odorant receptor transgenes. Nat Neurosci. 2000, 3 (7): 687-693. 10.1038/76641.
Ishii T, Serizawa S, Kohda A, Nakatani H, Shiroishi T, Okumura K, Iwakura Y, Nagawa F, Tsuboi A, Sakano H: Monoallelic expression of the odourant receptor gene and axonal projection of olfactory sensory neurones. Genes Cells. 2001, 6 (1): 71-78. 10.1046/j.1365-2443.2001.00398.x.
Hudson R: From molecule to mind: the role of experience in shaping olfactory function. J Comp Physiol [A]. 1999, 185 (4): 297-304. 10.1007/s003590050390.
Korsching SI: Odor maps in the brain: spatial aspects of odor representation in sensory surface and olfactory bulb. Cell Mol Life Sci. 2001, 58 (4): 520-530.
Alkasab TK, White J, Kauer JS: A computational system for simulating and analyzing arrays of biological and artificial chemical sensors. Chem Senses. 2002, 27 (3): 261-275. 10.1093/chemse/27.3.261.
Dickinson TA, Michael KL, Kauer JS, Walt DR: Convergent, self-encoded bead sensor arrays in the design of an artificial nose. Anal Chem. 1999, 71 (11): 2192-2198. 10.1021/ac981457i.
Mombaerts P: How smell develops. Nat Neurosci. 2001, 4 Suppl: 1192-1198. 10.1038/nn751.
Ressler KJ, Sullivan SL, Buck LB: Information coding in the olfactory system: evidence for a stereotyped and highly organized epitope map in the olfactory bulb. Cell. 1994, 79 (7): 1245-1255. 10.1016/0092-8674(94)90015-9.
Mombaerts P, Wang F, Dulac C, Chao SK, Nemes A, Mendelsohn M, Edmondson J, Axel R: Visualizing an olfactory sensory map. Cell. 1996, 87 (4): 675-686. 10.1016/S0092-8674(00)81387-2.
Linster C, Johnson BA, Yue E, Morse A, Xu Z, Hingco EE, Choi Y, Choi M, Messiha A, Leon M: Perceptual correlates of neural representations evoked by odorant enantiomers. J Neurosci. 2001, 21 (24): 9837-9843.
Cleland TA, Morse A, Yue EL, Linster C: Behavioral models of odor similarity. Behav Neurosci. 2002, 116 (2): 222-231. 10.1037/0735-7044.116.2.222.
Mouly AM, Gervais R, Holley A: Evidence for the involvement of rat olfactory bulb in processes supporting long-term olfactory memory. Eur J Neurosci. 1990, 2 (11): 978-984.
Kohonen T: Self-organized formation of topology correct feature maps. Biol Cybern. 1982, 43: 59-69. 10.1007/BF00337288.
Kohonen T, Hari R: Where the abstract feature maps of the brain might come from. Trends Neurosci. 1999, 22 (3): 135-139. 10.1016/S0166-2236(98)01342-3.
Penfield W, Rasmussen T: The cerebral cortex of man: a clinical study of localization of function. 1950, New York , MacMillan
Kell CA, von Kriegstein K, Roesler A, Kleinschmidt A, Laufs H: The sensory cortical representation of the human penis: revisiting somatotopy in the male homunculus. J Neurosci. 2005, 25 (25): 5984-5987. 10.1523/JNEUROSCI.0712-05.2005.
Igarashi KM, Mori K: Spatial representation of hydrocarbon odorants in the ventrolateral zones of the rat olfactory bulb. J Neurophysiol. 2005, 93 (2): 1007-1019. 10.1152/jn.00873.2004.
Johnson BA, Woo CC, Hingco EE, Pham KL, Leon M: Multidimensional chemotopic responses to n-aliphatic acid odorants in the rat olfactory bulb. J Comp Neurol. 1999, 409 (4): 529-548. 10.1002/(SICI)1096-9861(19990712)409:4<529::AID-CNE2>3.0.CO;2-N.
Johnson BA, Leon M: Odorant molecular length: one aspect of the olfactory code. J Comp Neurol. 2000, 426 (2): 330-338. 10.1002/1096-9861(20001016)426:2<330::AID-CNE12>3.0.CO;2-5.
Schoenfeld TA, Cleland TA: The anatomical logic of smell. Trends Neurosci. 2005, 28 (11): 620-627. 10.1016/j.tins.2005.09.005.
Bozza T, McGann JP, Mombaerts P, Wachowiak M: In vivo imaging of neuronal activity by targeted expression of a genetically encoded probe in the mouse. Neuron. 2004, 42 (1): 9-21. 10.1016/S0896-6273(04)00144-8.
Meister M, Bonhoeffer T: Tuning and topography in an odor map on the rat olfactory bulb. J Neurosci. 2001, 21 (4): 1351-1360.
Belluscio L, Katz LC: Symmetry, stereotypy, and topography of odorant representations in mouse olfactory bulbs. J Neurosci. 2001, 21 (6): 2113-2122.
Leon M, Johnson BA: Olfactory coding in the mammalian olfactory bulb. Brain Res Brain Res Rev. 2003, 42 (1): 23-32. 10.1016/S0165-0173(03)00142-5.
Johnson BA, Ho SL, Xu Z, Yihan JS, Yip S, Hingco EE, Leon M: Functional mapping of the rat olfactory bulb using diverse odorants reveals modular responses to functional groups and hydrocarbon structural features. J Comp Neurol. 2002, 449 (2): 180-194. 10.1002/cne.10284.
Uchida N, Takahashi YK, Tanifuji M, Mori K: Odor maps in the mammalian olfactory bulb: domain organization and odorant structural features. Nat Neurosci. 2000, 3 (10): 1035-1043. 10.1038/79857.
Takahashi YK, Kurosaki M, Hirono S, Mori K: Topographic representation of odorant molecular features in the rat olfactory bulb. J Neurophysiol. 2004
Friedrich RW, Korsching SI: Combinatorial and chemotopic odorant coding in the zebrafish olfactory bulb visualized by optical imaging. Neuron. 1997, 18 (5): 737-752. 10.1016/S0896-6273(00)80314-1.
Wachowiak M, Cohen LB: Representation of odorants by receptor neuron input to the mouse olfactory bulb. Neuron. 2001, 32 (4): 723-735. 10.1016/S0896-6273(01)00506-2.
Shepherd GM, Greer CA: Olfactory bulb. The synaptic organization of the brain. Edited by: Shepherd GM. 1998, New York , Oxford University Press, 159-204. 4th
Pinching AJ, Powell TP: The neuropil of the glomeruli of the olfactory bulb. J Cell Sci. 1971, 9 (2): 347-377.
White EL: Synaptic organization in the olfactory glomerulus of the mouse. Brain Res. 1972, 37 (1): 69-80. 10.1016/0006-8993(72)90346-0.
Hines ML, Carnevale NT: The NEURON simulation environment. Neural Comput. 1997, 9 (6): 1179-1209.
Hines ML, Carnevale NT: NEURON: a tool for neuroscientists. Neuroscientist. 2001, 7 (2): 123-135.
Chen WR, Shen GY, Shepherd GM, Hines ML, Midtgaard J: Multiple modes of action potential initiation and propagation in mitral cell primary dendrite. J Neurophysiol. 2002, 88 (5): 2755-2764.
Shen GY, Chen WR, Midtgaard J, Shepherd GM, Hines ML: Computational analysis of action potential initiation in mitral cell soma and dendrites based on dual patch recordings. J Neurophysiol. 1999, 82 (6): 3006-3020.
Ezeh PI, Wellis DP, Scott JW: Organization of inhibition in the rat olfactory bulb external plexiform layer. J Neurophysiol. 1993, 70 (1): 263-274.
McQuiston AR, Katz LC: Electrophysiology of interneurons in the glomerular layer of the rat olfactory bulb. J Neurophysiol. 2001, 86 (4): 1899-1907.
Pinching AJ, Powell TP: The neuron types of the glomerular layer of the olfactory bulb. J Cell Sci. 1971, 9 (2): 305-345.
Pinching AJ, Powell TP: The neuropil of the periglomerular region of the olfactory bulb. J Cell Sci. 1971, 9 (2): 379-409.
Wellis DP, Scott JW: Intracellular responses of identified rat olfactory bulb interneurons to electrical and odor stimulation. J Neurophysiol. 1990, 64 (3): 932-947.
Mel BW, Schiller J: On the fight between excitation and inhibition: location is everything. Sci STKE. 2004, 2004 (250): PE44-10.1126/stke.2502004pe44.
Liu G: Local structural balance and functional interaction of excitatory and inhibitory synapses in hippocampal dendrites. Nat Neurosci. 2004, 7 (4): 373-379. 10.1038/nn1206.
Koch C, Poggio T, Torre V: Nonlinear interactions in a dendritic tree: localization, timing, and role in information processing. Proc Natl Acad Sci U S A. 1983, 80 (9): 2799-2802.
Chen WR, Shepherd GM: Membrane and synaptic properties of mitral cells in slices of rat olfactory bulb. Brain Res. 1997, 745: 189-196. 10.1016/S0006-8993(96)01150-X.
Hamilton KA, Kauer JS: Patterns of intracellular potentials in salamander mitral/tufted cells in response to odor stimulation. J Neurophysiol. 1989, 62 (3): 609-625.
Kauer JS, Hamilton KA, Neff SR, Cinelli AR: Chemosensory information processing. Edited by: Schild D. 1990, Berlin , Springer-Verlag, H 39: 305-314. Temporal patterns of membrane potential in the olfactory bulb observed with intracellular recording and voltage-sensitive dye imaging: early hyperpolarization, NATO ASI Series.
Carandini M, Heeger DJ: Summation and division by neurons in primate visual cortex. Science. 1994, 264 (5163): 1333-1336.
Holt GR, Koch C: Shunting inhibition does not have a divisive effect on firing rates. Neural Comput. 1997, 9 (5): 1001-1013.
Spratling MW, Johnson MH: Preintegration lateral inhibition enhances unsupervised learning. Neural Comput. 2002, 14 (9): 2157-2179. 10.1162/089976602320264033.
Spratling MW, Johnson MH: Dendritic inhibition enhances neural coding properties. Cereb Cortex. 2001, 11 (12): 1144-1149. 10.1093/cercor/11.12.1144.
Johnson BA, Leon M: Modular representations of odorants in the glomerular layer of the rat olfactory bulb and the effects of stimulus concentration. J Comp Neurol. 2000, 422 (4): 496-509. 10.1002/1096-9861(20000710)422:4<496::AID-CNE2>3.0.CO;2-4.
Wachowiak M, Zochowski M, Cohen LB, Falk CX: The spatial representation of odors by olfactory receptor neuron input to the olfactory bulb is concentration invariant. Biol Bull. 2000, 199 (2): 162-163.
Rubin BD, Katz LC: Optical imaging of odorant representations in the mammalian olfactory bulb. Neuron. 1999, 23 (3): 499-511. 10.1016/S0896-6273(00)80803-X.
Harrison TA, Scott JW: Olfactory bulb responses to odor stimulation: analysis of response pattern and intensity relationships. J Neurophysiol. 1986, 56 (6): 1571-1589.
Wellis DP, Scott JW, Harrison TA: Discrimination among odorants by single neurons of the rat olfactory bulb. J Neurophysiol. 1989, 61 (6): 1161-1177.
Meredith M: Patterned response to odor in mammalian olfactory bulb: the influence of intensity. J Neurophysiol. 1986, 56 (3): 572-597.
Cang J, Isaacson JS: In vivo whole-cell recording of odor-evoked synaptic transmission in the rat olfactory bulb. J Neurosci. 2003, 23 (10): 4108-4116.
Cleland TA, Linster C: How synchronization properties among second-order sensory neurons can mediate stimulus salience. Behav Neurosci. 2002, 116 (2): 212-221. 10.1037/0735-7044.116.2.212.
Nagayama S, Takahashi YK, Yoshihara Y, Mori K: Mitral and tufted cells differ in the decoding manner of odor maps in the rat olfactory bulb. J Neurophysiol. 2004, 91 (6): 2532-2540. 10.1152/jn.01266.2003.
Schoppa NE, Urban NN: Dendritic processing within olfactory bulb circuits. Trends Neurosci. 2003, 26 (9): 501-506. 10.1016/S0166-2236(03)00228-5.
Mori K, Nagao H, Yoshihara Y: The olfactory bulb: coding and processing of odor molecule information. Science. 1999, 286 (5440): 711-715. 10.1126/science.286.5440.711.
Margrie TW, Sakmann B, Urban NN: Action potential propagation in mitral cell lateral dendrites is decremental and controls recurrent and lateral inhibition in the mammalian olfactory bulb. Proc Natl Acad Sci U S A. 2001, 98 (1): 319-324. 10.1073/pnas.011523098.
Rall W, Shepherd GM, Reese TS, Brightman MW: Dendrodendritic synaptic pathway for inhibition in the olfactory bulb. Exp Neurol. 1966, 14 (1): 44-56. 10.1016/0014-4886(66)90023-9.
Linster C, 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. 1996, 3 (3): 225-246. 10.1007/BF00161133.
Linster C, Hasselmo M: Modulation of inhibition in a model of olfactory bulb reduces overlap in the neural representation of olfactory stimuli. Behav Brain Res. 1997, 84 (1-2): 117-127. 10.1016/S0166-4328(97)83331-1.
Shepherd GM, Brayton RK: Computer simulation of a dendrodendritic synaptic circuit for self- and lateral-inhibition in the olfactory bulb. Brain Res. 1979, 175 (2): 377-382. 10.1016/0006-8993(79)91020-5.
Davison AP, Feng J, Brown D: Dendrodendritic inhibition and simulated odor responses in a detailed olfactory bulb network model. J Neurophysiol. 2003, 90 (3): 1921-1935.
DeVries SH, Baylor DA: Synaptic circuitry of the retina and olfactory bulb. Cell. 1993, 72 Suppl: 139-149. 10.1016/S0092-8674(05)80033-9.
Sachse S, Galizia CG: The coding of odour-intensity in the honeybee antennal lobe: local computation optimizes odour representation. Eur J Neurosci. 2003, 18 (8): 2119-2132. 10.1046/j.1460-9568.2003.02931.x.
Egana JI, Aylwin ML, Maldonado PE: Odor response properties of neighboring mitral/tufted cells in the rat olfactory bulb. Neuroscience. 2005
Buonviso N, Chaput MA: Response similarity to odors in olfactory bulb output cells presumed to be connected to the same glomerulus: electrophysiological study using simultaneous single-unit recordings. J Neurophysiol. 1990, 63 (3): 447-454.
Linster C, Hasselmo ME: Behavioral responses to aliphatic aldehydes can be predicted from known electrophysiological responses of mitral cells in the olfactory bulb. Physiol Behav. 1999, 66 (3): 497-502. 10.1016/S0031-9384(98)00324-2.
Laska M, Teubner P: Odor structure-activity relationships of carboxylic acids correspond between squirrel monkeys and humans. Am J Physiol. 1998, 274 (6 Pt 2): R1639-45.
Luo M, Katz LC: Response correlation maps of neurons in the mammalian olfactory bulb. Neuron. 2001, 32 (6): 1165-1179. 10.1016/S0896-6273(01)00537-2.
Debarbieux F, Audinat E, Charpak S: Action potential propagation in dendrites of rat mitral cells in vivo. J Neurosci. 2003, 23 (13): 5553-5560.
Lagier S, Carleton A, Lledo PM: Interplay between local GABAergic interneurons and relay neurons generates gamma oscillations in the rat olfactory bulb. J Neurosci. 2004, 24 (18): 4382-4392. 10.1523/JNEUROSCI.5570-03.2004.
Djurisic M, Antic S, Chen WR, Zecevic D: Voltage imaging from dendrites of mitral cells: EPSP attenuation and spike trigger zones. J Neurosci. 2004, 24 (30): 6703-6714. 10.1523/JNEUROSCI.0307-04.2004.
Cecchi GA, Petreanu LT, Alvarez-Buylla A, Magnasco MO: Unsupervised learning and adaptation in a model of adult neurogenesis. J Comput Neurosci. 2001, 11 (2): 175-182. 10.1023/A:1012849801892.
Rochefort C, Gheusi G, Vincent JD, Lledo PM: Enriched odor exposure increases the number of newborn neurons in the adult olfactory bulb and improves odor memory. J Neurosci. 2002, 22 (7): 2679-2689.
Petreanu L, Alvarez-Buylla A: Maturation and death of adult-born olfactory bulb granule neurons: role of olfaction. J Neurosci. 2002, 22 (14): 6106-6113.
Mandairon N, Jourdan F, Didier A: Deprivation of sensory inputs to the olfactory bulb up-regulates cell death and proliferation in the subventricular zone of adult mice. Neuroscience. 2003, 119 (2): 507-516. 10.1016/S0306-4522(03)00172-6.
Kashiwadani H, Sasaki YF, Uchida N, Mori K: Synchronized oscillatory discharges of mitral/tufted cells with different molecular receptive ranges in the rabbit olfactory bulb. J Neurophysiol. 1999, 82 (4): 1786-1792.
Nusser Z, Kay LM, Laurent G, Homanics GE, Mody I: Disruption of GABA(A) receptors on GABAergic interneurons leads to increased oscillatory power in the olfactory bulb network. J Neurophysiol. 2001, 86 (6): 2823-2833.
Martin C, Gervais R, Hugues E, Messaoudi B, Ravel N: Learning modulation of odor-induced oscillatory responses in the rat olfactory bulb: a correlate of odor recognition?. J Neurosci. 2004, 24 (2): 389-397. 10.1523/JNEUROSCI.3433-03.2004.
Kay LM: Two species of gamma oscillations in the olfactory bulb: dependence on behavioral state and synaptic interactions. J Integr Neurosci. 2003, 2 (1): 31-44. 10.1142/S0219635203000196.
Linster C, Cleland TA: How spike synchronization among olfactory neurons can contribute to sensory discrimination. J Comput Neurosci. 2001, 10 (2): 187-193. 10.1023/A:1011221131212.
Koster NL, Norman AB, Richtand NM, Nickell WT, Puche AC, Pixley SK, Shipley MT: Olfactory receptor neurons express D2 dopamine receptors. J Comp Neurol. 1999, 411 (4): 666-673. 10.1002/(SICI)1096-9861(19990906)411:4<666::AID-CNE10>3.0.CO;2-S.
Nickell WT, Norman AB, Wyatt LM, Shipley MT: Olfactory bulb DA receptors may be located on terminals of the olfactory nerve. Neuroreport. 1991, 2 (1): 9-12.
Bonino M, Cantino D, Sassoe-Pognetto M: Cellular and subcellular localization of gamma-aminobutyric acidB receptors in the rat olfactory bulb. Neurosci Lett. 1999, 274 (3): 195-198. 10.1016/S0304-3940(99)00697-7.
Palouzier-Paulignan B, Duchamp-Viret P, Hardy AB, Duchamp A: GABA(B) receptor-mediated inhibition of mitral/tufted cell activity in the rat olfactory bulb: a whole-cell patch-clamp study in vitro. Neuroscience. 2002, 111 (2): 241-250. 10.1016/S0306-4522(02)00003-9.
Ennis M, Zhou FM, Ciombor KJ, Aroniadou-Anderjaska V, Hayar A, Borrelli E, Zimmer LA, Margolis F, Shipley MT: Dopamine D2 receptor-mediated presynaptic inhibition of olfactory nerve terminals. J Neurophysiol. 2001, 86 (6): 2986-2997.
Aroniadou-Anderjaska V, Ennis M, Shipley MT: Dendrodendritic recurrent excitation in mitral cells of the rat olfactory bulb. J Neurophysiol. 1999, 82 (1): 489-494.
Salin PA, Lledo PM, Vincent JD, Charpak S: Dendritic glutamate autoreceptors modulate signal processing in rat mitral cells. J Neurophysiol. 2001, 85 (3): 1275-1282.
Schoppa NE, Westbrook GL: AMPA autoreceptors drive correlated spiking in olfactory bulb glomeruli. Nat Neurosci. 2002, 5 (11): 1194-1202. 10.1038/nn953.
Urban NN, Sakmann B: Reciprocal intraglomerular excitation and intra- and interglomerular lateral inhibition between mouse olfactory bulb mitral cells. J Physiol. 2002, 542 (Pt 2): 355-367. 10.1113/jphysiol.2001.013491.
Smith TC, Jahr CE: Self-inhibition of olfactory bulb neurons. Nat Neurosci. 2002, 5 (8): 760-766.
Eisthen HL: Why are olfactory systems of different animals so similar?. Brain Behav Evol. 2002, 59 (5-6): 273-293. 10.1159/000063564.
Cleland TA, Linster C: Concentration tuning mediated by spare receptor capacity in olfactory sensory neurons: A theoretical study. Neural Comput. 1999, 11 (7): 1673-1690. 10.1162/089976699300016188.
Kang K, Shelley M, Sompolinsky H: Mexican hats and pinwheels in visual cortex. Proc Natl Acad Sci U S A. 2003, 100 (5): 2848-2853. 10.1073/pnas.0138051100.
Anton PS, Granger R, Lynch G: Simulated dendritic spines influence reciprocal synaptic strengths and lateral inhibition in the olfactory bulb. Brain Res. 1993, 628 (1-2): 157-165. 10.1016/0006-8993(93)90951-I.
Hines ML, Carnevale NT: Expanding NEURON's repertoire of mechanisms with NMODL. Neural Comput. 2000, 12 (5): 995-1007. 10.1162/089976600300015475.
Destexhe A, Neubig M, Ulrich D, Huguenard J: Dendritic low-threshold calcium currents in thalamic relay cells. J Neurosci. 1998, 18 (10): 3574-3588.
Destexhe A, Contreras D, Steriade M, Sejnowski TJ, Huguenard JR: In vivo, in vitro, and computational analysis of dendritic calcium currents in thalamic reticular neurons. J Neurosci. 1996, 16 (1): 169-185.
Cadetti L, Belluzzi O: Hyperpolarisation-activated current in glomerular cells of the rat olfactory bulb. Neuroreport. 2001, 12 (14): 3117-3120. 10.1097/00001756-200110080-00027.
Migliore M, Cook EP, Jaffe DB, Turner DA, Johnston D: Computer simulations of morphologically reconstructed CA3 hippocampal neurons. J Neurophysiol. 1995, 73 (3): 1157-1168.
Coenen OJMD, Sejnowski TJ, Jabri MA: Parallel fiber coding in the cerebellum for life-long learning. Autonomous Robots. 2001, 11: 291-297. 10.1023/A:1012403510221.
Urban NN: Lateral inhibition in the olfactory bulb and in olfaction. Physiol Behav. 2002, 77 (4-5): 607-612. 10.1016/S0031-9384(02)00895-8.
We gratefully acknowledge Christiane Linster and Mary Lou Zeeman for technical discussions and comments, Michael Hines and Ted Carnevale for technical advice regarding NEURON, and Laura Hurley, Michael Leon, Ann Marie McNamara, Adam Puche, Merri Rosen, Thomas Schoenfeld, and Michael Shipley for comments on earlier versions of the manuscript. This work was supported by NIH grant DC005727 to ThAC.
ThAC conceived and designed the study, coded and analyzed models in Matlab and NEURON, and drafted the manuscript. PS coded and analyzed models in NEURON.