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.

Von Economo Neurons and Fork Cells: A Neurochemical Signature Linked to Monoaminergic Function

Single cell transcriptomics of vomeronasal neuroepithelium reveals a differential endoplasmic reticulum environment amongst neuronal subtypes

Transcriptional profiling aligned with in situ expression image analysis reveals mosaically expressed molecular markers for GABA neuron sub‐groups in the ventral tegmental area

Proportion and distribution of neurotransmitter-defined cell types in the ventral tegmental area and substantia nigra pars compacta

Neuron-type specificity of dorsal raphe projections to ventral tegmental area

Serotonin Neurons in the Dorsal and Medial Raphe Nuclei: from Single-Cell Transcriptomes to Whole-Brain Projections

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

Genetic Single Neuron Anatomy reveals fine granularity of cortical interneuron subtypes

Molecular interrogation of hypothalamic organization reveals distinct dopamine neuronal subtypes

The discovery of central monoamine neurons gave volume transmission to the wired brain

Molecular Profiling of VGluT1 AND VGluT2 Ventral Subiculum to Nucleus Accumbens Shell Projections

Genetic Single Neuron Anatomy Reveals Fine Granularity of Cortical Axo-Axonic Cells.

The proteomic profile of Von Economo neurons supports their selective vulnerability in bvFTD

Unique functional responses differentially map onto genetic subtypes of dopamine neurons

Transcriptomic diversity of amygdalar subdivisions across humans and nonhuman primates

nNOS-Expressing Neurons in the Ventral Tegmental Area and Substantia Nigra Pars Compacta

Shades of gray in human white matter

Spatial Distribution of Neurons Expressing Single, Double, and Triple Molecular Characteristics of Glutamatergic, Dopaminergic, or GABAergic Neurons in the Mouse Ventral Tegmental Area

Human neocortical expansion involves glutamatergic neuron diversification

Epigenomic Diversity of Cortical Projection Neurons in the Mouse Brain

Whole-Brain Monosynaptic Afferents to Rostromedial Tegmental Nucleus Gamma-Aminobutyric Acid-Releasing Neurons in Mice