Higher-Order Cumulants Drive Neuronal Activity Patterns, Inducing UP-DOWN States in Neural Populations

A major challenge in neuroscience is to understand the role of the higher-order correlations structure of neuronal populations. The dichotomized Gaussian model (DG) generates spike trains by means of thresholding a multivariate Gaussian random variable. The DG inputs are Gaussian distributed, and thus have no interactions beyond the second order in their inputs; however, they can induce higher-order correlations in the outputs. We propose a combination of analytical and numerical techniques to estimate higher-order, above the second, cumulants of the firing probability distributions. Our findings show that a large amount of pairwise interactions in the inputs can induce the system into two possible regimes, one with low activity (“DOWN state”) and another one with high activity (“UP state”), and the appearance of these states is due to a combination between the third- and fourth-order cumulant. This could be part of a mechanism that would help the neural code to upgrade specific information about the stimuli, motivating us to examine the behavior of the critical fluctuations through the Binder cumulant close to the critical point. We show, using the Binder cumulant, that higher-order correlations in the outputs generate a critical neural system that portrays a second-order phase transition.

[1]  Kurt Binder,et al.  Critical Properties from Monte Carlo Coarse Graining and Renormalization , 1981 .

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

[3]  Fernando Montani,et al.  Quantifying higher-order correlations in a neuronal pool , 2015 .

[4]  Shun-ichi Amari,et al.  Methods of information geometry , 2000 .

[5]  Eric Shea-Brown,et al.  A Moment-Based Maximum Entropy Model for Fitting Higher-Order Interactions in Neural Data , 2018, bioRxiv.

[6]  B. Roy Frieden,et al.  Science from Fisher Information: A Unification , 2004 .

[7]  Silvia Scarpetta,et al.  Alternation of up and down states at a dynamical phase-transition of a neural network with spatiotemporal attractors , 2014, Front. Syst. Neurosci..

[8]  W. Bialek,et al.  Are Biological Systems Poised at Criticality? , 2010, 1012.2242.

[9]  Ifije E. Ohiorhenuan,et al.  Sparse coding and high-order correlations in fine-scale cortical networks , 2010, Nature.

[10]  Jaime de la Rocha,et al.  Supplementary Information for the article ‘ Correlation between neural spike trains increases with firing rate ’ , 2007 .

[11]  X. R. Wang,et al.  Fisher transfer entropy: quantifying the gain in transient sensitivity , 2015, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[12]  Yutaka Sakai,et al.  Synchronous Firing and Higher-Order Interactions in Neuron Pool , 2003, Neural Computation.

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

[14]  Sheila Nirenberg,et al.  Decoding neuronal spike trains: How important are correlations? , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[15]  G. Sumbre,et al.  Whole-Brain Neuronal Activity Displays Crackling Noise Dynamics , 2018, Neuron.

[16]  Fernando Montani,et al.  Common inputs in subthreshold membrane potential: The role of quiescent states in neuronal activity. , 2018, Physical review. E.

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

[18]  Kurt Binder,et al.  Finite size scaling analysis of ising model block distribution functions , 1981 .

[19]  Stefano Panzeri,et al.  The impact of high-order interactions on the rate of synchronous discharge and information transmission in somatosensory cortex , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[20]  Peter E. Latham,et al.  Pairwise Maximum Entropy Models for Studying Large Biological Systems: When They Can Work and When They Can't , 2008, PLoS Comput. Biol..

[21]  Stefan Rotter,et al.  Higher-order correlations in non-stationary parallel spike trains : statistical modeling and inference , 2022 .

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

[23]  Adam Kohn,et al.  Laminar dependence of neuronal correlations in visual cortex. , 2013, Journal of neurophysiology.

[24]  C. Wilson,et al.  Spontaneous firing patterns and axonal projections of single corticostriatal neurons in the rat medial agranular cortex. , 1994, Journal of neurophysiology.

[25]  Shun-ichi Amari,et al.  Information-Geometric Measure for Neural Spikes , 2002, Neural Computation.

[26]  K. Binder,et al.  Wetting transition in the two-dimensional Blume-Capel model: a Monte Carlo study. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  M. Bethge,et al.  Common input explains higher-order correlations and entropy in a simple model of neural population activity. , 2011, Physical review letters.

[28]  D Contreras,et al.  Mechanisms of long‐lasting hyperpolarizations underlying slow sleep oscillations in cat corticothalamic networks. , 1996, The Journal of physiology.

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

[30]  Shan Yu,et al.  Higher-Order Interactions Characterized in Cortical Activity , 2011, The Journal of Neuroscience.

[31]  Serena Bradde,et al.  PCA Meets RG , 2016, Journal of Statistical Physics.

[32]  Luciano Fadiga,et al.  Superlinear Summation of Information in Premotor Neuron Pairs , 2017, Int. J. Neural Syst..

[33]  Sang Joon Kim,et al.  A Mathematical Theory of Communication , 2006 .

[34]  Yuji Ikegaya,et al.  Simultaneous silence organizes structured higher-order interactions in neural populations , 2015, Scientific Reports.

[35]  Fernando Montani,et al.  Higher-order correlations in common input shapes the output spiking activity of a neural population , 2017 .

[36]  Simon R. Schultz,et al.  Statistical modelling of higher-order correlations in pools of neural activity , 2012, 1211.6348.

[37]  A. Destexhe,et al.  Are corticothalamic ‘up’ states fragments of wakefulness? , 2007, Trends in Neurosciences.

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

[39]  Michael J. Berry,et al.  Ising models for networks of real neurons , 2006, q-bio/0611072.

[40]  M. A. Smith,et al.  The Role of Correlations in Direction and Contrast Coding in the Primary Visual Cortex , 2007, The Journal of Neuroscience.

[41]  R. Segev,et al.  Sparse low-order interaction network underlies a highly correlated and learnable neural population code , 2011, Proceedings of the National Academy of Sciences.

[42]  Charles J. Wilson,et al.  The origins of two-state spontaneous membrane potential fluctuations of neostriatal spiny neurons , 1996, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[43]  X. R. Wang,et al.  Relating Fisher information to order parameters. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[44]  R. Fisher,et al.  On the Mathematical Foundations of Theoretical Statistics , 1922 .

[45]  A. Destexhe,et al.  The high-conductance state of neocortical neurons in vivo , 2003, Nature Reviews Neuroscience.

[46]  Eric Shea-Brown,et al.  Triplet correlations among similarly tuned cells impact population coding , 2015, Front. Comput. Neurosci..

[47]  Jonathan D. Victor,et al.  Information-geometric measure of 3-neuron firing patterns characterizes scale-dependence in cortical networks , 2011, Journal of Computational Neuroscience.

[48]  Jorge Stolfi,et al.  Phase transitions and self-organized criticality in networks of stochastic spiking neurons , 2016, Scientific Reports.