# Fast and accurate single-cell RNA-seq analysis by clustering of transcript-compatibility counts

- Vasilis Ntranos†
^{1}, - Govinda M. Kamath†
^{2}, - Jesse M. Zhang†
^{2}, - Lior Pachter
^{3}Email author and - David N. Tse
^{2, 1}Email author

**17**:112

**DOI: **10.1186/s13059-016-0970-8

© Ntranos et al. 2016

**Received: **24 February 2016

**Accepted: **29 April 2016

**Published: **26 May 2016

## Abstract

Current approaches to single-cell transcriptomic analysis are computationally intensive and require assay-specific modeling, which limits their scope and generality. We propose a novel method that compares and clusters cells based on their transcript-compatibility read counts rather than on the transcript or gene quantifications used in standard analysis pipelines. In the reanalysis of two landmark yet disparate single-cell RNA-seq datasets, we show that our method is up to two orders of magnitude faster than previous approaches, provides accurate and in some cases improved results, and is directly applicable to data from a wide variety of assays.

## Introduction

Single-cell RNA-seq (scRNA-seq) has proved to be a powerful tool for probing cell states [1–5], defining cell types [6–9], and describing cell lineages [10–13]. These applications of scRNA-seq all rely on two computational steps: quantification of gene or transcript abundances in each cell and clustering of the data in the resulting abundance × cell expression matrix [14, 15]. There are a number of challenges in both of these steps that are specific to scRNA-seq analysis. While methods for transcript/gene abundance estimation from bulk RNA-seq have been extensively tested and benchmarked [16], the wide variety of assay types in scRNA-seq [17–25] have required a plethora of customized solutions [2, 6, 7, 9, 11–13, 24, 26–37] that are difficult to compare to each other. Furthermore, the quantification methods used all rely on read alignment to transcriptomes or genomes, a time-consuming step that will not scale well with the increasing numbers of reads predicted for scRNA-seq [15, 38]. Clustering based on scRNA-seq expression matrices can also require domain-specific information, e.g., temporal information [33] or functional constraints [37], so that in some cases hand curation of clusters is performed after unsupervised clustering [7].

*transcript-compatibility counts*(TCCs), correspond to the sufficient statistics for a standard RNA-seq model [42]. In other words, the use of transcript-compatibility counts was an intermediate computation step towards quantifying transcript abundances. In this paper we instead consider the direct use of such counts for the comparison and clustering of scRNA-seq cells. Figure 2 shows an outline of a method we have developed for clustering and analyzing scRNA-seq data; the key idea is to base clustering not on the quantification of transcripts or genes but on the transcript-compatibility counts for each cell. We note that equivalence classes have also been used in [43, 44] to define similarity scores between de novo assembled transcripts.

To better understand the relevance of transcript-compatibility counts, consider their relationship to the “gene-level” counts used in many RNA-seq analyses. In the same way that “genes” represent groupings of transcripts [45], equivalence classes as introduced by [39] are also groups of transcripts. However, while the former is a biologically motivated construction, the latter is technical, consisting of groupings that capture the extent of ambiguous multiple mappings among reads. The lack of direct biological interpretation of equivalence classes makes transcript-compatibility counts less intuitive. However, as we will show, there are two significant advantages to working with them: (1) unlike transcript or gene-level quantifications, transcript-compatibility counts can be computed without a read-generating model, and hence a single clustering pipeline based on transcript-compatibility counts can be used across a wide range of scRNA-seq assays; (2) transcript-compatibility counts can be computed by pseudoalignment, a process that does not require read alignment and can be done extremely efficiently [41].

To demonstrate both the general applicability of our method as well as its accuracy, we reanalyzed data from the topics of two recently published scRNA-seq papers: the pseudotemporal ordering of primary human myoblasts [12] and the cell classification in the mouse cortex and hippocampus [7]. We show that not only are we able to recapitulate the analyses of the papers two orders of magnitude faster than previously possible, but we also provide a refinement of the published results, suggesting that our approach is both fast and accurate. The speedup of our method makes single-cell RNA-seq analysis interactive for the first time: sensitivity of results to parameters and annotations can be easily explored and analyses can be easily reproduced by individuals without access to significant computing resources. Furthermore, the efficiency of our methods will take on increasing significance as single-cell RNA sequencing scales to experiments with hundreds of thousands of cells and improved technologies make the acquisition of single-cell data easier and faster (see, for example, [46]). We also illustrate the advantages of the broad applicability of our approach via its suitability to a multitude of assays. Existing pipelines must be tailored to specific assays, making it difficult to perform meta-analyses and to compare results across experiments.

## Results

### Transcript-compatibility counts from pseudoalignments

To demonstrate the effectiveness of transcript-compatibility counts for scRNA-seq analysis, we first examined how efficiently they can be computed. While transcript-compatibility counts can be extracted from read alignments (e.g., in SAM/BAM format), they do not require the full information contained in alignments. Instead, we examined the speedup possible with pseudoalignment [41], which obtains for each read the set of transcripts it is compatible with and therefore can be directly used to obtain transcript-compatibility counts.

### Pseudotime for differentiating human myoblasts

A central idea in the pseudo-temporal ordering of cells relies upon the construction of a minimum spanning tree (MST) over the pairwise distances of their corresponding gene expression vectors [48]. This attempts to capture the trajectory of a hypothetical cell that gradually “moves” through different cellular states or differentiation stages in a high-dimensional gene expression space. Our results show that the same concept can be applied to transcript-compatibility counts. A key step in Monocle is to *first* reduce the dimensionality of the data by independent component analysis (ICA) and *then* compute the MST based on Euclidean distances on the plane. Here we take a different approach (see Fig. 4
a–d) and compute the MST on “cluster centers” in high dimensions (see Methods). Both approaches aim to battle the biological and technical noise that is inevitably introduced in scRNA-seq experiments. Even though we could have used Monocle directly on transcript-compatibility counts, the design and comparison of specialized tools is beyond the scope of this paper.

Figure 4
d validates the three primary clusters and the pseudo-temporal ordering obtained by our method based on three key myoblast differentiation markers, *MYOG*, *CDK1*, and *PDGFRA* (see Additional file 1: Figure S3 for an additional set of genes taken from [12]). Interestingly, the expression of these genes gradually evolves over the pseudo-temporally ordered clusters, capturing both the underlying differentiation trajectory of proliferating cells to myoblasts and the corresponding branching towards mesenchymal cells, as was observed in [12].

### Cell classification in the mouse cortex and hippocampus

The reanalysis of [12] shows that clustering of transcript-compatibility counts can be useful on a single dataset, but we believe that the true power of our approach lies in its broad applicability to multiple single-cell assays. In contrast to the standard quantification pipeline, obtaining transcript-compatibility counts does not require a read-generating model; our method can be directly applied to a wide range of scRNA-seq datasets, and transcript-compatibility counts can be used to analyze sequenced reads without any assay-specific information. To make this point, we reanalyzed a recent large scRNA-seq experiment published earlier this year [7] that uses an assay based on unique molecular identifiers (UMIs). In contrast to [12], where paired-end reads were sampled from fragments covering the entire length of the transcripts, [7] used single-end reads that were only obtained from the 3’-end of the transcripts.

Zeisel et al. [7] examined a very diverse population of 3005 cells obtained from the cortical and hippocampal regions of the mouse brain. In order to analyze this complex dataset, the authors developed a state-of-the-art hierarchical biclustering method called BackSPIN (based on SPIN [30]) and were able to identify 47 distinct subpopulations of cells within nine major brain cell types. This fine-grained analysis also revealed a previously unknown post-mitotic oligodendrocyte subclass, referred to as Oligo1 in [7].

Quite remarkably, our method (via affinity propagation on all cells) was further able to recover the Oligo1 cluster of cells, showing that transcript-compatibility counts can indeed capture distinct cell signatures without actually quantifying their gene expression (Fig. 6, Methods). Overall, in our experiments we observed that unsupervised clustering of transcript-compatibility counts typically yielded more than 47 clusters, which was also the case in [7]. Some of our clusters were very small, probably capturing outlier cells, while others seemed to be further splitting the 47 cell subtypes identified in [7].

## Discussion

In this paper we introduced a novel method that uses transcript-compatibility counts — instead of gene expression profiles — as distinct cell signatures for clustering single-cell data. Note, however, that the main focus of our method is not about *how* to cluster (i.e., the particular choice of clustering algorithms) but rather *what* to cluster on. To emphasize this point, we considered simple, “off-the-shelf” clustering methods that directly use the corresponding TCCs as their input. Interestingly, while these methods may not be able to recover accurate clusters when applied to gene expression vectors (see Additional file 1: Figure S4 or ([7], Figure S3), for example), our results showed that TCCs maintain all the necessary information to recover the analyses of [12] and [7].

Even though clustering alone can reveal important information about a single-cell RNA-seq experiment, further biological interpretation of the results (marker gene identification or differential expression) requires some form of quantification of expression profiles within and in between clusters. So it is natural for one to think that eventually the quantification bottleneck will still manifest itself in single-cell analysis. A key observation, however, is that given an accurate clustering of the cells, each and every individual cell’s gene expression profile is no longer needed; one can extract an accurate *statistical* representation of the gene expression within each cluster — without having to quantify all cells separately. In particular, one can quantify the aggregate gene expression in each cluster (cluster centers) by pooling single-cell TCCs together and further estimate the corresponding gene variability by subsampling and quantifying only a few cells per cluster. For example, in Additional file 1: Figure S5 we used kallisto to quantify subsampled cells and the corresponding cluster centers (after clustering on TCCs) for Trapnell et al.’s dataset and generated results that are very similar to the ones obtained in Fig. 4 (where the corresponding gene expression profiles were obtained from [12]). Our method can therefore be used to effectively *reverse* the quantification and clustering steps in the conventional pipeline and potentially provide further end-to-end processing gains, depending on the needs/goals of each scRNA-seq experiment. Overall, we believe that clustering before quantifying is a promising future direction for scRNA-seq analysis which may lead to more robust and accurate quantification algorithms.

## Conclusions

The extraordinary developments in single-cell RNA-seq technology over the past few years have demonstrated that “single-cell resolution” is not just a gimmick but an unprecedented tool for probing transcriptomes that can reveal the inner workings of developmental programs and their resulting tissues. However, the computational challenges of scRNA-seq analysis, already very high due to the large number of cells to analyze, have been further exacerbated by the smorgasbord of assays, each of which introduces unique technical challenges.

The new method we have proposed and evaluated in this paper, namely analysis of scRNA-seq based on transcript-compatibility counts, offers an efficient, accurate, and broadly applicable solution for extracting information from scRNA-seq experiments. In the same way that single-cell analysis can be viewed as the ultimate resolution for transcriptomics, transcript-compatibility counts are the most direct way to “count” reads. While we have focused on clustering of cells in this paper, we believe that transcript-compatibility counts may have applications in many other sequencing-based assays, and that further development of methods based on such counts offers a fruitful avenue of exploration.

## Methods

The code used to generate the results presented in this paper is available online on GitHub [49]. The *Mus musculus* transcriptome assembly used was GRCm38. The *Homo sapiens* transcriptome assembly used was GRCh38. The reference genome used for HISAT was build 10 of the mouse genome (mm10) from the UCSC genome browser. The genome annotation used in the analysis is the release 83 of Ensembl.

### Computation of transcript-compatibility counts

In our implementation of the method, we use kallisto to compute transcript-compatibility counts via pseudoalignment (avoiding the quantification step that is usually performed when running kallisto altogether). In particular, we utilized the “pseudo” option of the kallisto RNA-seq program, which computes equivalence classes of reads after pseudoalignment. We used kallisto version 0.42.3 with *k* set to kallisto’s default value of 31. Even though kallisto pseudoalignment is a natural approach to obtain transcript-compatibility counts, one can in principle extract the same information from exact read alignments. To explore this alternative, we used HISAT (with the no-spliced-alignment option enabled) to align reads on the mouse transcriptome (GRCm38) in the case of Zeisel’s dataset and the human transcriptome (GRCh38) in the case of Trapnell’s dataset. Then, we generated the corresponding “alignment-based TCCs” by directly counting the number of multi-mapped reads aligned to each set of transcripts, and evaluated their performance in Additional file 1: Figure S7.

### Transcript-compatibility counts based on UMI information

The dataset of [7] has reads with unique molecular identifiers (UMIs). UMIs are typically used in scRNA-seq to correct for PCR bias; biological copies of a transcript (distinct molecules) can be identified based on their UMIs. This information can be utilized in generating the transcript-compatibility counts from equivalence classes. Instead of counting all the reads in each equivalence class, we only count the reads with distinct UMIs. Transcript-compatibility counts with UMIs are shown in Figs. 6 b and 8 a (represented as “TCC with UMI” in the figures).

### Clustering methodology

On obtaining the transcript-compatibility counts for each cell, we normalize by the total number of mapped reads to obtain a probability distribution called the transcript-compatibility count distribution or TCC distribution. We then compute the square root of the Jensen-Shannon divergence [50] between the TCC distributions for each pair of cells. As a distance metric which satisfies the triangle [51] inequality, the square root of the Jensen-Shannon divergence is a natural choice for computing pairwise distances between two probability distributions. However, the results obtained here are not contingent on using the square root of Jensen-Shannon divergences as the measure of distances, and quite similar results are obtained when we use other distances between probability distributions, such as the *ℓ*
_{1} distance to compute the pairwise distance matrix (see Additional file 1: Figure S6). The *ℓ*
_{1} distance (which is just twice the total-variation distance) in fact seems to perform better than the Jensen-Shannon distance for low coverage (Additional file 1: Figure S6b). In contrast, Euclidean distance (*ℓ*
_{2} distance) seems to perform much worse (see Additional file 1: Figure S6). The fact that Euclidean distance is not a good distance metric to measure distances between probability distributions is widely documented (see, for instance, [52]).

All clustering carried out in this paper was done using off-the-shelf clustering methods.

We used spectral clustering using the pairwise distance matrices when we knew the number of clusters in the data. This includes Figs. 6 b, 8 a, and Additional file 1: Figure S6b with two clusters for the pairwise distance matrix (from TCC distributions) obtained for the data from [7].

The clustering method used when the number of clusters is not known is affinity propagation [53]. This is an unsupervised clustering algorithm based on message passing, which needs a pairwise similarity matrix as input. The pairwise similarity matrix is computed as the negative of the pairwise distance matrix that was computed.

To evaluate the clustering accuracy of our method in Fig. 6 b, we performed binary classification tests using the labels reported in [7] as the ground truth. In particular, we randomly subsampled two different types of cells and evaluated the ability of each pipeline to separate them into two clusters via spectral clustering. We performed these binary classification tests between (1) the subclasses Oligo1 (45 cells) and Oligo4 (106 cells), (2) the cell types Astrocytes (198 cells) and Interneurons (290 cells), and (3) the more general cell types neurons (1628 cells) and non-neurons (1377 cells). The error rates for each test were obtained by randomly sampling 22, 99, and 200 cells from each of the two labels respectively, averaged over 10 Monte Carlo iterations.

For clustering the dataset of [7], we used affinity propagation with the default parameters which set the preference value equal to the median of the similarity scores and the damping parameter equal to 0.5. On doing this, we obtained 89 clusters. Of the 89 clusters obtained, cluster number 22 had the largest match with the set of cells the authors labeled as Oligo1 (which was the new type of cells discovered in [7]). A total of 24 out of the 28 cells in the cluster were labeled Oligo1 by [7]. There were a total of 45 cells labeled Oligo1 in [7] out of the total of 3005 considered. This is investigated in Fig. 6 c.

Also, affinity propagation with different parameters seems to split the class labeled Oligo3 in [7] into two classes. This is investigated in Fig. 7, where the two classes considered were classes obtained with parameters set as before.

For clustering the dataset of [12], we used affinity propagation with preference parameter set to 1.3 and damping parameter set to 0.95 to obtain three clusters in Fig. 4. To obtain eight clusters on the dataset of [12], we used affinity propagation with preference parameter set to 0.6 and damping parameter set to 0.95 to obtain eight clusters. Then after collapsing any cluster with less than five cells into the cluster closest to it, we obtain the seven clusters investigated in Fig. 5. More details regarding our parameter choices for affinity propagation in this dataset are provided in Additional file 1: Figure S8.

### Partial order on clusters

On the [12] dataset, for the seven clusters obtained, we first find the centroid TCC distribution of each cluster as the mean TCC distribution of all cells in the cluster. Then, we compute the pairwise Jensen-Shannon distances between the centroid TCC distributions (cluster centers). We then run a minimum weight spanning tree on the complete graph between the cluster centers with weights given by the computed pairwise distances. This gives us a partial order on the clusters, which is investigated in Fig. 4 and Additional file 1: Figure S2.

### Quantification

In this paper, we used kallisto and eXpress as representative methods for model-based quantification as they demonstrate similar accuracy to other available quantification tools [41]. In particular, for Zeisel et al.’s dataset (Figs. 6 b and 8 a) we used these tools as a “negative control” to demonstrate the importance of the read-generating model when using quantified transcript abundances to obtain gene expression profiles and cluster single-cell data. Although the significant mismatch between the assumed model (i.e., full transcript length coverage) and the 3’-end bias in this dataset is expected to affect the quantification accuracy for each individual cell, our results show that such a mismatch further impacts the accuracy of clustering and cell-type classification. To evaluate the cell-type classification performance of kallisto and eXpress, we took the minimum of the error rates obtained with bias modeling turned on and off.

For the Trapnell et al. dataset, we used kallisto to quantify transcript abundances and obtain the gene expression profiles within the clusters obtained from TCCs. Note that the read-generating model in this dataset is similar to the standard RNA-seq model that kallisto uses for quantification. More specifically, in Additional file 1: Figure S5a we quantified the corresponding cluster centers by running kallisto’s EM algorithm on the pooled TCCs of each cluster. Using kallisto in this setting resembles bulk RNA-seq quantification applied to the pooled reads coming from each individual cluster (instead of the entire population of cells). In Additional file 1: Figure S5b we further quantified randomly subsampled cells (20 cells per cluster) to obtain an accurate estimate of the gene expression variability within each cluster.

### Visualization of cells and clusters

We used t-SNE [54] to visualize the cells and clusters in Figs. 6 a, 8 b, and Additional file 1: Figure S6a.

The left panel of Figs. 4 d, 5 a, b and Additional file 1: Figure S3a was created using an implementation [55] of the diffusion map algorithm of [56].

## Availability of data and materials

The code used to generate the results presented in this paper is available online on GitHub [49]. All sequencing reads for the Zeisel et al. dataset [7] are available through Gene Expression Omnibus [GEO:GSE60361] and for the Trapnell et al. dataset [12] through [GEO:GSE52529]. The method is publically available on GitHub (https://github.com/govinda-kamath/clustering_on_transcript_compatibility_counts) under the MIT license.

## Ethics

No ethics approval was required for this study.

## Declarations

### Acknowledgements

We thank Páll Melsted for implementing the pseudo command in kallisto. This is the command that allows for direct output of transcript-compatibility counts via pseudoalignment. We would also like to thank Bo Li, Allon Wagner, and Nir Yosef for useful discussions about single-cell RNA-seq assays and their biases.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Patel AP, Tirosh I, Trombetta JJ, Shalek AK, Gillespie SM, Wakimoto H, Cahill DP, Nahed BV, Curry WT, Martuza RL, Louis DN, Rozenblatt-Rosen O, Suvà ML, Regev A, Bernstein BE. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science. 2014; 344(6190):1396–1401. doi:10.1126/science.1254257.View ArticlePubMedPubMed CentralGoogle Scholar
- Pollen AA, Nowakowski TJ, Chen J, Retallack H, Sandoval-Espinosa C, Nicholas CR, Shuga J, Liu SJ, Oldham MC, Diaz A, Lim DA, Leyrat AA, West JA, Kriegstein AR. Molecular identity of human outer radial glia during cortical development. Cell; 163(1):55–67. doi:10.1016/j.cell.2015.09.004.
- Gaublomme JT, Yosef N, Lee Y, Gertner RS, Yang LV, Wu C, Pandolfi PP, Mak T, Satija R, Shalek AK, Kuchroo VK, Park H, Regev A. Single-cell genomics unveils critical regulators of Th17 cell pathogenicity. Cell; 163(6):1400–12. doi:10.1016/j.cell.2015.11.009.
- Kowalczyk MS, Tirosh I, Heckl D, Rao TN, Dixit A, Haas BJ, Schneider RK, Wagers AJ, Ebert BL, Regev A. Single-cell RNA-seq reveals changes in cell cycle and differentiation programs upon aging of hematopoietic stem cells. Genome Res. 2015; 25(12):1860–1872. doi:10.1101/gr.192237.115. http://genome.cshlp.org/content/25/12/1860.full.pdf+html.View ArticlePubMedGoogle Scholar
- Lande-Diner L, Stewart-Ornstein J, Weitz CJ, Lahav G. Single-cell analysis of circadian dynamics in tissue explants. Mol Biol Cell. 2015; 26(22):3940–945. doi:10.1091/mbc.E15-06-0403.View ArticlePubMedPubMed CentralGoogle Scholar
- Usoskin D, Furlan A, Islam S, Abdo H, Lonnerberg P, Lou D, Hjerling-Leffler J, Haeggstrom J, Kharchenko O, Kharchenko PV, Linnarsson S, Ernfors P. Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing. Nat Neurosci. 2015; 18(1):145–53. doi:10.1038/nn.3881.View ArticlePubMedGoogle Scholar
- Zeisel A, Muñoz-Manchado AB, Codeluppi S, Lönnerberg P, La Manno G, Jurèus A, Marques S, Munguba H, He L, Betsholtz C, Rolny C, Castelo-Branco G, Hjerling-Leffler J, Linnarsson S. Cell types in the mouse cortex and hippocampus revealed by single-cell RNA-seq. Science. 2015; 347(6226):1138–42. doi:10.1126/science.aaa1934. http://www.sciencemag.org/content/347/6226/1138.full.pdf.View ArticlePubMedGoogle Scholar
- Burns JC, Kelly MC, Hoa M, Morell RJ, Kelley MW. Single-cell RNA-Seq resolves cellular complexity in sensory organs from the neonatal inner ear. Nat Commun. 2015; 6. doi:10.1038/ncomms9557.
- Grün D, Lyubimova A, Kester L, Wiebrands K, Basak O, Sasaki N, Clevers H, van Oudenaarden A. Single-cell messenger RNA sequencing reveals rare intestinal cell types. Nature. 2015; 525(7568):251–5.View ArticlePubMedGoogle Scholar
- Kafri R, Levy J, Ginzberg MB, Oh S, Lahav G, Kirschner MW. Dynamics extracted from fixed cells reveal feedback linking cell growth to cell cycle. Nature. 2013; 494(7438):480–3. doi:10.1038/nature11897.View ArticlePubMedPubMed CentralGoogle Scholar
- Bendall SC, Davis KL, Amir E-aD, Tadmor MD, Simonds EF, Chen TJ, Shenfeld DK, Nolan GP, Pe’er D. Single-cell trajectory detection uncovers progression and regulatory coordination in human B cell development. Cell. 2014; 157(3):714–25. doi:10.1016/j.cell.2014.04.005.View ArticlePubMedPubMed CentralGoogle Scholar
- Trapnell C, Cacchiarelli D, Grimsby J, Pokharel P, Li S, Morse M, Lennon NJ, Livak KJ, Mikkelsen TS, Rinn JL. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat Biotech. 2014; 32(4):381–6. doi:10.1038/nbt.2859.View ArticleGoogle Scholar
- Buettner F, Natarajan KN, Casale FP, Proserpio V, Scialdone A, Theis FJ, Teichmann SA, Marioni JC, Stegle O. Computational analysis of cell-to-cell heterogeneity in single-cell RNA-sequencing data reveals hidden subpopulations of cells. Nat Biotech. 2015; 33(2):155–60. doi:10.1038/nbt.3102.View ArticleGoogle Scholar
- Shapiro E, Biezuner T, Linnarsson S. Single-cell sequencing-based technologies will revolutionize whole-organism science. Nat Rev Genet. 2013; 14(9):618–30. doi:10.1038/nrg3542.View ArticlePubMedGoogle Scholar
- Stegle O, Teichmann SA, Marioni JC. Computational and analytical challenges in single-cell transcriptomics. Nat Rev Genet. 2015; 16(3):133–45. doi:10.1038/nrg3833.View ArticlePubMedGoogle Scholar
- Oshlack A, Robinson M, Young M. From RNA-seq reads to differential expression results. Genome Biol. 2010; 11(12):220. doi:10.1186/gb-2010-11-12-220.View ArticlePubMedPubMed CentralGoogle Scholar
- Tang F, Barbacioru C, Wang Y, Nordman E, Lee C, Xu N, Wang X, Bodeau J, Tuch BB, Siddiqui A, Lao K, Surani MA. mRNA-Seq whole-transcriptome analysis of a single cell. Nat Meth. 2009; 6(5):377–82. doi:10.1038/nmeth.1315.View ArticleGoogle Scholar
- Islam S, Kjällquist U, Moliner A, Zajac P, Fan JB, Lönnerberg P, Linnarsson S. Characterization of the single-cell transcriptional landscape by highly multiplex RNA-seq. Genome Res. 2011; 21(7):1160–1167. doi:10.1101/gr.110882.110.View ArticlePubMedPubMed CentralGoogle Scholar
- Ramsköld D, Luo S, Wang YC, Li R, Deng Q, Faridani OR, Daniels GA, Khrebtukova I, Loring JF, Laurent LC, Schroth GP, Sandberg R. Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells. Nat Biotech. 2012; 30(8):777–82. doi:10.1038/nbt.2282.View ArticleGoogle Scholar
- Hashimshony T, Wagner F, Sher N, Yanai I. CEL-Seq: single-cell RNA-Seq by multiplexed linear amplification. Cell Rep. 2012; 2(3):666–73. doi:10.1016/j.celrep.2012.08.003.View ArticlePubMedGoogle Scholar
- Picelli S, Bjorklund AK, Faridani OR, Sagasser S, Winberg G, Sandberg R. Smart-seq2 for sensitive full-length transcriptome profiling in single cells. Nat Meth. 2013; 10(11):1096–1098. doi:10.1038/nmeth.2639.View ArticleGoogle Scholar
- Sasagawa Y, Nikaido I, Hayashi T, Danno H, Uno K, Imai T, Ueda H. Quartz-Seq: a highly reproducible and sensitive single-cell RNA sequencing method, reveals non-genetic gene-expression heterogeneity. Genome Biol. 2013; 14(4):31. doi:10.1186/gb-2013-14-4-r31.View ArticleGoogle Scholar
- Jaitin DA, Kenigsberg E, Keren-Shaul H, Elefant N, Paul F, Zaretsky I, Mildner A, Cohen N, Jung S, Tanay A, Amit I. Massively parallel single-cell RNA-Seq for marker-free decomposition of tissues into cell types. Science. 2014; 343(6172):776–9. doi:10.1126/science.1247651.View ArticlePubMedPubMed CentralGoogle Scholar
- Macosko EZ, Basu A, Satija R, Nemesh J, Shekhar K, Goldman M, Tirosh I, Bialas AR, Kamitaki N, Martersteck EM, et al. Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell. 2015; 161(5):1202–1214. doi:10.1016/j.cell.2015.05.002.View ArticlePubMedPubMed CentralGoogle Scholar
- Klein AM, Mazutis L, Akartuna I, Tallapragada N, Veres A, Li V, Peshkin L, Weitz DA, Kirschner MW. Droplet barcoding for single-cell transcriptomics applied to embryonic stem cells. Cell; 161(5):1187–201. doi:10.1016/j.cell.2015.04.044.
- Amir E-aD, Davis KL, Tadmor MD, Simonds EF, Levine JH, Bendall SC, Shenfeld DK, Krishnaswamy S, Nolan GP, Pe’er D. viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia. Nat Biotech. 2013; 31(6):545–52. doi:10.1038/nbt.2594.View ArticleGoogle Scholar
- Mahfouz A, van de Giessen M, van der Maaten L, Huisman S, Reinders M, Hawrylycz MJ, Lelieveldt BPF. Visualizing the spatial gene expression organization in the brain through non-linear similarity embeddings. Methods. 2015; 73:79–89. doi:10.1016/j.ymeth.2014.10.004.View ArticlePubMedGoogle Scholar
- Shekhar K, Brodin P, Davis MM, Chakraborty AK. Automatic Classification of Cellular Expression by Nonlinear Stochastic Embedding (ACCENSE). Proc Natl Acad Sci U S A. 2014; 111(1):202–7. doi:10.1073/pnas.1321405111.View ArticlePubMedPubMed CentralGoogle Scholar
- Islam S, Zeisel A, Joost S, La Manno G, Zajac P, Kasper M, Lonnerberg P, Linnarsson S. Quantitative single-cell RNA-seq with unique molecular identifiers. Nat Meth. 2014; 11(2):163–6. doi:10.1038/nmeth.2772.View ArticleGoogle Scholar
- Tsafrir D, Tsafrir I, Ein-Dor L, Zuk O, Notterman DA, Domany E. Sorting points into neighborhoods (SPIN): data analysis and visualization by ordering distance matrices. Bioinformatics. 2005; 21(10):2301–308. doi:10.1093/bioinformatics.View ArticlePubMedGoogle Scholar
- Qiu P, Simonds EF, Bendall SC, Gibbs Jr KD, Bruggner RV, Linderman MD, Sachs K, Nolan GP, Plevritis SK. Extracting a cellular hierarchy from high-dimensional cytometry data with SPADE. Nat Biotech. 2011; 29(10):886–91. doi:10.1038/nbt.1991.View ArticleGoogle Scholar
- Levine JH, Simonds EF, Bendall SC, Davis KL, Amir E-aD, Tadmor MD, Litvin O, Fienberg HG, Jager A, Zunder ER, Finck R, Gedman AL, Radtke I, Downing JR, Pe’er D, Nolan GP. Data-driven phenotypic dissection of AML reveals progenitor-like cells that correlate with prognosis. Cell. 2015; 162(1):184–97. doi:10.1016/j.cell.2015.05.047.View ArticlePubMedGoogle Scholar
- Marco E, Karp RL, Guo G, Robson P, Hart AH, Trippa L, Yuan GC. Bifurcation analysis of single-cell gene expression data reveals epigenetic landscape. Proc Natl Acad Sci. 2014; 111(52):5643–650. doi:10.1073/pnas.1408993111.View ArticleGoogle Scholar
- Shin J, Berg DA, Zhu Y, Shin JY, Song J, Bonaguidi MA, Enikolopov G, Nauen DW, Christian KM, Ming G-l, Song H. Single-cell RNA-seq with Waterfall reveals molecular cascades underlying adult Neurogenesis. Cell Stem Cell; 17(3):360–72. doi:10.1016/j.stem.2015.07.013.
- Haghverdi L, Buettner F, Theis FJ. Diffusion maps for high-dimensional single-cell analysis of differentiation data. Bioinformatics. 2015. doi:10.1093/bioinformatics.
- Moignard V, Woodhouse S, Haghverdi L, Lilly AJ, Tanaka Y, Wilkinson AC, Buettner F, Macaulay IC, Jawaid W, Diamanti E, Nishikawa SI, Piterman N, Kouskoff V, Theis FJ, Fisher J, Gottgens B. Decoding the regulatory network of early blood development from single-cell gene expression measurements. Nat Biotech. 2015; 33(3):269–76. doi:10.1038/nbt.3154.View ArticleGoogle Scholar
- Fan J, Salathia N, Liu R, Kaeser G, Yung Y, Herman JL, Kaper F, Fan JB, Zhang K, Chun J, Kharchenko P. Characterizing transcriptional heterogeneity through pathway and gene set overdispersion analysis. bioRxiv. 2015. doi:10.1101/026948.
- Saliba AE, Westermann AJ, Gorski SA, Vogel J. Single-cell RNA-seq: advances and future challenges. Nucleic Acids Res. 2014. doi:10.1093/nar.
- Nicolae M, Mangul S, Mandoiu II, Zelikovsky A. Estimation of alternative splicing isoform frequencies from RNA-seq data. Algorithm Mol Biol. 2011; 6(1):9.View ArticleGoogle Scholar
- Patro R, Mount SM, Kingsford C. Sailfish enables alignment-free isoform quantification from RNA-seq reads using lightweight algorithms. Nat Biotechnol. 2014; 32(5):462–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat Biotech. 2016; advance online publication. doi:10.1038/nbt.3519.
- Pachter L. Models for transcript quantification from RNA-Seq. arXiv preprint arXiv:1104.3889. 2011. https://arxiv.org/pdf/1104.3889v2.pdf.
- Davidson NM, Oshlack A. Corset: enabling differential gene expression analysis for de novo assembled transcriptomes. Genome Biol. 2014; 15(7):410.PubMedPubMed CentralGoogle Scholar
- Srivastava A, Sarkar H, Patro R. RapMap: a rapid, sensitive and accurate tool for mapping RNA-seq reads to transcriptomes. bioRxiv. 2015. 029652. http://www.biorxiv.org/content/biorxiv/early/2015/10/22/029652.full.pdf.
- Gerstein MB, Bruce C, Rozowsky JS, Zheng D, Du J, Korbel JO, Emanuelsson O, Zhang ZD, Weissman S, Snyder M. What is a gene, post-ENCODE? History and updated definition. Genome Res. 2007; 17(6):669–81. http://genome.cshlp.org/content/17/6/669.abstract.View ArticlePubMedGoogle Scholar
- 10x Genomics to unveil new single-cell genetic analysis product. 2016. http://www.10xgenomics.com/news/10x-genomics-to-unveil-new-single-cell-genetic-analysis-product. Accessed 16 Feb 2016.
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, Van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010; 28(5):511–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Magwene PM, Lizardi P, Kim J. Reconstructing the temporal ordering of biological samples using microarray data. Bioinformatics. 2003; 19(7):842–50. doi:10.1093/bioinformatics.View ArticlePubMedGoogle Scholar
- Clustering on transcript compatibility counts. 2016. https://github.com/govinda-kamath/clustering_on_transcript_compatibility_counts. Github repository (2016). Accessed 15 May 2016.
- Lin J. Divergence measures based on the Shannon entropy. IEEE Trans Inf Theory. 1991; 37(1):145–51.View ArticleGoogle Scholar
- Fuglede B, Topsoe F. Jensen-Shannon divergence and Hilbert space embedding. In: IEEE International Symposium on Information Theory. Chicago: ISIT: 2004. p. 31–1.Google Scholar
- Batu T, Fortnow L, Rubinfeld R, Smith WD, White P. Testing that distributions are close. In: Proceedings 41st Annual Symposium on Foundations of Computer Science. Redondo Beach: IEEE: 2000. p. 259–69.Google Scholar
- Frey BJ, Dueck D. Clustering by passing messages between data points. science. 2007; 315(5814):972–6.View ArticlePubMedGoogle Scholar
- Van der Maaten L, Hinton G. Visualizing data using t-SNE. J Mach Learn Res. 2008; 9(2579–2605):85.Google Scholar
- Mühlbacher P. A python implementation of the diffusion maps algorithm introduced by Lafon. GitHub. 2015.
- Lafon SS. Diffusion maps and geometric harmonics. Ph.D. thesis: Yale University; 2004.