Abstract:BackgroundNext-generation sequencing (NGS) technologies are arguably the most revolutionary technical development to join the list of tools available to molecular biologists since PCR. For researchers working with nonconventional model organisms one major problem with the currently dominant NGS platform (Illumina) stems from the obligatory fragmentation of nucleic acid material that occurs prior to sequencing during library preparation. This step creates a significant bioinformatic challenge for accurate de novo assembly of novel transcriptome data. This challenge becomes apparent when a variety of modern assembly tools (of which there is no shortage) are applied to the same raw NGS dataset. With the same assembly parameters these tools can generate markedly different assembly outputs.ResultsIn this study we present an approach that generates an optimized consensus de novo assembly of eukaryotic coding transcriptomes. This approach does not represent a new assembler, rather it combines the outputs of a variety of established assembly packages, and removes redundancy via a series of clustering steps. We test and validate our approach using Illumina datasets from six phylogenetically diverse eukaryotes (three metazoans, two plants and a yeast) and two simulated datasets derived from metazoan reference genome annotations. All of these datasets were assembled using three currently popular assembly packages (CLC, Trinity and IDBA-tran). In addition, we experimentally demonstrate that transcripts unique to one particular assembly package are likely to be bioinformatic artefacts. For all eight datasets our pipeline generates more concise transcriptomes that in fact possess more unique annotatable protein domains than any of the three individual assemblers we employed. Another measure of assembly completeness (using the purpose built BUSCO databases) also confirmed that our approach yields more information.ConclusionsOur approach yields coding transcriptome assemblies that are more likely to be closer to biological reality than any of the three individual assembly packages we investigated. This approach (freely available as a simple perl script) will be of use to researchers working with species for which there is little or no reference data against which the assembly of a transcriptome can be performed.

BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs

BUSCO: Assessing Genomic Data Quality and Beyond

BUSCO Update: Novel and Streamlined Workflows along with Broader and Deeper Phylogenetic Coverage for Scoring of Eukaryotic, Prokaryotic, and Viral Genomes

OrthoDB and BUSCO update: annotation of orthologs with wider sampling of genomes

GenomeQC: a quality assessment tool for genome assemblies and gene structure annotations

Versatile genome assembly evaluation with QUAST-LG

Gencube: Efficient retrieval, download, and unification of genomic data from leading biodiversity databases

Benchmarking of bioinformatics tools for the hybrid de novo assembly of human whole-genome sequencing data

BAGS: An automated Barcode, Audit & Grade System for DNA barcode reference libraries

Benchmarking of Hi-C tools for scaffolding de novo genome assemblies

Exploring Genome Characteristics and Sequence Quality Without a Reference

Genome size and chromosome number are critical metrics for accurate genome assembly assessment in Eukaryota

BusyBee Web: metagenomic data analysis by bootstrapped supervised binning and annotation

Combining independent de novo assemblies optimizes the coding transcriptome for nonconventional model eukaryotic organisms

Quest for Orthologs in the Era of Biodiversity Genomics

Baa.pl: A tool to evaluate de novo genome assemblies with RNA transcripts

Benchmarking of Long-Read Sequencing, Assemblers and Polishers for Yeast Genome

Identifying the Causes and Consequences of Assembly Gaps Using a Multiplatform Genome Assembly of a Bird-of-paradise.

Benchmarking Recent Computational Tools for DNA-binding Protein Identification

UnigeneFinder: An automated pipeline for gene calling from transcriptome assemblies without a reference genome

Estimation of genomic characteristics by analyzing k-mer frequency in de novo genome projects