Encoding whisker deflection velocity within the rodent barrel cortex using phase-delayed inhibition

Mechanism of carbachol-induced 40 Hz gamma oscillations and the effects of NMDA activation on oscillatory dynamics in a model of the CA3 subfield of the hippocampus

Gamma oscillations are a prominent feature of various neural systems, including the CA3 subfield of the hippocampus. In CA3, in vitro carbachol application induces ∼40 Hz gamma oscillations in the network of glutamatergic excitatory pyramidal neurons (PNs) and local GABAergic inhibitory neurons (INs). Activation of NMDA receptors within CA3 leads to an increase in the frequency of carbachol-induced oscillations to ∼60 Hz, a broadening of the distribution of individual oscillation cycle frequencies, and a decrease in the time lag between PN and IN spike bursts. In this work, we develop a biophysical integrate-and-fire model of the CA3 subfield, we show that the dynamics of our model are in concordance with physiological observations, and we provide computational support for the hypothesis that the 'E-I' mechanism is responsible for the emergence of ∼40 Hz gamma oscillations in the absence of NMDA activation. We then incorporate NMDA receptors into our CA3 model, and we show that our model exhibits the increase in gamma oscillation frequency, broadening of the cycle frequency distribution, and decrease in the time lag between PN and IN spike bursts observed experimentally. Remarkably, we find an inverse relationship in our model between the net NMDA current delivered to PNs and INs in an oscillation cycle and cycle frequency. Furthermore, we find a disparate effect of NMDA receptors on PNs versus INs - we show that NMDA receptors on INs tend to increase oscillation frequency, while NMDA receptors on PNs tend to slightly decrease or not affect oscillation frequency. We find that these observations can be explained if NMDA activity above a threshold level causes a shift in the mechanism underlying gamma oscillations; in the absence of NMDA receptors, the 'E-I' mechanism is primarily responsible for the generation of gamma oscillations (at 40 Hz), while when NMDA receptors are active, the mechanism of gamma oscillations shifts to the 'I-I' mechanism, and we argue that within the 'I-I' regime (which displays a higher baseline oscillation frequency of ∼60 Hz), slight changes in the level of NMDA activity are inversely related to cycle frequency.

A model of lateral interactions as the origin of multiwhisker receptive fields in rat barrel cortex

While cells within barrel cortex respond primarily to deflections of their principal whisker (PW), they also exhibit responses to non-principal, or adjacent, whiskers (AWs), albeit responses with diminished amplitudes and longer latencies. The origin of multiwhisker receptive fields of barrel cells remains a point of controversy within the experimental literature, with three contending possibilities: (i) barrel cells inherit their AW responses from the AW responses of thalamocortical (TC) cells within their aligned barreloid; (ii) the axons of TC cells within a barreloid ramify to innervate multiple barrels, rather than only terminating within their aligned barrel; (iii) lateral intracortical transmission between barrels conveys AW responsivity to barrel cells. In this work, we develop a detailed, biologically plausible model of multiple barrels in order to examine possibility (iii); in order to isolate the dynamics that possibility (iii) entails, we incorporate lateral connections between barrels while assuming that TC cells respond only to their PW and that TC cell axons are confined to their home barrel. We show that our model is capable of capturing a broad swath of experimental observations on multiwhisker receptive field dynamics within barrels, and we compare and contrast the dynamics of this model with model dynamics from prior work in which employ a similar general modeling strategy to examine possibility (i).

Olfactory encoding within the insect antennal lobe: The emergence and role of higher order temporal correlations in the dynamics of antennal lobe spiking activity

In this review, we focus on the antennal lobe (AL) of three insect species - the fruit fly, sphinx moth, and locust. We first review the experimentally elucidated anatomy and physiology of the early olfactory system of each species; empirical studies of AL activity, however, often focus on assessing firing rates (averaged over time scales of about 100 ms), and hence the AL odor code is often analyzed in terms of a temporally evolving vector of firing rates. However, such a perspective necessarily misses the possibility of higher order temporal correlations in spiking activity within a single cell and across multiple cells over shorter time scales (of about 10 ms). Hence, we then review our prior theoretical work, where we constructed biophysically detailed, species-specific AL models within the fly, moth, and locust, finding that in each case higher order temporal correlations in spiking naturally emerge from model dynamics (i.e., without a prioriincorporation of elements designed to produce correlated activity). We therefore use our theoretical work to argue the perspective that temporal correlations in spiking over short time scales, which have received little experimental attention to-date, may provide valuable coding dimensions (complementing the coding dimensions provided by the vector of firing rates) that nature has exploited in the encoding of odors within the AL. We further argue that, if the AL does indeed utilize temporally correlated activity to represent odor information, such an odor code could be naturally and easily deciphered within the Mushroom Body.

Analysis of feedforward mechanisms of multiwhisker receptive field generation in a model of the rat barrel cortex

There is substantial anatomical segregation in the organization of the rodent barrel system - each whisker on the mystacial pad sends input to TC cells within a dedicated thalamic barreloid, which in turn innervates a corresponding cortical barrel, and RS cells within a barrel respond primarily to deflections of the corresponding whisker at the beginning of the dedicated transmission line (the principal whisker, PW). However, it is also well-established that barrel cells exhibit multiwhisker receptive fields (RFs), and display lower amplitude, longer latency responses to deflections of non-PWs (or adjacent whiskers, AWs). There is considerable controversy regarding the origin of such multiwhisker RFs; three possibilities include: (i) TC cells within a barreloid respond to multiple whiskers, and barrel RS cells simply inherit multiwhisker responses from their aligned barreloid; (ii) TC cells respond only to the PW, but individual barreloids innervate multiple barrels; (iii) multiwhisker responses of barrel cells arise from lateral corticocortical (barrel-to-barrel) synaptic transmission. Ablation studies attempting to pinpoint the source of RS cell AW responses are often contradictory (though experimental work tends to suggest possibilities (i) or (iii) to be most plausible), and hence it is important to carefully evaluate these hypotheses in terms of available physiological data on barreloid and barrel response dynamics. In this work, I employ a biologically detailed model of the rat barrel cortex to evaluate possibility (i), and I show that, within the model, hypothesis (i) is capable of explaining a broad range of the available physiological data on responses to single (PW or AW) deflections and paired whisker deflections (AW deflection followed by PW deflection), as well as the dependence of such responses on the angular direction of whisker deflection. In particular, the model shows that barrel RS cells can exhibit AW direction tuning despite the fact that barreloid to barrel wiring has no systematic dependence on the AW direction preference of TC cells. Future modeling work will examine the other possibilities for the generation of multiwhisker RS cell RFs, and compare and contrast the different possible mechanisms within the context of available experimental data.

Effects of Adaptation on Discrimination of Whisker Deflection Velocity and Angular Direction in a Model of the Barrel Cortex

Two important stimulus features represented within the rodent barrel cortex are velocity and angular direction of whisker deflection. Each cortical barrel receives information from thalamocortical (TC) cells that relay information from a single whisker, and TC input is decoded by barrel regular-spiking (RS) cells through a feedforward inhibitory architecture (with inhibition delivered by cortical fast-spiking or FS cells). TC cells encode deflection velocity through population synchrony, while deflection direction is encoded through the distribution of spike counts across the TC population. Barrel RS cells encode both deflection direction and velocity with spike rate, and are divided into functional domains by direction preference. Following repetitive whisker stimulation, system adaptation causes a weakening of synaptic inputs to RS cells and diminishes RS cell spike responses, though evidence suggests that stimulus discrimination may improve following adaptation. In this work, I construct a model of the TC, FS, and RS cells comprising a single barrel system—the model incorporates realistic synaptic connectivity and dynamics and simulates both angular direction (through the spatial pattern of TC activation) and velocity (through synchrony of the TC population spikes) of a deflection of the primary whisker, and I use the model to examine direction and velocity selectivity of barrel RS cells before and after adaptation. I find that velocity and direction selectivity of individual RS cells (measured over multiple trials) sharpens following adaptation, but stimulus discrimination using a simple linear classifier by the RS population response during a single trial (a more biologically meaningful measure than single cell discrimination over multiple trials) exhibits strikingly different behavior—velocity discrimination is similar both before and after adaptation, while direction classification improves substantially following adaptation. This is the first model, to my knowledge, that simulates both whisker deflection velocity and angular direction and examines the ability of the RS population response to pinpoint both stimulus features within the context of adaptation.

Spiking and Excitatory/Inhibitory Input Dynamics of Barrel Cells in Response to Whisker Deflections of Varying Velocity and Angular Direction

The spiking of barrel regular-spiking (RS) cells is tuned for both whisker deflection direction and velocity. Velocity tuning arises due to thalamocortical (TC) synchrony (but not spike quantity) varying with deflection velocity, coupled with feedforward inhibition, while direction selectivity is not fully understood, though may be due partly to direction tuning of TC spiking. Data show that as deflection direction deviates from the preferred direction of an RS cell, excitatory input to the RS cell diminishes minimally, but temporally shifts to coincide with the time-lagged inhibitory input. This work constructs a realistic large-scale model of a barrel; model RS cells exhibit velocity and direction selectivity due to TC input dynamics, with the experimentally observed sharpening of direction tuning with decreasing velocity. The model puts forth the novel proposal that RS→RS synapses can naturally and simply account for the unexplained direction dependence of RS cell inputs - as deflection direction deviates from the preferred direction of an RS cell, and TC input declines, RS→RS synaptic transmission buffers the decline in total excitatory input and causes a shift in timing of the excitatory input peak from the peak in TC input to the delayed peak in RS input. The model also provides several experimentally testable predictions on the velocity dependence of RS cell inputs. This model is the first, to my knowledge, to study the interaction of direction and velocity and propose physiological mechanisms for the stimulus dependence in the timing and amplitude of RS cell inputs.