Functional connectomics from neural dynamics: probabilistic graphical models for neuronal network of Caenorhabditis elegans

We propose an approach to represent neuronal network dynamics as a probabilistic graphical model (PGM). To construct the PGM, we collect time series of neuronal responses produced by the neuronal network and use singular value decomposition to obtain a low-dimensional projection of the time-series data. We then extract dominant patterns from the projections to get pairwise dependency information and create a graphical model for the full network. The outcome model is a functional connectome that captures how stimuli propagate through the network and thus represents causal dependencies between neurons and stimuli. We apply our methodology to a model of the Caenorhabditis elegans somatic nervous system to validate and show an example of our approach. The structure and dynamics of the C. elegans nervous system are well studied and a model that generates neuronal responses is available. The resulting PGM enables us to obtain and verify underlying neuronal pathways for known behavioural scenarios and detect possible pathways for novel scenarios. This article is part of a discussion meeting issue ‘Connectome to behaviour: modelling C. elegans at cellular resolution’.

[1]  S. R. Wicks,et al.  A Dynamic Network Simulation of the Nematode Tap Withdrawal Circuit: Predictions Concerning Synaptic Function Using Behavioral Criteria , 1996, The Journal of Neuroscience.

[2]  Junwei Han,et al.  Inferring functional interaction and transition patterns via dynamic bayesian variable partition models , 2013, Human brain mapping.

[3]  Eli Shlizerman,et al.  Classification of Fixed Point Network Dynamics from Multiple Node Timeseries Data , 2017, Front. Neuroinform..

[4]  Theodore H. Lindsay,et al.  Global Brain Dynamics Embed the Motor Command Sequence of Caenorhabditis elegans , 2015, Cell.

[5]  S. Lockery,et al.  Active Currents Regulate Sensitivity and Dynamic Range in C. elegans Neurons , 1998, Neuron.

[6]  Selin Aviyente,et al.  Recursive Tensor Subspace Tracking for Dynamic Brain Network Analysis , 2017, IEEE Transactions on Signal and Information Processing over Networks.

[7]  Le Song,et al.  Estimating time-varying networks , 2008, ISMB 2008.

[8]  Adam A. Margolin,et al.  Reverse engineering of regulatory networks in human B cells , 2005, Nature Genetics.

[9]  Nir Friedman,et al.  Probabilistic Graphical Models - Principles and Techniques , 2009 .

[10]  Lav R. Varshney,et al.  Structural Properties of the Caenorhabditis elegans Neuronal Network , 2009, PLoS Comput. Biol..

[11]  Andrew M. Leifer,et al.  Monoaminergic Orchestration of Motor Programs in a Complex C. elegans Behavior , 2013, PLoS biology.

[12]  Gal Haspel,et al.  A connectivity model for the locomotor network of Caenorhabditis elegans , 2012, Worm.

[13]  Michal Linial,et al.  Using Bayesian Networks to Analyze Expression Data , 2000, J. Comput. Biol..

[14]  Zhaoyang Feng,et al.  Light-sensitive neurons and channels mediate phototaxis in C. elegans , 2008, Nature Neuroscience.

[15]  Cornelia I Bargmann,et al.  A Distributed Chemosensory Circuit for Oxygen Preference in C. elegans , 2006, PLoS biology.

[16]  A. Kullyev,et al.  Regulation of Extrasynaptic 5-HT by Serotonin Reuptake Transporter Function in 5-HT-Absorbing Neurons Underscores Adaptation Behavior in Caenorhabditis elegans , 2011, The Journal of Neuroscience.

[17]  Eli Shlizerman,et al.  Neural Interactome: Interactive Simulation of a Neuronal System , 2017 .

[18]  Joshua W Shaevitz,et al.  Whole-brain calcium imaging with cellular resolution in freely behaving Caenorhabditis elegans , 2015, Proceedings of the National Academy of Sciences.

[19]  G. Stewart Perturbation theory for the singular value decomposition , 1990 .

[20]  Trevor Hastie,et al.  The Elements of Statistical Learning , 2001 .

[21]  I S Kohane,et al.  Mutual information relevance networks: functional genomic clustering using pairwise entropy measurements. , 1999, Pacific Symposium on Biocomputing. Pacific Symposium on Biocomputing.

[22]  Yuichi Nakamura,et al.  Approximation of dynamical systems by continuous time recurrent neural networks , 1993, Neural Networks.

[23]  J. Nathan Kutz,et al.  Multistability and Long-Timescale Transients Encoded by Network Structure in a Model of C. elegans Connectome Dynamics , 2017, Front. Comput. Neurosci..

[24]  J. Nathan Kutz,et al.  Neural Activity Measures and Their Dynamics , 2012, SIAM J. Appl. Math..

[25]  Marina Vannucci,et al.  Bayesian models for functional magnetic resonance imaging data analysis , 2014 .

[26]  Greg J. Stephens,et al.  Dimensionality and Dynamics in the Behavior of C. elegans , 2007, PLoS Comput. Biol..

[27]  Ben Taskar,et al.  Graphical Models in a Nutshell , 2007 .

[28]  Karl Deisseroth,et al.  Optogenetics in Neural Systems , 2011, Neuron.

[29]  Navin Pokala,et al.  Neurons Detect Increases and Decreases in Oxygen Levels Using Distinct Guanylate Cyclases , 2009, Neuron.

[30]  I. Mori,et al.  Neural regulation of thermotaxis in Caenorhabditis elegans , 1995, Nature.

[31]  Cori Bargmann,et al.  A circuit for navigation in Caenorhabditis elegans , 2005 .

[32]  Kevin P. Murphy,et al.  Machine learning - a probabilistic perspective , 2012, Adaptive computation and machine learning series.

[33]  S. Brenner,et al.  The neural circuit for touch sensitivity in Caenorhabditis elegans , 1985, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[34]  Eli Shlizerman,et al.  Low-dimensional functionality of complex network dynamics: neurosensory integration in the Caenorhabditis Elegans connectome. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[35]  Olaf Sporns,et al.  Network structure of cerebral cortex shapes functional connectivity on multiple time scales , 2007, Proceedings of the National Academy of Sciences.

[36]  Amr Ahmed,et al.  Recovering time-varying networks of dependencies in social and biological studies , 2009, Proceedings of the National Academy of Sciences.