Compacting and correcting Trinity and Oases RNA-Seq de novo assemblies

PeerJ. 2017 Feb 16:5:e2988. doi: 10.7717/peerj.2988. eCollection 2017.

Abstract

Background: De novo transcriptome assembly of short reads is now a common step in expression analysis of organisms lacking a reference genome sequence. Several software packages are available to perform this task. Even if their results are of good quality it is still possible to improve them in several ways including redundancy reduction or error correction. Trinity and Oases are two commonly used de novo transcriptome assemblers. The contig sets they produce are of good quality. Still, their compaction (number of contigs needed to represent the transcriptome) and their quality (chimera and nucleotide error rates) can be improved.

Results: We built a de novo RNA-Seq Assembly Pipeline (DRAP) which wraps these two assemblers (Trinity and Oases) in order to improve their results regarding the above-mentioned criteria. DRAP reduces from 1.3 to 15 fold the number of resulting contigs of the assemblies depending on the read set and the assembler used. This article presents seven assembly comparisons showing in some cases drastic improvements when using DRAP. DRAP does not significantly impair assembly quality metrics such are read realignment rate or protein reconstruction counts.

Conclusion: Transcriptome assembly is a challenging computational task even if good solutions are already available to end-users, these solutions can still be improved while conserving the overall representation and quality of the assembly. The de novo RNA-Seq Assembly Pipeline (DRAP) is an easy to use software package to produce compact and corrected transcript set. DRAP is free, open-source and available under GPL V3 license at http://www.sigenae.org/drap.

Keywords: Compaction; Correction; De novo assembly; Quality assessment; RNA-Seq.

Grants and funding

This work was supported by funds from the ANR (http://www.agence-nationale-recherche.fr/) in the frame of the PhyloFish (ANR-10-GENM-017) and the France Génomique (ANR10-INBS-09-08) projects. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.