State-space analysis of time-varying higher-order spike correlation for multiple neural spike train data

Precise spike coordination between the spiking activities of multiple neurons is suggested as an indication of coordinated network activity in active cell assemblies. Spike correlation analysis aims to identify such cooperative network activity by detecting excess spike synchrony in simultaneously recorded multiple neural spike sequences. Cooperative activity is expected to organize dynamically during behavior and cognition; therefore currently available analysis techniques must be extended to enable the estimation of multiple time-varying spike interactions between neurons simultaneously. In particular, new methods must take advantage of the simultaneous observations of multiple neurons by addressing their higher-order dependencies, which cannot be revealed by pairwise analyses alone. In this paper, we develop a method for estimating time-varying spike interactions by means of a state-space analysis. Discretized parallel spike sequences are modeled as multi-variate binary processes using a log-linear model that provides a well-defined measure of higher-order spike correlation in an information geometry framework. We construct a recursive Bayesian filter/smoother for the extraction of spike interaction parameters. This method can simultaneously estimate the dynamic pairwise spike interactions of multiple single neurons, thereby extending the Ising/spin-glass model analysis of multiple neural spike train data to a nonstationary analysis. Furthermore, the method can estimate dynamic higher-order spike interactions. To validate the inclusion of the higher-order terms in the model, we construct an approximation method to assess the goodness-of-fit to spike data. In addition, we formulate a test method for the presence of higher-order spike correlation even in nonstationary spike data, e.g., data from awake behaving animals. The utility of the proposed methods is tested using simulated spike data with known underlying correlation dynamics. Finally, we apply the methods to neural spike data simultaneously recorded from the motor cortex of an awake monkey and demonstrate that the higher-order spike correlation organizes dynamically in relation to a behavioral demand.     

Fourier analysis of spike train data

A formal representation of nerve spike trains in the form of a sum of rectangular functions is proposed. This formal instantaneous frequency function can be Fourier analyzed. The resulting algorithm has the useful properties of spike by spike calculations and an insensitivity to the mean (carrier) spike rate. The technique is also useful for producing a smooth (filtered) reconstruction of a spike train.     

Analytical reconstruction of the neuronal input current from spike train data

The time course of the current driving action potential generation at a neuron investigated experimentally is in general not measurable directly. In this paper an indirect method is introduced that allows estimation of this unknown current time course using only spike train data. Assuming the leaky integrator model as valid for the action potential encoding site of the investigated neuron, the unknown input current is obtained by determining (analytically) a current time course that upon injection into the leaky integrator model evokes action potential sequences identical to those observed experimentally. Applications of this current-reconstruction procedure to neuronal output data obtained from a leaky integrator model showed that the procedure allows a good estimation of the underlying input current even if the membrane time constant of the investigated neuron is not known exactly. Additionally, an application of current reconstruction to experimental data obtained from a cat muscle spindle primary afferent subject to repeated -stimuli is demonstrated.     

Nonlinear principal components analysis of neuronal spike train data

Many recent approaches to decoding neural spike trains depend critically on the assumption that for low-pass filtered spike trains, the temporal structure is optimally represented by a small number of linear projections onto the data. We therefore tested this assumption of linearity by comparing a linear factor analysis technique (principal components analysis) with a nonlinear neural network based method. It is first shown that the nonlinear technique can reliably identify a neuronally plausible nonlinearity in synthetic spike trains. However, when applied to the outputs from primary visual cortical neurons, this method shows no evidence for significant temporal nonlinearities. The implications of this are discussed. Received: 29 November 1996 / Accepted in revised form: 1 July 1997     

An information-geometric framework for statistical inferences in the neural spike train space

Statistical inferences are essentially important in analyzing neural spike trains in computational neuroscience. Current approaches have followed a general inference paradigm where a parametric probability model is often used to characterize the temporal evolution of the underlying stochastic processes. To directly capture the overall variability and distribution in the space of the spike trains, we focus on a data-driven approach where statistics are defined and computed in the function space in which spike trains are viewed as individual points. To this end, we at first develop a parametrized family of metrics that takes into account different warpings in the time domain and generalizes several currently used spike train distances. These new metrics are essentially penalized L ^p norms, involving appropriate functions of spike trains, with penalties associated with time-warping. The notions of means and variances of spike trains are then defined based on the new metrics when p = 2 (corresponding to the “Euclidean distance”). Using some restrictive conditions, we present an efficient recursive algorithm, termed Matching-Minimization algorithm, to compute the sample mean of a set of spike trains with arbitrary numbers of spikes. The proposed metrics as well as the mean spike trains are demonstrated using simulations as well as an experimental recording from the motor cortex. It is found that all these methods achieve desirable performance and the results support the success of this novel framework.     

Quantification and statistical verification of neuronal stimulus responses from noisy spike train data

Usually neuronal responses to short-lasting stimuli are displayed as peri-stimulus time histogram. The function estimated by such a histogram allows to obtain informations about stimulus-induced postsynaptic events as long as the interpretation is restricted to the first response component after the stimulus. The interpretation of secondary response components is much more difficult, as they may be either due to stimulus effects or represent an echo of the primary response. In the present paper two output functions are developed that do not show such an echoing of responses. The first one, the interspike interval change function, represents an ideal way to quantify a neuronal stimulus response as its amplitude was found to be almost independent of the stimulation strategy used during acquisition of the spike train data. The other function, the displaced impulses function, allows to verify the statistical significance of an observed response component. Both functions may be estimated from stimulus-correlated spike train data, even if the neuron under investigation shows considerable interspike-interval variability in the absence of stimulation. The concepts underlying these neuronal output functions are developed on simulated responses of a Hodgkin-Huxley-type model for a mammalian neuron at body temperature that is exposed to a transient excitatory conductance increase. Additionally, estimation of these output functions is also demonstrated on responses of human soleus motoneurons that were exposed to electrical stimuli of the tibial nerve in the popliteal fossa.     

Signal-to-signal neural networks for improved spike estimation from calcium imaging data

Spiking information of individual neurons is essential for functional and behavioral analysis in neuroscience research. Calcium imaging techniques are generally employed to obtain activities of neuronal populations. However, these techniques result in slowly-varying fluorescence signals with low temporal resolution. Estimating the temporal positions of the neuronal action potentials from these signals is a challenging problem. In the literature, several generative model-based and data-driven algorithms have been studied with varied levels of success. This article proposes a neural network-based signal-to-signal conversion approach, where it takes as input raw-fluorescence signal and learns to estimate the spike information in an end-to-end fashion. Theoretically, the proposed approach formulates the spike estimation as a single channel source separation problem with unknown mixing conditions. The source corresponding to the action potentials at a lower resolution is estimated at the output. Experimental studies on the spikefinder challenge dataset show that the proposed signal-to-signal conversion approach significantly outperforms state-of-the-art-methods in terms of Pearson’s correlation coefficient, Spearman’s rank correlation coefficient and yields comparable performance for the area under the receiver operating characteristics measure. We also show that the resulting system: (a) has low complexity with respect to existing supervised approaches and is reproducible; (b) is layer-wise interpretable, and (c) has the capability to generalize across different calcium indicators.     

Simulations of a ventrolateral medullary neural network for respiratory rhythmogenesis inferred from spike train cross-correlation

Connections among ventrolateral medullary respiratory neurons inferred from spike train analysis were incorporated into a model and simulated with the program SYSTM11 (MacGregor 1987). Inspiratory (I) and expiratory (E) neurons with augmenting (AUG) and decrementing (DEC) discharge patterns and rostral I-E/I neurons exhibited varying degrees of adaptation, but no endogenous bursting properties. Simulation parameters were adjusted so that respiratory phase durations, neuronal discharge patterns, and short-time scale correlations were similar to corresponding measurements from anesthetized, vagotomized, adult cats. Rhythmogenesis persisted when the strength of each set of connections was increased 100% over a smaller effective value. Changes in phase durations and discharge patterns caused by manipulation of connection strengths or population activity led to several predictions. (a) Excitation of the I-E/I population prolongs the inspiratory phase. (b) Rhythmic activity can be reestablished in the absence of I-E/I activity by unpatterned excitation of I-DEC and I-AUG neurons. (c) An increase in I-DEC neuron activity can cause an apneustic respiratory pattern. (d) A decrease in I-DEC neuron activity increases the slope of the inspiratory ramp and shortens inspiration. (e) Excitation of the E-DEC population prolongs the expiratory phase or produces apnea; inhibition of E-DEC neurons reduces expiratory time. (f) Excitation of E-AUG cells causes I-AUG neurons to exhibit a step rather than a ramp increase in firing rate at the onset of their active phase. The results suggest mechanisms by which the duration of each phase of breathing and neuronal discharge patterns may be regulated. Received: 24 February 1993/Accepted in revised form: 8 September 1993     

Estimating the directed information to infer causal relationships in ensemble neural spike train recordings

Advances in recording technologies have given neuroscience researchers access to large amounts of data, in particular, simultaneous, individual recordings of large groups of neurons in different parts of the brain. A variety of quantitative techniques have been utilized to analyze the spiking activities of the neurons to elucidate the functional connectivity of the recorded neurons. In the past, researchers have used correlative measures. More recently, to better capture the dynamic, complex relationships present in the data, neuroscientists have employed causal measures—most of which are variants of Granger causality—with limited success. This paper motivates the directed information, an information and control theoretic concept, as a modality-independent embodiment of Granger’s original notion of causality. Key properties include: (a) it is nonzero if and only if one process causally influences another, and (b) its specific value can be interpreted as the strength of a causal relationship. We next describe how the causally conditioned directed information between two processes given knowledge of others provides a network version of causality: it is nonzero if and only if, in the presence of the present and past of other processes, one process causally influences another. This notion is shown to be able to differentiate between true direct causal influences, common inputs, and cascade effects in more two processes. We next describe a procedure to estimate the directed information on neural spike trains using point process generalized linear models, maximum likelihood estimation and information-theoretic model order selection. We demonstrate that on a simulated network of neurons, it (a) correctly identifies all pairwise causal relationships and (b) correctly identifies network causal relationships. This procedure is then used to analyze ensemble spike train recordings in primary motor cortex of an awake monkey while performing target reaching tasks, uncovering causal relationships whose directionality are consistent with predictions made from the wave propagation of simultaneously recorded local field potentials.     

Sensitivity of different stimulus-timing strategies for the detection of small excitations in noisy spike train data

There are several different strategies to control the timing of a stimulus with respect to the ongoing discharge during the recording of neuronal stimulus-response characteristics. One possible strategy consists of delivering stimuli in such a way that a constant pre-stimulus spike density is reached. Another strategy enforces spike application with a constant stimulus latency after a spontaneous discharge. In this paper the sensitivity of these different strategies for statistical verification of small excitatory response components was investigated. It was found that the difference between observed poststimulus spike distribution and expected spike distribution under the null hypothesis of no stimulus effect was larger using a constant-stimulus-latency (CSL) strategy with an appropriate value for the stimulus latency. Thus, the statistical verification of neuronal response components is clearly facilitated if a CSL strategy is used. This superiority of the CSL strategy is marked, especially for small excitations at neurons discharging slowly with low discharge variability.     

Neural spike train reconstruction from calcium imaging via a signal-shape composition model

Science China Life Sciences -     

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     

The correlation grid: analysis of synchronous spiking in multi-dimensional spike train data and identification of feasible connection architectures

This paper presents a visualization technique specifically designed to support the analysis of synchronous firings in multiple, simultaneously recorded, spike trains. This technique, called the correlation grid, enables investigators to identify groups of spike trains, where each pair of spike trains has a high probability of generating spikes approximately simultaneously or within a constant time shift. Moreover, the correlation grid was developed to help solve the following reverse problem: identification of the connection architecture between spike train generating units, which may produce a spike train dataset similar to the one under analysis. To demonstrate the efficacy of this approach, results are presented from a study of three simulated, noisy, spike train datasets. The parameters of the simulated neurons were chosen to reflect the typical characteristics of cortical pyramidal neurons. The schemes of neuronal connections were not known to the analysts. Nevertheless, the correlation grid enabled the analysts to find the correct connection architecture for each of these three data sets.     

Balanced synaptic input shapes the correlation between neural spike trains

Stimulus properties, attention, and behavioral context influence correlations between the spike times produced by a pair of neurons. However, the biophysical mechanisms that modulate these correlations are poorly understood. With a combined theoretical and experimental approach, we show that the rate of balanced excitatory and inhibitory synaptic input modulates the magnitude and timescale of pairwise spike train correlation. High rate synaptic inputs promote spike time synchrony rather than long timescale spike rate correlations, while low rate synaptic inputs produce opposite results. This correlation shaping is due to a combination of enhanced high frequency input transfer and reduced firing rate gain in the high input rate state compared to the low state. Our study extends neural modulation from single neuron responses to population activity, a necessary step in understanding how the dynamics and processing of neural activity change across distinct brain states.     

Accurate intraocular pressure prediction from applanation response data using genetic algorithm and neural networks

The fact that Goldmann applanation tonometry does not accurately account for individual corneal elastic stiffness often leads to inaccuracy in the measurement of intraocular pressure (IOP). IOP should account not only for the effect of central corneal thickness (CCT) but should also account for other corneal biomechanical factors. A computational method for accurate and reliable determination of IOP is investigated with a modified applanation tonometer in this paper. The proposed method uses a combined genetic algorithm/neural network procedure to match the clinically measured applanation force-displacement history with that obtained from a nonlinear finite element simulation of applanation. An additional advantage of the proposed method is that it also provides the ability to determine CCT and material properties of the cornea from the same applanation response data. The performance of the proposed method has been demonstrated through a parametric study and via comparison with a well known clinical case. The proposed method is also shown to be computationally efficient, which is an important practical consideration for clinical application.