Previous Article |
Table of Contents
| Next Article
BIOLOGICAL SCIENCES / NEUROSCIENCE
Relational representation in the olfactory system

,
*Department of Neurobiology and Behavior, Cornell University, Ithaca, NY 14853; and
Department of Neurobiology and Behavior, University of California, Irvine, CA 92697
Edited by Linda M. Bartoshuk, University of Florida, Gainesville, FL, and approved December 13, 2006 (received for review October 4, 2006)
| Abstract |
|---|
|
|
|---|
coding | computational neuroscience | glomerulus | mitral cell | short | axon cell
Relational representations depend on the broad patterns of primary sensory activation that typically are observed in response to elemental stimuli. For example, the movement of a single vibrissa evokes responses across a broad region of somatosensory cortex (1), a single point of light evokes widespread activity in visual cortex (2), and a single pure tone elicits a response across large areas of auditory cortex (3). Odorants evoke comparably broad responses in the mammalian olfactory system. Primary olfactory sensory neurons (OSNs) express a single olfactory receptor (OR) species in rodents, yet are broadly selective for odorant stimuli (46), such that even monomolecular odorants can evoke responses across much of the epithelium (7) and the glomerular layer of the olfactory bulb (OB) (8). Because of the convergence of OSNs expressing the same OR species onto common glomeruli, the pattern of activated glomeruli on the surface of the OB reflects the pattern of activated ORs, thus clearly identifying the constellation of chemical qualities that constitute the presented odor (811).
We investigated the possibility that relational representations enable the maintenance of odor quality over large differences in odorant concentration and the mechanism by which the olfactory system can maintain these relational representations. It has been well established that the degree of overlap among different glomerular activity maps reflects odorant structural commonalities and corresponds to the perceptual similarity of odors (1215). However, this relationship breaks down when odorant concentration is included as a variable. OSN activation profiles and the corresponding glomerular activity maps universally become broader and more intense when odors are presented at higher concentrations (Fig. 1A) (1620), simply because higher odorant concentrations recruit additional OSN populations with progressively lower affinities for the presented agonist(s).
|
The computation of intensity-independent stimulus representations in sensory systems depends on the existence of associated inhibitory systems that are activated more broadly by sensory stimuli than are the principal neurons receiving direct sensory inputs (2224). Specifically, inhibition that is uniformly distributed across the input field, perhaps representing an average or sum of all input activity, will tend to preserve the pattern of relative activation levels across the field irrespective of total intensity, whereas inhibition that is delivered in proportion to local activation levels, or in an otherwise biased manner, will retain concentration-dependent distortions in the resulting output patterns. Recent work in the OB (25, 26) has identified an interconnected network of excitatory juxtaglomerular interneurons, external tufted (ET) and short-axon (SA) cells, that delivers inhibition onto OB principal neurons (mitral cells) via local inhibitory neurons [periglomerular (PG) cells], and it has been suggested that this recurrent excitatory network mediates the uniform inhibition required for concentration-invariant secondary representations and other odor-processing mechanisms (22). We show here that this interglomerular excitatory network, via activation of local inhibitory neurons and with estimates of synaptic densities and connectivity derived from experimental data, is not only capable of computing the required normalization but also seems optimized to perform this operation. Using a full-scale, 2,200-glomerulus model to process and analyze 2-deoxyglucose imaging data from the rat OB, along with behavioral data, we demonstrate the feasibility and potential mechanism of perceptual concentration invariance in olfaction.
| Results |
|---|
|
|
|---|
x]. To measure the similarity between glomerular activation patterns, we calculated pairwise indices of dissimilarity (normalized Euclidean distance) between the raw and normalized patterns of activation evoked by four odorants (2-hexanone, methyl valerate, n-pentanal, and n-pentanol) at vapor-phase concentrations of 25, 75, and 250 ppm. The indices of dissimilarity between the normalized patterns evoked by any given odorant across concentrations were reduced significantly as compared with those of the corresponding raw patterns (P < 0.05 in all cases; Fig. 1 C and D), indicating that normalized patterns represent concentration-invariant quality information significantly better than do raw response patterns. In contrast, dissimilarity indices between pairs of different odorants at the same concentration were not reduced consistently by normalization; effects varied from a marginal reduction to a substantial increase in dissimilarity (Fig. 1E). If normalized glomerular activity patterns represent the output of bulbar computations performed on raw input patterns, then they should be better predictors of perceptual similarity than the raw patterns. Indeed, quantitative comparisons between glomerular activation patterns and olfactory perception generally have been based on normalized data (1214,27,28). To test this hypothesis, we selected odorant pairs in which the raw glomerular patterns evoked by the two odorants at a given concentration are more similar to one another than are those evoked by different concentrations of the same odorant (Fig. 2A). In contrast, when the activity patterns of these odorants are normalized, the representations of each odorant are more similar to one another, irrespective of concentration, than they are to those evoked by the other odorant at any concentration (Fig. 2B). We tested which of these predictions of perceptual similarity was supported by behavioral data by using a spontaneous, nonassociative odor-discrimination test (1214). Briefly, we first habituated rats to the presentation of one odorant of each pair (2-hexanone or n-pentanal) at 25 ppm and then measured their novelty responses to presentation of (i) the same odorant at 75 ppm and (ii) the paired odorant (methyl valerate or n-pentanol, respectively) at 25 ppm (Fig. 2C). The results show that rats perceive the two concentrations of either odorant as more similar to one another than to the different odorant (Fig. 2D), demonstrating that normalized glomerular representations predict perceptual similarity whereas raw patterns do not.
|
We constructed a detailed, full-scale, 2,200-glomerulus model of the glomerular layer of the OB to study how normalization could be computed. OSNs and PG, ET, SA, and mitral cells were represented and synaptically connected as described experimentally (Fig. 3A) (25, 26, 34). Briefly, each population of OSNs converging on a given glomerulus was represented by a single aggregate OSN; hence, the model contained 2,200 such OSNs. For any given odor, the activation level of each OSN was drawn from a corresponding data point collected in 2-deoxyglucose imaging experiments, enabling the model to be presented with receptor activation profiles corresponding to those evoked by real odorant stimuli (absent any temporal component). Within each glomerulus, OSNs synaptically excited PG, ET, and mitral cells, whereas PG cells in turn directly inhibited mitral cells. ET cells excited both PG cells associated with the same glomerulus and SA cells within a 2-glomerulus radius (25, 26). SA cells did not receive monosynaptic input from OSNs but were activated by ET cells; they in turn excited ET and PG cells in several neighboring glomeruli. Importantly, the mutual excitation between ET and SA neurons constituted a recurrent excitatory network extending laterally across the glomerular layer, the output of which was delivered to PG cells and, subsequently, to mitral cells as inhibition (26). The strength, density, and physical distribution of these recurrent connections among glomeruli were variables of interest (Fig. 3B).
|
Optimization of Local Inhibitory Synaptic Strength. Normalization with respect to the mean requires subtraction of the average input strength from the representations of individual stimulus features. Hence, the network must both compute an estimate of this average input strength and scale this average appropriately for delivery onto mitral cells as inhibition. We tested the sufficiency of the interglomerular network to meet this criterion. Presentation of an odorant stimulus to the model evoked a characteristic pattern of activation among OSNs (Fig. 4Ai). OSN activity in turn excited ET, SA, and PG interneurons (directly or indirectly) as well as mitral cells. Critically, lateral excitatory projections mediated by ET/SA cells acted to reduce the variance in activation levels among the interneurons associated with different glomeruli (Fig. 4Aii).
|
Optimization of Lateral Excitatory Network Properties.
The SA cell is the main contributor to interglomerular connectivity (26). Interglomerular projections from SA interneurons radiate in all directions from the soma. The density of these processes declines with distance, indicating a center-surround topological organization. Specifically, based on retrograde dye transport from focal injections into individual glomeruli, it has been estimated that half of the interneurons projecting to a given glomerulus are located at distances within 57 glomerular diameters (
350 µm) of the injection site, whereas only 10% are located at distances farther than 1518 glomerular diameters (
850 µm). Functional imaging has shown that the lateral spread of excitation from focal glomerular stimulation extends up to 900 ± 93 µm from the stimulus site, covering nearly the entirety of the visible OB. The functional effect of this lateral excitatory network on mitral cells is inhibitory; stimulation of interglomerular connections abruptly (6.5 ± 1.2-ms latency) terminates long-lasting depolarizations in mitral cells located 150525 µm away and reduces the amplitudes of newly generated long-lasting depolarizations via a GABAAergic mechanism (26).
The ability of this interglomerular network to mediate the uniform global inhibition of mitral cells relies critically on the connectivity of the ET/SA recurrent excitatory network. Recurrent excitatory networks, if activity is bounded or balanced so as to avoid runaway positive feedback, have the property of equalizing the level of activity among interconnected neurons. Consequently, if the lateral projections of SA cells and their mutually excitatory synapses with ET and other SA cells are sufficiently extensive to approximate a fully connected network, then their inhibition of mitral cells via PG cells, in principle, would be uniform across the entire glomerular layer, reflecting the average activity of all OSNs. We asked whether experimental estimates of ET/SA lateral connectivity among glomeruli could support this operation.
Initial parameters in the model were set to replicate experimental estimates of ET/SA connectivity. Each SA cell had a baseline P = 0.1 probability to synapse on ET, SA, and PG cells in other glomeruli, modulated as a function of distance d along the trajectory of a Gaussian distribution [fgauss(R, d)] with a standard deviation R of 7 glomerular diameters (Fig. 3B). These parameters generated appropriate connectivity maps in which
50% of all connections to a given cell originated from SA cells located within 57 glomerular diameters and 10% originated from SA cells located 1518 glomeruli away. The strength of the PG-to-mitral inhibitory synapse was set to its optimum (Fig. 4B). We then varied the baseline synaptic probability (P) and projection distance (R) parameters of the ET/SA network systematically in the model. Increases in either of these parameter values increased the average number of synapses that each cell received from SA cells [SA synaptic density; Fig. 5A; also see supporting information (SI) Fig. 6]. To ensure that only the distribution, and not the overall strength, of the interglomerular excitatory network was varied, individual synaptic weights were reduced as synaptic probabilities, and projection distances were increased so that the average total weight of excitatory synapses received by a given cell type was held constant. The SA synaptic density thereby directly reflected the level of interglomerular connectivity. The standard deviation of SA cell activity across the glomerular layer was large in the absence of lateral connectivity, reflecting the diverse odor-activation levels of different glomeruli, but converged rapidly toward an asymptotic minimum when the average SA synaptic density exceeded roughly four SA synapses per postsynaptic neuron (Fig. 5B) irrespective of stimulus concentration (SI Fig. 7). The model output more closely approximated a z score transformation of input activity as ET/SA activity became more uniform across the glomerular layer (Fig. 5C). Critically, the experimentally estimated values for interglomerular connectivity correspond precisely with the minimum level of connectivity required to achieve asymptotically uniform global feed-forward inhibition (Fig. 5B and SI Fig. 7, dashed vertical lines). This finding has two major implications. First, it suggests an optimization for uniform global inhibition because the minimum synaptic density necessary to achieve this effect appears to have been deployed in the OB. Second, it demonstrates that a center-surround projection architecture, such as that exhibited by SA neurons (26), does not necessarily imply center-surround functionality.
|
| Discussion |
|---|
|
|
|---|
The OB recurrent excitatory network, comprising ET and SA neurons, is activated by odor presentation that excites ET cells along with mitral cells and a minority of PG cells (25, 26). The spread of excitation through the ET/SA network equalizes the activity of these cells across the glomerular layer at a level corresponding to the total input strength. Importantly, although the lateral projections of SA interneurons are distributed more densely among neighboring than among distant glomeruli, demonstrating a type of center-surround architecture in the glomerular layer (26), our results show that this localized architecture does not necessarily imply center-surround functionality. Rather, the recurrent excitatory network generates a spatially uniform level of activity across the glomerular layer, perhaps corresponding to the broad, diffuse activity observed in intrinsic imaging studies of the OB (18). Moreover, these results suggest that this network deploys the minimum number of synaptic connections and the narrowest distribution of projection distances sufficient to produce such uniform activity. The net effect of this global feed-forward inhibition is to reduce significantly the impact of stimulus concentration on the activity level of the secondary olfactory representation, constructing a relational odor representation at the level of mitral cells that enables the perceptual constancy of odor quality across the ranges of concentrations inevitably encountered in natural environments.
| Materials and Methods |
|---|
|
|
|---|
A z score-normalized counterpart for each raw glomerular activation map was calculated by subtracting the mean value of all data points from each data point and then dividing each difference by the standard deviation of all data points, resulting in a transformed data set with a mean of zero and a standard deviation of unity. Raw and normalized glomerular activation maps, as well as model output maps, corresponding to four different odorants presented at three different concentrations (see SI Table 1) were compared by calculating the normalized Euclidean distances between all possible pairs of maps.
Perceptual Comparisons. We measured the perceptual similarity of odor pairs in rats (n = 10) by using a cross-habituation test (12, 13). Odorants were dissolved in mineral oil and presented by pipetting 60 µl of the odor stimulus onto a filter paper disk in a weighing dish that was placed on top of the wire cage cover. Test sessions consisted of one 50-sec presentation of vehicle followed by four 50-sec presentations of the habituation odorant at 5-min intervals and, finally, one 50-sec presentation of the test odorant. The time that the rat spent investigating the odorant (active sniffing within 1 cm) during each trial was recorded. Perceptual dissimilarity indices were calculated as the root-mean-square of the differences between investigation times during the final habituation trial (H4) and the test trial, normalized by the root-mean-square of the differences between H4 and the initial presentation of the habituation odorant (H1; Fig. 2C).
Computational Modeling. Our OB model (27) includes five neuron types: OSNs and PG, ET, SA, and mitral cells (Fig. 3A). Because these simulations were concerned exclusively with computations in the glomerular layer, granule cells and mitral cell lateral dendrites were omitted. All neurons were represented as single compartments, each characterized by a membrane time constant that can be regarded as the mean product of the membrane capacitance and input resistance. Odor inputs to the simulations were drawn directly from raw 2-deoxyglucose activity maps of the OB glomerular layer; i.e., real primary odorant representations were used as input.
The model simulated interactions among the major cell types described in the glomerular layer (Fig. 3A), with estimates of interglomerular connectivity derived from the literature (26). Synapses received from OSNs as well as those between PG and mitral cells were local to each glomerulus, whereas synapses received from SA cells, as well as those from ET to SA cells, were drawn from different glomeruli as a probabilistic function of interglomerular distance (Fig. 3B). For each pair of interacting cells, the probability of a synaptic connection was computed according to a probability density function defined by a baseline probability of connection and the breadth of the distribution of possible interactions (Fig. 3B). These distributions were modeled as Gaussian with a standard deviation equal to the breadth term R and a peak scaled to the baseline probability P of synaptic interactions. Modeling equations can be found in SI Text, and parameters are outlined in SI Table 2.
| Acknowledgements |
|---|
|
|
|---|
| Footnotes |
|---|
Abbreviations: OSN, olfactory sensory neuron; OB, olfactory bulb; OR, olfactory receptor; ET, external tufted; SA, short-axon; PG, periglomerular; Hn, habituation trial n.
To whom correspondence should be addressed. E-mail: mleon{at}uci.edu
Author contributions: T.A.C., B.A.J., M.L., and C.L. designed research, performed research, analyzed data, and wrote the paper.
The authors declare no conflict of interest.
This article contains supporting information online at www.pnas.org/cgi/content/full/0608564104/DC1
This article is a PNAS direct submission.
© 2007 by The National Academy of Sciences of the USA
| References |
|---|
|
|
|---|
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
C. A. Lowry and L. M. Kay Chemical Factors Determine Olfactory System Beta Oscillations in Waking Rats J Neurophysiol, July 1, 2007; 98(1): 394 - 404. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||