Integrative analysis of in vivo recording with single-cell RNA-seq data reveals molecular properties of light-sensitive neurons in mouse V1

Abstract

Vision formation is classically based on projections from retinal ganglion cells (RGC) to the lateral geniculate nucleus (LGN) and the primary visual cortex (V1). Neurons in the mouse V1 are tuned to light stimuli. Although the cellular information of the retina and the LGN has been widely studied, the transcriptome profiles of single light-stimulated neuron in V1 remain unknown. In our study, in vivo calcium imaging and whole-cell electrophysiological patch-clamp recording were utilized to identify 53 individual cells from layer 2/3 of V1 as lightsensitive (LS) or non-light-sensitive (NS) by single-cell light-evoked calcium evaluation and action potential spiking. The contents of each cell after functional tests were aspirated in vivo through a patch-clamp pipette for mRNA sequencing. Moreover, the three-dimensional (3-D) morphological characterizations of the neurons were reconstructed in a live mouse after the whole-cell recordings. Our sequencing results indicated that V1 neurons with a high expression of genes related to transmission regulation, such as Rtn4r and Rgs7, and genes involved in membrane transport, such as Na⁺/K⁺ ATPase and NMDA-type glutamatergic receptors, preferentially responded to light stimulation. Furthermore, an antagonist that blocks Rtn4r signals could inactivate the neuronal responses to light stimulation in live mice. In conclusion, our findings of the vivo-seq analysis indicate the key role of the strength of synaptic transmission possesses neurons in V1 of light sensory.

Keywords

light sensitivity / vivo-seq / patch-seq / calcium imaging in vivo / whole cell recording in vivo

References

[1] Aperia A, Akkuratov EE, Fontana JM, Brismar H (2016) Na+-K+- ATPase, a new class of plasma membrane receptors. Am J Physiol Cell Physiol 310:C491–495

[2] Arganda S, Guantes R, de Polavieja GG (2007) Sodium pumps adapt spike bursting to stimulus statistics. Nat Neurosci 10:1467–1473

[3] Bardy C, van den Hurk M, Kakaradov B, Erwin JA, Jaeger BN, Hernandez RV, Eames T, Paucar AA, Gorris M, Marchand C (2016) Predicting the functional states of human iPSC-derived neurons with single-cell RNA-seq and electrophysiology. Mol Psychiatry 21:1573–1588

[4] Bjartmar L, Huberman AD, Ullian EM, Renteria RC, Liu X, Xu W, Prezioso J, Susman MW, Stellwagen D, Stokes CC (2006) Neuronal pentraxins mediate synaptic refinement in the devel- oping visual system. J Neurosci 26:6269–6281

[5] Bock DD, Lee WC, Kerlin AM, Andermann ML, Hood G, Wetzel AW, Yurgenson S, Soucy ER, Kim HS, Reid RC (2011) Network anatomy and in vivo physiology of visual cortical neurons. Nature 471:177–182

[6] Butler A, Hoffman P, Smibert P, Papalexi E, Satija R (2018) Integrating single-cell transcriptomic data across different condi- tions, technologies, and species. Nat Biotechnol 36:411

[7] Cadwell CR, Palasantza A, Jiang X, Berens P, Deng Q, Yilmaz M, Reimer J, Shen S, Bethge M, Tolias KF (2016) Electro- physiological, transcriptomic and morphologic profiling of single neurons using Patch-seq. Nat Biotechnol 34:199–203

[8] Cadwell CR, Scala F, Li S, Livrizzi G, Shen S, Sandberg R, Jiang X, Tolias AS (2017) Multimodal profiling of single-cell morphology, electrophysiology, and gene expression using Patch-seq. Nat Protoc 12:2531–2553

[9] Cang J, Renteria RC, Kaneko M, Liu X, Copenhagen DR, Stryker MP (2005) Development of precise maps in visual cortex requires patterned spontaneous activity in the retina. Neuron 48:797–809

[10] Chalupa LMW, Werner JS (2003) The visual neurosciences. MIT Press, Cambridge

[11] Chen S, Huang T, Zhou Y, Han Y, Xu M, Gu J (2017) AfterQC: automatic filtering, trimming, error removing and quality control for fastq data. BMC Bioinform 18:80

[12] Chen X, Zhang K, Zhou L, Gao X, Wang J, Yao Y, He F, Luo Y, Yu Y, Li S (2016) Coupled electrophysiological recording and single cell transcriptome analyses revealed molecular mecha- nisms underlying neuronal maturation. Protein Cell 7:175–186

[13] Cruz-Martin A, El-Danaf RN, Osakada F, Sriram B, Dhande OS, Nguyen PL, Callaway EM, Ghosh A, Huberman AD (2014) A dedicated circuit links direction-selective retinal ganglion cells to the primary visual cortex. Nature 507:358–361

[14] Cubelos B, Briz CG, Esteban-Ortega GM, Nieto M (2015) Cux1 and Cux2 selectively target basal and apical dendritic compartments of layer II-III cortical neurons. Dev Neurobiol 75:163–172

[15] da Huang W, Sherman BT, Lempicki RA (2009a) Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 37:1–13

[16] da Huang W, Sherman BT, Lempicki RA (2009b) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4:44–57

[17] Ding Y, Zheng Y, Liu T, Chen T, Wang C, Sun Q, Hua M, Hua T (2017) Changes in GABAergic markers accompany degradation of neuronal function in the primary visual cortex of senescent rats. Sci Rep 7:14897

[18] Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR (2013) STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29:15–21

[19] Fan X, Dong J, Zhong S, Wei Y, Wu Q, Yan L, Yong J, Sun L, Wang X, Zhao Y (2018) Spatial transcriptomic survey of human embryonic cerebral cortex by single-cell RNA-seq analysis. Cell Res 28:730–745

[20] Flavell SW, Greenberg ME (2008) Signaling mechanisms linking neuronal activity to gene expression and plasticity of the nervous system. Annu Rev Neurosci 31:563–590

[21] Fuzik J, Zeisel A, Mate Z, Calvigioni D, Yanagawa Y, Szabo G, Linnarsson S, Harkany T (2016) Integration of electrophysiolog- ical recordings with single-cell RNA-seq data identifies neuronal subtypes. Nat Biotechnol 34:175–183

[22] Gao G, Fernandez CS, Stapleton D, Auster AS, Widmer J, Dyck JR, Kemp BE, Witters LA (1996) Non-catalytic beta- and gamma- subunit isoforms of the 5’-AMP-activated protein kinase. J Biol Chem 271:8675–8681

[23] Gerber KJ, Squires KE, Hepler JR (2016) Roles for Regulator of G Protein Signaling Proteins in Synaptic Signaling and Plasticity. Mol Pharmacol 89:273–286

[24] Gray LT, Yao Z (2017) Layer-specific chromatin accessibility landscapes reveal regulatory networks in adult mouse visual cortex. Elife 6:e21883

[25] Hagihara KM, Ohki K (2013) Long-term down-regulation of GABA decreases orientation selectivity without affecting direction selec- tivity in mouse primary visual cortex. Front Neural Circuits 7:28

[26] Herrero JL, Gieselmann MA, Thiele A (2017) Muscarinic and Nicotinic Contribution to Contrast Sensitivity of Macaque Area V1 Neurons. Front Neural Circuits 11:106

[27] Hille B (2001) Ion channels of excitable membranes, 3rd edn. Sinauer, Sunderland

[28] Hrvatin S, Hochbaum DR, Nagy MA, Cicconet M, Robertson K, Cheadle L, Zilionis R, Ratner A, Borges-Monroy R, Klein AM (2018) Single-cell analysis of experience-dependent transcrip- tomic states in the mouse visual cortex. Nat Neurosci 21:120–129

[29] Hucho F (1993) Neurotransmitter receptors. Elsevier, Amsterdam

[30] Juliandi B, Abematsu M, Sanosaka T, Tsujimura K, Smith A, Nakashima K (2012) Induction of superficial cortical layer neurons from mouse embryonic stem cells by valproic acid. Neurosci Res 72:23–31

[31] Kaczmarek L, Chaudhuri A (1997) Sensory regulation of immediate– early gene expression in mammalian visual cortex: implications for functional mapping and neural plasticity. Brain Res Brain Res Rev 23:237–256

[32] Kersey PJ, Allen JE, Allot A, Barba M, Boddu S, Bolt BJ, Carvalho- Silva D, Christensen M, Davis P, Grabmueller C (2018) Ensembl Genomes 2018: an integrated omics infrastructure for non-vertebrate species. Nucleic Acids Res 46:D802–D808

[33] Kirkwood A, Rioult MC, Bear MF (1996) Experience-dependent modification of synaptic plasticity in visual cortex. Nature 381:526–528

[34] Lam SK, Yoda N, Schekman R (2010) A vesicle carrier that mediates peroxisome protein traffic from the endoplasmic reticulum. Proc Natl Acad Sci USA 107:21523–21528

[35] Langfelder P, Horvath S (2008) WGCNA: an R package for weighted correlation network analysis. BMC Bioinform 9:559

[36] Langfelder P, Horvath S (2012) Fast R functions for robust correlations and hierarchical clustering. J Stat Softw 46

[37] Lee SH, Kwan AC, Zhang S, Phoumthipphavong V, Flannery JG, Masmanidis SC, Taniguchi H, Huang ZJ, Zhang F, Boyden ES (2012) Activation of specific interneurons improves V1 feature selectivity and visual perception. Nature 488:379–383

[38] Lee SJ, Wei M, Zhang C, Maxeiner S, Pak C, Calado Botelho S, Trotter J, Sterky FH, Sudhof TC (2017) Presynaptic neuronal pentraxin receptor organizes excitatory and inhibitory synapses. J Neurosci 37:1062–1080

[39] Lee WC, Bonin V, Reed M, Graham BJ, Hood G, Glattfelder K, Reid RC (2016) Anatomy and function of an excitatory network in the visual cortex. Nature 532:370–374

[40] Leifer D, Krainc D, Yu YT, McDermott J, Breitbart RE, Heng J, Neve RL, Kosofsky B, Nadal-Ginard B, Lipton SA (1993) MEF2C, a MADS/MEF2-family transcription factor expressed in a laminar distribution in cerebral cortex. Proc Natl Acad Sci USA 90:1546–1550

[41] Li CL, Li KC, Wu D, Chen Y, Luo H, Zhao JR, Wang SS, Sun MM, Lu YJ, Zhong YQ (2016) Somatosensory neuron types identified by high-coverage single-cell RNA-sequencing and functional heterogeneity. Cell Res 26:967

[42] Li J, Zhang J, Wang M, Pan J, Chen X, Liao X (2017a) Functional imaging of neuronal activity of auditory cortex by using Cal-520 in anesthetized and awake mice. Biomed Opt Express 8:2599–2610

[43] Li S, Wang L, Tie X, Sohya K, Lin X, Kirkwood A, Jiang B (2017b) Brief novel visual experience fundamentally changes synaptic plasticity in the mouse visual cortex. J Neurosci 37:9353–9360

[44] Liao Y, Smyth GK, Shi W (2014) featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30:923–930

[45] Liu J, Liu W, Yang L, Wu Q, Zhang H, Fang A, Li L, Xu X, Sun L, Zhang J (2017) The primate-specific gene TMEM14B marks outer radial glia cells and promotes cortical expansion and folding. Cell Stem Cell 21(635–649):e638

[46] Love MI, Huber W, Anders S (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15:550

[47] Madabhushi R, Gao F, Pfenning AR, Pan L, Yamakawa S, Seo J, Rueda R, Phan T, Yamakawa H, Pao PC (2015) Activity- induced DNA breaks govern the expression of neuronal early- response genes. Cell 161:1592–1605

[48] Malagon G, Miki T, Llano I, Neher E, Marty A (2016) Counting vesicular release events reveals binomial release statistics at single glutamatergic synapses. J Neurosci 36:4010–4025

[49] Mangini NJ, Pearlman AL (1980) Laminar distribution of receptive field properties in the primary visual cortex of the mouse. J Comp Neurol 193:203–222

[50] McCarthy DJ, Campbell KR, Lun AT, Wills QF (2017) Scater: pre- processing, quality control, normalization and visualization of single-cell RNA-seq data in R. Bioinformatics 33:1179–1186

[51] Meisami E, Timiras PS (1974) Influence of early visual deprivation on regional activity of brain ATPases in developing rats. J Neu- rochem 22:725–729

[52] Moroni RF, Inverardi F, Regondi MC, Watakabe A, Yamamori T, Spreafico R, Frassoni C (2009) Expression of layer-specific markers in the adult neocortex of BCNU-Treated rat, a model of cortical dysplasia. Neuroscience 159:682–691

[53] Niell CM, Stryker MP (2008) Highly selective receptive fields in mouse visual cortex. J Neurosci 28:7520–7536

[54] O’Brien RJ, Kamboj S, Ehlers MD, Rosen KR, Fischbach GD, Huganir RL (1998) Activity-dependent modulation of synaptic AMPA receptor accumulation. Neuron 21:1067–1078

[55] Okuno H (2011) Regulation and function of immediate-early genes in the brain: beyond neuronal activity markers. Neurosci Res 69:175–186

[56] Paoletti P, Bellone C, Zhou Q (2013) NMDA receptor subunit diversity: impact on receptor properties, synaptic plasticity and disease. Nat Rev Neurosci 14:383–400

[57] Pei X, Vidyasagar TR, Volgushev M, Creutzfeldt OD (1994) Receptive field analysis and orientation selectivity of postsynaptic potentials of simple cells in cat visual cortex. J Neurosci 14:7130–7140

[58] Pelkey KA, Barksdale E, Craig MT, Yuan X, Sukumaran M, Vargish GA, Mitchell RM, Wyeth MS, Petralia RS, Chittajallu R (2015) Pentraxins coordinate excitatory synapse maturation and circuit integration of parvalbumin interneurons. Neuron 85:1257–1272

[59] Picelli S, Faridani OR, Bjorklund AK, Winberg G, Sagasser S, Sandberg R (2014) Full-length RNA-seq from single cells using Smart-seq2. Nat Protoc 9:171–181

[60] Plossl K, Royer M, Bernklau S, Tavraz NN, Friedrich T, Wild J, Weber BHF, Friedrich U (2017) Retinoschisin is linked to retinal Na/K-ATPase signaling and localization. Mol Biol Cell 28:2178–2189

[61] Sarria I, Orlandi C, McCall MA, Gregg RG, Martemyanov KA (2016) Intermolecular interaction 2between anchoring subunits specify subcellular targeting and function of RGS proteins in retina ON- bipolar neurons. J Neurosci 36:2915–2925

[62] Sato H, Hata Y, Masui H, Tsumoto T (1987) A functional role of cholinergic innervation to neurons in the cat visual cortex. J Neurophysiol 58:765–780

[63] Sia GM, Beique JC, Rumbaugh G, Cho R, Worley PF, Huganir RL (2007) Interaction of the N-terminal domain of the AMPA receptor GluR4 subunit with the neuronal pentraxin NP1 mediates GluR4 synaptic recruitment. Neuron 55:87–102

[64] Sillito AM, Kemp JA (1983) Cholinergic modulation of the functional organization of the cat visual cortex. Brain Res 289:143–155

[65] Stephany CE, Chan LL, Parivash SN, Dorton HM, Piechowicz M, Qiu S, McGee AW (2014) Plasticity of binocularity and visual acuity are differentially limited by nogo receptor. J Neurosci 34:11631–11640

[66] Stephany CE, Frantz MG, McGee AW (2016a) Multiple roles for nogo receptor 1 in visual system plasticity. Neuroscientist 22:653–666

[67] Stephany CE, Ikrar T, Nguyen C, Xu X, McGee AW (2016b) Nogo receptor 1 confines a disinhibitory microcircuit to the critical period in visual cortex. J Neurosci 36:11006–11012

[68] Stosiek C, Garaschuk O, Holthoff K, Konnerth A (2003) In vivo two- photon calcium imaging of neuronal networks. Proc Natl Acad Sci USA 100:7319–7324

[69] Thomas RA, Gibon J, Chen CXQ, Chierzi S, Soubannier VG, Baulac S, Seguela P, Murai K, Barker PA (2018) The nogo receptor ligand LGI1 regulates synapse number and synaptic activity in hippocampal and cortical neurons. eNeuro. https://doi.org/10.1523/ENEURO.0185-18.2018

[70] Traynelis SF, Wollmuth LP, McBain CJ, Menniti FS, Vance KM, Ogden KK, Hansen KB, Yuan H, Myers SJ, Dingledine R (2010) Glutamate receptor ion channels: structure, regulation, and function. Pharmacol Rev 62:405–496

[71] Wu QF, Yang L, Li S, Wang Q, Yuan XB, Gao X, Bao L, Zhang X (2012)Fibroblast growth factor 13 is a microtubule-stabilizing protein regulating neuronal polarization and migration. Cell 149:1549–1564

[72] Zariwala HA, Madisen L, Ahrens KF, Bernard A, Lein ES, Jones AR, Zeng H (2011) Visual tuning properties of genetically identified layer 2/3 neuronal types in the primary visual cortex of cre- transgenic mice. Front Syst Neurosci 4:162

[73] Zhong S, Zhang S, Fan X, Wu Q, Yan L, Dong J, Zhang H, Li L, Sun L, Pan N (2018) A single-cell RNA-seq survey of the developmental landscape of the human prefrontal cortex. Nature 555:524

RIGHTS & PERMISSIONS

2020 The Author(s)