Improved functional connectivity network estimation for brain networks using multivariate partial coherence

OBJECTIVE Graphical networks and network metrics are widely used to understand and characterise brain networks and brain function. These methods can be applied to a range of electrophysiological data including electroencephalography, local field potential and single unit recordings. Functional networks are often constructed using pair-wise correlation between variables. The objective of this study is to demonstrate that functional networks can be more accurately estimated using partial correlation than with pair-wise correlation. APPROACH We compared network metrics derived from unconditional and conditional graphical networks, obtained using coherence and multivariate partial coherence (MVPC), respectively. Graphical networks were constructed using coherence and MVPC estimates, and binary and weighted network metrics derived from these: node degree, path length, clustering coefficients and small-world index. MAIN RESULTS Network metrics were applied to simulated and experimental single unit spike train data. Simulated data used a 10x10 grid of simulated cortical neurons with centre-surround connectivity. Conditional network metrics gave a more accurate representation of the known connectivity: Numbers of excitatory connections had range 3-11, unconditional binary node degree had range 6-80, conditional node degree had range 2-13. Experimental data used multi-electrode array recording with 19 single-units from left and right hippocampal brain areas in a rat model for epilepsy. Conditional network analysis showed similar trends to simulated data, with lower binary node degree and longer binary path lengths compared to unconditional networks. SIGNIFICANCE We conclude that conditional networks, where common dependencies are removed through partial coherence analysis, give a more accurate representation of the interactions in a graphical network model. These results have important implications for graphical network analyses of brain networks and suggest that functional networks should be derived using partial correlation, based on MVPC estimates, as opposed to the common approach of pair-wise correlation.

[1]  Marcus Kaiser,et al.  A tutorial in connectome analysis: Topological and spatial features of brain networks , 2011, NeuroImage.

[2]  Javier M. Buldú,et al.  Beware of the Small-World Neuroscientist! , 2016, Front. Hum. Neurosci..

[3]  J. R. Rosenberg,et al.  Identification of directed interactions in networks , 2011, Biological Cybernetics.

[4]  J. R. Rosenberg,et al.  The Fourier approach to the identification of functional coupling between neuronal spike trains. , 1989, Progress in biophysics and molecular biology.

[5]  Michael Eichler,et al.  A graphical approach for evaluating effective connectivity in neural systems , 2005, Philosophical Transactions of the Royal Society B: Biological Sciences.

[6]  Lubos Buzna,et al.  Robustness of trans-European gas networks. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[7]  Lin Song,et al.  Comparison of co-expression measures: mutual information, correlation, and model based indices , 2012, BMC Bioinformatics.

[8]  Duncan J. Watts,et al.  Collective dynamics of ‘small-world’ networks , 1998, Nature.

[9]  Andrew T. Walden,et al.  Partial Coherence Estimation via Spectral Matrix Shrinkage under Quadratic Loss , 2015, IEEE Transactions on Signal Processing.

[10]  Guy Karlebach,et al.  Modelling and analysis of gene regulatory networks , 2008, Nature Reviews Molecular Cell Biology.

[11]  Carsten O. Daub,et al.  The mutual information: Detecting and evaluating dependencies between variables , 2002, ECCB.

[12]  Hiroaki Kitano,et al.  Gene regulatory network analysis supports inflammation as a key neurodegeneration process in prion disease , 2012, BMC Systems Biology.

[13]  A. T. Walden,et al.  Graphical modelling for brain connectivity via partial coherence , 2009, Journal of Neuroscience Methods.

[14]  A M Amjad,et al.  A framework for the analysis of mixed time series/point process data--theory and application to the study of physiological tremor, single motor unit discharges and electromyograms. , 1995, Progress in biophysics and molecular biology.

[15]  E. Bullmore,et al.  Undirected graphs of frequency-dependent functional connectivity in whole brain networks , 2005, Philosophical Transactions of the Royal Society B: Biological Sciences.

[16]  K. Kaski,et al.  Intensity and coherence of motifs in weighted complex networks. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  O. Sporns,et al.  Motifs in Brain Networks , 2004, PLoS biology.

[18]  R. Mason,et al.  Inhibition of endocannabinoid metabolism attenuates enhanced hippocampal neuronal activity induced by kainic acid , 2008, Synapse.

[19]  Vince D. Calhoun,et al.  Application of Graph Theory to Assess Static and Dynamic Brain Connectivity: Approaches for Building Brain Graphs , 2018, Proceedings of the IEEE.

[20]  Bin He,et al.  A weighted small world network measure for assessing functional connectivity , 2013, Journal of Neuroscience Methods.

[21]  David R. Brillinger,et al.  Mutual information in the frequency domain , 2007 .

[22]  Andrew T. Walden,et al.  Constructing brain connectivity group graphs from EEG time series , 2018, Journal of Applied Statistics.

[23]  Vadas Gintautas,et al.  Identification of functional information subgraphs in complex networks. , 2007, Physical review letters.

[24]  Edward T. Bullmore,et al.  On the use of correlation as a measure of network connectivity , 2012, NeuroImage.

[25]  R. Dahlhaus Graphical interaction models for multivariate time series1 , 2000 .

[26]  Guillaume A. Rousselet,et al.  A statistical framework for neuroimaging data analysis based on mutual information estimated via a gaussian copula , 2016, bioRxiv.

[27]  K. Gurney,et al.  Network ‘Small-World-Ness’: A Quantitative Method for Determining Canonical Network Equivalence , 2008, PloS one.

[28]  P. Welch The use of fast Fourier transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms , 1967 .

[29]  Partha P. Mitra,et al.  Sampling Properties of the Spectrum and Coherency of Sequences of Action Potentials , 2000, Neural Computation.

[30]  Jens Timmer,et al.  Networks: On the relation of bi- and multivariate measures , 2015, Scientific Reports.

[31]  Olaf Sporns,et al.  Complex network measures of brain connectivity: Uses and interpretations , 2010, NeuroImage.

[32]  O. Sporns,et al.  Complex brain networks: graph theoretical analysis of structural and functional systems , 2009, Nature Reviews Neuroscience.

[33]  Rainer Dahlhaus,et al.  Partial correlation analysis for the identification of synaptic connections , 2003, Biological Cybernetics.

[34]  Huafu Chen,et al.  Altered Functional Connectivity and Small-World in Mesial Temporal Lobe Epilepsy , 2010, PloS one.

[35]  J. R. Rosenberg,et al.  Identification of patterns of neuronal connectivity—partial spectra, partial coherence, and neuronal interactions , 1998, Journal of Neuroscience Methods.

[36]  David M. Halliday,et al.  Multivariate partial coherence analysis for identification of neuronal connectivity from multiple electrode array recordings , 2014, 2014 IEEE Conference on Biomedical Engineering and Sciences (IECBES).