Molecular Architecture of the Mouse Nervous System.

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.

Molecular Diversity of Midbrain Development in Mouse, Human, and Stem Cells.

Understanding human embryonic ventral midbrain is of major interest for Parkinson's disease. However, the cell types, their gene expression dynamics, and their relationship to commonly used rodent models remain to be defined. We performed single-cell RNA sequencing to examine ventral midbrain development in human and mouse. We found 25 molecularly defined human cell types, including five subtypes of radial glia-like cells and four progenitors. In the mouse, two mature fetal dopaminergic neuron subtypes diversified into five adult classes during postnatal development. Cell types and gene expression were generally conserved across species, but with clear differences in cell proliferation, developmental timing, and dopaminergic neuron development. Additionally, we developed a method to quantitatively assess the fidelity of dopaminergic neurons derived from human pluripotent stem cells, at a single-cell level. Thus, our study provides insight into the molecular programs controlling human midbrain development and provides a foundation for the development of cell replacement therapies.

Single-Cell Profiling Uncovers Evolutionary Divergence of Hypocretin/Orexin Neuronal Subpopulations.

Brain nuclei are traditionally defined by their anatomy, activity, and expression of specific markers. The hypothalamus contains discrete neuronal populations that coordinate fundamental behavioral functions, including sleep and wakefulness, in all vertebrates. Particularly, the diverse roles of hypocretin/orexin (Hcrt)-releasing neurons suggest functional heterogeneity among Hcrt neurons. Using single-cell RNA sequencing (scRNA-seq) and high-resolution imaging of the adult male and female zebrafish hypothalamic periventricular zone, we identified 21 glutamatergic and 28 GABAergic cell types. Integration of zebrafish and mouse scRNA-seq revealed evolutionary conserved and divergent hypothalamic cell types. The expression of specific genes, including npvf, which encodes a sleep-regulating neuropeptide, was enriched in subsets of glutamatergic Hcrt neurons in both larval and adult zebrafish. The genetic profile, activity, and neurite processing of the neuronal subpopulation that coexpresses both Hcrt and Npvf (Hcrt+Npvf+) differ from other Hcrt neurons. These interspecies findings provide a unified annotation of hypothalamic cell types and suggest that the heterogeneity of Hcrt neurons enables multifunctionality, such as consolidation of both wake and sleep by the Hcrt- and Npvf-releasing neuronal subpopulation.

RNA velocity of single cells.

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.

Spatial organization of the somatosensory cortex revealed by osmFISH.

Global efforts to create a molecular census of the brain using single-cell transcriptomics are producing a large catalog of molecularly defined cell types. However, spatial information is lacking and new methods are needed to map a large number of cell type-specific markers simultaneously on large tissue areas. Here, we describe a cyclic single-molecule fluorescence in situ hybridization methodology and define the cellular organization of the somatosensory cortex.

A secretagogin locus of the mammalian hypothalamus controls stress hormone release.

A hierarchical hormonal cascade along the hypothalamic-pituitary-adrenal axis orchestrates bodily responses to stress. Although corticotropin-releasing hormone (CRH), produced by parvocellular neurons of the hypothalamic paraventricular nucleus (PVN) and released into the portal circulation at the median eminence, is known to prime downstream hormone release, the molecular mechanism regulating phasic CRH release remains poorly understood. Here, we find a cohort of parvocellular cells interspersed with magnocellular PVN neurons expressing secretagogin. Single-cell transcriptome analysis combined with protein interactome profiling identifies secretagogin neurons as a distinct CRH-releasing neuron population reliant on secretagogin's Ca(2+) sensor properties and protein interactions with the vesicular traffic and exocytosis release machineries to liberate this key hypothalamic releasing hormone. Pharmacological tools combined with RNA interference demonstrate that secretagogin's loss of function occludes adrenocorticotropic hormone release from the pituitary and lowers peripheral corticosterone levels in response to acute stress. Cumulatively, these data define a novel secretagogin neuronal locus and molecular axis underpinning stress responsiveness.

Quantitative single-cell RNA-seq with unique molecular identifiers.

Single-cell RNA sequencing (RNA-seq) is a powerful tool to reveal cellular heterogeneity, discover new cell types and characterize tumor microevolution. However, losses in cDNA synthesis and bias in cDNA amplification lead to severe quantitative errors. We show that molecular labels--random sequences that label individual molecules--can nearly eliminate amplification noise, and that microfluidic sample preparation and optimized reagents produce a fivefold improvement in mRNA capture efficiency.

High-resolution microbial community reconstruction by integrating short reads from multiple 16S rRNA regions.

The emergence of massively parallel sequencing technology has revolutionized microbial profiling, allowing the unprecedented comparison of microbial diversity across time and space in a wide range of host-associated and environmental ecosystems. Although the high-throughput nature of such methods enables the detection of low-frequency bacteria, these advances come at the cost of sequencing read length, limiting the phylogenetic resolution possible by current methods. Here, we present a generic approach for integrating short reads from large genomic regions, thus enabling phylogenetic resolution far exceeding current methods. The approach is based on a mapping to a statistical model that is later solved as a constrained optimization problem. We demonstrate the utility of this method by analyzing human saliva and Drosophila samples, using Illumina single-end sequencing of a 750 bp amplicon of the 16S rRNA gene. Phylogenetic resolution is significantly extended while reducing the number of falsely detected bacteria, as compared with standard single-region Roche 454 Pyrosequencing. Our approach can be seamlessly applied to simultaneous sequencing of multiple genes providing a higher resolution view of the composition and activity of complex microbial communities.

An accessible database for mouse and human whole transcriptome qPCR primers.

MOTIVATION: Real time quantitative polymerase chain reaction (qPCR) is an important tool in quantitative studies of DNA and RNA molecules; especially in transcriptome studies, where different primer combinations allow identification of specific transcripts such as splice variants or precursor messenger RNA. Several softwares that implement various rules for optimal primer design are available. Nevertheless, as designing qPCR primers needs to be done manually, the repeated task is tedious, time consuming and prone to errors. RESULTS: We used a set of rules to automatically design all possible exon-exon and intron-exon junctions in the human and mouse transcriptomes. The resulting database is included as a track in the UCSC genome browser, making it widely accessible and easy to use. AVAILABILITY: The database is available from the UCSC genome browser (http://genome.ucsc.edu/), track name 'Whole Transcriptome qPCR Primers' for the hg19 (Human) and mm10 (Mouse) genome versions. Batch query is available in the following: http://www.weizmann.ac.il/complex/compphys/software/Amit/primers/batch_query_qpcr_primers.htm CONTACT: amit.zeisel@weizmann.ac.il or eytan.domany@weizmann.ac.il SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Blastocyst implantation failure relates to impaired translational machinery gene expression.

Oocyte quality is a well-established determinant of embryonic fate. However, the molecular participants and biological markers that affect and may predict adequate embryonic development are largely elusive. Our aim was to identify the components of the oocyte molecular machinery that part take in the production of a healthy embryo. For this purpose, we used an animal model, generated by us previously, the oocytes of which do not express Cx43 (Cx43(del/del)). In these mice, oogenesis appears normal, fertilisation does occur, early embryonic development is successful but implantation fails. We used magnetic resonance imaging analysis combined with histological examination to characterise the embryonic developmental incompetence. Reciprocal embryo transfer confirmed that the blastocyst evolved from the Cx43(del/del) oocyte is responsible for the implantation disorder. In order to unveil the genes, the impaired expression of which brings about the development of defective embryos, we carried out a genomic screening of both the oocytes and the resulting blastocysts. This microarray analysis revealed a low expression of Egr1, Rpl21 and Eif4a1 in Cx43(del/del) oocytes and downregulation of Rpl15 and Eif4g2 in the resulting blastocysts. We propose that global deficiencies in genes related to the expression of ribosomal proteins and translation initiation factors in apparently normal oocytes bring about accumulation of defects, which significantly compromise their developmental capacity. The blastocysts resulting from such oocytes, which grow within a confined space until implantation, may be unable to generate enough biological mass to allow their expansion. This information could be implicated to diagnosis and treatment of infertility, particularly to IVF.