The role of individual spikes and spike patterns in population coding of stimulus location in rat somatosensory cortex

This report addresses the nature of population coding in sensory cortex by applying information theoretic analysis to data recorded simultaneously from neuron pairs located in primary somatosensory cortex of anaesthetised rats. We studied how cortical spike trains code for the location of a whisker stimulus on the rat's snout. We found that substantially more information was conveyed by 10 ms precision spike timing compared with that conveyed by the number of spikes counted over a 40 ms response interval. Most of this information was accounted for by the timing of individual spikes. In particular, it was the first post-stimulus spikes that were crucial. Spike patterns within individual cells played a smaller role; spike patterns across cells were negligible. This pattern of results was robust both to the exact nature of the stimulus set and to the precision at which spikes were binned.     

Analysis of Slow (Theta) Oscillations as a Potential Temporal Reference Frame for Information Coding in Sensory Cortices

While sensory neurons carry behaviorally relevant information in responses that often extend over hundreds of milliseconds, the key units of neural information likely consist of much shorter and temporally precise spike patterns. The mechanisms and temporal reference frames by which sensory networks partition responses into these shorter units of information remain unknown. One hypothesis holds that slow oscillations provide a network-intrinsic reference to temporally partitioned spike trains without exploiting the millisecond-precise alignment of spikes to sensory stimuli. We tested this hypothesis on neural responses recorded in visual and auditory cortices of macaque monkeys in response to natural stimuli. Comparing different schemes for response partitioning revealed that theta band oscillations provide a temporal reference that permits extracting significantly more information than can be obtained from spike counts, and sometimes almost as much information as obtained by partitioning spike trains using precisely stimulus-locked time bins. We further tested the robustness of these partitioning schemes to temporal uncertainty in the decoding process and to noise in the sensory input. This revealed that partitioning using an oscillatory reference provides greater robustness than partitioning using precisely stimulus-locked time bins. Overall, these results provide a computational proof of concept for the hypothesis that slow rhythmic network activity may serve as internal reference frame for information coding in sensory cortices and they foster the notion that slow oscillations serve as key elements for the computations underlying perception.     

Spike-based population coding and working memory

Compelling behavioral evidence suggests that humans can make optimal decisions despite the uncertainty inherent in perceptual or motor tasks. A key question in neuroscience is how populations of spiking neurons can implement such probabilistic computations. In this article, we develop a comprehensive framework for optimal, spike-based sensory integration and working memory in a dynamic environment. We propose that probability distributions are inferred spike-per-spike in recurrently connected networks of integrate-and-fire neurons. As a result, these networks can combine sensory cues optimally, track the state of a time-varying stimulus and memorize accumulated evidence over periods much longer than the time constant of single neurons. Importantly, we propose that population responses and persistent working memory states represent entire probability distributions and not only single stimulus values. These memories are reflected by sustained, asynchronous patterns of activity which make relevant information available to downstream neurons within their short time window of integration. Model neurons act as predictive encoders, only firing spikes which account for new information that has not yet been signaled. Thus, spike times signal deterministically a prediction error, contrary to rate codes in which spike times are considered to be random samples of an underlying firing rate. As a consequence of this coding scheme, a multitude of spike patterns can reliably encode the same information. This results in weakly correlated, Poisson-like spike trains that are sensitive to initial conditions but robust to even high levels of external neural noise. This spike train variability reproduces the one observed in cortical sensory spike trains, but cannot be equated to noise. On the contrary, it is a consequence of optimal spike-based inference. In contrast, we show that rate-based models perform poorly when implemented with stochastically spiking neurons.     

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

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

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.     

Extra spike formation in sensory neurons and the disruption of afferent spike patterning

The peculiar pseudounipolar geometry of primary sensory neurons can lead to ectopic generation of "extra spikes" in the region of the dorsal root ganglion potentially disrupting the fidelity of afferent signaling. We have used an explicit model of myelinated vertebrate sensory neurons to investigate the location and mechanism of extra spike formation, and its consequences for distortion of afferent impulse patterning. Extra spikes originate in the initial segment axon under conditions in which the soma spike becomes delayed and broadened. The broadened soma spike then re-excites membrane it has just passed over, initiating an extra spike which propagates outwards into the main conducting axon. Extra spike formation depends on cell geometry, electrical excitability, and the recent history of impulse activity. Extra spikes add to the impulse barrage traveling toward the spinal cord, but they also travel antidromically in the peripheral nerve colliding with and occluding normal orthodromic spikes. As a result there is no net increase in afferent spike number. However, extra spikes render firing more staccato by increasing the number of short and long interspike intervals in the train at the expense of intermediate intervals. There may also be more complex changes in the pattern of afferent spike trains, and hence in afferent signaling.     

Neuronal precision and the limits for acoustic signal recognition in a small neuronal network

Recognition of acoustic signals may be impeded by two factors: extrinsic noise, which degrades sounds before they arrive at the receiver’s ears, and intrinsic neuronal noise, which reveals itself in the trial-to-trial variability of the responses to identical sounds. Here we analyzed how these two noise sources affect the recognition of acoustic signals from potential mates in grasshoppers. By progressively corrupting the envelope of a female song, we determined the critical degradation level at which males failed to recognize a courtship call in behavioral experiments. Using the same stimuli, we recorded intracellularly from auditory neurons at three different processing levels, and quantified the corresponding changes in spike train patterns by a spike train metric, which assigns a distance between spike trains. Unexpectedly, for most neurons, intrinsic variability accounted for the main part of the metric distance between spike trains, even at the strongest degradation levels. At consecutive levels of processing, intrinsic variability increased, while the sensitivity to external noise decreased. We followed two approaches to determine critical degradation levels from spike train dissimilarities, and compared the results with the limits of signal recognition measured in behaving animals.     

Processing of sensory signals by a non-spiking neuron in the leech

The non-spiking neurons 151 are present as bilateral pairs in each midbody ganglion of the leech nervous system and they are electrically coupled to several motorneurons. Intracellular recordings were used to investigate how these neurons process input from the mechanosensory P neurons in isolated ganglia. Induction of spike trains (15 Hz) in single P cells evoked responses that combined depolarizing and hyperpolarizing phases in cells 151. The phasic depolarizations, transmitted through spiking interneurons, reversed at around -20 mV. The hyperpolarization had two components, both reversing at around -65 mV, and which were inhibited by strychnine (10 micromol l(-1)). The faster component was transmitted through spiking interneurons and the slower component through a direct P-151 interaction. Short trains (<400 ms) of P cell spikes (15 Hz) evoked the phasic depolarizations superimposed on the hyperpolarization, while long spike trains (>500 ms) produced a succession of depolarizations that masked the hyperpolarizing phase. The amplitude and duration of the hyperpolarization reached their maximum at the initial spikes in a train, while the depolarizations persisted throughout the duration of the stimulus train. Both phases of the response were relatively unaffected by the spike frequency (5-25 Hz). The non-spiking neurons 151 processed the sensory signals in the temporal rather than in the amplitude domain.     

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.     

Regulation of spike timing in visual cortical circuits

A train of action potentials (a spike train) can carry information in both the average firing rate and the pattern of spikes in the train. But can such a spike-pattern code be supported by cortical circuits? Neurons in vitro produce a spike pattern in response to the injection of a fluctuating current. However, cortical neurons in vivo are modulated by local oscillatory neuronal activity and by top-down inputs. In a cortical circuit, precise spike patterns thus reflect the interaction between internally generated activity and sensory information encoded by input spike trains. We review the evidence for precise and reliable spike timing in the cortex and discuss its computational role.     

Determination of the precision of spike timing in the visual cortex of anaesthetised cats

Neuronal cortical spike trains contain precisely replicating patterns whose presence cannot be accounted for by chance production. A comparison of the number of triplets of spikes present two times with the number of doublets replicated three times in the same window duration gives a frequency-insensitive measure of this type of fine temporal organisation. By varying the tolerance with which such precisely replicating patterns are detected, one can evaluate the accuracy of spike timing in spike trains. In the sample of data here analysed, it was found that replicating patterns were best seen in the precision range 0.4–1.4 ms (a result evidently at variance with a simple ‘integrate and fire’ model of neurons). Surprisingly, the fine temporal structure of spike trains thus evidenced was stronger at relatively low firing rate discharges and was present in both the ‘spontaneous’ and ‘evoked’ responses.     

Computing Complex Visual Features with Retinal Spike Times

Neurons in sensory systems can represent information not only by their firing rate, but also by the precise timing of individual spikes. For example, certain retinal ganglion cells, first identified in the salamander, encode the spatial structure of a new image by their first-spike latencies. Here we explore how this temporal code can be used by downstream neural circuits for computing complex features of the image that are not available from the signals of individual ganglion cells. To this end, we feed the experimentally observed spike trains from a population of retinal ganglion cells to an integrate-and-fire model of post-synaptic integration. The synaptic weights of this integration are tuned according to the recently introduced tempotron learning rule. We find that this model neuron can perform complex visual detection tasks in a single synaptic stage that would require multiple stages for neurons operating instead on neural spike counts. Furthermore, the model computes rapidly, using only a single spike per afferent, and can signal its decision in turn by just a single spike. Extending these analyses to large ensembles of simulated retinal signals, we show that the model can detect the orientation of a visual pattern independent of its phase, an operation thought to be one of the primitives in early visual processing. We analyze how these computations work and compare the performance of this model to other schemes for reading out spike-timing information. These results demonstrate that the retina formats spatial information into temporal spike sequences in a way that favors computation in the time domain. Moreover, complex image analysis can be achieved already by a simple integrate-and-fire model neuron, emphasizing the power and plausibility of rapid neural computing with spike times.     

Predictive Coding of Dynamical Variables in Balanced Spiking Networks

Two observations about the cortex have puzzled neuroscientists for a long time. First, neural responses are highly variable. Second, the level of excitation and inhibition received by each neuron is tightly balanced at all times. Here, we demonstrate that both properties are necessary consequences of neural networks that represent information efficiently in their spikes. We illustrate this insight with spiking networks that represent dynamical variables. Our approach is based on two assumptions: We assume that information about dynamical variables can be read out linearly from neural spike trains, and we assume that neurons only fire a spike if that improves the representation of the dynamical variables. Based on these assumptions, we derive a network of leaky integrate-and-fire neurons that is able to implement arbitrary linear dynamical systems. We show that the membrane voltage of the neurons is equivalent to a prediction error about a common population-level signal. Among other things, our approach allows us to construct an integrator network of spiking neurons that is robust against many perturbations. Most importantly, neural variability in our networks cannot be equated to noise. Despite exhibiting the same single unit properties as widely used population code models (e.g. tuning curves, Poisson distributed spike trains), balanced networks are orders of magnitudes more reliable. Our approach suggests that spikes do matter when considering how the brain computes, and that the reliability of cortical representations could have been strongly underestimated.     

Modeling the insect feeding decision system: A new approach utilizing fuzzy system theory

Conclusions Preliminary tests of this model have been successful using as inputs responses to pure chemicals. When the trained model is presented with computer-synthesized spike trains as inputs, it predicts behavioral results that are plausible. Future tests will employ more complex stimuli, such as binary and trinary mixtures of compounds and natural stimuli from host and non-host plants.The use of this operational computer model represents a powerful new approach to the currently intractable problem of deciphering complex, acrossfiber sensory codes. This approach provides a more comprehensive and objective assessment of the data than do the current alternatives, such as analysis of individual records by trained experts. In addition, the model can assimilate data from any number of inputs, which is well beyond the abilities of individual investigators compiling data manually. Finally, the model is dynamic and can be updated with new data, and is sufficiently general to be used for any species or comparison across species.A model trained with responses to plants ranging from highly acceptable host plants to highly deterrent non-host plants would be capable of simulating the insect's entire repertoire of feeding decisions. Because the model can divulge the rules by which it made a decision, inspection of these rules will permit identification of which sensory inputs and input patterns were the most important for that decision, thereby providing a crucial insight into complex across-fiber sensory coding. In addition, probing the model with synthesized spike trains as inputs will permit understanding of input/output relationships over wider ranges or in more combinations than would be possible to test experimentally, and thus would be helpful in predicting behavioral responses to multi-component mixtures. Modeling acceptance or rejection, and thereby determining the operative rules associated with these behaviors, may supply a rational basis for the development of baits and antifeedants and provide direction for genetic engineering of semiochemical in plants.     

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.     

Inhibition of the HIV-1 Spike by Single-PG9/16-Antibody Binding Suggests a Coordinated-Activation Model for Its Three Protomeric Units

The HIV-1 spike is composed of three protomeric units, each containing a peripheral gp120 and a transmembrane gp41 subunit. Binding to the CD4 and the chemokine receptors triggers them to mediate virus entry into cells by membrane fusion. The spikes also represent the major target for neutralizing antibodies (Abs) against the virus. We have studied how two related broadly neutralizing Abs, PG9 and PG16, react with the spike. Unexpectedly, this also suggested how the functions of the individual protomers in the spike depend on each other. The Abs have been shown to bind the V1/V2 loops of gp120, located at the top of the spike. Using blue native-polyacrylamide gel electrophoresis (BN-PAGE), we show that only single Abs or antigen-binding fragments could bind to the spikes of HIV-1 virus-like particles. Apparently, binding to one gp120 sterically interferes with binding to the other two subunits in the spike top. Despite this constraint, all of the protomers of the spike became resistant to CD4 binding and subsequent formation of the coreceptor binding site. These activities were measured by monitoring the sequential complex formation of the spike first with Abs and then with soluble 2d- or 4d-CD4 or with soluble CD4 and the CD4 inducible coreceptor binding site Ab 17b in BN-PAGE. The inhibition of the spike by single-Ab binding suggested that the activation reactions of the individual protomeric units are linked to each other in a coordinated activation process.     

Critical parameters of the spike trains in a cell assembly: coding of turn direction by the giant interneurons of the cockroach

Cockroaches (Periplaneta americana) respond to air displacement produced by an approaching predator by turning and running away. A set of 4 bilateral pairs of ventral giant interneurons is important in determining turn direction. Wind from a given side is known to produce more spikes, an earlier onset of the spike trains, and different fine temporal patterning, in the ipsilateral vs the contralateral set of these interneurons. Here we investigate which of these spike train parameters the cockroach actually uses to determine the direction it will turn.We delivered controlled wind puffs from the right front, together with intracellular injection of spike trains in a left ventral giant interneuron, under conditions where the animal could make normally directed turning movements of the legs and body. In trials where our stimuli caused the left side to give both the first spike and more total spikes than the right, but where our injected spike train included none of the normal fine temporal patterning, 92% of the evoked turns were to the rightopposite of normal (Figs. 4–6). In trials where the left side gave the first spike, but the right side gave more spikes, 100% of the turns were to the left-the normal direction (Figs. 8, 9). Comparable results were obtained when each of the left giant interneurons 1, 2 or 3 were electrically stimulated, and when either weak or stronger wind puffs were used. Stimulating a left giant interneuron electrically in the absence of a wind puff evoked an escape-like turn on 9% of the trials, and these were all to the right (Fig. 9).These results indicate that fine temporal patterning in the spike trains is not necessary, and information about which side gives the first spike is not sufficient, to determine turn direction. Rather, the key parameter appears to be relative numbers of action potentials in the left vs the right group of cells. These conclusions were supported by similar experiments in which extracellular stimulation of several left giant interneurons was paired with right wind (Figs. 11, 12).Abbreviations GI giant interneuron - vGI ventral giant interneuron - dGI dorsal giant interneuron - LY Lucifer yellow - CF carboxyfluorescein     

Population coding of stimulus location in rat somatosensory cortex.

This study explores the nature of population coding in sensory cortex by applying information theoretic analyses to neuron pairs recorded simultaneously from rat barrel cortex. We quantified the roles of individual spikes and spike patterns in encoding whisker stimulus location. 82%-85% of the total information was contained in the timing of individual spikes: first spike time was particularly crucial. Spike patterns within neurons accounted for the remaining 15%-18%. Neuron pairs located in the same barrel column coded redundantly, whereas pairs in neighboring barrel columns coded independently. The barrel cortical population code for stimulus location appears to be the time of single neurons' first poststimulus spikes-a fast, robust coding mechanism that does not rely on "synergy" in crossneuronal spike patterns.     

Repetitive firing: a quantitative study of feedback in model encoders

Recognition of nonlinearities in the neuronal encoding of repetitive spike trains has generated a number of models to explain this behavior. Here we develop the mathematics and a set of tests for two such models: the leaky integrator and the variable-gamma model. Both of these are nearly sufficient to explain the dynamic behavior of a number of repetitively firing, sensory neurons. Model parameters can be related to possible underlying basic mechanisms. Summed and nonsummed, spike- locked negative feedback are examined in conjunction with the models. Transfer functions are formulated to predict responses to steady state, steps, and sinusoidally varying stimuli in which output data are the times of spike-train events only. An electrical analog model for the leaky integrator is tested to verify predicted responses. Curve fitting and parameter variation techniques are explored for the purpose of extracting basic model parameters from spike train data. Sinusoidal analysis of spike trains appear to be a very accurate method for determining spike-locked feedback parameters, and it is to a large extent a model independent method that may also be applied to neuronal responses.