The diversity of endogenous viral elements in insects

Unveiling the hidden viromes across the animal tree of life: insights from a taxonomic classification pipeline applied to invertebrates of 31 metazoan phyla

Invertebrates constitute the majority of animal species on Earth, including most disease-causing agents or vectors, with more diverse viromes when compared to vertebrates. Recent advancements in high-throughput sequencing have significantly expanded our understanding of invertebrate viruses, yet this knowledge remains biased toward a few well-studied animal lineages. In this study, we analyze invertebrate DNA and RNA viromes for 31 phyla using 417 publicly available RNA-Seq data sets from diverse environments in the marine-terrestrial and marine-freshwater gradients. This study aims to (i) estimate virome compositions at the family level for the first time across the animal tree of life, including the first exploration of the virome in several phyla, (ii) quantify the diversity of invertebrate viromes and characterize the structure of invertebrate-virus infection networks, and (iii) investigate host phylum and habitat influence on virome differences. Results showed that a set of few viral families of eukaryotes, comprising Retroviridae, Flaviviridae, and several families of giant DNA viruses, were ubiquitous and highly abundant. Nevertheless, some differences emerged between phyla, revealing for instance a less diverse virome in Ctenophora compared to the other animal phyla. Compositional analysis of the viromes showed that the host phylum explained over five times more variance in composition than its habitat. Moreover, significant similarities were observed between the viromes of some phylogenetically related phyla, which could highlight the influence of co-evolution in shaping invertebrate viromes.IMPORTANCEThis study significantly enhances our understanding of the global animal virome by characterizing the viromes of previously unexamined invertebrate lineages from a large number of animal phyla. It showcases the great diversity of viromes within each phylum and investigates the role of habitat shaping animal viral communities. Furthermore, our research identifies dominant virus families in invertebrates and distinguishes phyla with analogous viromes. This study sets the road toward a deeper understanding of the virome across the animal tree of life.

A tale of caution: How endogenous viral elements affect virus discovery in transcriptomic data

Large-scale metagenomic and -transcriptomic studies have revolutionized our understanding of viral diversity and abundance. In contrast, endogenous viral elements (EVEs), remnants of viral sequences integrated into host genomes, have received limited attention in the context of virus discovery, especially in RNA-Seq data. EVEs resemble their original viruses, a challenge that makes distinguishing between active infections and integrated remnants difficult, affecting virus classification and biases downstream analyses. Here, we systematically assess the effects of EVEs on a prototypical virus discovery pipeline, evaluate their impact on data integrity and classification accuracy, and provide some recommendations for better practices. We examined EVEs and exogenous viral sequences linked to Orthomyxoviridae, a diverse family of negative-sense segmented RNA viruses, in 13 genomic and 538 transcriptomic datasets of Culicinae mosquitoes. Our analysis revealed a substantial number of viral sequences in transcriptomic datasets. However, a significant portion appeared not to be exogenous viruses but transcripts derived from EVEs. Distinguishing between transcribed EVEs and exogenous virus sequences was especially difficult in samples with low viral abundance. For example, three transcribed EVEs showed full-length segments, devoid of frameshift and nonsense mutations, exhibiting sufficient mean read depths that qualify them as exogenous virus hits. Mapping reads on a host genome containing EVEs before assembly somewhat alleviated the EVE burden, but it led to a drastic reduction of viral hits and reduced quality of assemblies, especially in regions of the viral genome relatively similar to EVEs. Our study highlights that our knowledge of the genetic diversity of viruses can be altered by the underestimated presence of EVEs in transcriptomic datasets, leading to false positives and altered or missing sequence information. Thus, recognizing and addressing the influence of EVEs in virus discovery pipelines will be key in enhancing our ability to capture the full spectrum of viral diversity.

Co-option of a non-retroviral endogenous viral element in planthoppers

Non-retroviral endogenous viral elements (nrEVEs) are widely dispersed throughout the genomes of eukaryotes. Although nrEVEs are known to be involved in host antiviral immunity, it remains an open question whether they can be domesticated as functional proteins to serve cellular innovations in arthropods. In this study, we found that endogenous toti-like viral elements (ToEVEs) are ubiquitously integrated into the genomes of three planthopper species, with highly variable distributions and polymorphism levels in planthopper populations. Three ToEVEs display exon‒intron structures and active transcription, suggesting that they might have been domesticated by planthoppers. CRISPR/Cas9 experiments revealed that one ToEVE in Nilaparvata lugens, NlToEVE14, has been co-opted by its host and plays essential roles in planthopper development and fecundity. Large-scale analysis of ToEVEs in arthropod genomes indicated that the number of arthropod nrEVEs is currently underestimated and that they may contribute to the functional diversity of arthropod genes.

Endogenous Viral Elements in Ixodid Tick Genomes

The documentation of endogenous viral elements (EVEs; virus-derived genetic material integrated into the genome of a nonviral host) has offered insights into how arthropods respond to viral infection via RNA interference pathways. Small non-coding RNAs derived from EVE loci serve to direct RNAi pathways in limiting replication and infection from cognate viruses, thus benefiting the host's fitness and, potentially, vectorial capacity. Here we use informatic approaches to analyze nine available genome sequences of hard ticks (Acari: Ixodidae; Rhipicephalus sanguineus, R. microplus, R. annulatus, Ixodes ricinus, I. persulcatus, I. scapularis, Hyalomma asiaticum, Haemaphysalis longicornis, and Dermacentor silvarum) to identify endogenous viral elements and to illustrate the shared ancestry of all elements identified. Our results highlight a broad diversity of viral taxa as having given rise to 1234 identified EVEs in ticks, with Mononegavirales (specifically Rhabdoviridae) well-represented in this subset of hard ticks. Further investigation revealed extensive adintovirus integrations in several Ixodes species, the prevalence of Bunyavirales EVEs (notably not observed in mosquitoes), and the presence of several elements similar to known emerging human and veterinary pathogens. These results will inform subsequent work on current and past associations with tick species with regard to the viruses from which their "viral fossils" are derived and may serve as a reference for quality control of various tick-omics data that may suffer from misidentification of EVEs as viral genetic material.

Untangling an insect's virome from its endogenous viral elements

BACKGROUND

Insects are an important reservoir of viral biodiversity, but the vast majority of viruses associated with insects have not been discovered. Recent studies have employed high-throughput RNA sequencing, which has led to rapid advances in our understanding of insect viral diversity. However, insect genomes frequently contain transcribed endogenous viral elements (EVEs) with significant homology to exogenous viruses, complicating the use of RNAseq for viral discovery.

METHODS

In this study, we used a multi-pronged sequencing approach to study the virome of an important agricultural pest and prolific vector of plant pathogens, the potato aphid Macrosiphum euphorbiae. We first used rRNA-depleted RNAseq to characterize the microbes found in individual insects. We then used PCR screening to measure the frequency of two heritable viruses in a local aphid population. Lastly, we generated a quality draft genome assembly for M. euphorbiae using Illumina-corrected Nanopore sequencing to identify transcriptionally active EVEs in the host genome.

RESULTS

We found reads from two insect-specific viruses (a Flavivirus and an Ambidensovirus) in our RNAseq data, as well as a parasitoid virus (Bracovirus), a plant pathogenic virus (Tombusvirus), and two phages (Acinetobacter and APSE). However, our genome assembly showed that part of the 'virome' of this insect can be attributed to EVEs in the host genome.

CONCLUSION

Our work shows that EVEs have led to the misidentification of aphid viruses from RNAseq data, and we argue that this is a widespread challenge for the study of viral diversity in insects.

Endoparasitoid lifestyle promotes endogenization and domestication of dsDNA viruses

The accidental endogenization of viral elements within eukaryotic genomes can occasionally provide significant evolutionary benefits, giving rise to their long-term retention, that is, to viral domestication. For instance, in some endoparasitoid wasps (whose immature stages develop inside their hosts), the membrane-fusion property of double-stranded DNA viruses have been repeatedly domesticated following ancestral endogenizations. The endogenized genes provide female wasps with a delivery tool to inject virulence factors that are essential to the developmental success of their offspring. Because all known cases of viral domestication involve endoparasitic wasps, we hypothesized that this lifestyle, relying on a close interaction between individuals, may have promoted the endogenization and domestication of viruses. By analyzing the composition of 124 Hymenoptera genomes, spread over the diversity of this clade and including free-living, ecto, and endoparasitoid species, we tested this hypothesis. Our analysis first revealed that double-stranded DNA viruses, in comparison with other viral genomic structures (ssDNA, dsRNA, ssRNA), are more often endogenized and domesticated (that is, retained by selection) than expected from their estimated abundance in insect viral communities. Second, our analysis indicates that the rate at which dsDNA viruses are endogenized is higher in endoparasitoids than in ectoparasitoids or free-living hymenopterans, which also translates into more frequent events of domestication. Hence, these results are consistent with the hypothesis that the endoparasitoid lifestyle has facilitated the endogenization of dsDNA viruses, in turn, increasing the opportunities of domestications that now play a central role in the biology of many endoparasitoid lineages.

Detection of the Xanthi Chryso-like Virus in New Geographical Area and a Novel Arthropod Carrier

Here, we report on a serendipitous finding of a chryso-like virus associated with Culex pipiens mosquitos in the course of study aimed to detect and characterize West Nile virus (WNV) circulating in mosquitos in Serbia, Southern Europe. Upon initial detection of unexpected product in a PCR protocol for partial WNV NS5 gene amplification, further confirmation and identification was obtained through additional PCR and Sanger sequencing experiments. Bioinformatic and phylogenetic analysis identified the obtained sequences as Xanthi chryso-like virus (XCLV). The finding is particular for the fact that it associates XCLV with a new potential vector species and documents a novel geographical area of its distribution.

Identification of RNA Virus-Derived RdRp Sequences in Publicly Available Transcriptomic Data Sets

RNA viruses are abundant and highly diverse and infect all or most eukaryotic organisms. However, only a tiny fraction of the number and diversity of RNA virus species have been catalogued. To cost-effectively expand the diversity of known RNA virus sequences, we mined publicly available transcriptomic data sets. We developed 77 family-level Hidden Markov Model profiles for the viral RNA-dependent RNA polymerase (RdRp)-the only universal "hallmark" gene of RNA viruses. By using these to search the National Center for Biotechnology Information Transcriptome Shotgun Assembly database, we identified 5,867 contigs encoding RNA virus RdRps or fragments thereof and analyzed their diversity, taxonomic classification, phylogeny, and host associations. Our study expands the known diversity of RNA viruses, and the 77 curated RdRp Profile Hidden Markov Models provide a useful resource for the virus discovery community.

The Metagenomic Analysis of Viral Diversity in Colorado Potato Beetle Public NGS Data

The Colorado potato beetle (CPB) is one of the most serious insect pests due to its high ecological plasticity and ability to rapidly develop resistance to insecticides. The use of biological insecticides based on viruses is a promising approach to control insect pests, but the information on viruses which infect leaf feeding beetles is scarce. We performed a metagenomic analysis of 297 CPB genomic and transcriptomic samples from the public National Center for Biotechnology Information Sequence Read Archive (NCBI SRA) database. The reads that were not aligned to the reference genome were assembled with metaSPAdes, and 13314 selected contigs were analyzed with BLAST tools. The contigs and non-aligned reads were also analyzed with Kraken2 software. A total of 3137 virus-positive contigs were attributed to different viruses belonging to 6 types, 17 orders, and 32 families, matching over 97 viral species. The annotated sequences can be divided into several groups: those that are homologous to genetic sequences of insect viruses (Adintoviridae, Ascoviridae, Baculoviridae, Dicistroviridae, Chuviridae, Hytrosaviridae, Iflaviridae, Iridoviridae, Nimaviridae, Nudiviridae, Phasmaviridae, Picornaviridae, Polydnaviriformidae, Xinmoviridae etc.), plant viruses (Betaflexiviridae, Bromoviridae, Kitaviridae, Potyviridae), and endogenous retroviral elements (Retroviridae, Metaviridae). Additionally, the full-length genomes and near-full length genome sequences of several viruses were assembled. We also found sequences belonging to Bracoviriform viruses and, for the first time, experimentally validated the presence of bracoviral genetic fragments in the CPB genome. Our work represents the first attempt to discover the viral genetic material in CPB samples, and we hope that further studies will help to identify new viruses to extend the arsenal of biopesticides against CPB.

Putative Mitoviruses without In-Frame UGA(W) Codons: Evolutionary Implications

Mitoviruses are small vertically transmitted RNA viruses found in fungi, plants and animals. Taxonomically, a total of 105 species and 4 genera have been formally recognized by ICTV, and recently, 18 new putative species have been included in a new proposed genus. Transcriptomic and metatranscriptomic studies are a major source of countless new virus-like sequences that are continually being added to open databases and these may be good sources for identifying new putative mitoviruses. The search for mitovirus-like sequences in the NCBI databases resulted in the discovery of more than one hundred new putative mitoviruses, with important implications for taxonomy and also for the evolutionary scenario. Here, we propose the inclusion of four new putative members to the genus Kvaramitovirus, and the existence of a new large basally divergent lineage composed of 144 members that lack internal UGA codons (subfamily "Arkeomitovirinae"), a feature not shared by the vast majority of mitoviruses. Finally, a taxonomic categorization proposal and a detailed description of the evolutionary history of mitoviruses were carried out. This in silico study supports the hypothesis of the existence of a basally divergent lineage that could have had an impact on the early evolutionary history of mitoviruses.

Ancient origin of Jingchuvirales derived glycoproteins integrated in arthropod genomes

Endogenous virus elements (EVEs) are viral-derived sequences integrated into their host genomes. EVEs of the Jingchuvirales order were detected in a wide range of insect genomes covering several distantly related families. Moreover, Jingchuvirales-derived glycoproteins were recently associated by our group with the origin of a putative new retrovirus based on a glycoprotein captured by a mosquito retrotransposon. But, except for mosquitoes, there is a lack of a more detailed understanding of the endogenization mechanism, timing, and frequency per Jingchuvirales viral lineages. Here we screened Jingchuvirales glycoprotein-derived EVEs (Jg-EVEs) in eukaryotic genomes. We found six distinct endogenization events of Jg-EVEs, that belong to two out of five known Jingchuvirales families (Chuviridae and Natareviridae). For seven arthropod families bearing Jg-EVEs there is no register of bona fide circulating chuvirus infection. Hence, our results show that Jingchuvirales viruses infected or still infect these host families. Although we found abundant evidence of LTR-Gypsy retrotransposons fragments associated with the glycoprotein in Hymenoptera and other insect orders, our results show that the widespread distribution of Jingchuvirales glycoproteins in extant Arhtropods is a result of multiple ancient endogenization events and that these virus fossils are being vertically inherited in Arthropods genomes for millions of years.

An in silico analysis revealed a novel evolutionary lineage of putative mitoviruses

Mitoviruses (family Mitoviridae) are small capsid-less RNA viruses that replicate in the mitochondria of fungi and plants. However, to date, the only authentic animal mitovirus infecting an insect was identified as Lutzomyia longipalpis mitovirus 1 (LulMV1). Public databases of transcriptomic studies from several animals may be a good source for identifying the often missed mitoviruses. Consequently, a search of mitovirus-like transcripts at the NCBI transcriptome shotgun assembly (TSA) library, and a search for the mitoviruses previously recorded at the NCBI non-redundant (nr) protein sequences library, were performed in order to identify new mitovirus-like sequences associated with animals. In total, 10 new putative mitoviruses were identified in the TSA database and 8 putative mitoviruses in the nr protein database. To our knowledge, these results represent the first evidence of putative mitoviruses associated with poriferan, cnidarians, echinoderms, crustaceans, myriapods and arachnids. According to different phylogenetic inferences using the maximum likelihood method, these 18 putative mitoviruses form a robust monophyletic lineage with LulMV1, the only known animal-infecting mitovirus. These findings based on in silico procedures provide strong evidence for the existence of a clade of putative mitoviruses associated with animals, which has been provisionally named 'kvinmitovirus'.

First determination of DNA virus and some additional bacteria from Melophagus ovinus (sheep ked) in Tibet, China

Melophagus ovinus (sheep ked) is one of the common ectoparasites in sheep. In addition to causing direct damage to the host through biting and sucking blood, sheep ked is a potential vector of helminths, protozoa, bacteria, and viruses. Sheep M. ovinus samples from three regions in Tibet were selected for DNA extraction. The 16S rDNA V3-V4 hypervariable region was amplified, after genomic DNA fragmentation, Illumina Hiseq libraries were constructed. The 16S rRNA sequencing and viral metagenomics sequencing were separately conducted on the Illumina Novaseq 6000 platform and molecular biology software and platforms were employed to analyze the sequencing data. Illumina PE250 sequencing results demonstrated that the dominant bacteria phylum in M. ovinus from Tibet, China was Proteobacteria, where 29 bacteria genera were annotated. The dominant bacterial genera were Bartonella, Wolbachia, and Arsenophonus; Bartonella chomelii, Wolbachia spp., and Arsenophonus spp. were the dominant bacterial species in M. ovinus from Tibet, China. We also detected Kluyvera intermedia, Corynebacterium maris DSM 45190, Planomicrobium okeanokoites, and Rhodococcus erythropolis, of which the relative abundance of Kluyvera intermedia was high. Illumina Hiseq sequencing results demonstrated that 4 virus orders were detected in M. ovinus from Tibet, China, and 3 samples were annotated into 29 families, 30 families, and 28 families of viruses, respectively. Virus families related to vertebrates and insects mainly included Mimiviridae, Marseilleviridae, Poxviridae, Ascoviridae, Iridoviridae, Baculoviridae, Hytrosaviridae, Nudiviridae, Polydnaviridae, Adomaviridae, Asfarviridae, Hepeviridae, Herpesviridae, and Retroviridae; at the species level, the relative abundance of Tupanvirus_soda_lake, Klosneuvirus_KNV1, and Indivirus_ILV1 was higher. African swine fever virus and many poxviruses from the family Poxviridae were detected, albeit their relative abundance was low. The dominant bacterial phylum of M. ovinus from Tibet, China was Proteobacteria, and the dominant bacterial genera were Bartonella, Wolbachia, and Arsenophonus, where 23 out of 29 annotated bacteria genera were first reported in M. ovinus. Kluyvera intermedia, Corynebacterium maris DSM 45190, Planomicrobium okeanokoites, and Rhodococcus erythropolis were detected for the first time. All DNA viruses detected in this study have been reported in M. ovinus for the first time.

First Evidence of Past and Present Interactions between Viruses and the Black Soldier Fly, Hermetia illucens

Black soldier flies (BSFs, Hermetia illucens) are becoming a prominent research model encouraged by the insect as food and feed and waste bioconversion industries. Insect mass-rearing facilities are at risk from the spread of viruses, but so far, none have been described in BSFs. To fill this knowledge gap, a bioinformatic approach was undertaken to discover viruses specifically associated with BSFs. First, BSF genomes were screened for the presence of endogenous viral elements (EVEs). This led to the discovery and mapping of seven orthologous EVEs integrated into three BSF genomes originating from five viral families. Secondly, a virus discovery pipeline was used to screen BSF transcriptomes. This led to detecting a new exogenous totivirus that we named hermetia illucens totivirus 1 (HiTV1). Phylogenetic analyses showed this virus belongs to a clade of insect-specific totiviruses and is closely related to the largest EVE located on chromosome 1 of the BSF genome. Lastly, this EVE was found to express a small transcript in some BSFs infected by HiTV1. Altogether, this data mining study showed that far from being unscathed from viruses, BSFs bear traces of past interactions with several viral families and of present interactions with the exogenous HiTV1.

A Satellite dsRNA Attenuates the Induction of Helper Virus-Mediated Symptoms in Aspergillus flavus

Aspergillus flavus is an important fungal pathogen of animals and plants. Previously, we reported a novel partitivirus, Aspergillus flavus partitivirus 1 (AfPV1), infecting A. flavus. In this study, we obtained a small double-stranded (ds) RNA segment (734 bp), which is a satellite RNA of the helper virus, AfPV1. The presence of AfPV1 altered the colony morphology, decreased the number of conidiophores, created significantly larger vacuoles, and caused more sensitivity to osmotic, oxidative, and UV stresses in A. flavus, but the small RNA segment could attenuate the above symptoms caused by the helper virus AfPV1 in A. flavus. Moreover, AfPV1 infection reduced the pathogenicity of A. flavus in corn (Zea mays), honeycomb moth (Galleria mellonella), mice (Mus musculus), and the adhesion of conidia to host epithelial cells, and increased conidial death by macrophages. However, the small RNA segment could also attenuate the above symptoms caused by the helper virus AfPV1 in A. flavus, perhaps by reducing the genomic accumulation of the helper virus AfPV1 in A. flavus. We used this model to investigate transcriptional genes regulated by AfPV1 and the small RNA segment in A. flavus, and their role in generating different phenotypes. We found that the pathways of the genes regulated by AfPV1 in its host were similar to those of retroviral viruses. Therefore, some pathways may be of benefit to non-retroviral viral integration or endogenization into the genomes of its host. Moreover, some potential antiviral substances were also found in A. flavus using this system.