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

Deep Sequencing Analysis of Virome Components, Viral Gene Expression and Antiviral RNAi Responses in Myzus persicae Aphids

The green peach aphid (Myzus persicae) is a generalist pest damaging crops and transmitting viral pathogens. Using Illumina sequencing of small (s)RNAs and poly(A)-enriched long RNAs, we analyzed aphid virome components, viral gene expression and antiviral RNA interference (RNAi) responses. Myzus persicae densovirus (family Parvoviridae), a single-stranded (ss)DNA virus persisting in the aphid population, produced 22 nucleotide sRNAs from both strands of the entire genome, including 5'- and 3'-inverted terminal repeats. These sRNAs likely represent Dicer-dependent small interfering (si)RNAs, whose double-stranded RNA precursors are produced by readthrough transcription beyond poly(A) signals of the converging leftward and rightward transcription units, mapped here with Illumina reads. Additionally, the densovirus produced 26-28 nucleotide sRNAs, comprising those enriched in 5'-terminal uridine and mostly derived from readthrough transcripts and those enriched in adenosine at position 10 from their 5'-end and mostly derived from viral mRNAs. These sRNAs likely represent PIWI-interacting RNAs generated by a ping-pong mechanism. A novel ssRNA virus, reconstructed from sRNAs and classified into the family Flaviviridae, co-persisted with the densovirus and produced 22 nucleotide siRNAs from the entire genome. Aphids fed on plants versus artificial diets exhibited distinct RNAi responses affecting densovirus transcription and flavivirus subgenomic RNA production. In aphids vectoring turnip yellows virus (family Solemoviridae), a complete virus genome was reconstituted from 21, 22 and 24 nucleotide viral siRNAs likely acquired with plant phloem sap. Collectively, deep-sequencing analysis allowed for the identification and de novo reconstruction of M. persicae virome components and uncovered RNAi mechanisms regulating viral gene expression and replication.

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.

Diversity of RNA viruses in agricultural insects

Recent advancements in next-generation sequencing (NGS) technology and bioinformatics tools have revealed a vast array of viral diversity in insects, particularly RNA viruses. However, our current understanding of insect RNA viruses has primarily focused on hematophagous insects due to their medical importance, while research on the viromes of agriculturally relevant insects remains limited. This comprehensive review aims to address the gap by providing an overview of the diversity of RNA viruses in agricultural pests and beneficial insects within the agricultural ecosystem. Based on the NCBI Virus Database, over eight hundred RNA viruses belonging to 39 viral families have been reported in more than three hundred agricultural insect species. These viruses are predominantly found in the insect orders of Hymenoptera, Hemiptera, Thysanoptera, Lepidoptera, Diptera, Coleoptera, and Orthoptera. These findings have significantly enriched our understanding of RNA viral diversity in agricultural insects. While further virome investigations are necessary to expand our knowledge to more insect species, it is crucial to explore the biological roles of these identified RNA viruses within insects in future studies. This review also highlights the limitations and challenges for the effective virus discovery through NGS and their potential solutions, which might facilitate for the development of innovative bioinformatic tools in the future.

The Characterization of Three Novel Insect-Specific Viruses Discovered in the Bean Bug, Riptortus pedestris

Insect-specific virus (ISV) is one of the most promising agents for the biological control of insects, which is abundantly distributed in hematophagous insects. However, few ISVs have been reported in Riptortus pedestris (Fabricius), one of the major pests threatening soybeans and causing great losses in yield and quality. In this work, field Riptortus pedestris was collected from six soybean-producing regions in China, and their virome was analyzed with the metatranscriptomic approach. Altogether, seven new insect RNA viruses were identified, three of which had complete RNA-dependent RNA polymerase (RdRp) and nearly full-length genome sequences, which were named Riptortus pedestris alphadrosrha-like virus 1 (RpALv1), Riptortus pedestris alphadrosrha-like virus 2 (RpALv2) and Riptortus pedestris almendra-like virus (RiALv). The three identified novel ISVs belonged to the family Rhabdoviridae, and phylogenetic tree analysis indicated that they were clustered into new distinct clades. Interestingly, the analysis of virus-derived small-interfering RNAs (vsiRNAs) indicated that only RiALv-derived siRNAs exhibited 22 nt length preference, whereas no clear 21 or 22 nt peaks were observed for RpALv1 and RpALv2, suggesting the complexity of siRNA-based antiviral immunity in R. pedestris. In conclusion, this study contributes to a better understanding of the microenvironment in R. pedestris and provides viral information for the development of potential soybean insect-specific biocontrol agents.

Discovery of a widespread presence bunyavirus that may have symbiont-like relationships with different species of aphids

Aphids are important agricultural pests, vectors of many plant viruses and have sophisticated relationships with symbiotic microorganisms. Abundant asymptomatic RNA viruses have been reported in aphids due to the application of RNA-seq, but aphid-virus interactions remain unclear. Bunyavirales is the most abundant RNA virus order, which can infect mammals, arthropods, and plants. However, many bunyaviruses have specific hosts, such as insects. Here, we discovered 18 viruses from 10 aphid species by RNA-seq. Importantly, a widespread presence bunyavirus, Aphid bunyavirus 1 (ABV-1), was determined to have a wide host range, infecting and replicating in all 10 tested aphid species. ABV-1 may be transmitted horizontally during feeding on plant leaves and vertically through reproduction. In a comparison of the physiological parameters of ABV-1^high and ABV-1^low strains of pea aphid, higher ABV-1 titers reduced the total nymphal duration and induced the reproduction. Moreover, viral titer significantly affected the lipid and protein contents in pea aphids. In summary, we proposed that ABV-1 may have stable symbiont-like relationships with aphids, and these observations may provide a new direction for studying bunyaviruses in aphids and establishing a model for virus-aphid interactions.

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.

Complete genome sequence of a novel nege-like virus in aphids (genus Indomegoura)

Background

Aphids are important vectors of numerous plant viruses. Besides plant viruses, a number of insect specific viruses (ISVs), such as nege/nege-like viruses, have been recently discovered in aphids of the genera Aphis, Rhopalosiphum, and Sitobion.

Findings

In this study, the complete genome sequence of a novel nege-like virus, tentatively named “Indomegoura nege-like virus 1” (INLV1), was identified in aphids of the genus Indomegoura. INLV1 possessed a single positive-stranded RNA genome with 8945 nucleotides, which was predicted to contain three typical open reading frames (ORFs) of negeviruses (including ORF1, ORF2, and ORF3), a 44-nt 5′ untranslated region (UTR) and a 98-nt 3′ UTR. Five conserved domains were predicted for INLV1, including an Alphavirus-like methyltransferase domain, a RNA virus helicase core domain, and a RNA-dependent RNA polymerase domain (RdRP) in ORF1, a DISB-ORF2_chro domain in ORF2, and a SP24 domain in ORF3. According to the maximum likelihood phylogenetic tree based on RdRP, INLV1 was grouped with barley aphid RNA virus 1 and Hubei virga-like virus 4, together with another two invertebrate viruses, which formed a distinct clade in the proposed group Centivirus. The alignment of RdRP domains for INLV1 and other nege/kita-like viruses suggested that RdRP of INLV1 contained the permuted C (GDD)- A [DX(4–5)D] –B [GX(2–3)TX(3)N] motifs, which were conserved in the Centivirus and Sandewavirus groups. Furthermore, the high abundance and typical characteristics of INLV1 derived small interfering RNAs clearly showed the active replication of INLV1 in the aphid Indomegoura.

Conclusion

INLV1 is the first nege-like virus infecting aphids of the genus Indomegoura. As far as we know, it is also the first ISV revealed in this aphid genus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-021-01552-w.

Viromics of extant insect orders unveil the evolution of the flavi-like superfamily

Insects are the most diversified and species-rich group of animals and harbor an immense diversity of viruses. Several taxa in the flavi-like superfamily, such as the genus Flavivirus, are associated with insects; however, systematic studies on insect virus genetic diversity are lacking, limiting our understanding of the evolution of the flavi-like superfamily. Here, we examined the diversity of flavi-like viruses within the most complete and up-to-date insect transcriptome collection comprising 1,243 insect species by employing a Flaviviridae RdRp profile hidden Markov model search. We identified seventy-six viral sequences in sixty-one species belonging to seventeen insect, one entognathan, and one arachnidan orders. Phylogenetic analyses revealed that twenty-seven sequences fell within the Flaviviridae phylogeny but did not group with established genera. Despite the large diversity of insect hosts studied, we only detected one virus in a blood-feeding insect, which branched within the genus Flavivirus, indicating that this genus likely diversified only in hematophagous arthropods. Nine new jingmenviruses with novel host associations were identified. One of the jingmenviruses established a deep rooting lineage additional to the insect- and tick-associated clades. Segment co-segregation phylogenies support the separation of tick- and insect-associated groups within jingmenviruses, with evidence for segment reassortment. In addition, fourteen viruses grouped with unclassified flaviviruses encompassing genome length of up to 20 kb. Species-specific clades for Hymenopteran- and Orthopteran-associated viruses were identified. Forty-nine viruses populated three highly diversified clades in distant relationship to Tombusviridae, a plant-infecting virus family, suggesting the detection of three previously unknown insect-associated families that contributed to tombusvirus evolution.

Structures of diverse poxin cGAMP nucleases reveal a widespread role for cGAS-STING evasion in host-pathogen conflict

DNA viruses in the family Poxviridae encode poxin enzymes that degrade the immune second messenger 2'3'-cGAMP to inhibit cGAS-STING immunity in mammalian cells. The closest homologs of poxin exist in the genomes of insect viruses suggesting a key mechanism of cGAS-STING evasion may have evolved outside of mammalian biology. Here we use a biochemical and structural approach to discover a broad family of 369 poxins encoded in diverse viral and animal genomes and define a prominent role for 2'3'-cGAMP cleavage in metazoan host-pathogen conflict. Structures of insect poxins reveal unexpected homology to flavivirus proteases and enable identification of functional self-cleaving poxins in RNA-virus polyproteins. Our data suggest widespread 2'3'-cGAMP signaling in insect antiviral immunity and explain how a family of cGAS-STING evasion enzymes evolved from viral proteases through gain of secondary nuclease activity. Poxin acquisition by poxviruses demonstrates the importance of environmental connections in shaping evolution of mammalian pathogens.

Deep Roots and Splendid Boughs of the Global Plant Virome

Land plants host a vast and diverse virome that is dominated by RNA viruses, with major additional contributions from reverse-transcribing and single-stranded (ss) DNA viruses. Here, we introduce the recently adopted comprehensive taxonomy of viruses based on phylogenomic analyses, as applied to the plant virome. We further trace the evolutionary ancestry of distinct plant virus lineages to primordial genetic mobile elements. We discuss the growing evidence of the pivotal role of horizontal virus transfer from invertebrates to plants during the terrestrialization of these organisms, which was enabled by the evolution of close ecological associations between these diverse organisms. It is our hope that the emerging big picture of the formation and global architecture of the plant virome will be of broad interest to plant biologists and virologists alike and will stimulate ever deeper inquiry into the fascinating field of virus-plant coevolution.

Virome Analysis of Aphid Populations That Infest the Barley Field: The Discovery of Two Novel Groups of Nege/Kita-Like Viruses and Other Novel RNA Viruses

Aphids (order Hemiptera) are important insect pests of crops and are also vectors of many plant viruses. However, little is known about aphid-infecting viruses, particularly their diversity and relationship to plant viruses. To investigate the aphid viromes, we performed deep sequencing analyses of the aphid transcriptomes from infested barley plants in a field in Japan. We discovered virus-like sequences related to nege/kita-, flavi-, tombus-, phenui-, mononega-, narna-, chryso-, partiti-, and luteoviruses. Using RT-PCR and sequence analyses, we determined almost complete sequences of seven nege/kitavirus-like virus genomes; one of which was a variant of the Wuhan house centipede virus (WHCV-1). The other six seem to belong to four novel viruses distantly related to Wuhan insect virus 9 (WhIV-9) or Hubei nege-like virus 4 (HVLV-4). We designated the four viruses as barley aphid RNA virus 1 to 4 (BARV-1 to -4). Moreover, some nege/kitavirus-like sequences were found by searches on the transcriptome shotgun assembly (TSA) libraries of arthropods and plants. Phylogenetic analyses showed that BARV-1 forms a clade with WHCV-1 and HVLV-4, whereas BARV-2 to -4 clustered with WhIV-9 and an aphid virus, Aphis glycines virus 3. Both virus groups (tentatively designated as Centivirus and Aphiglyvirus, respectively), together with arthropod virus-like TSAs, fill the phylogenetic gaps between the negeviruses and kitaviruses lineages. We also characterized the flavi/jingmen-like and tombus-like virus sequences as well as other RNA viruses, including six putative novel viruses, designated as barley aphid RNA viruses 5 to 10. Interestingly, we also discovered that some aphid-associated viruses, including nege/kita-like viruses, were present in different aphid species, raising a speculation that these viruses might be distributed across different aphid species with plants being the reservoirs. This study provides novel information on the diversity and spread of nege/kitavirus-related viruses and other RNA viruses that are associated with aphids.

Proteomics Computational Analyses Suggest That the Envelope Glycoproteins of Segmented Jingmen Flavi-Like Viruses are Class II Viral Fusion Proteins (b-Penetrenes) with Mucin-Like Domains

Jingmen viruses are newly described segmented flavi-like viruses that have a worldwide distribution in ticks and have been associated with febrile illnesses in humans. Computational analyses were used to predict that Jingmen flavi-like virus glycoproteins have structural features of class II viral fusion proteins, including an ectodomain consisting of beta-sheets and short alpha-helices, a fusion peptide with interfacial hydrophobicity and a three-domain architecture. Jingmen flavi-like virus glycoproteins have a sequence enriched in serine, threonine, and proline at the amino terminus, which is a feature of mucin-like domains. Several of the serines and threonines are predicted be modified by the addition of O-linked glycans. Some of the glycoproteins are predicted to have an additional mucin-like domain located prior to the transmembrane anchor, whereas others are predicted to have a stem consisting of two alpha-helices. The flavivirus envelope protein and Jingmen flavi-virus glycoproteins may have diverged from a common class II precursor glycoprotein with a mucin-like domain or domains acquired after divergence.

A Symbiotic Virus Facilitates Aphid Adaptation to Host Plants by Suppressing Jasmonic Acid Responses

Symbiotic viruses exist in many insects; however, their functions in host insects are not well understood. In this study, we explored the role of acyrthosiphon pisum virus (APV) in the interaction of its host aphid Acyrthosiphon pisum with plants. APV is primarily located in aphid salivary glands and gut and propagated in the insect. APV is horizontally transmitted to host plants during aphid feeding, but the virus does not replicate in the host plant. When the pea host race of aphids colonized two low-fitness plants, Medicago truncatula and Vicia villosa, the virus titers in both the aphids and plants significantly increased. Furthermore, APV infection strongly promoted the survival rate of the pea host race on V. villosa. Transcriptomic analysis showed that only 0.85% of aphid genes responded to APV infection when aphids fed on V. villosa, with a fold change in transcript levels of no more than fourfold. The improved survival due to APV infection was apparently related to the inhibitory effect of the virus on levels of phytohormone jasmonic acid (JA) and JA-isoleucine. Our data suggest a benefit of the symbiotic virus to its aphid host and demonstrate a novel case of symbiotic virus-mediated three-species interaction.

Origins and Evolution of the Global RNA Virome

Viruses with RNA genomes dominate the eukaryotic virome, reaching enormous diversity in animals and plants. The recent advances of metaviromics prompted us to perform a detailed phylogenomic reconstruction of the evolution of the dramatically expanded global RNA virome. The only universal gene among RNA viruses is the gene encoding the RNA-dependent RNA polymerase (RdRp). We developed an iterative computational procedure that alternates the RdRp phylogenetic tree construction with refinement of the underlying multiple-sequence alignments. The resulting tree encompasses 4,617 RNA virus RdRps and consists of 5 major branches; 2 of the branches include positive-sense RNA viruses, 1 is a mix of positive-sense (+) RNA and double-stranded RNA (dsRNA) viruses, and 2 consist of dsRNA and negative-sense (-) RNA viruses, respectively. This tree topology implies that dsRNA viruses evolved from +RNA viruses on at least two independent occasions, whereas -RNA viruses evolved from dsRNA viruses. Reconstruction of RNA virus evolution using the RdRp tree as the scaffold suggests that the last common ancestors of the major branches of +RNA viruses encoded only the RdRp and a single jelly-roll capsid protein. Subsequent evolution involved independent capture of additional genes, in particular, those encoding distinct RNA helicases, enabling replication of larger RNA genomes and facilitating virus genome expression and virus-host interactions. Phylogenomic analysis reveals extensive gene module exchange among diverse viruses and horizontal virus transfer between distantly related hosts. Although the network of evolutionary relationships within the RNA virome is bound to further expand, the present results call for a thorough reevaluation of the RNA virus taxonomy.IMPORTANCE The majority of the diverse viruses infecting eukaryotes have RNA genomes, including numerous human, animal, and plant pathogens. Recent advances of metagenomics have led to the discovery of many new groups of RNA viruses in a wide range of hosts. These findings enable a far more complete reconstruction of the evolution of RNA viruses than was attainable previously. This reconstruction reveals the relationships between different Baltimore classes of viruses and indicates extensive transfer of viruses between distantly related hosts, such as plants and animals. These results call for a major revision of the existing taxonomy of RNA viruses.

Sequence analysis of the potato aphid Macrosiphum euphorbiae transcriptome identified two new viruses

The potato aphid, Macrosiphum euphorbiae, is an important agricultural pest that causes economic losses to potato and tomato production. To establish the transcriptome for this aphid, RNA-Seq libraries constructed from aphids maintained on tomato plants were used in Illumina sequencing generating 52.6 million 75-105 bp paired-end reads. The reads were assembled using Velvet/Oases software with SEED preprocessing resulting in 22,137 contigs with an N50 value of 2,003bp. After removal of contigs from tomato host origin, 20,254 contigs were annotated using BLASTx searches against the non-redundant protein database from the National Center for Biotechnology Information (NCBI) as well as IntereProScan. This identified matches for 74% of the potato aphid contigs. The highest ranking hits for over 12,700 contigs were against the related pea aphid, Acyrthosiphon pisum. Gene Ontology (GO) was used to classify the identified M. euphorbiae contigs into biological process, cellular component and molecular function. Among the contigs, sequences of microbial origin were identified. Sixty five contigs were from the aphid bacterial obligate endosymbiont Buchnera aphidicola origin and two contigs had amino acid similarities to viruses. The latter two were named Macrosiphum euphorbiae virus 2 (MeV-2) and Macrosiphum euphorbiae virus 3 (MeV-3). The highest sequence identity to MeV-2 had the Dysaphis plantaginea densovirus, while to MeV-3 is the Hubei sobemo-like virus 49. Characterization of MeV-2 and MeV-3 indicated that both are transmitted vertically from adult aphids to nymphs. MeV-2 peptides were detected in the aphid saliva and only MeV-2 and not MeV-3 nucleic acids were detected inside tomato leaves exposed to virus-infected aphids. However, MeV-2 nucleic acids did not persist in tomato leaf tissues, after clearing the plants from aphids, indicating that MeV-2 is likely an aphid virus.

Metagenomics reshapes the concepts of RNA virus evolution by revealing extensive horizontal virus transfer

Virus metagenomics is a young research filed but it has already transformed our understanding of virus diversity and evolution, and illuminated at a new level the connections between virus evolution and the evolution and ecology of the hosts. In this review article, we examine the new picture of the evolution of RNA viruses, the dominant component of the eukaryotic virome, that is emerging from metagenomic data analysis. The major expansion of many groups of RNA viruses through metagenomics allowed the construction of substantially improved phylogenetic trees for the conserved virus genes, primarily, the RNA-dependent RNA polymerases (RdRp). In particular, a new superfamily of widespread, small positive-strand RNA viruses was delineated that unites tombus-like and noda-like viruses. Comparison of the genome architectures of RNA viruses discovered by metagenomics and by traditional methods reveals an extent of gene module shuffling among diverse virus genomes that far exceeds the previous appreciation of this evolutionary phenomenon. Most dramatically, inclusion of the metagenomic data in phylogenetic analyses of the RdRp resulted in the identification of numerous, strongly supported groups that encompass RNA viruses from diverse hosts including different groups of protists, animals and plants. Notwithstanding potential caveats, in particular, incomplete and uneven sampling of eukaryotic taxa, these highly unexpected findings reveal horizontal virus transfer (HVT) between diverse hosts as the central aspect of RNA virus evolution. The vast and diverse virome of invertebrates, particularly nematodes and arthropods, appears to be the reservoir, from which the viromes of plants and vertebrates evolved via multiple HVT events.

Complete Genome Sequence of the Largest Known Flavi-Like Virus, Diaphorina citri flavi-like virus, a Novel Virus of the Asian Citrus Psyllid, Diaphorina citri

A novel flavi-like virus tentatively named Diaphorina citri flavi-like virus (DcFLV) was identified in field populations of Diaphorina citri through small RNA and transcriptome sequencing followed by reverse transcription (RT)-PCR. We report here the complete nucleotide sequence and genome organization of DcFLV, the largest flavi-like virus identified to date.