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.

RNA-Seq Mapping Errors When Using Incomplete Reference Transcriptomes of Vertebrates

Improving transcriptome assembly through error correction of high-throughput sequence reads

GeneScissors: a comprehensive approach to detecting and correcting spurious transcriptome inference owing to RNA-seq reads misalignment.

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

Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences

Enhancing transcriptome expression quantification through accurate assignment of long RNA sequencing reads with TranSigner

Effects of transcriptional noise on estimates of gene and transcript expression in RNA sequencing experiments

Estimating and correcting index hopping misassignments in single-cell RNA-seq data

Comparative analysis of RNA-Seq alignment algorithms and the RNA-Seq unified mapper (RUM)

Characterizing and Annotating the Genome Using RNA-seq Data

Importance of Transcript Variants in Transcriptome Analyses

Knowledge-Based Reconstruction of Mrna Transcripts with Short Sequencing Reads for Transcriptome Research

Illuminating the dark side of the human transcriptome with long read transcript sequencing

Revealing the missing expressed genes beyond the human reference genome by RNA-Seq

Multimapping confounds ribosome profiling analysis: A case‐study of the Hsp90 molecular chaperone

Modeling non-uniformity in short-read rates in RNA-Seq data

High-Resolution Transcriptome Analysis with Long-Read RNA Sequencing

TransRef Enables Accurate Transcriptome Assembly by Redefining Accurate Neo-Splicing Graphs.

RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome

The Difficult Calls in RNA Editing