For more than 100 years, the fruit fly Drosophila melanogaster has been one of the most studied model organisms. Here, we present a single-cell atlas of the adult fly, Tabula Drosophilae , that includes 580,000 nuclei from 15 individually dissected sexed tissues as well as the entire head and body, annotated to >250 distinct cell types. We provide an in-depth analysis of cell type–related gene signatures and transcription factor markers, as well as sexual dimorphism, across the whole animal. Analysis of common cell types between tissues, such as blood and muscle cells, reveals rare cell types and tissue-specific subtypes. This atlas provides a valuable resource for the Drosophila community and serves as a reference to study genetic perturbations and disease models at single-cell resolution.

Tyrosine kinase inhibitors were found to be clinically effective for treatment of patients with certain subsets of cancers carrying somatic mutations in receptor tyrosine kinases. However, the duration of clinical response is often limited, and patients ultimately develop drug resistance. Here, we use single-cell RNA sequencing to demonstrate the existence of multiple cancer cell subpopulations within cell lines, xenograft tumors and patient tumors. These subpopulations exhibit epigenetic changes and differential therapeutic sensitivity. Recurrently overrepresented ontologies in genes that are differentially expressed between drug tolerant cell populations and drug sensitive cells include epithelial-to-mesenchymal transition, epithelium development, vesicle mediated transport, drug metabolism and cholesterol homeostasis. We show analysis of identified markers using the LINCS database to predict and functionally validate small molecules that target selected drug tolerant cell populations. In combination with EGFR inhibitors, crizotinib inhibits the emergence of a defined subset of EGFR inhibitor-tolerant clones. In this study, we describe the spectrum of changes associated with drug tolerance and inhibition of specific tolerant cell subpopulations with combination agents.

The function of Retinoblastoma tumor suppressor (pRB) is greatly influenced by the cellular context, therefore the consequences of pRB inactivation are cell-type-specific. Here we employ single cell RNA-sequencing (scRNA-seq) to profile the impact of an Rbf mutation during Drosophila eye development. First, we build a catalogue of 11,500 wild type eye disc cells containing major known cell types. We find a transcriptional switch occurring in differentiating photoreceptors at the time of axonogenesis. Next, we map a cell landscape of Rbf mutant and identify a mutant-specific cell population that shows intracellular acidification due to increase in glycolytic activity. Genetic experiments demonstrate that such metabolic changes, restricted to this unique Rbf mutant population, sensitize cells to apoptosis and define the pattern of cell death in Rbf mutant eye disc. Thus, these results illustrate how scRNA-seq can be applied to dissect mutant phenotypes.

In Drosophila, the wing disc-associated muscle precursor cells give rise to the fibrillar indirect flight muscles (IFM) and the tubular direct flight muscles (DFM). To understand early transcriptional events underlying this muscle diversification, we performed singlecell RNA-sequencing experiments and built a cell atlas of myoblasts associated with third instar larval wing disc. Our analysis identified distinct transcriptional signatures for IFM and DFM myoblasts that underlie the molecular basis of their divergence. The atlas further revealed various states of differentiation of myoblasts, thus illustrating previously unappreciated spatial and temporal heterogeneity among them. We identified and validated novel markers for both IFM and DFM myoblasts at various states of differentiation by immunofluorescence and genetic cell-tracing experiments. Finally, we performed a systematic genetic screen using a panel of markers from the reference cell atlas as an entry point and found a novel gene, Amalgam which is functionally important in muscle development. Our work provides a framework for leveraging scRNA-seq for gene discovery and details a strategy that can be applied to other scRNA-seq datasets.