Comparison of de novo and reference genome-based transcriptome assembly pipelines for differential expression analysis of RNA sequencing data

Abstract

ObjectiveAs sequencing technologies become more accessible and bioinformatic tools improve, genomic resources are increasingly available for non-model species. Using a draft genome to guide transcriptome assembly from RNA sequencing data, rather than performing assembly de novo, affects downstream analyses. Yet, direct comparisons of these approaches are rare. Here, we compare the results of the standard de novo assembly pipeline ( Trinity) and two reference genome-based pipelines ( Tuxedo and the new Tuxedo) for differential expression and gene ontology enrichment analysis of a companion study on Atlantic cod (Gadus morhua). ResultsThe new Tuxedo pipeline produced a higher quality assembly than the Tuxedo suite. However, greater enrichment of Trinity-identified differentially expressed genes suggests that a higher proportion of them represent biologically meaningful differences in transcription, as opposed to transcriptional noise or false positives. Coupled with the ability to annotate novel loci, the increased sensitivity of the Trinity pipeline might make it preferable over the reference genome-based approaches for studies aimed at broadly characterizing variation in the magnitude of expression differences and biological processes. However, the new Tuxedo pipeline might be appropriate when a more conservative approach is warranted, such as for the identification of candidate genes.