Abstract:Central respiratory chemoreceptors are specialized neurons with intrinsic sensitivity to hypercapnia and/or acidosis that couple breathing and pH/CO2 levels. Prior data indicate that a sub-population of brainstem serotonin (5-HT) neurons are likely central respiratory chemoreceptors. However, it remains unclear which 5-HT neurons develop chemosensitivity and what molecular markers may identify this unique sub-population of 5-HT neurons. Here we employed the “patch-to-seq” technique to measure action potential firing rate responses of 5-HT neuron to hypercapnic acidosis using cell-attached patch-clamp electrophysiology followed by isolation of that neurons’ intracellular RNA content for single cell RNA Sequencing. Acute brainstem slices (200 µm) from young (P18-23) transgenic rats expressing eGFP in all 5-HT neurons (SSeGFP) were recorded while superfused with artificial CSF (aCSF) containing inhibitors for synaptic blockade (10 mM CNQX, 50 mM D-AP5, 20 mM Gabazine) bubbled with either 5% CO2 (bal. O2; pH= 7.36; 5 min) or 15% CO2 (bal. O2; pH= 7.10). The Chemosensitivity Index (C.I.) determined cell phenotypes. pH-sensitive eGFP+ neurons (n=48) had an average C.I. of 169.8± 8.25 SEM whereas pH-insensitive eGFP+ neurons (n=44) had an average C.I. of 100.3 ± 1.4 (P < 0.0001; t-test). Cell extracts from 11 chemosensitive and 10 insensitive recorded cells were subjected to single cell RNA Sequencing (scRNA-Seq), from which an average of 48.3M reads were generated with a quality score of ~32.7 and ~75% mapping rate per sample with an average of ~8,000 genes detected. Purity of the isolated samples was confirmed by significantly more expression (Log2(FPKM+1)) of 6 serotonergic vs 7 neuronal and 8 glial gene markers, and more expression of 7 neuronal vs 8 known glial gene markers. There were 166differentially expressed genes (q < 0.05) between pH-sensitive and -insensitive 5-HT neuron populations none of which have a known contribution to pH regulation or sensing. Predicted upstream regulators of 352 DEGs (p<0.01) point to the importance of Egr2 activation, a developmental transcription factor previously shown to give rise to pH sensitive 5-HT neurons in mice. Machine learning algorithms identified unique markers of CO2/pH sensitive 5-HT neurons in rats, where supervised analyses (using Support Vector Machine learning and Recursive Feature Elimination (SVM-RFE)) yielded several genes in which expression levels significantly (P<0.001) and accurately (%) predict the predefined phenotype from the DEGs (CD46, 68.7%; Iba57, 71.2%; Chd5, 72.0%; and Maz, 74.0%) and among the 1,000s of detected genes (Hnmpk, 81.0%; Diablo, 87.9%) but the biological/function significance of these genes in CO2/pH sensing remains unclear. RNAScope confirmed 25% of Tph2+ neurons in the raphe midline co-express 4 or more RNA molecules of CD46 or Iba57 RNA. Together, these data suggest: 1) 5-HT pH sensitive neurons are transcriptionally distinct from insensitive 5-HT neurons, and 2) novel molecular markers of 5-HT pH sensitive rat neurons may be useful identifiers for this unique subpopulation.

5-HT3 Receptors in Rat Dorsal Root Ganglion Neurons: Ca2+ Entry and Modulation of Neurotransmitter Release

P2X7 Receptors in Satellite Glial Cells Mediate High Functional Expression of P2X3 Receptors in Immature Dorsal Root Ganglion Neurons

Patch‐to‐Seq Reveals Unique Transcriptomic Profiles of Chemosensitive Serotonergic Neurons

Developmental Expression Of Ph-Sensitive Kir4.1 And Kir5.1 Channels In Serotonergic (5-Ht) Neurons

5-Ht2c Receptors Modulate the Discharge Activities of Inspiratory Neurons in the Medial Region of Nucleus Retrofacialis of Neonatal Rats in Vitro

Identifying Candidate Genes that Underlie Cellular pH Sensitivity in Serotonin Neurons Using Transcriptomics: A Potential Role for Kir5.1 Channels.

5-HT potentiation of the GABA(A) response in the rat sacral dorsal commissural neurones

Comparison of metabotropic glutamate receptor responses at segmental and descending inputs to motoneurons in neonatal rat spinal cord

Electrophysiological Properties of Rat Spinal Dorsal Horn Neurones in Vitro: Calcium‐dependent Action Potentials.

5-HT3A receptor subunit is expressed in a subpopulation of GABAergic and enkephalinergic neurons in the mouse dorsal spinal cord

G Protein-Coupled Inwardly Rectifying Potassium Channels in Dorsal Root Ganglion Neurons

Functional properties of ionotropic glutamate receptor channels in rat sacral dorsal commissural neurons

Potentiation of 5-HT3 Receptor Function by the Activation of Coexistent 5-HT2 Receptors in Trigeminal Ganglion Neurons of Rats

Propagation of neuronal excitation in the rat central nucleus of amygdala revealed by the optical imaging

Activation of 5-HT3 receptors in the rat and mouse intestinal tract: a comparative study

Expression Pattern of Type 3 Adenylyl Cyclase in Rodent Dorsal Root Ganglion and Its Primary Afferent Terminals

Small Sensory Neurons in the Rat Dorsal Root Ganglia Express Functional NK‐1 Tachykinin Receptor

Serotonin Enhances Depolarizing Spontaneous Fluctuations, Excitability, and Ongoing Activity in Isolated Rat DRG Neurons Via 5-HT4 Receptors and Camp-Dependent Mechanisms

Histamine H3 Heteroreceptors Suppress Glutamatergic and GABAergic Synaptic Transmission in the Rat Insular Cortex

GABAergic Neurons in the Dorsal Raphe Nucleus that Express 5-HT3A Receptors Participate in Responses to Stress Hormones

Satellite glia cells in dorsal root ganglia express functional NMDA receptors