Statistical technique for analysing functional connectivity of multiple spike trains

A new statistical technique, the Cox method, used for analysing functional connectivity of simultaneously recorded multiple spike trains is presented. This method is based on the theory of modulated renewal processes and it estimates a vector of influence strengths from multiple spike trains (called reference trains) to the selected (target) spike train. Selecting another target spike train and repeating the calculation of the influence strengths from the reference spike trains enables researchers to find all functional connections among multiple spike trains. In order to study functional connectivity an "influence function" is identified. This function recognises the specificity of neuronal interactions and reflects the dynamics of postsynaptic potential. In comparison to existing techniques, the Cox method has the following advantages: it does not use bins (binless method); it is applicable to cases where the sample size is small; it is sufficiently sensitive such that it estimates weak influences; it supports the simultaneous analysis of multiple influences; it is able to identify a correct connectivity scheme in difficult cases of "common source" or "indirect" connectivity. The Cox method has been thoroughly tested using multiple sets of data generated by the neural network model of the leaky integrate and fire neurons with a prescribed architecture of connections. The results suggest that this method is highly successful for analysing functional connectivity of simultaneously recorded multiple spike trains.

[1]  V. I. Kryukov,et al.  A new statistical method for identifying interconnections between neuronal network elements , 1985, Biological Cybernetics.

[2]  Petr Lánský,et al.  The parameters of the stochastic leaky integrate-and-fire neuronal model , 2006, Journal of Computational Neuroscience.

[3]  G. P. Moore,et al.  Neuronal spike trains and stochastic point processes. II. Simultaneous spike trains. , 1967, Biophysical journal.

[4]  Kathryn B. Laskey,et al.  Neural Coding: Higher-Order Temporal Patterns in the Neurostatistics of Cell Assemblies , 2000, Neural Computation.

[5]  H. L. Le Roy,et al.  Proceedings of the Fifth Berkeley Symposium on Mathematical Statistics and Probability; Vol. IV , 1969 .

[6]  D. Brillinger.  Maximum likelihood analysis of spike trains of interacting nerve cells , 2004, Biological Cybernetics.

[7]  L. Paninski.  Maximum likelihood estimation of cascade point-process neural encoding models , 2004, Network.

[8]  Fabio Rigat,et al.  Bayesian modelling and analysis of spatio-temporal neuronal networks , 2006 .

[9]  J. Kiefer,et al.  An Introduction to Stochastic Processes. , 1956 .

[10]  H C Tuckwell,et al.  Neuronal interspike time distributions and the estimation of neurophysiological and neuroanatomical parameters. , 1978, Journal of theoretical biology.

[11]  R. Kass,et al.  Multiple neural spike train data analysis: state-of-the-art and future challenges , 2004, Nature Neuroscience.

[12]  M. Wilson,et al.  Analyzing Functional Connectivity Using a Network Likelihood Model of Ensemble Neural Spiking Activity , 2005, Neural Computation.

[13]  Sonja Grün,et al.  Significance of joint-spike events based on trial-shuffling by efficient combinatorial methods , 2003, Complex..

[14]  Sonja Grün,et al.  Spatially organized spike correlation in cat visual cortex , 2007, Neurocomputing.

[15]  S. K. Srinivasan,et al.  Stochastic Point Processes , 1978 .

[16]  L Stuart,et al.  The correlation grid: analysis of synchronous spiking in multi-dimensional spike train data and identification of feasible connection architectures. , 2005, Bio Systems.

[17]  Christian Borgelt,et al.  Generation and Selection of Surrogate Methods for Correlation Analysis , 2010 .

[18]  David R. Brillinger,et al.  Nerve Cell Spike Train Data Analysis: A Progression of Technique , 1992 .

[19]  Lee E Miller,et al.  Inferring functional connections between neurons , 2008, Current Opinion in Neurobiology.

[20]  Wulfram Gerstner,et al.  SPIKING NEURON MODELS Single Neurons , Populations , Plasticity , 2002 .

[21]  Rong Jin,et al.  Identifying Functional Connectivity in Large-Scale Neural Ensemble Recordings: A Multiscale Data Mining Approach , 2009, Neural Computation.

[22]  H. L. Bryant,et al.  Identification of synaptic interactions , 1976, Biological Cybernetics.

[23]  R. Quian Quiroga,et al.  Unsupervised Spike Detection and Sorting with Wavelets and Superparamagnetic Clustering , 2004, Neural Computation.

[24]  E. S. Chornoboy,et al.  Maximum likelihood identification of neural point process systems , 1988, Biological Cybernetics.

[25]  Rajesh P. N. Rao,et al.  Bayesian brain : probabilistic approaches to neural coding , 2006 .

[26]  David R. Brillinger,et al.  Measuring The Association Of Point Processes: A Case History , 1976 .

[27]  Christian Borgelt,et al.  Efficient Identification of Assembly Neurons within Massively Parallel Spike Trains , 2009, Comput. Intell. Neurosci..

[28]  Daryl J. Daley,et al.  An Introduction to the Theory of Point Processes , 2013 .

[29]  Liam Paninski,et al.  Statistical models for neural encoding, decoding, and optimal stimulus design. , 2007, Progress in brain research.

[30]  Uri T Eden,et al.  A point process framework for relating neural spiking activity to spiking history, neural ensemble, and extrinsic covariate effects. , 2005, Journal of neurophysiology.

[31]  P. A. W. Lewis,et al.  Multivariate point processes , 2018, Point Processes.

[32]  Matthias Bethge,et al.  Bayesian Inference for Spiking Neuron Models with a Sparsity Prior , 2007, NIPS.

[33]  D. Brillinger.  The Identification of Point Process Systems , 1975 .

[34]  Jonathan W. Pillow,et al.  Likelihood-based approaches to modeling the neural code , 2007 .

[35]  S. Grün,et al.  Higher-Order Correlations and Cumulants , 2010 .

[36]  G L Gerstein,et al.  Mutual temporal relationships among neuronal spike trains. Statistical techniques for display and analysis. , 1972, Biophysical journal.

[37]  D. Perkel,et al.  Simultaneously Recorded Trains of Action Potentials: Analysis and Functional Interpretation , 1969, Science.

[38]  M S Smith,et al.  A computer program for simulating a network of interacting neurons. II. Programming aspects. , 1976, Computers and biomedical research, an international journal.

[39]  Ian H. Stevenson,et al.  Bayesian Inference of Functional Connectivity and Network Structure From Spikes , 2009, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[40]  G. P. Moore,et al.  Neuronal spike trains and stochastic point processes. I. The single spike train. , 1967, Biophysical journal.

[41]  Roman Borisyuk,et al.  Oscillatory activity in the neural networks of spiking elements. , 2002, Bio Systems.

[42]  D. Brillinger.  Estimation of the Second-Order Intensities of a Bivariate Stationary Point Process , 1976 .

[43]  D H Perkel,et al.  A computer program for simulating a network of interacting neurons. I. Organization and physiological assumptions. , 1976, Computers and biomedical research, an international journal.

[44]  José Carlos Príncipe,et al.  A comparison of binless spike train measures , 2010, Neural Computing and Applications.

[45]  Satish Iyengar,et al.  Parameter estimation for a leaky integrate-and-fire neuronal model from ISI data , 2008, Journal of Computational Neuroscience.

[46]  Sonja Grün,et al.  Frontiers in Computational Neuroscience , 2022 .

[47]  M K Habib,et al.  Dynamics of neuronal firing correlation: modulation of "effective connectivity". , 1989, Journal of neurophysiology.

[48]  Jon H. Kaas,et al.  Statistical analysis of large-scale neuronal recording data , 2010, Neural Networks.

[49]  Eero P. Simoncelli,et al.  Spatio-temporal correlations and visual signalling in a complete neuronal population , 2008, Nature.

[50]  Günther Palm,et al.  Detecting higher-order interactions among the spiking events in a group of neurons , 1995, Biological Cybernetics.

[51]  Sonja Grün,et al.  CuBIC: cumulant based inference of higher-order correlations in massively parallel spike trains , 2009, Journal of Computational Neuroscience.