Characterizing the transcriptome of individual cells is fundamental to understanding complex biological systems. We describe a droplet-based system that enables 3′ mRNA counting of tens of thousands of single cells per sample. Cell encapsulation, of up to 8 samples at a time, takes place in ∼6 min, with ∼50% cell capture efficiency. To demonstrate the system's technical performance, we collected transcriptome data from ∼250k single cells across 29 samples. We validated the sensitivity of the system and its ability to detect rare populations using cell lines and synthetic RNAs. We profiled 68k peripheral blood mononuclear cells to demonstrate the system's ability to characterize large immune populations. Finally, we used sequence variation in the transcriptome data to determine host and donor chimerism at single-cell resolution from bone marrow mononuclear cells isolated from transplant patients.

40Characterizing the transcriptome of individual cells is fundamental to understanding complex 41 biological systems. We describe a droplet-based system that enables 3' mRNA counting of up 56peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.The copyright holder for this preprint (which was not . http://dx.doi.org/10.1101/065912 doi: bioRxiv preprint first posted online 84 RESULTS 86Droplet-based platform enables barcoding of tens of thousands of cells 88The scRNA-seq microfluidics platform builds on the GemCode ® technology, which has 89 been used for genome haplotyping, structural variant analysis and de novo assembly of a 90human genome [10][11][12] . The core of the technology is a Gel bead in Emulsion (GEM). GEM 91 peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.The copyright holder for this preprint (which was not . http://dx.doi.org/10.1101/065912 doi: bioRxiv preprint first posted online 5 generation takes place in an 8-channel microfluidic chip that encapsulates single gel beads at ~80% fill rate (Fig. 1a-c). Each gel bead is functionalized with barcoded oligonucleotides that 93 consist of: i) sequencing adapters and primers, ii) a 14bp barcode drawn from ~750,000 94 designed sequences to index GEMs, iii) a 10bp randomer to index molecules (unique molecular 95 identifier, UMI), and iv) an anchored 30bp oligo-dT to prime poly-adenylated RNA transcripts 96 (Fig. 1d). Within each microfluidic channel, ~100,000 GEMs are formed per ~6-min run, 97encapsulating thousands of cells in GEMs. Cells are loaded at a limiting dilution to minimize co- 98occurrence of multiple cells in the same GEM. 100Cell lysis begins immediately after encapsulation. Gel beads automatically dissolve to 101 release their oligonucleotides for reverse transcription of poly-adenylated RNAs. Each cDNA 102 molecule contains a UMI and shared barcode per GEM, and ends with a template switching 103 oligo at the 3' end (Fig. 1e). Next, the emulsion is broken and barcoded cDNA is pooled for 104PCR amplification, using primers complementary to the switch oligos and sequencing adapters. Methods, Fig. 1f). Briefly, 98-nt of Read1s were aligned against the union of human (hg19) 123peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.The copyright holder for this preprint (which was not . http://dx.doi.org/10.1101/065912 doi: bioRxiv preprint first posted online 6 Based on the distribution of total UMI counts for each barcode (Online Methods), we 124 estimated that 1,012 GEMs contained cells, of which 482 and 538 contained reads that mapped 125 primarily to the human and mouse transcriptome, respectively (and will be referred to as human 126 and mouse GEMs) (Fig. 2a). >83% of UMI counts were associated with cell barcodes, 127indicating low background of cell-free RNA. Eight cell-containing GEMs had a substantial 128 fraction of human and mouse UMI counts (the UMI count is >1% of each species' UMI...

No abstract

SUMMARY The canonical Wnt/β-catenin signaling pathway governs diverse developmental, homeostatic and pathologic processes. Palmitoylated Wnt ligands engage cell surface Frizzled (Fzd) receptors and Lrp5/6 co-receptors enabling β-catenin nuclear translocation and Tcf/Lef-dependent gene transactivation1–3. Mutations in Wnt downstream signaling components have revealed diverse functions presumptively attributed to Wnt ligands themselves, although direct attribution remains elusive, as complicated by redundancy between 19 mammalian Wnts and 10 Fzds1 and Wnt hydrophobicity2,3. For example, individual Wnt ligand mutations have not revealed homeostatic phenotypes in the intestinal epithelium4, an archetypal canonical Wnt pathway-dependent rapidly self-renewing tissue whose regeneration is fueled by proliferative crypt Lgr5+ intestinal stem cells (ISCs)5–9. R-spondin ligands (Rspo1–4) engage distinct Lgr4-6 and Rnf43/Znrf3 receptor classes10–13, markedly potentiate canonical Wnt/β-catenin signaling and induce intestinal organoid growth in vitro and Lgr5+ ISCs in vivo8,14–17. However, the interchangeability, functional cooperation and relative contributions of Wnt versus Rspo ligands to in vivo canonical Wnt signaling and ISC biology remain unknown. Here, we deconstructed functional roles of Wnt versus Rspo ligands in the intestinal crypt stem cell niche. We demonstrate that the default fate of Lgr5+ ISCs is lineage commitment, escape from which requires both Rspo and Wnt ligands. However, gain-of-function studies using Rspo versus a novel non-lipidated Wnt analog reveal qualitatively distinct, non-interchangeable roles for these ligands in ISCs. Wnts are insufficient to induce Lgr5+ ISC self-renewal, but rather confer a basal competency by maintaining Rspo receptor expression that enables Rspo to actively drive and specify the extent of stem cell expansion. This functionally non-equivalent yet cooperative interplay between Wnt and Rspo ligands establishes a molecular precedent for regulation of mammalian stem cells by distinct priming and self-renewal factors, with broad implications for precision control of tissue regeneration.

SUMMARY Several cell populations have been reported to possess intestinal stem cell (ISC) activity during homeostasis and injury-induced regeneration. Here, we explored inter-relationships between putative mouse ISC populations by comparative RNA-sequencing (RNA-seq). The transcriptomes of multiple cycling ISC populations closely resembled Lgr5+ ISCs, the most well-defined ISC pool, but Bmi1-GFP+ cells were distinct and enriched for enteroendocrine (EE) markers including Prox1. Prox1-GFP+ cells exhibited sustained clonogenic growth in vitro, and lineage-tracing of Prox1+ cells revealed long-lived clones during homeostasis and after radiation-induced injury in vivo. Single-cell mRNA-seq revealed two subsets of Prox1-GFP+ cells, one of which resembled mature EE cells while the other displayed low level EE gene expression but co-expressed tuft cell markers, Lgr5 and Ascl2, reminiscent of label-retaining secretory progenitors. Our data suggest that the EE lineage, including mature EE cells, comprise a reservoir of homeostatic and injury-inducible ISCs, extending our understanding of cellular plasticity and stemness.