A coupled-oscillator model of olfactory bulb gamma oscillations

Heterogeneous quantization regularizes spiking neural network activity

The learning and recognition of object features from unregulated input has been a longstanding challenge for artificial intelligence systems. Brains, on the other hand, are adept at learning stable sensory representations given noisy observations, a capacity mediated by a cascade of signal conditioning steps informed by domain knowledge. The olfactory system, in particular, solves a source separation and denoising problem compounded by concentration variability, environmental interference, and unpredictably correlated sensor affinities using a plastic network that requires statistically well-behaved input. We present a data-blind neuromorphic signal conditioning strategy, based on the biological system architecture, that normalizes and quantizes analog data into spike-phase representations, thereby transforming uncontrolled sensory input into a regular form with minimal information loss. Normalized input is delivered to a column of spiking principal neurons via heterogeneous synaptic weights; this gain diversification strategy regularizes neuronal utilization, yoking total activity to the network's operating range and rendering internal representations robust to uncontrolled open-set stimulus variance. To dynamically optimize resource utilization while balancing activity regularization and resolution, we supplement this mechanism with a data-aware calibration strategy in which the range and density of the quantization weights adapt to accumulated input statistics.

Sex differences in olfactory behavior and neurophysiology in Long Evans rats

In many species, olfactory abilities in females are more acute than those in males. Studies in humans show that women have lower olfactory thresholds and are better able to discriminate and identify odors than men. In mice, odorants elicit faster activation from a larger number of olfactory bulb glomeruli in females than in males. Our study explores sex differences in olfaction in Long Evans rats from a behavioral and electrophysiological perspective. Local field potentials (LFPs) in the olfactory bulb (OB) represent the coordinated activity of bulbar neurons. Olfactory gamma (65-120 Hz) and beta (15-30 Hz) oscillations have been functionally linked to odor perception. Spontaneous and odor-evoked OB LFPs were recorded from awake rats at the same time for 12 days. Odors used included urine of both sexes and monomolecular odorants characterized previously for correlation of volatility with behavior and OB oscillations. Sampling duration in a habituation context, baseline gamma and beta power, and odor-elicited beta and gamma power were analyzed. We find that females sample odorants for a shorter duration than males (just over 1-s difference). Although baseline gamma and beta power do not show significant differences between the two sexes, odor-elicited gamma and beta power in females is significantly lower than in males. Neither sampling duration nor beta and gamma power in females varied systematically with day of estrus. We further verify that variance of these behavioral and physiological measures is not different across sexes, adding to growing evidence that researchers need not be concerned about often-claimed additional variance in female subjects.NEW & NOTEWORTHY Olfaction plays a large role in evolutionary processes. However, we know little about sex differences in olfactory bulb neurophysiology, and many scientists believe that females are more variable because of estrus. We show that female rats sniff odors for shorter durations than males and have lower power in neural oscillations related to cognition. Estrus was not related to variance in any measures. Finally, males and females show equal variance on these behavioral and physiological processes.

Fast-spiking interneuron detonation drives high-fidelity inhibition in the olfactory bulb

Inhibitory circuits in the mammalian olfactory bulb (OB) dynamically reformat olfactory information as it propagates from peripheral receptors to downstream cortex. To gain mechanistic insight into how specific OB interneuron types support this sensory processing, we examine unitary synaptic interactions between excitatory mitral and tufted cells (MTCs), the OB projection neurons, and a conserved population of anaxonic external plexiform layer interneurons (EPL-INs) using pair and quartet whole-cell recordings in acute mouse brain slices. Physiological, morphological, neurochemical, and synaptic analyses divide EPL-INs into distinct subtypes and reveal that parvalbumin-expressing fast-spiking EPL-INs (FSIs) perisomatically innervate MTCs with release-competent dendrites and synaptically detonate to mediate fast, short-latency recurrent and lateral inhibition. Sparse MTC synchronization supralinearly increases this high-fidelity inhibition, while sensory afferent activation combined with single-cell silencing reveals that individual FSIs account for a substantial fraction of total network-driven MTC lateral inhibition. OB output is thus powerfully shaped by detonation-driven high-fidelity perisomatic inhibition.

Fast-spiking interneuron detonation drives high-fidelity inhibition in the olfactory bulb

Inhibitory circuits in the mammalian olfactory bulb (OB) dynamically reformat olfactory information as it propagates from peripheral receptors to downstream cortex. To gain mechanistic insight into how specific OB interneuron types support this sensory processing, we examine unitary synaptic interactions between excitatory mitral and tufted cells (MTCs), the OB projection cells, and a conserved population of anaxonic external plexiform layer interneurons (EPL-INs) using pair and quartet whole-cell recordings in acute mouse brain slices. Physiological, morphological, neurochemical, and synaptic analyses divide EPL-INs into distinct subtypes and reveal that parvalbumin-expressing fast-spiking EPL-INs (FSIs) perisomatically innervate MTCs with release-competent dendrites and synaptically detonate to mediate fast, short-latency recurrent and lateral inhibition. Sparse MTC synchronization supralinearly increases this high-fidelity inhibition, while sensory afferent activation combined with single-cell silencing reveals that individual FSIs account for a substantial fraction of total network-driven MTC lateral inhibition. OB output is thus powerfully shaped by detonation-driven high-fidelity perisomatic inhibition.

Coherent olfactory bulb gamma oscillations arise from coupling independent columnar oscillators

Spike timing-based representations of sensory information depend on embedded dynamical frameworks within neuronal networks that establish the rules of local computation and interareal communication. Here, we investigated the dynamical properties of olfactory bulb circuitry in mice of both sexes using microelectrode array recordings from slice and in vivo preparations. Neurochemical activation or optogenetic stimulation of sensory afferents evoked persistent gamma oscillations in the local field potential. These oscillations arose from slower, GABA(A) receptor-independent intracolumnar oscillators coupled by GABA(A)-ergic synapses into a faster, broadly coherent network oscillation. Consistent with the theoretical properties of coupled-oscillator networks, the spatial extent of zero-phase coherence was bounded in slices by the reduced density of lateral interactions. The intact in vivo network, however, exhibited long-range lateral interactions that suffice in simulation to enable zero-phase gamma coherence across the olfactory bulb. The timing of action potentials in a subset of principal neurons was phase-constrained with respect to evoked gamma oscillations. Coupled-oscillator dynamics in olfactory bulb thereby enable a common clock, robust to biological heterogeneities, that is capable of supporting gamma-band spike synchronization and phase coding across the ensemble of activated principal neurons.NEW & NOTEWORTHY Odor stimulation evokes rhythmic gamma oscillations in the field potential of the olfactory bulb, but the dynamical mechanisms governing these oscillations have remained unclear. Establishing these mechanisms is important as they determine the biophysical capacities of the bulbar circuit to, for example, maintain zero-phase coherence across a spatially extended network, or coordinate the timing of action potentials in principal neurons. These properties in turn constrain and suggest hypotheses of sensory coding.

Stimulus-Induced Theta-Band LFP Oscillations Format Neuronal Representations of Social Chemosignals in the Mouse Accessory Olfactory Bulb

Social communication is crucial for the survival of many species. In most vertebrates, a dedicated chemosensory system, the vomeronasal system (VNS), evolved to process ethologically relevant chemosensory cues. The first central processing stage of the VNS is the accessory olfactory bulb (AOB), which sends information to downstream brain regions via AOB mitral cells (AMCs). Recent studies provided important insights about the functional properties of AMCs, but little is known about the principles that govern their coordinated activity. Here, we recorded local field potentials (LFPs) and single-unit activity in the AOB of adult male and female mice during presentation of natural stimuli. Our recordings reveal prominent LFP theta-band oscillatory episodes with a characteristic spatial pattern across the AOB. Throughout an experiment, the AOB network shows varying degrees of similarity to this pattern, in a manner that depends on the sensory stimulus. Analysis of LFP signal polarity and single-unit activity indicates that oscillatory episodes are generated locally within the AOB, likely representing a reciprocal interaction between AMCs and granule cells. Notably, spike times of many AMCs are constrained to the negative LFP oscillation phase in a manner that can drastically affect integration by downstream processing stages. Based on these observations, we propose that LFP oscillations may gate, bind, and organize outgoing signals from individual AOB neurons to downstream processing stages. Our findings suggest that, as in other neuronal systems and brain regions, population-level oscillations play a key role in organizing and enhancing transmission of socially relevant chemosensory information.SIGNIFICANCE STATEMENT The accessory olfactory bulb (AOB) is the first central stage of the vomeronasal system, a chemosensory system dedicated to processing cues from other organisms. Information from the AOB is conveyed to other brain regions via activity of its principal neurons, AOB mitral cells (AMCs). Here, we show that socially relevant sensory stimulation of the mouse vomeronasal system leads not only to changes in AMC activity, but also to distinct theta-band (∼5 Hz) oscillatory episodes in the local field potential. Notably AMCs favor the negative phase of these oscillatory events. Our findings suggest a novel mechanism for the temporal coordination of distributed patterns of neuronal activity, which can serve to efficiently activate downstream processing stages.

Odor modality is transmitted to cortical brain regions from the olfactory bulb

Odor perception is the impetus for important animal behaviors with two predominate modes of processing: odors pass through the front of the nose (orthonasal) while inhaling and sniffing, or through the rear (retronasal) during exhalation and while eating. Despite the importance of olfaction for an animal's well-being and that ortho and retro naturally occur, it is unknown how the modality (ortho vs. retro) is even transmitted to cortical brain regions, which could significantly affect how odors are processed and perceived. Using multielectrode array recordings in tracheotomized anesthetized rats, which decouples ortho-retro modality from breathing, we show that mitral cells in rat olfactory bulb can reliably and directly transmit orthonasal versus retronasal modality with ethyl butyrate, a common food odor. Drug manipulations affecting synaptic inhibition via GABAA lead to worse decoding of ortho versus retro, independent of whether overall inhibition increases or decreases, suggesting that the olfactory bulb circuit may naturally favor encoding this important aspect of odors. Detailed data analysis paired with a firing rate model that captures population trends in spiking statistics shows how this circuit can encode odor modality. We have not only demonstrated that ortho/retro information is encoded to downstream brain regions but also used modeling to demonstrate a plausible mechanism for this encoding; due to synaptic adaptation, it is the slower time course of the retronasal stimulation that causes retronasal responses to be stronger and less sensitive to inhibitory drug manipulations than orthonasal responses.NEW & NOTEWORTHY Whether ortho (sniffing odors) versus retro (exhalation and eating) is encoded from the olfactory bulb to other brain areas is not completely known. Using multielectrode array recordings in anesthetized rats, we show that the olfactory bulb transmits this information downstream via spikes. Altering inhibition degrades ortho/retro information on average. We use theory and computation to explain our results, which should have implications on cortical processing considering that only food odors occur retronasally.

The Effects of Background Noise on a Biophysical Model of Olfactory Bulb Mitral Cells

The spiking activity of mitral cells (MC) in the olfactory bulb is a key attribute in olfactory sensory information processing to downstream cortical areas. A more detailed understanding of the modulation of MC spike statistics could shed light on mechanistic studies of olfactory bulb circuits and olfactory coding. We study the spike response of a recently developed single-compartment biophysical MC model containing seven known ionic currents and calcium dynamics subject to constant current input with background white noise. We observe rich spiking dynamics even with constant current input, including multimodal peaks in the interspike interval distribution (ISI). Although weak-to-moderate background noise for a fixed current input does not change the firing rate much, the spike dynamics can change dramatically, exhibiting non-monotonic spike variability not commonly observed in standard neuron models. We explain these dynamics with a phenomenological model of the ISI probability density function. Our study clarifies some of the complexities of MC spiking dynamics.

An ErbB4-Positive Neuronal Network in the Olfactory Bulb for Olfaction

Olfactory information is relayed and processed in the olfactory bulb (OB). Mitral cells, the principal output excitatory neurons of the OB, are controlled by multiple types of interneurons. However, mechanisms that regulate the activity of OB interneurons are not well understood. We provide evidence that the transmembrane tyrosine kinase ErbB4 is selectively expressed in subsets of OB inhibitory neurons in both male and female mice. ErbB4-positive (ErbB4+) neurons are mainly located in the glomerular layer (GL) and granule cell layer (GCL) and do not express previously defined markers. Optogenetic activation of GL-ErbB4+ neurons promotes theta oscillation, whereas activation of those in the GCL generates γ oscillations. Stimulation of OB slices with NRG1, a ligand that activates ErbB4, increases GABA transmission onto mitral cells, suggesting a role of OB NRG1-ErbB4 signaling in olfaction. In accord, ErbB4 mutant mice or acute inhibition of ErbB4 by a chemical genetic approach diminishes GABA transmission, reduces bulbar local field potential power, increases the threshold of olfactory sensitivity, and impairs odor discrimination. Together, these results identified a bulbar inhibitory network of ErbB4+ neurons for olfaction. Considering that both Nrg1 and Erbb4 are susceptibility genes for neuropsychiatric disorders, our study provides insight into pathologic mechanisms of olfactory malfunctions in these disorders.SIGNIFICANCE STATEMENT This study demonstrates that ErbB4+ neurons are a new subset of olfactory bulb inhibitory neurons in the glomerular layer and granule cell layer that innervate mitral cells and ErbB4- cells. They regulate olfaction by controlling local synchrony and distinct oscillations. ErbB4 inhibition diminishes GABA transmission, reduces bulbar local field potential power, increases the threshold of olfactory sensitivity, and impairs odor discrimination. Our results provide insight into pathophysiological mechanism of olfaction deficits in brain disorders associated with Nrg1 or Erbb4 mutations.

The voltage and spiking responses of subthreshold resonant neurons to structured and fluctuating inputs: persistence and loss of resonance and variability

We systematically investigate the response of neurons to oscillatory currents and synaptic-like inputs and we extend our investigation to non-structured synaptic-like spiking inputs with more realistic distributions of presynaptic spike times. We use two types of chirp-like inputs consisting of (i) a sequence of cycles with discretely increasing frequencies over time, and (ii) a sequence having the same cycles arranged in an arbitrary order. We develop and use a number of frequency-dependent voltage response metrics to capture the different aspects of the voltage response, including the standard impedance (Z) and the peak-to-trough amplitude envelope ([Formula: see text]) profiles. We show that Z-resonant cells (cells that exhibit subthreshold resonance in response to sinusoidal inputs) also show [Formula: see text]-resonance in response to sinusoidal inputs, but generally do not (or do it very mildly) in response to square-wave and synaptic-like inputs. In the latter cases the resonant response using Z is not predictive of the preferred frequencies at which the neurons spike when the input amplitude is increased above subthreshold levels. We also show that responses to conductance-based synaptic-like inputs are attenuated as compared to the response to current-based synaptic-like inputs, thus providing an explanation to previous experimental results. These response patterns were strongly dependent on the intrinsic properties of the participating neurons, in particular whether the unperturbed Z-resonant cells had a stable node or a focus. In addition, we show that variability emerges in response to chirp-like inputs with arbitrarily ordered patterns where all signals (trials) in a given protocol have the same frequency content and the only source of uncertainty is the subset of all possible permutations of cycles chosen for a given protocol. This variability is the result of the multiple different ways in which the autonomous transient dynamics is activated across cycles in each signal (different cycle orderings) and across trials. We extend our results to include high-rate Poisson distributed current- and conductance-based synaptic inputs and compare them with similar results using additive Gaussian white noise. We show that the responses to both Poisson-distributed synaptic inputs are attenuated with respect to the responses to Gaussian white noise. For cells that exhibit oscillatory responses to Gaussian white noise (band-pass filters), the response to conductance-based synaptic inputs are low-pass filters, while the response to current-based synaptic inputs may remain band-pass filters, consistent with experimental findings. Our results shed light on the mechanisms of communication of oscillatory activity among neurons in a network via subthreshold oscillations and resonance and the generation of network resonance.

Connectivity and dynamics in the olfactory bulb

Dendrodendritic interactions between excitatory mitral cells and inhibitory granule cells in the olfactory bulb create a dense interaction network, reorganizing sensory representations of odors and, consequently, perception. Large-scale computational models are needed for revealing how the collective behavior of this network emerges from its global architecture. We propose an approach where we summarize anatomical information through dendritic geometry and density distributions which we use to calculate the connection probability between mitral and granule cells, while capturing activity patterns of each cell type in the neural dynamical systems theory of Izhikevich. In this way, we generate an efficient, anatomically and physiologically realistic large-scale model of the olfactory bulb network. Our model reproduces known connectivity between sister vs. non-sister mitral cells; measured patterns of lateral inhibition; and theta, beta, and gamma oscillations. The model in turn predicts testable relationships between network structure and several functional properties, including lateral inhibition, odor pattern decorrelation, and LFP oscillation frequency. We use the model to explore the influence of cortex on the olfactory bulb, demonstrating possible mechanisms by which cortical feedback to mitral cells or granule cells can influence bulbar activity, as well as how neurogenesis can improve bulbar decorrelation without requiring cell death. Our methodology provides a tractable tool for other researchers.

Cell and circuit origins of fast network oscillations in the mammalian main olfactory bulb

Neural synchrony generates fast network oscillations throughout the brain, including the main olfactory bulb (MOB), the first processing station of the olfactory system. Identifying the mechanisms synchronizing neurons in the MOB will be key to understanding how network oscillations support the coding of a high-dimensional sensory space. Here, using paired recordings and optogenetic activation of glomerular sensory inputs in MOB slices, we uncovered profound differences in principal mitral cell (MC) vs. tufted cell (TC) spike-time synchrony: TCs robustly synchronized across fast- and slow-gamma frequencies, while MC synchrony was weaker and concentrated in slow-gamma frequencies. Synchrony among both cell types was enhanced by shared glomerular input but was independent of intraglomerular lateral excitation. Cell-type differences in synchrony could also not be traced to any difference in the synchronization of synaptic inhibition. Instead, greater TC than MC synchrony paralleled the more periodic firing among resonant TCs than MCs and emerged in patterns consistent with densely synchronous network oscillations. Collectively, our results thus reveal a mechanism for parallel processing of sensory information in the MOB via differential TC vs. MC synchrony, and further contrast mechanisms driving fast network oscillations in the MOB from those driving the sparse synchronization of irregularly firing principal cells throughout cortex.

Odor-evoked increases in olfactory bulb mitral cell spiking variability

The spiking variability of neural networks has important implications for how information is encoded to higher brain regions. It has been well documented by numerous labs in many cortical and motor regions that spiking variability decreases with stimulus onset, yet whether this principle holds in the OB has not been tested. In stark contrast to this common view, we demonstrate that the onset of sensory input can cause an increase in the variability of neural activity in the mammalian OB. We show this in both anesthetized and awake rodents. Furthermore, we use computational models to describe the mechanisms of this phenomenon. Our findings establish sensory evoked increases in spiking variability as a viable alternative coding strategy.

Differences in olfactory bulb mitral cell spiking with ortho- and retronasal stimulation revealed by data-driven models

The majority of olfaction studies focus on orthonasal stimulation where odors enter via the front nasal cavity, while retronasal olfaction, where odors enter the rear of the nasal cavity during feeding, is understudied. The coding of retronasal odors via coordinated spiking of neurons in the olfactory bulb (OB) is largely unknown despite evidence that higher level processing is different than orthonasal. To this end, we use multi-electrode array in vivo recordings of rat OB mitral cells (MC) in response to a food odor with both modes of stimulation, and find significant differences in evoked firing rates and spike count covariances (i.e., noise correlations). Differences in spiking activity often have implications for sensory coding, thus we develop a single-compartment biophysical OB model that is able to reproduce key properties of important OB cell types. Prior experiments in olfactory receptor neurons (ORN) showed retro stimulation yields slower and spatially smaller ORN inputs than with ortho, yet whether this is consequential for OB activity remains unknown. Indeed with these specifications for ORN inputs, our OB model captures the salient trends in our OB data. We also analyze how first and second order ORN input statistics dynamically transfer to MC spiking statistics with a phenomenological linear-nonlinear filter model, and find that retro inputs result in larger linear filters than ortho inputs. Finally, our models show that the temporal profile of ORN is crucial for capturing our data and is thus a distinguishing feature between ortho and retro stimulation, even at the OB. Using data-driven modeling, we detail how ORN inputs result in differences in OB dynamics and MC spiking statistics. These differences may ultimately shape how ortho and retro odors are coded.

Olfaction is a key sense for many cognitive and behavioral tasks, and is particularly unique because odors can naturally enter the nasal cavity from the front or rear, i.e., ortho- and retro-nasal, respectively. Yet little is known about the differences in coordinated spiking in the olfactory bulb with ortho versus retro stimulation, let alone how these different modes of olfaction may alter coding of odors. We simultaneously record many cells in rat olfactory bulb to assess the differences in spiking statistics, and develop a biophysical olfactory bulb network model to study the reasons for these differences. Using theoretical and computational methods, we find that the olfactory bulb transfers input statistics differently for retro stimulation relative to ortho stimulation. Furthermore, our models show that the temporal profile of inputs is crucial for capturing our data and is thus a distinguishing feature between ortho and retro stimulation, even at the olfactory bulb. Understanding the spiking dynamics of the olfactory bulb with both ortho and retro stimulation is a key step for ultimately understanding how the brain codes odors with different modes of olfaction.

Increased oscillatory power in a computational model of the olfactory bulb due to synaptic degeneration

Several neurodegenerative diseases impact the olfactory system, and in particular the olfactory bulb, early in disease progression. One mechanism by which damage occurs is via synaptic dysfunction. Here, we implement a computational model of the olfactory bulb and investigate the effect of weakened connection weights on network oscillatory behavior. Olfactory bulb network activity can be modeled by a system of equations that describes a set of coupled nonlinear oscillators. In this modeling framework, we propagate damage to synaptic weights using several strategies, varying from localized to global. Damage propagated in a dispersed or spreading manner leads to greater oscillatory power at moderate levels of damage. This increase arises from a higher average level of mitral cell activity due to a shift in the balance between excitation and inhibition. That this shift leads to greater oscillations depends critically on the nonlinearity of the activation function. Linearized analysis of the network dynamics predicts when this shift leads to loss of oscillatory activity. We thus demonstrate one potential mechanism involved in the increased gamma oscillations seen in some animal models of Alzheimer's disease, and we highlight the potential that pathological olfactory bulb behavior presents as an early biomarker of disease.

Computational Evidence for a Competitive Thalamocortical Model of Spikes and Spindle Activity in Rolandic Epilepsy

Rolandic epilepsy (RE) is the most common idiopathic focal childhood epilepsy syndrome, characterized by sleep-activated epileptiform spikes and seizures and cognitive deficits in school age children. Recent evidence suggests that this disease may be caused by disruptions to the Rolandic thalamocortical circuit, resulting in both an abundance of epileptiform spikes and a paucity of sleep spindles in the Rolandic cortex during non-rapid eye movement sleep (NREM); electrographic features linked to seizures and cognitive symptoms, respectively. The neuronal mechanisms that support the competitive shared thalamocortical circuitry between pathological epileptiform spikes and physiological sleep spindles are not well-understood. In this study we introduce a computational thalamocortical model for the sleep-activated epileptiform spikes observed in RE. The cellular and neuronal circuits of this model incorporate recent experimental observations in RE, and replicate the electrophysiological features of RE. Using this model, we demonstrate that: (1) epileptiform spikes can be triggered and promoted by either a reduced NMDA current or h-type current; and (2) changes in inhibitory transmission in the thalamic reticular nucleus mediates an antagonistic dynamic between epileptiform spikes and spindles. This work provides the first computational model that both recapitulates electrophysiological features and provides a mechanistic explanation for the thalamocortical switch between the pathological and physiological electrophysiological rhythms observed during NREM sleep in this common epileptic encephalopathy.

A Systematic Framework for Olfactory Bulb Signal Transformations

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.

Predicting brain organization with a computational model: 50-year perspective on lateral inhibition and oscillatory gating by dendrodendritic synapses

The first compartmental computer models of brain neurons using the Rall method predicted novel and unexpected dendrodendritic interactions between mitral and granule cells in the olfactory bulb. We review the models from a 50-year perspective on the work that has challenged, supported, and extended the original proposal that these interactions mediate both lateral inhibition and oscillatory activity, essential steps in the neural basis of olfactory processing and perception. We highlight strategies behind the neurophysiological experiments and the Rall methods that enhance the ability of detailed compartmental modeling to give counterintuitive predictions that lead to deeper insights into neural organization at the synaptic and circuit level. The application of these methods to mechanisms of neurogenesis and plasticity are exciting challenges for the future.

Top-Down Control of Inhibitory Granule Cells in the Main Olfactory Bulb Reshapes Neural Dynamics Giving Rise to a Diversity of Computations

Growing evidence shows that top-down projections from excitatory neurons in piriform cortex selectively synapse onto local inhibitory granule cells in the main olfactory bulb, effectively gating their own inputs by controlling inhibition. An open question in olfaction is the role this feedback plays in shaping the dynamics of local circuits, and the resultant computational benefits it provides. Using rate models of neuronal firing in a network consisting of excitatory mitral and tufted cells, inhibitory granule cells and top-down piriform cortical neurons, we found that changes in the weight of feedback to inhibitory neurons generated diverse network dynamics and complex transitions between these dynamics. Changes in the weight of top-down feedback supported a number of computations, including both pattern separation and oscillatory synchrony. Additionally, the network could generate gamma oscillations though a mechanism we termed Top-down control of Inhibitory Neuron Gamma (TING). Collectively, these functions arose from a codimension-2 bifurcation in the dynamical system. Our results highlight a key role for this top-down feedback, gating inhibition to facilitate often diametrically different computations.

Rapid online learning and robust recall in a neuromorphic olfactory circuit

We present a neural algorithm for the rapid online learning and identification of odourant samples under noise, based on the architecture of the mammalian olfactory bulb and implemented on the Intel Loihi neuromorphic system. As with biological olfaction, the spike timing-based algorithm utilizes distributed, event-driven computations and rapid (one-shot) online learning. Spike timing-dependent plasticity rules operate iteratively over sequential gamma-frequency packets to construct odour representations from the activity of chemosensor arrays mounted in a wind tunnel. Learned odourants then are reliably identified despite strong destructive interference. Noise resistance is further enhanced by neuromodulation and contextual priming. Lifelong learning capabilities are enabled by adult neurogenesis. The algorithm is applicable to any signal identification problem in which high-dimensional signals are embedded in unknown backgrounds.

Learning improves decoding of odor identity with phase-referenced oscillations in the olfactory bulb

Local field potential oscillations reflect temporally coordinated neuronal ensembles-coupling distant brain regions, gating processing windows, and providing a reference for spike timing-based codes. In phase amplitude coupling (PAC), the amplitude of the envelope of a faster oscillation is larger within a phase window of a slower carrier wave. Here, we characterized PAC, and the related theta phase-referenced high gamma and beta power (PRP), in the olfactory bulb of mice learning to discriminate odorants. PAC changes throughout learning, and odorant-elicited changes in PRP increase for rewarded and decrease for unrewarded odorants. Contextual odorant identity (is the odorant rewarded?) can be decoded from peak PRP in animals proficient in odorant discrimination, but not in naïve mice. As the animal learns to discriminate the odorants the dimensionality of PRP decreases. Therefore, modulation of phase-referenced chunking of information in the course of learning plays a role in early sensory processing in olfaction.

A Spike Time-Dependent Online Learning Algorithm Derived From Biological Olfaction

We have developed a spiking neural network (SNN) algorithm for signal restoration and identification based on principles extracted from the mammalian olfactory system and broadly applicable to input from arbitrary sensor arrays. For interpretability and development purposes, we here examine the properties of its initial feedforward projection. Like the full algorithm, this feedforward component is fully spike timing-based, and utilizes online learning based on local synaptic rules such as spike timing-dependent plasticity (STDP). Using an intermediate metric to assess the properties of this initial projection, the feedforward network exhibits high classification performance after few-shot learning without catastrophic forgetting, and includes a none of the above outcome to reflect classifier confidence. We demonstrate online learning performance using a publicly available machine olfaction dataset with challenges including relatively small training sets, variable stimulus concentrations, and 3 years of sensor drift.

Rhythmic modulation of thalamic oscillations depends on intrinsic cellular dynamics

OBJECTIVE

Rhythmic brain stimulation has emerged as a powerful tool to modulate cognition and to target pathological oscillations related to neurological and psychiatric disorders. However, we lack a systematic understanding of how periodic stimulation interacts with endogenous neural activity as a function of the brain state and target.

APPROACH

To address this critical issue, we applied periodic stimulation to a unified biophysical thalamic network model that generates multiple distinct oscillations, and examined thoroughly the impact of rhythmic stimulation on different oscillatory states.

MAIN RESULTS

We found that rhythmic perturbation induces four basic response mechanisms: entrainment, acceleration, resonance and suppression. Importantly, the appearance and expression of these mechanisms depend highly on the intrinsic cellular dynamics in each state. Specifically, the low-threshold bursting of thalamocortical cells (TCs) in delta (δ) oscillation renders the network relatively insensitive to entrainment; the high-threshold bursting of TCs in alpha (α) oscillation leads to widespread oscillation suppression while the tonic spiking of TC cells in gamma (γ) oscillation results in prominent entrainment and resonance. In addition, we observed entrainment discontinuity during α oscillation that is mediated by firing pattern switching of high-threshold bursting TC cells. Furthermore, we demonstrate that direct excitatory stimulation of the lateral geniculate nucleus (LGN) entrains thalamic oscillations via an asymmetric Arnold tongue that favors higher frequency entrainment and resonance, while stimulation of the inhibitory circuit, the reticular nucleus, induces much weaker and more symmetric entrainment and resonance. These results support the notion that rhythmic stimulation engages brain oscillations in a state- and target-dependent manner.

SIGNIFICANCE

Overall, our study provides, for the first time, insights into how the biophysics of thalamic oscillations guide the emergence of complex, state-dependent mechanisms of target engagement, which can be leveraged for the future rational design of novel therapeutic stimulation modalities.

Generalized Cross-Frequency Decomposition: A Method for the Extraction of Neuronal Components Coupled at Different Frequencies

Perceptual, motor and cognitive processes are based on rich interactions between remote regions in the human brain. Such interactions can be carried out through phase synchronization of oscillatory signals. Neuronal synchronization has been primarily studied within the same frequency range, e.g., within alpha or beta frequency bands. Yet, recent research shows that neuronal populations can also demonstrate phase synchronization between different frequency ranges. An extraction of such cross-frequency interactions in EEG/MEG recordings remains, however, methodologically challenging. Here we present a new method for the robust extraction of cross-frequency phase-to-phase synchronized components. Generalized Cross-Frequency Decomposition (GCFD) reconstructs the time courses of synchronized neuronal components, their spatial filters and patterns. Our method extends the previous state of the art, Cross-Frequency Decomposition (CFD), to the whole range of frequencies: it works for any f₁ and f₂ whenever f₁:f₂ is a rational number. GCFD gives a compact description of non-linearly interacting neuronal sources on the basis of their cross-frequency phase coupling. We successfully validated the new method in simulations and tested it with real EEG recordings including resting state data and steady state visually evoked potentials (SSVEP).

Pulse-coupled mixed-mode oscillators: Cluster states and extreme noise sensitivity

Motivated by rhythms in the olfactory system of the brain, we investigate the synchronization of all-to-all pulse-coupled neuronal oscillators exhibiting various types of mixed-mode oscillations (MMOs) composed of sub-threshold oscillations (STOs) and action potentials ("spikes"). We focus particularly on the impact of the delay in the interaction. In the weak-coupling regime, we reduce the system to a Kuramoto-type equation with non-sinusoidal phase coupling and the associated Fokker-Planck equation. Its linear stability analysis identifies the appearance of various cluster states. Their type depends sensitively on the delay and the width of the pulses. Interestingly, long delays do not imply slow population rhythms, and the number of emerging clusters only loosely depends on the number of STOs. Direct simulations of the oscillator equations reveal that for quantitative agreement of the weak-coupling theory the coupling strength and the noise have to be extremely small. Even moderate noise leads to significant skipping of STO cycles, which can enhance the diffusion coefficient in the Fokker-Planck equation by two orders of magnitude. Introducing an effective diffusion coefficient extends the range of agreement significantly. Numerical simulations of the Fokker-Planck equation reveal bistability and solutions with oscillatory order parameters that result from nonlinear mode interactions. These are confirmed in simulations of the full spiking model.

Generative Biophysical Modeling of Dynamical Networks in the Olfactory System

Generative models are computational models designed to generate appropriate values for all of their embedded variables, thereby simulating the response properties of a complex system based on the coordinated interactions of a multitude of physical mechanisms. In systems neuroscience, generative models are generally biophysically based compartmental models of neurons and networks that are explicitly multiscale, being constrained by experimental data at multiple levels of organization from cellular membrane properties to large-scale network dynamics. As such, they are able to explain the origins of emergent properties in complex systems, and serve as tests of sufficiency and as quantitative instantiations of working hypotheses that may be too complex to simply intuit. Moreover, when adequately constrained, generative biophysical models are able to predict novel experimental outcomes, and consequently are powerful tools for experimental design. We here outline a general strategy for the iterative design and implementation of generative, multiscale biophysical models of neural systems. We illustrate this process using our ongoing, iteratively developing model of the mammalian olfactory bulb. Because the olfactory bulb exhibits diverse and interesting properties at multiple scales of organization, it is an attractive system in which to illustrate the value of generative modeling across scales.