Threshold fatigue and information transfer

A bifurcation integrates information from many noisy ion channels and allows for milli-Kelvin thermal sensitivity in the snake pit organ

In various biological systems, information from many noisy molecular receptors must be integrated into a collective response. A striking example is the thermal imaging organ of pit vipers. Single nerve fibers in the organ reliably respond to milli-Kelvin (mK) temperature increases, a thousand times more sensitive than their molecular sensors, thermo-transient receptor potential (TRP) ion channels. Here, we propose a mechanism for the integration of this molecular information. In our model, amplification arises due to proximity to a dynamical bifurcation, separating a regime with frequent and regular action potentials (APs), from a regime where APs are irregular and infrequent. Near the transition, AP frequency can have an extremely sharp dependence on temperature, naturally accounting for the thousand-fold amplification. Furthermore, close to the bifurcation, most of the information about temperature available in the TRP channels' kinetics can be read out from the times between consecutive APs even in the presence of readout noise. A key model prediction is that the coefficient of variation in the distribution of interspike times decreases with AP frequency, and quantitative comparison with experiments indeed suggests that nerve fibers of snakes are located very close to the bifurcation. While proximity to such bifurcation points typically requires fine-tuning of parameters, we propose that having feedback act from the order parameter (AP frequency) onto the control parameter robustly maintains the system in the vicinity of the bifurcation. This robustness suggests that similar feedback mechanisms might be found in other sensory systems which also need to detect tiny signals in a varying environment.

Short-term neuronal and synaptic plasticity act in synergy for deviance detection in spiking networks

Sensory areas of cortex respond more strongly to infrequent stimuli when these violate previously established regularities, a phenomenon known as deviance detection (DD). Previous modeling work has mainly attempted to explain DD on the basis of synaptic plasticity. However, a large fraction of cortical neurons also exhibit firing rate adaptation, an underexplored potential mechanism. Here, we investigate DD in a spiking neuronal network model with two types of short-term plasticity, fast synaptic short-term depression (STD) and slower threshold adaptation (TA). We probe the model with an oddball stimulation paradigm and assess DD by evaluating the network responses. We find that TA is sufficient to elicit DD. It achieves this by habituating neurons near the stimulation site that respond earliest to the frequently presented standard stimulus (local fatigue), which diminishes the response and promotes the recovery (global fatigue) of the wider network. Further, we find a synergy effect between STD and TA, where they interact with each other to achieve greater DD than the sum of their individual effects. We show that this synergy is caused by the local fatigue added by STD, which inhibits the global response to the frequently presented stimulus, allowing greater recovery of TA-mediated global fatigue and making the network more responsive to the deviant stimulus. Finally, we show that the magnitude of DD strongly depends on the timescale of stimulation. We conclude that highly predictable information can be encoded in strong local fatigue, which allows greater global recovery and subsequent heightened sensitivity for DD.

A dynamic spike threshold with correlated noise predicts observed patterns of negative interval correlations in neuronal spike trains

Negative correlations in the sequential evolution of interspike intervals (ISIs) are a signature of memory in neuronal spike-trains. They provide coding benefits including firing-rate stabilization, improved detectability of weak sensory signals, and enhanced transmission of information by improving signal-to-noise ratio. Primary electrosensory afferent spike-trains in weakly electric fish fall into two categories based on the pattern of ISI correlations: non-bursting units have negative correlations which remain negative but decay to zero with increasing lags (Type I ISI correlations), and bursting units have oscillatory (alternating sign) correlation which damp to zero with increasing lags (Type II ISI correlations). Here, we predict and match observed ISI correlations in these afferents using a stochastic dynamic threshold model. We determine the ISI correlation function as a function of an arbitrary discrete noise correlation function [Formula: see text], where k is a multiple of the mean ISI. The function permits forward and inverse calculations of the correlation function. Both types of correlation functions can be generated by adding colored noise to the spike threshold with Type I correlations generated with slow noise and Type II correlations generated with fast noise. A first-order autoregressive (AR) process with a single parameter is sufficient to predict and accurately match both types of afferent ISI correlation functions, with the type being determined by the sign of the AR parameter. The predicted and experimentally observed correlations are in geometric progression. The theory predicts that the limiting sum of ISI correlations is [Formula: see text] yielding a perfect DC-block in the power spectrum of the spike train. Observed ISI correlations from afferents have a limiting sum that is slightly larger at [Formula: see text] ([Formula: see text]). We conclude that the underlying process for generating ISIs may be a simple combination of low-order AR and moving average processes and discuss the results from the perspective of optimal coding.

Cellular and Network Mechanisms May Generate Sparse Coding of Sequential Object Encounters in Hippocampal-Like Circuits

The localization of distinct landmarks plays a crucial role in encoding new spatial memories. In mammals, this function is performed by hippocampal neurons that sparsely encode an animal’s location relative to surrounding objects. Similarly, the dorsolateral pallium (DL) is essential for spatial learning in teleost fish. The DL of weakly electric gymnotiform fish receives both electrosensory and visual input from the preglomerular nucleus (PG), which has been hypothesized to encode the temporal sequence of electrosensory or visual landmark/food encounters. Here, we show that DL neurons in the Apteronotid fish and in the Carassius auratus (goldfish) have a hyperpolarized resting membrane potential (RMP) combined with a high and dynamic spike threshold that increases following each spike. Current-evoked spikes in DL cells are followed by a strong small-conductance calcium-activated potassium channel (SK)-mediated after-hyperpolarizing potential (AHP). Together, these properties prevent high frequency and continuous spiking. The resulting sparseness of discharge and dynamic threshold suggest that DL neurons meet theoretical requirements for generating spatial memory engrams by decoding the landmark/food encounter sequences encoded by PG neurons. Thus, DL neurons in teleost fish may provide a promising, simple system to study the core cell and network mechanisms underlying spatial memory.

Statistics of inverse interspike intervals: The instantaneous firing rate revisited

The rate coding hypothesis is the oldest and still one of the most accepted and investigated scenarios in neuronal activity analyses. However, the actual neuronal firing rate, while informally understood, can be mathematically defined in several different ways. These definitions yield distinct results; even their average values may differ dramatically for the simplest neuronal models. Such an inconsistency, together with the importance of "firing rate," motivates us to revisit the classical concept of the instantaneous firing rate. We confirm that different notions of firing rate can in fact be compatible, at least in terms of their averages, by carefully discerning the time instant at which the neuronal activity is observed. Two general cases are distinguished: either the inspection time is synchronised with a reference time or with the neuronal spiking. The statistical properties of the instantaneous firing rate, including parameter estimation, are analyzed, and compatibility with the intuitively understood concept is demonstrated.

The effect of positive interspike interval correlations on neuronal information transmission

Experimentally it is known that some neurons encode preferentially information about low-frequency (slow) components of a time-dependent stimulus while others prefer intermediate or high-frequency (fast) components. Accordingly, neurons can be categorized as low-pass, band-pass or high-pass information filters. Mechanisms of information filtering at the cellular and the network levels have been suggested. Here we propose yet another mechanism, based on noise shaping due to spontaneous non-renewal spiking statistics. We compare two integrate-and-fire models with threshold noise that differ solely in their interspike interval (ISI) correlations: the renewal model generates independent ISIs, whereas the non-renewal model exhibits positive correlations between adjacent ISIs. For these simplified neuron models we analytically calculate ISI density and power spectrum of the spontaneous spike train as well as approximations for input-output cross-spectrum and spike-train power spectrum in the presence of a broad-band Gaussian stimulus. This yields the spectral coherence, an approximate frequency-resolved measure of information transmission. We demonstrate that for low spiking variability the renewal model acts as a low-pass filter of information (coherence has a global maximum at zero frequency), whereas the non-renewal model displays a pronounced maximum of the coherence at non-vanishing frequency and thus can be regarded as a band-pass filter of information.

Impact of slow K(+) currents on spike generation can be described by an adaptive threshold model

A neuron that is stimulated by rectangular current injections initially responds with a high firing rate, followed by a decrease in the firing rate. This phenomenon is called spike-frequency adaptation and is usually mediated by slow K(+) currents, such as the M-type K(+) current (I M ) or the Ca(2+)-activated K(+) current (I AHP ). It is not clear how the detailed biophysical mechanisms regulate spike generation in a cortical neuron. In this study, we investigated the impact of slow K(+) currents on spike generation mechanism by reducing a detailed conductance-based neuron model. We showed that the detailed model can be reduced to a multi-timescale adaptive threshold model, and derived the formulae that describe the relationship between slow K(+) current parameters and reduced model parameters. Our analysis of the reduced model suggests that slow K(+) currents have a differential effect on the noise tolerance in neural coding.

Enhanced Sensitivity to Rapid Input Fluctuations by Nonlinear Threshold Dynamics in Neocortical Pyramidal Neurons

The way in which single neurons transform input into output spike trains has fundamental consequences for network coding. Theories and modeling studies based on standard Integrate-and-Fire models implicitly assume that, in response to increasingly strong inputs, neurons modify their coding strategy by progressively reducing their selective sensitivity to rapid input fluctuations. Combining mathematical modeling with in vitro experiments, we demonstrate that, in L5 pyramidal neurons, the firing threshold dynamics adaptively adjust the effective timescale of somatic integration in order to preserve sensitivity to rapid signals over a broad range of input statistics. For that, a new Generalized Integrate-and-Fire model featuring nonlinear firing threshold dynamics and conductance-based adaptation is introduced that outperforms state-of-the-art neuron models in predicting the spiking activity of neurons responding to a variety of in vivo-like fluctuating currents. Our model allows for efficient parameter extraction and can be analytically mapped to a Generalized Linear Model in which both the input filter—describing somatic integration—and the spike-history filter—accounting for spike-frequency adaptation—dynamically adapt to the input statistics, as experimentally observed. Overall, our results provide new insights on the computational role of different biophysical processes known to underlie adaptive coding in single neurons and support previous theoretical findings indicating that the nonlinear dynamics of the firing threshold due to Na⁺-channel inactivation regulate the sensitivity to rapid input fluctuations.

Over the last decades, a variety of simplified spiking models have been shown to achieve a surprisingly high performance in predicting the neuronal responses to in vitro somatic current injections. Because of the complex adaptive behavior featured by cortical neurons, this success is however restricted to limited stimulus ranges: model parameters optimized for a specific input regime are often inappropriate to describe the response to input currents with different statistical properties. In the present study, a new spiking neuron model is introduced that captures single-neuron computation over a wide range of input statistics and explains different aspects of the neuronal dynamics within a single framework. Our results indicate that complex forms of single neuron adaptation are mediated by the nonlinear dynamics of the firing threshold and that the input-output transformation performed by cortical pyramidal neurons can be intuitively understood in terms of an enhanced Generalized Linear Model in which both the input filter and the spike-history filter adapt to the input statistics.

Balanced ionotropic receptor dynamics support signal estimation via voltage-dependent membrane noise

Encoding behaviorally relevant stimuli in a noisy background is critical for animals to survive in their natural environment. We identify core biophysical and synaptic mechanisms that permit the encoding of low-frequency signals in pyramidal neurons of the weakly electric fish Apteronotus leptorhynchus, an animal that can accurately encode even miniscule amplitude modulations of its self-generated electric field. We demonstrate that slow NMDA receptor (NMDA-R)-mediated excitatory postsynaptic potentials (EPSPs) are able to summate over many interspike intervals (ISIs) of the primary electrosensory afferents (EAs), effectively eliminating the baseline EA ISI correlations from the pyramidal cell input. Together with a dynamic balance of NMDA-R and GABA-A-R currents, this permits stimulus-evoked changes in EA spiking to be transmitted efficiently to target electrosensory lobe (ELL) pyramidal cells, for encoding low-frequency signals. Interestingly, AMPA-R activity is depressed and appears to play a negligible role in the generation of action potentials. Instead, we hypothesize that cell-intrinsic voltage-dependent membrane noise supports the encoding of perithreshold sensory input; this noise drives a significant proportion of pyramidal cell spikes. Together, these mechanisms may be sufficient for the ELL to encode signals near the threshold of behavioral detection.

The interplay between cooperativity and diversity in model threshold ensembles

The interplay between cooperativity and diversity is crucial for biological ensembles because single molecule experiments show a significant degree of heterogeneity and also for artificial nanostructures because of the high individual variability characteristic of nanoscale units. We study the cross-effects between cooperativity and diversity in model threshold ensembles composed of individually different units that show a cooperative behaviour. The units are modelled as statistical distributions of parameters (the individual threshold potentials here) characterized by central and width distribution values. The simulations show that the interplay between cooperativity and diversity results in ensemble-averaged responses of interest for the understanding of electrical transduction in cell membranes, the experimental characterization of heterogeneous groups of biomolecules and the development of biologically inspired engineering designs with individually different building blocks.

Optimal decoding and information transmission in Hodgkin-Huxley neurons under metabolic cost constraints

Information theory quantifies the ultimate limits on reliable information transfer by means of the channel capacity. However, the channel capacity is known to be an asymptotic quantity, assuming unlimited metabolic cost and computational power. We investigate a single-compartment Hodgkin-Huxley type neuronal model under the spike-rate coding scheme and address how the metabolic cost and the decoding complexity affects the optimal information transmission. We find that the sub-threshold stimulation regime, although attaining the smallest capacity, allows for the most efficient balance between the information transmission and the metabolic cost. Furthermore, we determine post-synaptic firing rate histograms that are optimal from the information-theoretic point of view, which enables the comparison of our results with experimental data.

Colored noise and memory effects on formal spiking neuron models

Simplified neuronal models capture the essence of the electrical activity of a generic neuron, besides being more interesting from the computational point of view when compared to higher-dimensional models such as the Hodgkin-Huxley one. In this work, we propose a generalized resonate-and-fire model described by a generalized Langevin equation that takes into account memory effects and colored noise. We perform a comprehensive numerical analysis to study the dynamics and the point process statistics of the proposed model, highlighting interesting new features such as (i) nonmonotonic behavior (emergence of peak structures, enhanced by the choice of colored noise characteristic time scale) of the coefficient of variation (CV) as a function of memory characteristic time scale, (ii) colored noise-induced shift in the CV, and (iii) emergence and suppression of multimodality in the interspike interval (ISI) distribution due to memory-induced subthreshold oscillations. Moreover, in the noise-induced spike regime, we study how memory and colored noise affect the coherence resonance (CR) phenomenon. We found that for sufficiently long memory, not only is CR suppressed but also the minimum of the CV-versus-noise intensity curve that characterizes the presence of CR may be replaced by a maximum. The aforementioned features allow to interpret the interplay between memory and colored noise as an effective control mechanism to neuronal variability. Since both variability and nontrivial temporal patterns in the ISI distribution are ubiquitous in biological cells, we hope the present model can be useful in modeling real aspects of neurons.

First-passage times in integrate-and-fire neurons with stochastic thresholds

We consider a leaky integrate-and-fire neuron with deterministic subthreshold dynamics and a firing threshold that evolves as an Ornstein-Uhlenbeck process. The formulation of this minimal model is motivated by the experimentally observed widespread variation of neural firing thresholds. We show numerically that the mean first-passage time can depend nonmonotonically on the noise amplitude. For sufficiently large values of the correlation time of the stochastic threshold the mean first-passage time is maximal for nonvanishing noise. We provide an explanation for this effect by analytically transforming the original model into a first-passage-time problem for Brownian motion. This transformation also allows for a perturbative calculation of the first-passage-time histograms. In turn this provides quantitative insights into the mechanisms that lead to the nonmonotonic behavior of the mean first-passage time. The perturbation expansion is in excellent agreement with direct numerical simulations. The approach developed here can be applied to any deterministic subthreshold dynamics and any Gauss-Markov processes for the firing threshold. This opens up the possibility to incorporate biophysically detailed components into the subthreshold dynamics, rendering our approach a powerful framework that sits between traditional integrate-and-fire models and complex mechanistic descriptions of neural dynamics.

A frequency-resolved mutual information rate and its application to neural systems

The encoding and processing of time-dependent signals into sequences of action potentials of sensory neurons is still a challenging theoretical problem. Although, with some effort, it is possible to quantify the flow of information in the model-free framework of Shannon's information theory, this yields just a single number, the mutual information rate. This rate does not indicate which aspects of the stimulus are encoded. Several studies have identified mechanisms at the cellular and network level leading to low- or high-pass filtering of information, i.e., the selective coding of slow or fast stimulus components. However, these findings rely on an approximation, specifically, on the qualitative behavior of the coherence function, an approximate frequency-resolved measure of information flow, whose quality is generally unknown. Here, we develop an assumption-free method to measure a frequency-resolved information rate about a time-dependent Gaussian stimulus. We demonstrate its application for three paradigmatic descriptions of neural firing: an inhomogeneous Poisson process that carries a signal in its instantaneous firing rate; an integrator neuron (stochastic integrate-and-fire model) driven by a time-dependent stimulus; and the synchronous spikes fired by two commonly driven integrator neurons. In agreement with previous coherence-based estimates, we find that Poisson and integrate-and-fire neurons are broadband and low-pass filters of information, respectively. The band-pass information filtering observed in the coherence of synchronous spikes is confirmed by our frequency-resolved information measure in some but not all parameter configurations. Our results also explicitly show how the response-response coherence can fail as an upper bound on the information rate.

Emergence of band-pass filtering through adaptive spiking in the owl's cochlear nucleus

In the visual, auditory, and electrosensory modalities, stimuli are defined by first- and second-order attributes. The fast time-pressure signal of a sound, a first-order attribute, is important, for instance, in sound localization and pitch perception, while its slow amplitude-modulated envelope, a second-order attribute, can be used for sound recognition. Ascending the auditory pathway from ear to midbrain, neurons increasingly show a preference for the envelope and are most sensitive to particular envelope modulation frequencies, a tuning considered important for encoding sound identity. The level at which this tuning property emerges along the pathway varies across species, and the mechanism of how this occurs is a matter of debate. In this paper, we target the transition between auditory nerve fibers and the cochlear nucleus angularis (NA). While the owl's auditory nerve fibers simultaneously encode the fast and slow attributes of a sound, one synapse further, NA neurons encode the envelope more efficiently than the auditory nerve. Using in vivo and in vitro electrophysiology and computational analysis, we show that a single-cell mechanism inducing spike threshold adaptation can explain the difference in neural filtering between the two areas. We show that spike threshold adaptation can explain the increased selectivity to modulation frequency, as input level increases in NA. These results demonstrate that a spike generation nonlinearity can modulate the tuning to second-order stimulus features, without invoking network or synaptic mechanisms.

Event-triggered feedback in noise-driven phase oscillators

Using a stochastic nonlinear phase oscillator model, we study the effect of event-triggered feedback on the statistics of interevent intervals. Events are associated with the entering of a new cycle. The feedback is modeled by an instantaneous increase (positive feedback) or decrease (negative feedback) of the oscillator frequency whenever an event occurs followed by an exponential decay on a slow time scale. In addition to the known excitable and oscillatory regimes, which are separated by a saddle node on invariant circle bifurcation, positive feedback can lead to bistable dynamics and a change of the system's excitability. The feedback has also a strong effect on noise-induced phenomena like coherence resonance or anticoherence resonance. Both positive and negative feedback can lead to more regular output for particular noise strengths. Finally, we investigate serial correlations in the sequence of interevent intervals that occur due to the additional slow dynamics. We derive approximations for the serial correlation coefficient and show that positive feedback results in extended positive interval correlations, whereas negative feedback yields short-ranging negative correlations. Investigating the interplay of feedback and the nonlinear phase dynamics close to the bifurcation, we find that correlations are most pronounced for optimal feedback strengths.

A simple Markov model of sodium channels with a dynamic threshold

Characteristics of action potential generation are important to understanding brain functioning and, thus, must be understood and modeled. It is still an open question what model can describe concurrently the phenomena of sharp spike shape, the spike threshold variability, and the divisive effect of shunting on the gain of frequency-current dependence. We reproduced these three effects experimentally by patch-clamp recordings in cortical slices, but we failed to simulate them by any of 11 known neuron models, including one- and multi-compartment, with Hodgkin-Huxley and Markov equation-based sodium channel approximations, and those taking into account sodium channel subtype heterogeneity. Basing on our voltage-clamp data characterizing the dependence of sodium channel activation threshold on history of depolarization, we propose a 3-state Markov model with a closed-to-open state transition threshold dependent on slow inactivation. This model reproduces the all three phenomena. As a reduction of this model, a leaky integrate-and-fire model with a dynamic threshold also shows the effect of gain reduction by shunt. These results argue for the mechanism of gain reduction through threshold dynamics determined by the slow inactivation of sodium channels.

Predicting spike timing in highly synchronous auditory neurons at different sound levels

A challenge for sensory systems is to encode natural signals that vary in amplitude by orders of magnitude. The spike trains of neurons in the auditory system must represent the fine temporal structure of sounds despite a tremendous variation in sound level in natural environments. It has been shown in vitro that the transformation from dynamic signals into precise spike trains can be accurately captured by simple integrate-and-fire models. In this work, we show that the in vivo responses of cochlear nucleus bushy cells to sounds across a wide range of levels can be precisely predicted by deterministic integrate-and-fire models with adaptive spike threshold. Our model can predict both the spike timings and the firing rate in response to novel sounds, across a large input level range. A noisy version of the model accounts for the statistical structure of spike trains, including the reliability and temporal precision of responses. Spike threshold adaptation was critical to ensure that predictions remain accurate at different levels. These results confirm that simple integrate-and-fire models provide an accurate phenomenological account of spike train statistics and emphasize the functional relevance of spike threshold adaptation.

Channel noise from both slow adaptation currents and fast currents is required to explain spike-response variability in a sensory neuron

Spike-timing variability has a large effect on neural information processing. However, for many systems little is known about the noise sources causing the spike-response variability. Here we investigate potential sources of spike-response variability in auditory receptor neurons of locusts, a classic insect model system. At low-spike frequencies, our data show negative interspike-interval (ISI) correlations and ISI distributions that match the inverse Gaussian distribution. These findings can be explained by a white-noise source that interacts with an adaptation current. At higher spike frequencies, more strongly peaked distributions and positive ISI correlations appear, as expected from a canonical model of suprathreshold firing driven by temporally correlated (i.e., colored) noise. Simulations of a minimal conductance-based model of the auditory receptor neuron with stochastic ion channels exclude the delayed rectifier as a possible noise source. Our analysis suggests channel noise from an adaptation current and the receptor or sodium current as main sources for the colored and white noise, respectively. By comparing the ISI statistics with generic models, we find strong evidence for two distinct noise sources. Our approach does not involve any dendritic or somatic recordings that may harm the delicate workings of many sensory systems. It could be applied to various other types of neurons, in which channel noise dominates the fluctuations that shape the neuron's spike statistics.

Enhancing the resolution of a sensor via negative correlation: a biologically inspired approach

We demonstrate that a neuronal system, underpinned by "fire-then-reset" dynamics, can display an enhanced resolution R~T(ob)(-1) where T(ob) is the observation time of the measurement; this occurs when the interspike intervals are negatively correlated and T(ob)<Δ/ε, where ε is a parameter characterizing the level of correlation between interspike intervals and Δ is the average interspike interval. We also show that by introducing negative correlations into the time domain response of a nonlinear dynamical sensor it is possible to replicate this enhanced scaling of the resolution. Thus, we demonstrate the potential for designing a novel class of biomimetic sensors that afford improved signal resolution by functionally utilizing negative correlations.

Information filtering by synchronous spikes in a neural population

Information about time-dependent sensory stimuli is encoded by the spike trains of neurons. Here we consider a population of uncoupled but noisy neurons (each subject to some intrinsic noise) that are driven by a common broadband signal. We ask specifically how much information is encoded in the synchronous activity of the population and how this information transfer is distributed with respect to frequency bands. In order to obtain some insight into the mechanism of information filtering effects found previously in the literature, we develop a mathematical framework to calculate the coherence of the synchronous output with the common stimulus for populations of simple neuron models. Within this frame, the synchronous activity is treated as the product of filtered versions of the spike trains of a subset of neurons. We compare our results for the simple cases of (1) a Poisson neuron with a rate modulation and (2) an LIF neuron with intrinsic white current noise and a current stimulus. For the Poisson neuron, formulas are particularly simple but show only a low-pass behavior of the coherence of synchronous activity. For the LIF model, in contrast, the coherence function of the synchronous activity shows a clear peak at high frequencies, comparable to recent experimental findings. We uncover the mechanism for this shift in the maximum of the coherence and discuss some biological implications of our findings.

Bursts and isolated spikes code for opposite movement directions in midbrain electrosensory neurons

Directional selectivity, in which neurons respond strongly to an object moving in a given direction but weakly or not at all to the same object moving in the opposite direction, is a crucial computation that is thought to provide a neural correlate of motion perception. However, directional selectivity has been traditionally quantified by using the full spike train, which does not take into account particular action potential patterns. We investigated how different action potential patterns, namely bursts (i.e. packets of action potentials followed by quiescence) and isolated spikes, contribute to movement direction coding in a mathematical model of midbrain electrosensory neurons. We found that bursts and isolated spikes could be selectively elicited when the same object moved in opposite directions. In particular, it was possible to find parameter values for which our model neuron did not display directional selectivity when the full spike train was considered but displayed strong directional selectivity when bursts or isolated spikes were instead considered. Further analysis of our model revealed that an intrinsic burst mechanism based on subthreshold T-type calcium channels was not required to observe parameter regimes for which bursts and isolated spikes code for opposite movement directions. However, this burst mechanism enhanced the range of parameter values for which such regimes were observed. Experimental recordings from midbrain neurons confirmed our modeling prediction that bursts and isolated spikes can indeed code for opposite movement directions. Finally, we quantified the performance of a plausible neural circuit and found that it could respond more or less selectively to isolated spikes for a wide range of parameter values when compared with an interspike interval threshold. Our results thus show for the first time that different action potential patterns can differentially encode movement and that traditional measures of directional selectivity need to be revised in such cases.

Parallel coding of first- and second-order stimulus attributes by midbrain electrosensory neurons

Natural stimuli often have time-varying first-order (i.e., mean) and second-order (i.e., variance) attributes that each carry critical information for perception and can vary independently over orders of magnitude. Experiments have shown that sensory systems continuously adapt their responses based on changes in each of these attributes. This adaptation creates ambiguity in the neural code as multiple stimuli may elicit the same neural response. While parallel processing of first- and second-order attributes by separate neural pathways is sufficient to remove this ambiguity, the existence of such pathways and the neural circuits that mediate their emergence have not been uncovered to date. We recorded the responses of midbrain electrosensory neurons in the weakly electric fish Apteronotus leptorhynchus to stimuli with first- and second-order attributes that varied independently in time. We found three distinct groups of midbrain neurons: the first group responded to both first- and second-order attributes, the second group responded selectively to first-order attributes, and the last group responded selectively to second-order attributes. In contrast, all afferent hindbrain neurons responded to both first- and second-order attributes. Using computational analyses, we show how inputs from a heterogeneous population of ON- and OFF-type afferent neurons are combined to give rise to response selectivity to either first- or second-order stimulus attributes in midbrain neurons. Our study thus uncovers, for the first time, generic and widely applicable mechanisms by which parallel processing of first- and second-order stimulus attributes emerges in the brain.

Nonlinear dynamic modeling of neuron action potential threshold during synaptically driven broadband intracellular activity

Activity-dependent variation of neuronal thresholds for action potential (AP) generation is one of the key determinants of spike-train temporal-pattern transformations from presynaptic to postsynaptic spike trains. In this study, we model the nonlinear dynamics of the threshold variation during synaptically driven broadband intracellular activity. First, membrane potentials of single CA1 pyramidal cells were recorded under physiologically plausible broadband stimulation conditions. Second, a method was developed to measure AP thresholds from the continuous recordings of membrane potentials. It involves measuring the turning points of APs by analyzing the third-order derivatives of the membrane potentials. Four stimulation paradigms with different temporal patterns were applied to validate this method by comparing the measured AP turning points and the actual AP thresholds estimated with varying stimulation intensities. Results show that the AP turning points provide consistent measurement of the AP thresholds, except for a constant offset. It indicates that 1) the variation of AP turning points represents the nonlinearities of threshold dynamics; and 2) an optimization of the constant offset is required to achieve accurate spike prediction. Third, a nonlinear dynamical third-order Volterra model was built to describe the relations between the threshold dynamics and the AP activities. Results show that the model can predict threshold accurately based on the preceding APs. Finally, the dynamic threshold model was integrated into a previously developed single neuron model and resulted in a 33% improvement in spike prediction.

Coherence depression in stochastic excitable systems with two-frequency forcing

We study the response of two generic neuron models, the leaky integrate-and-fire (LIF) model and the leaky integrate-and-fire model with dynamic threshold (LIFDT) (i.e., with memory) to a stimulus consisting of two sinusoidal drives with incommensurate frequency, an amplitude modulation ("envelope") noise and a relatively weak additive noise. Spectral and coherence analysis of responses to such naturalistic stimuli reveals how the LIFDT model exhibits better correlation between modulation and spike train even in the presence of both noises. However, a resonance-induced synchrony, occurring when the beat frequency between the sinusoids is close to the intrinsic neuronal firing rate, decreases the coherence in the dynamic threshold case. Under suprathreshold conditions, the modulation noise simultaneously decreases the linear spectral coherence between the spikes and the whole stimulus, as well as between spikes and the stimulus envelope. Our study shows that the coefficient of variation of the envelope fluctuations is positively correlated with the degree of coherence depression. As the coherence function quantifies the linear information transmission, our findings indicate that under certain conditions, a transmission loss results when an excitable system with adaptive properties encodes a beat with frequency in the vicinity of its mean firing rate.

Impact of neural noise on a sensory-motor pathway signaling impending collision

Noise is a major concern in circuits processing electrical signals, including neural circuits. There are many factors that influence how noise propagates through neural circuits, and there are few systems in which noise levels have been studied throughout a processing pathway. We recorded intracellularly from multiple stages of a sensory-motor pathway in the locust that detects approaching objects. We found that responses are more variable and that signal-to-noise ratios (SNRs) are lower further from the sensory periphery. SNRs remain low even with the use of stimuli for which the pathway is most selective and for which the neuron representing its final sensory level must integrate many synaptic inputs. Modeling of this neuron shows that variability in the strength of individual synaptic inputs within a large population has little effect on the variability of the spiking output. In contrast, jitter in the timing of individual inputs and spontaneous variability is important for shaping the responses to preferred stimuli. These results suggest that neural noise is inherent to the processing of visual stimuli signaling impending collision and contributes to shaping neural responses along this sensory-motor pathway.

Coding movement direction by burst firing in electrosensory neurons

Directional selectivity, in which neurons respond strongly to an object moving in a given direction ("preferred") but respond weakly or not at all to an object moving in the opposite direction ("null"), is a critical computation achieved in brain circuits. Previous measures of direction selectivity have compared the numbers of action potentials elicited by each direction of movement, but most sensory neurons display patterning, such as bursting, in their spike trains. To examine the contribution of patterned responses to direction selectivity, we recorded from midbrain neurons in weakly electric fish and found that most neurons responded with a combination of both bursts and isolated spikes to moving object stimuli. In these neurons, we separated bursts and isolated spikes using an interspike interval (ISI) threshold. The directional bias of bursts was significantly higher than that of either the full spike train or the isolated spike train. To examine the encoding and decoding of bursts, we built biologically plausible models that examine 1) the upstream mechanisms that generate these spiking patterns and 2) downstream decoders of bursts. Our model of upstream mechanisms uses an interaction between afferent input and subthreshold calcium channels to give rise to burst firing that occurs preferentially for one direction of movement. We tested this model in vivo by application of calcium antagonists, which reduced burst firing and eliminated the differences in direction selectivity between bursts, isolated spikes, and the full spike train. Our model of downstream decoders used strong synaptic facilitation to achieve qualitatively similar results to those obtained using the ISI threshold criterion. This model shows that direction selective information carried by bursts can be decoded by downstream neurons using biophysically plausible mechanisms.

Cell assembly sequences arising from spike threshold adaptation keep track of time in the hippocampus

Hippocampal neurons can display reliable and long-lasting sequences of transient firing patterns, even in the absence of changing external stimuli. We suggest that time-keeping is an important function of these sequences, and propose a network mechanism for their generation. We show that sequences of neuronal assemblies recorded from rat hippocampal CA1 pyramidal cells can reliably predict elapsed time (15-20 s) during wheel running with a precision of 0.5 s. In addition, we demonstrate the generation of multiple reliable, long-lasting sequences in a recurrent network model. These sequences are generated in the presence of noisy, unstructured inputs to the network, mimicking stationary sensory input. Identical initial conditions generate similar sequences, whereas different initial conditions give rise to distinct sequences. The key ingredients responsible for sequence generation in the model are threshold-adaptation and a Mexican-hat-like pattern of connectivity among pyramidal cells. This pattern may arise from recurrent systems such as the hippocampal CA3 region or the entorhinal cortex. We hypothesize that mechanisms that evolved for spatial navigation also support tracking of elapsed time in behaviorally relevant contexts.

Nonrenewal spike train statistics: causes and functional consequences on neural coding

Many neurons display significant patterning in their spike trains (e.g. oscillations, bursting), and there is accumulating evidence that information is contained in these patterns. In many cases, this patterning is caused by intrinsic mechanisms rather than external signals. In this review, we focus on spiking activity that displays nonrenewal statistics (i.e. memory that persists from one firing to the next). Such statistics are seen in both peripheral and central neurons and appear to be ubiquitous in the CNS. We review the principal mechanisms that can give rise to nonrenewal spike train statistics. These are separated into intrinsic mechanisms such as relative refractoriness and network mechanisms such as coupling with delayed inhibitory feedback. Next, we focus on the functional roles for nonrenewal spike train statistics. These can either increase or decrease information transmission. We also focus on how such statistics can give rise to an optimal integration timescale at which spike train variability is minimal and how this might be exploited by sensory systems to maximize the detection of weak signals. We finish by pointing out some interesting future directions for research in this area. In particular, we explore the interesting possibility that synaptic dynamics might be matched with the nonrenewal spiking statistics of presynaptic spike trains in order to further improve information transmission.

Linear versus nonlinear signal transmission in neuron models with adaptation currents or dynamic thresholds

Spike-frequency adaptation is a prominent aspect of neuronal dynamics that shapes a neuron's signal processing properties on timescales ranging from about 10 ms to >1 s. For integrate-and-fire model neurons spike-frequency adaptation is incorporated either as an adaptation current or as a dynamic firing threshold. Whether a physiologically observed adaptation mechanism should be modeled as an adaptation current or a dynamic threshold, however, is not known. Here we show that a dynamic threshold has a divisive effect on the onset f-I curve (the initial maximal firing rate following a step increase in an input current) measured at increasing mean threshold levels, i.e., adaptation states. In contrast, an adaptation current subtractively shifts this f-I curve to higher inputs without affecting its slope. As a consequence, an adaptation current acts essentially linearly, resulting in a high-pass filter component of the neuron's transfer function for current stimuli. With a dynamic threshold, however, the transfer function strongly depends on the input range because of the multiplicative effect on the f-I curves. Simulations of conductance-based spiking models with adaptation currents, such as afterhyperpolarization (AHP)-type, M-type, and sodium-activated potassium currents, do not show the divisive effects of a dynamic threshold, but agree with the properties of integrate-and-fire neurons with adaptation current. Notably, the effects of slow inactivation of sodium currents cannot be reproduced by either model. Our results suggest that, when lateral shifts of the onset f-I curve are seen in response to adapting inputs, adaptation should be modeled with adaptation currents and not with a dynamic threshold. In contrast, when the slope of onset f-I curves depends on the adaptation state, then adaptation should be modeled with a dynamic threshold. Further, the observation of divisively altered onset f-I curves in adapted neurons with notable variability of their spike threshold could hint to yet known biophysical mechanisms directly affecting the threshold.

In vivo conditions influence the coding of stimulus features by bursts of action potentials

The functional role of burst firing (i.e. the firing of packets of action potentials followed by quiescence) in sensory processing is still under debate. Should bursts be considered as unitary events that signal the presence of a particular feature in the sensory environment or is information about stimulus attributes contained within their temporal structure? We compared the coding of stimulus attributes by bursts in vivo and in vitro of electrosensory pyramidal neurons in weakly electric fish by computing correlations between burst and stimulus attributes. Our results show that, while these correlations were strong in magnitude and significant in vitro, they were actually much weaker in magnitude if at all significant in vivo. We used a mathematical model of pyramidal neuron activity in vivo and showed that such a model could reproduce the correlations seen in vitro, thereby suggesting that differences in burst coding were not due to differences in bursting seen in vivo and in vitro. We next tested whether variability in the baseline (i.e. without stimulation) activity of ELL pyramidal neurons could account for these differences. To do so, we injected noise into our model whose intensity was calibrated to mimic baseline activity variability as quantified by the coefficient of variation. We found that this noise caused significant decreases in the magnitude of correlations between burst and stimulus attributes and could account for differences between in vitro and in vivo conditions. We then tested this prediction experimentally by directly injecting noise in vitro through the recording electrode. Our results show that this caused a lowering in magnitude of the correlations between burst and stimulus attributes in vitro and gave rise to values that were quantitatively similar to those seen under in vivo conditions. While it is expected that noise in the form of baseline activity variability will lower correlations between burst and stimulus attributes, our results show that such variability can account for differences seen in vivo. Thus, the high variability seen under in vivo conditions has profound consequences on the coding of information by bursts in ELL pyramidal neurons. In particular, our results support the viewpoint that bursts serve as a detector of particular stimulus features but do not carry detailed information about such features in their structure.

Spontaneous dynamics and response properties of a Hodgkin-Huxley-type neuron model driven by harmonic synaptic noise

We study statistical properties, response dynamics, and information transmission in a Hodgkin-Huxley-type neuron system, modeling peripheral electroreceptors in paddlefish. In addition to sodium and potassium currents, the neuron model includes fast calcium and slow afterhyperpolarization (AHP) potassium currents. The synaptic transmission from sensory epithelium is modeled by a Poission process with a rate modulated by narrow-band noise, mimicking stochastic epithelial oscillations observed experimentally. We study how the interplay of parameters of AHP current and synaptic noise affects the statistics of spontaneous dynamics and response properties of the system. In particular, we confirm predictions made earlier with perfect integrate and fire and phase neuron models that epithelial oscillations enhance stimulus-response coherence and thus information transmission in electroreceptor system. In addition, we consider a strong stimulus regime and show that coherent epithelial oscillations may reduce variability of electroreceptor responses to time-varying stimuli.

Noise Shaping in Neural Populations with Global Delayed Feedback

The interplay between intrinsic and network dynamics has been the focus of many investigations. Here we use a combination of theoretical and numerical approaches to study the effects of delayed global feedback on the information transmission properties of neural networks. Specifically, we compare networks of neurons that display intrinsic interspike interval correlations (nonrenewal) to networks that do not (renewal). We find that excitatory and inhibitory delays can tune information transmission by single neurons but not by the entire network. Most surprisingly, addition of a delay can change the dependence of the information on the coupling strength for renewal neurons and not for nonrenewal neurons. Our results show that intrinsic ISI correlations can have nontrivial interactions with network-induced phenomena.

SK channels gate information processing in vivo by regulating an intrinsic bursting mechanism seen in vitro

Understanding the mechanistic substrates of neural computations that lead to behavior remains a fundamental problem in neuroscience. In particular, the contributions of intrinsic neural properties such as burst firing and dendritic morphology to the processing of behaviorally relevant sensory input have received much interest recently. Pyramidal cells within the electrosensory lateral line lobe of weakly electric fish display an intrinsic bursting mechanism that relies on somato-dendritic interactions when recorded in vitro: backpropagating somatic action potentials trigger dendritic action potentials that lead to a depolarizing afterpotential (DAP) at the soma. We recorded intracellularly from these neurons in vivo and found firing patterns that were quite different from those seen in vitro: we found no evidence for DAPs as each somatic action potential was followed by a pronounced afterhyperpolarization (AHP). Calcium chelators injected in vivo reduced the AHP, thereby unmasking the DAP and inducing in vitro-like bursting in pyramidal cells. These bursting dynamics significantly reduced the cell's ability to encode the detailed time course of sensory input. We performed additional in vivo pharmacological manipulations and mathematical modeling to show that calcium influx through N-methyl-d-aspartate (NMDA) receptors activate dendritic small conductance (SK) calcium-activated potassium channels, which causes an AHP that counteracts the DAP and leads to early termination of the burst. Our results show that ion channels located in dendrites can have a profound influence on the processing of sensory input by neurons in vivo through the modulation of an intrinsic bursting mechanism.

Serial correlation in neural spike trains: experimental evidence, stochastic modeling, and single neuron variability

The activity of spiking neurons is frequently described by renewal point process models that assume the statistical independence and identical distribution of the intervals between action potentials. However, the assumption of independent intervals must be questioned for many different types of neurons. We review experimental studies that reported the feature of a negative serial correlation of neighboring intervals, commonly observed in neurons in the sensory periphery as well as in central neurons, notably in the mammalian cortex. In our experiments we observed the same short-lived negative serial dependence of intervals in the spontaneous activity of mushroom body extrinsic neurons in the honeybee. To model serial interval correlations of arbitrary lags, we suggest a family of autoregressive point processes. Its marginal interval distribution is described by the generalized gamma model, which includes as special cases the log-normal and gamma distributions, which have been widely used to characterize regular spiking neurons. In numeric simulations we investigated how serial correlation affects the variance of the neural spike count. We show that the experimentally confirmed negative correlation reduces single-neuron variability, as quantified by the Fano factor, by up to 50%, which favors the transmission of a rate code. We argue that the feature of a negative serial correlation is likely to be common to the class of spike-frequency-adapting neurons and that it might have been largely overlooked in extracellular single-unit recordings due to spike sorting errors.

Spike-rate coding and spike-time coding are affected oppositely by different adaptation mechanisms

Spike-frequency adaptation causes reduced spiking during prolonged stimulation, but the full impact of adaptation on neural coding is far more complex, especially if one takes into account the diversity of biophysical mechanisms mediating adaptation and the different ways in which neural information can be encoded. Here, we show that adaptation has opposite effects depending on the neural coding strategy and the biophysical mechanism responsible for adaptation. Under noisy conditions, calcium-activated K(+) current (I(AHP)) improved efficient spike-rate coding at the expense of spike-time coding by regularizing the spike train elicited by slow or constant inputs; noise power was increased at high frequencies but reduced at low frequencies, consistent with noise shaping that improves coding of low- frequency signals. In contrast, voltage-activated M-type K(+) current (I(M)) improved spike-time coding at the expense of spike-rate coding by stopping the neuron from spiking repetitively to slow inputs so that it could generate isolated, well timed spikes in response to fast inputs. Using dynamical systems analysis, we demonstrate how I(AHP) minimizes perturbation of the interspike interval caused by high- frequency noise, whereas I(M) minimizes disruption of spike-timing accuracy caused by repetitive spiking. The dichotomous outcomes are related directly to the distinct activation requirements for I(AHP) and I(M), which in turn dictate whether those currents mediate negative feedback onto spiking or membrane potential. Thus, based on their distinct activation properties, I(AHP) implements noise shaping that improves spike-rate coding of low-frequency signals, whereas I(M) implements high-pass filtering that improves spike-time coding of high- frequency signals.

Noise shaping in neural populations

Many neurons display intrinsic interspike interval correlations in their spike trains. However, the effects of such correlations on information transmission in neural populations are not well understood. We quantified signal processing using linear response theory supported by numerical simulations in networks composed of two different models: One model generates a renewal process where interspike intervals are not correlated while the other generates a nonrenewal process where subsequent interspike intervals are negatively correlated. Our results show that the fractional rate of increase in information rate as a function of network size and stimulus intensity is lower for the nonrenewal model than for the renewal one. We show that this is mostly due to the lower amount of effective noise in the nonrenewal model. We also show the surprising result that coupling has opposite effects in renewal and nonrenewal networks: Excitatory (inhibitory coupling) will decrease (increase) the information rate in renewal networks while inhibitory (excitatory coupling) will decrease (increase) the information rate in nonrenewal networks. We discuss these results and their applicability to other classes of excitable systems.

The response of cortical neurons to in vivo-like input current: theory and experiment : I. Noisy inputs with stationary statistics

The study of several aspects of the collective dynamics of interacting neurons can be highly simplified if one assumes that the statistics of the synaptic input is the same for a large population of similarly behaving neurons (mean field approach). In particular, under such an assumption, it is possible to determine and study all the equilibrium points of the network dynamics when the neuronal response to noisy, in vivo-like, synaptic currents is known. The response function can be computed analytically for simple integrate-and-fire neuron models and it can be measured directly in experiments in vitro. Here we review theoretical and experimental results about the neural response to noisy inputs with stationary statistics. These response functions are important to characterize the collective neural dynamics that are proposed to be the neural substrate of working memory, decision making and other cognitive functions. Applications to the case of time-varying inputs are reviewed in a companion paper (Giugliano et al. in Biol Cybern, 2008). We conclude that modified integrate-and-fire neuron models are good enough to reproduce faithfully many of the relevant dynamical aspects of the neuronal response measured in experiments on real neurons in vitro.

Spontaneous firing statistics and information transfer in electroreceptors of paddlefish

We study information processing in a peripheral sensory receptor system which possesses spontaneous dynamics with two distinct rhythms. Such organization was found in the electrosensory system of paddlefish and is represented by two distinct and unidirectionally coupled oscillators, resulting in biperiodic spontaneous firing patterns of sensory neurons. We use computational modeling to elucidate the functional role of spontaneous oscillations in conveying information from sensory periphery to the brain. We show that biperiodic organization resulting in nonrenewal statistics of background neuronal activity leads to significant improvement in information transfer through the system as compared to an equivalent renewal model.

Dynamic I-V curves are reliable predictors of naturalistic pyramidal-neuron voltage traces

Neuronal response properties are typically probed by intracellular measurements of current-voltage (I-V) relationships during application of current or voltage steps. Here we demonstrate the measurement of a novel I-V curve measured while the neuron exhibits a fluctuating voltage and emits spikes. This dynamic I-V curve requires only a few tens of seconds of experimental time and so lends itself readily to the rapid classification of cell type, quantification of heterogeneities in cell populations, and generation of reduced analytical models. We apply this technique to layer-5 pyramidal cells and show that their dynamic I-V curve comprises linear and exponential components, providing experimental evidence for a recently proposed theoretical model. The approach also allows us to determine the change of neuronal response properties after a spike, millisecond by millisecond, so that postspike refractoriness of pyramidal cells can be quantified. Observations of I-V curves during and in absence of refractoriness are cast into a model that is used to predict both the subthreshold response and spiking activity of the neuron to novel stimuli. The predictions of the resulting model are in excellent agreement with experimental data and close to the intrinsic neuronal reproducibility to repeated stimuli.