Inferring causal connectivity from pairwise recordings and optogenetics

To study how the brain works, it is crucial to identify causal interactions between neurons, which is thought to require perturbations. However, when using optogenetics we typically perturb multiple neurons, producing a confound - any of the stimulated neurons can have affected the postsynaptic neuron. Here we show how this produces large biases, and how they can be reduced using the instrumental variable (IV) technique from econometrics. The interaction between stimulation and the absolute refractory period produces a weak, approximately random signal which can be exploited to estimate causal connectivity. When simulating integrate-and-fire neurons, we find that estimates from IV are better than naïve techniques (R2 = 0.77 vs R2 = 0.01). The difference is important as the estimates disagree when applied to experimental data from stimulated neurons with recorded spiking activity. Presented is a robust analysis framework for mapping out network connectivity based on causal neuron interactions.

[1]  Feng Zhang,et al.  An optical neural interface: in vivo control of rodent motor cortex with integrated fiberoptic and optogenetic technology , 2007, Journal of neural engineering.

[2]  Konrad Paul Kording,et al.  Spiking allows neurons to estimate their causal effect , 2018, bioRxiv.

[3]  G. Ermentrout,et al.  Reliability, synchrony and noise , 2008, Trends in Neurosciences.

[4]  L. Wasserman All of Nonparametric Statistics , 2005 .

[5]  Konrad P. Kording,et al.  Neural spiking for causal inference , 2019 .

[6]  Laurence Aitchison,et al.  With or without you: predictive coding and Bayesian inference in the brain , 2017, Current Opinion in Neurobiology.

[7]  Ryan R. Wick,et al.  Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads , 2016, bioRxiv.

[8]  Marc-Oliver Gewaltig,et al.  NEST (NEural Simulation Tool) , 2007, Scholarpedia.

[9]  Konrad Paul Kording,et al.  Could a Neuroscientist Understand a Microprocessor? , 2016, bioRxiv.

[10]  Arvind Kumar,et al.  The High-Conductance State of Cortical Networks , 2008, Neural Computation.

[11]  Illtyd Trethowan Causality , 1938 .

[12]  S. Wright,et al.  The Theory of Path Coefficients a Reply to Niles's Criticism. , 1923, Genetics.

[13]  A. Wald The Fitting of Straight Lines if Both Variables are Subject to Error , 1940 .

[14]  M. J. Friedlander,et al.  The time course and amplitude of EPSPs evoked at synapses between pairs of CA3/CA1 neurons in the hippocampal slice , 1990, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[15]  J. Hertz,et al.  Ising model for neural data: model quality and approximate methods for extracting functional connectivity. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[16]  Bernhard Schölkopf,et al.  Elements of Causal Inference: Foundations and Learning Algorithms , 2017 .

[17]  Alberto Abadie Semiparametric Difference-in-Differences Estimators , 2005 .

[18]  A. Destexhe,et al.  Impact of spontaneous synaptic activity on the resting properties of cat neocortical pyramidal neurons In vivo. , 1998, Journal of neurophysiology.

[19]  W. C. Hall,et al.  High-speed mapping of synaptic connectivity using photostimulation in Channelrhodopsin-2 transgenic mice , 2007, Proceedings of the National Academy of Sciences.

[20]  D. A. Kenny,et al.  Correlation and Causation. , 1982 .

[21]  B. Zemelman,et al.  Selective Photostimulation of Genetically ChARGed Neurons , 2002, Neuron.

[22]  D. Pinault,et al.  A novel single-cell staining procedure performed in vivo under electrophysiological control: morpho-functional features of juxtacellularly labeled thalamic cells and other central neurons with biocytin or Neurobiotin , 1996, Journal of Neuroscience Methods.

[23]  D. A. Kenny,et al.  Correlation and Causation , 1937, Wilmott.

[24]  Skipper Seabold,et al.  Statsmodels: Econometric and Statistical Modeling with Python , 2010, SciPy.

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

[26]  G. Buzsáki,et al.  Pyramidal Cell-Interneuron Circuit Architecture and Dynamics in Hippocampal Networks , 2017, Neuron.

[27]  K. Stratford,et al.  Synaptic transmission between individual pyramidal neurons of the rat visual cortex in vitro , 1991, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[28]  Frederico A. C. Azevedo,et al.  Equal numbers of neuronal and nonneuronal cells make the human brain an isometrically scaled‐up primate brain , 2009, The Journal of comparative neurology.

[29]  Eran Stark,et al.  Unbiased estimation of precise temporal correlations between spike trains , 2009, Journal of Neuroscience Methods.

[30]  Alain Destexhe,et al.  On the Use of Analytical Expressions for the Voltage Distribution to Analyze Intracellular Recordings , 2006, Neural Computation.

[31]  K. Deisseroth,et al.  Millisecond-timescale, genetically targeted optical control of neural activity , 2005, Nature Neuroscience.

[32]  Elizabeth A Stuart,et al.  Matching methods for causal inference: A review and a look forward. , 2010, Statistical science : a review journal of the Institute of Mathematical Statistics.

[33]  Bernhard Schölkopf,et al.  Inferring deterministic causal relations , 2010, UAI.

[34]  Richard A. Nielsen,et al.  Why Propensity Scores Should Not Be Used for Matching , 2019, Political Analysis.

[35]  Gaël Varoquaux,et al.  Scikit-learn: Machine Learning in Python , 2011, J. Mach. Learn. Res..

[36]  Aapo Hyvärinen,et al.  A Linear Non-Gaussian Acyclic Model for Causal Discovery , 2006, J. Mach. Learn. Res..

[37]  Rafael Yuste,et al.  Two-photon photostimulation and imaging of neural circuits , 2007, Nature Methods.

[38]  Philip G. Wright,et al.  The tariff on animal and vegetable oils , 1928 .

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

[40]  Jonathan V. Gill,et al.  Two Photon Holographic Stimulation System for Cellular-Resolution Interrogation of Olfactory Coding , 2017 .

[41]  Stergios B. Fotopoulos,et al.  All of Nonparametric Statistics , 2007, Technometrics.

[42]  M. Häusser,et al.  All-Optical Interrogation of Neural Circuits , 2015, The Journal of Neuroscience.

[43]  Matteo Carandini,et al.  Kilosort: realtime spike-sorting for extracellular electrophysiology with hundreds of channels , 2016, bioRxiv.

[44]  Joshua D. Angrist,et al.  Mostly Harmless Econometrics: An Empiricist's Companion , 2008 .

[45]  M. Abeles Quantification, smoothing, and confidence limits for single-units' histograms , 1982, Journal of Neuroscience Methods.

[46]  David Pfau,et al.  Robust learning of low-dimensional dynamics from large neural ensembles , 2013, NIPS.

[47]  Donghyun Kim,et al.  Handbook of Photonics for Biomedical Engineering , 2016 .

[48]  M. Drton,et al.  Global identifiability of linear structural equation models , 2010, 1003.1146.

[49]  O Sporns,et al.  Predicting human resting-state functional connectivity from structural connectivity , 2009, Proceedings of the National Academy of Sciences.

[50]  A. Destexhe,et al.  Impact of network activity on the integrative properties of neocortical pyramidal neurons in vivo. , 1999, Journal of neurophysiology.

[51]  David Card,et al.  Using Geographic Variation in College Proximity to Estimate the Return to Schooling , 1993 .