Summary It has long been the dream of biologists to map gene expression at the single cell level. With such data one might track heterogeneous cell sub-populations, and infer regulatory relationships between genes and pathways. Recently, RNA sequencing has achieved single cell resolution. What is limiting is an effective way to routinely isolate and process large numbers of individual cells for quantitative in-depth sequencing. We have developed a high-throughput droplet-microfluidic approach for barcoding the RNA from thousands of individual cells for subsequent analysis by next-generation sequencing. The method shows a surprisingly low noise profile and is readily adaptable to other sequencing-based assays. We analyzed mouse embryonic stem cells, revealing in detail the population structure and the heterogeneous onset of differentiation after LIF withdrawal. The reproducibility of these high-throughput single cell data allowed us to deconstruct cell populations and infer gene expression relationships.

Intestinal stem cells, characterized by high Lgr5 expression, reside between Paneth cells at the small intestinal crypt base and divide every day. We have carried out fate mapping of individual stem cells by generating a multicolor Cre-reporter. As a population, Lgr5(hi) stem cells persist life-long, yet crypts drift toward clonality within a period of 1-6 months. We have collected short- and long-term clonal tracing data of individual Lgr5(hi) cells. These reveal that most Lgr5(hi) cell divisions occur symmetrically and do not support a model in which two daughter cells resulting from an Lgr5(hi) cell division adopt divergent fates (i.e., one Lgr5(hi) cell and one transit-amplifying [TA] cell per division). The cellular dynamics are consistent with a model in which the resident stem cells double their numbers each day and stochastically adopt stem or TA fates. Quantitative analysis shows that stem cell turnover follows a pattern of neutral drift dynamics.

Summary Although the function of the mammalian pancreas hinges on complex interactions of distinct cell types, gene expression profiles have primarily been described with bulk mixtures. Here we implemented a droplet-based, single-cell RNA-seq method to determine the transcriptomes of over 12,000 individual pancreatic cells from four human donors and two mouse strains. Cells could be divided into 15 clusters that matched previously characterized cell types: all endocrine cell types, including rare epsilon-cells; exocrine cell types; vascular cells; Schwann cells; quiescent and activated stellate cells; and four types of immune cells. We detected subpopulations of ductal cells with distinct expression profiles and validated their existence with immuno-histochemistry stains. Moreover, among human beta- cells, we detected heterogeneity in the regulation of genes relating to functional maturation and levels of ER stress. Finally, we deconvolved bulk gene expression samples using the single-cell data to detect disease-associated differential expression. Our dataset provides a resource for the discovery of novel cell type-specific transcription factors, signaling receptors, and medically relevant genes.

Highlights d Human dendritic cell and monocyte subsets show one-toone equivalence in mouse d Neutrophils exhibit tumor-associated phenotypes that are conserved across species d Myeloid subsets in patient blood only partially overlap with those in their tumors d Unique markers define myeloid cell subsets and associate with clinical prognosis

According to the current model of adult epidermal homeostasis, skin tissue is maintained by two discrete populations of progenitor cells: self-renewing stem cells; and their progeny, known as transit amplifying cells, which differentiate after several rounds of cell division. By making use of inducible genetic labelling, we have tracked the fate of a representative sample of progenitor cells in mouse tail epidermis at single-cell resolution in vivo at time intervals up to one year. Here we show that clone-size distributions are consistent with a new model of homeostasis involving only one type of progenitor cell. These cells are found to undergo both symmetric and asymmetric division at rates that ensure epidermal homeostasis. The results raise important questions about the potential role of stem cells on tissue maintenance in vivo.

The functions of epithelial tissues are dictated by the types, abundance and distribution of the differentiated cells they contain. Attempts to restore tissue function after damage require knowledge of how physiological tasks are distributed among cell types, and how cell states vary between homeostasis, injury-repair and disease. In the conducting airway, a heterogeneous basal cell population gives rise to specialized luminal cells that perform mucociliary clearance. Here we perform single-cell profiling of human bronchial epithelial cells and mouse tracheal epithelial cells to obtain a comprehensive census of cell types in the conducting airway and their behaviour in homeostasis and regeneration. Our analysis reveals cell states that represent known and novel cell populations, delineates their heterogeneity and identifies distinct differentiation trajectories during homeostasis and tissue repair. Finally, we identified a novel, rare cell type that we call the 'pulmonary ionocyte', which co-expresses FOXI1, multiple subunits of the vacuolar-type H-ATPase (V-ATPase) and CFTR, the gene that is mutated in cystic fibrosis. Using immunofluorescence, modulation of signalling pathways and electrophysiology, we show that Notch signalling is necessary and FOXI1 expression is sufficient to drive the production of the pulmonary ionocyte, and that the pulmonary ionocyte is a major source of CFTR activity in the conducting airway epithelium.

Highlights d We define two multiplet errors in single-cell RNA-seq data: ''embedded'' and ''neotypic'' d Neotypic errors can lead to misidentification of cell types or transitional states d Scrublet code identifies neotypic doublets and predicts the overall doublet rate d The algorithm is tested against several experimental methods for labeling multiplets

It is currently thought that life-long blood cell production is driven by the action of a small number of multipotent haematopoietic stem cells. Evidence supporting this view has been largely acquired through the use of functional assays involving transplantation. However, whether these mechanisms also govern native non-transplant haematopoiesis is entirely unclear. Here we have established a novel experimental model in mice where cells can be uniquely and genetically labelled in situ to address this question. Using this approach, we have performed longitudinal analyses of clonal dynamics in adult mice that reveal unprecedented features of native haematopoiesis. In contrast to what occurs following transplantation, steady-state blood production is maintained by the successive recruitment of thousands of clones, each with a minimal contribution to mature progeny. Our results demonstrate that a large number of long-lived progenitors, rather than classically defined haematopoietic stem cells, are the main drivers of steady-state haematopoiesis during most of adulthood. Our results also have implications for understanding the cellular origin of haematopoietic disease.