The chronotron: a neuron that learns to fire temporally precise spike patterns

In many cases, neurons process information carried by the precise timings of spikes. Here we show how neurons can learn to generate specific temporally precise output spikes in response to input patterns of spikes having precise timings, thus processing and memorizing information that is entirely temporally coded, both as input and as output. We introduce two new supervised learning rules for spiking neurons with temporal coding of information (chronotrons), one that provides high memory capacity (E-learning), and one that has a higher biological plausibility (I-learning). With I-learning, the neuron learns to fire the target spike trains through synaptic changes that are proportional to the synaptic currents at the timings of real and target output spikes. We study these learning rules in computer simulations where we train integrate-and-fire neurons. Both learning rules allow neurons to fire at the desired timings, with sub-millisecond precision. We show how chronotrons can learn to classify their inputs, by firing identical, temporally precise spike trains for different inputs belonging to the same class. When the input is noisy, the classification also leads to noise reduction. We compute lower bounds for the memory capacity of chronotrons and explore the influence of various parameters on chronotrons' performance. The chronotrons can model neurons that encode information in the time of the first spike relative to the onset of salient stimuli or neurons in oscillatory networks that encode information in the phases of spikes relative to the background oscillation. Our results show that firing one spike per cycle optimizes memory capacity in neurons encoding information in the phase of firing relative to a background rhythm.     

Spatiotemporal burst coding for extracting features of spatiotemporally varying stimuli

Encoding features of spatiotemporally varying stimuli is quite important for understanding the neural mechanisms of various sensory coding. Temporal coding can encode features of time-varying stimulus, and population coding with temporal coding is adequate for encoding spatiotemporal correlation of stimulus features into spatiotemporal activity of neurons. However, little is known about how spatiotemporal features of stimulus are encoded by spatiotemporal property of neural activity. To address this issue, we propose here a population coding with burst spikes, called here spatiotemporal burst (STB) coding. In STB coding, the temporal variation of stimuli is encoded by the precise onset timing of burst spike, and the spatiotemporal correlation of stimuli is emphasized by one specific aspect of burst firing, or spike packet followed by silent interval. To show concretely the role of STB coding, we study the electrosensory system of a weakly electric fish. Weakly electric fish must perceive the information about an object nearby by analyzing spatiotemporal modulations of electric field around it. On the basis of well-characterized circuitry, we constructed a neural network model of the electrosensory system. Here we show that STB coding encodes well the information of object distance and size by extracting the spatiotemporal correlation of the distorted electric field. The burst activity of electrosensory neurons is also affected by feedback signals through synaptic plasticity. We show that the control of burst activity caused by the synaptic plasticity leads to extracting the stimulus features depending on the stimulus context. Our results suggest that sensory systems use burst spikes as a unit of sensory coding in order to extract spatiotemporal features of stimuli from spatially distributed stimuli.     

大鼠初级传入纤维与脊髓背角神经元间的动作电序列的突触传递

外周感觉神经元通过动作电位序列对信号进行编码,这些动作电位序列经过突触传递最终到达脑部。但是各种脉冲序列如何通过神经元之间的化学突触进行传递依然是一个悬而未决的问题。研究了初级传入A6纤维与背角神经元之间各种动作电位序列的突触传递过程。用于刺激的规则,周期、随机脉冲序列由短簇脉冲或单个脉冲构成。定义“事件”(event)为峰峰问期(intefspike interval)小于或等于规定阈值的最长动作电位串,然后从脉冲序列中提取事件间间期(interevent interval,IEI)。用时间,IEI图与回归映射的方法分析IEI序列,结果表明在突触后输出脉冲序列中可以检测到突触前脉冲序列的主要时间结构特征,特别是在短簇脉冲作为刺激单位时。通过计算输入与输出脉冲序列的互信息,发现短簇脉冲可以更可靠地跨突触传递由输入序列携带的神经信息。这些结果表明外周输入脉冲序列的主要时间结构特征可以跨突触传递,在突触传递神经信息的过程中短簇脉冲更为有效。这一研究在从突触传递角度探索神经信息编码方面迈出了一步。     

Recurrent Network Models for Perfect Temporal Integration of Fluctuating Correlated Inputs

Temporal integration of input is essential to the accumulation of information in various cognitive and behavioral processes, and gradually increasing neuronal activity, typically occurring within a range of seconds, is considered to reflect such computation by the brain. Some psychological evidence suggests that temporal integration by the brain is nearly perfect, that is, the integration is non-leaky, and the output of a neural integrator is accurately proportional to the strength of input. Neural mechanisms of perfect temporal integration, however, remain largely unknown. Here, we propose a recurrent network model of cortical neurons that perfectly integrates partially correlated, irregular input spike trains. We demonstrate that the rate of this temporal integration changes proportionately to the probability of spike coincidences in synaptic inputs. We analytically prove that this highly accurate integration of synaptic inputs emerges from integration of the variance of the fluctuating synaptic inputs, when their mean component is kept constant. Highly irregular neuronal firing and spike coincidences are the major features of cortical activity, but they have been separately addressed so far. Our results suggest that the efficient protocol of information integration by cortical networks essentially requires both features and hence is heterotic.     

Neural Coding with Graded Membrane Potential Changes and Spikes

The neural encoding of sensory stimuli is usually investigated for spike responses, although many neurons are known to convey information by graded membrane potential changes. We compare by model simulations how well different dynamical stimuli can be discriminated on the basis of spiking or graded responses. Although a continuously varying membrane potential contains more information than binary spike trains, we find situations where different stimuli can be better discriminated on the basis of spike responses than on the basis of graded responses. Spikes can be superior to graded membrane potential fluctuations if spikes sharpen the temporal structure of neuronal responses by amplifying fast transients of the membrane potential. Such fast membrane potential changes can be induced deterministically by the stimulus or can be due to membrane potential noise that is influenced in its statistical properties by the stimulus. The graded response mode is superior for discrimination between stimuli on a fine time scale.     

Emergence of resonances in neural systems: the interplay between adaptive threshold and short-term synaptic plasticity

In this work we study the detection of weak stimuli by spiking (integrate-and-fire) neurons in the presence of certain level of noisy background neural activity. Our study has focused in the realistic assumption that the synapses in the network present activity-dependent processes, such as short-term synaptic depression and facilitation. Employing mean-field techniques as well as numerical simulations, we found that there are two possible noise levels which optimize signal transmission. This new finding is in contrast with the classical theory of stochastic resonance which is able to predict only one optimal level of noise. We found that the complex interplay between adaptive neuron threshold and activity-dependent synaptic mechanisms is responsible for this new phenomenology. Our main results are confirmed by employing a more realistic FitzHugh-Nagumo neuron model, which displays threshold variability, as well as by considering more realistic stochastic synaptic models and realistic signals such as poissonian spike trains.     

Temporal sub units in dendritic trees

This simulation study examines the possibility that dendritic sub units can be defined according to temporal aspects in the timing of populations of synaptic inputs. A two cell model with passive dendritic trees is used, which is subject to both common and independent synaptic inputs, the presence of common synaptic input results in a tendency for correlated firing in the two cell model. The strength of this correlation is used to measure the efficacy of the common synaptic inputs in modulating the output discharge of each neurone. Our results suggest that a small fraction of the total synaptic input can effectively modulate the timing of output spikes, this phenomenon is not dependent on the physical location of the inputs on the dendritic tree. This phenomenon depends on the presence of temporal correlation between the pre-synaptic spike trains that provide the common input. We propose to refer to these as temporal sub units.     

Studying spike trains using a van Rossum metric with a synapse-like filter

Spike trains are unreliable. For example, in the primary sensory areas, spike patterns and precise spike times will vary between responses to the same stimulus. Nonetheless, information about sensory inputs is communicated in the form of spike trains. A challenge in understanding spike trains is to assess the significance of individual spikes in encoding information. One approach is to define a spike train metric, allowing a distance to be calculated between pairs of spike trains. In a good metric, this distance will depend on the information the spike trains encode. This method has been used previously to calculate the timescale over which the precision of spike times is significant. Here, a new metric is constructed based on a simple model of synaptic conductances which includes binding site depletion. Including binding site depletion in the metric means that a given individual spike has a smaller effect on the distance if it occurs soon after other spikes. The metric proves effective at classifying neuronal responses by stimuli in the sample data set of electro-physiological recordings from the primary auditory area of the zebra finch fore-brain. This shows that this is an effective metric for these spike trains suggesting that in these spike trains the significance of a spike is modulated by its proximity to previous spikes. This modulation is a putative information-coding property of spike trains.     

STOCHASTIC RESONANCE MEASURED BY THEMUTUAL INFORMATION IN THE NEURONSTHAT TRANSMIT CHAOTIC SPIKE TRAINS

The spike trains generated by a neuron model are studied by the methods of nonlinear time series analysis. The results show that the spike trains are chaotic. To investigate effect of noise on transmission of chaotic spike trains, this chaotic spike trains are used as a discrete subthreshold input signal to the integrate-and-fire neuronal model and the FitzHugh-Nagumo(FHN) neuronal model working in noisy environment. The mutual information between the input spike trains and the output spike trains is calculated, the result shows that the transformation of information encoded by the chaotic spike trains is optimized by some level of noise, and stochastic resonance(SR) measured by mutual information is a property available for neurons to transmit chaotic spike trains.     

STOCHASTIC RESONANCE MEASURED BY THEMUTUAL INFORMATION IN THE NEURONSTHAT TRANSMIT CHAOTIC SPIKE TRAINS

１ＩｎｔｒｏｄｕｃｔｉｏｎＩｔｉｓｗｅｌｌｋｎｏｗｎｔｈａｔｎｅｒｖｅｃｅｌｌｓｗｏｒｋｉｎｎｏｉｓｙｅｎｖｉｒｏｎｍｅｎｔ,ａｎｄｎｏｉｓｅｓｏｕｒｃｅｓｒａｎｇｉｎｇｆｒｏｍｉｎｔｅｒｎａｌｔｈｅｒｍａｌｎｏｉｓｅｔｏｅｘｔｅｒｎａｌｐｅｒｔｕｒｂａｔｉｏｎ．Ｏｎｅｐｕｚｚｌｉｎｇｐｒｏｂｌｅｍｉｓｈｏｗｄｏｎｅｒｖｅｃｅｌｌｓａｃｃｏｍｍｏｄａｔｅｎｏｉｓｅｉｎｃｏｄｉｎｇａｎｄｔｒａｎｓｆｏｒｍｉｎｇｉｎｆｏｒｍａｔｉｏｎ,ｒｅｃｅｎｔｒｅｓｅａｒｃｈｓｈｏｗｓｔｈａｔｎｏｉｓｅｍａｙｐ…     

Dynamic Alignment Models for Neural Coding

Recently, there have been remarkable advances in modeling the relationships between the sensory environment, neuronal responses, and behavior. However, most models cannot encompass variable stimulus-response relationships such as varying response latencies and state or context dependence of the neural code. Here, we consider response modeling as a dynamic alignment problem and model stimulus and response jointly by a mixed pair hidden Markov model (MPH). In MPHs, multiple stimulus-response relationships (e.g., receptive fields) are represented by different states or groups of states in a Markov chain. Each stimulus-response relationship features temporal flexibility, allowing modeling of variable response latencies, including noisy ones. We derive algorithms for learning of MPH parameters and for inference of spike response probabilities. We show that some linear-nonlinear Poisson cascade (LNP) models are a special case of MPHs. We demonstrate the efficiency and usefulness of MPHs in simulations of both jittered and switching spike responses to white noise and natural stimuli. Furthermore, we apply MPHs to extracellular single and multi-unit data recorded in cortical brain areas of singing birds to showcase a novel method for estimating response lag distributions. MPHs allow simultaneous estimation of receptive fields, latency statistics, and hidden state dynamics and so can help to uncover complex stimulus response relationships that are subject to variable timing and involve diverse neural codes.     

Spike timing precision of neuronal circuits

Spike timing is believed to be a key factor in sensory information encoding and computations performed by the neurons and neuronal circuits. However, the considerable noise and variability, arising from the inherently stochastic mechanisms that exist in the neurons and the synapses, degrade spike timing precision. Computational modeling can help decipher the mechanisms utilized by the neuronal circuits in order to regulate timing precision. In this paper, we utilize semi-analytical techniques, which were adapted from previously developed methods for electronic circuits, for the stochastic characterization of neuronal circuits. These techniques, which are orders of magnitude faster than traditional Monte Carlo type simulations, can be used to directly compute the spike timing jitter variance, power spectral densities, correlation functions, and other stochastic characterizations of neuronal circuit operation. We consider three distinct neuronal circuit motifs: Feedback inhibition, synaptic integration, and synaptic coupling. First, we show that both the spike timing precision and the energy efficiency of a spiking neuron are improved with feedback inhibition. We unveil the underlying mechanism through which this is achieved. Then, we demonstrate that a neuron can improve on the timing precision of its synaptic inputs, coming from multiple sources, via synaptic integration: The phase of the output spikes of the integrator neuron has the same variance as that of the sample average of the phases of its inputs. Finally, we reveal that weak synaptic coupling among neurons, in a fully connected network, enables them to behave like a single neuron with a larger membrane area, resulting in an improvement in the timing precision through cooperation.     

A computational model of cellular mechanisms of temporal coding in the medial geniculate body (MGB)

Acoustic stimuli are often represented in the early auditory pathway as patterns of neural activity synchronized to time-varying features. This phase-locking predominates until the level of the medial geniculate body (MGB), where previous studies have identified two main, largely segregated response types: Stimulus-synchronized responses faithfully preserve the temporal coding from its afferent inputs, and Non-synchronized responses, which are not phase locked to the inputs, represent changes in temporal modulation by a rate code. The cellular mechanisms underlying this transformation from phase-locked to rate code are not well understood. We use a computational model of a MGB thalamocortical neuron to test the hypothesis that these response classes arise from inferior colliculus (IC) excitatory afferents with divergent properties similar to those observed in brain slice studies. Large-conductance inputs exhibiting synaptic depression preserved input synchrony as short as 12.5 ms interclick intervals, while maintaining low firing rates and low-pass filtering responses. By contrast, small-conductance inputs with Mixed plasticity (depression of AMPA-receptor component and facilitation of NMDA-receptor component) desynchronized afferent inputs, generated a click-rate dependent increase in firing rate, and high-pass filtered the inputs. Synaptic inputs with facilitation often permitted band-pass synchrony along with band-pass rate tuning. These responses could be tuned by changes in membrane potential, strength of the NMDA component, and characteristics of synaptic plasticity. These results demonstrate how the same synchronized input spike trains from the inferior colliculus can be transformed into different representations of temporal modulation by divergent synaptic properties.     

Modeling the impact of common noise inputs on the network activity of retinal ganglion cells

Synchronized spontaneous firing among retinal ganglion cells (RGCs), on timescales faster than visual responses, has been reported in many studies. Two candidate mechanisms of synchronized firing include direct coupling and shared noisy inputs. In neighboring parasol cells of primate retina, which exhibit rapid synchronized firing that has been studied extensively, recent experimental work indicates that direct electrical or synaptic coupling is weak, but shared synaptic input in the absence of modulated stimuli is strong. However, previous modeling efforts have not accounted for this aspect of firing in the parasol cell population. Here we develop a new model that incorporates the effects of common noise, and apply it to analyze the light responses and synchronized firing of a large, densely-sampled network of over 250 simultaneously recorded parasol cells. We use a generalized linear model in which the spike rate in each cell is determined by the linear combination of the spatio-temporally filtered visual input, the temporally filtered prior spikes of that cell, and unobserved sources representing common noise. The model accurately captures the statistical structure of the spike trains and the encoding of the visual stimulus, without the direct coupling assumption present in previous modeling work. Finally, we examined the problem of decoding the visual stimulus from the spike train given the estimated parameters. The common-noise model produces Bayesian decoding performance as accurate as that of a model with direct coupling, but with significantly more robustness to spike timing perturbations.     

Differential coding of humoral stimuli by timing and amplitude of intracellular calcium spike trains

The ubiquitous Ca2(+)-phosphoinositide pathway transduces extracellular signals to cellular effectors. Using a mathematical model, we simulated intracellular Ca2+ fluctuations in hepatocytes upon humoral stimulation. We estimated the information encoded about random humoral stimuli in these Ca2+ spike trains using an information-theoretic approach based on stimulus estimation methods. We demonstrate accurate transfer of information about random humoral signals with low temporal cutoff frequencies. In contrast, our results suggest that high-frequency stimuli are poorly transduced by the transmembrane machinery. We found that humoral signals are encoded in both the timing and amplitude of intracellular Ca2+ spikes. The information transmitted per spike is similar to that of sensory neuronal systems, in spite of several orders of magnitude difference in firing rate.     

Integrate-and-fire vs Poisson models of LGN input to V1 cortex: noisier inputs reduce orientation selectivity

One of the reasons the visual cortex has attracted the interest of computational neuroscience is that it has well-defined inputs. The lateral geniculate nucleus (LGN) of the thalamus is the source of visual signals to the primary visual cortex (V1). Most large-scale cortical network models approximate the spike trains of LGN neurons as simple Poisson point processes. However, many studies have shown that neurons in the early visual pathway are capable of spiking with high temporal precision and their discharges are not Poisson-like. To gain an understanding of how response variability in the LGN influences the behavior of V1, we study response properties of model V1 neurons that receive purely feedforward inputs from LGN cells modeled either as noisy leaky integrate-and-fire (NLIF) neurons or as inhomogeneous Poisson processes. We first demonstrate that the NLIF model is capable of reproducing many experimentally observed statistical properties of LGN neurons. Then we show that a V1 model in which the LGN input to a V1 neuron is modeled as a group of NLIF neurons produces higher orientation selectivity than the one with Poisson LGN input. The second result implies that statistical characteristics of LGN spike trains are important for V1’s function. We conclude that physiologically motivated models of V1 need to include more realistic LGN spike trains that are less noisy than inhomogeneous Poisson processes.     

Reconstruction of Underlying Nonlinear Deterministic Dynamics Embedded in Noisy Spike Trains

An experimentally recorded time series formed by the exact times of occurrence of the neuronal spikes (spike train) is likely to be affected by observational noise that provokes events mistakenly confused with neuronal discharges, as well as missed detection of genuine neuronal discharges. The points of the spike train may also suffer a slight jitter in time due to stochastic processes in synaptic transmission and to delays in the detecting devices. This study presents a procedure aimed at filtering the embedded noise (denoising the spike trains) the spike trains based on the hypothesis that recurrent temporal patterns of spikes are likely to represent the robust expression of a dynamic process associated with the information carried by the spike train. The rationale of this approach is tested on simulated spike trains generated by several nonlinear deterministic dynamical systems with embedded observational noise. The application of the pattern grouping algorithm (PGA) to the noisy time series allows us to extract a set of points that form the reconstructed time series. Three new indices are defined for assessment of the performance of the denoising procedure. The results show that this procedure may indeed retrieve the most relevant temporal features of the original dynamics. Moreover, we observe that additional spurious events affect the performance to a larger extent than the missing of original points. Thus, a strict criterion for the detection of spikes under experimental conditions, thus reducing the number of spurious spikes, may raise the possibility to apply PGA to detect endogenous deterministic dynamics in the spike train otherwise masked by the observational noise.     

Predicting spike occurrence and neuronal responsiveness from LFPs in primary somatosensory cortex

Local Field Potentials (LFPs) integrate multiple neuronal events like synaptic inputs and intracellular potentials. LFP spatiotemporal features are particularly relevant in view of their applications both in research (e.g. for understanding brain rhythms, inter-areal neural communication and neuronal coding) and in the clinics (e.g. for improving invasive Brain-Machine Interface devices). However the relation between LFPs and spikes is complex and not fully understood. As spikes represent the fundamental currency of neuronal communication this gap in knowledge strongly limits our comprehension of neuronal phenomena underlying LFPs. We investigated the LFP-spike relation during tactile stimulation in primary somatosensory (S-I) cortex in the rat. First we quantified how reliably LFPs and spikes code for a stimulus occurrence. Then we used the information obtained from our analyses to design a predictive model for spike occurrence based on LFP inputs. The model was endowed with a flexible meta-structure whose exact form, both in parameters and structure, was estimated by using a multi-objective optimization strategy. Our method provided a set of nonlinear simple equations that maximized the match between models and true neurons in terms of spike timings and Peri Stimulus Time Histograms. We found that both LFPs and spikes can code for stimulus occurrence with millisecond precision, showing, however, high variability. Spike patterns were predicted significantly above chance for 75% of the neurons analysed. Crucially, the level of prediction accuracy depended on the reliability in coding for the stimulus occurrence. The best predictions were obtained when both spikes and LFPs were highly responsive to the stimuli. Spike reliability is known to depend on neuron intrinsic properties (i.e. on channel noise) and on spontaneous local network fluctuations. Our results suggest that the latter, measured through the LFP response variability, play a dominant role.     

Network variability limits stimulus-evoked spike timing precision in retinal ganglion cells

Visual, auditory, somatosensory, and olfactory stimuli generate temporally precise patterns of action potentials (spikes). It is unclear, however, how the precision of spike generation relates to the pattern and variability of synaptic input elicited by physiological stimuli. We determined how synaptic conductances evoked by light stimuli that activate the rod bipolar pathway control spike generation in three identified types of mouse retinal ganglion cells (RGCs). The relative amplitude, timing, and impact of excitatory and inhibitory input differed dramatically between On and Off RGCs. Spikes evoked by repeated somatic injection of identical light-evoked synaptic conductances were more temporally precise than those evoked by light. However, the precision of spikes evoked by conductances that varied from trial to trial was similar to that of light-evoked spikes. Thus, the rod bipolar pathway modulates different RGCs via unique combinations of synaptic input, and RGC temporal variability reflects variability in the input this circuit provides.     

Representation of dynamical stimuli in populations of threshold neurons

Many sensory or cognitive events are associated with dynamic current modulations in cortical neurons. This raises an urgent demand for tractable model approaches addressing the merits and limits of potential encoding strategies. Yet, current theoretical approaches addressing the response to mean- and variance-encoded stimuli rarely provide complete response functions for both modes of encoding in the presence of correlated noise. Here, we investigate the neuronal population response to dynamical modifications of the mean or variance of the synaptic bombardment using an alternative threshold model framework. In the variance and mean channel, we provide explicit expressions for the linear and non-linear frequency response functions in the presence of correlated noise and use them to derive population rate response to step-like stimuli. For mean-encoded signals, we find that the complete response function depends only on the temporal width of the input correlation function, but not on other functional specifics. Furthermore, we show that both mean- and variance-encoded signals can relay high-frequency inputs, and in both schemes step-like changes can be detected instantaneously. Finally, we obtain the pairwise spike correlation function and the spike triggered average from the linear mean-evoked response function. These results provide a maximally tractable limiting case that complements and extends previous results obtained in the integrate and fire framework.