Emergence of the erythroid lineage from multipotent hematopoiesis [preprint]

Red cell formation begins with the hematopoietic stem cell, but the manner by which it gives rise to erythroid progenitors, and their subsequent developmental path, remain unclear. Here we combined single-cell transcriptomics of murine hematopoietic tissues with fate potential assays to infer a continuous yet hierarchical structure for the hematopoietic network. We define the erythroid differentiation trajectory as it emerges from multipotency and diverges from 6 other blood lineages. With the aid of a new flow-cytometric sorting strategy, we validated predicted cell fate potentials at the single cell level, revealing a coupling between erythroid and basophil/mast cell fates. We uncovered novel growth factor receptor regulators of the erythroid trajectory, including the proinflammatory IL- 17RA, found to be a strong erythroid stimulator; and identified a global hematopoietic response to stress erythropoiesis. We further identified transcriptional and high-purity FACS gates for the complete isolation of all classically-defined erythroid burst-forming (BFU-e) and colony-forming progenitors (CFU-e), finding that they express a dedicated transcriptional program, distinct from that of terminally-differentiating erythroblasts. Intriguingly, profound remodeling of the cell cycle is intimately entwined with CFU-e developmental progression and with a sharp transcriptional switch that extinguishes the CFU-e stage and activates terminal differentiation. Underlying these results, our work showcases the utility of theoretic approaches linking transcriptomic data to predictive fate models, providing key insights into lineage development in vivo

End Sequence Analysis Toolkit (ESAT) expands the extractable information from single-cell RNA-seq data

RNA-seq protocols that focus on transcript termini are well suited for applications in which template quantity is limiting. Here we show that, when applied to end-sequencing data, analytical methods designed for global RNA-seq produce computational artifacts. To remedy this, we created the End Sequence Analysis Toolkit (ESAT). As a test, we first compared end-sequencing and bulk RNA-seq using RNA from dendritic cells stimulated with lipopolysaccharide (LPS). As predicted by the telescripting model for transcriptional bursts, ESAT detected an LPS-stimulated shift to shorter 3\u27-isoforms that was not evident by conventional computational methods. Then, droplet-based microfluidics was used to generate 1000 cDNA libraries, each from an individual pancreatic islet cell. ESAT identified nine distinct cell types, three distinct beta-cell types, and a complex interplay between hormone secretion and vascularization. ESAT, then, offers a much-needed and generally applicable computational pipeline for either bulk or single-cell RNA end-sequencing

Single-Cell Analysis of Experience-Dependent Transcriptomic States in Mouse Visual Cortex

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

Identification of a targetable KRAS-mutant epithelial population in non-small cell lung cancer

Lung cancer is the leading cause of cancer deaths. Tumor heterogeneity, which hampers development of targeted therapies, was herein deconvoluted via single cell RNA sequencingin aggressive human adenocarcinomas (carrying Kras-mutations) and comparable murine model. We identified a tumor-specific, mutant-KRAS-associated subpopulation which is conserved in both human and murine lung cancer. We previously reported a key role for the oncogene BMI-1 in adenocarcinomas. We therefore investigated the effects of in vivo PTC596 treatment, which affects BMI-1 activity, in our murine model. Post-treatment, MRI analysis showed decreased tumor size, while single cell transcriptomics concomitantly detected near complete ablation of the mutant-KRAS-associated subpopulation, signifying the presence of a pharmacologically targetable, tumor-associated subpopulation. Our findings therefore hold promise for the development of a targeted therapy for KRAS-mutant adenocarcinomas

Quantitative Clonal Analysis and Single-Cell Transcriptomics Reveal Division Kinetics, Hierarchy, and Fate of Oral Epithelial Progenitor Cells

The oral mucosa is one of the most rapidly dividing tissues in the body and serves as a barrier to physical and chemical insults from mastication, food, and microorganisms. Breakdown of this barrier can lead to significant morbidity and potentially life-threatening infections for patients. Determining the identity and organization of oral epithelial progenitor cells (OEPCs) is therefore paramount to understanding their roles in homeostasis and disease. Using lineage tracing and label retention experiments, we show that rapidly dividing OEPCs are located broadly within the basal layer of the mucosa throughout the oral cavity. Quantitative clonal analysis demonstrated that OEPCs undergo population-asymmetrical divisions following neutral drift dynamics and that they respond to chemotherapy-induced damage by altering daughter cell fates. Finally, using single-cell RNA-seq, we establish the basal layer population structure and propose a model that defines the organization of cells within the basal layer

A single-cell atlas of the airway epithelium reveals the CFTR-rich pulmonary ionocyte

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 clearance1. 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

A large pool of actively cycling progenitors orchestrates self-renewal and injury repair of an ectodermal appendage.

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

Single-cell transcriptomics of human and mouse lung cancers reveals conserved myeloid populations across individuals and species

Tumor-infiltrating myeloid cells (TIMs) comprise monocytes, macrophages, dendritic cells, and neutrophils, and have emerged as key regulators of cancer growth. These cells can diversify into a spectrum of states, which might promote or limit tumor outgrowth but remain poorly understood. Here, we used single-cell RNA sequencing (scRNA-seq) to map TIMs in non-small-cell lung cancer patients. We uncovered 25 TIM states, most of which were reproducibly found across patients. To facilitate translational research of these populations, we also profiled TIMs in mice. In comparing TIMs across species, we identified a near-complete congruence of population structures among dendritic cells and monocytes; conserved neutrophil subsets; and species differences among macrophages. By contrast, myeloid cell population structures in patients' blood showed limited overlap with those of TIMs. This study determines the lung TIM landscape and sets the stage for future investigations into the potential of TIMs as immunotherapy targets