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.

PRAM: a novel pooling approach for discovering intergenic transcripts from large-scale RNA sequencing experiments

A Novel Analytical Strategy To Identify Fusion Transcripts Between Repetitive Elements And Protein Coding-Exons Using Rna-Seq

PRAPI: Post-Transcriptional Regulation Analysis Pipeline for Iso-Seq

Enhancing Transcriptome Mapping with Rapid PRO-seq Profiling of Nascent RNA

Genome-Wide Mapping of RNA-Protein Associations via Sequencing

PRADA: pipeline for RNA sequencing data analysis

MPRAnalyze: Statistical Framework for Massively Parallel Reporter Assays

QuaPra: Efficient transcript assembly and quantification using quadratic programming with Apriori algorithm

RADPRE: A Computational Program for Identification of Differential Mrna Processing Including Alternative Polyadenylation

P_RNA_scaffolder: a Fast and Accurate Genome Scaffolder Using Paired-End RNA-sequencing Reads

PrimerPooler: Automated Primer Pooling to Prepare Library for Targeted Sequencing.

Haplotype-aware pantranscriptome analyses using spliced pangenome graphs

Agtar: A Novel Approach For Transcriptome Assembly And Abundance Estimation Using An Adapted Genetic Algorithm From Rna-Seq Data

RAG-seq: A NSR Primed and Transposase Tagmentation Mediated Strand-specific Total RNA Sequencing in Single Cell

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

BinPacker: Packing-Based De Novo Transcriptome Assembly from RNA-seq Data

InPACT: a computational method for accurate characterization of intronic polyadenylation from RNA sequencing data

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

Deep annotation of long noncoding RNAs by assembling RNA-seq and small RNA-seq data

APAtrap: Identification and Quantification of Alternative Polyadenylation Sites from RNA-seq Data.

scAPAtrap: identification and quantification of alternative polyadenylation sites from single-cell RNA-seq data