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.

Mapping short reads to a genome without using hash look-up table algorithm and Burrows Wheeler Transformation

CUDA-FRESCO: Frequency-Based RE-Sequencing Tool Based on CO-clustering Segmentation by GPU

A Re-Sequencing Tool For High Mismatch-Tolerant Short Read Alignment Based On Burrows-Wheeler Transform

An efficient Burrows-Wheeler transform-based aligner for short read mapping

ReadsMap: a new tool for high precision mapping of DNAseq and RNAseq read sequences

Minimap and miniasm: fast mapping and de novo assembly for noisy long sequences

Perm: Efficient Mapping of Short Sequencing Reads with Periodic Full Sensitive Spaced Seeds

FANSe: an accurate algorithm for quantitative mapping of large scale sequencing reads.

FANSe2: a robust and cost-efficient alignment tool for quantitative next-generation sequencing applications.

Cloud Based Short Read Mapping Service

Short read mapping for exome sequencing.

The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote

A mapping-free NLP-based technique for sequence search in Nanopore long-reads

Fast and Accurate Read Alignment for Resequencing.

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

featureCounts: an efficient general purpose program for assigning sequence reads to genomic features

GRIM-Filter: Fast Seed Location Filtering in DNA Read Mapping Using Processing-in-Memory Technologies

Seed-and-Vote based In-Memory Accelerator for DNA Read Mapping

Accelerating massive short reads mapping for next generation sequencing (abstract only).

Fast and efficient Rmap assembly using the Bi-labelled de Bruijn graph

RNA-Skim: a Rapid Method for RNA-Seq Quantification at Transcript Level.