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.

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

TransLiG: a de novo transcriptome assembler that uses line graph iteration

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

TransBorrow: Genome-Guided Transcriptome Assembly by Borrowing Assemblies from Different Assemblers.

MultiTrans: an Algorithm for Path Extraction Through Mixed Integer Linear Programming for Transcriptome Assembly

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

TransComb: Genome-Guided Transcriptome Assembly Via Combing Junctions in Splicing Graphs

Scallop2 Enables Accurate Assembly of Multiple-End RNA-seq Data

A novel min-cost flow method for estimating transcript expression with RNA-Seq

IsoTree: De Novo Transcriptome Assembly from RNA-Seq Reads

TransMeta Simultaneously Assembles Multisample RNA-seq Reads

Piecing the Puzzle Together: a Revisit to Transcript Reconstruction Problem in RNA-seq

StringTie enables improved reconstruction of a transcriptome from RNA-seq reads

TransCoord: Genome-guided Transcripts Assembly by Coordinating Candidate Paths into Two-phased Linear Programming

Bridger: A New Framework for De Novo Transcriptome Assembly Using Rna-Seq Data

IsoTree: A New Framework for De Novo Transcriptome Assembly from RNA-seq Reads

Scrnass: a Single-Cell RNA-seq Assembler Via Imputing Dropouts and Combing Junctions

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

Improved Transcript Isoform Discovery Using ORF Graphs.

A Novel Combinatorial Method for Estimating Transcript Expression with RNA-Seq: Bounding the Number of Paths

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