The mammalian nervous system executes complex behaviors controlled by specialized, precisely positioned, and interacting cell types. Here, we used RNA sequencing of half a million single cells to create a detailed census of cell types in the mouse nervous system. We mapped cell types spatially and derived a hierarchical, data-driven taxonomy. Neurons were the most diverse and were grouped by developmental anatomical units and by the expression of neurotransmitters and neuropeptides. Neuronal diversity was driven by genes encoding cell identity, synaptic connectivity, neurotransmission, and membrane conductance. We discovered seven distinct, regionally restricted astrocyte types that obeyed developmental boundaries and correlated with the spatial distribution of key glutamate and glycine neurotransmitters. In contrast, oligodendrocytes showed a loss of regional identity followed by a secondary diversification. The resource presented here lays a solid foundation for understanding the molecular architecture of the mammalian nervous system and enables genetic manipulation of specific cell types.

RNA abundance is a powerful indicator of the state of individual cells. Single-cell RNA sequencing can reveal RNA abundance with high quantitative accuracy, sensitivity and throughput1. However, this approach captures only a static snapshot at a point in time, posing a challenge for the analysis of time-resolved phenomena such as embryogenesis or tissue regeneration. Here we show that RNA velocity-the time derivative of the gene expression state-can be directly estimated by distinguishing between unspliced and spliced mRNAs in common single-cell RNA sequencing protocols. RNA velocity is a high-dimensional vector that predicts the future state of individual cells on a timescale of hours. We validate its accuracy in the neural crest lineage, demonstrate its use on multiple published datasets and technical platforms, reveal the branching lineage tree of the developing mouse hippocampus, and examine the kinetics of transcription in human embryonic brain. We expect RNA velocity to greatly aid the analysis of developmental lineages and cellular dynamics, particularly in humans.

Mosquito blood cells are immune cells that help control infection by vector-borne pathogens. Despite their importance, little is known about mosquito blood cell biology beyond morphological and functional criteria used for their classification. Here, we combined the power of single-cell RNA sequencing, high-content imaging flow cytometry, and single-molecule RNA hybridization to analyze a subset of blood cells of the malaria mosquito Anopheles gambiae By demonstrating that blood cells express nearly half of the mosquito transcriptome, our dataset represents an unprecedented view into their transcriptional program. Analyses of differentially expressed genes identified transcriptional signatures of two cell types and provide insights into the current classification of these cells. We further demonstrate the active transfer of a cellular marker between blood cells that may confound their identification. We propose that cell-to-cell exchange may contribute to cellular diversity and functional plasticity seen across biological systems.

The mammary gland consists of cells with gene expression patterns reflecting their cellular origins, function, and spatiotemporal context. However, knowledge of developmental kinetics and mechanisms of lineage specification is lacking. We address this significant knowledge gap by generating a single-cell transcriptome atlas encompassing embryonic, postnatal, and adult mouse mammary development. From these data, we map the chronology of transcriptionally and epigenetically distinct cell states and distinguish fetal mammary stem cells (fMaSCs) from their precursors and progeny. fMaSCs show balanced co-expression of factors associated with discrete adult lineages and a metabolic gene signature that subsides during maturation but reemerges in some human breast cancers and metastases. These data provide a useful resource for illuminating mammary cell heterogeneity, the kinetics of differentiation, and developmental correlates of tumorigenesis.

Recent sequencing analyses have shed light on heterogeneous patterns of genomic aberrations in human gastric cancers (GCs). To explore how individual genetic events translate into cancer phenotypes, we established a biological library consisting of genetically engineered gastric organoids carrying various GC mutations and 37 patient-derived organoid lines, including rare genomically stable GCs. Phenotype analyses of GC organoids revealed divergent genetic and epigenetic routes to gain Wnt and R-spondin niche independency. An unbiased phenotype-based genetic screening identified a significant association between CDH1/TP53 compound mutations and the R-spondin independency that was functionally validated by CRISPR-based knockout. Xenografting of GC organoids further established the feasibility of Wnt-targeting therapy for Wnt-dependent GCs. Our results collectively demonstrate that multifaceted genetic abnormalities render human GCs independent of the stem cell niche and highlight the validity of the genotype-phenotype screening strategy in gaining deeper understanding of human cancers.

The mammalian brain is composed of diverse, specialized cell populations. To systematically ascertain and learn from these cellular specializations, we used Drop-seq to profile RNA expression in 690,000 individual cells sampled from 9 regions of the adult mouse brain. We identified 565 transcriptionally distinct groups of cells using computational approaches developed to distinguish biological from technical signals. Cross-region analysis of these 565 cell populations revealed features of brain organization, including a gene-expression module for synthesizing axonal and presynaptic components, patterns in the co-deployment of voltage-gated ion channels, functional distinctions among the cells of the vasculature and specialization of glutamatergic neurons across cortical regions. Systematic neuronal classifications for two complex basal ganglia nuclei and the striatum revealed a rare population of spiny projection neurons. This adult mouse brain cell atlas, accessible through interactive online software (DropViz), serves as a reference for development, disease, and evolution.

Interleukin-22 (IL-22)-producing group 3 innate lymphoid cells (ILC3) maintains gut homeostasis but can also promote inflammatory bowel disease (IBD). The regulation of ILC3-dependent colitis remains to be elucidated. Here we show that Foxp3 + regulatory T cells (Treg cells) prevented ILC3-mediated colitis in an IL-10-independent manner. Treg cells inhibited IL-23 and IL-1β production from intestinal-resident CX3CR1 + macrophages but not CD103 + dendritic cells. Moreover, Treg cells restrained ILC3 production of IL-22 through suppression of CX3CR1 + macrophage production of IL-23 and IL-1β. This suppression was contact dependent and was mediated by latent activation gene-3 (LAG-3)—an immune checkpoint receptor—expressed on Treg cells. Engagement of LAG-3 on MHC class II drove profound immunosuppression of CX3CR1+ tissue-resident macrophages. Our study reveals that the health of the intestinal mucosa is maintained by an axis driven by Treg cells communication with resident macrophages that withhold inflammatory stimuli required for ILC3 function.

Most intronic RNAs are degraded within seconds or minutes after their excision from newly formed transcripts. However, stable intronic sequence RNAs (sisRNAs) have been described from oocytes of the frog Xenopus, from Drosophila embryos, and from human cell lines. In Xenopus oocytes, sisRNAs are abundant in both the nucleus and cytoplasm, they occur in the form of lariats, and they are stable for days. In this study we demonstrate that cytoplasmic sisRNAs are also found in human, mouse, chicken, and zebrafish cells. They exist as circular (lariat) molecules, mostly 100-500 nucleotides in length, and are derived from many housekeeping genes. They tend to have an unusual cytosine branchpoint (with the exception of those from the frog). Stable lariats are exported from the nucleus to the cytoplasm by the NXF1/NXT1 system, demonstrating that their presence in the cytoplasm is not due to passive diffusion. Lariats in the cytoplasm are not associated with transcripts of the genes from which they are derived. The biological significance of cytoplasmic sisRNAs remains obscure.

Our appreciation for the extent of Epstein Barr virus (EBV) transcriptome complexity continues to grow through findings of EBV encoded microRNAs, new long non-coding RNAs as well as the more recent discovery of over a hundred new polyadenylated lytic transcripts. Here we report an additional layer to the EBV transcriptome through the identification of a repertoire of latent and lytic viral circular RNAs. Utilizing RNase R-sequencing with cell models representing latency types I, II, and III, we identified EBV encoded circular RNAs expressed from the latency Cp promoter involving backsplicing from the W1 and W2 exons to the C1 exon, from the EBNA BamHI U fragment exon, and from the latency long non-coding RPMS1 locus. In addition, we identified circular RNAs expressed during reactivation including backsplicing from exon 8 to exon 2 of the LMP2 gene and a highly expressed circular RNA derived from intra-exonic backsplicing within the BHLF1 gene. While expression of most of these circular RNAs was found to depend on the EBV transcriptional program utilized and the transcription levels of the associated loci, expression of LMP2 exon 8 to exon 2 circular RNA was found to be cell model specific. Altogether we identified over 30 unique EBV circRNAs candidates and we validated and determined the structural features, expression profiles and nuclear/cytoplasmic distributions of several predominant and notable viral circRNAs. Further, we show that two of the EBV circular RNAs derived from the RPMS1 locus are detected in EBV positive clinical stomach cancer specimens. This study increases the known EBV latency and lytic transcriptome repertoires to include viral circular RNAs and it provides an essential foundation and resource for investigations into the functions and roles of this new class of EBV transcripts in EBV biology and diseases.

Telomerase is reactivated in most cancers and is possibly an early driver event in melanoma. Our aim was to test a novel in situ hybridization technique, RNAscope, for the detection of human telomerase reverse transcriptase (hTERT) mRNA in archival formalin-fixed, paraffin-embedded (FFPE) tissue and to compare the mRNA expression of melanomas and benign naevi. Furthermore, we wanted to see if hTERT mRNA could be a diagnostic or prognostic marker of melanoma. In situ hybridization for the detection of hTERT mRNA was performed on FFPE tissue of 17 melanomas and 13 benign naevi. We found a significant difference in the expression of hTERT mRNA between melanomas and benign naevi (P<0.001) and the expression of hTERT mRNA correlated with Breslow thickness (ρ=0.56, P=0.0205) and the Ki67 proliferation index (ρ=0.72, P=0.001). This study showed that RNAscope was a reliable in situ hybridization method for the detection of hTERT mRNA in FFPE tissue of melanomas and benign naevi. hTERT mRNA was more abundantly expressed in melanomas compared with benign naevi, but cannot be used solely as a diagnostic marker due to an overlap in expression. The hTERT mRNA expression in melanomas correlated with the prognostic markers Breslow thickness and the Ki67 index indicating a prognostic potential of hTERT mRNA.This is an open-access article distributed under the terms of the Creative Commons Attribution-Non Commercial-No Derivatives License 4.0 (CCBY-NC-ND), where it is permissible to download and share the work provided it is properly cited.

Chronic itch is a highly debilitating condition affecting about 10% of the general population. The relay of itch signals is under tight control by inhibitory circuits of the spinal dorsal horn, which may offer a hitherto unexploited therapeutic opportunity. Here, we found that specific pharmacological targeting of inhibitory α2 and α3GABAA receptors reduces acute histaminergic and non-histaminergic itch in mice. Systemic treatment with an α2/α3GABAA receptor selective modulator alleviates also chronic itch in a mouse model of atopic dermatitis and in dogs sensitized to house dust mites, without inducing sedation, motor dysfunction, or loss of antipruritic activity after prolonged treatment. Transsynaptic circuit tracing, immunofluorescence, and electrophysiological experiments identify spinal α2 and α3GABAA receptors as likely molecular targets underlying the antipruritic effect. Our results indicate that drugs targeting α2 and α3GABAA receptors are well-suited to alleviate itch, including non-histaminergic chronic itch for which currently no approved treatment exists.

The Hyperpolarization activated, cyclic nucleotide gated (HCN) channel is a candidate mediator of neuroendocrine influence over detrusor tonus during filling. In other tissues, HCN loss with aging is linked to declines in rhythmicity and function. We hypothesized that HCN has an age-sensitive expression profile and functional role in adrenergic bladder relaxation. HCN was examined in bladders from young (2–6 months) and old (18–24 months) C57BL/6 female mice, using qRT-PCR, RNAScope, and Western blots. Isometric tension studies were conducted using bladder strips from young wild-type (YWT), old wild-type (OWT), and young HCN1 knock-out (YKO) female mice to test the role HCN in effects of β-adrenergic stimulation. Hcn1 is the dominant HCN isoform RNA in the mouse bladder wall, and is diminished with age. Location of Hcn RNA within the mouse bladder wall is isoform-specific, with HCN1 limited to the detrusor layer. Passively-tensioned YWT bladder strips are relaxed by isoproterenol in the presence of HCN function, where OWT strips are relaxed only in the presence of HCN blockade. HCN has an age-specific expression and function in adrenergic detrusor relaxation in mouse bladder strips.