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     

Quantification of Clustering in Joint Interspike Interval Scattergrams of Spike Trains

Joint interval scattergrams are usually employed in determining serial correlations between events of spike trains. However, any inherent structures in such scattergrams that are often seen in experimental records are not quantifiable by serial correlation coefficients. Here, we develop a method to quantify clustered structures in any two-dimensional scattergram of pairs of interspike intervals. The method gives a cluster coefficient as well as clustering density function that could be used to quantify clustering in scattergrams obtained from first- or higher-order interval return maps of single spike trains, or interspike interval pairs drawn from simultaneously recorded spike trains. The method is illustrated using numerical spike trains as well as in vitro pairwise recordings of rat striatal tonically active neurons.     

Temporal patterns recognized by a network of coordinated time delays and coincidence detectors

A computational model of a neuronal network is described which performs a fundamental task of general perception: recognition of temporal patterns in continuous and uncued neuronal spike trains. The presented network is able to recognize each pattern element (100 ms interval composed of sets of 10, 20, 30 and 40 ms interspike intervals combined in linear order) as it arrives. Its operation is based upon biologically plausible filtering mechanisms and population neurodynamics.     

Multiscale spike train variability in primary electrosensory afferents.

Spike train variability is of fundamental importance for understanding how information is encoded and processed in the nervous system. Most studies in this area have focused on short-term variability, as characterized by the coefficient of variation of the interspike interval distribution. Here we discuss the importance of extending the analysis of spike train variability to longer time scales that span multiple interspike intervals. Recent experimental and modeling studies of probability coding (P-type) electrosensory afferent nerve fibers in weakly electric fish have provided new insights into the functional importance of multiscale spike train variability. P-type afferent spike trains are moderately irregular on short time scales of a few milliseconds, but show significantly enhanced regularity on time scales of a few hundred milliseconds. This increased regularity is beyond what would be expected for a renewal process model in which successive intervals are uncorrelated. Modeling studies suggest that the correlation structure that underlies spike train regularization arises from relative refractory effects associated with a dynamic spike threshold. Spike train regularization in P-type afferents has been shown to significantly enhance signal detectability and information transmission on time scales that are functionally relevant for electrolocation behavior.     

Interspike interval dependency from arterial chemoreceptors

The carotid body impulse generator has been previously characterized as a Poisson-type random process. We examined the validity of this characterization by analyzing sinus nerve spike trains for interspike interval dependency. Fifteen single chemoreceptive afferents were recorded in vivo under hypoxic-hypercapnic conditions, and approximately 1,000 consecutive interspike intervals for each fiber were timed and analyzed for serial dependence. The same set of intervals placed in shuffled order served as a control series without serial dependence. The original spike interval trains showed significantly negative first-order serial correlation coefficients and less variability in joint interval distributions than did the shuffled interval trains. These results suggest that the chemoreceptor afferent train is not random and may reflect a negative feedback system operating within the carotid body that limits variation about a mean frequency.     

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

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

Influence of temporal correlation of synaptic input on the rate and variability of firing in neurons

The spike trains that transmit information between neurons are stochastic. We used the theory of random point processes and simulation methods to investigate the influence of temporal correlation of synaptic input current on firing statistics. The theory accounts for two sources for temporal correlation: synchrony between spikes in presynaptic input trains and the unitary synaptic current time course. Simulations show that slow temporal correlation of synaptic input leads to high variability in firing. In a leaky integrate-and-fire neuron model with spike afterhyperpolarization the theory accurately predicts the firing rate when the spike threshold is higher than two standard deviations of the membrane potential fluctuations. For lower thresholds the spike afterhyperpolarization reduces the firing rate below the theory's predicted level when the synaptic correlation decays rapidly. If the synaptic correlation decays slower than the spike afterhyperpolarization, spike bursts can occur during single broad peaks of input fluctuations, increasing the firing rate over the prediction. Spike bursts lead to a coefficient of variation for the interspike intervals that can exceed one, suggesting an explanation of high coefficient of variation for interspike intervals observed in vivo.     

Transmission of spike trains by nerve fibers after near-threshold stimulation

Spike trains in group A nerve fibers were studied in anesthetized cats in response to stimulation of the hind-limb nerves by random (Poisson) and regular sequences of stimuli. In response to above-threshold stimulation spike trains in the nerve fibers were shown to differ from the stimulating trains purely in the absence of intervals less than 1–1.5 msec in duration, as a result of the presence of a refractory period. With near-threshold stimulation with an average frequency of over 10 per second, spike trains differed significantly from the stimulating trains, as reflected in histograms of interspike intervals, the shape of the intensity function, and the magnitude of the coefficient of correlation for successive intervals. It is postulated that changes in the structure of the spike trains conveyed by a nerve fiber are attributable to the presence of after-activity.A. A. Bogomolets Institute of Physiology, Academy of Sciences of the Ukrainian SSR, Kiev. Translated from Neirofiziologiya, Vol. 8, No. 1, pp. 91–98, January–February, 1976.     

A Pseudo-Markov Model for Series of Neuronal Spike Events

Spike trains of spontaneous neuronal activity in the rabbit brain are submitted to statistical analyses based on the following pseudo-Markov model. The nerve cell is supposed to alternate between a bursting and a resting state. The numbers of consecutive spikes within each state are assumed to be independent integer-valued random variables with discrete probability distributions. Given the state, the interspike intervals are independent real-valued random variables. The two state semi-Markov model is obtained as a special case when the discrete distributions are geometrical. Statistical second-order properties of recorded spike trains are compared with those predicted by the model on the basis of known first-order properties. For that purpose, serial correlation coefficients and intensity functions for spike trains produced by the model are computed. A comparison between observed and predicted results for the spontaneous activity of 17 brain cells yields a good fit in eight cells and discloses some salient features of the statistical structure in the activity of six other cells. By making it feasible to compute theoretical correlograms, the model may advance the understanding of empirical correlograms. The possibilities for integrating this statistical model of spike trains with a model of the mechanism of spike train production are discussed.     

Statistical inference on spontaneous neuronal discharge patterns

A statistical analysis was performed on extracellularly recorded spike trains of spontaneously active mesencephalic reticular neurons of rats. Only stationary records were used for detailed examination. The moments of interspike intervals were computed, hypothesis of renewal process and its specific forms was tested. Implications for statistical methodology are considered on the basis of the results. The main emphasis is laid on the connection between experimental results and stochastic neuronal models.     

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.     

Recurrence plots of neuronal spike trains

The recently developed qualitative method of diagnosis of dynamical systems — recurrence plots has been applied to the analysis of dynamics of neuronal spike trains recorded from cerebellum and red nucleus of anesthetized cats. Recurrence plots revealed robust and common changes in the similarity structure of interspike interval sequences as well as significant deviations from randomness in serial ordering of intervals. Recurring episodes of alike, quasi-deterministic firing patterns suggest the spontaneous modulation of the dynamical complexity of the trajectories of observed neurons. These modulations are associated with changing dynamical properties of a neuronal spike-train-generating system. Their existence is compatible with the information processing paradigm of attractor neural networks.     

Multifractal analysis of wide dynamic range neuron discharge profiles in normal rats and in rats with sciatic nerve constriction

Using normal rats and rats with a chronic constriction injury of the sciatic nerve (injured, model of Bennett-Xie) we investigated the possibility of classifying, by statistical tools, the temporal sequences of neuronal discharges during different noxious and non-noxious stimuli. An analysis was made of both the distribution of the inter-spike intervals and the temporal density of spike trains, the latter being studied within the framework of stochastic universal multifractals, to allow the identification of different random processes involved in the discharge distributions through the Levy index alpha. The statistical analysis shows that the parametrization based on the Levy index seems able to discriminate between different noxious stimuli (mechanical pinching and thermal), both in normal and injured animals. Furthermore, comparing normal and injured animals, although the spontaneous basal and non-noxious stimuli (brushing) evoked activities presented different frequencies, these seem to have the same multifractal structure, while the corresponding statistics of the inter-spike intervals are quite different. This information might be relevant to the understanding of a code of neuronal firing and to the modelling of temporal patterns in acute and chronic noxious signals.     

神经放电序列模式的一种识别方法及应用

提出了神经放电序列模式识别的一种新方法。首先,把放电序列用阶梯状的响应函数来表示,然后定义了其一阶、二阶形式导数以及形式积分。这三个特征量均有着不同的几何和物理意义,因此采用这三个特征量来刻画神经放电序列的模式,就可以较全面地表示其特征。对神经放电序列的重构也表明通过这几个特征量可以很好地反映序列中所包含的信息。作为应用例子,这种量化方法用来研究冷热感受器模型所产生的放电模式,结果表明它能够识别在不同温度条件下的放电模式。     

Discharge patterns of tonically firing human motoneurones

We have attempted to reconcile the different patterns of distribution of interspike intervals that are found in motoneurones made to discharge by intracellular injection of constant current in reduced animal preparations and by voluntary control in human subjects. We recorded long spike trains from single motor units in three human muscles made to discharge at constant mean frequencies with the help of auditory and visual feedback. The distribution of interspike intervals in each spike train was analysed quantitatively. We found that the different pattern of discharge of the human motor units could be accounted for when due allowance was made for the variability of the drive to the human motoneurone which arose because of the feedback process used to maintain the target frequency. A model testing this hypothesis gave results that were qualitatively consistent with the human data.     

The stochastic properties of input spike trains control neuronal arithmetic

In the nervous system, the representation of signals is based predominantly on the rate and timing of neuronal discharges. In most everyday tasks, the brain has to carry out a variety of mathematical operations on the discharge patterns. Recent findings show that even single neurons are capable of performing basic arithmetic on the sequences of spikes. However, the interaction of the two spike trains, and thus the resulting arithmetic operation may be influenced by the stochastic properties of the interacting spike trains. If we represent the individual discharges as events of a random point process, then an arithmetical operation is given by the interaction of two point processes. Employing a probabilistic model based on detection of coincidence of random events and complementary computer simulations, we show that the point process statistics control the arithmetical operation being performed and, particularly, that it is possible to switch from subtraction to division solely by changing the distribution of the inter-event intervals of the processes. Consequences of the model for evaluation of binaural information in the auditory brainstem are demonstrated. The results accentuate the importance of the stochastic properties of neuronal discharge patterns for information processing in the brain; further studies related to neuronal arithmetic should therefore consider the statistics of the interacting spike trains.     

A Superposition Model of the Spontaneous Activity of Cerebellar Purkinje Cells

Based on physiological evidence for multiple firing zones in the dendritic arborizations of cerebellar Purkinje cells, a superposition model is proposed for spike triggering in these cells. Spike trains from 10 Purkinje cells were analyzed in terms of independence of interspike intervals and the properties of their variance-time curves. The results of this analysis were found consistent with the hypothesis that the spike train of a cerebellar Purkinje cell is the pooled output of a relatively large number of independent component processes. Simplifying assumptions as to the statistical nature of these processes lead to a very rough estimate of the number of firing zones.     

A method of neuronal spike train study based on autoregressive analysis

The method of autoregressive (AR) analysis for neuronal spike trains (NST) is proposed in the paper. The AR model and the Green's function as well as the power spectral density function are used to process and analyse the neuronal interspike interval (ISI) sequences of cat's first somatosensory area of cortex (SI area) under various situations. With these methods the characteristics of the ISI sequence such as the AR order and parameters, memory property, correlativity and periodicity etc. can be extracted.     

Background neuronal activity in thenucleus supraopticus and its modifications under the influence of parathyroid hormone

In experiments on anesthetized cats, we found that i. v. injection of 5.0 U/kg of parathyroid hormone (PTH) results in modifications of the statistical parameters of the neuronal impulse activity in thenucleus supraopticus (SO) of the hypothalamus. Sliding frequency graphs, histograms of interspike intervals, autocorrelograms, and serial correlation coefficients were plotted and calculated before and after PTH injections; their comparison demonstrates that the hormone significantly modulates the temporal organization of spike trains generated by the neurons of this nucleus. We observed that PTH mostly activated SO neurons and diminished the level of spike grouping in their activity. The effect of PTH to a certain level depended on the initial frequency of background activity: an increase in the spiking frequency was typical of primarily dominating “low-frequency” neurons, while “high-frequency” units were mostly inhibited. The possible mechanisms of the observed modifications are discussed.     

Neuronal spike trains and stochastic point processes. I. The single spike train

In a growing class of neurophysiological experiments, the train of impulses (“spikes”) produced by a nerve cell is subjected to statistical treatment involving the time intervals between spikes. The statistical techniques available for the analysis of single spike trains are described and related to the underlying mathematical theory, that of stochastic point processes, i.e., of stochastic processes whose realizations may be described as series of point events occurring in time, separated by random intervals. For single stationary spike trains, several orders of complexity of statistical treatment are described; the major distinction is that between statistical measures that depend in an essential way on the serial order of interspike intervals and those that are order-independent. The interrelations among the several types of calculations are shown, and an attempt is made to ameliorate the current nomenclatural confusion in this field. Applications, interpretations, and potential difficulties of the statistical techniques are discussed, with special reference to types of spike trains encountered experimentally. Next, the related types of analysis are described for experiments which involve repeated presentations of a brief, isolated stimulus. Finally, the effects of nonstationarity, e.g. long-term changes in firing rate, on the various statistical measures are discussed. Several commonly observed patterns of spike activity are shown to be differentially sensitive to such changes. A companion paper covers the analysis of simultaneously observed spike trains.