Reverse engineering environmental metatranscriptomes clarifies best practices for eukaryotic assembly

Abstract

Diverse communities of microbial eukaryotes in the global ocean provide a variety of essential ecosystem services, from primary production and carbon flow through trophic transfer to cooperation via symbioses. Increasingly, these communities are being understood through the lens of omics tools, which enable high-throughput processing of communities of multiple species. Metatranscriptomics in particular offers an understanding of real-time gene expression in microbial eukaryotic communities, providing a window into community metabolic activity. However, these approaches are rarely validated. A systematic validation of metatranscriptome assembly and annotation methods is necessary to appropriately assess (1) the fidelity of our community composition measurements and (2) the confidence of novel taxonomic and functional content discovered with metatranscriptomics. Here, we present euk rhythmic, an automated and customizable multi-assembler workflow for eukaryotic metatranscriptome assembly, annotation, and analysis, and validate the ability of the pipeline to recapitulate simulated eukaryotic community-level expression data treated as a gold standard. Finally, we apply this pipeline to two previously published metatranscriptomic datasets, showing both recapitulation of previous results as well as enhanced taxonomic and functional gene discovery.