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

Objective As 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 ). Results The 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. ### Competing Interest Statement The authors have declared no competing interest. * List of abbreviations : DE : differential expression dph : days post hatch FPKM : Fragments Per Kilobase Mapped CPM : Counts Per Million GO : Gene Ontology