Differential Covariance: A New Class of Methods to Estimate Sparse Connectivity from Neural Recordings

With our ability to record more neurons simultaneously, making sense of these data is a challenge. Functional connectivity is one popular way to study the relationship of multiple neural signals. Correlation-based methods are a set of currently well-used techniques for functional connectivity estimation. However, due to explaining away and unobserved common inputs (Stevenson, Rebesco, Miller, & Körding, 2008), they produce spurious connections. The general linear model (GLM), which models spike trains as Poisson processes (Okatan, Wilson, & Brown, 2005; Truccolo, Eden, Fellows, Donoghue, & Brown, 2005; Pillow et al., 2008), avoids these confounds. We develop here a new class of methods by using differential signals based on simulated intracellular voltage recordings. It is equivalent to a regularized AR(2) model. We also expand the method to simulated local field potential recordings and calcium imaging. In all of our simulated data, the differential covariance-based methods achieved performance better than or similar to the GLM method and required fewer data samples. This new class of methods provides alternative ways to analyze neural signals.

[1]  A. McIntosh,et al.  Structural modeling of functional neural pathways mapped with 2-deoxyglucose: effects of acoustic startle habituation on the auditory system , 1991, Brain Research.

[2]  Karl J. Friston Functional and Effective Connectivity: A Review , 2011, Brain Connect..

[3]  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.

[4]  Maxime Bonjean,et al.  Corticothalamic Feedback Controls Sleep Spindle Duration In Vivo , 2011, The Journal of Neuroscience.

[5]  Haiping Huang Sparse Hopfield network reconstruction with ℓ1 regularization , 2012 .

[6]  E. Kandel,et al.  Neuroscience thinks big (and collaboratively) , 2013, Nature Reviews Neuroscience.

[7]  Xiuming Shan,et al.  Covariances of linear stochastic differential equations for analyzing computer networks , 2011 .

[8]  Alexander S. Ecker,et al.  Improved Estimation and Interpretation of Correlations in Neural Circuits , 2015, PLoS Comput. Biol..

[9]  L. Ricciardi,et al.  The Ornstein-Uhlenbeck process as a model for neuronal activity , 1979, Biological Cybernetics.

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

[11]  Karl J. Friston,et al.  Multivariate Autoregressive Modelling of fMRI time series , 2003 .

[12]  Theo Geisel,et al.  Model-Free Reconstruction of Excitatory Neuronal Connectivity from Calcium Imaging Signals , 2012, PLoS Comput. Biol..

[13]  Nanny Wermuth,et al.  Multivariate Dependencies: Models, Analysis and Interpretation , 1996 .

[14]  Brendon O. Watson,et al.  Spike inference from calcium imaging using sequential Monte Carlo methods. , 2009, Biophysical journal.

[15]  Konrad P Kording,et al.  How advances in neural recording affect data analysis , 2011, Nature Neuroscience.

[16]  Pablo A. Parrilo,et al.  Rank-Sparsity Incoherence for Matrix Decomposition , 2009, SIAM J. Optim..

[17]  T. Sejnowski,et al.  Origin of slow cortical oscillations in deafferented cortical slabs. , 2000, Cerebral cortex.

[18]  Petr Lánský,et al.  Ornstein-Uhlenbeck model neuron revisited , 1995, Biological Cybernetics.

[19]  R. Mazo On the theory of brownian motion , 1973 .

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

[21]  T. Sejnowski,et al.  Model of Thalamocortical Slow-Wave Sleep Oscillations and Transitions to Activated States , 2002, The Journal of Neuroscience.

[22]  John Hertz,et al.  Belief-Propagation and replicas for inference and learning in a kinetic Ising model with hidden spins , 2014, 1412.1727.

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

[24]  I. Timofeev,et al.  Interneuron‐mediated inhibition synchronizes neuronal activity during slow oscillation , 2012, The Journal of physiology.

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

[26]  Shaoqun Zeng,et al.  Method to reconstruct neuronal action potential train from two-photon calcium imaging. , 2010, Journal of biomedical optics.

[27]  Stefan Rotter,et al.  Reconstruction of sparse connectivity in neural networks from spike train covariances , 2013 .

[28]  Federico Ricci-Tersenghi,et al.  Inferring Synaptic Structure in Presence of Neural Interaction Time Scales , 2014, PloS one.

[29]  Yasser Roudi,et al.  Learning and inference in a nonequilibrium Ising model with hidden nodes. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[30]  R. Tibshirani,et al.  Sparse inverse covariance estimation with the graphical lasso. , 2008, Biostatistics.

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

[32]  Yasser Roudi,et al.  Multi-neuronal activity and functional connectivity in cell assemblies , 2015, Current Opinion in Neurobiology.

[33]  Alexandre d'Aspremont,et al.  Convex optimization techniques for fitting sparse Gaussian graphical models , 2006, ICML.

[34]  Karin Schwab,et al.  Comparison of linear signal processing techniques to infer directed interactions in multivariate neural systems , 2005, Signal Process..

[35]  Zhixun Su,et al.  Linearized Alternating Direction Method with Adaptive Penalty for Low-Rank Representation , 2011, NIPS.

[36]  J. Hertz,et al.  Mean field theory for nonequilibrium network reconstruction. , 2010, Physical review letters.

[37]  P. Nunez,et al.  Electric fields of the brain , 1981 .

[38]  Stefan J. Kiebel,et al.  Inferring Neuronal Dynamics from Calcium Imaging Data Using Biophysical Models and Bayesian Inference , 2016, PLoS Comput. Biol..

[39]  Anthony N. Burkitt,et al.  A Review of the Integrate-and-fire Neuron Model: I. Homogeneous Synaptic Input , 2006, Biological Cybernetics.

[40]  A. Destexhe Spike-and-Wave Oscillations Based on the Properties of GABAB Receptors , 1998, The Journal of Neuroscience.