Our understanding of the development and maintenance of tissues has been greatly aided by large-scale gene expression analysis. However, tissues are invariably complex, and expression analysis of a tissue confounds the true expression patterns of its constituent cell types. Here we describe a novel strategy to access such complex samples. Single-cell RNA-seq expression profiles were generated, and clustered to form a two-dimensional cell map onto which expression data were projected. The resulting cell map integrates three levels of organization: the whole population of cells, the functionally distinct subpopulations it contains, and the single cells themselves—all without need for known markers to classify cell types. The feasibility of the strategy was demonstrated by analyzing the transcriptomes of 85 single cells of two distinct types. We believe this strategy will enable the unbiased discovery and analysis of naturally occurring cell types during development, adult physiology, and disease.

Single-cell analysis of gene expression is increasingly important for the analysis of complex tissues, including cancer, developing organs and adult stem cell niches. Here we present a detailed protocol for quantitative gene expression analysis in single cells, by the sequencing of mRNA 5' ends. In all, 96 cells are lysed, and their mRNA is converted to cDNA. By using a template-switching mechanism, a bar code and an upstream primer-binding sequence are introduced simultaneously with reverse transcription. All cDNA is pooled and then prepared for 5' end sequencing, including fragmentation, adapter ligation and PCR amplification. The chief advantage of this approach is the great reduction in cost and time, afforded by the early bar-coding strategy. Compared with previous methods, it is more suitable for large-scale quantitative analysis, as well as for the characterization of transcription start sites, but it is unsuitable for the detection of alternatively spliced transcripts. Sample preparation takes 3 d, and two sets of 96 cells can be prepared in parallel. Finally, the sequencing and data analysis can take an additional 4 d altogether.

Evidence for distinct human cancer stem cells (CSCs) remains contentious and the degree to which different cancer cells contribute to propagating malignancies in patients remains unexplored. In low- to intermediate-risk myelodysplastic syndromes (MDS), we establish the existence of rare multipotent MDS stem cells (MDS-SCs), and their hierarchical relationship to lineage-restricted MDS progenitors. All identified somatically acquired genetic lesions were backtracked to distinct MDS-SCs, establishing their distinct MDS-propagating function in vivo. In isolated del(5q)-MDS, acquisition of del(5q) preceded diverse recurrent driver mutations. Sequential analysis in del(5q)-MDS revealed genetic evolution in MDS-SCs and MDS-progenitors prior to leukemic transformation. These findings provide definitive evidence for rare human MDS-SCs in vivo, with extensive implications for the targeting of the cells required and sufficient for MDS-propagation.

Metastatic breast cancers are still incurable. Characterizing the evolutionary landscape of these cancers, including the role of metastatic axillary lymph nodes (ALNs) in seeding distant organ metastasis, can provide a rational basis for effective treatments. Here, we have described the genomic analyses of the primary tumors and metastatic lesions from 99 samples obtained from 20 patients with breast cancer. Our evolutionary analyses revealed diverse spreading and seeding patterns that govern tumor progression. Although linear evolution to successive metastatic sites was common, parallel evolution from the primary tumor to multiple distant sites was also evident. Metastatic spreading was frequently coupled with polyclonal seeding, in which multiple metastatic subclones originated from the primary tumor and/or other distant metastases. Synchronous ALN metastasis, a well-established prognosticator of breast cancer, was not involved in seeding the distant metastasis, suggesting a hematogenous route for cancer dissemination. Clonal evolution coincided frequently with emerging driver alterations and evolving mutational processes, notably an increase in apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like-associated (APOBEC-associated) mutagenesis. Our data provide genomic evidence for a role of ALN metastasis in seeding distant organ metastasis and elucidate the evolving mutational landscape during cancer progression.

Inhibitory signaling during natural killer (NK) cell education translates into increased responsiveness to activation; however, the intracellular mechanism for functional tuning by inhibitory receptors remains unclear. Secretory lysosomes are part of the acidic lysosomal compartment that mediates intracellular signalling in several cell types. Here we show that educated NK cells expressing self-MHC specific inhibitory killer cell immunoglobulin-like receptors (KIR) accumulate granzyme B in dense-core secretory lysosomes that converge close to the centrosome. This discrete morphological phenotype is independent of transcriptional programs that regulate effector function, metabolism and lysosomal biogenesis. Meanwhile, interference of signaling from acidic Ca2+ stores in primary NK cells reduces target-specific Ca2+-flux, degranulation and cytokine production. Furthermore, inhibition of PI(3,5)P2 synthesis, or genetic silencing of the PI(3,5)P2-regulated lysosomal Ca2+-channel TRPML1, leads to increased granzyme B and enhanced functional potential, thereby mimicking the educated state. These results indicate an intrinsic role for lysosomal remodeling in NK cell education.