Evaluating Characteristics of De Novo Assembly Software on 454 Transcriptome Data: A Simulation Approach

Background: The quantity of transcriptome data is rapidly increasing for non-model organisms. As sequencing technology advances, focus shifts towards solving bioinformatic challenges, of which sequence read assembly is the first task. Recent studies have compared the performance of different softwar...

Verfasser: Mundry, Marvin
Bornberg-Bauer, Erich
Sammeth, Michael
Feulner, Philine G. D.
FB/Einrichtung:FB 13: Biologie
Dokumenttypen:Artikel
Medientypen:Text
Erscheinungsdatum:2012
Publikation in MIAMI:21.02.2013
Datum der letzten Änderung:26.10.2023
Angaben zur Ausgabe:[Electronic ed.]
Quelle:PLoS ONE 7 (2012) 2, e31410
Fachgebiet (DDC):570: Biowissenschaften; Biologie
Lizenz:CC BY 2.5
Sprache:English
Anmerkungen:Finanziert durch den Open-Access-Publikationsfonds 2011/2012 der Deutschen Forschungsgemeinschaft (DFG) und der Westfälischen Wilhelms-Universität Münster (WWU Münster).
Format:PDF-Dokument
URN:urn:nbn:de:hbz:6-77379446593
Weitere Identifikatoren:DOI: doi:10.1371/journal.pone.0031410
Permalink:https://nbn-resolving.de/urn:nbn:de:hbz:6-77379446593
Onlinezugriff:journal.pone.0031410.pdf

Background: The quantity of transcriptome data is rapidly increasing for non-model organisms. As sequencing technology advances, focus shifts towards solving bioinformatic challenges, of which sequence read assembly is the first task. Recent studies have compared the performance of different software to establish a best practice for transcriptome assembly. Here, we adapted a simulation approach to evaluate specific features of assembly programs on 454 data. The novelty of our study is that the simulation allows us to calculate a model assembly as reference point for comparison. Findings: The simulation approach allows us to compare basic metrics of assemblies computed by different software applications (CAP3, MIRA, Newbler, and Oases) to a known optimal solution. We found MIRA and CAP3 are conservative in merging reads. This resulted in comparably high number of short contigs. In contrast, Newbler more readily merged reads into longer contigs, while Oases produced the overall shortest assembly. Due to the simulation approach, reads could be traced back to their correct placement within the transcriptome. Together with mapping reads onto the assembled contigs, we were able to evaluate ambiguity in the assemblies. This analysis further supported the conservative nature of MIRA and CAP3, which resulted in low proportions of chimeric contigs, but high redundancy. Newbler produced less redundancy, but the proportion of chimeric contigs was higher. Conclusion: Our evaluation of four assemblers suggested that MIRA and Newbler slightly outperformed the other programs, while showing contrasting characteristics. Oases did not perform very well on the 454 reads. Our evaluation indicated that the software was either conservative (MIRA) or liberal (Newbler) about merging reads into contigs. This suggested that in choosing an assembly program researchers should carefully consider their follow up analysis and consequences of the chosen approach to gain an assembly.