Next Article in Journal
Coagulation System Activation for Targeting of COVID-19: Insights into Anticoagulants, Vaccine-Loaded Nanoparticles, and Hypercoagulability in COVID-19 Vaccines
Next Article in Special Issue
Densovirus Oil Suspension Significantly Improves the Efficacy and Duration of Larvicidal Activity against Aedes albopictus
Previous Article in Journal
Spike Gene Evolution and Immune Escape Mutations in Patients with Mild or Moderate Forms of COVID-19 and Treated with Monoclonal Antibodies Therapies
Previous Article in Special Issue
Actin Contributes to the Hyperexpression of Baculovirus Polyhedrin (polh) and p10 as a Component of Transcription Initiation Complex (TIC)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 14
Issue 2
10.3390/v14020227
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Virome Analysis Reveals Diverse and Divergent RNA Viruses in Wild Insect Pollinators in Beijing, China

by
Nannan Li
1,
Yizhao Huang
2,
Wei Li
2,*,† and
Shufa Xu
1,*,†
1
Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing 100093, China
2
College of Food Science and Pharmaceutical Engineering, Zaozhuang University, Zaozhuang 277160, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Viruses 2022, 14(2), 227; https://doi.org/10.3390/v14020227
Submission received: 31 December 2021 / Revised: 16 January 2022 / Accepted: 19 January 2022 / Published: 24 January 2022
(This article belongs to the Special Issue State-of-the-Art Insect Viruses Research in China)
Browse Figures
Versions Notes

Abstract

:
Insect pollinators provide major pollination services for wild plants and crops. Honeybee viruses can cause serious damage to honeybee colonies. However, viruses of other wild pollinating insects have yet to be fully explored. In the present study, we used RNA sequencing to investigate the viral diversity of 50 species of wild pollinating insects. A total of 3 pathogenic honeybee viruses, 8 previously reported viruses, and 26 novel viruses were identified in sequenced samples. Among these, 7 novel viruses were shown to be closely related to honeybee pathogenic viruses, and 4 were determined to have potential pathogenicity for their hosts. The viruses detected in wild insect pollinators were mainly from the order Picornavirales and the families Orthomyxoviridae, Sinhaliviridae, Rhabdoviridae, and Flaviviridae. Our study expanded the species range of known insect pollinator viruses, contributing to future efforts to protect economic honeybees and wild pollinating insects.

1. Introduction

Insect pollinators are necessary for most flowering plants, playing a key role in both wild plant reproduction and food production security [1,2]. Pollinating insects have a large species range, mainly belonging to three orders (i.e., Hymenoptera, Diptera, and Lepidoptera [3]. Pollinator biodiversity is critical for pollination quality in agricultural productivity and the conservation of the ecosystem [4,5], but the number of insect pollinators, including honeybees, is declining worldwide [6,7]. Several factors underlying pollinator decline are: heavy use of pesticides; worldwide spread of parasites, especially Varroa destructor; diseases caused by pathogenic viruses and other pathogens; monoculture cropping and plant biodiversity reduction; and competition between native and invasive species [6,8,9]. Additionally, many of these factors act simultaneously on insect pollinators and, due to their concurrent relationship, can exert additive or even synergistic harmful effects [6]. In the case of honeybees, several factors leading to the outbreak of viral diseases have already been identified: V. destructor parasitism [10], low-quality food [11], exposure to pesticides, and so forth [12,13,14,15,16,17].
Commercial pollination services have significantly contributed to the growth of agricultural production [18,19,20]. To protect the population of pollinators and avoid the spread of viruses, virus detection can help in designing a better management program to prevent the transmission of viruses from commercial insect pollinators to wild insect pollinators. The widespread application of next-generation sequencing technology in virus research has revealed a large number of novel viruses and contributed to the understanding of both the distribution and the evolutionary relationship of viruses [21,22]. For instance, the recent discovery of rhabdoviruses and flaviviruses in honeybees is based on sequencing technology [23]. Forty-two novel viruses in the order Picornavirales associated with Apis mellifera have been detected by RNA sequencing (RNA-seq) [24], and an orthomyxovirus has been detected in A. mellifera and V. destructor by using transcriptome sequencing. Detection of positive-sense RNA strands subsequently confirmed that this virus can infect both V. destructor and A. mellifera [25].
An RNA virus of an invertebrate was redefined by profiling the viromics of over 220 invertebrate species sampled across nine animal phyla; identified viruses fill major gaps in the RNA virus phylogeny and reveal an evolutionary history that is characterized by both host switching and codivergence [21]. Bees are ancient insects, and the evolution of viruses carried by them may provide us with a new perspective to understand the evolution of insect viruses [23]. Though many studies have focused on the diversity of honeybee viruses, their evolutionary relationships still need further research. For example, while the Lake Sinai virus (LSV) has been classified as Nodamuvirales, the LSV-containing clade has been shown to be distinct from the Nodaviridae family in an unrooted phylogenetic tree [26,27]. Interestingly, although Hypera postica-associated sinaivirus in the Sinhaliviridae family has been detected in Hypera postica, this virus is most closely related to the LSV [28]. Therefore, a study of viruses in host-related species has the potential to clarify as-yet-undefined viral evolutionary relationships.
Previous studies of insect pollinator viruses have primarily focused on individual viruses in honeybees [13,29,30]. Honeybee viruses transmitting to wild pollinating insects is a common phenomenon [31]. Deformed wing virus (DWV) has been proved to be able to infect wild bumblebees and cause wing deformities, and acute bee paralysis virus (ABPV), black queen cell virus (BQCV), and chronic bee paralysis virus (CBPV) were also detected in bumblebee species [32,33]. Understanding the diversity and composition of viruses in the category of pollinators as a whole can enable the development of more effective practices for the monitoring and control of the spread of viral diseases in pollinators [25,34,35]. In the present study, we employed next-generation sequencing technology for the first time to investigate viruses in wild insect pollinators and increase understanding of the diversity and evolutionary relationships of viruses in wild pollinating insects.

2. Materials and Methods

2.1. Sample Collection and Processing

Fifty specimens of wild insect pollinators (Hymenoptera: 38, Diptera: 10, and Amphiesmenoptera: 2) presenting no overt symptoms of viral diseases in the field were collected by an insect net on or next to flowering plants in Xiangshan, Beijing, from August to October 2020. The species of the specimens were identified according to their morphology, and all specimens were identified to the family level. These specimens were immediately stored in 5 mL centrifuge tubes in an ice box after capture and transferred to a refrigerator and maintained at −80 °C until processing. Specimens were then individually ground into powder in liquid nitrogen and divided into Pool A and Pool B for subsequent sequencing. Pool A contained Apidae: 10, Halictidae: 5, Syrphidae: 5, Ichneumonidae: 2, Megachilidae: 2, Asilidae: 1, Bombyliidae: 1, Calliphoridae: 1, Lycaenidae: 1, Sphingidae: 1, Tabanidae: 1, Tachinidae: 1, and Tenthredinidae: 1. Pool B contained Vespidae: 9, Sphecidae: 3, Scoliidae: 2, Apidae: 1, Halictidae: 1, Ichneumonidae: 1, and Pompilidae: 1. A complete list of specimens can be found in Supplementary Table S1.

2.2. Total RNA Extraction and Sequencing

The total RNA of the two mixed pools was extracted using the RNApure Total RNA Kit (RN0302; Aidlab Biotechnologies Co., Ltd., Beijing, China), according to the manufacturer’s recommendations. Sequencing libraries were generated using the VAHTS mRNA-seq v2 Library Prep Kit for Illumina (NR601-01; Vazyme, Nanjing, China), following the manufacturer’s recommendations, and index codes were added to attribute sequences to each sample. The libraries were sequenced on an Illumina NovaSeq platform to generate 150-bp paired-end reads, according to the manufacturer’s recommendations. We obtained a total of 37 Gb clean reads.

2.3. Sequence Assembly and Virus Detection

Sequencing reads were de novo assembled using Trinity (version 2.12.0) [34]. The contigs were compared with the GenBank nonredundant (NR) database with Diamond (version 2.0.8.146) [35]. We used a Python script to filter out nonvirus contigs, while viruslike contigs were mapped using Bowtie 2 (version 2.4.2) [36] to evaluate the contig quantity. Viral genome assembly results were manually corrected, and consensus sequences were generated by Geneious Primer (version 2020.0.4). Virus open reading frames (ORFs) were annotated based on the results of ORFfinder (version 0.4.3) [37] and the structure of the most closely related viral genome. Conserved domains of new viruses were identified using NCBI conserved domain database BLAST searches (version 2.11.0) [38].

2.4. Phylogenetic Analyses

Protein sequences of representatives from all viral species, as proposed by the International Committee on Taxonomy of Viruses (ICTV), were collected from GenBank and used for phylogenetic analyses. The RdRp region of novel viruses was used to retrieve related ICTV-classified viruses by BlastP. The RdRp domains of novel viruses and homologous viral proteins were aligned using MUSCLE (version 3.8.31) [39,40]. All ambiguously aligned regions were subsequently removed using Trimal (version 1.2rev59) [41]. The best-fit model of amino acid substitution in each dataset was determined using ModelTest [42]. Phylogenetic trees were inferred using the maximum likelihood method implemented in RAxML with 1000 bootstrap replicates [43]. Phylogenetic trees were displayed and annotated using FigTree (version 1.4.4). Pairwise sequence comparisons were performed by using MUSCLE to align the pairwise combinations of viral RdRp regions in each group and calculate the identity value of each alignment. The numerical matrix of each group was displayed by a heat map using the pheatmap package (version 1.0.12), and the pairwise sequence comparison results were used to verify the classification results of novel viruses based on their evolutionary trees.

2.5. Identification of Virus Hosts

We extracted RNA from each specimen powder sample before mixing using the RNApure Total RNA Kit (Aidlab). cDNA was prepared using the TaKaRa PrimeScript RT Reagent Kit (Perfect Real Time RR037A; Takara Biomedical Technology Co., Ltd., Beijing, China). The presence of a virus was determined using reverse transcription-polymerase chain reaction (RT-PCR), which was carried out with the first-strand cDNA products using the 2 × TSINGKE Master Mix (blue) (TSE004; Beijing TsingKe Biotech Co., Ltd., Beijing, China) in 25 µL reactions and specific primers. The specific primers were designed based on the assembled viral genome sequences (Supplementary Table S2). PCR products were confirmed for the presence of a target fragment by agarose gel electrophoresis and DNA sequencing.

2.6. Strand-Specific Detection of Novel Viruses

For identifying the reverse complementary strand of the novel viruses, the strand-specific primers (Supplementary Table S2) were designed in conformity with the method of a previous study, with slight modifications [44]. Specific primers were designed based on the assembled viral genome sequences, and a tag sequence (TCATGGTGGCGAATAA) was added to the forward primer. In reverse transcription, RNA was reverse-transcribed with tagged forward primers, and the control group with normal, random primers. Then the tag sequence as the forward primer and normal reverse primer was used for PCR detection of cDNA products. PCR products were confirmed for the presence of a target fragment by agarose gel electrophoresis and DNA sequencing.

3. Results

3.1. Identification of Known Viruses in Wild Pollinators

The Trinity software was used to assemble the sequencing data, and the Diamond program was used to annotate contigs before querying the GenBank NR database. In the present study, in total, eight known viruses were identified in the annotated results. The related contigs were further assembled to obtain consensus sequences (Table 1). Among these eight viruses, DWV, CBPV, and ABPV were found to be widely distributed honeybee pathogenic viruses; the Mayfield virus 1, previously detected in Bombus terrestris in Lebanon and the UK; the Vespa velutina-associated acypi-like virus, previously detected in Vespa velutina nigrithorax in France; the Scaldis River bee virus, previously detected in Osmia cornuta in Belgium; the Arboretum almendravirus, previously detected in the mosquito Psorophora albigenu in Peru; and the Hubei Diptera virus 6, previously detected in Diptera in China.

3.2. Identification of Novel RNA Viruses

A total of 26 novel RNA viruses were identified in this study (Table 2). ORFs of novel viruses were identified by the ORFfinder software and referenced to the genome structure of closely related viruses (Figure 1). Compared with the related viruses, ORFs of all novel viruses were complete or nearly complete. RdRp regions of the novel viruses were checked with the online BLASTp tool. The identity value with the most closely related viruses ranged from 34.8% to 73.9%, indicating that these viruses were significantly different from known viruses. The phylogenetic tree constructed using the RdRp region (Figure 2) and the pairwise identity comparison (Supplementary Figure S1) in each group showed that 21 novel viruses could be divided into 9 virus families: Iflaviridae: 5 (Figure 2A), Polycipiviridae: 1 (Figure 2A), Flaviviridae: 2 (Figure 2B), Tymoviridae: 1 (Figure 2C), Rhabdoviridae: 5 (Figure 2D), Nyamiviridae: 1 (Figure 2D), Kitaviridae: 4 (Figure 2E), Orthomyxoviridae: 1 (Figure 2F), Sinhaliviridae: 1 (Figure 2G).
Phylogenetic analysis showed that seven novel viruses have closer relationships with honeybee pathogenic viruses than other viruses within the tree, determined using multiple sequence alignment comparison by the MUSCLE program. Of the seven novel viruses, four novel viruses in the order Picornavirales were similar to honeybee pathogenic viruses. XPLV2 and XPLV4 share 45.94% and 46.39% amino acid identity with DWV (NC_004830.2), respectively. XPLV5 share 41.4% amino acid identity with slow bee paralysis virus (SBPV, NC_014137.1). XPLV1 share 47.26% amino acid identity with Sacbrood virus (SBV, NC_002066.1). A novel virus in Orthomyxoviridae, the Xiangshan orthomyxo-like virus (XOLV), share 34.98% amino acid identity in polymerase subunit PA with Varroa orthomyxovirus-1 (VOV-1, MK032467.1). The Xiangshan insect virus (XIV) and the Xiangshan sinhali-like virus (XSLV) share 35.1% and 38.26% amino acid identity in RdRp with CBPV (NC_010711.1) and LSV (NC_032433.1), respectively.

3.3. Identification of Complementary Strand of Novel Viruses

We selected seven novel viruses (XOLV, XTLV, XPLV2, XIV, XSLV, XPLV1, and XPLV4) to identify their reverse complementary chain to the host and used RT-PCR-tested insect samples to confirm the hosts of these seven viruses (Table 3). We further confirmed that these segments came from one sample by designing specific primers for each segment of three RNA viruses with segmented genomes (Figure 3A). Strand-specific RT-PCRs were employed to investigate whether the virus could replicate in the hosts. We screened for the presence of the positive-sense RNA strand of negative-sense genome viruses and the negative-sense RNA strand of positive-sense genome viruses. Specific primer information is listed in Supplementary Table S2. Results showed reverse complementary strands of XOLV, XPLV1, XSLV, and XPLV4 in the corresponding host samples, indicating that these viruses can replicate within and might be pathogenic to the hosts (Figure 3B).

4. Discussion

RNA-seq is a powerful tool that has been widely used to analyze the composition of known viruses, discover novel RNA viruses, and detect pathogenic viruses [22,45,46]. Many studies have focused on the viromics of honeybees using sequencing methods, significantly improving the scientific understanding of bee viruses [23,26,47,48,49,50]. Far fewer studies have been performed on wild insect pollinators, although their survival is also threatened by viruses [7,51,52,53]. By using the RNA-seq method in the present study, we performed the first preliminary exploration of the viromic characteristics of wild insect pollinators collected from flowering plants in Xiangshan, Beijing, China.
In the present study, DWV, ABPV, and CBPV were detected in pool B (which did not include Apis species), indicating that the three honeybee pathogenic viruses can also inhabit non-Apis wild insect pollinators. Previous studies have indicated that DWV and BQCV are able to replicate their genomes in bumblebees [53,54,55]. Therefore, it is confirmed that honeybee viruses can spread to other pollinators of the same family [56].
In this study, 5 previously reported insect viruses, the Mayfield virus 1 (host B. terrestris), the Vespa velutina-associated acypi-like virus (host V. velutina nigrithorax), the Scaldis River bee virus (host O. cornuta), the Hubei Diptera virus 6 (host Diptera), and the Arboretum almendravirus (host P. albigenu) were also detected in our samples, indicating that they may be prevalent and have been distributed in the corresponding hosts [21,57,58,59]. The Scaldis River bee virus and Arboretum almendravirus detected in this study shared 74.88% and 77.56% nucleotide identity with the reference sequences, but they shared higher amino acid identity in L protein (81% and 90.7%). In terms of similarity, they might not be classified as a new species. In the current study, they were classified as an “isolate” of the reference viruses.
We confirmed the existence of 26 novel RNA viruses by assembling sequencing data. Previous studies reported that bee-infecting viruses were primarily positive-sense single-stranded ssRNA (+ssRNA) viruses of the order Picornavirales [24,60]. Our study showed that, of the 26 identified novel viruses, 7 can be classified in the order Picornavirales, suggesting that Picornavirales viruses are ubiquitous in insect pollinators. Previous studies also suggested that viruses in the families Rhabdoviridae, Flaviviridae, Orthomyxoviridae, and Sinhaliviridae can infect honeybees [23,61]. Our study here also showed that novel viruses in wild pollinators belong to these families, suggesting that categories of viruses in wild insect pollinators are similar to categories of viruses found in honeybees. Moreover, 5 rhabdoviruses were detected in this study, clearly indicating that rhabdoviruses are more common in wild insect pollinators. Previous studies reported that the Kitaviridae family consists of plant viruses [62]. Four novel viruses belonging to this family were detected in this study, and we think these 4 viruses might originate from the plant’s pollen from which we collected our specimens. Seven novel viruses (XPLV2, XPLV4, XPLV5, XPLV1, XSLV, XIV, and XOLV) showed high amino acid sequence similarity with honeybee pathogenic viruses (DWV, SBPV, SBV, LSV, CBPV, and VOV-1). Single-stranded RNA (ssRNA) viruses produce reverse complementary strands during their replication in hosts. Several studies have already proved that an ssRNA virus can replicate in the host when the reverse complementary strand of the ssRNA virus is present [25,34,54,63]. By detecting the complementary strand of seven novel viral positive- or negative-sense RNA, 4 novel viruses (XOLV, XPLV1, XSLV, and XPLV4) could replicate in their corresponding host, which suggests that these viruses present the risk of virulence to their hosts and can closely interact with honeybees in a shared ecosystem [64].

5. Conclusions

We detected 26 novel RNA viruses in wild insect pollinators using RNA-seq, largely expanding previous understanding of the species range of insect pollinator viruses and their interspecies virulence risk and revealing the diversity of the virus composition of wild pollinating insects. These findings have applications for efforts to protect pollinator populations and ecosystems going forward, adding to knowledge of the risks faced by not only honeybees but all insect pollinators. Furthermore, this work suggests that there are still many undiscovered viruses in wild insect pollinators, indicating that more research is called for in this field.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v14020227/s1, Supplementary Figure S1: Pairwise sequence comparisons of novel viruses; Supplementary Table S1: Sample information of collected insects; Supplementary Table S2: Primer information of RT-PCR.

Author Contributions

Conceptualization, S.X. and W.L.; methodology, N.L.; software, N.L.; validation, S.X.; formal analysis, N.L.; investigation, Y.H.; resources, S.X.; data curation, N.L.; writing—original draft preparation, W.L. and N.L.; writing—review and editing, S.X.; visualization, N.L.; supervision, S.X. and W.L.; project administration, S.X.; funding acquisition, S.X. and W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Scientific Research Foundation for Ph.D. Programs of Zaozhuang University (1020708), the Agricultural Science and Technology Innovation Program (CAAS-ASTIP-2021-IAR), and the Modern Agro-Industry Technology Research System (CARS-44-KXJ6) in China.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw sequence data used for analysis are available in NCBI under the Sequence Read Archive (SRA), BioProject No. PRJNA728541 and SRA accession numbers SRR14554108–SRR14554111. SRR14554109 and SRR14554108 are from the same library of pool B. SRR14554109 and SRR14554110 are from the same library of pool A. The sequences of detected viruses were submitted to GenBank and given accession numbers OK491477–OK491521.

Acknowledgments

We thank W. Shi and T. Hu from Shandong First Medical University for their technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fijen, T.P.M.; Scheper, J.A.; Boom, T.M.; Janssen, N.; Raemakers, I.; Kleijn, D. Insect pollination is at least as important for marketable crop yield as plant quality in a seed crop. Ecol. Lett. 2018, 21, 1704–1713. [Google Scholar] [CrossRef] [PubMed]
  2. Ollerton, J.; Winfree, R.; Tarrant, S. How many flowering plants are pollinated by animals? Oikos 2011, 120, 321–326. [Google Scholar] [CrossRef]
  3. Corlett, R.T. Flower visitors and pollination in the Oriental (Indomalayan) Region. Biol. Rev. 2004, 79, 497–532. [Google Scholar] [CrossRef] [PubMed]
  4. Garibaldi, L.A.; Steffan-Dewenter, I.; Winfree, R.; Aizen, M.A.; Bommarco, R.; Cunningham, S.A.; Kremen, C.; Carvalheiro, L.G.; Harder, L.D.; Afik, O.; et al. Wild pollinators enhance fruit set of crops regardless of honey bee abundance. Science 2013, 339, 1608–1611. [Google Scholar] [CrossRef]
  5. Vasiliev, D.; Greenwood, S. Pollinator biodiversity and crop pollination in temperate ecosystems, implications for national pollinator conservation strategies: Mini review. Sci. Total Environ. 2020, 744, 140880. [Google Scholar] [CrossRef] [PubMed]
  6. Potts, S.G.; Biesmeijer, J.C.; Kremen, C.; Neumann, P.; Schweiger, O.; Kunin, W.E. Global pollinator declines: Trends, impacts and drivers. Trends Ecol. Evol. 2010, 25, 345–353. [Google Scholar] [CrossRef] [PubMed]
  7. Powney, G.D.; Carvell, C.; Edwards, M.; Morris, R.K.A.; Roy, H.E.; Woodcock, B.A.; Isaac, N.J.B. Widespread losses of pollinating insects in Britain. Nat. Commun. 2019, 10, 1018. [Google Scholar] [CrossRef]
  8. Cameron, S.A.; Sadd, B.M. Global trends in bumble bee health. Annu. Rev. Entomol. 2020, 65, 209–232. [Google Scholar] [CrossRef] [Green Version]
  9. Mitchell, E.A.D.; Mulhauser, B.; Mulot, M.; Mutabazi, A.; Glauser, G.; Aebi, A. A worldwide survey of neonicotinoids in honey. Science 2017, 358, 109–111. [Google Scholar] [CrossRef] [Green Version]
  10. Weaver, D.B.; Cantarel, B.L.; Elsik, C.G.; Boncristiani, D.L.; Evans, J.D. Multi-tiered analyses of honey bees that resist or succumb to parasitic mites and viruses. BMC Genom. 2021, 22, 720. [Google Scholar] [CrossRef]
  11. Dolezal, A.G.; Carrillo-Tripp, J.; Judd, T.M.; Allen Miller, W.; Bonning, B.C.; Toth, A.L. Interacting stressors matter: Diet quality and virus infection in honeybee health. R. Soc. Open Sci. 2019, 6, 181803. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Santamaria, J.; Villalobos, E.M.; Brettell, L.E.; Nikaido, S.; Graham, J.R.; Martin, S. Evidence of Varroa-mediated deformed wing virus spillover in Hawaii. J. Invertebr. Pathol. 2018, 151, 126–130. [Google Scholar] [CrossRef] [PubMed]
  13. Ratti, V.; Kevan, P.G.; Eberl, H.J. A mathematical model of the honeybee-varroa destructor-acute bee paralysis virus system with seasonal effects. Bull. Math. Biol. 2015, 77, 1493–1520. [Google Scholar] [CrossRef] [PubMed]
  14. Harwood, G.P.; Dolezal, A.G. Pesticide-virus interactions in honey bees: Challenges and opportunities for understanding drivers of bee declines. Viruses 2020, 12, 566. [Google Scholar] [CrossRef]
  15. Giacobino, A.; Molineri, A.I.; Pacini, A.; Fondevila, N.; Pietronave, H.; Rodriguez, G.; Palacio, A.; Bulacio Cagnolo, N.; Orellano, E.; Salto, C.E.; et al. Varroa destructor and viruses association in honey bee colonies under different climatic conditions. Environ. Microbiol. Rep. 2016, 8, 407–412. [Google Scholar] [CrossRef]
  16. Di Prisco, G.; Pennacchio, F.; Caprio, E.; Boncristiani, H.F., Jr.; Evans, J.D.; Chen, Y. Varroa destructor is an effective vector of Israeli acute paralysis virus in the honeybee, Apis mellifera. J. Gen. Virol. 2011, 92 Pt 1, 151–155. [Google Scholar] [CrossRef]
  17. Ryabov, E.V.; Wood, G.R.; Fannon, J.M.; Moore, J.D.; Bull, J.C.; Chandler, D.; Mead, A.; Burroughs, N.; Evans, D.J. A virulent strain of deformed wing virus (DWV) of honeybees (Apis mellifera) prevails after Varroa destructor-mediated, or in vitro, transmission. PLoS Pathog. 2014, 10, e1004230. [Google Scholar] [CrossRef] [Green Version]
  18. Klatt, B.K.; Holzschuh, A.; Westphal, C.; Clough, Y.; Smit, I.; Pawelzik, E.; Tscharntke, T. Bee pollination improves crop quality, shelf life and commercial value. Proc. Biol. Sci. 2014, 281, 20132440. [Google Scholar] [CrossRef]
  19. Saez, A.; Morales, J.M.; Morales, C.L.; Harder, L.D.; Aizen, M.A. The costs and benefits of pollinator dependence: Empirically based simulations predict raspberry fruit quality. Ecol. Appl. 2018, 28, 1215–1222. [Google Scholar] [CrossRef] [Green Version]
  20. Martin, C.D.; Fountain, M.T.; Brown, M.J.F. Varietal and seasonal differences in the effects of commercial bumblebees on fruit quality in strawberry crops. Agric. Ecosyst. Environ. 2019, 281, 124–133. [Google Scholar] [CrossRef]
  21. Shi, M.; Lin, X.D.; Tian, J.H.; Chen, L.J.; Chen, X.; Li, C.X.; Qin, X.C.; Li, J.; Cao, J.P.; Eden, J.S.; et al. Redefining the invertebrate RNA virosphere. Nature 2016, 540, 539–543. [Google Scholar] [CrossRef] [PubMed]
  22. Shi, M.; Lin, X.D.; Chen, X.; Tian, J.H.; Chen, L.J.; Li, K.; Wang, W.; Eden, J.S.; Shen, J.J.; Liu, L.; et al. The evolutionary history of vertebrate RNA viruses. Nature 2018, 556, 197–202. [Google Scholar] [CrossRef] [PubMed]
  23. Remnant, E.J.; Shi, M.; Buchmann, G.; Blacquiere, T.; Holmes, E.C.; Beekman, M.; Ashe, A. A diverse range of novel RNA viruses in geographically distinct honey bee populations. J. Virol. 2017, 91, e00158-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Roberts, J.M.K.; Anderson, D.L.; Durr, P.A. Metagenomic analysis of Varroa-free Australian honey bees (Apis mellifera) shows a diverse Picornavirales virome. J. Gen. Virol. 2018, 99, 818–826. [Google Scholar] [CrossRef] [PubMed]
  25. Levin, S.; Sela, N.; Erez, T.; Nestel, D.; Pettis, J.; Neumann, P.; Chejanovsky, N. New viruses from the ectoparasite mite varroa destructor infesting Apis mellifera and Apis cerana. Viruses 2019, 11, 94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Cornman, R.S. Relative abundance and molecular evolution of Lake Sinai Virus (Sinaivirus) clades. PeerJ 2019, 7, e6305. [Google Scholar] [CrossRef] [Green Version]
  27. Daughenbaugh, K.; Martin, M.; Brutscher, L.; Cavigli, I.; Garcia, E.; Lavin, M.; Flenniken, M. Honey bee infecting lake sinai viruses. Viruses 2015, 7, 3285–3309. [Google Scholar] [CrossRef] [Green Version]
  28. Francois, S.; Antoine-Lorquin, A.; Kulikowski, M.; Frayssinet, M.; Filloux, D.; Fernandez, E.; Roumagnac, P.; Froissart, R.; Ogliastro, M. Characterisation of the viral community associated with the Alfalfa Weevil (Hypera postica) and its host plant, Alfalfa (Medicago sativa). Viruses 2021, 13, 791. [Google Scholar] [CrossRef]
  29. Kevill, J.L.; de Souza, F.S.; Sharples, C.; Oliver, R.; Schroeder, D.C.; Martin, S.J. DWV-A lethal to honey bees (Apis mellifera): A colony level survey of DWV variants (A, B, and C) in England, Wales, and 32 states across the US. Viruses 2019, 11, 426. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, J.Q.; Zhang, Y.; Han, R.C. The high-throughput production of dsRNA against sacbrood virus for use in the honey bee Apis cerana (Hymenoptera: Apidae). Virus Genes 2016, 52, 698–705. [Google Scholar] [CrossRef]
  31. McMahon, D.P.; Furst, M.A.; Caspar, J.; Theodorou, P.; Brown, M.J.F.; Paxton, R.J. A sting in the spit: Widespread cross-infection of multiple RNA viruses across wild and managed bees. J. Anim. Ecol. 2015, 84, 615–624. [Google Scholar] [CrossRef] [PubMed]
  32. Genersch, E.; Yue, C.; Fries, I.; de Miranda, J.R. Detection of Deformed wing virus, a honey bee viral pathogen, in bumble bees (Bombus terrestris and Bombus pascuorum) with wing deformities. J. Invertebr. Pathol. 2006, 91, 61–63. [Google Scholar] [CrossRef] [PubMed]
  33. Tlak Gajger, I.; Simenc, L.; Toplak, I. The first detection and genetic characterization of four different honeybee viruses in wild bumblebees from Croatia. Pathogens 2021, 10, 808. [Google Scholar] [CrossRef] [PubMed]
  34. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.D.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef]
  36. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [Green Version]
  37. Rombel, I.T.; Sykes, K.F.; Rayner, S.; Johnston, S.A. ORF-FINDER: A vector for high-throughput gene identification. Gene 2002, 282, 33–41. [Google Scholar] [CrossRef]
  38. Marchler-Bauer, A.; Derbyshire, M.K.; Gonzales, N.R.; Lu, S.N.; Chitsaz, F.; Geer, L.Y.; Geer, R.C.; He, J.; Gwadz, M.; Hurwitz, D.I.; et al. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 2015, 43, D222–D226. [Google Scholar] [CrossRef] [Green Version]
  39. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  40. Edgar, R.C. MUSCLE: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004, 5, 113. [Google Scholar] [CrossRef] [Green Version]
  41. Capella-Gutierrez, S.; Silla-Martinez, J.M.; Gabaldon, T. trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 2009, 25, 1972–1973. [Google Scholar] [CrossRef] [PubMed]
  42. Darriba, D.; Posada, D.; Kozlov, A.M.; Stamatakis, A.; Morel, B.; Flouri, T. ModelTest-NG: A new and scalable tool for the selection of DNA and protein evolutionary models. Mol. Biol. Evol. 2020, 37, 291–294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Kozlov, A.M.; Darriba, D.; Flouri, T.; Morel, B.; Stamatakis, A. RAxML-NG: A fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 2019, 35, 4453–4455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Craggs, J.K.; Ball, J.K.; Thomson, B.J.; Irving, W.L.; Grabowska, A.M. Development of a strand-specific RT-PCR based assay to detect the replicative form of hepatitis C virus RNA. J. Virol. Methods 2001, 94, 111–120. [Google Scholar] [CrossRef]
  45. Wolf, Y.I.; Silas, S.; Wang, Y.; Wu, S.; Bocek, M.; Kazlauskas, D.; Krupovic, M.; Fire, A.; Dolja, V.V.; Koonin, E.V. Doubling of the known set of RNA viruses by metagenomic analysis of an aquatic virome. Nat. Microbiol. 2020, 5, 1262–1270. [Google Scholar] [CrossRef]
  46. Charlebois, R.L.; Sathiamoorthy, S.; Logvinoff, C.; Gisonni-Lex, L.; Mallet, L.; Ng, S.H.S. Sensitivity and breadth of detection of high-throughput sequencing for adventitious virus detection. NPJ Vaccines 2020, 5, 61. [Google Scholar] [CrossRef]
  47. Traniello, I.M.; Bukhari, S.A.; Kevill, J.; Ahmed, A.C.; Hamilton, A.R.; Naeger, N.L.; Schroeder, D.C.; Robinson, G.E. Meta-analysis of honey bee neurogenomic response links Deformed wing virus type A to precocious behavioral maturation. Sci. Rep. 2020, 10, 3101. [Google Scholar] [CrossRef] [Green Version]
  48. Brettell, L.E.; Schroeder, D.C.; Martin, S.J. RNAseq of Deformed wing virus and other honey bee-associated viruses in eight insect taxa with or without Varroa infestation. Viruses 2020, 12, 1229. [Google Scholar] [CrossRef]
  49. Caesar, L.; Cibulski, S.P.; Canal, C.W.; Blochtein, B.; Sattler, A.; Haag, K.L. The virome of an endangered stingless bee suffering from annual mortality in southern Brazil. J. Gen. Virol. 2019, 100, 1153–1164. [Google Scholar] [CrossRef]
  50. Chen, G.; Wang, S.; Jia, S.; Feng, Y.; Hu, F.; Chen, Y.; Zheng, H. A new strain of virus discovered in china specific to the parasitic mite varroa destructor poses a potential threat to honey bees. Viruses 2021, 13, 679. [Google Scholar] [CrossRef]
  51. Bartomeus, I.; Stavert, J.R.; Ward, D.; Aguado, O. Historical collections as a tool for assessing the global pollination crisis. Philos. Trans. R. Soc. B Biol. Sci. 2018, 374, 20170389. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Yoshioka, A.; Mishima, Y.; Fukasawa, K. Pollinators and other flying insects inside and outside the fukushima evacuation zone. PLoS ONE 2015, 10, e0140957. [Google Scholar] [CrossRef]
  53. Singh, R.; Levitt, A.L.; Rajotte, E.G.; Holmes, E.C.; Ostiguy, N.; Vanengelsdorp, D.; Lipkin, W.I.; Depamphilis, C.W.; Toth, A.L.; Cox-Foster, D.L. RNA viruses in hymenopteran pollinators: Evidence of inter-Taxa virus transmission via pollen and potential impact on non-Apis hymenopteran species. PLoS ONE 2010, 5, e14357. [Google Scholar] [CrossRef] [PubMed]
  54. Radzevičiūtė, R.; Theodorou, P.; Husemann, M.; Japoshvili, G.; Kirkitadze, G.; Zhusupbaeva, A.; Paxton, R.J. Replication of honey bee-associated RNA viruses across multiple bee species in apple orchards of Georgia, Germany and Kyrgyzstan. J. Invertebr. Pathol. 2017, 146, 14–23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Manley, R.; Boots, M.; Wilfert, L. Emerging viral disease risk to pollinating insects: Ecological, evolutionary and anthropogenic factors. J. Appl. Ecol. 2015, 52, 331–340. [Google Scholar] [CrossRef]
  56. Gisder, S.; Genersch, E. Viruses of commercialized insect pollinators. J. Invertebr. Pathol. 2017, 147, 51–59. [Google Scholar] [CrossRef] [PubMed]
  57. Dalmon, A.; Gayral, P.; Decante, D.; Klopp, C.; Bigot, D.; Thomasson, M.; Herniou, E.A.; Alaux, C.; Le Conte, Y. Viruses in the invasive hornet Vespa velutina. Viruses 2019, 11, 1041. [Google Scholar] [CrossRef] [Green Version]
  58. Schoonvaere, K.; De Smet, L.; Smagghe, G.; Vierstraete, A.; Braeckman, B.P.; de Graaf, D.C. Unbiased RNA shotgun metagenomics in social and solitary wild bees detects associations with eukaryote parasites and new viruses. PLoS ONE 2016, 11, e0168456. [Google Scholar] [CrossRef]
  59. Vasilakis, N.; Castro-Llanos, F.; Widen, S.G.; Aguilar, P.V.; Guzman, H.; Guevara, C.; Fernandez, R.; Auguste, A.J.; Wood, T.G.; Popov, V.; et al. Arboretum and Puerto Almendras viruses: Two novel rhabdoviruses isolated from mosquitoes in Peru. J. Gen. Virol. 2014, 95 Pt 4, 787–792. [Google Scholar] [CrossRef]
  60. Chen, Y.P.; Siede, R. Honey bee viruses. In Advances in Virus Research; Elsevier: Amsterdam, The Netherlands, 2007; Volume 70, pp. 33–80. [Google Scholar]
  61. McMenamin, A.J.; Flenniken, M.L. Recently identified bee viruses and their impact on bee pollinators. Curr. Opin. Insect Sci. 2018, 26, 120–129. [Google Scholar] [CrossRef]
  62. Adams, M.J.; Antoniw, J.F.; Kreuze, J. Virgaviridae: A new family of rod-shaped plant viruses. Arch. Virol. 2009, 154, 1967–1972. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Levitt, A.L.; Singh, R.; Cox-Foster, D.L.; Rajotte, E.; Hoover, K.; Ostiguy, N.; Holmes, E.C. Cross-species transmission of honey bee viruses in associated arthropods. Virus Res. 2013, 176, 232–240. [Google Scholar] [CrossRef] [PubMed]
  64. Martin, S.J.; Brettell, L.E. Deformed wing virus in honeybees and other insects. Annu. Rev. Virol. 2019, 6, 49–69. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Viruses 14 00227 g001 550
Figure 1. Genome structures of novel viruses. Conservative domains (represented by brown rectangles) were annotated by CDD, while ORFs (represented by gray rectangles) were annotated based on the results of ORFfinder and the structure of the most closely related viral genome. The Xiangshan tymo-like virus may have an unassembled ORF encode capsid protein, which is represented by a dotted line.
Figure 1. Genome structures of novel viruses. Conservative domains (represented by brown rectangles) were annotated by CDD, while ORFs (represented by gray rectangles) were annotated based on the results of ORFfinder and the structure of the most closely related viral genome. The Xiangshan tymo-like virus may have an unassembled ORF encode capsid protein, which is represented by a dotted line.
Viruses 14 00227 g001
Viruses 14 00227 g002 550
Figure 2. Phylogenetic relationship of novel viruses. The phylogenetic relationship of novel viruses was inferred from conserved RdRp amino acid sequences. Related viruses only include those accurately classified by ICTV. RdRp-conserved sequences of novel viruses and related viruses in the trees were identified by CDD. Midpoint rooting phylogenetic trees were built using the maximum likelihood method with a bootstrap value of 1000, and bootstrap values (>80%) were shown on the branches. The amino acid substitution model was annotated below each tree. The tree of Orthomyxoviridae was based on polymerase subunit PA amino acid sequences, and the tree of narnavirus was based on complete RdRp amino acid sequences.
Figure 2. Phylogenetic relationship of novel viruses. The phylogenetic relationship of novel viruses was inferred from conserved RdRp amino acid sequences. Related viruses only include those accurately classified by ICTV. RdRp-conserved sequences of novel viruses and related viruses in the trees were identified by CDD. Midpoint rooting phylogenetic trees were built using the maximum likelihood method with a bootstrap value of 1000, and bootstrap values (>80%) were shown on the branches. The amino acid substitution model was annotated below each tree. The tree of Orthomyxoviridae was based on polymerase subunit PA amino acid sequences, and the tree of narnavirus was based on complete RdRp amino acid sequences.
Viruses 14 00227 g002
Viruses 14 00227 g003 550
Figure 3. (A) Detection of each segment of segmented RNA viruses. (a) Xiangshan tombus-like virus segments 1 and 2 correspond to lanes 1 and 2. (b) Xiangshan insect virus segments 1 and 2 correspond to lanes 1 and 2. (c) Xiangshan orthomyxo-like virus segments 1–6 correspond to lanes 1–6. All PCR products were further verified by DNA sequencing. (B) Detection of complementary strands of several novel viruses. XOLV (lanes 3, 4), XTLV (lanes 5, 6), XPLV1 (lanes 7, 8), XPLV2 (lanes 9, 10), XSLV (lanes 11, 12), XIV (lanes 13, 14), and XPLV4 (lanes 15, 16). Lanes 3, 5, 7, 9, 11, 13, and 15 were the respective controls. DNA sequencing of PCR products of lanes 4, 8, 12, and 16 shows these bands in the specific amplification of virus sequences.
Figure 3. (A) Detection of each segment of segmented RNA viruses. (a) Xiangshan tombus-like virus segments 1 and 2 correspond to lanes 1 and 2. (b) Xiangshan insect virus segments 1 and 2 correspond to lanes 1 and 2. (c) Xiangshan orthomyxo-like virus segments 1–6 correspond to lanes 1–6. All PCR products were further verified by DNA sequencing. (B) Detection of complementary strands of several novel viruses. XOLV (lanes 3, 4), XTLV (lanes 5, 6), XPLV1 (lanes 7, 8), XPLV2 (lanes 9, 10), XSLV (lanes 11, 12), XIV (lanes 13, 14), and XPLV4 (lanes 15, 16). Lanes 3, 5, 7, 9, 11, 13, and 15 were the respective controls. DNA sequencing of PCR products of lanes 4, 8, 12, and 16 shows these bands in the specific amplification of virus sequences.
Viruses 14 00227 g003
Table
Table 1. Identified known viruses in wild pollinators.
Table 1. Identified known viruses in wild pollinators.
Virus NamePoolContig LengthRead NumberRead CoverageClosest Nucleotide AccessionQuery Coverage (%)Nucleotide Identity (%)
DWVA894833,345100AB070959.110096.65
B10,08658,750100MH267695.110096.49
Mayfield virus 1A902413,680100MH614304.19694.94
B7668554100MH614304.19694.4
Vespa velutina-associated acypi-like virusA9919249,864100MN565043.19894.64
B987791,996100MN565043.19895
Scaldis River bee virusA7842694100KY053857.19474.88
Hubei Diptera virus 6 RNA1A64068923100KX884805.19992.54
Hubei Diptera virus 6 RNA2A2089251096.7KX884806.19692.32
Arboretum almendravirusA11,49230,383100KC994644.110077.56
ABPVB94812598100MN565031.19996.78
CBPV RNA1B3592332100KX168412.110097.88
CBPV RNA2B2349228100MF175174.19598.44
Read number and mapping coverage were obtained using Bowtie 2, and mapped sequence reads were used to assemble contigs. The online BLASTn tool was used to compare the assembled viral contig sequences with sequences in the database. Closest nucleotide accessions were came from GenBank.
Table
Table 2. Identified novel viruses in wild pollinators.
Table 2. Identified novel viruses in wild pollinators.
Virus NamePoolFamilyGenusContig Size (bp)Query CoverageSubject AccessionClosest Relative (RdRp Amino Acid Identity)
Xiangshan martelli-like virus 1AKitaviridaeUnclassified884898.82QOJ43136Sandewavirus dungfly (67.1)
Xiangshan martelli-like virus 2AKitaviridaeUnclassified9205100QTW97796Riboviria sp. (51.8)
Xiangshan martelli-like virus 3BKitaviridaeUnclassified944054.5YP_009337423Hubei virga-like virus 1 (54.5)
Xiangshan martelli-like virus 4BKitaviridaeUnclassified11,26999.81YP_009337693Hubei virga-like virus 15 (52.1)
Xiangshan orthomyxo-like virusAOrthomyxoviridaeUnclassified2336 + 2233 + 1994 + 1604 + 1443 + 96799.84QOQ34681Dhori thogotovirus (subunit PA: 36.2)
Xiangshan tombus-like virusAUnclassifiedUnclassified2261 + 1722100QED21532Cushing virus (64.7)
Xiangshan insect virusBUnclassifiedUnclassified3665 + 2045100YP_009011225Anopheline-associated C virus (49.7)
Xiangshan sinhali-like virusBSinhaliviridaeUnclassified5770100ASS83305Lake Sinai virus (49.7)
Xiangshan picorna-like virus 1AIflaviridaeIflavirus937999.46QQX28927Soybean thrips ifla-like virus 10 (60.7)
Xiangshan picorna-like virus 2AIflaviridaeIflavirus10,05098.92QKW94218PNG bee virus 13 (67.8)
Xiangshan picorna-like virus 3BIflaviridaeIflavirus908797.28QPI13031Iflavirus IricIV-4 (50.1)
Xiangshan picorna-like virus 4BIflaviridaeIflavirus10,08799.47AWK77848Darwin bee virus 3 (56.3)
Xiangshan picorna-like virus 5BIflaviridaeIflavirus9918100YP_009337760Hubei odonate virus 4 (69.0)
Xiangshan picorna-like virus 6APolycipiviridaeSopolycivirus11,38799.28AXA52568Linepithema humile polycipivirus 1 (51.2)
Xiangshan picorna-like virus 7AUnclassifiedUnclassified10,20085.95QIT20099Diabrotica undecimpunctata virus 1 (38.5)
Xiangshan rhabdo-like virus 1ARhabdoviridaeUnclassified10,79898.74YP_009094394Puerto Almendras virus (47.1)
Xiangshan rhabdo-like virus 2ARhabdoviridaeUnclassified13,24192.08AJG39108Jingshan fly virus 2 (34.8)
Xiangshan rhabdo-like virus 3ARhabdoviridaeAlphahymrhavirus11,77499.06QMP82144Hymenopteran rhabdo-related virus OKIAV109 (51.5)
Xiangshan rhabdo-like virus 4BRhabdoviridaeBetahymrhavirus12,64599.77QPB73983Hymenopteran rhabdo-related virus OKIAV24 (58.9)
Xiangshan rhabdo-like virus 5BRhabdoviridaeUnclassified11,51799.11QMP82217Lepidopteran rhabdo-related virus OKIAV3 (47.8)
Xiangshan nyami-like virusBNyamiviridaeFormivirus970799.15QPB73978Hymenopteran orino-related virus OKIAV85 (63.5)
Xiangshan narna-like virusAUnclassifiedUnclassified319183.7APG77263Wenling narna-like virus 8 (35.3)
Xiangshan tymo-like virusBTymoviridaeMarafivirus9087100QQG34658Nasturtium officinale macula-like virus 1 (73.9)
Xiangshan toli-like virusBUnclassifiedUnclassified600990.87YP_009143313Camponotus yamaokai virus (57.3)
Xiangshan flavi-like virus 1BFlaviviridaeUnclassified14,732100QTJ63564Hymenopteran flavi-related virus (54.9)
Xiangshan flavi-like virus 2BFlaviviridaeUnclassified17,417100QTJ63570Hymenopteran flavi-related virus (58.1)
RdRp regions of novel viruses were identified using CDD, and the BLASTp online tool was used to search for the closest relative virus. The closest relative of Xiangshan orthomyxo-like virus and Xiangshan narna-like virus was searched using polymerase subunit PA aa and RdRp aa, respectively. Only the results supported both by evolutionary tree and pairwise sequence comparison were shown in classification information. Subject accessions came from GenBank. “+” in the column “Contig size” means the separator of each segment size of segmented RNA viruses. Numbers in the brackets in the column “Closest relative” mean the identity between the RdRp region of novel virus and subject sequence.
Table
Table 3. Complementary strand detection of seven novel viruses.
Table 3. Complementary strand detection of seven novel viruses.
Virus NameGenome TypeHostComplementary Strand Detection
Xiangshan orthomyxo-like virusssRNA (−)Eristalisi tenax (L.)+
Xiangshan tombus-like virusssRNA (+)Sphaerophoria indiana Bigot −
Xiangshan picorna-like virus 2ssRNA (+)Amegilla zonata L.−
Xiangshan insect virusssRNA (+)Vespidae−
Xiangshan sinhali-like virusssRNA (+)Sphecidae+
Xiangshan picorna-like virus 1ssRNA (+)Tabanidae+
Xiangshan picorna-like virus 4ssRNA (+)Scolia sinensis Saussure et Siehel +
Positive results of the reverse complementary strand are denoted by ‘+’, and negative results are denoted by ‘−’.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, N.; Huang, Y.; Li, W.; Xu, S. Virome Analysis Reveals Diverse and Divergent RNA Viruses in Wild Insect Pollinators in Beijing, China. Viruses 2022, 14, 227. https://doi.org/10.3390/v14020227

AMA Style

Li N, Huang Y, Li W, Xu S. Virome Analysis Reveals Diverse and Divergent RNA Viruses in Wild Insect Pollinators in Beijing, China. Viruses. 2022; 14(2):227. https://doi.org/10.3390/v14020227

Chicago/Turabian Style

Li, Nannan, Yizhao Huang, Wei Li, and Shufa Xu. 2022. "Virome Analysis Reveals Diverse and Divergent RNA Viruses in Wild Insect Pollinators in Beijing, China" Viruses 14, no. 2: 227. https://doi.org/10.3390/v14020227

APA Style

Li, N., Huang, Y., Li, W., & Xu, S. (2022). Virome Analysis Reveals Diverse and Divergent RNA Viruses in Wild Insect Pollinators in Beijing, China. Viruses, 14(2), 227. https://doi.org/10.3390/v14020227

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop