Functional two-way analysis of variance and bootstrap methods for neural synchrony analysis

BMC Neuroscience

Jorge Mariño¹  Javier Cudeiro¹  Nelson Espinosa¹  Ricardo Cao²  Aldana M González Montoro² 

[1] Neuroscience and Motor Control Group (NEUROcom), Department of Medicine, Facultad de Ciencias de la Salud, Universidade da Coruña, Campus de Oza s/n, 15006 A Coruña, Spain;Department of Mathematics, Facultad de Informática, Universidade da Coruña, Campus de Elviña s/n, 15071 A Coruña, Spain

关键词: Functional data;    Low firing-rate;    Dependence;    Spike-trains;    Bootstrap;    Cross-correlation analysis;   

Others  :  1091434
DOI  :  10.1186/1471-2202-15-96

 received in 2013-12-18, accepted in 2014-07-31,  发布年份 2014

【 摘 要 】

Background

Pairwise association between neurons is a key feature in understanding neural coding. Statistical neuroscience provides tools to estimate and assess these associations. In the mammalian brain, activating ascending pathways arise from neuronal nuclei located at the brainstem and at the basal forebrain that regulate the transition between sleep and awake neuronal firing modes in extensive regions of the cerebral cortex, including the primary visual cortex, where neurons are known to be selective for the orientation of a given stimulus. In this paper, the estimation of neural synchrony as a function of time is studied in data obtained from anesthetized cats. A functional data analysis of variance model is proposed. Bootstrap statistical tests are introduced in this context; they are useful tools for the study of differences in synchrony strength regarding 1) transition between different states (anesthesia and awake), and 2) affinity given by orientation selectivity.

Results

An analysis of variance model for functional data is proposed for neural synchrony curves, estimated with a cross-correlation based method. Dependence arising from the experimental setting needs to be accounted for. Bootstrap tests allow the identification of differences between experimental conditions (modes of activity) and between pairs of neurons formed by cells with different affinities given by their preferred orientations. In our test case, interactions between experimental conditions and preferred orientations are not statistically significant.

Conclusions

The results reflect the effect of different experimental conditions, as well as the affinity regarding orientation selectivity in neural synchrony and, therefore, in neural coding. A cross-correlation based method is proposed that works well under low firing activity. Functional data statistical tools produce results that are useful in this context. Dependence is shown to be necessary to account for, and bootstrap tests are an appropriate method with which to do so.

【 授权许可】

   
2014 González Montoro et al.; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

【 参考文献 】

• [1]Kuhn A, Aertsen A, Rotter S: Higher-order statistics of input ensembles and the response of simple model neurons. Neural Comput 2003, 15:67-101.
• [2]Steinmetz PN, Roy A, Fitzgerald PJ, Hsiao SS, Johnson KO, Niebur E: Attention modulates synchronized neuronal firing in primate somatosensory cortex. Nature 2000, 404:187-190.
• [3]Stopfer M, Bhagavan S, Smith BH, Laurent G: Impaired odour discrimination on desynchronization of odour-encoding neural assemblies. Nature 1997, 390:70-74.
• [4]Vaadia E, Haalman I, Abeles M, Bergman H, Prut Y, Slovin H, Aertsen A: Dynamics of neuronal interactions in monkey cortex in relation to behavioral events. Nature 1995, 373:515-518.
• [5]Zohary E, Shadlen M, Newsome W: Correlated neuronal discharge rate and its implications for psychophysical performance. Nature 1994, 370:140-143.
• [6]Gerstein G, Perkel D: Simultaneously recorded trains of action potentials - Analysis and functional interpretation. Science 1969, 164:828-830.
• [7]Grün S, Diesmann M, Aertsen A: Unitary events in multiple single-neuron spiking activity: I. Detection and significance. Neural Comput 2002, 14:43-80.
• [8]Grün S, Diesmann M, Aertsen A: Unitary events in multiple single-neuron spiking activity: II. Nonstationary Data. Neural Comput 2002, 14:81-119.
• [9]Faes C, Geys H, Molenberghs G, Aerts M, Cadarso-Suárez C, Acuña C, Cano M: A flexible method to measure synchrony in neuronal firing. J Am Stat Assoc 2008, 103:149-261.
• [10]González-Montoro AM, Cao R, Faes C, Molenberghs G, Espinosa N, Cudeiro J, Mariño J: Cross nearest-spike interval based method to measure synchrony dynamics. Math Biosci Eng 2014. doi:10.3934/mbe.2014.11.27
• [11]Grün S: Data-driven significance estimation for precise spike correlation. J Neurophysiol 2009, 101:1126-1140.
• [12]Ventura V, Cai C, Kass RE: Statistical assessment of time-varying dependence between two neurons. J Neurophysiol 2005, 94:2940-2947.
• [13]Steriade M: Sleep oscillations and their blockage by activating systems. J Psychiat Neurosci 1994, 19:354-358.
• [14]Steriade M, Gloor P, Llinás RR, Lopes da Silva FH, Mesulam MM: Basic mechanisms of cerebral rhythmic activities. Electroencephalogr Clin Neurophysiol 1990, 76:481-508.
• [15]Brown EN, Lydic R, Schiff D: General anesthesia, sleep, and coma. N Engl J Med 2010, 363:2638-2650.
• [16]Franks NP: General anaesthesia: from molecular targets to neuronal pathways of sleep and arousal. Net Rev Neurosci 2008, 9:370-386.
• [17]Sleigh JW, Voss L, Steyn-Ross ML, Steyn-Ross DA, Wilson MT: Modelling sleep and general anaesthesia. Springer Ser Comp Neurosci 2011, 15:21-41.
• [18]Mariño J, Cudeiro J: Nitric oxide-mediated cortical activation: A diffuse wake-up system. J Neurosci 2003, 23:4299-4307.
• [19]Hubel DH, Wiesel TN: Receptive fields, binocular interaction and functional architecture in the cat’ s visual cortex. J Physiol 1962, 10:106-154.
• [20]Cuesta-Albertos J, Febrero-Bande M: A simple multiway ANOVA for functional data. Test 2010, 19:537-557.
• [21]Perkel DH, Gerstein GL, Moore GP: Neuronal spike trains and stochastic processes. II. Simultaneous spike trains. Biophys J 1967, 7:419-440.
• [22]Ramsay JO, Silverman BW: Applied Functional Data Analysis: Methods and Case Studies. New York: Springer; 2002.
• [23]Ramsay JO, Silverman BW: Functional Data Analysis. New York: Springer Series in Statistics; 2005.
• [24]Ferraty F, Vieu P: Nonparametric Functional Data Analysis: Theory and Practice. New York: Springer; 2006.
• [25]Cuesta-Albertos JA, Fraiman R, Ransford T: A sharp form of the Cramer-Wold theorem. J Theor Prob 2007, 20:201-209.
• [26]Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B 1995, 57:289-300.
• [27]Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Ann Stat 2001, 29:1165-1188.
• [28]Provost SB, Ridiuk EM: The exact distribution of indefinite quadratic forms iin noncentral normal vectors. Ann Inst Statist Math 1996, 48:381-394.
• [29]Fuller PM, Sherman D, Pedersen NP, Saper CB, Lu J: Reassessment of the structural basis of the ascending arousal system. J Comp Neurol 2011, 519:933-956.
• [30]Steriade M: Arousal: revisiting the reticular activating system. Science 1996, 272:225-226.
• [31]Henny P, Jones BE: Projections from basal forebrain to prefrontal cortex comprise cholinergic, GABAergic and glutamatergic inputs to pyramidal cells or interneurons. Eur J Neurosci 2008, 27:654-670.
• [32]Schummers J, Mariño J, Sur M: Synaptic integration by V1 neurons depends on location within the orientation map. Neuron 2002, 36:969-978.
• [33]Mariño J, Schummers J, Lyon DC, Schwabe L, Beck O, Wiesing P, Obermayer K, Sur M: Invariant computations in local cortical networks with balanced excitation and inhibition. Nat Neurocsi 2005, 8:194-201.
• [34]Stimberg M, Wimmer K, Martin R, Schwabe L, Mariño J, Schummers J, Lyon DC, Sur M, Obermayer K: The operating regime of local computations in primary visual cortex. Cereb Cortex 2009, 19:2166-2180.