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

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.

Single-cell RNA sequencing has recently emerged as a powerful tool for mapping cellular heterogeneity in diseased and healthy tissues, yet high-throughput methods are needed for capturing the unbiased diversity of cells. Droplet microfluidics is among the most promising candidates for capturing and processing thousands of individual cells for whole-transcriptome or genomic analysis in a massively parallel manner with minimal reagent use. We recently established a method called inDrops, which has the capability to index >15,000 cells in an hour. A suspension of cells is first encapsulated into nanoliter droplets with hydrogel beads (HBs) bearing barcoding DNA primers. Cells are then lysed and mRNA is barcoded (indexed) by a reverse transcription (RT) reaction. Here we provide details for (i) establishing an inDrops platform (1 d); (ii) performing hydrogel bead synthesis (4 d); (iii) encapsulating and barcoding cells (1 d); and (iv) RNA-seq library preparation (2 d). inDrops is a robust and scalable platform, and it is unique in its ability to capture and profile >75% of cells in even very small samples, on a scale of thousands or tens of thousands of cells.

Activity-dependent transcriptional responses shape cortical function. However, we lack a comprehensive understanding of the diversity of these responses across the full range of cortical cell types, and how these changes contribute to neuronal plasticity and disease. Here we applied high-throughput single-cell RNA-sequencing to investigate the breadth of transcriptional changes that occur across cell types in mouse visual cortex following exposure to light. We identified significant and divergent transcriptional responses to stimulation in each of the 30 cell types characterized, revealing 611 stimulus-responsive genes. Excitatory pyramidal neurons exhibit inter- and intra-laminar heterogeneity in the induction of stimulus responsive genes. Non-neuronal cells demonstrated clear transcriptional responses that may regulate experience-dependent changes in neurovascular coupling and myelination. Together, these results reveal the dynamic landscape of stimulus-dependent transcriptional changes that occur across cell types in visual cortex, which are likely critical for cortical function and may be sites of de-regulation in developmental brain disorders.

Red cell formation begins with the differentiation of multipotent hematopoietic progenitors. Reconstructing the steps of differentiation represents a stereotypical challenge in stem cell biology. Combining single-cell transcriptomics, fate assays, and theory for predicting fate from population snapshots, we inferred a continuous, hierarchical structure of murine hematopoietic progenitors committing to seven blood lineages. We uncovered coupling between erythroid and basophil/mast cell fates, a global hematopoietic response to erythroid stress, and novel growth factor receptor regulators of erythropoiesis. We also defined a new flow-cytometric sorting strategy to purify progressive early stages of erythroid differentiation, completely isolating classically-defined burst-forming (BFU-e) and colony-forming progenitors (CFU-e). Intriguingly, profound remodeling of the cell cycle is intimately entwined with erythroid development and with a sharp transcriptional switch that extinguishes the CFU-e stage and activates terminal differentiation. Our work showcases the utility of theory linking transcriptomic data to predictive fate models, providing insights into lineage development in vivo.

CD4+ T regulatory (Treg) cells are central to immune homeostasis, their phenotypic heterogeneity reflecting the diverse environments and target cells they regulate. To understand this heterogeneity, we combined single-cell RNAseq, activation reporter and TCR analysis to profile thousands of Tregs or Tconvs from mouse lymphoid organs or human blood. Treg and Tconv pools showed areas of overlap, as resting “furtive” Tregs with overall similarity to Tconv, or as a convergence of activated states. All Tregs express a small core of FoxP3-dependent transcripts, onto which additional programs are added less uniformly. Among suppressive functions, Il2ra and Ctla4 were quasi-constant, inhibitory cytokines being more sparsely distributed. TCR signal intensity didn’t affect resting/activated Treg proportions, but molded activated Treg programs. The main lines of Treg heterogeneity in mice were strikingly conserved in human blood. These results reveal unexpected TCR-shaped states of activation, providing a framework to synthesize previous observations about Treg heterogeneity.

Bone marrow-derived myeloid cells can accumulate within tumors and foster cancer outgrowth. Local immune-neoplastic interactions have been intensively investigated, but the contribution of the systemic host environment to tumor growth remains poorly understood. Here, we show in mice and cancer patients ( = 70) that lung adenocarcinomas increase bone stromal activity in the absence of bone metastasis. Animal studies reveal that the cancer-induced bone phenotype involves bone-resident osteocalcin-expressing (Ocn) osteoblastic cells. These cells promote cancer by remotely supplying a distinct subset of tumor-infiltrating SiglecF neutrophils, which exhibit cancer-promoting properties. Experimentally reducing Ocn cell numbers suppresses the neutrophil response and lung tumor outgrowth. These observations posit osteoblasts as remote regulators of lung cancer and identify SiglecF neutrophils as myeloid cell effectors of the osteoblast-driven protumoral response.

The classical model of tissue renewal posits that small numbers of quiescent stem cells (SCs) give rise to proliferating transit-amplifying cells before terminal differentiation. However, many organs house pools of SCs with proliferative and differentiation potentials that diverge from this template. Resolving SC identity and organization is therefore central to understanding tissue renewal. Here, using a combination of single-cell RNA sequencing (scRNA-seq), mouse genetics and tissue injury approaches, we uncover cellular hierarchies and mechanisms that underlie the maintenance and repair of the continuously growing mouse incisor. Our results reveal that, during homeostasis, a group of actively cycling epithelial progenitors generates enamel-producing ameloblasts and adjacent layers of non-ameloblast cells. After injury, tissue repair was achieved through transient increases in progenitor-cell proliferation and through direct conversion of Notch1-expressing cells to ameloblasts. We elucidate epithelial SC identity, position and function, providing a mechanistic basis for the homeostasis and repair of a fast-turnover ectodermal appendage.