Comprehensive benchmarking of single cell RNA sequencing technologies for characterizing cellular perturbation

Technological advances in transcriptome sequencing of single cells continues to provide an unprecedented view on tissue composition and cellular heterogeneity. While several studies have compared different single cell RNA-seq methods with respect to data quality and their ability to distinguish cell subpopulations, none of these studies investigated the heterogeneity of the cellular transcriptional response upon a chemical perturbation. In this study, we evaluated the transcriptional response of NGP neuroblastoma cells upon nutlin-3 treatment using the C1, ddSeq and Chromium single cell systems. These devices and library preparation methods are representative for the wide variety of platforms, ranging from microfluid chips to droplet-based systems and from full transcript sequencing to 3-prime end sequencing. In parallel, we used bulk RNA-seq for molecular characterization of the transcriptional response. Two complementary metrics to evaluate performance were applied: the first is the number and identity of differentially expressed genes as defined in consensus by two statistical models, and the second is the enrichment analysis of biological signals. Where relevant, to make the data more comparable, we downsampled sequencing library size, selected cell subpopulations based on specific RNA abundance features, or created pseudobulk samples. While the C1 detects the highest number of genes per cell and better resembles bulk RNA-seq, the Chromium identifies most differentially expressed genes, albeit still substantially fewer than bulk RNA-seq. Gene set enrichment analyses reveals that detection of a limited set of the most abundant genes in single cell RNA-seq experiments is sufficient for molecular phenotyping. Finally, single cell RNA-seq reveals a heterogeneous response of NGP neuroblastoma cells upon nutlin-3 treatment, revealing putative late-responder or resistant cells, both undetected in bulk RNA-seq experiments.