RNA virus discovery in insects

Identification and detection of known and new viruses in larvae of laboratory-reared fall armyworm, Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is a significant pest that causes economic losses worldwide. Many companies and research centers rear FAW to produce their microbiological-based products, such as viruses that target FAW. Nevertheless, colonies are vulnerable to collapsing, mainly due to uncontrolled and unexpected viral infections. In this work, dead FAWs exhibiting signs of viral infection were collected from unsuccessful attempts to propagate a baculovirus at a baculovirus production facility in Brazil. Total RNA was extracted and used to construct a cDNA library that was sequenced. The results showed the presence of five viruses, including three RNA viruses (alphanodavirus, rhabdovirus, and iflavirus) and two DNA viruses (densovirus and alphabaculovirus). To confirm the presence of the identified viruses in laboratory-reared FAWs, ten individual larvae from four accredited laboratories in Brazil were analyzed by RT-PCR with specific primers for each virus identified by sequencing, except the alphabaculovirus. Alphanodavirus and rhabdovirus were not detected in any of the four tested colonies, whereas the iflavirus was detected in two laboratories. A putative new densovirus was found in all samples. Accurate identification and timely detection of viruses that could disrupt the health of laboratory-reared insect colonies are crucial to ensure the production of high-quality biological products.

Novel viruses in the families Iflaviridae and Partitiviridae associated with the common eastern firefly Photinus pyralis

Fireflies are iconic insects that are under threat from environmental change. Knowledge of the viral diversity associated with natural firefly populations is important to our understanding of the basic biology of these insects and could be relevant to firefly conservation. We performed metatranscriptome sequencing of the Common Eastern Firefly ( Photinus pyralis) and assembled genomes for two new species of virus in the families Iflaviridae and Partitiviridae. We surveyed multiple individuals for these viruses using PCR, and we showed that both viruses are found at intermediate frequences in a natural population.

Diversity and transmissibility of RNA viruses in the small brown planthopper, Laodelphax striatellus

While a considerable number of viruses have been recently discovered in hematophagous insects, there remains insufficient research on virus diversity and their association with phytophagous insect hosts. In this study, we conducted a systematic investigation of the RNA virome in the small brown planthopper (SBPH), Laodelphax striatellus, an important vector of plant viruses. We successfully identified a total of 22 RNA viruses, including 17 novel viruses, from various families. These viruses were prevalent and abundant in SBPH populations, originating from the lab or field, with +ssRNA viruses composing the core SBPH viruses. Subsequent analysis revealed that the overall abundance of RNA viruses in SBPH remained relatively consistent across different developmental stages of the insects, although the titers of individual viruses varied among different insect tissues. This indicates a delicate balance between the viruses and their insect hosts. Interestingly, cross-species experiments confidently indicated that certain SBPH viruses could successfully infect and replicate in two other rice planthopper species (the brown planthopper and the white-backed planthopper) through microinjection. In conclusion, this study provides valuable insights into the RNA virome and its adaptability in a phytophagous insect, contributing to a better understanding of the intimate relationship between viruses and host insects.

IMPORTANCE

In the last decade, advances in the next-generation sequencing technology have unveiled a vast diversity of viruses in insects, particularly RNA viruses in hematophagous insects. However, research on virus diversity and their association with phytophagous insect hosts remains insufficient. This study presents a comprehensive analysis of the RNA virome in the small brown planthopper (SBPH), Laodelphax striatellus, a critical vector of plant viruses. The results indicated that the +ssRNA viruses, especially picorna-like viruses, comprised the core RNA viruses of SBPH that were prevalent in both laboratory and field populations. Moreover, a delicate balance was observed between the viruses and insect hosts. Significantly, some RNA viruses of SBPH could successfully infect and replicate in two other rice planthopper species belonging to different genera. This study provides valuable insights into the RNA virome and its adaptability in a phytophagous insect.

The multispecies stinkbug iflavirus Halyomorpha halys virus detected in the multispecies stinkbug egg parasitoid microwasp, Telenomus podisi (Ashmead) (Hymenoptera: Platygastridae)

Wasps are important parasitoids of stinkbugs and frequently exposed to various types of microorganisms through environmental contact and fecal-oral transmission route. Many parasitize stinkbug eggs and are commercially used in the field to control insect population. The parasitoid T. podisi is known for its high parasitism capacity and ability to target multiple species of stinkbugs. In this study we asked whether T. podisi exposed to eggs infected by a multispecies asymptomatic stinkbug virus, the Halyomorpha halys virus (HhV) would get infected. HhV is a geographically distributed multispecies iflavirus previously found to infect four stinkbug hosts, including three Brazilian species, Chinavia ubica, Euschistus heros and Diceraeus melacanthus, and T. posidi can parasitize all of them. As results, RT-PCR screening revealed positive samples for the HhV genome in two out of four tested pools of T. podisi, whereas the antigenome, indicative of replicative activity, was not detected. The wasps were raised in E. heros eggs that presented both the genome and the antigenome forms of the HhV genome. Subsequent RNA-deep sequencing of HhV positive T. podisi RNA pools yielded a complete genome of HhV with high coverage. Phylogenetic analysis positioned the isolate HhV-Tp (isolate Telenomus podisi) alongside with the stinkbug HhV. Analysis of transcriptomes from several hymenopteran species revealed HhV-Tp reads in four species. However, the transmission mechanism and the ecological significance of HhV remain elusive, warranting further studies to illuminate both the transmission process and its capacity for environmental propagation using T. podisi as a potential vector.

Insect-specific RNA virus affects the stylet penetration activity of brown citrus aphid (Aphis citricidus) to facilitate its transmission

Sap-sucking insects often transmit plant viruses but also carry insect viruses, which infect insects but not plants. The impact of such insect viruses on insect host biology and ecology is largely unknown. Here, we identified a novel insect-specific virus carried by brown citrus aphid (Aphis citricidus), which we tentatively named Aphis citricidus picornavirus (AcPV). Phylogenetic analysis discovered a monophyletic cluster with AcPV and other unassigned viruses, suggesting that these viruses represent a new family in order Picornavirales. Systemic infection with AcPV triggered aphid antiviral immunity mediated by RNA interference, resulting in asymptomatic tolerance. Importantly, we found that AcPV was transmitted horizontally by secretion of the salivary gland into the feeding sites of plants. AcPV influenced aphid stylet behavior during feeding and increased the time required for intercellular penetration, thus promoting its transmission among aphids with plants as an intermediate site. The gene expression results suggested that this mechanism was linked with transcription of salivary protein genes and plant defense hormone signaling. Together, our results show that the horizontal transmission of AcPV in brown citrus aphids evolved in a manner similar to that of the circulative transmission of plant viruses by insect vectors, thus providing a new ecological perspective on the activity of insect-specific viruses found in aphids and improving the understanding of insect virus ecology.

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.

RNA Viruses Are Prevalent and Active Tenants of the Predatory Mite Phytoseiulus persimilis (Acari: Phytoseiidae)

Many arthropod species harbor a diverse range of viruses. While much is known about pathogenic viruses of some economically important insects and arthropods involved in disease transmission, viruses associated with mites have rarely been studied. The main objective of this study was to characterize the virome of Phytoseiulus persimilis (Phytoseiidae), a predatory mite commercially used worldwide for the biological control of the key pest Tetranychus urticae (Tetranichidae). A combination of de novo transcriptome assembly and virion sequencing, revealed that RNA viruses are highly prevalent and active tenants of commercial populations of P. persimilis, comprising on average 9% of the mite's total mRNA. Seventeen RNA viruses dominated the mite's virome (i.e., were highly transcribed) with over half (n = 10) belonging to the order Picornavirales, + ssRNA viruses that infect a large range of hosts, including arthropods. Screening of the 17 dominant virus sequences in P. persimilis and T. urticae revealed that three viruses (two Picornavirales of the families Iflaviridae and Dicistroviridae, and one unclassified Riboviria) are unique to P. persimilis and three others (two unclassified Picornavirales and one unclassified Riboviria) are present in both mite species. Most of the sequences were related to viruses previously documented in economically important arthropods, while others have rarely been documented before in arthropods. These findings demonstrate that P. persimilis, like many other arthropods, harbors a diverse RNA virome, which might affect the mite's physiology and consequently its efficiency as a biological control agent.

Characterization of Iflavirus in the Red Flour Beetle, Tribolium castaneum (Coleoptera; Tenebrionidae)

Iflavirus is a group of viruses distributed mainly in arthropod species. We surveyed Tribolium castaneum iflavirus (TcIV) in different laboratory strains and in Sequence Read Archives (SRA) in GenBank. TcIV is highly specific to only T. castaneum and is not found in seven other Tenebrionid species, including the closely related species T. freemani. The same strains from different laboratories and different strains displayed largely different degrees of infections in the examination of 50 different lines by using Taqman-based quantitative PCR. We found that ~63% (27 out of 43 strains) of T. castaneum strains in different laboratories are positive for TcIV PCR with large degrees of variation, in the range of seven orders of magnitude, indicating that the TcIV is highly fluctuating depending on the rearing conditions. The TcIV was prevalent in the nervous system with low levels found in the gonad and gut. The transovarial transmission was supported in the experiment with surface-sterilized eggs. Interestingly, TcIV infection did not show observable pathogenicity. TcIV offers an opportunity to study the interaction between the virus and the immune system of this model beetle species.

Meta-transcriptomics reveals potential virus transfer between Aedes communis mosquitoes and their parasitic water mites

Arthropods harbor a largely undocumented diversity of RNA viruses. Some arthropods, like mosquitoes, can transmit viruses to vertebrates but are themselves parasitized by other arthropod species, such as mites. Very little is known about the viruses of these ectoparasites and how they move through the host-parasite relationship. To address this, we determined the virome of both mosquitoes and the mites that feed on them. The mosquito Aedes communis is an abundant and widely distributed species in Sweden, in northern Europe. These dipterans are commonly parasitized by water mite larvae (Trombidiformes: Mideopsidae) that are hypothesized to impose negative selection pressures on the mosquito by reducing fitness. In turn, viruses are dual-host agents in the mosquito-mite interaction. We determined the RNA virus diversity of mite-free and mite-detached mosquitoes, as well as their parasitic mites, using meta-transcriptomic sequencing. Our results revealed an extensive RNA virus diversity in both mites and mosquitoes, including thirty-seven putative novel RNA viruses that cover a wide taxonomic range. Notably, a high proportion of viruses (20/37) were shared between mites and mosquitoes, while a limited number of viruses were present in a single host. Comparisons of virus composition and abundance suggest potential virus transfer between mosquitoes and mites during their symbiotic interaction. These findings shed light on virome diversity and ecology in the context of arthropod host-parasite-virus relationships.

Identification of a Novel Densovirus in Aphid, and Uncovering the Possible Antiviral Process During Its Infection

Densoviruses (DVs) are single-stranded DNA viruses and exclusively happen in invertebrates. Most of DVs reported in insects are pathogenic to their native hosts, however, no pathogenic effect of them has been examined in vertebrates. Hence, DVs are the potential agents used in pest managements. Aphids are the primary vectors of plant viruses. In this study, we identified a novel DV in Chinese Sitobion miscanthi population, provisionally named “Sitobion miscanthi densovirus” (SmDV). Taxonomically, SmDV belongs to genus Hemiambidensovirus. In S. miscanthi, SmDV is hosted in diverse cells and can be horizontally transmitted via wheat feeding. Subject to SmDV, aphids activate their intrinsic antiviral autophagy pathway. Grouped with ascorbate and aldarate metabolism, chlorophyll metabolism, p450 related drug metabolism, and retinoid metabolism, aphids form a complex immune network response to the infection of SmDV. Obviously, it works as elder aphids still alive even they contain the highest examined concentration of SmDV. This study provides a foundation for the identifications of novel DVs, and further improves the understanding of the molecular interactions between insects and DVs.

The RNA virome of echinoderms

Echinoderms are a phylum of marine invertebrates that include model organisms, keystone species, and animals commercially harvested for seafood. Despite their scientific, ecological, and economic importance, there is little known about the diversity of RNA viruses that infect echinoderms compared to other invertebrates. We screened over 900 transcriptomes and viral metagenomes to characterize the RNA virome of 38 echinoderm species from all five classes (Crinoidea, Holothuroidea, Asteroidea, Ophiuroidea and Echinoidea). We identified 347 viral genome fragments that were classified to genera and families within nine viral orders - Picornavirales, Durnavirales, Martellivirales, Nodamuvirales, Reovirales, Amarillovirales, Ghabrivirales, Mononegavirales, and Hepelivirales. We compared the relative viral representation across three life stages (embryo, larvae, adult) and characterized the gene content of contigs which encoded complete or near-complete genomes. The proportion of viral reads in a given transcriptome was not found to significantly differ between life stages though the majority of viral contigs were discovered from transcriptomes of adult tissue. This study illuminates the biodiversity of RNA viruses from echinoderms, revealing the occurrence of viral groups in natural populations.

Novel viral RNA genomes of the vine mealybug Planococcus ficus

The vine mealybug, Planococcus ficus (Signoret, 1875), is the most important insect pest in growing areas of the grapevine Vitis vinifera L. in several countries, including Mexico. In Mexico, Baja California (B.C.) is the region with the highest production of V. vinifera L. grapes for industrial purposes. Recently, the diversity of viruses infecting insects only (insect-specific viruses) has been broadly explored to elucidate further ecological viral-host interactions in many insect species, which in some cases has resulted in the application of virus-based biological control agents for insect pests. However, a survey of the Pl. ficus virome has not been done yet. In the present study, we pooled Pl. ficus individuals collected through different vineyards of Ensenada, B.C., Mexico and analysed them by meta-transcriptomics. Novel nearly complete genomes of five RNA viruses were retrieved. These viruses were related to the Iflaviridae and Reoviridae families, and to the Picornavirales and Tolivirales orders. A new isolate belonging to the Dicistroviridae family was also found. Phylogenetic analyses showed that these putative viral genomes group with viruses having hemipteran (including a mealybug species) or other insect hosts, or with viruses associated with insects. Our results suggest that the identified novel RNA viruses could be insect-specific viruses of Pl. ficus. This work is the first insight into the Pl. ficus virome; it guarantees further studies aimed to characterize those viruses with potential for application in biological control of this economically important insect.

Diagnosis of Viral Diseases Using Deep Sequencing and Metagenomics Analyses

Viruses are ubiquitous in nature and exist in a variety of habitats. The advancement in sequencing technologies has revolutionized the understanding of viral biodiversity associated with plant diseases. Deep sequencing combined with metagenomics is a powerful approach that has proven to be revolutionary in the last decade and involves the direct analysis of viral genomes present in a diseased tissue sample. This protocol describes the details of RNA extraction and purification from wild rice plant and their yield, RNA purity, and integrity assessment. As a final step, bioinformatics data analysis including demultiplexing, quality control, de novo transcriptome assembly, taxonomic allocation and read mapping following Illumina HiSeq small and total RNA sequencing are described. Furthermore, the total RNAs extraction protocol and an additional ribosomal rRNAs depletion step which are significantly important for viral genomes construction are provided.

The Novel Agrotis ipsilon Nora Virus Confers Deleterious Effects to the Fitness of Spodoptera frugiperda (Lepidoptera: Noctuidae)

In the present study, we identified a novel, positive-sense single-stranded RNA virus in the Chinese black cutworm, Agrotis ipsilon. It has a genome length of 11,312 nucleotides, excluding the poly(A) tails, and contains five open reading frames. The ORF2 encodes the conserved domains of RNA helicase and RNA-dependent RNA polymerase, while ORF4 and 5 encode three viral proteins. Herein, the A. ipsilon virus was clustered with a Helicoverpa armigera Nora virus and was thus provisionally named “Agrotis ipsilon Nora virus” (AINV). AINV was successfully transmitted into a novel host, Spodoptera frugiperda, through injection, causing a stable infection. This found the possibility of horizontal AINV transmission among moths belonging to the same taxonomic family. Nonetheless, AINV infection was deleterious to S. frugiperda and mainly mediated by antiviral and amino acid metabolism-related pathways. Furthermore, the infection significantly increased the S. frugiperda larval period but significantly reduced its moth eclosion rate. It suggests that AINV is probably to be a parasitic virus of S. frugiperda.

Novel Dicistroviruses in an Unexpected Wide Range of Invertebrates

Dicistroviruses are members of a rapidly growing family of small RNA viruses. Related sequences have been discovered in many environmental samples, indicating that our knowledge about dicistrovirus diversity and host range is still limited. In this study, we performed a systematic search against the publicly available transcriptome database, and identified large numbers of dicistrovirus-like sequences in a wide variety of eukaryotic species. The origins of these sequences were 108 invertebrates (including 77 insect species belonging to 18 orders) and 11 plants, revealing new associations between dicistroviruses and hosts. Finally, 83 transcripts corresponding to nearly-complete viral genomes were retrieved from the RNA-seq data, of which most sequences showed limited similarity to known dicistroviruses and might present previously unreported virus species. Phylogenetic analysis suggested that horizontal virus transfer has occurred between diverse hosts and has important implications for dicistrovirus evolution. The results will provide new insight into the hidden diversity of the Dicistroviridae, and help us to better understand the viral evolution, host range and the possible way of transmission.

Novel RNA Viruses from the Transcriptome of Pheromone Glands in the Pink Bollworm Moth, Pectinophora gossypiella

Simple Summary

The pink bollworm, Pectinophora gossypiella (Lepidoptera: Gelechiidae), is a major pest of cotton. In this study, we analyzed the mRNA from pheromone glands of two populations in Israel. We found several virus sequences that were the same in these populations. We identified these viruses based on high-throughput sequencing data and analysis of the assembled transcripts. Through analysis of the sequences, we identified several unique viral sequences representing possible novel viral species. Two of the viral sequences were found in relatively high abundance in pheromone glands. One of the virus sequences was also found through analysis of previous transcriptome sequencing data from the midgut of pink bollworm larvae. This is the first report of these unique viral sequences found in the pink bollworm, and these viruses could be developed to help control this pest around the world, but more research is needed to determine their utility as biological control agents.

Abstract

In this study, we analyzed the transcriptome obtained from the pheromone gland isolated from two Israeli populations of the pink bollworm Pectinophora gossypiella to identify viral sequences. The lab population and the field samples carried the same viral sequences. We discovered four novel viruses: two positive-sense single-stranded RNA viruses, Pectinophora gossypiella virus 1 (PecgV1, a virus of Iflaviridae) and Pectinophora gossypiella virus 4 (PecgV4, unclassified), and two negative-sense single-stranded RNA viruses, Pectinophora gossypiella virus 2 (PecgV2, a virus of Phasmaviridae) and Pectinophora gossypiella virus 3 (PecgV3, a virus of Phenuiviridae). In addition, sequences derived from two negative-sense single-stranded RNA viruses that belong to Mononegavirales were found in the data. Analysis of previous transcriptome sequencing data derived from the midgut of pink bollworm larvae of a USA population only identified PecgV1, but no other viruses. High viral sequence coverages of PecgV1 and PecgV4 were observed in both field and lab populations. This is the first report of viral sequences discovered from the pink bollworm. Results from this investigation suggest that the pink bollworm harbors multiple viruses. Further investigation of the viral pathogens may help to develop novel pest management strategies for control of the pink bollworm.

Diverse RNA Viruses Discovered in Three Parasitoid Wasps of the Rice Weevil Sitophilus oryzae

The enormous diversity of RNA viruses in insects is continuously validated. Parasitoid wasps, as biocontrol insects which are widely used against insect pests in agroecosystems, may also carry many “good” RNA viruses.

In this study, many virus-like fragments were obtained from transcriptomes of three wasp species, including Anisopteromalus calandrae (8), Lariophagus distinguendus (3), and Theocolax elegans (18), which can parasitize and control rice weevil Sitophilus oryzae, a serious insect pest of farm-stored grains. By further bioinformatic analysis and sequencing, we identified six novel RNA viruses with complete genomes and named them WWPSRV-1, WWPSRV-2, AcPSRV-1, AcNSRV-1, AcNSRV-2, and LdNSRV-1. PCR-based detection revealed that WWPSRV-1 and WWPSRV-2 had the possibility of interspecies virus transmission, especially WWPSRV-2, which was also present in the rice weevil adults. Phylogenetically, three out of these six viruses appeared to be members of order Picornavirales: WWPSRV-1 belonged to unassigned virus families of this order, whereas WWPSRV-2 and AcPSRV-1 belonged to families Iflaviridae and Dicistroviridae, respectively. The conserved picornavirus-typical domains helicase, protease, and RNA-dependent RNA polymerase could be found in the nonstructural protein encoded by the three viruses, whose genomes consisted of the different numbers of open reading frames (ORFs). The other three RNA viruses could be classified to order Mononegavirales: AcNSRV-1 and AcNSRV-2 belonged to family Lispiviridae, whereas LdNSRV-1 belonged to a big family Rhabdoviridae. The genomes of the three viruses contained at least five ORFs, encoding deduced proteins in the following order: 3′-N-P-M-G-L-5′. All the ORFs were separated by conserved intergenic sequences which likely regulated the transcription termination and initiation. Our findings enhance the understanding of RNA viruses in weevil wasps and set the foundation for the future study of the association among weevils, weevil wasps, and RNA viruses.

IMPORTANCE The enormous diversity of RNA viruses in insects is continuously validated. Parasitoid wasps, as biocontrol insects which are widely used against insect pests in agroecosystems, may also carry many “good” RNA viruses. Some RNA viruses in parasitoid wasps have been reported to affect the host wasps or the wasps’ host. Here, six novel RNA viruses with complete genomes were identified in three parasitoid wasps of the rice weevil. One of these viruses was also detected in the rice weevil adults. Phylogenetically, WWPSRV-1 was the first unambiguous detection of Nora-like virus in insect parasitoids. WWPSRV-2 and AcPSRV-1 belong to families Iflaviridae and Dicistroviridae, some viruses of which can result in lethal infections in silkworms and honeybees. The other three RNA viruses belong to order Mononegavirales, which comprises many well-known insect-associated viruses.

Sequence diversity and evolution of a group of iflaviruses associated with ticks

We studied a group of tick-associated viruses with characteristics of members of the family Iflaviridae, a family of viruses frequently found in arthropods. Our aim was to gain insight into the evolutionary dynamics of this group of viruses, which may be linked to the biology of ticks. We explored assembled RNA-Seq data sets for different species of ticks. We identified members of five different iflavirus species, four of them novel, and discovered nine new genome sequences, including variants. Five variants represented a virus species associated with Ixodes ricinus. Unexpectedly, a sequence found in the Ixodes scapularis cell line ISE6 was nearly identical to the sequences of I. ricinus variants, suggesting a contamination of this cell line by I. ricinus material. Analysing patterns of substitutions between these variants, we detected a strong excess of synonymous mutations, suggesting evolution under strong positive selection. The phylogenies of the viruses and of their tick hosts were not congruent, suggesting recurrent host changes across tick genera during their evolution. Overall, our work constitutes a step in the understanding of the interactions between this family of viruses and ticks.

Evaluating assembly and variant calling software for strain-resolved analysis of large DNA viruses

Infection with human cytomegalovirus (HCMV) can cause severe complications in immunocompromised individuals and congenitally infected children. Characterizing heterogeneous viral populations and their evolution by high-throughput sequencing of clinical specimens requires the accurate assembly of individual strains or sequence variants and suitable variant calling methods. However, the performance of most methods has not been assessed for populations composed of low divergent viral strains with large genomes, such as HCMV. In an extensive benchmarking study, we evaluated 15 assemblers and 6 variant callers on 10 lab-generated benchmark data sets created with two different library preparation protocols, to identify best practices and challenges for analyzing such data. Most assemblers, especially metaSPAdes and IVA, performed well across a range of metrics in recovering abundant strains. However, only one, Savage, recovered low abundant strains and in a highly fragmented manner. Two variant callers, LoFreq and VarScan2, excelled across all strain abundances. Both shared a large fraction of false positive variant calls, which were strongly enriched in T to G changes in a 'G.G' context. The magnitude of this context-dependent systematic error is linked to the experimental protocol. We provide all benchmarking data, results and the entire benchmarking workflow named QuasiModo, Quasispecies Metric determination on omics, under the GNU General Public License v3.0 (https://github.com/hzi-bifo/Quasimodo), to enable full reproducibility and further benchmarking on these and other data.

Exploration of the virome of the European brown shrimp (Crangon crangon)

Crangon crangon is economically a very important species. Recently, promising culture attempts have been made, but a major problem is the uncontrollable mortality during the grow-out phase. As of yet, the life cycle of C. crangon is not closed in captivity so wild-caught individuals are used for further rearing. Therefore, it is important to investigate the virome of C. crangon both in wild-caught animals as in cultured animals. In recent years, next-generation-sequencing (NGS) technologies have been very important in the unravelling of the virome of a wide range of environments and matrices, such as soil, sea, potable water, but also of a wide range of animal species. This will be the first report of a virome study in C. crangon using NGS in combination with the NetoVIR protocol. The near complete genomes of 16 novel viruses were described, most of which were rather distantly related to unclassified viruses or viruses belonging to the Picornavirales, Bunyavirales Nudiviridae, Parvoviridae, Flaviviridae, Hepeviridae, Tombusviridae, Narnaviridae, Nodaviridae, Sobemovirus. A difference in virome composition was observed between muscle and hepatopancreatic tissue, suggesting a distinct tissue tropism of several of these viruses. Some differences in the viral composition were noted between the cultured and wild shrimp, which could indicate that in sub-optimal aquaculture conditions some viruses become more abundant. This research showed that a plethora of unknown viruses is present in C. crangon and that more research is needed to determine which virus is potentially dangerous for the culture of C. crangon.

Iflavirus Covert Infection Increases Susceptibility to Nucleopolyhedrovirus Disease in Spodoptera exigua

Naturally occurring covert infections in lepidopteran populations can involve multiple viruses with potentially different transmission strategies. In this study, we characterized covert infection by two RNA viruses, Spodoptera exigua iflavirus 1 (SeIV-1) and Spodoptera exigua iflavirus 2 (SeIV-2) (family Iflaviridae) that naturally infect populations of Spodoptera exigua, and examined their influence on susceptibility to patent disease by the nucleopolyhedrovirus Spodoptera exigua multiple nucleopolyhedrovirus (SeMNPV) (family Baculoviridae). The abundance of SeIV-1 genomes increased up to ten-thousand-fold across insect developmental stages after surface contamination of host eggs with a mixture of SeIV-1 and SeIV-2 particles, whereas the abundance of SeIV-2 remained constant across all developmental stages. Low levels of SeIV-2 infection were detected in all groups of insects, including those that hatched from surface-decontaminated egg masses. SeIV-1 infection resulted in reduced larval weight gain, and an unbalanced sex ratio, whereas larval developmental time, pupal weight, and adult emergence and fecundity were not significantly affected in infected adults. The inoculation of S. exigua egg masses with iflavirus, followed by a subsequent infection with SeMNPV, resulted in an additive effect on larval mortality. The 50% lethal concentration (LC50) of SeMNPV was reduced nearly 4-fold and the mean time to death was faster by 12 h in iflavirus-treated insects. These results suggest that inapparent iflavirus infections may be able to modulate the host response to a new pathogen, a finding that has particular relevance to the use of SeMNPV as the basis for biological pest control products.

First Report of Aleurocanthus spiniferus on Ailanthus altissima: Profiling of the Insect Microbiome and MicroRNAs

We report the first occurrence of the orange spiny whitefly (Aleurocanthus spiniferus; OSW) on the tree of heaven (Ailanthus altissima) in Bari, Apulia region, Italy. After our first observation in 2016, the infestation recurred regularly during the following years and expanded to the neighboring trees. Since then, we have also found the insect on numerous patches of the tree of heaven and other plant species in the Bari province. Nevertheless, the tree of heaven was not particularly threatened by the insect, so that a possible contribution by OSW for the control of such an invasive plant cannot be hypothesized hitherto. This work was also aimed at profiling the microbiome of OSW feeding on A. altissima. For this purpose, we used the denaturing gradient gel electrophoresis (DGGE) and the deep sequencing of small RNAs (sRNAs). Both techniques unveiled the presence of "Candidatus Portiera" (primary endosymbiont), Wolbachia sp. and Rickettsia sp., endosymbionts already reported for other Aleyrodidae. Deep sequencing data were analyzed by four computational pipelines in order to understand the reliability of the detection of fungi, bacteria, and viruses: Kraken, Kaiju, Velvet, and VelvetOptimiser. Some contigs assembled by Velvet or VelvetOptimiser were associated with insects, but not necessarily in the Aleurocanthus genus or Aleyrodidae family, suggesting the non-specificity of sRNAs or possible traces of parasitoids in the sample (e.g., Eretmocerus sp.). Finally, deep sequencing data were used to describe the microtranscriptome of OSW: 56 canonical and at least four high-confidence novel microRNAs (miRNAs) were identified. The overall miRNA abundance in OSW was in agreement with previous works on Bemisia tabaci, and bantam-3p, miR-276a-3p, miR-317-3p, miR-750-3p, and mir-8-3p were the most represented miRNAs.

Analysis of the RNA virome of basal hexapods

The diversity and evolution of RNA viruses has been well studied in arthropods and especially in insects. However, the diversity of RNA viruses in the basal hexapods has not been analysed yet. To better understand their diversity, evolutionary histories and genome organizations, we searched for RNA viruses in transcriptome and genome databases of basal hexapods. We discovered  40 novel RNA viruses, some of which are also present as endogenous viral elements derived from RNA viruses. Here, we demonstrated that basal hexapods host 14 RNA viral clades that have been recently identified in invertebrates. The following RNA viral clades are associated with basal hexapods: Reo, Partiti-Picobirna, Toti-Chryso, Mono-Chu, Bunya-Arena, Orthomyxo, Qinvirus, Picorna-Calici, Hepe-Virga, Narna-Levi, Tombus-Noda, Luteo-Sobemo, Permutotetra and Flavi. We have found representatives of the nine RNA viral clades that are present as endogenous genomic copies in the genomes of Machilis (Monocondylia) and Catajapyx (Diplura). Our study provided a first insight into the diversity of RNA viruses in basal hexapods and demonstrated that the basal hexapods possess quite high diversity of RNA viral clades.

A Novel RNA Virus in the Parasitoid Wasp Lysiphlebus fabarum: Genomic Structure, Prevalence, and Transmission

We report on a novel RNA virus infecting the wasp Lysiphlebus fabarum, a parasitoid of aphids. This virus, tentatively named "Lysiphlebus fabarum virus" (LysV), was discovered in transcriptome sequences of wasps from an experimental evolution study in which the parasitoids were allowed to adapt to aphid hosts (Aphis fabae) with or without resistance-conferring endosymbionts. Based on phylogenetic analyses of the viral RNA-dependent RNA polymerase (RdRp), LysV belongs to the Iflaviridae family in the order of the Picornavirales, with the closest known relatives all being parasitoid wasp-infecting viruses. We developed an endpoint PCR and a more sensitive qPCR assay to screen for LysV in field samples and laboratory lines. These screens verified the occurrence of LysV in wild parasitoids and identified the likely wild-source population for lab infections in Western Switzerland. Three viral haplotypes could be distinguished in wild populations, of which two were found in the laboratory. Both vertical and horizontal transmission of LysV were demonstrated experimentally, and repeated sampling of laboratory populations suggests that the virus can form persistent infections without obvious symptoms in infected wasps.

Codon Usage in the Iflaviridae Family Is Not Diverse Though the Family Members Are Isolated from Diverse Host Taxa

All iflavirus members belong to the unique genus, Iflavirus, of the family, Iflaviridae. The host taxa and sequence identities of these viruses are diverse. A codon usage bias, maintained by a balance between selection, mutation, and genetic drift, exists in a wide variety of organisms. We characterized the codon usage patterns of 44 iflavirus genomes that were isolated from the classes, Insecta, Arachnida, Mammalia, and Malacostraca. Iflaviruses lack a strong codon usage bias when they are evaluated using an effective number of codons. The odds ratios of the majority of dinucleotides are within the normal range. However, the dinucleotides at the 1st-2nd codon positions are more biased than those at the 2nd-3rd codon positions. Plots of effective numbers of codons, relative neutrality analysis, and PR2 bias analysis all indicate that selection pressure dominates mutations in shaping codon usage patterns in the family, Iflaviridae. When these viruses were grouped into their host taxa, we found that the indices, including the nucleotide composition, effective number of codons, relative synonymous codon usage, and the influencing factors behind the codon usage patterns, all show that there are non-significant differences between the six host-taxa-groups. Our results disagree with our assumption that diverse viruses should possess diverse codon usage patterns, suggesting that the nucleotide composition and codon usage in the family, Iflaviridae, are not host taxa-specific signatures.

Genome organization and molecular characterization of the three Formica exsecta viruses-FeV1, FeV2 and FeV4

We present the genome organization and molecular characterization of the three Formica exsecta viruses, along with ORF predictions, and functional annotation of genes. The Formica exsecta virus-4 (FeV4; GenBank ID: MF287670) is a newly discovered negative-sense single-stranded RNA virus representing the first identified member of order Mononegavirales in ants, whereas the Formica exsecta virus-1 (FeV1; GenBank ID: KF500001), and the Formica exsecta virus-2 (FeV2; GenBank ID: KF500002) are positive single-stranded RNA viruses initially identified (but not characterized) in our earlier study. The new virus FeV4 was found by re-analyzing data from a study published earlier. The Formica exsecta virus-4 genome is 9,866 bp in size, with an overall G + C content of 44.92%, and containing five predicted open reading frames (ORFs). Our bioinformatics analysis indicates that gaps are absent and the ORFs are complete, which based on our comparative genomics analysis suggests that the genomes are complete. Following the characterization, we validate virus infection for FeV1, FeV2 and FeV4 for the first time in field-collected worker ants. Some colonies were infected by multiple viruses, and the viruses were observed to infect all castes, and multiple life stages of workers and queens. Finally, highly similar viruses were expressed in adult workers and queens of six other Formica species: F. fusca, F. pressilabris, F. pratensis, F. aquilonia, F. truncorum and F. cinerea. This research indicates that viruses can be shared between ant species, but further studies on viral transmission are needed to understand viral infection pathways.

Insect-specific viruses: from discovery to potential translational applications

Over the past decade the scientific community has experienced a new age of virus discovery in arthropods in general, and in insects in particular. Next generation sequencing and advanced bioinformatics tools have provided new insights about insect viromes and viral evolution. In this review, we discuss some high-throughput sequencing technologies used to discover viruses in insects and the challenges raised in data interpretations. Additionally, the discovery of these novel viruses that are considered as insect-specific viruses (ISVs) has gained increasing attention in their potential use as biological agents. As example, we show how the ISV Nhumirim virus was used to reduce West Nile virus transmission when co-infecting the mosquito vector. We also discuss new translational opportunities of using ISVs to limit insect vector competence by using them to interfere with pathogen acquisition, to directly target the insect vector or to confer pathogen resistance by the insect vector.

Insect Virus Discovery by Metagenomic and Cell Culture-Based Approaches

Insects are the most abundant and diverse group of animals on earth, but our knowledge of their viruses is biased toward insect-borne viruses that cause disease in plants, animals, or humans. Recent metagenomic studies and systematic surveys of viruses in wild-caught insects have identified an unanticipated large repertoire of novel viruses and viral sequences. These include new members of existing clades, new clades, and even entirely new virus families. These studies greatly expand the known virosphere in insects, provide opportunities to study virus-host interactions, and generate new insights into virus evolution. In this chapter, we discuss the methods used to identify novel viruses in insects and highlight some notable surprises arising from these studies.

Discovery of Culex pipiens associated tunisia virus: a new ssRNA(+) virus representing a new insect associated virus family

In the global context of arboviral emergence, deep sequencing unlocks the discovery of new mosquito-borne viruses. Mosquitoes of the species Culex pipiens, C. torrentium, and C. hortensis were sampled from 22 locations worldwide for transcriptomic analyses. A virus discovery pipeline was used to analyze the dataset of 0.7 billion reads comprising 22 individual transcriptomes. Two closely related 6.8 kb viral genomes were identified in C. pipiens and named as Culex pipiens associated tunisia virus (CpATV) strains Ayed and Jedaida. The CpATV genome contained four ORFs. ORF1 possessed helicase and RNA-dependent RNA polymerase (RdRp) domains related to new viral sequences recently found mainly in dipterans. ORF2 and 4 contained a capsid protein domain showing strong homology with Virgaviridae plant viruses. ORF3 displayed similarities with eukaryotic Rhoptry domain and a merozoite surface protein (MSP7) domain only found in mosquito-transmitted Plasmodium, suggesting possible interactions between CpATV and vertebrate cells. Estimation of a strong purifying selection exerted on each ORFs and the presence of a polymorphism maintained in the coding region of ORF3 suggested that both CpATV sequences are genuine functional viruses. CpATV is part of an entirely new and highly diversified group of viruses recently found in insects, and that bears the genomic hallmarks of a new viral family.

Covert Infection of Insects by Baculoviruses

Baculoviruses (Baculoviridae) are occluded DNA viruses that are lethal pathogens of the larval stages of some lepidopterans, mosquitoes, and sawflies (phytophagous Hymenoptera). These viruses have been developed as biological insecticides for control of insect pests and as expression vectors in biotechnological applications. Natural and laboratory populations frequently harbor covert infections by baculoviruses, often at a prevalence exceeding 50%. Covert infection can comprise either non-productive latency or sublethal infection involving low level production of virus progeny. Latency in cell culture systems involves the expression of a small subset of viral genes. In contrast, covert infection in lepidopterans is associated with differential infection of cell types, modulation of virus gene expression and avoidance of immune system clearance. The molecular basis for covert infection may reside in the regulation of host-virus interactions through the action of microRNAs (miRNA). Initial findings suggest that insect nudiviruses and vertebrate herpesviruses may provide useful analogous models for exploring the mechanisms of covert infection by baculoviruses. These pathogens adopt mixed-mode transmission strategies that depend on the relative fitness gains that accrue through vertical and horizontal transmission. This facilitates virus persistence when opportunities for horizontal transmission are limited and ensures virus dispersal in migratory host species. However, when host survival is threatened by environmental or physiological stressors, latent or persistent infections can be activated to produce lethal disease, followed by horizontal transmission. Covert infection has also been implicated in population level effects on host-pathogen dynamics due to the reduced reproductive capacity of infected females. We conclude that covert infections provide many opportunities to examine the complexity of insect-virus pathosystems at the organismal level and to explore the evolutionary and ecological relationships of these pathogens with major crop and forest pests.

Viral Delivery of dsRNA for Control of Insect Agricultural Pests and Vectors of Human Disease: Prospects and Challenges

RNAi is applied as a new and safe method for pest control in agriculture but efficiency and specificity of delivery of dsRNA trigger remains a critical issue. Various agents have been proposed to augment dsRNA delivery, such as engineered micro-organisms and synthetic nanoparticles, but the use of viruses has received relatively little attention. Here we present a critical view of the potential of the use of recombinant viruses for efficient and specific delivery of dsRNA. First of all, it requires the availability of plasmid-based reverse genetics systems for virus production, of which an overview is presented. For RNA viruses, their application seems to be straightforward since dsRNA is produced as an intermediate molecule during viral replication, but DNA viruses also have potential through the production of RNA hairpins after transcription. However, application of recombinant virus for dsRNA delivery may not be straightforward in many cases, since viruses can encode RNAi suppressors, and virus-induced silencing effects can be determined by the properties of the encoded RNAi suppressor. An alternative is virus-like particles that retain the efficiency and specificity determinants of natural virions but have encapsidated non-replicating RNA. Finally, the use of viruses raises important safety issues which need to be addressed before application can proceed.

Viruses of insects reared for food and feed

The use of insects as food for humans or as feed for animals is an alternative for the increasing high demand for meat and has various environmental and social advantages over the traditional intensive production of livestock. Mass rearing of insects, under insect farming conditions or even in industrial settings, can be the key for a change in the way natural resources are utilized in order to produce meat, animal protein and a list of other valuable animal products. However, because insect mass rearing technology is relatively new, little is known about the different factors that determine the quality and yield of the production process. Obtaining such knowledge is crucial for the success of insect-based product development. One of the issues that is likely to compromise the success of insect rearing is the outbreak of insect diseases. In particular, viral diseases can be devastating for the productivity and the quality of mass rearing systems. Prevention and management of viral diseases imply the understanding of the different factors that interact in insect mass rearing. This publication provides an overview of the known viruses in insects most commonly reared for food and feed. Nowadays with large-scale sequencing techniques, new viruses are rapidly being discovered. We discuss factors affecting the emergence of viruses in mass rearing systems, along with virus transmission routes. Finally we provide an overview of the wide range of measures available to prevent and manage virus outbreaks in mass rearing systems, ranging from simple sanitation methods to highly sophisticated methods including RNAi and transgenics.

Aphis Glycines Virus 2, a Novel Insect Virus with a Unique Genome Structure

The invasive soybean aphid, Aphis glycines, is a major pest in soybeans, resulting in substantial economic loss. We analyzed the A. glycines transcriptome to identify sequences derived from viruses of A. glycines. We identified sequences derived from a novel virus named Aphis glycines virus 2 (ApGlV2). The assembled virus genome sequence was confirmed by reverse transcription polymerase chain reaction (RT-PCR) and Sanger sequencing, conserved domains were characterized, and distribution, and transmission examined. This virus has a positive sense, single-stranded RNA genome of ~4850 nt that encodes three proteins. The RNA-dependent RNA polymerase (RdRp) of ApGlV2 is a permuted RdRp similar to those of some tetraviruses, while the capsid protein is structurally similar to the capsid proteins of plant sobemoviruses. ApGlV2 also encodes a larger minor capsid protein, which is translated by a readthrough mechanism. ApGlV2 appears to be widespread in A. glycines populations and to persistently infect aphids with a 100% vertical transmission rate. ApGlV2 is susceptible to the antiviral RNA interference (RNAi) pathway. This virus, with its unique genome structure with both plant- and insect-virus characteristics, is of particular interest from an evolutionary standpoint.

Moku virus; a new Iflavirus found in wasps, honey bees and Varroa

There is an increasing global trend of emerging infectious diseases (EIDs) affecting a wide range of species, including honey bees. The global epidemic of the single stranded RNA Deformed wing virus (DWV), driven by the spread of Varroa destructor has been well documented. However, DWV is just one of many insect RNA viruses which infect a wide range of hosts. Here we report the full genome sequence of a novel Iflavirus named Moku virus (MV), discovered in the social wasp Vespula pensylvanica collected in Hawaii. The novel genome is 10,056 nucleotides long and encodes a polyprotein of 3050 amino acids. Phylogenetic analysis showed that MV is most closely related to Slow bee paralysis virus (SBPV), which is highly virulent in honey bees but rarely detected. Worryingly, MV sequences were also detected in honey bees and Varroa from the same location, suggesting that MV can also infect other hymenopteran and Acari hosts.

Application of Genomics for Understanding Plant Virus-Insect Vector Interactions and Insect Vector Control

The relationships between plant viruses and their vectors have evolved over the millennia, and yet, studies on viruses began <150 years ago and investigations into the virus and vector interactions even more recently. The advent of next generation sequencing, including rapid genome and transcriptome analysis, methods for evaluation of small RNAs, and the related disciplines of proteomics and metabolomics offer a significant shift in the ability to elucidate molecular mechanisms involved in virus infection and transmission by insect vectors. Genomic technologies offer an unprecedented opportunity to examine the response of insect vectors to the presence of ingested viruses through gene expression changes and altered biochemical pathways. This review focuses on the interactions between viruses and their whitefly or thrips vectors and on potential applications of genomics-driven control of the insect vectors. Recent studies have evaluated gene expression in vectors during feeding on plants infected with begomoviruses, criniviruses, and tospoviruses, which exhibit very different types of virus-vector interactions. These studies demonstrate the advantages of genomics and the potential complementary studies that rapidly advance our understanding of the biology of virus transmission by insect vectors and offer additional opportunities to design novel genetic strategies to manage insect vectors and the viruses they transmit.

Discovery of parvovirus-related sequences in an unexpected broad range of animals

Our knowledge of the genetic diversity and host ranges of viruses is fragmentary. This is particularly true for the Parvoviridae family. Genetic diversity studies of single stranded DNA viruses within this family have been largely focused on arthropod- and vertebrate-infecting species that cause diseases of humans and our domesticated animals: a focus that has biased our perception of parvovirus diversity. While metagenomics approaches could help rectify this bias, so too could transcriptomics studies. Large amounts of transcriptomic data are available for a diverse array of animal species and whenever this data has inadvertently been gathered from virus-infected individuals, it could contain detectable viral transcripts. We therefore performed a systematic search for parvovirus-related sequences (PRSs) within publicly available transcript, genome and protein databases and eleven new transcriptome datasets. This revealed 463 PRSs in the transcript databases of 118 animals. At least 41 of these PRSs are likely integrated within animal genomes in that they were also found within genomic sequence databases. Besides illuminating the ubiquity of parvoviruses, the number of parvoviral sequences discovered within public databases revealed numerous previously unknown parvovirus-host combinations; particularly in invertebrates. Our findings suggest that the host-ranges of extant parvoviruses might span the entire animal kingdom.

Trade-offs and mixed infections in an obligate-killing insect pathogen

Natural populations of pathogens are frequently composed of numerous interacting strains. Understanding what maintains this diversity remains a key focus of research in disease ecology. In addition, within-host pathogen dynamics can have a strong impact on both infection outcome and the evolution of pathogen virulence, and thus, understanding the impact of pathogen diversity is important for disease management. We compared eight genetically distinguishable variants from Spodoptera exempta nucleopolyhedrovirus (SpexNPV) isolated from the African armyworm, Spodoptera exempta. NPVs are obligate killers, and the vast majority of transmission stages are not released until after the host has died. The NPV variants differed significantly in their virulence and could be clustered into two groups based on their dose-response curves. They also differed in their speed of kill and productivity (transmission potential) for S. exempta. The mixed-genotype wild-type (WT) SpexNPV, from which each variant was isolated, was significantly more virulent than any individual variant and its mean mortality rate was within the fastest group of individual variants. However, the WT virus produced fewer new infectious stages than any single variant, which might reflect competition among the variants. A survival analysis, combining the mortality and speed of kill data, confirmed the superiority of the genetically mixed WT virus over any single variant. Spodoptera exempta larvae infected with WT SpexNPV were predicted to die 2·7 and 1·9 times faster than insects infected with isolates from either of the two clusters of genotypes. Theory suggests that there are likely to be trade-offs between pathogen fitness traits. Across all larvae, there was a negative linear relationship between virus yield and speed of kill, such that more rapid host death carried the cost of producing fewer transmission stages. We also found a near-significant relationship for the same trend at the intervariant level. However, there was no evidence for a significant relationship between the induced level of mortality and transmission potential (virus yield) or speed of kill.

A novel virus from Macrosiphum euphorbiae with similarities to members of the family Flaviviridae

A virus with a large genome was identified in the transcriptome of the potato aphid (Macrosiphum euphorbiae) and was named Macrosiphum euphorbiae virus 1 (MeV-1). The MeV-1 genome is 22 780 nt in size, including 3' and 5' non-coding regions, with a single large ORF encoding a putative polyprotein of 7333 aa. The C-terminal region of the predicted MeV-1 polyprotein contained sequences with similarities to helicase, methyltransferase and RNA-dependent RNA polymerase (RdRp) motifs, while the N-terminal region lacked any motifs including structural proteins. Phylogenetic analysis of the helicase placed MeV-1 close to pestiviruses, while the RdRp region placed it close to pestiviruses and flaviviruses, suggesting MeV-1 has a positive-polarity ssRNA genome and is a member of the family Flaviviridae. Since the MeV-1 genome is predicted to contain a methyltransferase, a gene present typically in flaviviruses but not pestiviruses, MeV-1 is likely a member of the genus Flavivirus. MeV-1 was present in nymphal and adult stages of the aphid, aphid saliva and plant tissues fed upon by aphids. However, the virus was unable to multiply and spread in tomato plants. In addition, dsRNA, the replication intermediate of RNA viruses, was isolated from virus-infected M. euphorbiae and not from tomato plants infested with the aphid. Furthermore, nymphs laid without exposure to infected plants harboured the virus, indicating that MeV-1 is an aphid-infecting virus likely transmitted transovarially. The virus was present in M. euphorbiae populations from Europe but not from North America and was absent in all other aphid species tested.

Diverse Array of New Viral Sequences Identified in Worldwide Populations of the Asian Citrus Psyllid (Diaphorina citri) Using Viral Metagenomics

The Asian citrus psyllid, Diaphorina citri, is the natural vector of the causal agent of Huanglongbing (HLB), or citrus greening disease. Together; HLB and D. citri represent a major threat to world citrus production. As there is no cure for HLB, insect vector management is considered one strategy to help control the disease, and D. citri viruses might be useful. In this study, we used a metagenomic approach to analyze viral sequences associated with the global population of D. citri. By sequencing small RNAs and the transcriptome coupled with bioinformatics analysis, we showed that the virus-like sequences of D. citri are diverse. We identified novel viral sequences belonging to the picornavirus superfamily, the Reoviridae, Parvoviridae, and Bunyaviridae families, and an unclassified positive-sense single-stranded RNA virus. Moreover, a Wolbachia prophage-related sequence was identified. This is the first comprehensive survey to assess the viral community from worldwide populations of an agricultural insect pest. Our results provide valuable information on new putative viruses, some of which may have the potential to be used as biocontrol agents.

IMPORTANCE

Insects have the most species of all animals, and are hosts to, and vectors of, a great variety of known and unknown viruses. Some of these most likely have the potential to be important fundamental and/or practical resources. In this study, we used high-throughput next-generation sequencing (NGS) technology and bioinformatics analysis to identify putative viruses associated with Diaphorina citri, the Asian citrus psyllid. D. citri is the vector of the bacterium causing Huanglongbing (HLB), currently the most serious threat to citrus worldwide. Here, we report several novel viral sequences associated with D. citri.