RNA virus EVEs in insect genomes

Highly divergent and diverse viral community infecting sylvatic mosquitoes from Northeast Brazil

Mosquitoes can transmit several pathogenic viruses to humans, but their natural viral community is also composed of a myriad of other viruses such as insect-specific viruses (ISVs) and those that infect symbiotic microorganisms. Besides a growing number of studies investigating the mosquito virome, the majority are focused on few urban species, and relatively little is known about the virome of sylvatic mosquitoes, particularly in high biodiverse biomes such as the Brazilian biomes. Here, we characterized the RNA virome of 10 sylvatic mosquito species from Atlantic forest remains at a sylvatic-urban interface in Northeast Brazil employing a metatranscriptomic approach. A total of 16 viral families were detected. The phylogenetic reconstructions of 14 viral families revealed that the majority of the sequences are putative ISVs. The phylogenetic positioning and, in most cases, the association with a high RNA-dependent RNA polymerase amino acid divergence from other known viruses suggests that the viruses characterized here represent at least 34 new viral species. Therefore, the sylvatic mosquito viral community is predominantly composed of highly divergent viruses highlighting the limited knowledge we still have about the natural virome of mosquitoes in general. Moreover, we found that none of the viruses recovered were shared between the species investigated, and only one showed high identity to a virus detected in a mosquito sampled in Peru, South America. These findings add further in-depth understanding about the interactions and coevolution between mosquitoes and viruses in natural environments.

IMPORTANCE

Mosquitoes are medically important insects as they transmit pathogenic viruses to humans and animals during blood feeding. However, their natural microbiota is also composed of a diverse set of viruses that cause no harm to the insect and other hosts, such as insect-specific viruses. In this study, we characterized the RNA virome of sylvatic mosquitoes from Northeast Brazil using unbiased metatranscriptomic sequencing and in-depth bioinformatic approaches. Our analysis revealed that these mosquitoes species harbor a diverse set of highly divergent viruses, and the majority comprises new viral species. Our findings revealed many new virus lineages characterized for the first time broadening our understanding about the natural interaction between mosquitoes and viruses. Finally, it also provided several complete genomes that warrant further assessment for mosquito and vertebrate host pathogenicity and their potential interference with pathogenic arboviruses.

Endogenous Bornavirus-like Elements in Bats: Evolutionary Insights from the Conserved Riboviral L-Gene in Microbats and Its Antisense Transcription in Myotis daubentonii

Bats are ecologically diverse vertebrates characterized by their ability to host a wide range of viruses without apparent illness and the presence of numerous endogenous viral elements (EVEs). EVEs are well preserved, expressed, and may affect host biology and immunity, but their role in bat immune system evolution remains unclear. Among EVEs, endogenous bornavirus-like elements (EBLs) are bornavirus sequences integrated into animal genomes. Here, we identified a novel EBL in the microbat Myotis daubentonii, EBLL-Cultervirus.10-MyoDau (short name is CV.10-MyoDau) that shows protein-level conservation with the L-protein of a Cultervirus (Wuhan sharpbelly bornavirus). Surprisingly, we discovered a transcript on the antisense strand comprising three exons, which we named AMCR-MyoDau. The active transcription in Myotis daubentonii tissues of AMCR-MyoDau, confirmed by RNA-Seq analysis and RT-PCR, highlights its potential role during viral infections. Using comparative genomics comprising 63 bat genomes, we demonstrate nucleotide-level conservation of CV.10-MyoDau and AMCR-MyoDau across various bat species and its detection in 22 Yangochiropera and 12 Yinpterochiroptera species. To the best of our knowledge, this marks the first occurrence of a conserved EVE shared among diverse bat species, which is accompanied by a conserved antisense transcript. This highlights the need for future research to explore the role of EVEs in shaping the evolution of bat immunity.

Rhabdoviral Endogenous Sequences Identified in the Leishmaniasis Vector Lutzomyia longipalpis Are Widespread in Sandflies from South America

Sandflies are known vectors of leishmaniasis. In the Old World, sandflies are also vectors of viruses while little is known about the capacity of New World insects to transmit viruses to humans. Here, we relate the identification of RNA sequences with homology to rhabdovirus nucleocapsids (NcPs) genes, initially in the Lutzomyia longipalpis LL5 cell lineage, named NcP1.1 and NcP2. The Rhabdoviridae family never retrotranscribes its RNA genome to DNA. The sequences here described were identified in cDNA and DNA from LL-5 cells and in adult insects indicating that they are transcribed endogenous viral elements (EVEs). The presence of NcP1.1 and NcP2 in the L. longipalpis genome was confirmed in silico. In addition to showing the genomic location of NcP1.1 and NcP2, we identified another rhabdoviral insertion named NcP1.2. Analysis of small RNA molecules derived from these sequences showed that NcP1.1 and NcP1.2 present a profile consistent with elements targeted by primary piRNAs, while NcP2 was restricted to the degradation profile. The presence of NcP1.1 and NcP2 was investigated in sandfly populations from South America and the Old World. These EVEs are shared by different sandfly populations in South America while none of the Old World species studied presented the insertions.

In vitro characterization of cell-fusing agent virus DNA forms in Aedes aegypti mosquitoes

How non-retroviral endogenous viral elements (EVEs) are established is a long-standing question. Viral DNA (vDNA) forms of RNA viruses are likely to be EVE precursors. Cell-fusing agent virus (CFAV) is a major insect-specific virus (ISV) in the Aedes aegypti mosquitoes and one of the few existing non-retroviral RNA viruses found as EVEs. We characterized CFAV-derived vDNA in the cell line to understand the mechanism of why current viruses are rarely endogenized. vDNA production was affected by cell culture media independent of CFAV replication. vDNAs that correspond to different regions covering the entire viral genome were detected, implying multiple initiation sites exist. A considerable proportion of vDNAs corresponded to ssDNA. Higher vDNA copies were detected in the cytoplasm than the nucleus. Our findings provide valuable insights into the intracellular characteristics of ISV-derived vDNAs, which will aid in understanding the underlying mechanisms of non-retroviral EVE formation.

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.

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.

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.

Integrated Jingmenvirus Polymerase Gene in Ixodes ricinus Genome

Members of the jingmenviruses group have been found in arthropods and mammals on all continents except Australia and Antarctica. Two viruses of this group were isolated from patients with fever after a tick bite. Using a nested RT-PCR assay targeting a jingmenvirus polymerase gene fragment, we screened ticks collected in seven regions of Russia and found that the abundant jingmenvirus-positive were of Ixodes ricinus species, with the prevalence ranging from 19.8% to 34.3%. In all cases, DNase/RNase treatment suggested that the detected molecule was DNA and subsequent next generation sequencing (NGS) proved that the viral polymerase gene was integrated in the I. ricinus genome. The copy number of the integrated polymerase gene was quantified by qPCR relative to the ITS2 gene and estimated as 1.32 copies per cell. At least three different genetic variants of the integrated polymerase gene were found in the territory of Russia. Phylogenetic analysis of the integrated jingmenvirus polymerase gene showed the highest similarity with the sequence of the correspondent gene obtained in Serbia from I. ricinus.

Local-scale virome depiction in Medellín, Colombia, supports significant differences between Aedes aegypti and Aedes albopictus

Aedes spp. comprise the primary group of mosquitoes that transmit arboviruses such as dengue, Zika, and chikungunya viruses to humans, and thus these insects pose a significant burden on public health worldwide. Advancements in next-generation sequencing and metagenomics have expanded our knowledge on the richness of RNA viruses harbored by arthropods such as Ae. aegypti and Ae. albopictus. Increasing evidence suggests that vector competence can be modified by the microbiome (comprising both bacteriome and virome) of mosquitoes present in endemic zones. Using an RNA-seq-based metataxonomic approach, this study determined the virome structure, Wolbachia presence and mitochondrial diversity of field-caught Ae. aegypti and Ae. albopictus mosquitoes in Medellín, Colombia, a municipality with a high incidence of mosquito-transmitted arboviruses. The two species are sympatric, but their core viromes differed considerably in richness, diversity, and abundance; although the community of viral species identified was large and complex, the viromes were dominated by few virus species. BLAST searches of assembled contigs suggested that at least 17 virus species (16 of which are insect-specific viruses [ISVs]) infect the Ae. aegypti population. Dengue virus 3 was detected in one sample and it was the only pathogenic virus detected. In Ae. albopictus, up to 11 ISVs and one plant virus were detected. Therefore, the virome composition appears to be species-specific. The bacterial endosymbiont Wolbachia was identified in all Ae. albopictus samples and in some Ae. aegypti samples collected after 2017. The presence of Wolbachia sp. in Ae. aegypti was not related to significant changes in the richness, diversity, or abundance of this mosquito’s virome, although it was related to an increase in the abundance of Aedes aegypti To virus 2 (Metaviridae). The mitochondrial diversity of these mosquitoes suggested that the Ae. aegypti population underwent a change that started in the second half of 2017, which coincides with the release of Wolbachia-infected mosquitoes in Medellín, indicating that the population of wMel-infected mosquitoes released has introduced new alleles into the wild Ae. aegypti population of Medellín. However, additional studies are required on the dispersal speed and intergenerational stability of wMel in Medellín and nearby areas as well as on the introgression of genetic variants in the native mosquito population.

Aphid Viruses: A Brief View of a Long History

Aphids are common agricultural pests with a wide range of hosts from agriculture to forestry plants. As known, aphids also serve as the major vectors to transmit plant viruses. Although numerous studies have focused on interactions between aphids and plant viruses, little is known about the aphid viruses, i.e., the insect viruses that are infectious to aphids. In the past four decades, several aphid viruses have been identified in diverse aphid species. In this review, we present a brief view of the aphid pathogenic viruses from several aspects, including classification of aphid viruses and characters of the viral genome, integration of viral sequences in host genomes, infection symptoms and influence on aphids, as well as host range and transmission modes. Taken together, these studies have increased our understanding of the rarely known aphid viruses, and will potentially contribute to the development of new strategies for controlling aphid populations.

Viral Metagenomic Data Analyses of Five New World Bat Species from Argentina: Identification of 35 Novel DNA Viruses

Bats are natural reservoirs of a variety of zoonotic viruses, many of which cause severe human diseases. Characterizing viruses of bats inhabiting different geographical regions is important for understanding their viral diversity and for detecting viral spillovers between animal species. Herein, the diversity of DNA viruses of five arthropodophagous bat species from Argentina was investigated using metagenomics. Fecal samples of 29 individuals from five species (Tadarida brasiliensis, Molossus molossus, Eumops bonariensis, Eumops patagonicus, and Eptesicus diminutus) living at two different geographical locations, were investigated. Enriched viral DNA was sequenced using Illumina MiSeq, and the reads were trimmed and filtered using several bioinformatic approaches. The resulting nucleotide sequences were subjected to viral taxonomic classification. In total, 4,520,370 read pairs were sequestered by sequencing, and 21.1% of them mapped to viral taxa. Circoviridae and Genomoviridae were the most prevalent among vertebrate viral families in all bat species included in this study. Samples from the T. brasiliensis colony exhibited lower viral diversity than samples from other species of New World bats. We characterized 35 complete genome sequences of novel viruses. These findings provide new insights into the global diversity of bat viruses in poorly studied species, contributing to prevention of emerging zoonotic diseases and to conservation policies for endangered species.