Generation of Spike Trains with Controlled Auto- and Cross-Correlation Functions

Emerging evidence indicates that information processing, as well as learning and memory processes, in both the network and single-neuron levels are highly dependent on the correlation structure of multiple spike trains. Contemporary experimental as well as theoretical studies that involve quasi-realistic neuronal stimulation thus require a method for controlling spike train correlations. This letter introduces a general new strategy for generating multiple spike trains with exactly controlled mean firing rates and correlation structure (defined in terms of auto- and cross-correlation functions). Our approach nonlinearly transforms random gaussian-distributed processes with a predistorted correlation structure into nonnegative rate processes, which are then used to generate doubly stochastic Poisson point processes with the required correlation structure. We show how this approach can be used to generate stationary or nonstationary spike trains from small or large groups of neurons with diverse auto- and cross-correlation structures. We analyze and derive analytical formulas for the high-order correlation structure of generated spike trains and discuss the limitations of this approach.

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

[2]  G. E. Johnson.  Constructions of particular random processes , 1994, Proc. IEEE.

[3]  Roberto Araya,et al.  Dendritic spines linearize the summation of excitatory potentials , 2006, Proceedings of the National Academy of Sciences.

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

[5]  L. Isserlis.  ON A FORMULA FOR THE PRODUCT-MOMENT COEFFICIENT OF ANY ORDER OF A NORMAL FREQUENCY DISTRIBUTION IN ANY NUMBER OF VARIABLES , 1918 .

[6]  Ernst Niebur,et al.  Generation of Synthetic Spike Trains with Defined Pairwise Correlations , 2007, Neural Computation.

[7]  J. Magee,et al.  State-Dependent Dendritic Computation in Hippocampal CA1 Pyramidal Neurons , 2006, The Journal of Neuroscience.

[8]  Stefan Rotter,et al.  Dependence of Neuronal Correlations on Filter Characteristics and Marginal Spike Train Statistics , 2008, Neural Computation.

[9]  S. Rosenbaum.  Moments of a Truncated Bivariate Normal Distribution , 1961 .

[10]  L. Paninski,et al.  Superlinear Population Encoding of Dynamic Hand Trajectory in Primary Motor Cortex , 2004, The Journal of Neuroscience.

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

[12]  D. Cox.  Some Statistical Methods Connected with Series of Events , 1955 .

[13]  Sonja Grün,et al.  Can Spike Coordination Be Differentiated from Rate Covariation? , 2008, Neural Computation.

[14]  C. K. Knox,et al.  Signal transmission in random spike trains with applications to the statocyst neurons of the lobster , 1970, Kybernetik.

[15]  Moshe Abeles,et al.  Corticonics: Neural Circuits of Cerebral Cortex , 1991 .

[16]  Romain Brette,et al.  Generation of Correlated Spike Trains , 2009, Neural Computation.

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

[18]  Peter Dayan,et al.  Theoretical Neuroscience: Computational and Mathematical Modeling of Neural Systems , 2001 .

[19]  A. Aertsen,et al.  Controlling synaptic input patterns in vitro by dynamic photo stimulation. , 2005, Journal of neurophysiology.

[20]  Bartlett W. Mel,et al.  Information Processing in Dendritic Trees , 1994, Neural Computation.

[21]  Donald L. Snyder,et al.  Random Point Processes in Time and Space , 1991 .

[22]  Don H. Johnson,et al.  Point process models of single-neuron discharges , 1996, Journal of Computational Neuroscience.

[23]  I Segev,et al.  Untangling dendrites with quantitative models. , 2000, Science.

[24]  Jonathon Shlens,et al.  The Structure of Multi-Neuron Firing Patterns in Primate Retina , 2006, The Journal of Neuroscience.

[25]  Bartlett W. Mel,et al.  Dendrites: bug or feature? , 2003, Current Opinion in Neurobiology.

[26]  Rafael Yuste,et al.  Two-photon photostimulation and imaging of neural circuits , 2007, Nature Methods.

[27]  Alexander S. Ecker,et al.  Generating Spike Trains with Specified Correlation Coefficients , 2009, Neural Computation.

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

[29]  S. Shoham,et al.  Patterned Optical Activation of Retinal Ganglion Cells , 2007, 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[30]  M. London,et al.  Dendritic computation. , 2005, Annual review of neuroscience.

[31]  Yasushi Miyashita,et al.  Dendritic spine geometry is critical for AMPA receptor expression in hippocampal CA1 pyramidal neurons , 2001, Nature Neuroscience.

[32]  Sonja Grün,et al.  Impact of Higher-Order Correlations on Coincidence Distributions of Massively Parallel Data , 2007, Summer School on Neural Networks.

[33]  Stefan Rotter,et al.  Higher-Order Statistics of Input Ensembles and the Response of Simple Model Neurons , 2003, Neural Computation.

[34]  Barry Simon,et al.  The P(φ)[2] Euclidean (quantum) field theory , 1974 .

[35]  Matthias Bethge,et al.  Near-Maximum Entropy Models for Binary Neural Representations of Natural Images , 2007, NIPS.

[36]  Matthew Fellows,et al.  Statistical encoding model for a primary motor cortical brain-machine interface , 2005, IEEE Transactions on Biomedical Engineering.

[37]  J. Magee,et al.  Integrative Properties of Radial Oblique Dendrites in Hippocampal CA1 Pyramidal Neurons , 2006, Neuron.

[38]  Shy Shoham,et al.  Correlation-distortion based identification of Linear-Nonlinear-Poisson models , 2009, Journal of Computational Neuroscience.

[39]  Michael J. Berry,et al.  Weak pairwise correlations imply strongly correlated network states in a neural population , 2005, Nature.

[40]  Shy Shoham,et al.  Multivariate Autoregressive Modeling and Granger Causality Analysis of Multiple Spike Trains , 2010, Comput. Intell. Neurosci..

[41]  Shy Shoham,et al.  Rapid neurotransmitter uncaging in spatially defined patterns , 2005, Nature Methods.