Single-neuron models linking electrophysiology, morphology, and transcriptomics across cortical cell types

Identifying the cell types constituting brain circuits is a fundamental question in neuroscience and motivates the generation of taxonomies based on electrophysiological, morphological and molecular single cell properties. Establishing the correspondence across data modalities and understanding the underlying principles has proven challenging. Bio-realistic computational models offer the ability to probe cause-and-effect and have historically been used to explore phenomena at the single-neuron level. Here we introduce a computational optimization workflow used for the generation and evaluation of more than 130 million single neuron models with active conductances. These models were based on 230 in vitro electrophysiological experiments followed by morphological reconstruction from the mouse visual cortex. We show that distinct ion channel conductance vectors exist that distinguish between major cortical classes with passive and h-channel conductances emerging as particularly important for classification. Next, using models of genetically defined classes, we show that differences in specific conductances predicted from the models reflect differences in gene expression in excitatory and inhibitory cell types as experimentally validated by single-cell RNA-sequencing. The differences in these conductances, in turn, explain many of the electrophysiological differences observed between cell types. Finally, we show the robustness of the herein generated single-cell models as representations and realizations of specific cell types in face of biological variability and optimization complexity. Our computational effort generated models that reconcile major single-cell data modalities that define cell types allowing for causal relationships to be examined. Highlights Generation and evaluation of more than 130 million single-cell models with active conductances along the reconstructed morphology faithfully recapitulate the electrophysiology of 230 in vitro experiments. Optimized ion channel conductances along the cellular morphology (‘all-active’) are characteristic of model complexity and offer enhanced biophysical realism. Ion channel conductance vectors of all-active models classify transcriptomically defined cell-types. Cell type differences in ion channel conductances predicted by the models correlate with experimentally measured single-cell gene expression differences in inhibitory (Pvalb, Sst, Htr3a) and excitatory (Nr5a1, Rbp4) classes. A set of ion channel conductances identified by comparing between cell type model populations explain electrophysiology differences between these types in simulations and brain slice experiments. All-active models recapitulate multimodal properties of excitatory and inhibitory cell types offering a systematic and causal way of linking differences between them.

[1]  B. Hille Ionic channels in excitable membranes. Current problems and biophysical approaches. , 1978, Biophysical journal.

[2]  Patrick R. Hof,et al.  Automated evolutionary optimization of ion channel conductances and kinetics in models of young and aged rhesus monkey pyramidal neurons , 2016, Journal of Computational Neuroscience.

[3]  H. Markram,et al.  Organizing principles for a diversity of GABAergic interneurons and synapses in the neocortex. , 2000, Science.

[4]  T. Sejnowski,et al.  A model of spike initiation in neocortical pyramidal neurons , 1995, Neuron.

[5]  A. Erisir,et al.  Function of specific K(+) channels in sustained high-frequency firing of fast-spiking neocortical interneurons. , 1999, Journal of neurophysiology.

[6]  T. Hackett,et al.  A Distributed Network for Social Cognition Enriched for Oxytocin Receptors , 2016, The Journal of Neuroscience.

[7]  J. Christie,et al.  Distinct Kv Channel Subtypes Contribute to Differences in Spike Signaling Properties in the Axon Initial Segment and Presynaptic Boutons of Cerebellar Interneurons , 2014, The Journal of Neuroscience.

[8]  Henry Markram,et al.  BluePyOpt: Leveraging Open Source Software and Cloud Infrastructure to Optimise Model Parameters in Neuroscience , 2016, Front. Neuroinform..

[9]  G. Edelman,et al.  Degeneracy and complexity in biological systems , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[10]  P. Jonas,et al.  Kv3 Potassium Conductance is Necessary and Kinetically Optimized for High-Frequency Action Potential Generation in Hippocampal Interneurons , 2003, The Journal of Neuroscience.

[11]  Henry Markram,et al.  A Novel Multiple Objective Optimization Framework for Constraining Conductance-Based Neuron Models by Experimental Data , 2007, Front. Neurosci..

[12]  S. Linnarsson,et al.  Cell types in the mouse cortex and hippocampus revealed by single-cell RNA-seq , 2015, Science.

[13]  G A Ascoli,et al.  Generation, description and storage of dendritic morphology data. , 2001, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[14]  M. Daly,et al.  A recurrent de novo mutation in KCNC1 causes progressive myoclonus epilepsy , 2014, Nature Genetics.

[15]  D. McCormick,et al.  H-Current Properties of a Neuronal and Network Pacemaker , 1998, Neuron.

[16]  Brian R. Lee,et al.  Classification of electrophysiological and morphological neuron types in the mouse visual cortex , 2019, Nature Neuroscience.

[17]  Henry Markram,et al.  Models of Neocortical Layer 5b Pyramidal Cells Capturing a Wide Range of Dendritic and Perisomatic Active Properties , 2011, PLoS Comput. Biol..

[18]  Allan R. Jones,et al.  Conserved cell types with divergent features in human versus mouse cortex , 2019, Nature.

[19]  Walther Akemann,et al.  Transgenic mice expressing a fluorescent in vivo label in a distinct subpopulation of neocortical layer 5 pyramidal cells , 2004, The Journal of comparative neurology.

[20]  David S. Greenberg,et al.  Training deep neural density estimators to identify mechanistic models of neural dynamics , 2019, bioRxiv.

[21]  E. Welker,et al.  K+ Channel Expression Distinguishes Subpopulations of Parvalbumin- and Somatostatin-Containing Neocortical Interneurons , 1999, The Journal of Neuroscience.

[22]  E. Marder,et al.  Similar network activity from disparate circuit parameters , 2004, Nature Neuroscience.

[23]  Henry Markram,et al.  Channelpedia: An Integrative and Interactive Database for Ion Channels , 2011, Front. Neuroinform..

[24]  Hyungwon Choi,et al.  Moving beyond P values: data analysis with estimation graphics , 2019, Nature Methods.

[25]  Nicholas T. Carnevale,et al.  The NEURON Simulation Environment , 1997, Neural Computation.

[26]  Arnd Roth,et al.  Automated optimization of a reduced layer 5 pyramidal cell model based on experimental data , 2012, Journal of Neuroscience Methods.

[27]  T. Perney,et al.  Allele‐dependent changes of olivocerebellar circuit properties in the absence of the voltage‐gated potassium channels Kv3.1 and Kv3.3 , 2004, The European journal of neuroscience.

[28]  C. McBain,et al.  Frequency‐dependent regulation of rat hippocampal somato‐dendritic excitability by the K+ channel subunit Kv2.1 , 2000, The Journal of physiology.

[29]  Trygve E Bakken,et al.  h-Channels Contribute to Divergent Intrinsic Membrane Properties of Supragranular Pyramidal Neurons in Human versus Mouse Cerebral Cortex , 2018, Neuron.

[30]  H. Markram,et al.  Interneurons of the neocortical inhibitory system , 2004, Nature Reviews Neuroscience.

[31]  Christof Koch,et al.  Neurodata Without Borders: Creating a Common Data Format for Neurophysiology , 2015, Neuron.

[32]  Staci A. Sorensen,et al.  Adult Mouse Cortical Cell Taxonomy Revealed by Single Cell Transcriptomics , 2016 .

[33]  Thomas K. Berger,et al.  Effective Stimuli for Constructing Reliable Neuron Models , 2011, PLoS Comput. Biol..

[34]  Christopher. Simons,et al.  Machine learning with Python , 2017 .

[35]  L. Jan,et al.  The distribution and targeting of neuronal voltage-gated ion channels , 2006, Nature Reviews Neuroscience.

[36]  Yuchio Yanagawa,et al.  Integration of electrophysiological recordings with single-cell RNA-seq data identifies novel neuronal subtypes , 2015, Nature Biotechnology.

[37]  C. Pedroarena BK and Kv3.1 Potassium Channels Control Different Aspects of Deep Cerebellar Nuclear Neurons Action Potentials and Spiking Activity , 2011, The Cerebellum.

[38]  Bernardo Rudy,et al.  Kv3 channels: voltage-gated K+ channels designed for high-frequency repetitive firing , 2001, Trends in Neurosciences.

[39]  Tjerk P. Straatsma,et al.  NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations , 2010, Comput. Phys. Commun..

[40]  B. Kampa,et al.  Action potential generation requires a high sodium channel density in the axon initial segment , 2008, Nature Neuroscience.

[41]  Marc Parizeau,et al.  DEAP: evolutionary algorithms made easy , 2012, J. Mach. Learn. Res..

[42]  H. Markram,et al.  Correlation maps allow neuronal electrical properties to be predicted from single-cell gene expression profiles in rat neocortex. , 2004, Cerebral cortex.

[43]  Hanchuan Peng,et al.  Extensible visualization and analysis for multidimensional images using Vaa3D , 2014, Nature Protocols.

[44]  Michael L. Hines,et al.  Neuroinformatics Original Research Article Neuron and Python , 2022 .

[45]  Eve Marder,et al.  Alternative to hand-tuning conductance-based models: construction and analysis of databases of model neurons. , 2003, Journal of neurophysiology.

[46]  Hongkui Zeng,et al.  Neuronal cell-type classification: challenges, opportunities and the path forward , 2017, Nature Reviews Neuroscience.

[47]  Properties of a Neuronal and Network Pacemaker , 1998 .

[48]  Kristofer E. Bouchard,et al.  Inferring neuronal ionic conductances from membrane potentials using CNNs , 2019, bioRxiv.

[49]  Christof Koch,et al.  Systematic generation of biophysically detailed models for diverse cortical neuron types , 2018, Nature Communications.

[50]  Christof Koch,et al.  Generalized leaky integrate-and-fire models classify multiple neuron types , 2017, Nature Communications.

[51]  Jens Hjerling-Leffler,et al.  Disentangling neural cell diversity using single-cell transcriptomics , 2016, Nature Neuroscience.

[52]  Leonard Maler,et al.  Intrinsic frequency tuning in ELL pyramidal cells varies across electrosensory maps. , 2008, Journal of neurophysiology.

[53]  X. Chen,et al.  Genomic organization, nucleotide sequence, biophysical properties, and localization of the voltage-gated K+ channel gene KCNA4/Kv1.4 to mouse chromosome 2/human 11p14 and mapping of KCNC1/Kv3.1 to mouse 7/human 11p14.3-p15.2 and KCNA1/Kv1.1 to human 12p13. , 1994, Genomics.

[54]  John R Huguenard,et al.  Resilient RTN fast spiking in Kv3.1 null mice suggests redundancy in the action potential repolarization mechanism. , 2002, Journal of neurophysiology.

[55]  Catalin C. Barbacioru,et al.  mRNA-Seq whole-transcriptome analysis of a single cell , 2009, Nature Methods.

[56]  Paul Pavlidis,et al.  Transcriptomic correlates of neuron electrophysiological diversity , 2017, bioRxiv.

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

[58]  Guy Eyal,et al.  Human Cortical Pyramidal Neurons: From Spines to Spikes via Models , 2018, bioRxiv.

[59]  D. Johnston,et al.  Active properties of neuronal dendrites. , 1996, Annual review of neuroscience.

[60]  Michael Z. Lin,et al.  A Suite of Transgenic Driver and Reporter Mouse Lines with Enhanced Brain-Cell-Type Targeting and Functionality , 2018, Cell.

[61]  Kaisa Miettinen,et al.  Nonlinear multiobjective optimization , 1998, International series in operations research and management science.

[62]  Gaute T. Einevoll,et al.  Uncertainpy: A Python Toolbox for Uncertainty Quantification and Sensitivity Analysis in Computational Neuroscience , 2018, bioRxiv.

[63]  Henry Markram,et al.  A Kinetic Map of the Homomeric Voltage-Gated Potassium Channel (Kv) Family , 2019, Front. Cell. Neurosci..

[64]  Iain Murray,et al.  Sequential Neural Likelihood: Fast Likelihood-free Inference with Autoregressive Flows , 2018, AISTATS.

[65]  Athanasia G. Palasantza,et al.  Electrophysiological, transcriptomic and morphologic profiling of single neurons using Patch-seq , 2015, Nature Biotechnology.

[66]  Allan R. Jones,et al.  Shared and distinct transcriptomic cell types across neocortical areas , 2018, Nature.

[67]  Shreejoy J. Tripathy,et al.  Transcriptomic correlates of electrophysiological and morphological diversity within and across excitatory and inhibitory neuron classes , 2019, PLoS Comput. Biol..

[68]  Robert Tibshirani,et al.  Estimating the number of clusters in a data set via the gap statistic , 2000 .

[69]  Q. Hogan,et al.  Hyperpolarization-activated current (I h) contributes to excitability of primary sensory neurons in rats , 2008, Brain Research.