A cross-interval spike train analysis: the correlation between spike generation and temporal integration of doublets

A stochastic spike train analysis technique is introduced to reveal the correlation between the firing of the next spike and the temporal integration period of two consecutive spikes (i.e., a doublet). Statistics of spike firing times between neurons are established to obtain the conditional probability of spike firing in relation to the integration period. The existence of a temporal integration period is deduced from the time interval between two consecutive spikes fired in a reference neuron as a precondition to the generation of the next spike in a compared neuron. This analysis can show whether the coupled spike firing in the compared neuron is correlated with the last or the second-to-last spike in the reference neuron. Analysis of simulated and experimentally recorded biological spike trains shows that the effects of excitatory and inhibitory temporal integration are extracted by this method without relying on any subthreshold potential recordings. The analysis also shows that, with temporal integration, a neuron driven by random firing patterns can produce fairly regular firing patterns under appropriate conditions. This regularity in firing can be enhanced by temporal integration of spikes in a chain of polysynaptically connected neurons. The bandpass filtering of spike firings by temporal integration is discussed. The results also reveal that signal transmission delays may be attributed not just to conduction and synaptic delays, but also to the delay time needed for temporal integration. Received: 3 March 1997 / Accepted in revised form: 6 November 1997     

A joint interspike interval difference stochastic spike train analysis: detecting local trends in the temporal firing patterns of single neurons

We introduce a stochastic spike train analysis method called joint interspike interval difference (JISID) analysis. By design, this method detects changes in firing interspike intervals (ISIs), called local trends, within a 4-spike pattern in a spike train. This analysis classifies 4-spike patterns that have similar incremental changes. It characterizes the higher-order serial dependence in spike firing relative to changes in the firing history. Mathematically, this spike train analysis describes the statistical joint distribution of consecutive changes in ISIs, from which the serial dependence of the changes in higher-order intervals can be determined. It is similar to the joint interspike interval (JISI) analysis, except that the joint distribution of consecutive ISI differences (ISIDs) is quantified. The graphical location of points in the JISID scatter plot reveals the local trends in firing (i.e., monotonically increasing, monotonically decreasing, or transitional firing). The trajectory of these points in the serial-JISID plot traces the time evolution of these trends represented by a 5-spike pattern, while points in the JISID scatter plot represent trends of a 4-spike pattern. We provide complete theoretical interpretations of the JISID analysis. We also demonstrate that this method indeed identifies firing trends in both simulated spike trains and spike trains recorded from cultured neurons. Received: 13 May 1997 / Accepted in revised form: 9 December 1998     

Visual pattern discrimination by population retinal ganglion cells’ activities during natural movie stimulation

In the visual system, neurons often fire in synchrony, and it is believed that synchronous activities of group neurons are more efficient than single cell response in transmitting neural signals to down-stream neurons. However, whether dynamic natural stimuli are encoded by dynamic spatiotemporal firing patterns of synchronous group neurons still needs to be investigated. In this paper we recorded the activities of population ganglion cells in bullfrog retina in response to time-varying natural images (natural scene movie) using multi-electrode arrays. In response to some different brief section pairs of the movie, synchronous groups of retinal ganglion cells (RGCs) fired with similar but different spike events. We attempted to discriminate the movie sections based on temporal firing patterns of single cells and spatiotemporal firing patterns of the synchronous groups of RGCs characterized by a measurement of subsequence distribution discrepancy. The discrimination performance was assessed by a classification method based on Support Vector Machines. Our results show that different movie sections of the natural movie elicited reliable dynamic spatiotemporal activity patterns of the synchronous RGCs, which are more efficient in discriminating different movie sections than the temporal patterns of the single cells’ spike events. These results suggest that, during natural vision, the down-stream neurons may decode the visual information from the dynamic spatiotemporal patterns of the synchronous group of RGCs’ activities.     

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.     

Why spikes? Hebbian learning and retrieval of time-resolved excitation patterns

Hebbian learning allows a network of spiking neurons to store and retrieve spatio-temporal patterns with a time resolution of 1 ms, despite the long postsynaptic and dendritic integration times. To show this, we introduce and analyze a model of spiking neurons, the spike response model, with a realistic distribution of axonal delays and with realistic postsynaptic potentials. Learning is performed by a local Hebbian rule which is based on the synchronism of presynaptic neurotransmitter release and some short-acting postsynaptic process. The time window of this synchronism determines the temporal resolution of pattern retrieval, which can be initiated by applying a short external stimulus pattern. Furthermore, a rate quantization is found in dependence upon the threshold value of the neurons, i.e., in a given time a pattern runsn times as often as learned, wheren is a positive integer (n 0). We show that all information about the spike pattern is lost if only mean firing rates (temporal average) or ensemble activities (spatial average) are considered. An average over several retrieval runs in order to generate a post-stimulus time histogram may also deteriorate the signal. The full information on a pattern is contained in the spike raster of a single run. Our results stress the importance, and advantage, of coding by spatio-temporal spike patterns instead of firing rates and average ensemble activity. The implications regarding modelling and experimental data analysis are discussed.     

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.     

Firing pattern and synchronization property analysis in a network model of the olfactory bulb

In the olfactory system, both the temporal spike structure and spatial distribution of neuronal activity are important for processing odor information. In this paper, a biophysically-detailed, spiking neuronal model is used to simulate the activity of olfactory bulb. It is shown that by varying some key parameters such as maximal conductances of Ks and Nap the spike train of single neuron can exhibit various firing patterns. Synchronization in coupled neurons is also investigated as the coupling strength varying in the situation of two neurons and network. It is illustrated that the coupled neurons can exhibit different types of pattern when the coupling strength varies. These results may be instructive to understand information transmission in olfactory system.     

Temporal Features of Spike Trains in the Moth Antennal Lobe Revealed by a Comparative Time-Frequency Analysis

The discrimination of complex sensory stimuli in a noisy environment is an immense computational task. Sensory systems often encode stimulus features in a spatiotemporal fashion through the complex firing patterns of individual neurons. To identify these temporal features, we have developed an analysis that allows the comparison of statistically significant features of spike trains localized over multiple scales of time-frequency resolution. Our approach provides an original way to utilize the discrete wavelet transform to process instantaneous rate functions derived from spike trains, and select relevant wavelet coefficients through statistical analysis. Our method uncovered localized features within olfactory projection neuron (PN) responses in the moth antennal lobe coding for the presence of an odor mixture and the concentration of single component odorants, but not for compound identities. We found that odor mixtures evoked earlier responses in biphasic response type PNs compared to single components, which led to differences in the instantaneous firing rate functions with their signal power spread across multiple frequency bands (ranging from 0 to 45.71 Hz) during a time window immediately preceding behavioral response latencies observed in insects. Odor concentrations were coded in excited response type PNs both in low frequency band differences (2.86 to 5.71 Hz) during the stimulus and in the odor trace after stimulus offset in low (0 to 2.86 Hz) and high (22.86 to 45.71 Hz) frequency bands. These high frequency differences in both types of PNs could have particular relevance for recruiting cellular activity in higher brain centers such as mushroom body Kenyon cells. In contrast, neurons in the specialized pheromone-responsive area of the moth antennal lobe exhibited few stimulus-dependent differences in temporal response features. These results provide interesting insights on early insect olfactory processing and introduce a novel comparative approach for spike train analysis applicable to a variety of neuronal data sets.     

Reconstruction of firing rate changes across neuronal populations by temporally deconvolved Ca2+ imaging

Methods to record action potential (AP) firing in many individual neurons are essential to unravel the function of complex neuronal circuits in the brain. A promising approach is bolus loading of Ca(2+) indicators combined with multiphoton microscopy. Currently, however, this technique lacks cell-type specificity, has low temporal resolution and cannot resolve complex temporal firing patterns. Here we present simple solutions to these problems. We identified neuron types by colocalizing Ca(2+) signals of a red-fluorescing indicator with genetically encoded markers. We reconstructed firing rate changes from Ca(2+) signals by temporal deconvolution. This technique is efficient, dramatically enhances temporal resolution, facilitates data interpretation and permits analysis of odor-response patterns across thousands of neurons in the zebrafish olfactory bulb. Hence, temporally deconvolved Ca(2+) imaging (TDCa imaging) resolves limitations of current optical recording techniques and is likely to be widely applicable because of its simplicity, robustness and generic principle.     

Analyzing multiple spike trains with nonparametric granger causality

Simultaneous recordings of spike trains from multiple single neurons are becoming commonplace. Understanding the interaction patterns among these spike trains remains a key research area. A question of interest is the evaluation of information flow between neurons through the analysis of whether one spike train exerts causal influence on another. For continuous-valued time series data, Granger causality has proven an effective method for this purpose. However, the basis for Granger causality estimation is autoregressive data modeling, which is not directly applicable to spike trains. Various filtering options distort the properties of spike trains as point processes. Here we propose a new nonparametric approach to estimate Granger causality directly from the Fourier transforms of spike train data. We validate the method on synthetic spike trains generated by model networks of neurons with known connectivity patterns and then apply it to neurons simultaneously recorded from the thalamus and the primary somatosensory cortex of a squirrel monkey undergoing tactile stimulation.     

Gap junction permeability modulated by dopamine exerts effects on spatial and temporal correlation of retinal ganglion cells’ firing activities

Synchronized activities among retinal ganglion cells (RGCs) via gap junctions can be increased by exogenous dopamine (DA). During DA application, single neurons’ firing activities become more synchronized with its adjacent neighbors. One intriguing question is how the enhanced spatial synchronization alters the temporal firing structure of single neurons. In the present study, firing activities of bullfrog’s dimming detectors in response to binary pseudo-random checker-board flickering were recorded via a multi-channel recording system. DA was applied in the retina to modulate synchronized activities between RGCs, and the effect of DA on firing activities of single neurons was examined. It was found that, during application of DA, synchronized activities between single neuron and its neighboring neurons was enhanced. At the meantime, the temporal structures of single neuron spike train changed significantly, and the temporal correlation in single neuron’s response decreased. The pharmacological study results indicated that the activation of D1 receptor might have effects on gap junction permeability between RGCs. Our results suggested that the dopaminergic pathway participated in the modulation of spatial and temporal correlation of RGCs’ firing activities, and may exert critical effects on visual information processing in the retina.     

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.     

Processing of Auditory Midbrain Interspike Intervals by Model Neurons

A question central to sensory processing is how signal information is encoded and processed by single neurons. Stimulus features can be represented through rate coding (via firing rate), temporal coding (via firing synchronization to temporal periodicities), or temporal encoding (via intricate patterns of spike trains). Of the three, examples of temporal encoding are the least documented. One region in which temporal encoding is currently being explored is the auditory midbrain. Midbrain neurons in the plainfin midshipman generate different interspike interval (ISI) distributions depending on the frequencies of the concurrent vocal signals. However, these distributions differ only along certain lengths of ISIs, so that any neurons trying to distinguish the distributions would have to respond selectively to specific ISI ranges. We used this empirical observation as a realistic challenge with which to explore the plausibility of ISI-tuned neurons that could validate this form of temporal encoding. The resulting modeled cells—point neurons optimized through multidimensional searching—were successfully tuned to discriminate patterns in specific ranges of ISIs. Achieving this task, particularly with simplified neurons, strengthens the credibility of ISI coding in the brain and lends credence to its role in auditory processing.     

STDP allows fast rate-modulated coding with Poisson-like spike trains

Spike timing-dependent plasticity (STDP) has been shown to enable single neurons to detect repeatedly presented spatiotemporal spike patterns. This holds even when such patterns are embedded in equally dense random spiking activity, that is, in the absence of external reference times such as a stimulus onset. Here we demonstrate, both analytically and numerically, that STDP can also learn repeating rate-modulated patterns, which have received more experimental evidence, for example, through post-stimulus time histograms (PSTHs). Each input spike train is generated from a rate function using a stochastic sampling mechanism, chosen to be an inhomogeneous Poisson process here. Learning is feasible provided significant covarying rate modulations occur within the typical timescale of STDP (~10-20 ms) for sufficiently many inputs (~100 among 1000 in our simulations), a condition that is met by many experimental PSTHs. Repeated pattern presentations induce spike-time correlations that are captured by STDP. Despite imprecise input spike times and even variable spike counts, a single trained neuron robustly detects the pattern just a few milliseconds after its presentation. Therefore, temporal imprecision and Poisson-like firing variability are not an obstacle to fast temporal coding. STDP provides an appealing mechanism to learn such rate patterns, which, beyond sensory processing, may also be involved in many cognitive tasks.     

Predicting the Responses of Repetitively Firing Neurons to Current Noise

We used phase resetting methods to predict firing patterns of rat subthalamic nucleus (STN) neurons when their rhythmic firing was densely perturbed by noise. We applied sequences of contiguous brief (0.5–2 ms) current pulses with amplitudes drawn from a Gaussian distribution (10–100 pA standard deviation) to autonomously firing STN neurons in slices. Current noise sequences increased the variability of spike times with little or no effect on the average firing rate. We measured the infinitesimal phase resetting curve (PRC) for each neuron using a noise-based method. A phase model consisting of only a firing rate and PRC was very accurate at predicting spike timing, accounting for more than 80% of spike time variance and reliably reproducing the spike-to-spike pattern of irregular firing. An approximation for the evolution of phase was used to predict the effect of firing rate and noise parameters on spike timing variability. It quantitatively predicted changes in variability of interspike intervals with variation in noise amplitude, pulse duration and firing rate over the normal range of STN spontaneous rates. When constant current was used to drive the cells to higher rates, the PRC was altered in size and shape and accurate predictions of the effects of noise relied on incorporating these changes into the prediction. Application of rate-neutral changes in conductance showed that changes in PRC shape arise from conductance changes known to accompany rate increases in STN neurons, rather than the rate increases themselves. Our results show that firing patterns of densely perturbed oscillators cannot readily be distinguished from those of neurons randomly excited to fire from the rest state. The spike timing of repetitively firing neurons may be quantitatively predicted from the input and their PRCs, even when they are so densely perturbed that they no longer fire rhythmically.     

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.     

The Possible Role of Spike Patterns in Cortical Information Processing

When the same visual stimulus is presented across many trials, neurons in the visual cortex receive stimulus-related synaptic inputs that are reproducible across trials (S) and inputs that are not (N). The variability of spike trains recorded in the visual cortex and their apparent lack of spike-to-spike correlations beyond that implied by firing rate fluctuations, has been taken as evidence for a low S/N ratio. A recent re-analysis of in vivo cortical data revealed evidence for spike-to-spike correlations in the form of spike patterns. We examine neural dynamics at a higher S/N in order to determine what possible role spike patterns could play in cortical information processing. In vivo-like spike patterns were obtained in model simulations. Superpositions of multiple sinusoidal driving currents were especially effective in producing stable long-lasting patterns. By applying current pulses that were either short and strong or long and weak, neurons could be made to switch from one pattern to another. Cortical neurons with similar stimulus preferences are located near each other, have similar biophysical properties and receive a large number of common synaptic inputs. Hence, recordings of a single neuron across multiple trials are usually interpreted as the response of an ensemble of these neurons during one trial. In the presence of distinct spike patterns across trials there is ambiguity in what would be the corresponding ensemble, it could consist of the same spike pattern for each neuron or a set of patterns across neurons. We found that the spiking response of a neuron receiving these ensemble inputs was determined by the spike-pattern composition, which, in turn, could be modulated dynamically as a means for cortical information processing.     

Synchronization and sensitivity enhancement of the Hodgkin-Huxley neurons due to inhibitory inputs

Recent experimental results imply that inhibitory postsynaptic potentials can play a functional role in realizing synchronization of neuronal firing in the brain. In order to examine the relation between inhibition and synchronous firing of neurons theoretically, we analyze possible effects of synchronization and sensitivity enhancement caused by inhibitory inputs to neurons with a biologically realistic model of the Hodgkin-Huxley equations. The result shows that, after an inhibitory spike, the firing probability of a single postsynaptic neuron exposed to random excitatory background activity oscillates with time. The oscillation of the firing probability can be related to synchronous firing of neurons receiving an inhibitory spike simultaneously. Further, we show that when an inhibitory spike input precedes an excitatory spike input, the presence of such preceding inhibition raises the firing probability peak of the neuron after the excitatory input. The result indicates that an inhibitory spike input can enhance the sensitivity of the postsynaptic neuron to the following excitatory spike input. Two neural network models based on these effects on postsynaptic neurons caused by inhibitory inputs are proposed to demonstrate possible mechanisms of detecting particular spatiotemporal spike patterns. Received: 15 April 1999 /Accepted in revised form: 25 November 1999     

Compromised axonal functionality after neurodegeneration,concussion and/or traumatic brain injury

Axonal swellings are almost universal in neurodegenerative diseases of the central nervous system, including Alzheimer’s and Parkinson’s disease. Concussions and traumatic brain injuries can also produce cognitive and behavioral deficits by compromising neuronal morphology. Using a spike metric analysis, we characterize computationally the effects of such axonal varicosities on spike train propagation by comparing Poisson spike train classes before and after propagation through a prototypical axonal enlargement, or focused axonal swelling. Misclassification of spike train classes and low-pass filtering of firing rate activity increases with more pronounced axonal injury. We show that confusion matrices and a calculation of the loss of transmitted information provide a very practical way to characterize how injured neurons compromise the signal processing and faithful conductance of spike trains. The method demonstrates that (i) neural codes encoded with low firing rates are more robust to injury than those encoded with high firing rates, (ii) classification depends upon the length of the spike train used to encode information, and (iii) axonal injuries reduce the variance of spike trains within a given stimulus class. The work introduces a novel theoretical and computational framework to quantify the interplay between electrophysiological dynamics with focused axonal swellings generated by injury or other neurodegenerative processes. It further suggests how pharmacology and plasticity may play a role in recovery of neural computation. Ultimately, the work bridges vast experimental observations of in vitro morphological pathologies with post-traumatic cognitive and behavioral dysfunction.     

Rhythm-induced spike-timing patterns characterized by 1D firing maps

We explore patterns in the spike timing of neurons receiving periodic inputs, with an emphasis on stable characteristics which are realized in both models and in-vitro whole-cell recordings. We report on whole-cell recordings of pyramidal CA1 cells from rat hippocampus and entorhinal cortex and compare this data to model simulations. Cells were injected with a constant current to induce a steady firing rate and then a modest rhythm was added which altered the spike times and their corresponding phases relative to the rhythm. For both experiment and theory the relationship between consecutive spike phases is characterized by a probability distribution with peaks concentrated near a one-dimensional firing map. As is well-known, stable fixed points of this map correspond to the neuron phase-locking to the rhythm. We show that the interaction between noise and sufficiently steep maps can also cause a new kind of spike-time organization, in which consecutive spike time pairs organize into discrete clusters, with transitions between these clusters proceeding in a fixed sequence. This structure is not just a vestige of the noise-free dynamics. This slow dynamics and temporal organization in the relationship between consecutive spike phases is not evident in either the neuron’s voltage traces or single phase or interspike interval histograms. Furthermore, the consecutive spike relationship is also evident in consecutive ISIs, and hence this ordering can be observed without detailed knowledge of the rhythm (e.g. without concurrent LFP recordings).