Two distinct mechanisms shape the reliability of neural responses

Dynamics of interaction between IH and IKLT currents to mediate double resonances of medial superior olive neurons related to sound localization

Neurons in the medial superior olive (MSO) exhibit high frequency responses such as subthreshold resonance, which is helpful to sensitively detect a small difference in the arrival time of sounds between two ears for precise sound localization. Recently, except for the high frequency depolarization resonance mediated by a low threshold potassium (IKLT) current, a low frequency hyperpolarization resonance mediated by a hyperpolarization-activated cation (IH) current is observed in experiments on the MSO neurons, forming double resonances. The complex dynamics underlying double resonances are studied in an MSO neuron model in the present paper. Firstly, double resonances similar to the experimental observations are simulated as the resting membrane potential is between half-activation voltages of IH and IKLT currents, and stimulation current (IZAP) with large amplitude and exponentially increasing frequency is applied. Secondly, multiple effective factors to modulate double resonances are obtained. Especially, the decrease of time constant of IKLT current and increase of conductance of IH and IKLT currents can enhance the depolarization resonance frequency for precise sound localization. Last, different frequency responses of slow IH and fast IKLT currents in formation of the resonances are acquired. A middle phase difference between IZAP and IKLT currents appears at a high frequency, and the interaction between the positive part of IZAP and the negative IKLT current forms the depolarization resonance. Interaction between the negative part of IZAP and positive IH current with a middle phase difference results in hyperpolarization resonance at a low frequency. Furthermore, the phase difference between IZAP and resonance current can well explain the increase of depolarization resonance frequency modulated by the increase of conductance of IH or IKLT currents. The results present the dynamical and biophysical mechanisms for the double resonances mediated by two currents in the MSO neurons, which is helpful to enhance the depolarization resonance frequency for precise sound localization.

Persistent Interruption in Parvalbumin Positive Inhibitory Interneurons: Biophysical and Mathematical Mechanisms

Persistent activity in principal cells is a putative mechanism for maintaining memory traces during working memory. We recently demonstrated persistent interruption of firing in fast-spiking parvalbumin-expressing interneurons (PV-INs), a phenomenon which could serve as a substrate for persistent activity in principal cells through disinhibition lasting hundreds of milliseconds. Here, we find that hippocampal CA1 PV-INs exhibit type 2 excitability, like striatal and neocortical PV-INs. Modelling and mathematical analysis showed that the slowly inactivating potassium current Kv1 contributes to type 2 excitability, enables the multiple firing regimes observed experimentally in PV-INs, and provides a mechanism for robust persistent interruption of firing. Using a fast/slow separation of times scales approach with the Kv1 inactivation variable as a bifurcation parameter shows that the initial inhibitory stimulus stops repetitive firing by moving the membrane potential trajectory onto a co-existing stable fixed point corresponding to a non-spiking quiescent state. As Kv1 inactivation decays, the trajectory follows the branch of stable fixed points until it crosses a subcritical Hopf bifurcation then spirals out into repetitive firing. In a model describing entorhinal cortical PV-INs without Kv1, interruption of firing could be achieved by taking advantage of the bistability inherent in type 2 excitability based on a subcritical Hopf bifurcation, but the interruption was not robust to noise. Persistent interruption of firing is therefore broadly applicable to PV-INs in different brain regions but is only made robust to noise in the presence of a slow variable.

Synchrony in Networks of Type 2 Interneurons Is More Robust to Noise with Hyperpolarizing Inhibition Compared to Shunting Inhibition in Both the Stochastic Population Oscillator and the Coupled Oscillator Regimes

Synchronization in the gamma band (25-150 Hz) is mediated by PV+ inhibitory interneurons, and evidence is accumulating for the essential role of gamma oscillations in cognition. Oscillations can arise in inhibitory networks via synaptic interactions between individual oscillatory neurons (mean-driven) or via strong recurrent inhibition that destabilizes the stationary background firing rate in the fluctuation-driven balanced state, causing an oscillation in the population firing rate. Previous theoretical work focused on model neurons with Hodgkin's Type 1 excitability (integrators) connected by current-based synapses. Here we show that networks comprised of simple Type 2 oscillators (resonators) exhibit a supercritical Hopf bifurcation between synchrony and asynchrony and a gradual transition via cycle skipping from coupled oscillators to stochastic population oscillator (SPO), as previously shown for Type 1. We extended our analysis to homogeneous networks with conductance rather than current based synapses and found that networks with hyperpolarizing inhibitory synapses were more robust to noise than those with shunting synapses, both in the coupled oscillator and SPO regime. Assuming that reversal potentials are uniformly distributed between shunting and hyperpolarized values, as observed in one experimental study, converting synapses to purely hyperpolarizing favored synchrony in all cases, whereas conversion to purely shunting synapses made synchrony less robust except at very high conductance strengths. In mature neurons the synaptic reversal potential is controlled by chloride cotransporters that control the intracellular concentrations of chloride and bicarbonate ions, suggesting these transporters as a potential therapeutic target to enhance gamma synchrony and cognition.

Striatal Neurons Are Recruited Dynamically into Collective Representations of Self-Initiated and Learned Actions in Freely Moving Mice

Striatal spiny projection neurons are hyperpolarized-at-rest (HaR) and driven to action potential threshold by a small number of powerful inputs-an input-output configuration that is detrimental to response reliability. Because the striatum is important for habitual behaviors and goal-directed learning, we conducted a microendoscopic imaging in freely moving mice that express a genetically encoded Ca2+ indicator sparsely in striatal HaR neurons to evaluate their response reliability during self-initiated movements and operant conditioning. The sparse expression was critical for longitudinal studies of response reliability, and for studying correlations among HaR neurons while minimizing spurious correlations arising from contamination by the background signal. We found that HaR neurons are recruited dynamically into action representation, with distinct neuronal subsets being engaged in a moment-by-moment fashion. While individual neurons respond with little reliability, the population response remained stable across days. Moreover, we found evidence for the temporal coupling between neuronal subsets during conditioned (but not innate) behaviors.

Synchrony in Networks of Type 2 Interneurons is More Robust to Noise with Hyperpolarizing Inhibition Compared to Shunting Inhibition in Both the Stochastic Population Oscillator and the Coupled Oscillator Regimes

Synchronization in the gamma band (30-80 Hz) is mediated by PV+ inhibitory interneurons, and evidence is accumulating for the essential role of gamma oscillations in cognition. Oscillations can arise in inhibitory networks via synaptic interactions between individual oscillatory neurons (mean-driven) or via strong recurrent inhibition that destabilizes the stationary background firing rate in the fluctuation-driven balanced state, causing an oscillation in the population firing rate. Previous theoretical work focused on model neurons with Hodgkin's type 1 excitability (integrators) connected by current-based synapses. Here we show that networks comprised of simple type 2 oscillators (resonators) exhibit a supercritical Hopf bifurcation between synchrony and asynchrony and a gradual transition via cycle skipping from coupled oscillators to stochastic population oscillator, as previously shown for type 1. We extended our analysis to homogeneous networks with conductance rather than current based synapses and found that networks with hyperpolarizing inhibitory synapses were more robust to noise than those with shunting synapses, both in the coupled oscillator and stochastic population oscillator regime. Assuming that reversal potentials are uniformly distributed between shunting and hyperpolarized values, as observed in one experimental study, converting synapses to purely hyperpolarizing favored synchrony in all cases, whereas conversion to purely shunting synapses made synchrony less robust except at very high conductance strengths. In mature neurons the synaptic reversal potential is controlled by chloride cotransporters that control the intracellular concentrations of chloride and bicarbonate ions, suggesting these transporters as a potential therapeutic target to enhance gamma synchrony and cognition.

Significance Statement

Brain rhythms in the gamma frequency band (30-80 Hz) depend on the activity of inhibitory interneurons and evidence for a causal role for gamma oscillations in cognitive functions is accumulating. Here we extend previous studies on synchronization mechanisms to interneurons that have an abrupt threshold frequency below which they cannot sustain firing. In addition to current based synapses, we examined inhibitory networks with conductance based synapses. We found that if the reversal potential for inhibition was below the average membrane potential (hyperpolarizing), synchrony was more robust to noise than if the reversal potential was very close to the average potential (shunting). These results have implications for therapies to ameliorate cognitive deficits.

Simulations predict differing phase responses to excitation vs. inhibition in theta-resonant pyramidal neurons

Rhythmic activity is ubiquitous in neural systems, with theta-resonant pyramidal neurons integrating rhythmic inputs in many cortical structures. Impedance analysis has been widely used to examine frequency-dependent responses of neuronal membranes to rhythmic inputs, but it assumes that the neuronal membrane is a linear system, requiring the use of small signals to stay in a near-linear regime. However, postsynaptic potentials are often large and trigger nonlinear mechanisms (voltage-gated ion channels). The goals of this work were to 1) develop an analysis method to evaluate membrane responses in this nonlinear domain and 2) explore phase relationships between rhythmic stimuli and subthreshold and spiking membrane potential (Vmemb) responses in models of theta-resonant pyramidal neurons. Responses in these output regimes were asymmetrical, with different phase shifts during hyperpolarizing and depolarizing half-cycles. Suprathreshold theta-rhythmic stimuli produced nonstationary Vmemb responses. Sinusoidal inputs produced "phase retreat": action potentials occurred progressively later in cycles of the input stimulus, resulting from adaptation. Sinusoidal current with increasing amplitude over cycles produced "phase advance": action potentials occurred progressively earlier. Phase retreat, phase advance, and subthreshold phase shifts were modulated by multiple ion channel conductances. Our results suggest differential responses of cortical neurons depending on the frequency of oscillatory input, which will play a role in neuronal responses to shifts in network state. We hypothesize that intrinsic cellular properties complement network properties and contribute to in vivo phase-shift phenomena such as phase precession, seen in place and grid cells, and phase roll, also observed in hippocampal CA1 neurons.NEW & NOTEWORTHY We augmented electrical impedance analysis to characterize phase shifts between large-amplitude current stimuli and nonlinear, asymmetric membrane potential responses. We predict different frequency-dependent phase shifts in response excitation vs. inhibition, as well as shifts in spike timing over multiple input cycles, in theta-resonant pyramidal neurons. We hypothesize that these effects contribute to navigation-related phenomena such as phase precession and phase roll. Our neuron-level hypothesis complements, rather than falsifies, prior network-level explanations of these phenomena.

Physiological noise facilitates multiplexed coding of vibrotactile-like signals in somatosensory cortex

Neurons can use different aspects of their spiking to simultaneously represent (multiplex) different features of a stimulus. For example, some pyramidal neurons in primary somatosensory cortex (S1) use the rate and timing of their spikes to, respectively, encode the intensity and frequency of vibrotactile stimuli. Doing so has several requirements. Because they fire at low rates, pyramidal neurons cannot entrain 1:1 with high-frequency (100 to 600 Hz) inputs and, instead, must skip (i.e., not respond to) some stimulus cycles. The proportion of skipped cycles must vary inversely with stimulus intensity for firing rate to encode stimulus intensity. Spikes must phase-lock to the stimulus for spike times (intervals) to encode stimulus frequency, but, in addition, skipping must occur irregularly to avoid aliasing. Using simulations and in vitro experiments in which mouse S1 pyramidal neurons were stimulated with inputs emulating those induced by vibrotactile stimuli, we show that fewer cycles are skipped as stimulus intensity increases, as required for rate coding, and that intrinsic or synaptic noise can induce irregular skipping without disrupting phase locking, as required for temporal coding. This occurs because noise can modulate the reliability without disrupting the precision of spikes evoked by small-amplitude, fast-onset signals. Specifically, in the fluctuation-driven regime associated with sparse spiking, rate and temporal coding are both paradoxically improved by the strong synaptic noise characteristic of the intact cortex. Our results demonstrate that multiplexed coding by S1 pyramidal neurons is not only feasible under in vivo conditions, but that background synaptic noise is actually beneficial.

Firing statistics in the bistable regime of neurons with homoclinic spike generation

Neuronal voltage dynamics of regularly firing neurons typically has one stable attractor: either a fixed point (like in the subthreshold regime) or a limit cycle that defines the tonic firing of action potentials (in the suprathreshold regime). In two of the three spike onset bifurcation sequences that are known to give rise to all-or-none type action potentials, however, the resting-state fixed point and limit cycle spiking can coexist in an intermediate regime, resulting in bistable dynamics. Here, noise can induce switches between the attractors, i.e., between rest and spiking, and thus increase the variability of the spike train compared to neurons with only one stable attractor. Qualitative features of the resulting spike statistics depend on the spike onset bifurcations. This paper focuses on the creation of the spiking limit cycle via the saddle-homoclinic orbit (HOM) bifurcation and derives interspike interval (ISI) densities for a conductance-based neuron model in the bistable regime. The ISI densities of bistable homoclinic neurons are found to be unimodal yet distinct from the inverse Gaussian distribution associated with the saddle-node-on-invariant-cycle bifurcation. It is demonstrated that for the HOM bifurcation the transition between rest and spiking is mainly determined along the downstroke of the action potential-a dynamical feature that is not captured by the commonly used reset neuron models. The deduced spike statistics can help to identify HOM dynamics in experimental data.

Clusters of cooperative ion channels enable a membrane-potential-based mechanism for short-term memory

Across biological systems, cooperativity between proteins enables fast actions, supra-linear responses, and long-lasting molecular switches. In the nervous system, however, the function of cooperative interactions between voltage-dependent ionic channels remains largely unknown. Based on mathematical modeling, we here demonstrate that clusters of strongly cooperative ion channels can plausibly form bistable conductances. Consequently, clusters are permanently switched on by neuronal spiking, switched off by strong hyperpolarization, and remain in their state for seconds after stimulation. The resulting short-term memory of the membrane potential allows to generate persistent firing when clusters of cooperative channels are present together with non-cooperative spike-generating conductances. Dynamic clamp experiments in rodent cortical neurons confirm that channel cooperativity can robustly induce graded persistent activity - a single-cell based, multistable mnemonic firing mode experimentally observed in several brain regions. We therefore propose that ion channel cooperativity constitutes an efficient cell-intrinsic implementation for short-term memories at the voltage level.

Asymmetrical voltage response in resonant neurons shaped by nonlinearities

The conventional impedance profile of a neuron can identify the presence of resonance and other properties of the neuronal response to oscillatory inputs, such as nonlinear response amplifications, but it cannot distinguish other nonlinear properties such as asymmetries in the shape of the voltage response envelope. Experimental observations have shown that the response of neurons to oscillatory inputs preferentially enhances either the upper or lower part of the voltage envelope in different frequency bands. These asymmetric voltage responses arise in a neuron model when it is submitted to high enough amplitude oscillatory currents of variable frequencies. We show how the nonlinearities associated to different ionic currents or present in the model as captured by its voltage equation lead to asymmetrical response and how high amplitude oscillatory currents emphasize this response. We propose a geometrical explanation for the phenomenon where asymmetries result not only from nonlinearities in their activation curves but also from nonlinearites captured by the nullclines in the phase-plane diagram and from the system's time-scale separation. In addition, we identify an unexpected frequency-dependent pattern which develops in the gating variables of these currents and is a product of strong nonlinearities in the system as we show by controlling such behavior by manipulating the activation curve parameters. The results reported in this paper shed light on the ionic mechanisms by which brain embedded neurons process oscillatory information.

How to correctly quantify neuronal phase-response curves from noisy recordings

At the level of individual neurons, various coding properties can be inferred from the input-output relationship of a cell. For small inputs, this relation is captured by the phase-response curve (PRC), which measures the effect of a small perturbation on the timing of the subsequent spike. Experimentally, however, an accurate experimental estimation of PRCs is challenging. Despite elaborate measurement efforts, experimental PRC estimates often cannot be related to those from modeling studies. In particular, experimental PRCs rarely resemble the characteristic theoretical PRC expected close to spike initiation, which is indicative of the underlying spike-onset bifurcation. Here, we show for conductance-based model neurons that the correspondence between theoretical and measured phase-response curve is lost when the stimuli used for the estimation are too large. In this case, the derived phase-response curve is distorted beyond recognition and takes on a generic shape that reflects the measurement protocol and masks the spike-onset bifurcation. We discuss how to identify appropriate stimulus strengths for perturbation and noise-stimulation methods, which permit to estimate PRCs that reliably reflect the spike-onset bifurcation - a task that is particularly difficult if a lower bound for the stimulus amplitude is dictated by prominent intrinsic neuronal noise.

Diversity and Function of Somatostatin-Expressing Interneurons in the Cerebral Cortex

Inhibitory interneurons make up around 10-20% of the total neuron population in the cerebral cortex. A hallmark of inhibitory interneurons is their remarkable diversity in terms of morphology, synaptic connectivity, electrophysiological and neurochemical properties. It is generally understood that there are three distinct and non-overlapping interneuron classes in the mouse neocortex, namely, parvalbumin-expressing, 5-HT3A receptor-expressing and somatostatin-expressing interneuron classes. Each class is, in turn, composed of a multitude of subclasses, resulting in a growing number of interneuron classes and subclasses. In this review, I will focus on the diversity of somatostatin-expressing interneurons (SOM+ INs) in the cerebral cortex and elucidate their function in cortical circuits. I will then discuss pathological consequences of a malfunctioning of SOM+ INs in neurological disorders such as major depressive disorder, and present future avenues in SOM research and brain pathologies.

Gad1-promotor-driven GFP expression in non-GABAergic neurons of the nucleus endopiriformis in a transgenic mouse line

Transgenic animals have become a widely used model to identify and study specific cell types in whole organs. Promotor-driven reporter gene labeling of the cells under investigation has promoted experimental efficacy to a large degree. However, rigorous assessment of transgene expression specificity in these animal models is highly recommended to validate cellular identity and to isolate potentially mislabeled cell populations. Here, we report on one such mislabeled neuron population in a widely used transgenic mouse line in which GABAergic somatostatin-expressing interneurons (SOM^pos INs) are labeled by eGFP (so-called GIN mouse, FVB-Tg(GadGFP)45704Swn/J). These neurons represent a subpopulation of all SOM^pos INs. However, we report here on GFP labeling of non-GABAergic neurons in the nucleus endopiriformis of this mouse line.

Population dynamics and entrainment of basal ganglia pacemakers are shaped by their dendritic arbors

The theory of phase oscillators is an essential tool for understanding population dynamics of pacemaking neurons. GABAergic pacemakers in the substantia nigra pars reticulata (SNr), a main basal ganglia (BG) output nucleus, receive inputs from the direct and indirect pathways at distal and proximal regions of their dendritic arbors, respectively. We combine theory, optogenetic stimulation and electrophysiological experiments in acute brain slices to ask how dendritic properties impact the propensity of the various inputs, arriving at different locations along the dendrite, to recruit or entrain SNr pacemakers. By combining cable theory with sinusoidally-modulated optogenetic activation of either proximal somatodendritic regions or the entire somatodendritic arbor of SNr neurons, we construct an analytical model that accurately fits the empirically measured somatic current response to inputs arising from illuminating the soma and various portions of the dendritic field. We show that the extent of the dendritic tree that is illuminated generates measurable and systematic differences in the pacemaker’s phase response curve (PRC), causing a shift in its peak. Finally, we show that the divergent PRCs correctly predict differences in two major features of the collective dynamics of SNr neurons: the fidelity of population responses to sudden step-like changes in inputs; and the phase latency at which SNr neurons are entrained by rhythmic stimulation, which can occur in the BG under both physiological and pathophysiological conditions. Our novel method generates measurable and physiologically meaningful spatial effects, and provides the first empirical demonstration of how the collective responses of SNr pacemakers are determined by the transmission properties of their dendrites. SNr dendrites may serve to delay distal striatal inputs so that they impinge on the spike initiation zone simultaneously with pallidal and subthalamic inputs in order to guarantee a fair competition between the influence of the monosynaptic direct- and polysynaptic indirect pathways.

The substantia nigra pars reticulata (SNr) is a main output nucleus of the basal ganglia (BG), where inputs from the competing direct and indirect pathways converge onto the same neurons. Interestingly, these inputs are differentially distributed with direct and indirect pathway projections arriving at distal and proximal regions of the dendritic arbor, respectively. We employ a novel method combining theory with electrophysiological experiments and optogenetics to study the distinct effects of inputs arriving at different locations along the dendrite. Our approach represents a useful compromise between complexity and reduction in modelling. Our work addresses the question of high fidelity encoding of inputs by networks of neurons in the new context of pacemaking neurons, which are driven to fire by their intrinsic dynamics rather than by a network state. We provide the first empirical demonstration that dendritic delays can introduce latencies in the responses of a population of neurons that are commensurate with synaptic delays, suggesting a new role for SNr dendrites with implications for BG function.

Using computational models to predict in vivo synaptic inputs to interneuron specific 3 (IS3) cells of CA1 hippocampus that also allow their recruitment during rhythmic states

Brain coding strategies are enabled by the balance of synaptic inputs that individual neurons receive as determined by the networks in which they reside. Inhibitory cell types contribute to brain function in distinct ways but recording from specific, inhibitory cell types during behaviour to determine their contributions is highly challenging. In particular, the in vivo activities of vasoactive intestinal peptide-expressing interneuron specific 3 (IS3) cells in the hippocampus that only target other inhibitory cells are unknown at present. We perform a massive, computational exploration of possible synaptic inputs to IS3 cells using multi-compartment models and optimized synaptic parameters. We find that asynchronous, in vivo-like states that are sensitive to additional theta-timed inputs (8 Hz) exist when excitatory and inhibitory synaptic conductances are approximately equally balanced and with low numbers of activated synapses receiving correlated inputs. Specifically, under these balanced conditions, the input resistance is larger with higher mean spike firing rates relative to other activated synaptic conditions investigated. Incoming theta-timed inputs result in strongly increased spectral power relative to baseline. Thus, using a generally applicable computational approach we predict the existence and features of background, balanced states in hippocampal circuits.

Oscillations and Spike Entrainment

Oscillatory input to networks, as indicated by field potentials, must entrain neuronal firing to be a causal agent in brain activity. Even when the oscillatory input is prominent, entrainment of firing is not a foregone conclusion but depends on the intrinsic dynamics of the postsynaptic neurons, including cell type-specific resonances, and background firing rates. Within any local network of neurons, only a subset of neurons may have their firing entrained by an oscillating synaptic input, and oscillations of different frequency may engage separate subsets of neurons.

Two types of somatostatin-expressing GABAergic interneurons in the superficial layers of the mouse cingulate cortex

Somatostatin-expressing (SOM⁺), inhibitory interneurons represent a heterogeneous group of cells and given their remarkable diversity, classification of SOM⁺ interneurons remains a challenging task. Electrophysiological, morphological and neurochemical classes of SOM⁺ interneurons have been proposed in the past but it remains unclear as to what extent these classes are congruent. We performed whole-cell patch-clamp recordings from 127 GFP-labeled SOM⁺ interneurons ('GIN') of the superficial cingulate cortex with subsequent biocytin-filling and immunocytochemical labeling. Principal component analysis followed by k-means clustering predicted two putative subtypes of SOM⁺ interneurons, which we designated as group I and group II GIN. A key finding of our study is the fact that these electrophysiologically and morphologically distinct groups of SOM⁺ interneurons can be correlated with two neurochemical subtypes of SOM⁺ interneurons described recently in our laboratory. In particular, all SOM⁺ interneurons expressing calbindin but no calretinin could be classified as group I GIN, whereas all but one neuropeptide Y- and calretinin-positive interneurons were found in group II.

Predicting the response of striatal spiny neurons to sinusoidal input

Spike-timing effects of small-amplitude sinusoidal currents were measured in mouse striatal spiny neurons firing repetitively. Spike-timing reliability varied with the stimulus frequency. For frequencies near the cell's firing rate, the cells altered firing rate to match the stimulus and became phase locked to it. The stimulus phase of firing during lock depended on the stimulus frequency relative to the cell's unperturbed firing rate. Interspike intervals during sinusoidal stimulation were predicted using an iterative map constructed from the cells' phase-resetting curve. Variability of interspike intervals was reduced by stimulation at all frequencies higher than about half the cell's unperturbed rate, and interspike intervals were accurately predicted by the map. Long sequences of spike times were predicted by iterating on the map. The accuracy of that prediction varied with frequency. Spike time predictability was highest near and during phase lock. The map predicted the phase of firing on the input and its dependence on stimulus frequency. Prediction errors, when they occurred, were of two kinds: unpredicted variation in interspike interval from intrinsic cell noise and accumulation of prediction errors from previous interspike intervals. Each type of prediction error arose from a different mechanism, and their impact was also predicted from the phase model. When two oscillatory input currents were presented simultaneously, striatal neurons responded selectively to only one of them, the one closest in frequency to the cell's unperturbed firing rate. Their spike times encoded the frequency and phase of that single oscillatory input.NEW & NOTEWORTHY During repetitive firing, the timing of action potentials is determined by the interaction between the input and voltage-sensitive currents throughout the interspike interval. This interaction is encapsulated in the neuron's phase-resetting curve. The phase-resetting curve predicted spike timing to small sinusoidal currents over a wide range of stimulus frequencies. Firing patterns were most sensitive to oscillatory components near the cell's own firing rate, even in the presence of noise and other inputs.

Qualitative changes in phase-response curve and synchronization at the saddle-node-loop bifurcation

Prominent changes in neuronal dynamics have previously been attributed to a specific switch in onset bifurcation, the Bogdanov-Takens (BT) point. This study unveils another, relevant and so far underestimated transition point: the saddle-node-loop bifurcation, which can be reached by several parameters, including capacitance, leak conductance, and temperature. This bifurcation turns out to induce even more drastic changes in synchronization than the BT transition. This result arises from a direct effect of the saddle-node-loop bifurcation on the limit cycle and hence spike dynamics. In contrast, the BT bifurcation exerts its immediate influence upon the subthreshold dynamics and hence only indirectly relates to spiking. We specifically demonstrate that the saddle-node-loop bifurcation (i) ubiquitously occurs in planar neuron models with a saddle node on invariant cycle onset bifurcation, and (ii) results in a symmetry breaking of the system's phase-response curve. The latter entails an increase in synchronization range in pulse-coupled oscillators, such as neurons. The derived bifurcation structure is of interest in any system for which a relaxation limit is admissible, such as Josephson junctions and chemical oscillators.

Resonant Interneurons Can Increase Robustness of Gamma Oscillations

Gamma oscillations are believed to play a critical role in in information processing, encoding, and retrieval. Inhibitory interneuronal network gamma (ING) oscillations may arise from a coupled oscillator mechanism in which individual neurons oscillate or from a population oscillator in which individual neurons fire sparsely and stochastically. All ING mechanisms, including the one proposed herein, rely on alternating waves of inhibition and windows of opportunity for spiking. The coupled oscillator model implemented with Wang-Buzsáki model neurons is not sufficiently robust to heterogeneity in excitatory drive, and therefore intrinsic frequency, to account for in vitro models of ING. Similarly, in a tightly synchronized regime, the stochastic population oscillator model is often characterized by sparse firing, whereas interneurons both in vivo and in vitro do not fire sparsely during gamma, but rather on average every other cycle. We substituted so-called resonator neural models, which exhibit class 2 excitability and postinhibitory rebound (PIR), for the integrators that are typically used. This results in much greater robustness to heterogeneity that actually increases as the average participation in spikes per cycle approximates physiological levels. Moreover, dynamic clamp experiments that show autapse-induced firing in entorhinal cortical interneurons support the idea that PIR can serve as a network gamma mechanism. Furthermore, parvalbumin-positive (PV(+)) cells were much more likely to display both PIR and autapse-induced firing than GAD2(+) cells, supporting the view that PV(+) fast-firing basket cells are more likely to exhibit class 2 excitability than other types of inhibitory interneurons.

SIGNIFICANCE STATEMENT

Gamma oscillations are believed to play a critical role in information processing, encoding, and retrieval. Networks of inhibitory interneurons are thought to be essential for these oscillations. We show that one class of interneurons with an abrupt onset of firing at a threshold frequency may allow more robust synchronization in the presence of noise and heterogeneity. The mechanism for this robustness depends on the intrinsic resonance at this threshold frequency. Moreover, we show experimentally the feasibility of the proposed mechanism and suggest a way to distinguish between this mechanism and another proposed mechanism: that of a stochastic population oscillator independent of the dynamics of individual neurons.

Anomalous neuronal responses to fluctuated inputs

The irregular firing of a cortical neuron is thought to result from a highly fluctuating drive that is generated by the balance of excitatory and inhibitory synaptic inputs. A previous study reported anomalous responses of the Hodgkin-Huxley neuron to the fluctuated inputs where an irregularity of spike trains is inversely proportional to an input irregularity. In the current study, we investigated the origin of these anomalous responses with the Hindmarsh-Rose neuron model, map-based models, and a simple mixture of interspike interval distributions. First, we specified the parameter regions for the bifurcations in the Hindmarsh-Rose model, and we confirmed that the model reproduced the anomalous responses in the dynamics of the saddle-node and subcritical Hopf bifurcations. For both bifurcations, the Hindmarsh-Rose model shows bistability in the resting state and the repetitive firing state, which indicated that the bistability was the origin of the anomalous input-output relationship. Similarly, the map-based model that contained bistability reproduced the anomalous responses, while the model without bistability did not. These results were supported by additional findings that the anomalous responses were reproduced by mimicking the bistable firing with a mixture of two different interspike interval distributions. Decorrelation of spike trains is important for neural information processing. For such spike train decorrelation, irregular firing is key. Our results indicated that irregular firing can emerge from fluctuating drives, even weak ones, under conditions involving bistability. The anomalous responses, therefore, contribute to efficient processing in the brain.

Firing-rate resonances in the peripheral auditory system of the cricket, Gryllus bimaculatus

In many communication systems, information is encoded in the temporal pattern of signals. For rhythmic signals that carry information in specific frequency bands, a neuronal system may profit from tuning its inherent filtering properties towards a peak sensitivity in the respective frequency range. The cricket Gryllus bimaculatus evaluates acoustic communication signals of both conspecifics and predators. The song signals of conspecifics exhibit a characteristic pulse pattern that contains only a narrow range of modulation frequencies. We examined individual neurons (AN1, AN2, ON1) in the peripheral auditory system of the cricket for tuning towards specific modulation frequencies by assessing their firing-rate resonance. Acoustic stimuli with a swept-frequency envelope allowed an efficient characterization of the cells' modulation transfer functions. Some of the examined cells exhibited tuned band-pass properties. Using simple computational models, we demonstrate how different, cell-intrinsic or network-based mechanisms such as subthreshold resonances, spike-triggered adaptation, as well as an interplay of excitation and inhibition can account for the experimentally observed firing-rate resonances. Therefore, basic neuronal mechanisms that share negative feedback as a common theme may contribute to selectivity in the peripheral auditory pathway of crickets that is designed towards mate recognition and predator avoidance.

Subthreshold amplitude and phase resonance in models of quadratic type: nonlinear effects generated by the interplay of resonant and amplifying currents

We investigate the biophysical and dynamic mechanisms of generation of subthreshold amplitude and phase resonance in response to sinusoidal input currents in two-dimensional models of quadratic type. These models feature a parabolic voltage nullcline and a linear nullcline for the recovery gating variable, capturing the interplay of the so-called resonant currents (e.g., hyperpolarization-activated mixed-cation inward and slow potassium) and amplifying currents (e.g., persistent sodium) in biophysically realistic parameter regimes. These currents underlie the generation of resonance in medial entorhinal cortex layer II stellate cells and CA1 pyramidal cells. We show that quadratic models exhibit nonlinear amplifications of the voltage response to sinusoidal inputs in the resonant frequency band. These are expressed as an increase in the impedance profile as the input amplitude increases. They are stronger for values positive than negative to resting potential and are accompanied by a shift in the phase profile, a decrease in the resonant and phase-resonant frequencies, and an increase in the sharpness of the voltage response. These effects are more prominent for smaller values of ∊ (larger levels of the time scale separation between the voltage and the resonant gating variable) and for values of the resting potential closer to threshold for spike generation. All other parameter fixed, as ∊ increases the voltage response becomes "more linear"; i.e., the nonlinearities are present, but "ignored". In addition, the nonlinear effects are strongly modulated by the curvature of the parabolic voltage nullcline (partially reflecting the effects of the amplifying current) and the slope of the resonant current activation curve. Following the effects of changes in the biophysical conductances of realistic conductance-based models through the parameters of the quadratic model, we characterize the qualitatively different effects that resonant and amplifying currents have on the nonlinear properties of the voltage response. We identify different classes of resonant currents, represented by h- and slow potassium, according to whether they enhance (h-) or attenuate (slow potassium) the nonlinear effects. Finally, we use dynamical systems tools to investigate the dynamic mechanisms of generation of resonance and phase-resonance. We show that the nonlinear effects on the voltage response (e.g., amplification of the voltage response in the resonant frequency band and shifts in the resonant and phase-resonant frequencies) result from the ability of limit cycle trajectories to follow the unstable (right) branch of the voltage nullcline for a significant amount of time. This is a canard-related mechanism that has been shown to underlie the generation of intrinsic subthreshold oscillations in quadratic type models such as medial entorhinal cortex stellate cells. Overall, our results highlight the complexity of the voltage response to oscillatory inputs in nonlinear models and the roles that resonant and amplifying currents have in shaping these responses.

Dynamics and Reliability of Bistable Neurons Driven with Time-Dependent Stimuli

The reliability of a spiking neuron depends on the frequency content of the driving input signal. Previous studies have shown that well above threshold, regularly firing neurons generate reliable responses when the input signal resonates with the firing frequency of the cell. Instead, well below threshold, reliable responses are obtained when the input frequency resonates with the subthreshold oscillations of the neuron. Previous theories, however, provide no clear prediction for the input frequency giving rise to maximally reliable spiking at threshold, which is probably the most relevant firing regime in mammalian cortex under physiological conditions. In particular, when the firing onset is governed by a subcritical Hopf bifurcation, the frequency of subthreshold oscillations often differs from the firing rate at threshold. The predictions of previous studies, hence, cannot be smoothly merged at threshold. Here we explore the behavior of reliability in bistable neurons near threshold using three types of driving stimuli: constant, periodic, and stochastic. We find that the two natural frequencies of the system, associated with the two coexisting attractors, provide a rich variety of possible locking modes with the external signal. Reliability is determined by the sensitivity to noise of each locking mode and by the transition probabilities between modes. Noise increases the amount of spike time jitter, and minimal jitter is obtained for input frequencies coinciding with the suprathreshold firing rate of the cell. In addition, noise may either enhance or inhibit transitions between the two attractors, depending on the input frequency. The dual role played by noise in bistable systems implies that reliability is determined by a delicate balance between spike time jitter and the rate of transitions between attractors.

Neuronal spike timing adaptation described with a fractional leaky integrate-and-fire model

The voltage trace of neuronal activities can follow multiple timescale dynamics that arise from correlated membrane conductances. Such processes can result in power-law behavior in which the membrane voltage cannot be characterized with a single time constant. The emergent effect of these membrane correlations is a non-Markovian process that can be modeled with a fractional derivative. A fractional derivative is a non-local process in which the value of the variable is determined by integrating a temporal weighted voltage trace, also called the memory trace. Here we developed and analyzed a fractional leaky integrate-and-fire model in which the exponent of the fractional derivative can vary from 0 to 1, with 1 representing the normal derivative. As the exponent of the fractional derivative decreases, the weights of the voltage trace increase. Thus, the value of the voltage is increasingly correlated with the trajectory of the voltage in the past. By varying only the fractional exponent, our model can reproduce upward and downward spike adaptations found experimentally in neocortical pyramidal cells and tectal neurons in vitro. The model also produces spikes with longer first-spike latency and high inter-spike variability with power-law distribution. We further analyze spike adaptation and the responses to noisy and oscillatory input. The fractional model generates reliable spike patterns in response to noisy input. Overall, the spiking activity of the fractional leaky integrate-and-fire model deviates from the spiking activity of the Markovian model and reflects the temporal accumulated intrinsic membrane dynamics that affect the response of the neuron to external stimulation.

Spike adaptation is a property of most neurons. When spike time adaptation occurs over multiple time scales, the dynamics can be described by a power-law. We study the computational properties of a leaky integrate-and-fire model with power-law adaptation. Instead of explicitly modeling the adaptation process by the contribution of slowly changing conductances, we use a fractional temporal derivative framework. The exponent of the fractional derivative represents the degree of adaptation of the membrane voltage, where 1 is the normal leaky integrator while values less than 1 produce increasing correlations in the voltage trace. The temporal correlation is interpreted as a memory trace that depends on the value of the fractional derivative. We identify the memory trace in the fractional model as the sum of the instantaneous differentiation weighted by a function that depends on the fractional exponent, and it provides non-local information to the incoming stimulus. The spiking dynamics of the fractional leaky integrate-and-fire model show memory dependence that can result in downward or upward spike adaptation. Our model provides a framework for understanding how long-range membrane voltage correlations affect spiking dynamics and information integration in neurons.

Network models predict that reduced excitatory fluctuations can give rise to hippocampal network hyper-excitability in MeCP2-null mice

Rett syndrome is a severe pediatric neurological disorder caused by loss of function mutations within the gene encoding methyl CpG-binding protein 2 (MeCP2). Although MeCP2 is expressed near ubiquitously, the primary pathophysiology of Rett syndrome stems from impairments of nervous system function. One alteration within different regions of the MeCP2-deficient brain is the presence of hyper-excitable network responses. In the hippocampus, such responses exist despite there being an overall decrease in spontaneous excitatory drive within the network. In this study, we generated and used mathematical, neuronal network models to resolve this apparent paradox. We did this by taking advantage of previous mathematical modelling insights that indicated that decreased excitatory fluctuations, but not mean excitatory drive, more critically explain observed changes in hippocampal network oscillations from MeCP2-null mouse slices. Importantly, reduced excitatory fluctuations could also bring about hyper-excitable responses in our network models. Therefore, these results indicate that diminished excitatory fluctuations may be responsible for the hyper-excitable state of MeCP2-deficient hippocampal circuitry.

Somatic versus dendritic resonance: differential filtering of inputs through non-uniform distributions of active conductances

Synaptic inputs to neurons are processed in a frequency-dependent manner, with either low-pass or resonant response characteristics. These types of filtering play a key role in the frequency-specific information flow in neuronal networks. While the generation of resonance by specific ionic conductances is well investigated, less attention has been paid to the spatial distribution of the resonance-generating conductances across a neuron. In pyramidal neurons – one of the major excitatory cell-types in the mammalian brain – a steep gradient of resonance-generating h-conductances with a 60-fold increase towards distal dendrites has been demonstrated experimentally. Because the dendritic trees of these cells are large, spatial compartmentalization of resonant properties can be expected. Here, we use mathematical descriptions of spatially extended neurons to investigate the consequences of such a distal, dendritic localization of h-conductances for signal processing. While neurons with short dendrites do not exhibit a pronounced compartmentalization of resonance, i.e. the filter properties of dendrites and soma are similar, we find that neurons with longer dendrites ( space constant) can show distinct filtering of dendritic and somatic inputs due to electrotonic segregation. Moreover, we show that for such neurons, experimental classification as resonant versus nonresonant can be misleading when based on somatic recordings, because for these morphologies a dendritic resonance could easily be undetectable when using somatic input. Nevertheless, noise-driven membrane-potential oscillations caused by dendritic resonance can propagate to the soma where they can be recorded, hence contrasting with the low-pass filtering at the soma. We conclude that non-uniform distributions of active conductances can underlie differential filtering of synaptic input in neurons with spatially extended dendrites, like pyramidal neurons, bearing relevance for the localization-dependent targeting of synaptic input pathways to these cells.

A computational study of spike time reliability in two types of threshold dynamics

Spike time reliability (STR) refers to the phenomenon in which repetitive applications of a frozen copy of one stochastic signal to a neuron trigger spikes with reliable timing while a constant signal fails to do so. Observed and explored in numerous experimental and theoretical studies, STR is a complex dynamic phenomenon depending on the nature of external inputs as well as intrinsic properties of a neuron. The neuron under consideration could be either quiescent or spontaneously spiking in the absence of the external stimulus. Focusing on the situation in which the unstimulated neuron is quiescent but close to a switching point to oscillations, we numerically analyze STR treating each spike occurrence as a time localized event in a model neuron. We study both the averaged properties as well as individual features of spike-evoking epochs (SEEs). The effects of interactions between spikes is minimized by selecting signals that generate spikes with relatively long interspike intervals (ISIs). Under these conditions, the frequency content of the input signal has little impact on STR. We study two distinct cases, Type I in which the f-I relation (f for frequency, I for applied current) is continuous and Type II where the f-I relation exhibits a jump. STR in the two types shows a number of similar features and differ in some others. SEEs that are capable of triggering spikes show great variety in amplitude and time profile. On average, reliable spike timing is associated with an accelerated increase in the "action" of the signal as a threshold for spike generation is approached. Here, "action" is defined as the average amount of current delivered during a fixed time interval. When individual SEEs are studied, however, their time profiles are found important for triggering more precisely timed spikes. The SEEs that have a more favorable time profile are capable of triggering spikes with higher precision even at lower action levels.

Reliability of spike and burst firing in thalamocortical relay cells

The reliability and precision of the timing of spikes in a spike train is an important aspect of neuronal coding. We investigated reliability in thalamocortical relay (TCR) cells in the acute slice and also in a Morris-Lecar model with several extensions. A frozen Gaussian noise current, superimposed on a DC current, was injected into the TCR cell soma. The neuron responded with spike trains that showed trial-to-trial variability, due to amongst others slow changes in its internal state and the experimental setup. The DC current allowed to bring the neuron in different states, characterized by a well defined membrane voltage (between -80 and -50 mV) and by a specific firing regime that on depolarization gradually shifted from a predominantly bursting regime to a tonic spiking regime. The filtered frozen white noise generated a spike pattern output with a broad spike interval distribution. The coincidence factor and the Hunter and Milton measure were used as reliability measures of the output spike train. In the experimental TCR cell as well as the Morris-Lecar model cell the reliability depends on the shape (steepness) of the current input versus spike frequency output curve. The model also allowed to study the contribution of three relevant ionic membrane currents to reliability: a T-type calcium current, a cation selective h-current and a calcium dependent potassium current in order to allow bursting, investigate the consequences of a more complex current-frequency relation and produce realistic firing rates. The reliability of the output of the TCR cell increases with depolarization. In hyperpolarized states bursts are more reliable than single spikes. The analytically derived relations were capable to predict several of the experimentally recorded spike features.

Biologically inspired information processing and synchronization in ensembles of non-identical threshold-potential nanostructures

Nanotechnology produces basic structures that show a significant variability in their individual physical properties. This experimental fact may constitute a serious limitation for most applications requiring nominally identical building blocks. On the other hand, biological diversity is found in most natural systems. We show that reliable information processing can be achieved with heterogeneous groups of non-identical nanostructures by using some conceptual schemes characteristic of biological networks (diversity, frequency-based signal processing, rate and rank order coding, and synchronization). To this end, we simulate the integrated response of an ensemble of single-electron transistors (SET) whose individual threshold potentials show a high variability. A particular experimental realization of a SET is a metal nanoparticle-based transistor that mimics biological spiking synapses and can be modeled as an integrate-and-fire oscillator. The different shape and size distributions of nanoparticles inherent to the nanoscale fabrication procedures result in a significant variability in the threshold potentials of the SET. The statistical distributions of the nanoparticle physical parameters are characterized by experimental average and distribution width values. We consider simple but general information processing schemes to draw conclusions that should be of relevance for other threshold-based nanostructures. Monte Carlo simulations show that ensembles of non-identical SET may show some advantages over ensembles of identical nanostructures concerning the processing of weak signals. The results obtained are also relevant for understanding the role of diversity in biophysical networks.

Neurodynamics, tonality, and the auditory brainstem response

Tonal relationships are foundational in music, providing the basis upon which musical structures, such as melodies, are constructed and perceived. A recent dynamic theory of musical tonality predicts that networks of auditory neurons resonate nonlinearly to musical stimuli. Nonlinear resonance leads to stability and attraction relationships among neural frequencies, and these neural dynamics give rise to the perception of relationships among tones that we collectively refer to as tonal cognition. Because this model describes the dynamics of neural populations, it makes specific predictions about human auditory neurophysiology. Here, we show how predictions about the auditory brainstem response (ABR) are derived from the model. To illustrate, we derive a prediction about population responses to musical intervals that has been observed in the human brainstem. Our modeled ABR shows qualitative agreement with important features of the human ABR. This provides a source of evidence that fundamental principles of auditory neurodynamics might underlie the perception of tonal relationships, and forces reevaluation of the role of learning and enculturation in tonal cognition.

Fluctuating inhibitory inputs promote reliable spiking at theta frequencies in hippocampal interneurons

Theta-frequency (4–12 Hz) rhythms in the hippocampus play important roles in learning and memory. CA1 interneurons located at the stratum lacunosum-moleculare and radiatum junction (LM/RAD) are thought to contribute to hippocampal theta population activities by rhythmically pacing pyramidal cells with inhibitory postsynaptic potentials. This implies that LM/RAD cells need to fire reliably at theta frequencies in vivo. To determine whether this could occur, we use biophysically based LM/RAD model cells and apply different cholinergic and synaptic inputs to simulate in vivo-like network environments. We assess spike reliabilities and spiking frequencies, identifying biophysical properties and network conditions that best promote reliable theta spiking. We find that synaptic background activities that feature large inhibitory, but not excitatory, fluctuations are essential. This suggests that strong inhibitory input to these cells is vital for them to be able to contribute to population theta activities. Furthermore, we find that Type I-like oscillator models produced by augmented persistent sodium currents (I_NaP) or diminished A-type potassium currents (I_A) enhance reliable spiking at lower theta frequencies. These Type I-like models are also the most responsive to large inhibitory fluctuations and can fire more reliably under such conditions. In previous work, we showed that I_NaP and I_A are largely responsible for establishing LM/RAD cells’ subthreshold activities. Taken together with this study, we see that while both these currents are important for subthreshold theta fluctuations and reliable theta spiking, they contribute in different ways – I_NaP to reliable theta spiking and subthreshold activity generation, and I_A to subthreshold activities at theta frequencies. This suggests that linking subthreshold and suprathreshold activities should be done with consideration of both in vivo contexts and biophysical specifics.

GAD67-GFP+ neurons in the Nucleus of Roller. II. Subthreshold and firing resonance properties

In the companion paper we show that GAD67-GFP+ (GFP+) inhibitory neurons located in the Nucleus of Roller of the mouse brain stem can be classified into two main groups (tonic and phasic) based on their firing patterns in responses to injected depolarizing current steps. In this study we examined the responses of GFP+ cells to fluctuating sinusoidal ("chirp") current stimuli. Membrane impedance profiles in response to chirp stimulation showed that nearly all phasic cells exhibited subthreshold resonance, whereas the majority of tonic GFP+ cells were nonresonant. In general, subthreshold resonance was associated with a relatively fast passive membrane time constant and low input resistance. In response to suprathreshold chirp current stimulation at a holding potential just below spike threshold the majority of tonic GFP+ cells fired multiple action potentials per cycle at low input frequencies (<5 Hz) and either stopped firing or were not entrained by the chirp at higher input frequencies (= tonic low-pass cells). A smaller group of phasic GFP+ cells did not fire at low input frequency but were able to phase-lock 1:1 at intermediate chirp frequencies (= band-pass cells). Spike timing reliability was tested with repeated chirp stimuli and our results show that phasic cells were able to reliably fire when they phase-locked 1:1 over a relatively broad range of input frequencies. Most tonic low-pass cells showed low reliability and poor phase-locking ability. Computer modeling suggested that these different firing resonance properties among GFP+ cells are due to differences in passive and active membrane properties and spiking mechanisms. This heterogeneity of resonance properties might serve to selectively activate subgroups of interneurons.

Rate maintenance and resonance in the entorhinal cortex

Throughout the brain, neurons encode information in fundamental units of spikes. Each spike represents the combined thresholding of synaptic inputs and intrinsic neuronal dynamics. Here, we address a basic question of spike train formation: how do perithreshold synaptic inputs perturb the output of a spiking neuron? We recorded from single entorhinal principal cells in vitro and drove them to spike steadily at ∼5 Hz (theta range) with direct current injection, then used a dynamic-clamp to superimpose strong excitatory conductance inputs at varying rates. Neurons spiked most reliably when the input rate matched the intrinsic neuronal firing rate. We also found a striking tendency of neurons to preserve their rates and coefficients of variation, independently of input rates. As mechanisms for this rate maintenance, we show that the efficacy of the conductance inputs varied with the relationship of input rate to neuronal firing rate, and with the arrival time of the input within the natural period. Using a novel method of spike classification, we developed a minimal Markov model that reproduced the measured statistics of the output spike trains and thus allowed us to identify and compare contributions to the rate maintenance and resonance. We suggest that the strength of rate maintenance may be used as a new categorization scheme for neuronal response and note that individual intrinsic spiking mechanisms may play a significant role in forming the rhythmic spike trains of activated neurons; in the entorhinal cortex, individual pacemakers may dominate production of the regional theta rhythm.

Characterization of voltage-gated K+ currents contributing to subthreshold membrane potential oscillations in hippocampal CA1 interneurons

CA1 inhibitory interneurons at the stratum lacunosum-moleculare and radiatum junction (LM/RAD-INs) display subthreshold membrane potential oscillations (MPOs) involving voltage-dependent Na(+) and A-type K(+) currents. LM/RAD-INs also express other voltage-gated K(+) currents, although their properties and role in MPOs remain unclear. Here, we characterized these voltage-gated K(+) currents and investigated their role in MPOs. Using outside-out patch recordings from LM/RAD-IN somata, we distinguished four voltage-gated K(+) currents based on their pharmacology and activation/inactivation properties: a fast delayed rectifier current (I(Kfast)), a slow delayed rectifier current (I(Kslow)), a rapidly inactivating A-type current (I(A)), and a slowly inactivating current (I(D)). Their relative contribution to the total K(+) current was I(A) > I(Kfast) > I(Kslow) = I(D). The presence of I(D) and the relative contributions of K(+) currents in LM/RAD-INs are different from those of other CA1 interneurons, suggesting the presence of differential complement of K(+) currents in subgroups of interneurons. We next determined whether these K(+) currents were sufficient for MPO generation using a single-compartment model of LM/RAD-INs. The model captured the subthreshold voltage dependence of MPOs. Moreover, all K(+) currents were active at subthreshold potentials but I(D), I(A), and the persistent sodium current (I(NaP)) were most active near threshold. Using impedance analysis, we found that I(A) and I(NaP) contribute to MPO generation by modulating peak spectral frequency during MPOs and governing the voltage range over which MPOs occur. Our findings uncover a differential expression of a complement of K(+) channels that underlies intrinsic rhythmic activity in inhibitory interneurons.

Mode-locked spike trains in responses of ventral cochlear nucleus chopper and onset neurons to periodic stimuli

We report evidence of mode-locking to the envelope of a periodic stimulus in chopper units of the ventral cochlear nucleus (VCN). Mode-locking is a generalized description of how responses in periodically forced nonlinear systems can be closely linked to the input envelope, while showing temporal patterns of higher order than seen during pure phase-locking. Re-analyzing a previously unpublished dataset in response to amplitude modulated tones, we find that of 55% of cells (6/11) demonstrated stochastic mode-locking in response to sinusoidally amplitude modulated (SAM) pure tones at 50% modulation depth. At 100% modulation depth SAM, most units (3/4) showed mode-locking. We use interspike interval (ISI) scattergrams to unravel the temporal structure present in chopper mode-locked responses. These responses compared well to a leaky integrate-and-fire model (LIF) model of chopper units. Thus the timing of spikes in chopper unit responses to periodic stimuli can be understood in terms of the complex dynamics of periodically forced nonlinear systems. A larger set of onset (33) and chopper units (24) of the VCN also shows mode-locked responses to steady-state vowels and cosine-phase harmonic complexes. However, while 80% of chopper responses to complex stimuli meet our criterion for the presence of mode-locking, only 40% of onset cells show similar complex-modes of spike patterns. We found a correlation between a unit's regularity and its tendency to display mode-locked spike trains as well as a correlation in the number of spikes per cycle and the presence of complex-modes of spike patterns. These spiking patterns are sensitive to the envelope as well as the fundamental frequency of complex sounds, suggesting that complex cell dynamics may play a role in encoding periodic stimuli and envelopes in the VCN.