Abstract
We describe an integrated theory of olfactory systems operation that incorporates experimental findings across scales, stages, and methods of analysis into a common framework. In particular, we consider the multiple stages of olfactory signal processing as a collective system, in which each stage samples selectively from its antecedents. We propose that, following the signal conditioning operations of the nasal epithelium and glomerular-layer circuitry, the plastic external plexiform layer of the olfactory bulb effects a process of category learning—the basis for extracting meaningful, quasi-discrete odor representations from the metric space of undifferentiated olfactory quality. Moreover, this early categorization process also resolves the foundational problem of how odors of interest can be recognized in the presence of strong competitive interference from simultaneously encountered background odorants. This problem is fundamentally constraining on early-stage olfactory encoding strategies and must be resolved if these strategies and their underlying mechanisms are to be understood. Multiscale general theories of olfactory systems operation are essential in order to leverage the analytical advantages of engineered approaches together with our expanding capacity to interrogate biological systems.
Introduction
Theoretical models of neural systems serve as proofs of concept, tests of sufficiency, and quantitative embodiments of working hypotheses (Cleland and Linster, ). Well-vetted models impose discipline on hypothesis generation, and become particularly important when the systems under study are too complex to simply intuit. However, such models can reflect only the sources of variance that are adequately appreciated by theoreticians during model design. It is axiomatic in engineering that new designs are likely to require revision once physical prototypes are built, often because neglected minor variables generate significant unforeseen effects. This same principle applies to the reverse engineering of biological systems. Biological systems are physical systems that operate autonomously within insecure and incompletely predictable natural environments; their continued existence implies their adequate management of all relevant variables and successful exploitation of pertinent information. Identifying the properties that enable biological systems to adapt successfully to unmitigated natural variance is a core challenge both for our understanding of biological operations and for the design of autonomous robotic systems.
Here, we describe a set of physically explicit systems models of the early olfactory system that constrain and elucidate the neuronal networks that sample chemosensory data from the environment and present information extracted from these data to the central nervous system. This approach emphasizes an integrated appreciation of systems-level operations in situ: i.e., given a limited physical toolset of sensors and network elements, what data can be reliably acquired, what useful information can be plausibly extracted therefrom, and what conditions must be met in order to extract this useful information successfully? To date, this reverse engineering strategy has produced predictive models that have stood the test of time—despite their initial construction on relatively sparse direct evidence, and despite substantial revisions, owing to subsequent experimental studies, to our understanding of the cellular and network elements that directly underlie these computations. In particular, physical, engineered instantiations of biological systems models and the utilization of physical sensor data can present designed networks with the full spectrum of underappreciated difficulties that natural systems encounter routinely, thereby probing their capacities to manage unregulated stimuli encountered under unpredictable circumstances.
Three Core Problems in Engineering Sensory Systems
In sensory systems, this capacity for managing unpredictable input statistics is particularly important, because the external physical environment does not respect the limitations of biological systems, and the process of data sampling lacks the shared assumptions between sender and receiver that underlie the efficient encoding and transmission of information. Sources of potentially informative variance in the physical world are wildly diverse, and biological sensory systems must adapt to acquire the information that they need using tools with imperfect and idiosyncratic properties. These principles strongly constrain and determine the architecture of sensory systems. As a rough analogy, much of the circuitry in a commercial amplifier also is designed to manage unregulated external variance and compensate for the physical limitations of circuit elements so that the core amplifier elements can function properly and predictably in the real world. Neurons also exhibit limited bandwidth, notoriously non-linear response properties, and relatively narrow response ranges to sensory or synaptic input. It is to be expected that many of the circuit elements and properties of early sensory systems exist for the unlovely but essential processes of signal conditioning.
The core problems addressed by sensory systems can be divided into three phases:
Sampling the environmental variance of interest adequately, using available sensors, by any means necessary. In particular, the range of intensities of environmental signals often exceeds the response ranges of biological sensor elements by orders of magnitude. This renders it difficult to even encode afferent information, much less to organize it internally once encoded.
Identifying the particular aspects of the sampled variance that are relevant to a given sensory or behavioral task, and segregating this information from the uninteresting or confusing aspects of that variance.
Transforming the sensory information into a usable format or formats. This includes the problem of rendering this information relatively portable—i.e., usable by other brain systems and integrable with other modalities. Sometimes this may involve reframing such information into a common, abstract metric space.
Sampling Strategies
All sensory systems diversify their primary sensors in order to adequately sample the relevant environmental variance. Retinal cone photoreceptor sensitivities are diversified along the axis of wavelength, such that together they are able to detect a wider chromatic range of signals and, once detected, identify and distinguish specific colors (Birch and Wright, ). Adaptation to a wide range of environmental light levels is achieved by a set of strategies including iris contraction and the reliance on different photoreceptor populations during light-adapted and dark-adapted vision. Similarly, somatosensory tactile mechanoreceptors are diversified along axes of depth (under the skin) and adaptation timescale, as well as physical location on the body, so that in combination a wide range of contact modes, textures and other constructed tactile stimuli can be distinguished.
Nasal chemoreceptors are notoriously diverse, with mammals exhibiting hundreds of different olfactory receptor types in their noses that selectively respond to specific ranges of molecular ligands (Araneda et al., ). However, ligand concentration also strongly affects receptor binding and olfactory sensory neuron (OSN) activation (Schaefer et al., 2001; Storace and Cohen, 2017), substantially disrupting the naïve odor quality signal based on combinatorial activation patterns. Accordingly, a great deal of attention has been paid to the broad problem of concentration “invariance” in olfaction. However, there has been relatively little attention paid to the closely related problem imposed by the sharply limited dynamic response ranges of individual neurons. How can different odorant concentrations spanning multiple orders of magnitude even be sampled given the Boltzmann limitations on ligand-receptor binding properties (cooperativity ≥1), which restrict the dynamic input range of primary sensory receptors to roughly 1.5 orders of magnitude concentration? And, if this can be solved, how can this signal diversity then be organized and compressed back into a range that can be physically represented by downstream neural activity?
The answer to this question has not been fully experimentally demonstrated, but theoretical work has demarcated a set of strategies that could underlie this capacity. Briefly, if a large population of sibling OSNs (i.e., those expressing the same odorant receptor) could exhibit diverse EC50 values—hence each responding optimally to a different concentration range—then, in combination, their population activity generates a Boltzmann function with a reduced cooperativity and a correspondingly extended dynamic range across concentrations. This diversity can be achieved without loss of specificity if primary OSNs exhibit a range of spare receptor capacities (Cleland and Linster, ). The population activity then can be extracted from the glomeruli of the main olfactory bulb (MOB), upon which the axons of sibling OSNs converge (Figure 1). Indeed, optical imaging studies have demonstrated that the dynamic ranges of glomerular activation (effectively presynaptic population recordings from convergent OSN axons) extend across many orders of magnitude of odorant concentration, even as the corresponding postsynaptic activity in the apical dendrites of second-order principal neurons (mitral and tufted cells; MTCs) remains relatively stable (Storace and Cohen, 2017). These results reveal that the primary OSN population can successfully encode an extensive intensity range of environmental variance, encompassing diverse qualities, and concentrations of some number of odorous ligands. A series of established neurophysiological mechanisms then can compress this wide-range population signal into the limited dynamic range exhibited by activated MTCs (reviewed in Cleland et al., ; Cleland, ).
Figure 1
Selecting Relevant Variance
The processes by which primary sensory representations are mined for useful information are hugely diverse and challenging to fully understand. Primary representations generally reflect information derived from multiple conflated causes and sources and are correspondingly difficult to parse. A two-dimensional retinal pattern includes multiple objects of potential interest, occluding one another and perhaps in motion, embedded within regions that may exhibit very different intensities of light and contrast. Odorant receptors are affected by concentration as well as quality; moreover, the ligands that comprise different sources of interest may compete for the same receptors—as agonists, partial agonists, or antagonists—and hence occlude source-specific signal patterns (Gronowitz et al., 2020; Xu et al., 2020; Zak et al., 2020). Indeed, even two source odorants in binary mixtures often strongly interfere with one another, generating activation patterns that do not resemble the diagnostic patterns of any individual source (Linster and Cleland, 2004; Riffell, 2012; Thomas-Danguin et al., 2014). Aside from the specialized mixture interactions literature, most contemporary studies of “odor coding” elide this fundamental problem, focusing instead on single odorants delivered under clean conditions. We argue here that the problem of inter-source interference is fundamentally constraining on early-stage olfactory encoding strategies and must be resolved if these strategies and their underlying mechanisms are to be understood.
We propose that the most effective metaphor for framing the complexities of sensory processing is that neuronal populations comprise computational stages, each of which selectively constructs a sensory representation from the activity of its antecedents. That is, after the most peripheral sensory circuits have sampled a wide diversity of information from the environment, secondary sensory circuits sample in turn from this ensemble, deploying specific sampling strategies to selectively extract desired information as if they were themselves primary sensory neurons. Critically, multiple different neuronal populations can extract different information from the same antecedent ensemble, and these divergent pathways then can serve different, if potentially interacting, goals (Milner and Goodale, 2008; Milner, 2017). In the olfactory system, for example, mitral cells, and projecting tufted cells exhibit different response profiles to odors and concentrations (Igarashi et al., 2012; Geramita and Urban, 2017), although both sample from the same primary sensory information presented by MOB glomeruli. Moreover, different subsets of these principal neuron populations project from the olfactory bulb to divergent follower structures—including a lateral pathway to the piriform cortex (PCx) and cortically organized olfactory tubercle, a ventromedial pathway to the striatally organized portion of the olfactory tubercle, and a rostromedial pathway to the tenia tecta and indusium griseum (Ojima et al., 1984; Cleland and Linster,
These selective sampling processes are inescapably intertwined with the transformations in the representational metric itself that occur in the MOB. Activity across the population of primary sensory neurons is temporally structured in the theta band (2–8 Hz; owing primarily to respiration), but there is no suggestion there of spike timing regulation on any faster timescale. Hence, when considered at faster timescales this information is effectively rate-coded (with high OSN convergence ratios onto MOB glomeruli greatly reducing the time required for a reliable estimate of population activity). In contrast, activity in MTCs is structured on the gamma timescale (40–110 Hz). Specifically, mitral cell spikes are phase-constrained with respect to stimulus-evoked gamma rhythmicity (Eeckman and Freeman, 1990; Kashiwadani et al., 1999; Bathellier et al.,
For example, current models of olfactory bulb information processing (Cleland and Sethupathy, 2006; Cleland,
Transforming Variance of Interest
A third essential determinant of sensory processing mechanics is their need for integration with other sensory modalities and brain systems. Both the information architecture that organizes communication among disparate brain systems and the underlying physical mechanisms that directly mediate this communication must be considered when hypothesizing encoding and transformation strategies for particular systems. Two related questions, in particular, should be posed: (1) is there physiological evidence for a common abstract metric space in which multiple modalities or antecedent ensembles can manipulate shared variables? (2) is multimodal integration a discrete, discernable stage of processing, or might it be intercalated deeply into early sensory processing? These questions both relate to the overall question of just how modular and functionally independent we should consider such neural systems to be. How much do local networks insulate one another from a dependence on their many essential local parameters? For example, primary motor cortex in primates appears to encode hand movements as population vectors in veridical space, insulating this common, abstract representational metric from the direct governance of the contractions of many obliquely oriented muscle fibers that control the angles of multiple arm joints and the need to adapt to their non-linear viscoelastic properties (Georgopoulos, 1996; Georgopoulos and Carpenter, 2015). Higher-level olfactory representation often is considered more categorical, hence unlikely to utilize an underlying low-dimensional map akin to these population vectors, but the underlying question remains. How is odor information embedded into useful multimodal representations, and where and how does multimodal information merge with the olfactory sensory cascade?
Odor identification is commonly considered the chief goal of the olfactory system, but it is not clear that this constitutes a singular representational end point. One might instead ask, what are the implications of the sampled olfactory information? A particular odor might have different implications in different contexts, or when paired with specific inputs from other modalities. Should the sensory system be expected to purify the odor identification signal first, and only then coordinate it with other modalities and memories (the classical sensory cortex/association cortex dichotomy)? Or do physiological odor representations depend upon non-olfactory factors at early stages of processing, such that they quickly cease to be strictly “odor representations?”
There is increasing evidence that the latter is the case. In particular, recent work in the anterior olfactory nucleus (AON) has revealed that it is, at least in part, an efferent structure that mediates contextual information from the hippocampus into the olfactory system—very likely into the olfactory bulb itself (Aqrabawi et al.,
The governance of interareal communication in the brain is often based upon transient periods of coherence in higher frequency bands including beta (15–40 Hz) and gamma (Baird and Eeckman,
Theoretical Capacities of Olfactory Bulb Circuits
Assessing the capabilities of a complex circuit requires a heavy reliance on theory, and specifically theory that can be challenged and validated across multiple scales and levels of organization. Our early network models examined how OSN populations are able to encode odorants across wide concentration ranges despite the sharply limited dynamic ranges of individual OSNs (Cleland and Linster,
In contrast, other olfactory systems models are genuinely multiscale, combining algorithmic function with Marr's lowest level of implementational detail. These compartmental models require copious constraining data for validation, but are uniquely valuable as existence proofs for the operations of complex physiological systems. For example, fast, synchronous network oscillations in the gamma band have long been recognized as endogenous to the activated olfactory bulb, emerging from the reciprocal synaptic interactions of the external plexiform layer (EPL; Figure 1; Lagier et al., 2007). However, the dynamical system presumed to underlie these oscillations, pyramidal-interneuron network gamma (PING; Traub et al., 1997; Borgers et al.,
One limitation of both types of purely computational models is that they are only able to study the problems that are appreciated by their creators. Sensory systems deployed “in the wild” have many additional, underappreciated sources of variance to contend with: e.g., interference among odorants from multiple sources (competition for receptors, unknown ligand-receptor efficacies), the adaptation, decay, drift, and/or dropout of receptors, unpredictable fluctuations in odorant concentrations, effects of temperature and humidity, and odor plume dynamics. Testing models with data from physical chemosensor arrays and instantiating biomimetic networks within field-deployable systems provides yet another level of insight into the inconspicuous but essential operations that such systems must perform. Notably, these engineered implementations have revealed a number of problems encountered by our algorithms when deployed “in the wild,” each of which has suggested resolutions that now can be analyzed and experimentally investigated in the biological system. We here summarize and review our common, internally consistent, and broadly data-constrained multiscale framework for sensory processing in the two computational layers of the olfactory bulb, as instantiated in and vetted by computational models and/or engineered neuromorphic systems.
Glomerular Layer Transformations
The axons of primary OSNs project to the glomerular layer of the MOB, with those of OSNs that express the same odorant receptor converging together to form receptor-specific glomeruli within that layer (Figure 1). The activity of these convergent OSNs is sampled by tens of principal neurons (MTCs) extending apical dendrites into each glomerulus; however, the resulting activation of these MTCs is substantially modified by glomerular-layer computations mediated by a diverse population of interneurons. These computations effect critical signal conditioning functions including contrast enhancement and concentration tolerance, though evidence of structural plasticity among glomerular interneurons further suggests that this layer also mediates additional operations of interest.
Concentration Tolerance
The need for substantial concentration tolerance transformations in the early olfactory cascade, prior to MTC activation, has long been clear (Cleland et al.,
A substantial portion of the absolute compression of environmental intensities is probably mediated by a series of mechanisms including axonal convergence, feedback inhibition, and homeostatic synaptic scaling within MOB columns (Cleland et al.,
Subsequent experimental work, however, made it clear that short-axon cells are in fact inhibitory interneurons, expressing GABA and dopamine (Banerjee et al.,
Non-topographical Contrast Enhancement
Non-topographical contrast enhancement (NTCE; Cleland and Sethupathy, 2006; Cleland,
The mechanism proposed for NTCE constituted a self-normalization scheme in which the glomerular circuitry of individual MOB columns generated an “inhibitory surround” effect that was directly dependent on overall afferent activation levels. Briefly, when afferent input was weak to moderate, feedforward inhibition would dominate and the MTC would be inhibited, whereas when the input to that column was stronger, excitation would overcome inhibition and drive MTCs to fire. A relational normalization circuit (see Concentration Tolerance section) was required in order to render the NTCE mechanism robust to changing odorant concentrations, thereby ensuring that “weak/moderate” and “stronger” inputs in this context would be assessed relative to one another. That is, across an activated ensemble of MOB columns, the MTCs of the most strongly activated columns would fire, those of the more weakly to moderately activated columns would be inhibited below baseline, and the most weakly activated columns would, naturally, generate little physiological response. This algorithm recapitulated the observed phenomenon that, when recording from one mitral cell and moving through its chemoreceptive field by presenting a series of sequentially similar odorants, the MTC would respond with excitation to its preferred odorant stimuli and with below-baseline inhibition specifically to those odorants most similar to its preferred range (i.e., its “surround”). That is, NTCE effected surround inhibition with respect to a high-dimensional odor similarity space that we now term R-space (Figure 2; Gronowitz et al., 2020). Moreover, NTCE predicted that this high-dimensional surround inhibition would be mediated by periglomerular interneurons, rather than granule cells, both because only periglomerular cells can effect feedforward intraglomerular inhibition (Figure 1) and because, even with that point aside, granule cells are not well-positioned to suppress MTC spiking activity (Mcintyre and Cleland, 2016). Subsequent experimental and computational work in which the stringency of this contrast enhancement could be regulated by established mechanisms of cholinergic neuromodulation (Mandairon et al., 2006; Chaudhury et al.,
Figure 2

Depiction of a three-dimensional R-space. The space of all possible instantaneous odorant representations at the receptor (glomerular input) level can be depicted as a unit N-cube, where each axis corresponds to the activation level of one type of odorant receptor, and N is the number of odorant receptor types expressed. Note that the numbers of receptor types expressed, as well as their precise structures, are species-specific; accordingly, odorant similarity is an organism-dependent measure. In the three-receptor system depicted, the activation levels of the three receptor types R1, R2, and R3 each independently range from zero at the origin (black dot) to unity (maximally activated). Proximity in this metric space roughly corresponds to odorant similarity (neglecting, for example, learning-dependent effects on perception). Learned odors admit ranges of tolerated variance in the activation levels of each receptor type, reflecting category learning and measurable by generalization gradients; these are depicted in R-space as odor source volumes (colored shapes). Importantly, realistic R-spaces exhibit much higher dimensionalities than can be depicted here (human R-space is ~400-dimensional, and the R-spaces of mice and rats are ~1,200-dimensional), such that the distribution of odor source volumes is (1) much sparser than depicted here, and (2) likely to lie along a surface of the unit hypercube, in that odors generally do not activate all of the receptor types in an animal's complement.
Like the concentration tolerance circuit discussed above (see Concentration Tolerance section), the NTCE mechanism also requires modification owing to important experimental updates in the structure of glomerular-layer circuitry. However, the essence of its mechanism has been supported by a broad range of subsequent experimental tests. It is now quite clear that mitral cells are gated, and delayed in their activation, by feedforward intraglomerular inhibition delivered by periglomerular cells (Gire and Schoppa, 2009; Geramita and Urban, 2017). Direct tests of the inhibitory control profile of MTCs by feedforward inhibition, including the activation-dependent masking of evoked excitation by evoked inhibition, specifically favored the NTCE mechanism (Fukunaga et al., 2014). The relational normalization process in the deep glomerular layer has been experimentally established (Banerjee et al.,
Indeed, the indirect mechanism presents advantages to the NTCE framework, alleviating some of its vulnerabilities. For example, NTCE works best when the effects of (disynaptic) feedforward inhibition onto MTCs reliably precede MTC excitation, which had been considered exclusively monosynaptic. This was plausible, in that the exceptionally high input resistance of periglomerular cells and the likelihood that feedforward inhibition was mediated through their gemmules (spines) together facilitated rapid, graded inhibition of MTCs that, if not fast enough to precede MTC excitation per se, could be fast enough to precede MTC spike initiation. However, the inhibition computed by the global normalization circuit mediated by short-axon cells would necessarily be slower, requiring the initial inhibitory response to be sustained until this laterally computed feedback inhibition could be incorporated. Mitral cell spike times were in fact sufficiently delayed so as to permit this, and the T-current plateaus observed in periglomerular cells (McQuiston and Katz, 2001) were a potential supporting mechanism for (briefly) sustained inhibition, but this was nevertheless an underspecified weakness in the model mechanism. Indirect MTC excitation via external tufted cells, however, resolves both of these issues, as it slows the process of MTC excitation and additionally fills the role of adjudicating the diverse excitatory and inhibitory inputs to the column, delivering the result of this summation to MTCs. Finally, the indirect framework also simplifies the problem of ensuring that stimulus-activated MTC spikes are robustly constrained by emergent fast oscillatory dynamics. Respiration-synchronized feedforward inhibition evokes a phase reset in mitral cell subthreshold oscillations (Desmaisons et al., 1999; Rubin and Cleland, 2006), which is theoretically important for the reliable recruitment of odorant-activated MTCs into a common gamma regime (Li and Cleland, 2017). Mitral cells that must directly adjudicate inputs from diverse excitatory and inhibitory inputs distributed in time might be less likely to enjoy a common, distinct reset phase across sibling MTCs, or even across the apical dendritic arbor within a single MTC. Whereas MTCs associated with the same glomerulus exhibit heterogeneous properties and non-identical responses, there is evidence that they are coordinated to a certain extent as a common functional unit (Schoppa and Westbrook, 2001).
External Plexiform Layer Transformations
As the previously purported functions of the EPL now were assigned to the glomerular layer, the question naturally arose: what does the EPL do? It clearly comprises a lateral inhibitory network, but does not respect physical proximity—distributions of functional lateral inhibitory efficacies are largely distance-independent (Fantana et al., 2008). Additionally, accumulating evidence indicated that the EPL exhibited sophisticated synaptic and structural plasticity mechanisms (Balu et al.,
From this, two particular hypotheses stood out. First, the construction of category representations for meaningful “odors” presumably must be initiated in the MOB, where the full spectrum of afferent variance remains computationally available. MOB output diverges into at least three separate pathways, as noted above, and recordings from anterior PCx (for example) suggest that odor representations there already exhibit categorical indicators such as experience-dependent configural responses (Wilson et al., 2006) and discontinuous receptive fields (Stettler and Axel, 2009). Second, receptor occlusion by competing odorant sources is a fundamental, catastrophic challenge to contemporary models of olfactory encoding and representation. Background odorants powerfully disrupt the integrity of olfactory ensemble representations, and are ubiquitous in natural environments. The ability to identify simultaneously-encountered odors of interest in strong and unpredictable backgrounds is essential, and must be resolved early in the sensory cascade. Interestingly, these two hypotheses are likely to interact, as category learning can in principle provide the information needed to identify a known signal presented against a background of strongly disruptive interference.
Neuromorphic Design
To instantiate and test these hypotheses, we constructed models of the EPL network and embedded some variants of these models in the experimental Intel Loihi neuromorphic system (Davies et al., 2018). From the engineering perspective, the diagnostic feature of neuromorphic networks is a reliance on local computation; no global error term is assessed, and the weight of a synapse is modifiable only by the activity experienced directly by that synapse. Neuromorphic algorithms also tend to rely on spike-based communication (i.e., they operate as pulsed networks; Maass and Bishop, 1998), though this is not necessarily diagnostic. Finally, they often rely upon task-specific network architectures, and hence can be much faster to train at the cost of generality (Imam and Cleland, 2020). This combination of features enables neuromorphic hardware platforms to execute appropriate neuromorphic algorithms rapidly and with very low energy expenditure, notably by colocalizing memory and compute resources on-chip to leverage the benefits of local computation. A still-underappreciated advantage of neuromorphic design is the ease with which the properties of nodes (cells) and edges (synapses) can be arbitrarily heterogeneous and elaborated, and networks can fork, merge, and loop datastreams. These features, however, are difficult to optimize, and hence place substantial burdens upon the designer. Biologically inspired network motifs fill this gap admirably, as the biological systems—to the extent that they are correctly understood—provide existence proofs for the functionality and robustness of the network architecture. Conversely, and with the same caveats, physically instantiated biomimetic architectures tested with real-world datasets can challenge our understandings of biological systems and present quantitatively vetted hypotheses for the next iteration of conceptual frameworks.
EPL Network Plasticity Enables Rapid Learning and Odor Identification Under Noise
We constructed a network model of the MOB external plexiform layer based on its theoretical capacities and key circuit motifs in order to address the problems of (1) how categorical odor representations are generated and (2) how these might then be leveraged to identify known odor sources within highly interfering and unpredictable backgrounds (Imam and Cleland, 2020). Briefly, if the plastic synaptic matrix of the olfactory bulb EPL (Figure 1) can be leveraged to construct categorical odor representations, sparsely distributed within a high-dimensional receptor-based similarity space (R-space; Figure 2; Gronowitz et al., 2020), then an attractor network (Hopfield, 1982; Kay et al., 1996) should be capable of identifying a known odor despite this disruptive interference. Notably, the high dimensionality of this olfactory R-space is of great benefit, both because the greater number of receptor types increases the likelihood that some critical fraction will remain free of catastrophic occlusion, and because higher degrees of freedom means that random disruptive interference is increasingly unlikely to resemble a different learned odor representation.
To construct this network model, we instantiated some core principles of MOB computation that we have hypothesized for the biological system (Cleland,
The importance of gamma-discretized spike timing-based computation in the EPL. This enables the direct deployment of asymmetric spike timing-dependent plasticity (STDP) rules, which in turn enable the construction of higher-order, configural receptive fields in granule cells (Linster and Cleland, 2010). PCx pyramidal neurons also are known to be selective for synchronized MTC spikes on this timescale (Luna and Schoppa, 2008).
The principle that that only a minority of principal neurons participate in gamma dynamics during any given stimulus presentation, based on their afferent activation (Li and Cleland, 2017). MTCs are known for their relatively high background spiking activity (a side effect of their high sensitivity to inputs), but this background activity does not correspond to MTCs' responsivity to odor stimuli (Kollo et al., 2014), nor do background spikes align to the common, coherent gamma rhythm (Werth and Cleland, 2019). The synchrony-selective mechanisms of point (1) suggest that MTC activity that is disorganized in time simply fails to effectively activate follower neurons.
The principle that GC inhibition of MTCs manifests as delays in MTC spike times on the gamma timescale (inhibition-imposed phase lag). This counteracts the phase lead exhibited by more strongly activated MTCs owing to their more rapid escape from inhibition on each gamma cycle, thereby recapitulating the direct opposition of cellular excitation and inhibition within a gamma-scale spike timing regime. This metric transformation neatly avoids the catastrophic limitation of rate coding assumptions in sparsely active networks like the MOB: the need for postsynaptic neurons to integrate over extended periods of time in order to accurately estimate presynaptic activity.
The principle that the topology of the EPL network is not reciprocal, despite observations that individual synaptic connections between GCs and MTC lateral dendrites often are bidirectional. Specifically, MTC lateral dendrites project broadly across the EPL (Mori et al., 1983; Orona et al., 1984) and support regenerative spiking, and hence can deliver excitation to GCs independent of proximity (Xiong and Chen, 2002). In contrast, inhibition does not propagate; for GCs to affect spike timing propagation in MTCs, they must deliver their shunting inhibition very close to the afferent axis (i.e., the MTC soma; Mcintyre and Cleland, 2016). This facilitates the abstraction of MOB columns, each consisting of a glomerulus, its associated principal neurons, and the GCs that deliver effective inhibition onto those principal neurons. Transsynaptic tracing studies suggest just such a discontinuous and column-centric topology of MTC-GC connectivity (Willhite et al., 2006; Kim et al., 2011).
The principle that, within the constraints of this asymmetry, the topology of the lateral inhibitory network of the EPL is dynamically acquired and learning-dependent. Specifically, GCs learn higher-order configural receptive fields, becoming more specifically diagnostic of particular odorants than is possible for principal neurons, and then deploy this information to regulate MTC spike timing. Consistent with this principle, GCs require either coincident activity by a large number of MTCs, or by a smaller number of MTCs coupled with excitation from PCx, in order to fire reliably (Pressler and Strowbridge, 2017). Moreover, task learning is known to shape MTC activation profiles (Doucette and Restrepo, 2008), and acquired odor-reward associations regulate MTC spike synchronization patterns (Doucette et al., 2011).
The corollary to point (5) that GCs are permanently differentiated by recruitment during the process of odor learning, and therefore require constitutive replacement by adult neurogenesis in order to enable future (lifelong) learning. Such a strategy for category learning may explain why GCs outnumber MTCs by roughly 25:1 (Meisami and Safari, 1981; Frazier and Brunjes, 1988; Royet et al., 1998), whereas most brain regions exhibit comparable numbers of excitatory and inhibitory neurons, or even a bias toward the excitatory (Shepherd, 2004).
The possibility that neuromodulation should be conceived as an optimization trajectory rather than as a stationary, behaviorally-dependent state.
This combination of elemental structures and properties enabled the construction of a neuromorphic network that, for the first time, successfully addressed the problem of identifying odor sources despite powerful interference by simultaneously encountered background odorants within a complex chemical environment (Imam and Cleland, 2020). Odor presentations were drawn from the responses of a 72-chemosensor array to stimulation with a complex odorant plume in an wind tunnel (Vergara et al., 2013). Background interference included both Bernoulli and quasi-Gaussian components. Specifically, Bernoulli noise was delivered by replacing the activity levels of some fraction of the sensors (typically 40–60%) with random numbers. This process models competition for receptors by other, unknown odorant species that may strongly activate, weakly occupy, or block the activation of any given receptor in an unpredictable manner (Araneda et al.,
Training using this algorithm is rapid—even one-shot learning suffices to learn a new odorant signature. Moreover, the network naturally exhibits online learning—the capacity to learn multiple inputs (odors) sequentially, without the new training causing the catastrophic forgetting of previous training. Both of these properties stand in sharp contrast to the properties of contemporary deep networks. Specifically, a similarly sized deep network required thousands of training trials, plus foreknowledge of the variance structure of the interference, in order to reach the same levels of task performance as this neuromorphic strategy, and also lacked the capacity for online learning (Imam and Cleland, 2020). These properties of the EPL network model parallel those of the biological system: rapid learning, robust identification under noise, and the potential for sustained performance across unpredictable or changing environments.
This implementation demonstrates that the circuitry of the mammalian MOB is, in principle, capable of supporting odor category learning (potentially inclusive of context and valence, as discussed above) and of using these categorical memories to facilitate odor recognition in complex milieus. The particular benefit is that category learning here is not simply a cognitive outcome, but also can be applied very early in the sensory cascade in order to achieve an extraordinarily difficult sensory goal. That said, there are critical simplifications in this initial implementation that must be superseded with further development. Behavioral studies of learned odor generalization (Cleland et al., 2009,
Learning in the Wild
The initial EPL model proved robust to the amount of variance associated with the odor plume together with a measure of background occlusion. However, as discussed above, biological systems are subject to many additional sources of uncontrolled variance. Concentrations vary widely, sensors decay, and the network must continually update representations according to newly acquired information. Moreover, the goals of odor learning include not simply classification into one of a set of discrete categories, but the recognition of similarities among diverse stimuli sufficient to form generalization gradients (Shepard, 1987; Cleland et al.,
Some aspects of learning in the wild simply require the inclusion of established glomerular-layer computations affording concentration tolerance and regulated contrast enhancement (Imam et al., 2012; Borthakur and Cleland,
Heterogeneity in cellular and synaptic properties is of course expected in biological systems, but there is increasing evidence that this is a cultivated property of neuronal networks rather than arising solely from inescapable random variability. Indeed, heterogeneity in cellular and synaptic properties has been credited with improving the sensitivity, efficiency, and information content of neuronal ensembles (Lengler et al., 2013; Zohar et al., 2013) and regulating the balance between network robustness and flexibility (Gu et al., 2019). The advantages of parameter heterogeneity observed in artificial systems models are indicating a broader emerging principle of biomimetic design: considerable computational power and robustness may be achievable by relaxing control over critical variables and subsequently reconstructing the required specificity through effective sampling strategies.
Conclusions
Understanding any information processing system requires insight into the problems that it faces, and the forms that possible solutions to these problems can take (Marr and Poggio, 1977). Hence, the role of theoretical neuroscience is not only to replicate and integrate established results, but to extend beyond contemporary experimental findings so as to inspire and guide subsequent investigations. Engineered systems models expand our capacity to identify complex theoretical problems and limitations that neither neurobiological experiments nor more narrowly targeted models are able to appreciate. Such models can be speculative, but are at their most powerful when they analytically cross-reference information from multiple sources, constraining the possible interpretations more sharply than any of the individual contributing studies are able to do. This is doubly true when such models incorporate complex, multi-stage systems within a single, internally consistent, quantitative framework. Some problems are only addressed by sequences of mechanisms coordinated across stages. Some other problems are best solved by relaxing the diagnostic effort at one stage so as to improve the capacity of a subsequent stage to achieve that diagnostic goal in a flexible and task-dependent manner. Biological networks exemplify this multistage, systems approach to sensory sampling and information acquisition. A corollary of this position is that narrower, stage-specific analyses are at risk of misinterpreting the computational and “coding” strategies of the system.
We propose a common, integrated theory of olfactory systems operation that incorporates experimental findings across scales, stages, and methods of analysis into a common framework. This constructive approach to systems analysis considers individual stages of processing in a broadened context as contributors to global goals, while still generating concrete, stage-specific hypotheses that can be falsified or updated by experimental results. By leveraging the test platforms and analytical capabilities of engineered approaches together with our expanding capacity to interrogate biological systems, contemporary theorists can design, implement, and validate increasingly sophisticated and integrated theoretical models of neural and cognitive systems.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.
Author contributions
TC developed the general theory and wrote the paper. AB developed and analyzed olfactory bulb models and contributed to the writing of the paper. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by NIH/NIDCD grants R01 DC014701 and R01 DC014367 to TC. AB was supported in part by an Intel Neuromorphic Research Community grant to TC.
Acknowledgments
The authors would like to thank Aparna Calambur for related modeling work on the effects of network heterogeneity, and Alec J. Mutti for generating Figure 2.
Conflict of interest
The neuromorphic algorithm described in the EPL Network Plasticity Enables Rapid Learning and Odor Identification Under Noise section is the subject of a Cornell University patent application on which both authors are listed as inventors. TC is a member of the Intel Neuromorphic Research Community and has received research funding from Intel for related work. The authors declare that this study received funding from Intel. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.
References
1
AqrabawiA. J.BrowneC. J.DargaeiZ.GarandD.KhademullahC. S.WoodinM. A.et al. (2016). Top-down modulation of olfactory-guided behaviours by the anterior olfactory nucleus pars medialis and ventral hippocampus. Nat. Commun.7:13721. 10.1038/ncomms13721
2
AqrabawiA. J.KimJ. C. (2020). Olfactory memory representations are stored in the anterior olfactory nucleus. Nat. Commun.11:1246. 10.1038/s41467-020-15032-2
3
AranedaR. C.KiniA. D.FiresteinS. (2000). The molecular receptive range of an odorant receptor. Nat. Neurosci.3, 1248–1255. 10.1038/81774
4
AungstJ. L.HeywardP. M.PucheA. C.KarnupS. V.HayarA.SzaboG.et al. (2003). Centre-surround inhibition among olfactory bulb glomeruli. Nature426, 623–629. 10.1038/nature02185
5
BairdB.EeckmanF. H. (1993). A hierarchical sensory-motor architecture for oscillating cortical area subnetworks, in Neural Systems: Analysis and Modeling, ed F. H. Eeckman (Boston, MA: Springer), 47–66.
6
BaluR.LarimerP.StrowbridgeB. W. (2004). Phasic stimuli evoke precisely timed spikes in intermittently discharging mitral cells. J. Neurophysiol.92, 743–753. 10.1152/jn.00016.2004
7
BaluR.PresslerR. T.StrowbridgeB. W. (2007). Multiple modes of synaptic excitation of olfactory bulb granule cells. J. Neurosci.27, 5621–5632. 10.1523/JNEUROSCI.4630-06.2007
8
BanerjeeA.MarbachF.AnselmiF.KohM. S.DavisM. B.Garcia Da SilvaP.et al. (2015). An interglomerular circuit gates glomerular output and implements gain control in the mouse olfactory bulb. Neuron87, 193–207. 10.1016/j.neuron.2015.06.019
9
BaoX.RaguetL. L.ColeS. M.HowardJ. D.GottfriedJ. (2016). The role of piriform associative connections in odor categorization. Elife5:e13732. 10.7554/eLife.13732.012
10
BathellierB.LagierS.FaureP.LledoP. M. (2006). Circuit properties generating gamma oscillations in a network model of the olfactory bulb. J. Neurophysiol.95, 2678–2691. 10.1152/jn.01141.2005
11
BirchJ.WrightW. D. (1961). Colour discrimination. Phys. Med. Biol.6, 3–24. 10.1088/0031-9155/6/1/201
12
BoldingK. A.FranksK. M. (2017). Complementary codes for odor identity and intensity in olfactory cortex. Elife6:e22630. 10.7554/eLife.22630.023
13
BoldingK. A.FranksK. M. (2018). Recurrent cortical circuits implement concentration-invariant odor coding. Science361:eaat6904. 10.1126/science.aat6904
14
BorgersC.EpsteinS.KopellN. J. (2005). Background gamma rhythmicity and attention in cortical local circuits: a computational study. Proc. Natl. Acad. Sci. U.S.A.102, 7002–7007. 10.1073/pnas.0502366102
15
BorthakurA.ClelandT. A. (2019a). Signal conditioning for learning in the wild, in Neuro-Inspired Computational Elements Workshop (Albany, NY: Association for Computing Machinery (ACM)). 10.1145/3320288.3320293
16
BorthakurA.ClelandT. A. (2019b). A spike time-dependent online learning algorithm derived from biological olfaction. Front. Neurosci.13:656. 10.3389/fnins.2019.00656
17
BreaJ. N.KayL. M.KopellN. J. (2009). Biophysical model for gamma rhythms in the olfactory bulb via subthreshold oscillations. Proc. Natl. Acad. Sci. U.S.A.106, 21954–21959. 10.1073/pnas.0910964106
18
BurtonS. D. (2017). Inhibitory circuits of the mammalian main olfactory bulb. J. Neurophysiol.118, 2034–2051. 10.1152/jn.00109.2017
19
CenierT.DavidF.LitaudonP.GarciaS.AmatC.BuonvisoN. (2009). Respiration-gated formation of gamma and beta neural assemblies in the mammalian olfactory bulb. Eur. J. Neurosci.29, 921–930. 10.1111/j.1460-9568.2009.06651.x
20
ChaudhuryD.EscanillaO.LinsterC. (2009). Bulbar acetylcholine enhances neural and perceptual odor discrimination. J. Neurosci.29, 52–60. 10.1523/JNEUROSCI.4036-08.2009
21
ClapperJ. P. (2019). Graded similarity in free categorization. Cognition190, 1–19. 10.1016/j.cognition.2019.04.009
22
ClelandT. A. (2010). Early transformations in odor representation. Trends Neurosci.33, 130–139. 10.1016/j.tins.2009.12.004
23
ClelandT. A. (2014). Construction of odor representations by olfactory bulb microcircuits. Prog. Brain Res.208, 177–203. 10.1016/B978-0-444-63350-7.00007-3
24
ClelandT. A.ChenS. Y.HozerK. W.UkatuH. N.WongK. J.ZhengF. (2012). Sequential mechanisms underlying concentration invariance in biological olfaction. Front. Neuroeng.4:21. 10.3389/fneng.2011.00021
25
ClelandT. A.JohnsonB. A.LeonM.LinsterC. (2007). Relational representation in the olfactory system. Proc. Natl. Acad. Sci. U.S.A.104, 1953–1958. 10.1073/pnas.0608564104
26
ClelandT. A.LinsterC. (1999). Concentration tuning mediated by spare receptor capacity in olfactory sensory neurons: a theoretical study. Neural Comput.11, 1673–1690. 10.1162/089976699300016188
27
ClelandT. A.LinsterC. (2005). Computation in the olfactory system. Chem. Senses30, 801–813. 10.1093/chemse/bji072
28
ClelandT. A.LinsterC. (2019). Central olfactory structures, in Handbook of Clinical Neurology, Volume 164: Smell and Taste, eds DotyR. L.Leon-SarmientoF. E. (New York, NY: Elsevier), 79–96.
29
ClelandT. A.MorseA.YueE. L.LinsterC. (2002). Behavioral models of odor similarity. Behav. Neurosci.116, 222–231. 10.1037/0735-7044.116.2.222
30
ClelandT. A.NarlaV. A. (2003). Intensity modulation of olfactory acuity. Behav. Neurosci.117, 1434–1440. 10.1037/0735-7044.117.6.1434
31
ClelandT. A.NarlaV. A.BoudadiK. (2009). Multiple learning parameters differentially regulate olfactory generalization. Behav. Neurosci.123, 26–35. 10.1037/a0013991
32
ClelandT. A.SethupathyP. (2006). Non-topographical contrast enhancement in the olfactory bulb. BMC Neurosci.7:7. 10.1186/1471-2202-7-7
33
DaviesM.SrinivasaN.LinT.-H.ChinyaG.CaoY.ChodayS. H.et al. (2018). Loihi: a neuromorphic manycore processor with on-chip learning. IEEE Micro38, 82–99. 10.1109/MM.2018.112130359
34
DesmaisonsD.VincentJ. D.LledoP. M. (1999). Control of action potential timing by intrinsic subthreshold oscillations in olfactory bulb output neurons. J. Neurosci.19, 10727–10737. 10.1523/JNEUROSCI.19-24-10727.1999
35
DoucetteW.GireD. H.WhitesellJ.CarmeanV.LuceroM. T.RestrepoD. (2011). Associative cortex features in the first olfactory brain relay station. Neuron69, 1176–1187. 10.1016/j.neuron.2011.02.024
36
DoucetteW.RestrepoD. (2008). Profound context-dependent plasticity of mitral cell responses in olfactory bulb. PLoS Biol.6:e258. 10.1371/journal.pbio.0060258
37
EdelmanS.ShahbaziR. (2012). Renewing the respect for similarity. Front. Comput. Neurosci.6:45. 10.3389/fncom.2012.00045
38
EeckmanF. H.FreemanW. J. (1990). Correlations between unit firing and EEG in the rat olfactory system. Brain Res.528, 238–244. 10.1016/0006-8993(90)91663-2
39
FantanaA. L.SoucyE. R.MeisterM. (2008). Rat olfactory bulb mitral cells receive sparse glomerular inputs. Neuron59, 802–814. 10.1016/j.neuron.2008.07.039
40
FrazierL. L.BrunjesP. C. (1988). Unilateral odor deprivation: early postnatal changes in olfactory bulb cell density and number. J. Comp. Neurol.269, 355–370. 10.1002/cne.902690304
41
FrederickD. E.BrownA.BrimE.MehtaN.VujovicM.KayL. M. (2016). Gamma and beta oscillations define a sequence of neurocognitive modes present in odor processing. J. Neurosci.36, 7750–7767. 10.1523/JNEUROSCI.0569-16.2016
42
FriesP. (2015). Rhythms for cognition: communication through coherence. Neuron88, 220–235. 10.1016/j.neuron.2015.09.034
43
FukunagaI.HerbJ. T.KolloM.BoydenE. S.SchaeferA. T. (2014). Independent control of gamma and theta activity by distinct interneuron networks in the olfactory bulb. Nat. Neurosci.17, 1208–1216. 10.1038/nn.3760
44
GallianoE.FranzoniE.BretonM.ChandA. N.ByrneD. J.MurthyV. N.et al. (2018). Embryonic and postnatal neurogenesis produce functionally distinct subclasses of dopaminergic neuron. Elife7:e32373. 10.7554/eLife.32373.022
45
GaoY.StrowbridgeB. W. (2009). Long-term plasticity of excitatory inputs to granule cells in the rat olfactory bulb. Nat. Neurosci.12, 731–733. 10.1038/nn.2319
46
GeorgopoulosA. P. (1996). Arm movements in monkeys: behavior and neurophysiology. J. Comp. Physiol. A179, 603–612. 10.1007/BF00216125
47
GeorgopoulosA. P.CarpenterA. F. (2015). Coding of movements in the motor cortex. Curr. Opin. Neurobiol.33, 34–39. 10.1016/j.conb.2015.01.012
48
GeramitaM.UrbanN. N. (2017). Differences in glomerular-layer-mediated feedforward inhibition onto mitral and tufted cells lead to distinct modes of intensity coding. J. Neurosci.37, 1428–1438. 10.1523/JNEUROSCI.2245-16.2016
49
GeramitaM. A.BurtonS. D.UrbanN. N. (2016). Distinct lateral inhibitory circuits drive parallel processing of sensory information in the mammalian olfactory bulb. Elife5:e16039. 10.7554/eLife.16039.022
50
GireD. H.FranksK. M.ZakJ. D.TanakaK. F.WhitesellJ. D.MulliganA. A.et al. (2012). Mitral cells in the olfactory bulb are mainly excited through a multistep signaling path. J. Neurosci.32, 2964–2975. 10.1523/JNEUROSCI.5580-11.2012
51
GireD. H.SchoppaN. E. (2009). Control of on/off glomerular signaling by a local GABAergic microcircuit in the olfactory bulb. J. Neurosci.29, 13454–13464. 10.1523/JNEUROSCI.2368-09.2009
52
GronowitzM. E.LiuA.ClelandT. A. (2020). A physicochemical model of odor sampling. bioRxiv [Preprint]10.1101/2020.06.16.154807
53
GuC. G.WangP.WengT. F.YangH. J.RohlingJ. (2019). Heterogeneity of neuronal properties determines the collective behavior of the neurons in the suprachiasmatic nucleus. Math. Biosci. Eng.16, 1893–1913. 10.3934/mbe.2019092
54
HackM. A.SaghatelyanA.De ChevignyA.PfeiferA.Ashery-PadanR.LledoP. M.et al. (2005). Neuronal fate determinants of adult olfactory bulb neurogenesis. Nat. Neurosci.8, 865–872. 10.1038/nn1479
55
HopfieldJ. J. (1982). Neural networks and physical systems with emergent collective computational abilities. Proc. Natl. Acad. Sci. U.S.A.79, 2554–2558. 10.1073/pnas.79.8.2554
56
IgarashiK. M.IekiN.AnM.YamaguchiY.NagayamaS.KobayakawaK.et al. (2012). Parallel mitral and tufted cell pathways route distinct odor information to different targets in the olfactory cortex. J. Neurosci.32, 7970–7985. 10.1523/JNEUROSCI.0154-12.2012
57
ImamN.ClelandT. A. (2020). Rapid learning and robust recall in a neuromorphic olfactory circuit. Nat. Mach. Intell.2, 181–191. 10.1038/s42256-020-0159-4
58
ImamN.ClelandT. A.ManoharR.MerollaP. A.ArthurJ. V.AkopyanF.et al. (2012). Implementation of olfactory bulb glomerular-layer computations in a digital neurosynaptic core. Front. Neurosci.6:83. 10.3389/fnins.2012.00083
59
KashiwadaniH.SasakiY. F.UchidaN.MoriK. (1999). Synchronized oscillatory discharges of mitral/tufted cells with different molecular receptive ranges in the rabbit olfactory bulb. J. Neurophysiol.82, 1786–1792. 10.1152/jn.1999.82.4.1786
60
KayL. M.BeshelJ. (2010). A beta oscillation network in the rat olfactory system during a 2-alternative choice odor discrimination task. J. Neurophysiol.104, 829–839. 10.1152/jn.00166.2010
61
KayL. M.LancasterL. R.FreemanW. J. (1996). Reafference and attractors in the olfactory system during odor recognition. Int. J. Neural. Syst.7, 489–495. 10.1142/S0129065796000476
62
KimD. H.PhillipsM. E.ChangA. Y.PatelH. K.NguyenK. T.WillhiteD. C. (2011). Lateral connectivity in the olfactory bulb is sparse and segregated. Front. Neural. Circuits5:5. 10.3389/fncir.2011.00005
63
KiyokageE.PanY. Z.ShaoZ.KobayashiK.SzaboG.YanagawaY.et al. (2010). Molecular identity of periglomerular and short axon cells. J. Neurosci.30, 1185–1196. 10.1523/JNEUROSCI.3497-09.2010
64
KolloM.SchmaltzA.AbdelhamidM.FukunagaI.SchaeferA. T. (2014). 'Silent' mitral cells dominate odor responses in the olfactory bulb of awake mice. Nat. Neurosci.17, 1313–1315. 10.1038/nn.3768
65
LagierS.PanzanelliP.RussoR. E.NissantA.BathellierB.Sassoe-PognettoM.et al. (2007). GABAergic inhibition at dendrodendritic synapses tunes gamma oscillations in the olfactory bulb. Proc. Natl. Acad. Sci. U.S.A.104, 7259–7264. 10.1073/pnas.0701846104
66
LenglerJ.JugF.StegerA. (2013). Reliable neuronal systems: the importance of heterogeneity. PLoS ONE8:e80694. 10.1371/journal.pone.0080694
67
LepousezG.ValleyM. T.LledoP. M. (2013). The impact of adult neurogenesis on olfactory bulb circuits and computations. Annu. Rev. Physiol.75, 339–363. 10.1146/annurev-physiol-030212-183731
68
LevinsonM.KolendaJ. P.AlexandrouG. J.EscanillaO.ClelandT. A.SmithD. M.et al. (2020). Context-dependent odor learning requires the anterior olfactory nucleus. Behav. Neurosci.134, 332–343. 10.1037/bne0000371
69
LiG.ClelandT. A. (2013). A two-layer biophysical model of cholinergic neuromodulation in olfactory bulb. J. Neurosci.33, 3037–3058. 10.1523/JNEUROSCI.2831-12.2013
70
LiG.ClelandT. A. (2017). A coupled-oscillator model of olfactory bulb gamma oscillations. PLoS Comput. Biol.13:e1005760. 10.1371/journal.pcbi.1005760
71
LinsterC.ClelandT. A. (2004). Configurational and elemental odor mixture perception can arise from local inhibition. J. Comput. Neurosci.16, 39–47. 10.1023/B:JCNS.0000004840.87570.2e
72
LinsterC.ClelandT. A. (2010). Decorrelation of odor representations via spike timing-dependent plasticity. Front. Comput. Neurosci.4:157. 10.3389/fncom.2010.00157
73
LitaudonP.AmatC.BertrandB.VigourouxM.BuonvisoN. (2003). Piriform cortex functional heterogeneity revealed by cellular responses to odours. Eur. J. Neurosci.17, 2457–2461. 10.1046/j.1460-9568.2003.02654.x
74
LitaudonP.GarciaS.BuonvisoN. (2008). Strong coupling between pyramidal cell activity and network oscillations in the olfactory cortex. Neuroscience156, 781–787. 10.1016/j.neuroscience.2008.07.077
75
LocatelliF. F.FernandezP. C.SmithB. H. (2016). Learning about natural variation of odor mixtures enhances categorization in early olfactory processing. J. Exp. Biol.219, 2752–2762. 10.1242/jeb.141465
76
LunaV. M.SchoppaN. E. (2008). GABAergic circuits control input-spike coupling in the piriform cortex. J. Neurosci.28, 8851–8859. 10.1523/JNEUROSCI.2385-08.2008
77
MaassW.BishopC. M. (1998). Pulsed Neural Networks.Cambridge, MA: MIT Press.
78
ManabeH.MoriK. (2013). Sniff rhythm-paced fast and slow gamma-oscillations in the olfactory bulb: relation to tufted and mitral cells and behavioral states. J. Neurophysiol.110, 1593–1599. 10.1152/jn.00379.2013
79
MandaironN.FerrettiC. J.StackC. M.RubinD. B.ClelandT. A.LinsterC. (2006). Cholinergic modulation in the olfactory bulb influences spontaneous olfactory discrimination in adult rats. Eur. J. Neurosci.24, 3234–3244. 10.1111/j.1460-9568.2006.05212.x
80
MarrD. (1982). Vision: a Computational Approach.San Francisco, CA: Freeman & Co.
81
MarrD.PoggioT. (1977). From understanding computation to understanding neural circuitry. Neurosci. Res. Prog. Bull.15, 470–488.
82
McintyreA. B.ClelandT. A. (2016). Biophysical constraints on lateral inhibition in the olfactory bulb. J. Neurophysiol.115, 2937–2949. 10.1152/jn.00671.2015
83
McQuistonA. R.KatzL. C. (2001). Electrophysiology of interneurons in the glomerular layer of the rat olfactory bulb. J. Neurophysiol.86, 1899–1907. 10.1152/jn.2001.86.4.1899
84
MeisamiE.SafariL. (1981). A quantitative study of the effects of early unilateral olfactory deprivation on the number and distribution of mitral and tufted cells and of glomeruli in the rat olfactory bulb. Brain Res.221, 81–107. 10.1016/0006-8993(81)91065-9
85
MilnerA. D. (2017). How do the two visual streams interact with each other?Exp. Brain Res.235, 1297–1308. 10.1007/s00221-017-4917-4
86
MilnerA. D.GoodaleM. A. (2008). Two visual systems re-viewed. Neuropsychologia46, 774–785. 10.1016/j.neuropsychologia.2007.10.005
87
MorenoM. M.LinsterC.EscanillaO.SacquetJ.DidierA.MandaironN. (2009). Olfactory perceptual learning requires adult neurogenesis. Proc. Natl. Acad. Sci. U.S.A.106, 17980–17985. 10.1073/pnas.0907063106
88
MoriK.KishiK.OjimaH. (1983). Distribution of dendrites of mitral, displaced mitral, tufted, and granule cells in the rabbit olfactory bulb. J. Comp. Neurol.219, 339–355. 10.1002/cne.902190308
89
NagayamaS.EnervaA.FletcherM. L.MasurkarA. V.IgarashiK. M.MoriK.et al. (2010). Differential axonal projection of mitral and tufted cells in the mouse main olfactory system. Front. Neural. Circuits4:120. 10.3389/fncir.2010.00120
90
NajacM.De Saint JanD.RegueroL.GrandesP.CharpakS. (2011). Monosynaptic and polysynaptic feed-forward inputs to mitral cells from olfactory sensory neurons. J. Neurosci.31, 8722–8729. 10.1523/JNEUROSCI.0527-11.2011
91
Nunez-ParraA.LiA.RestrepoD. (2014). Coding odor identity and odor value in awake rodents. Prog. Brain Res.208, 205–222. 10.1016/B978-0-444-63350-7.00008-5
92
OjimaH.MoriK.KishiK. (1984). The trajectory of mitral cell axons in the rabbit olfactory cortex revealed by intracellular HRP injection. J. Comp. Neurol.230, 77–87. 10.1002/cne.902300107
93
OronaE.RainerE. C.ScottJ. W. (1984). Dendritic and axonal organization of mitral and tufted cells in the rat olfactory bulb. J. Comp. Neurol.226, 346–356. 10.1002/cne.902260305
94
OronaE.ScottJ. W.RainerE. C. (1983). Different granule cell populations innervate superficial and deep regions of the external plexiform layer in rat olfactory bulb. J. Comp. Neurol.217, 227–237. 10.1002/cne.902170209
95
PadmanabhanK.OsakadaF.TarabrinaA.KizerE.CallawayE. M.GageF. H.et al. (2018). Centrifugal inputs to the main olfactory bulb revealed through whole brain circuit-mapping. Front. Neuroanat.12:115. 10.3389/fnana.2018.00115
96
PeaceS. T.JohnsonB. C.LiG.KaiserM. E.FukunagaI.SchaeferA. T.et al. (2018). Coherent olfactory bulb gamma oscillations arise from coupling independent columnar oscillators. bioRxiv [Preprint]10.1101/213827
97
PresslerR. T.StrowbridgeB. W. (2017). Direct recording of dendrodendritic excitation in the olfactory bulb: divergent properties of local and external glutamatergic inputs govern synaptic integration in granule cells. J. Neurosci.37, 11774–11788. 10.1523/JNEUROSCI.2033-17.2017
98
QuL. P.KahntT.ColeS. M.GottfriedJ. A. (2016). De novo emergence of odor category representations in the human brain. J. Neurosci.36, 468–478. 10.1523/JNEUROSCI.3248-15.2016
99
RiffellJ. A. (2012). Olfactory ecology and the processing of complex mixtures. Curr. Opin. Neurobiol.22, 236–242. 10.1016/j.conb.2012.02.013
100
RoyetJ. P.DistelH.HudsonR.GervaisR. (1998). A re-estimation of the number of glomeruli and mitral cells in the olfactory bulb of rabbit. Brain Res.788, 35–42. 10.1016/S0006-8993(97)01504-7
101
RubinD. B.ClelandT. A. (2006). Dynamical mechanisms of odor processing in olfactory bulb mitral cells. J. Neurophysiol.96, 555–568. 10.1152/jn.00264.2006
102
SchaeferM. L.YoungD. A.RestrepoD. (2001). Olfactory fingerprints for major histocompatibility complex-determined body odors. J. Neurosci.21, 2481–2487. 10.1523/JNEUROSCI.21-07-02481.2001
103
SchoppaN. E.WestbrookG. L. (2001). Glomerulus-specific synchronization of mitral cells in the olfactory bulb. Neuron31, 639–651. 10.1016/S0896-6273(01)00389-0
104
SethupathyP.RubinD. B.LiG.ClelandT. A. (2013). A model of electrophysiological heterogeneity in periglomerular cells. Front. Comput. Neurosci.7:49. 10.3389/fncom.2013.00049
105
ShepardR. N. (1987). Toward a universal law of generalization for psychological science. Science237, 1317–1323. 10.1126/science.3629243
106
ShepherdG. M. (2004). The Synaptic Organization of the Brain.New York, NY: Oxford University Press.
107
SoucyE. R.AlbeanuD. F.FantanaA. L.MurthyV. N.MeisterM. (2009). Precision and diversity in an odor map on the olfactory bulb. Nat. Neurosci.12, 210–220. 10.1038/nn.2262
108
StettlerD. D.AxelR. (2009). Representations of odor in the piriform cortex. Neuron63, 854–864. 10.1016/j.neuron.2009.09.005
109
StoraceD. A.CohenL. B. (2017). Measuring the olfactory bulb input-output transformation reveals a contribution to the perception of odorant concentration invariance. Nat. Commun.8:81. 10.1038/s41467-017-00036-2
110
Strube-BlossM. F.RosslerW. (2018). Multimodal integration and stimulus categorization in putative mushroom body output neurons of the honeybee. R. Soc. Open Sci.5:171785. 10.1098/rsos.171785
111
Thomas-DanguinT.SindingC.RomagnyS.El MountassirF.AtanasovaB.Le BerreE.et al. (2014). The perception of odor objects in everyday life: a review on the processing of odor mixtures. Front. Psychol.5:504. 10.3389/fpsyg.2014.00504
112
TraubR. D.JefferysJ. G.WhittingtonM. A. (1997). Simulation of gamma rhythms in networks of interneurons and pyramidal cells. J. Comput. Neurosci.4, 141–150. 10.1023/A:1008839312043
113
VergaraA.FonollosaJ.MahiquesJ.TrincavelliM.RulkovN.HuertaR. (2013). On the performance of gas sensor arrays in open sampling systems using Inhibitory support vector machines. Sens. Actuat. B185, 462–477. 10.1016/j.snb.2013.05.027
114
WellisD. P.ScottJ. W.HarrisonT. A. (1989). Discrimination among odorants by single neurons of the rat olfactory bulb. J. Neurophysiol.61, 1161–1177. 10.1152/jn.1989.61.6.1161
115
WerthJ. C.ClelandT. A. (2019). Dynamics of spike time encoding in the olfactory bulb. Soc. Neurosci. Abstr. 399.07. Available online at: https://www.abstractsonline.com/pp8/#!/7883/presentation/70228
116
WillhiteD. C.NguyenK. T.MasurkarA. V.GreerC. A.ShepherdG. M.ChenW. R. (2006). Viral tracing identifies distributed columnar organization in the olfactory bulb. Proc. Natl. Acad. Sci. U.S.A.103, 12592–12597. 10.1073/pnas.0602032103
117
WilsonD. A.KadohisaM.FletcherM. L. (2006). Cortical contributions to olfaction: plasticity and perception. Semin. Cell. Dev. Biol.17, 462–470. 10.1016/j.semcdb.2006.04.008
118
XiaC. Z.AdjeiS.WessonD. W. (2015). Coding of odor stimulus features among secondary olfactory structures. J. Neurophysiol.114, 736–745. 10.1152/jn.00902.2014
119
XiongW.ChenW. R. (2002). Dynamic gating of spike propagation in the mitral cell lateral dendrites. Neuron34, 115–126. 10.1016/S0896-6273(02)00628-1
120
XuL.LiW.VoletiV.ZouD. J.HillmanE. M. C.FiresteinS. (2020). Widespread receptor-driven modulation in peripheral olfactory coding. Science368:eaaz5390. 10.1126/science.aaz5390
121
YokoiM.MoriK.NakanishiS. (1995). Refinement of odor molecule tuning by dendrodendritic synaptic inhibition in the olfactory bulb. Proc. Natl. Acad. Sci. U.S.A.92, 3371–3375. 10.1073/pnas.92.8.3371
122
ZakJ. D.ReddyG.VergassolaM.MurthyV. N. (2020). Antagonistic odor interactions in olfactory sensory neurons are widespread in freely breathing mice. Nat. Commun.11:3350. 10.1038/s41467-020-17124-5
123
ZavitzD.YoungstromI. A.BorisyukA.WachowiakM. (2020). Effect of interglomerular inhibitory networks on olfactory bulb odor representations. J. Neurosci.40, 5954–5969. 10.1523/JNEUROSCI.0233-20.2020
124
ZoharO.ShackletonT. M.PalmerA. R.ShamirM. (2013). The effect of correlated neuronal firing and neuronal heterogeneity on population coding accuracy in guinea pig inferior colliculus. PLoS ONE8:e81660. 10.1371/journal.pone.0081660
Summary
Keywords
spike synchronization, category learning, perceptual learning, olfaction, plasticity, neuromorphic, computational modeling, learning in the wild
Citation
Cleland TA and Borthakur A (2020) A Systematic Framework for Olfactory Bulb Signal Transformations. Front. Comput. Neurosci. 14:579143. doi: 10.3389/fncom.2020.579143
Received
01 July 2020
Accepted
17 August 2020
Published
23 September 2020
Volume
14 - 2020
Edited by
Justus V. Verhagen, John B. Pierce Laboratory, United States
Reviewed by
Daniel W. Wesson, University of Florida, United States; Keiji Miura, Kwansei Gakuin University, Japan; Douglas Anthony Storace, Florida State University, United States
Updates

Check for updates
Copyright
© 2020 Cleland and Borthakur.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Thomas A. Cleland tac29@cornell.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.