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.

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

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

Self-supervised deep learning of gene–gene interactions for improved gene expression recovery

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

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

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

Distribution‐Agnostic Deep Learning Enables Accurate Single‐Cell Data Recovery and Transcriptional Regulation Interpretation

SAVER: gene expression recovery for single-cell RNA sequencing

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

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

VIPER: variability-preserving imputation for accurate gene expression recovery in single-cell RNA sequencing studies

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

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

High-throughput Single-Cell RNA-seq Data Imputation and Characterization with Surrogate-Assisted Automated Deep Learning

Missing Value Imputation With Low-Rank Matrix Completion in Single-Cell RNA-Seq Data by Considering Cell Heterogeneity

Zero-preserving Imputation of Scrna-Seq Data Using Low-Rank Approximation

Cellular Similarity based Imputation for Single cell RNA Sequencing Data

A Unified Statistical Framework for Single Cell and Bulk RNA Sequencing Data

Accurate and interpretable gene expression imputation on scRNA-seq data using IGSimpute

Data Smoothing Filling Method based on ScRNA-Seq Data Zero-Value Identification

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