Abstract:Abstract Droplet-based 3’ single-cell RNA-sequencing (scRNA-seq) methods have proved transformational in characterizing cellular diversity and generating valuable hypotheses throughout biology 1,2 . Here we outline a common problem with 3’ scRNA-seq datasets where genes that have been documented to be expressed with other methods, are either completely missing or are dramatically under-represented thereby compromising the discovery of cell types, states, and genetic mechanisms. We show that this problem stems from three main sources of sequencing read loss: (1) reads mapping immediately 3’ to known gene boundaries due to poor 3’ UTR annotation; (2) intronic reads stemming from unannotated exons or pre-mRNA; (3) discarded reads due to gene overlaps 3 . Each of these issues impacts the detection of thousands of genes even in well-characterized mouse and human genomes rendering downstream analysis either partially or fully blind to their expression. We outline a simple three-step solution to recover the missing gene expression data that entails compiling a hybrid pre-mRNA reference to retrieve intronic reads 4 , resolving gene collision derived read loss through removal of readthrough and premature start transcripts, and redefining 3’ gene boundaries to capture false intergenic reads. We demonstrate with mouse brain and human peripheral blood datasets that this approach dramatically increases the amount of sequencing data included in downstream analysis revealing 20 - 50% more genes per cell and incorporates 15-20% more sequencing reads than with standard solutions 5 . These improvements reveal previously missing biologically relevant cell types, states, and marker genes in the mouse brain and human blood profiling data. Finally, we provide scRNA-seq optimized transcriptomic references for human and mouse data as well as simple algorithmic implementation of these solutions that can be deployed to both thoroughly as well as poorly annotated genomes. Our results demonstrate that optimizing the sequencing read mapping step can significantly improve the analysis resolution as well as biological insight from scRNA-seq. Moreover, this approach warrants a fresh look at preceding analyses of this popular and scalable cellular profiling technology.

SAVER: gene expression recovery for single-cell RNA sequencing

De Novo Identification of Expressed Cancer Somatic Mutations from Single-Cell RNA Sequencing Data

Enhanced recovery of single-cell RNA-sequencing reads for missing gene expression data

A new bioinformatics tool to recover missing gene expression in single-cell RNA sequencing data

Recovery of missing single-cell RNA-sequencing data with optimized transcriptomic references

scTSSR: gene expression recovery for single-cell RNA sequencing using two-side sparse self-representation

Identification of cell types, states and programs by learning gene set representations

Issues arising from benchmarking single-cell RNA sequencing imputation methods

sc-SPLASH provides ultra-efficient reference-free discovery in barcoded single-cell sequencing

scRNASequest: an ecosystem of scRNA-seq analysis, visualization, and publishing

Leveraging data-driven self-consistency for high-fidelity gene expression recovery

SMURF: Embedding Single-Cell RNA-seq Data with Matrix Factorization Preserving Self-Consistency

scRecover: Discriminating true and false zeros in single-cell RNA-seq data for imputation

V-SVA: an R Shiny application for detecting and annotating hidden sources of variation in single-cell RNA-seq data

A Highly Sensitive, Accurate, and Automated Single-Cell RNA Sequencing Platform with Digital Microfluidics.

Deep Generative Autoencoder for Low-Dimensional Embeding Extraction from Single-Cell RNAseq Data

Recovering single-cell expression profiles from spatial transcriptomics with scResolve

Reply to "Issues arising from benchmarking single-cell RNA sequencing imputation methods"

Sctssr-D: Gene Expression Recovery by Two-side Self-Representation and Dropout Information for Scrna-Seq Data

scIALM: A method for sparse scRNA-seq expression matrix imputation using the Inexact Augmented Lagrange Multiplier with low error

Quantification of transcript isoforms at the single-cell level using SCALPEL