Abstract:Motivation: RNA-seq techniques provide an unparalleled means for exploring a transcriptome with deep coverage and base pair level resolution. Various analysis tools have been developed to align and assemble RNA-seq data, such as the widely used TopHat/Cufflinks pipeline. A common observation is that a sizable fraction of the fragments/reads align to multiple locations of the genome. These multiple alignments pose substantial challenges to existing RNA-seq analysis tools. Inappropriate treatment may result in reporting spurious expressed genes (false positives) and missing the real expressed genes (false negatives). Such errors impact the subsequent analysis, such as differential expression analysis. In our study, we observe that similar to 3.5% of transcripts reported by TopHat/Cufflinks pipeline correspond to annotated nonfunctional pseudogenes. Moreover, similar to 10.0% of reported transcripts are not annotated in the Ensembl database. These genes could be either novel expressed genes or false discoveries. Results: We examine the underlying genomic features that lead to multiple alignments and investigate how they generate systematic errors in RNA-seq analysis. We develop a general tool, GeneScissors, which exploits machine learning techniques guided by biological knowledge to detect and correct spurious transcriptome inference by existing RNA-seq analysis methods. In our simulated study, GeneScissors can predict spurious transcriptome calls owing to misalignment with an accuracy close to 90%. It provides substantial improvement over the widely used TopHat/Cufflinks or MapSplice/Cufflinks pipelines in both precision and F-measurement. On real data, GeneScissors reports 53.6% less pseudogenes and 0.97% more expressed and annotated transcripts, when compared with the TopHat/Cufflinks pipeline. In addition, among the 10.0% unannotated transcripts reported by TopHat/Cufflinks, GeneScissors finds that > 16.3% of them are false positives.

RES-Scanner: a Software Package for Genome-Wide Identification of RNA-editing Sites

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

RED: A Java-MySQL Software for Identifying and Visualizing RNA Editing Sites Using Rule-Based and Statistical Filters

SPRINT: an SNP-free Toolkit for Identifying RNA Editing Sites.

Full-length PacBio Amplicon Sequencing to Unveil RNA Editing Sites

Sequence based identification of RNA editing sites

Dbres: A Web-Oriented Database for Annotated Rna Editing Sites

RED-ML: a Novel, Effective RNA Editing Detection Method Based on Machine Learning

Identifying RNA Editing Sites Using RNA Sequencing Data Alone

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

Elucidating the editome: bioinformatics approaches for RNA editing detection

Genome sequence–independent identification of RNA editing sites

REIA: A Database for Cancer A-to-I RNA Editing with Interactive Analysis.

MultiEditR: An easy validation method for detecting and quantifying RNA editing from Sanger sequencing

Snpdetector: A Software Tool for Sensitive and Accurate Snp Detection

CRISPR-detector: fast and accurate detection, visualization, and annotation of genome-wide mutations induced by genome editing events

RealSeq2: a Software Integrated with UMI Identification, Error Correction, and Methylation Modifications Storing

RNASEQR—a Streamlined and Accurate RNA-seq Sequence Analysis Program

RADAR: a rigorously annotated database of A-to-I RNA editing

BS-RNA: An efficient mapping and annotation tool for RNA bisulfite sequencing data

PRES: a webserver for decoding the functional perturbations of RNA editing sites