Metatranscriptomics unmasks Mollusca virome with a remarkable presence of rhabdovirus in cephalopods

Abstract

[Introduction] Mollusks play a significant role in marine ecosystems and have economic value for aquaculture. Sometimes, unexpected and unexplained mortalities among mollusks have been described. The role of potential pathogens such as viruses remains unknown due to the lack of molluscan cell cultures, which is one of the major drawbacks to determining the viral role in such mortalities. Several oceanographic studies have suggested a high abundance of viruses in the oceans. Virus identification and understanding of viral interaction with organisms in marine ecosystems are in their infancy. Metatranscriptomics could become a useful tool to identify viruses using a shotgun approach and the growing number of viral genomes and sequences deposited in public databases. [Methods] In this work, several bioinformatics approaches were set up to screen Mollusca RNA sequences to find and confirm viral traces in their transcriptomes. This meta-analysis included an extensive search of SRA datasets belonging to mollusks available in the NCBI database, selecting a total of 55 SRA datasets that were further analyzed searching for viral sequences. [Results] Twenty-two bivalves, 19 cephalopods and 16 gastropods from 16 geographical origins and 17 different tissues were considered. The domain search approach was the most productive method to find viral sequences. This virus search showed that Cephalopoda samples (Idiosepius notoides and Amphioctopus fangsiao) exhibited the highest number of virus identifications. Some of the detected viral sequences were similar or identical to others previously identified. However, 33 putative new viruses were identified and analyzed phylogenetically when the RdRp domain was available. Specifically, Cephalopoda samples showed a considerable number of viruses belonging to the Rhabdoviridae family