Patch-to-Seq and Transcriptomic Analyses Yield Molecular Markers of Functionally Distinct Brainstem Serotonin Neurons

Acute regulation of CO₂ and pH homeostasis requires sensory feedback from peripheral (carotid body) and central (central) CO₂/pH sensitive cells - so called respiratory chemoreceptors. Subsets of brainstem serotonin (5-HT) neurons in the medullary raphe are CO₂ sensitive or insensitive based on differences in embryonic origin, suggesting these functionally distinct subpopulations may have unique transcriptional profiles. Here, we used Patch-to-Seq to determine if the CO₂ responses in brainstem 5-HT neurons could be correlated to unique transcriptional profiles and/or unique molecular markers and pathways. First, firing rate changes with hypercapnic acidosis were measured in fluorescently labeled 5-HT neurons in acute brainstem slices from transgenic, Dahl SS (SSMcwi) rats expressing T2/ePet-eGFP transgene in Pet-1 expressing (serotonin) neurons (SS ^ePet1-eGFP rats). Subsequently, the transcriptomic and pathway profiles of CO₂ sensitive and insensitive 5-HT neurons were determined and compared by single cell RNA (scRNAseq) and bioinformatic analyses. Low baseline firing rates were a distinguishing feature of CO₂ sensitive 5-HT neurons. scRNAseq of these recorded neurons revealed 166 differentially expressed genes among CO₂ sensitive and insensitive 5-HT neurons. Pathway analyses yielded novel predicted upstream regulators, including the transcription factor Egr2 and Leptin. Additional bioinformatic analyses identified 6 candidate gene markers of CO₂ sensitive 5-HT neurons, and 2 selected candidate genes (CD46 and Iba57) were both expressed in 5-HT neurons determined via in situ mRNA hybridization. Together, these data provide novel insights into the transcriptional control of cellular chemoreception and provide unbiased candidate gene markers of CO₂ sensitive 5-HT neurons. Methodologically, these data highlight the utility of the patch-to-seq technique in enabling the linkage of gene expression to specific functions, like CO₂ chemoreception, in a single cell to identify potential mechanisms underlying functional differences in otherwise similar cell types.

Identification of TPBG-Expressing Amacrine Cells in DAT-tdTomato Mouse

Purpose

Neurons are the bricks of the neuronal system and experimental access to certain neuron subtypes will be of great help to decipher neuronal circuits. Here, we identified trophoblast glycoprotein (TPBG)-expressing GABAergic amacrine cells (ACs) that were selectively labeled in DAT-tdTomato transgenic mice.

Methods

Retina and brain sections were prepared for immunostaining with antibodies against various biomarkers. Patch-sequencing was performed to obtain the transcriptomes of tdTomato-positive cells in DAT-tdTomato mice. Whole-cell recordings were conducted to identify responses to light stimulation.

Results

Tyrosine hydroxylase immunoreactive cells were colocalized with tdTomato-positive cells in substantia nigra pars compacta, but not in the retina. Transcriptomes collected from tdTomato-positive cells in retinas via Patch-sequencing exhibited the expression of marker genes of ACs (Pax6 and Slc32a1) and marker genes of GABAergic neurons (Gad1, Gad2, and Slc6a1). Immunostaining with antibodies against relevant proteins (GAD67, GAD65, and GABA) also confirmed transcriptomic results. Furthermore, tdTomato-positive cells in retinas selectively expressed Tpbg, a marker gene for distinct clusters molecularly defined, which was proved with TPBG immunoreactivity in fluorescently labeled cells. Finally, tdTomato-positive cells recorded showed ON-OFF responses to light stimulation.

Conclusions

Ectopic expression occurs in the retina but not in the substantia nigra pars compacta in the DAT-tdTomato mouse, and fluorescently labeled cells in the retina are TPBG-expressing GABAergic ACs. This type of transgenic mice has been proved as an ideal tool to achieve efficient labeling of a distinct subset of ACs that selectively express Tpbg.

Heterogenous impairment of α cell function in type 2 diabetes is linked to cell maturation state

In diabetes, glucagon secretion from pancreatic α cells is dysregulated. The underlying mechanisms, and whether dysfunction occurs uniformly among cells, remain unclear. We examined α cells from human donors and mice using electrophysiological, transcriptomic, and computational approaches. Rising glucose suppresses α cell exocytosis by reducing P/Q-type Ca2+ channel activity, and this is disrupted in type 2 diabetes (T2D). Upon high-fat feeding of mice, α cells shift toward a "β cell-like" electrophysiological profile in concert with indications of impaired identity. In human α cells we identified links between cell membrane properties and cell surface signaling receptors, mitochondrial respiratory chain complex assembly, and cell maturation. Cell-type classification using machine learning of electrophysiology data demonstrated a heterogenous loss of "electrophysiologic identity" in α cells from donors with type 2 diabetes. Indeed, a subset of α cells with impaired exocytosis is defined by an enrichment in progenitor and lineage markers and upregulation of an immature transcriptomic phenotype, suggesting important links between α cell maturation state and dysfunction.

Signature morpho-electric, transcriptomic, and dendritic properties of human layer 5 neocortical pyramidal neurons

In the neocortex, subcerebral axonal projections originate largely from layer 5 (L5) extratelencephalic-projecting (ET) neurons. The unique morpho-electric properties of these neurons have been mainly described in rodents, where retrograde tracers or transgenic lines can label them. Similar labeling strategies are infeasible in the human neocortex, rendering the translational relevance of findings in rodents unclear. We leveraged the recent discovery of a transcriptomically defined L5 ET neuron type to study the properties of human L5 ET neurons in neocortical brain slices derived from neurosurgeries. Patch-seq recordings, where transcriptome, physiology, and morphology were assayed from the same cell, revealed many conserved morpho-electric properties of human and rodent L5 ET neurons. Divergent properties were often subtler than differences between L5 cell types within these two species. These data suggest a conserved function of L5 ET neurons in the neocortical hierarchy but also highlight phenotypic divergence possibly related to functional specialization of human neocortex.

Scaled, high fidelity electrophysiological, morphological, and transcriptomic cell characterization

The Patch-seq approach is a powerful variation of the patch-clamp technique that allows for the combined electrophysiological, morphological, and transcriptomic characterization of individual neurons. To generate Patch-seq datasets at scale, we identified and refined key factors that contribute to the efficient collection of high-quality data. We developed patch-clamp electrophysiology software with analysis functions specifically designed to automate acquisition with online quality control. We recognized the importance of extracting the nucleus for transcriptomic success and maximizing membrane integrity during nucleus extraction for morphology success. The protocol is generalizable to different species and brain regions, as demonstrated by capturing multimodal data from human and macaque brain slices. The protocol, analysis and acquisition software are compiled at https://githubcom/AllenInstitute/patchseqtools. This resource can be used by individual labs to generate data across diverse mammalian species and that is compatible with large publicly available Patch-seq datasets.

Heterogeneity of neurons reprogrammed from spinal cord astrocytes by the proneural factors Ascl1 and Neurogenin2

Astrocytes are a viable source for generating new neurons via direct conversion. However, little is known about the neurogenic cascades triggered in astrocytes from different regions of the CNS. Here, we examine the transcriptome induced by the proneural factors Ascl1 and Neurog2 in spinal cord-derived astrocytes in vitro. Each factor initially elicits different neurogenic programs that later converge to a V2 interneuron-like state. Intriguingly, patch sequencing (patch-seq) shows no overall correlation between functional properties and the transcriptome of the heterogenous induced neurons, except for K-channels. For example, some neurons with fully mature electrophysiological properties still express astrocyte genes, thus calling for careful molecular and functional analysis. Comparing the transcriptomes of spinal cord- and cerebral-cortex-derived astrocytes reveals profound differences, including developmental patterning cues maintained in vitro. These relate to the distinct neuronal identity elicited by Ascl1 and Neurog2 reflecting their developmental functions in subtype specification of the respective CNS region.

Maturational Changes in Mouse Cutaneous Touch and Piezo2-Mediated Mechanotransduction

The age of studied animals has a profound impact on experimental outcomes in animal-based research. In mice, age influences molecular, morphological, physiological, and behavioral parameters, particularly during rapid postnatal growth and maturation until adulthood (at 12 weeks of age). Despite this knowledge, most biomedical studies use a wide-spanning age range from 4 to 12 weeks, raising concerns about reproducibility and potential masking of relevant age differences. Here, using mouse behavior and electrophysiology in cultured dorsal root ganglia (DRG), we reveal a decline in behavioral cutaneous touch sensitivity and Piezo2-mediated mechanotransduction in vitro during mouse maturation but not thereafter. In addition, we identify distinct transcript changes in individual Piezo2-expressing mechanosensitive DRG neurons by combining electrophysiology with single-cell RNA sequencing (patch-seq). Taken together, our study emphasizes the need for accurate age matching and uncovers hitherto unknown maturational plasticity in cutaneous touch at the level of behavior, mechanotransduction, and transcripts.

Patch-Seq Links Single-Cell Transcriptomes to Human Islet Dysfunction in Diabetes

Impaired function of pancreatic islet cells is a major cause of metabolic dysregulation and disease in humans. Despite this, it remains challenging to directly link physiological dysfunction in islet cells to precise changes in gene expression. Here we show that single-cell RNA sequencing combined with electrophysiological measurements of exocytosis and channel activity (patch-seq) can be used to link endocrine physiology and transcriptomes at the single-cell level. We collected 1,369 patch-seq cells from the pancreata of 34 human donors with and without diabetes. An analysis of function and gene expression networks identified a gene set associated with functional heterogeneity in β cells that can be used to predict electrophysiology. We also report transcriptional programs underlying dysfunction in type 2 diabetes and extend this approach to cryopreserved cells from donors with type 1 diabetes, generating a valuable resource for understanding islet cell heterogeneity in health and disease.

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

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 light-sensitive (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.

Corrigendum: Patch-Seq Protocol to Analyze the Electrophysiology, Morphology and Transcriptome of Whole Single Neurons Derived From Human Pluripotent Stem Cells

[This corrects the article DOI: 10.3389/fnmol.2018.00261.].

Patch-Seq Protocol to Analyze the Electrophysiology, Morphology and Transcriptome of Whole Single Neurons Derived From Human Pluripotent Stem Cells

The human brain is composed of a complex assembly of about 171 billion heterogeneous cellular units (86 billion neurons and 85 billion non-neuronal glia cells). A comprehensive description of brain cells is necessary to understand the nervous system in health and disease. Recently, advances in genomics have permitted the accurate analysis of the full transcriptome of single cells (scRNA-seq). We have built upon such technical progress to combine scRNA-seq with patch-clamping electrophysiological recording and morphological analysis of single human neurons in vitro. This new powerful method, referred to as Patch-seq, enables a thorough, multimodal profiling of neurons and permits us to expose the links between functional properties, morphology, and gene expression. Here, we present a detailed Patch-seq protocol for isolating single neurons from in vitro neuronal cultures. We have validated the Patch-seq whole-transcriptome profiling method with human neurons generated from embryonic and induced pluripotent stem cells (ESCs/iPSCs) derived from healthy subjects, but the procedure may be applied to any kind of cell type in vitro. Patch-seq may be used on neurons in vitro to profile cell types and states in depth to unravel the human molecular basis of neuronal diversity and investigate the cellular mechanisms underlying brain disorders.

Multimodal profiling of single-cell morphology, electrophysiology, and gene expression using Patch-seq

Neurons exhibit a rich diversity of morphological phenotypes, electrophysiological properties, and gene-expression patterns. Understanding how these different characteristics are interrelated at the single-cell level has been difficult because of the lack of techniques for multimodal profiling of individual cells. We recently developed Patch-seq, a technique that combines whole-cell patch-clamp recording, immunohistochemistry, and single-cell RNA-sequencing (scRNA-seq) to comprehensively profile single neurons from mouse brain slices. Here, we present a detailed step-by-step protocol, including modifications to the patching mechanics and recording procedure, reagents and recipes, procedures for immunohistochemistry, and other tips to assist researchers in obtaining high-quality morphological, electrophysiological, and transcriptomic data from single neurons. Successful implementation of Patch-seq allows researchers to explore the multidimensional phenotypic variability among neurons and to correlate gene expression with phenotype at the level of single cells. The entire procedure can be completed in ∼2 weeks through the combined efforts of a skilled electrophysiologist, molecular biologist, and biostatistician.