Next Article in Journal
When Dendritic Cells Go Viral: The Role of Siglec-1 in Host Defense and Dissemination of Enveloped Viruses
Next Article in Special Issue
Tracing Baculovirus AcMNPV Infection Using a Real-Time Method Based on ANCHORTM DNA Labeling Technology
Previous Article in Journal
TREX1 531C>T Polymorphism is Associated with High Proviral Load Levels in HTLV-1-Infected Persons
Previous Article in Special Issue
Genome Analysis of a Novel Clade II.b Alphabaculovirus Obtained from Artaxa digramma
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 1
10.3390/v12010006
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

Occurrence and Molecular Phylogeny of Honey Bee Viruses in Vespids

by
Sa Yang
1,2,†,
Philippe Gayral
3,†,
Hongxia Zhao
4,†,
Yaojun Wu
5,†,
Xuejian Jiang
5,†,
Yanyan Wu
1,2,
Diane Bigot
3,‡,
Xinling Wang
1,2,
Dahe Yang
1,2,
Elisabeth A. Herniou
3,
Shuai Deng
1,2,
Fei Li
1,2,
Qingyun Diao
1,2,
Eric Darrouzet
3,*,† and
Chunsheng Hou
1,2,*,†
1
Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing 100093, China
2
Key Laboratory of Pollinating Insect Biology, Ministry of Agriculture, Beijing 100093, China
3
Institut de Recherche sur la Biologie de l’Insecte, UMR 7261, CNRS—Université de Tours, F-37200 Tours, France
4
Guangdong Key Laboratory of Animal Conservation and Resource Utilization, Guangdong Public Laboratory of Wild Animal Conservation and Utilization, Guangdong Institute of Applied Biological Resources, Guangzhou 510260, China
5
Institute of Forestry Protection, Guangxi Zhuang Autonomous Region Forestry Research Institute, Nanning 530002, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
‡
Present address: General Zoology, Institute for Biology, Martin Luther University Halle-Wittenberg, Hoher Weg 8, 06120 Halle (Saale), Germany.
Viruses 2020, 12(1), 6; https://doi.org/10.3390/v12010006
Submission received: 21 November 2019 / Revised: 15 December 2019 / Accepted: 16 December 2019 / Published: 19 December 2019
(This article belongs to the Special Issue Insect Viruses and Pest Management)
Browse Figures
Versions Notes

Abstract

:
Since the discovery that honey bee viruses play a role in colony decline, researchers have made major breakthroughs in understanding viral pathology and infection processes in honey bees. Work on virus transmission patterns and virus vectors, such as the mite Varroa destructor, has prompted intense efforts to manage honey bee health. However, little is known about the occurrence of honey bee viruses in bee predators, such as vespids. In this study, we characterized the occurrence of 11 honey bee viruses in five vespid species and one wasp from four provinces in China and two vespid species from four locations in France. The results showed that all the species from China carried certain honey bee viruses, notably Apis mellifera filamentous virus (AmFV), Deformed wing virus (DWV), and Israeli acute paralysis virus (IAPV); furthermore, in some vespid colonies, more than three different viruses were identified. In France, DWV was the most common virus; Sacbrood virus (SBV) and Black queen cell virus (BQCV) were observed in one and two samples, respectively. Phylogenetic analyses of IAPV and BQCV sequences indicated that most of the IAPV sequences belonged to a single group, while the BQCV sequences belonged to several groups. Additionally, our study is the first to detect Lake Sinai virus (LSV) in a hornet from China. Our findings can guide further research into the origin and transmission of honey bee viruses in Vespidae, a taxon of ecological, and potentially epidemiological, relevance.

1. Introduction

Honey bees (Apis mellifera) make vital contributions to plant diversity and agricultural ecosystems [1], but they are also attacked by a variety of pathogens, including viruses, bacteria, and fungi [2]. Researchers have identified more than 20 honey bee viruses that impact bee health and survival [3]. Among them are eight common viruses thought to pose a lethal threat: Israeli acute paralysis virus (IAPV), Deformed wing virus (DWV), Black queen cell virus (BQCV), Sacbrood virus (SBV), Chronic bee paralysis virus (CBPV), Acute bee paralysis virus (ABPV), Kashmir bee virus (KBV), and Kakugo virus (DWV-A/KV) [2]. In some cases, a vector influences infection and transmission patterns. For example, the parasitic mite Varroa destructor boosts viral titers and expands the transmission range of DWV [4] and IAPV [5].
Virus transmission dynamics could also be affected by honey bee predators, a subject about which little is currently known. Hornets, for example, are major bee predators, and the yellow-legged hornet (Vespa velutina nigrithorax)—an invasive species in Europe [6]—has a major impact on honey bee populations. At present, public authorities and the beekeeping industry are deeply concerned about the rapid spread of V. velutina and its dramatic impact on honey bee survival in Europe [7], South Korea [8,9], and Japan [10]. The yellow-legged hornet was first observed in France in 2005 [11,12]. It quickly spread throughout the country and then moved into other parts of Europe [6,13,14,15]. The species is native to China, where it has also recently begun to significantly threaten honey bee survival [16], especially in the autumn and winter, when pollinating insects face a lack of food resources.
Apiculture is a major industry in China. Beekeepers manage more than 9 million hives and supply enough honey for the country’s consumers. However, the European honey bee (A. mellifera) lacks effective defensives against predators like the yellow-legged hornet. Beekeepers in Europe have reported that V. velutina preys on honey bee colonies—impacting populations [17]—but quantitative data of this type remain scarce for China (Figure S1). Since viruses can move from primary hosts to new hosts, infected hornets could influence the host dynamics of honey bee viruses. Hornets hunt across broad areas for bees, which means that hornets may come in contact with and carry viruses from multiple honey bee colonies. Several studies have shown that honey bee viruses occur in the yellow-legged hornet [18,19,20], but little research has been dedicated to exploring virus diversity. This study aims to provide a better foundation for understanding honey bee virus dynamics in vespids.
Vespids of various species were collected from four provinces in China and four locations in France to (1) characterize the occurrence of major honey bee viruses in vespids and (2) explore the phylogenetic relationships among virus sequences from vespids and honey bees.

2. Materials and Methods

2.1. Samples

Vespids and wasps were sampled in 2016 and 2017. In China, 12 colonies were collected in four provinces, including five colonies of Vespa bicolor, two colonies of Polistes rotheyi, two colonies of Vespa velutina, one colony of Vespa koreansis, one colony of Vespa affinia, one colony of Vespa sp. (Figure 1a). For each colony, fifty workers and pupa (if present) were randomly selected and tested for the presence of honey bee viruses. In France, 14 vespids colonies were collected from four locations near the city of Tours, which is in the administrative department of Indre-et-Loire (Figure 1b). The vespids were sampled live, directly transported to the laboratory in a small iron wire cage, quickly freeze killed at −80 °C, and kept frozen until analyses could take place. None of the vespids displayed symptoms of infection.

2.2. Identification of Vespid Species

We identified the vespid species collected in China using a DNA barcoding approach. Total genomic DNA was extracted from one individual per colony using the Genomic DNA Kit (CW Biotech, Beijing, China). A universal primer set (LCO-1490 and HCO-2198) was used to amplify a 709 bp fragment of cytochrome c oxidase subunit I (CO I) via polymerase chain reaction (PCR) [21]. The same CO I primers were used to carry out PCR on cDNA to verify that RNA from each colony was of good quality (Figure S2). The PCR conditions for amplifying the CO I gene were as follows: an initial denaturation step at 95 °C for 5 min; 40 cycles of 94 °C for 40 s, 54 °C for 40 s, and 72 °C for 1 min; and a final extension step at 72 °C for 10 min. PCR products were checked using electrophoresis (1.2% agarose gel). Successfully amplified PCR products were purified using a DNA Purification Kit (CW Biotech, Beijing, China); they were then Sanger sequenced employing primers described by Yue et al. [21] (Sango Biotech, Shanghai, China). The sequences were then assembled using Vector NTI Advance software and identified via a BLAST search in GenBank. Vespid species and castes were also identified using morphological characters (as described by Matsuura [22]), and these identifications were confirmed by an expert (Li Yun, Institute of Apicultural Research, CAAS). French vespid species were identified morphologically by a different expert (ED).

2.3. Extraction of DNA and RNA for Virus Identification

For the Chinese samples, DNA and RNA were extracted from a pool of approximately 20 workers or pupae per colony; queens were also used, if present. The samples were prepared using cold grinding (liquid nitrogen and a mortar/pestle). Total RNA was extracted using a TRIzol Kit (Invitrogen, Carlsbad, CA, USA) and TissuePrep homogenizer (Gering Scientific Instruments, Beijing, China) in accordance with the manufacturers’ instructions. The resulting sample was then eluted in 20 µL of sterile water and stored at −80 °C. The quantity and purity of the RNA were measured using a Nanodrop spectrophotometer (Thermo Scientific, Beijing, China). We synthesized cDNA from 4 µg of the total RNA using M-MLV reverse transcriptase (Promega, Madison, WI, USA) and the oligo dT primer in accordance with the manufacturer’s instructions. Total DNA was extracted using a Wizard Genomic DNA Purification Kit (Promega, Madison, WI, USA) and TissuePrep homogenizer (Gering Scientific Instruments, Beijing, China) following the manufacturers’ instructions. The resulting sample was then eluted into 20 μL of nuclease-free water and stored at −20 °C.
For the French samples, RNA extraction was performed as previously described [23]. Briefly, the head, abdomen, and thorax of individual vespids were placed in three separate tubes with a 5 mm stainless steel bead. The tissues were then mechanically homogenized for 90 sec (3×) at 30 Hz using a TissueLyser II (Qiagen, Courtaboeuf, France). RA1 lysis buffer (NucleoSpin RNA Isolation Kit, Qiagen, Courtaboeuf, France) was then added in accordance with the manufacturer’s instructions; the amount depended on the sample mass. The lysates of the three body sections were then pooled and centrifuged for 3 min at 14,000× g. Next, total RNA extraction was performed using 350 µL of the supernatant and a NucleoSpin RNA Isolation Kit (Macherey-Nagel, Bethlehem, PA, USA) in accordance with the manufacturer’s instructions. The resulting sample was eluted into 40 µL of RNAse free water; this process was repeated three times with the same elute to increase final RNA concentration. RNA quality was visually controlled via electrophoresis (1% agarose gel), and RNA quantity was measured using a Qubit fluorometer. Reverse transcription was performed using a RevertAid First Strand cDNA Synthesis Kit (Life Technologies, Carlsbad, CA, USA) and random hexamer primers in accordance with the manufacturer’s instructions; in most cases, 1 µg of total RNA was employed.

2.4. Identification of Honey Bee Viruses via PCR

The primer sequences are provided in Table S1. We used PCR to test the Chinese samples for the presence of 10 RNA viruses (IAPV, SBV, DWV-A, DWV-A/KV, CBPV, ABPV, BQCV, KBV, Deformed wing virus-B, Chinese Sacbrood virus, Aphid lethal paralysis virus, Lake Sinai virus, Solenopsis invicta virus, SINV) and 1 DNA virus, AmFV. The PCR reaction mixture contained 2×GoTaq reaction buffer (Promega, USA), 1 µM of the sense and antisense primers, and 1 µL of cDNA; the total reaction volume was increased to 20 µL using nuclease-free water. The PCR conditions were as follows: an initial denaturation step at 95 °C for 1 min; 33 cycles of 94 °C for 30 s, the suitable Tm for 30 s, and 72 °C for 1 min; and a final extension step at 72 °C for 10 min. The PCR products were separated using electrophoresis (1.2% agarose gel stained with GV II [Biomec, China]) and visualized using a FR-200A luminescent and fluorescent biological image analysis system (Furi, China). Product size was determined using a 100 bp molecular size marker. Because DWV-A/KV and AmFV have never been reported in France, we only looked for these two viruses in the Chinese samples.
We used multiplex PCR (as previously described [24]) to test the French samples for the presence of the six most common honey bee viruses: ABPV, IAPV, SBV, DWV, BQCV, and CBPV. The PCR reaction mixture contained 0.5 U of Diamond Taq® polymerase (Invitrogen, Milan, Italy), 2.5 µL of 10× buffer (Invitrogen, Milan, Italy), 1.5 mM of MgCl2, 0.2 mM of each dNTP, 0.4 µM of the forward and reverse primers, and 100 ng of cDNA; the total reaction volume was increased to 25 µL using nuclease-free water. The PCR conditions were as follows: an initial denaturation step at 95 °C for 2 min; four cycles of 95 °C for 1 min, 54 °C for 1 min, and 72 °C for 1 min; 30 cycles of 95 °C for 1 min, 49 °C for 1 min, and 72 °C for 1 min; and a final extension step at 72 °C for 10 min. The PCR products were separated using electrophoresis (2.5% agarose gel stained with Red DNA Stain) and visualized using a UV light.

2.5. Sequencing and Phylogenetic Analyses

The PCR products were purified and sequenced by Sango Biotech (Shanghai, China) and by GATC Biotech (Konstanz, Germany). The sequences from the Chinese samples were assembled and manually corrected using Lasergene 8 and Geneious (v. 9) software.
For each virus, sequences were then aligned using the MAFFT algorithm [25] implemented in Geneious (v. 9) [26] and all the available homologous sequences found in the NCBI Nucleotide database. Maximum likelihood phylogenies were built using PhyML 3.0 online (http://www.atgc-montpellier.fr/phyml/) [27]. In each case, the best evolutionary model was determined using smart model selection (SMS) and the Akaike information criterion (AIC) [28]. Node support values were estimated using an approximate likelihood-ratio test (aLRT) [29].
The sequences of the viruses were submitted to GenBank under the accession numbers MF092814 to MF092832 for BQCV from China and MH133351 and MH133352 for BQCV from France; MF092819 to MF092825 for IAPV from China; MF092827 to MF092832 for DWV-A from China; MF092817 for AmFV from China; MF092826 for DWV-B from China; MF092818 for DWV-A/KV from China; and MH133361 for SBV from France. The DWV-A sequences VC-F1, VV-L2, VC-S2 were shorter than 200 bp long and were thus not attributed accession numbers.

3. Results

3.1. Virus Occurrence

We screened the vespid colonies for the presence of 11 viruses (Table 1). Several viruses were detected in all the vespid species sampled in China, with the exception of those from the Jiangxi and Qinghai provinces. Four honey bee viruses (CBPV, KBV, ALPV, and SINV) were not found in any of the vespids. IAPV and DWV-B were the most common viruses: they were present in nine individuals, representing five vespid species, and all the provinces sampled (Table 1). The third most common virus was AmFV, which was detected in six individuals, followed by BQCV and DWV-A/KV, which were detected in five individuals. LSV was found in pupae from the B2 colony. Interestingly, the honey bee viruses were most commonly found in pupae, not adults (Table 1). In the vespids sampled in France, the most common virus was DWV-A, followed by BQCV and SBV. Finally, ABPV was not observed in any vespids from either country.
In China, the co-occurrence of more than one virus in vespid colonies was fairly common for both adults and pupae (Table 2). In adults, all co-occurrences involved DWV. More co-occurrences were found in pupae than in adults, with up to six viruses detected in a single colony. The main viruses involved were AmFV, DWV, and IAPV. No co-occurrences were detected in French vespids.

3.2. Phylogenetic Relationships Among the Virus Sequences Found in Vespids

Our phylogenetic analyses confirmed the species identity of the new AmFV sequence from China (B4/Beijing, Polistes rothneyi) (node support value = 0.95) (Figure 2a). The BQCV sequences from China and France formed different groups (Figure 2b). The French sequences clustered with European BQCV sequences (mainly from Turkey; node support value = 0.83), and the three Chinese BQCV sequences were found within a clade formed by Asian and US sequences (node support values = 0.84 and 0.82, respectively) (Figure 2b).
In the phylogeny for IAPV, six of the seven IAPV sequences from this study formed a subgroup (node support value = 0.88), and the sister group of this clade was the cluster of the IAPV sequences found in A. mellifera in South Korea (Figure 2c). The seventh sequence was most closely related to IAPV sequences in A. mellifera in China (node support value = 0.91). One SBV sequence was discovered in a hornet in France, and it was located within a clade formed by SBV sequences in A. mellifera in Australia and South Korea (node support value = 0.75) (Figure 2d).
The phylogeny shows that, as expected, the DWV complex comprises two genotypes, A and B (formerly V. destructor virus-1). Interestingly, both DWV genotypes were represented in the sequences obtained from vespids in China and France. Because different DWV loci were amplified for the Chinese and French samples, separate phylogenetic trees were reconstructed for each country (Figure S3), and we generated a combined phylogenetic tree using concatenated aligned sequences (Figure 3). Both approaches gave similar results.

4. Discussion

Vespids are widely distributed across the world, from tropical to temperate regions. While some species can help disperse plant seeds, such as those of Stemona tuberosa [30], other species are predators. For example, the yellow-legged hornet (V. velutina) is a specialized predator of honey bee foragers and thus represents a serious threat to the apiculture industry and pollination ecosystem services [31]. This species was introduced into Europe in 2004 and since then it has placed significant biological stress on honey bee colonies.
Our study aimed to take a first look at the occurrence of 11 honey bee viruses in vespid species collected in China and France; it also examined phylogenetic relationships among virus sequences. The results show that several viruses were commonly present and broadly geographically distributed in a variety of vespid species. In addition, we found more viruses in pupae than in adults, which is not surprising given that adults feed larvae large quantities of masticated bees. While we did not investigate infection dynamics, this finding might suggest that pupae are more suitable or more susceptible hosts than adults, if the viruses are capable of replicating in vespids. In addition, we also did not detect the pupae from France vespids due to limited samples. Such stage-dependent susceptibility has been described for several viruses that infect both A. mellifera and A. cerana [32,33]. Further experimental research is now needed to test whether (1) honey bee viruses can replicate in hornets (2) hornets can then transmit the viruses to honey bees and (3) pupae from France vespids carry more honey bee viruses than adults.
Based on our phylogenetic analyses, the viruses infecting vespids do not form specific clades. This result could suggest that the honey bee viruses are not transmitted among vespids but rather from honey bees to vespids. It might also indicate that, if replication occurs in vespids, it primarily results from recent spillovers from honey bees. A recent study showed that such spillovers could result in DWV-B and KBV infections in the yellow-legged hornet, including a case in which an individual was symptomatic (i.e., had deformed wings) [34]. In our study, DWV was the most common virus in both China and France, which is unsurprising given its frequency in honey bees worldwide [4]. Both genotypes—DWV-A and DWV-B—were observed. Future research might reveal more specific replication dynamics in vespids.
IAPV was first identified in Israel in 2004 but has since been found to have a worldwide distribution [35]. It is a virus that displays high genetic diversity [36]. IAPV is quite common in China [37], and its replicative form has previously been detected in the yellow-legged hornet [18]. Consequently, there may be interspecific transmission of IAPV, an issue that should be explored further in epidemiological studies. In contrast, BQCV sequences from vespids in China and France formed two different groups. The long branches in the phylogeny suggest that this virus displays pronounced genetic diversity. As a result, concerted efforts should be made to construct a more global phylogeny for BQCV to better identify its possible adaptive diversity and host associations.
Although AmFV was first identified in the honey bee in 1963 [38] and its complete dsDNA genome was recently published [39], we actually know little about the virus’s occurrence [40], especially in China. It is thus quite notable that we observed AmFV in several hornets and especially in pupae. There could be several explanations for this result. First, AmFV might be more common in China than in Europe, where it was first detected. Second, honey bees might be frequently infected with AmFV, but we could be unaware of this fact because there is no specific monitoring of this virus. If AmFV prevalence is high in honey bees, vespids—and notably vespid larvae/pupae—may be more likely to ingest and accumulate the virus, which might explain why AmFV is so common in vespids in China. While AmFV could also be infectious in vespids, our phylogenetic analyses revealed no evidence that AmFV sequences from vespids were distinct from AmFV sequences from honey bees, which could suggest that vespids do not carry a specific form of the virus. Complete genome sequencing might reveal more divergent molecular markers that would allow us to better discriminate among virus sequences. Overall, our findings for AmFV highlight the need to further investigate the dynamics of this virus in honey bees and vespids.

5. Conclusions

To summarize, we characterized the occurrence of the 11 most important honey bee viruses in vespid species in China and France. We also identified differences in virus co-occurrence among different vespid species in China. All the virus sequences that we obtained from vespids were closely related to sequences obtained from honey bees, likely indicating that the vespids acquired the viruses while feeding on bees. However, further research is needed to explore infection dynamics in vespids, and the potential for virus transmission between vespids and honey bees.

Supplementary Materials

The supplementary materials are available online at https://www.mdpi.com/1999-4915/12/1/6/s1. They include Table S1. PCR primers used to identify the honey bee viruses in vespids collected in China and France. Figure S1. Picture of a yellow-legged hornet (Vespa velutina) attacking a honey bee. Figure S2. The quality of the RNA extracted from the vespid samples. M, DNA ladder; the numbers 1–10 refer to the colony IDs B1, B3, B4, B5, Q1, GD1, GD2, GD3, GD4, and JX, respectively. Figure S3. Separate maximum likelihood phylogenetic trees for each Deformed wing virus (DWV) sequence region in China and France. (a) Region 1: 2170–2438 (201 sites, HKY85+I model); (b) Region 2: 2843–3182 (340 sites, HKY85+I model); (c) Region 3: 6111–6310 (152 sites, TN93+I model); and (d) Region 4: 8377–8753 (311 sites, HKY85+G model). The country names, species names, and GenBank accession numbers were also added to each branch of the tree. The scale bars represent the substitution rate per site, and the values above the nodes are the aLRT statistics.

Author Contributions

Conceptualization, C.H. and E.D.; methodology, S.Y., H.Z., D.Y., E.A.H., P.G., X.W., Y.W. (Yanyan Wu), S.D. and F.L.; software, P.G. and D.B.; validation, C.H., E.D., E.A.H., P.G. and D.B.; formal analysis, C.H.; investigation, H.Z., Y.D., S.Y., S.D., C.H., Y.W. (Yaojun Wu), X.J.; writing and editing the manuscript, C.H., E.A.H., D.B. and E.D.; project administration, C.H. and E.D.; and funding acquisition, C.H., Q.D. and E.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant number 31811530276,31572471), the National Project for Upgrading Overall Bee-Product Quality in the Beekeeping Industry of China (1803086), The Key Technology Research and Development Program of Guangxi province (AB16380094), and the Agricultural Science and Technology Innovation Program (CAAS-ASTIP-2018-IAR).

Acknowledgments

We thank the beekeepers at all our sampling locations for helping to collect the samples.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Klein, A.-M.; Vaissiere, B.E.; Cane, J.H.; Steffan-Dewenter, I.; Cunningham, S.A.; Kremen, C. Importance of pollinators in changing landscapes for world crops. Proc. R. Soc. B Biol. Sci. 2007, 274, 303–313. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Goulson, D.; Nicholls, E.; Botías, C.; Rotheray, E.L. Bee declines driven by combined stress from parasites, pesticides, and lack of flowers. Science 2015, 347, 1255957. [Google Scholar] [CrossRef]
  3. Brutscher, L.M.; McMenamin, A.J.; Flenniken, M.L. The buzz about honey bee viruses. PLoS Pathog. 2016, 12, e1005757. [Google Scholar] [CrossRef]
  4. Wilfert, L.; Long, G.; Leggett, H.C.; Schmid-Hempel, P.; Butlin, R.; Martin, S.J.M.; Boots, M. Deformed wing virus is a recent global epidemic in honey bees driven by Varroa mites. Science 2016, 351, 594–597. [Google Scholar] [CrossRef] [Green Version]
  5. Di Prisco, G.; Pennacchio, F.; Caprio, E.; Boncristiani, H.F.; Evans, J.D.; Chen, Y.P. Varroa destructor is an effective vector of Israeli acute paralysis virus in the honey bee, Apis mellifera. J. Gen. Virol. 2011, 92, 151–155. [Google Scholar] [CrossRef]
  6. Takeuchi, T.; Takahashi, R.; Kiyoshi, T.; Nakamura, M.; Minoshima, Y.; Takahahi, J. The origin and genetic diversity of the yellow-legged hornet, Vespa velutina introduced in Japan. Insectes Sociaux 2017, 3, 313–320. [Google Scholar] [CrossRef]
  7. Perrard, A.; Haxaire, J.; Rortais, A.; Villemant, C. Observations on the colony activity of the Asian hornet Vespa velutina Lepeletier 1836 (Hymenoptera: Vespidae: Vespinae) in France. Ann. Soc. Entomol. Fr. 2009, 45, 119–127. [Google Scholar] [CrossRef] [Green Version]
  8. Kim, J.K.; Choi, M.; Moon, T.Y. Occurrence of Vespa velutina Lepeletier from Korea, and a revised key for Korean Vespa species (Hymenoptera:Vespidae). Entomol. Res. 2006, 2, 112–115. [Google Scholar] [CrossRef]
  9. Choi, M.B.; Martin, S.J.; Lee, J.W. Distribution, spread, and impact of the invasive hornet Vespa velutina in South Korea. Entomol. Res. 2012, 15, 473–477. [Google Scholar] [CrossRef]
  10. Villemant, C.; Zuccon, D.; Rome, Q.; Muller, F.; Poinar, G.O., Jr.; Justine, J.L. Can parasites halt the invader? Mermithid nematodes parasitizing the yellow-legged Asian hornet in France. PeerJ 2015, 3, e947. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Haxaire, J.; Bouguet, J.P.; Tamisier, J.P. Vespa velutina Lepeletier, 1836, une redoutable nouveauté pour la faune de France (Hym., Vespidae). Bull. Soc. Entomol. Fr. 2006, 111, 194. [Google Scholar]
  12. Villemant, C.; Barbet-Massin, M.; Perrard, A.; Muller, F.; Gargominy, O.; Jiguet, F.; Rome, Q. Predicting the invasion risk by the alien bee-hawking Yellow-legged hornet Vespa velutina nigrithorax across Europe and other continents with niche models. Biol. Conserv. 2011, 144, 2142–2150. [Google Scholar] [CrossRef]
  13. Robinet, C.; Darrouzet, E.; Suppo, C. Spread modelling: A suitable tool to explore the role of human-mediated dispersal in the range expansion of the yellow-legged hornet in Europe. Int. J. Pest Manag. 2018, 1, 1–10. [Google Scholar] [CrossRef]
  14. Rome, Q.; Dambrine, L.; Onate, C.; Muller, F.; Villemant, C.; Pérez, A.L.G.; Maia, M.; Carvalho, P.; Bruneau, E. Spread of the invasive hornet Vespa velutina Lepeletier, 1836, in Europe in 2012 (Hym., Vespidae). Bull. Soc. Entomol. Fr. 2013, 118, 21–22. [Google Scholar]
  15. Demichelis, S.; Manino, A.; Minuto, G.; Mariotti, M.; Porporato, M. Social wasp trapping in north west Italy: Comparison of different bait-traps and first detection of Vespa velutina. Bull. Insectol. 2014, 67, 307–317. [Google Scholar]
  16. Tan, K.; Wu, J.; Liu, Q.; Zhou, D.Y. Reduction of foraging activity by A.cerana colonies attacked by Vespa Velutina. J. Bee 2004, 2, 7–8. [Google Scholar]
  17. Requier, F.; Rome, Q.; Chiron, G.; Decante, D.; Marion, S.; Menard, M.; Muller, F.; Villemant, C.; Henry, M. Predation of the invasive Asian hornet affects foraging activity and survival probability of honey bees in Western Europe. J. Pest Sci. 2019, 92, 567–578. [Google Scholar] [CrossRef]
  18. Yañez, O.; Zheng, H.Q.; Hu, F.L.; Neumann, P.; Dietemann, V. A scientific note on Israeli acute paralysis virus infection of Eastern honeybee Apis cerana and vespine predator Vespa velutina. Apidologie 2012, 5, 587–589. [Google Scholar] [CrossRef] [Green Version]
  19. Garigliany, M.; Taminiau, B.; El Agrebi, N.; Cadar, D.; Gilliaux, G.; Hue, M.; Desmecht, D.; Daube, G.; Linden, A.; Farnir, F. Moku virus in invasive Asian hornets, Belgium, 2016. Emerg. Infect. Dis. 2017, 12, 2109. [Google Scholar] [CrossRef] [Green Version]
  20. Mazzei, M.; Forzan, M.; Cilia, G.; Sagona, S.; Bortolotti, L.; Felicioli, A. First detection of replicative deformed wing virus (DWV) in Vespa velutina nigrithorax. Bull. Insectol. 2018, 2, 211–216. [Google Scholar]
  21. Yue, Q.Y.; Qiu, D.Y.; Huang, Y.W.; Liu, G.X. Application of DNA barcoding to species identification of an unknown insect larva. Chin. J. Health Lab. Technol. 2011, 3, 615–617. [Google Scholar]
  22. Matsuura, M. Vespa and Provespa. In The Social Biology of Wasps; Ross, K.G., Matthews, R.W., Eds.; Cornell University Press: Ithaca, NY, USA, 1991; p. 30. [Google Scholar]
  23. Bigot, D.; Dalmon, A.; Roy, B.; Hou, C.; Germain, M.; Romary, M. The discovery of Halictivirus resolves the Sinaivirus phylogeny. J. Gen. Virol. 2017, 98, 2864–2875. [Google Scholar] [CrossRef] [PubMed]
  24. Sguazza, G.H.; Reynaldi, F.J.; Galosi, C.M.; Pecoraro, M.R. Simultaneous detection of bee viruses by multiplex PCR. J. Virol. Methods. 2013, 194, 102–106. [Google Scholar] [CrossRef] [PubMed]
  25. Katoh, K.; Misawa, K.; Kuma, K.; Miyata, T. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef] [Green Version]
  26. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef]
  27. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 3, 307–321. [Google Scholar] [CrossRef] [Green Version]
  28. Lefort, V.; Longueville, J.E.; Gascuel, O. SMS: Smart model selection in PhyML. Mol. Biol. Evol. 2017, 34, 2422–2424. [Google Scholar] [CrossRef] [Green Version]
  29. Anisimova, M.; Gascuel, O. Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst. Biol. 2006, 55, 539–552. [Google Scholar] [CrossRef]
  30. Chen, G.; Wang, Z.W.; Qin, Y.; Sun, W.B. Seed dispersal by hornets: An unusual insect-plant mutualism. J. Integr. Plant Biol. 2017, 59, 792–796. [Google Scholar] [CrossRef] [Green Version]
  31. López, S.; González, M.; Goldarazena, A. Vespa velutina Lepeletier, 1836 (Hymenoptera: Vespidae): First records in Iberian Peninsula. Eppo Bull. 2011, 41, 439–441. [Google Scholar] [CrossRef]
  32. Boncristiani, H.F.; Evans, J.D.; Chen, Y.; Pettis, J.; Murphy, C.; Lopez, D.L.; Simone-Finstrom, M.; Strand, M.; Tarpy, D.R.; Rueppell, O. In vitro infection of pupae with Israeli acute paralysis virus suggests disturbance of transcriptional homeostasis in honey bees (Apis mellifera). PLoS ONE 2013, 9, e73429. [Google Scholar] [CrossRef] [PubMed]
  33. Ai, H.; Yan, X.; Han, R. Occurrence and prevalence of seven bee viruses in Apis mellifera and Apis cerana apiaries in China. J. Invertebr. Pathol. 2012, 1, 160–164. [Google Scholar] [CrossRef] [PubMed]
  34. Dalmon, A.; Gayral, P.; Decante, D.; Klopp, C.; Bigot, D.; Thomasson, M.; Herniou, E.A.; Alaux, C.; Conte, Y.L. Viruses in the invasive hornet Vespa velutina. Viruses 2019, 11, 1041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Maori, E.; Lavi, S.; Mozes-Koch, R.; Gantman, Y.; Peretz, Y.; Edelbaum, O.; Tanne, E.; Sela, I. Isolation and characterization of Israeli acute paralysis virus, a dicistrovirus affecting honey bees in Israel: Evidence for diversity due to intra-and inter-species recombination. J. Gen. Virol. 2008, 88, 3428–3438. [Google Scholar] [CrossRef] [PubMed]
  36. Palacios, G.; Hui, J.; Quan, P.L.; Kalkstein, A.; Honkavuori, K.S.; Bussetti, A.V.; Conlan, S.; Evans, J.; Chen, Y.P.; VanEngelsdorp, D.; et al. Genetic analysis of Israel acute paralysis virus: Distinct clusters are circulating in the United States. J. Virol. 2008, 13, 6209–6217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Diao, Q.; Li, B.; Zhao, H.; Wu, Y.; Guo, R.; Dai, P.; Chen, D.; Wang, Q.; Hou, C. Enhancement of chronic bee paralysis virus levels in honeybees acute exposed to imidacloprid: A Chinese case study. Sci. Total Environ. 2018, 630, 487–494. [Google Scholar] [CrossRef] [PubMed]
  38. Clark, T.B. A filamentous virus of the honey bee. J. Invertebr. Pathol. 1978, 32, 332–340. [Google Scholar] [CrossRef]
  39. Gauthier, L.; Cornman, S.; Hartmann, U.; Cousserans, F.; Evans, J.D.; De Miranda, J.R.; Neumann, P. The Apis mellifera filamentous virus genome. Viruses 2015, 7, 3798–3815. [Google Scholar] [CrossRef] [Green Version]
  40. Hartmann, U.; Forsgren, E.; Charrière, J.D.; Neumann, P.; Gauthier, L. Dynamics of Apis mellifera filamentous virus (AmFV) infections in honey bees and relationships with other parasites. Viruses 2015, 7, 2654–2667. [Google Scholar] [CrossRef] [Green Version]
Viruses 12 00006 g001 550
Figure 1. Sampling locations in (a) China and (b) France. Sampled species are indicated by names and photos.
Figure 1. Sampling locations in (a) China and (b) France. Sampled species are indicated by names and photos.
Viruses 12 00006 g001
Viruses 12 00006 g002 550
Figure 2. Maximum Likelihood phylogenetic trees of new AmFV, BQCV, IAPV, and SBV sequences found in China and France. (a) Apis mellifera filamentous virus phylogeny (524 sites, TN93+G model), (b) Black queen cell virus phylogeny (634 sites, GTR+I+G model), (c) Israeli acute bee paralysis virus (439 sites, TN93+I model) and (d) Sacbrood virus phylogeny (277 sites, GTR+G model). Names of countries and species for each known and new sequence were added in each leaf of the tree, joined with the corresponding GenBank accession number. The scale bars represent the substitution rate per site. The values correspond to the values of nodes supported by aLRT statistics.
Figure 2. Maximum Likelihood phylogenetic trees of new AmFV, BQCV, IAPV, and SBV sequences found in China and France. (a) Apis mellifera filamentous virus phylogeny (524 sites, TN93+G model), (b) Black queen cell virus phylogeny (634 sites, GTR+I+G model), (c) Israeli acute bee paralysis virus (439 sites, TN93+I model) and (d) Sacbrood virus phylogeny (277 sites, GTR+G model). Names of countries and species for each known and new sequence were added in each leaf of the tree, joined with the corresponding GenBank accession number. The scale bars represent the substitution rate per site. The values correspond to the values of nodes supported by aLRT statistics.
Viruses 12 00006 g002
Viruses 12 00006 g003 550
Figure 3. Concatenate Maximum Likelihood phylogenetic tree of Deformed wing virus sequences found in China and France. This phylogeny was built by concatenated four sequences’ alignment regions, due to multiple PCR detection primers used in protocols (Region 1: 2170–2438; region 2: 2843–3182; Region 3: 6111–6310; Region 4: 8377–8753) (1004 sites, GTR+G model). To tree topology of the concatenation was verified using each independent region in their respective phylogenetics tree (Supplementary Figure S2a–d). Country and species names for each known and new sequence were added in each leaf of the tree, joined with the corresponding GenBank accession number. The scale bars represent the substitution rate per site. The values correspond to the values of nodes supporting by aLRT statistics.
Figure 3. Concatenate Maximum Likelihood phylogenetic tree of Deformed wing virus sequences found in China and France. This phylogeny was built by concatenated four sequences’ alignment regions, due to multiple PCR detection primers used in protocols (Region 1: 2170–2438; region 2: 2843–3182; Region 3: 6111–6310; Region 4: 8377–8753) (1004 sites, GTR+G model). To tree topology of the concatenation was verified using each independent region in their respective phylogenetics tree (Supplementary Figure S2a–d). Country and species names for each known and new sequence were added in each leaf of the tree, joined with the corresponding GenBank accession number. The scale bars represent the substitution rate per site. The values correspond to the values of nodes supporting by aLRT statistics.
Viruses 12 00006 g003
Table
Table 1. Presence of major honey bee viruses in various vespid species.
Table 1. Presence of major honey bee viruses in various vespid species.
CountryDevelopmental StageColony ID/Province (Species)AmFVIAPVDWV-ADWV-BDWV-A/KVCSBVSBVBQCVLSV
ChinaAdultB1/Beijing (Vespa bicolor)ntnt+nt+ntntntnt
Pupa+++nt+ntntntnt
AdultB2/Beijing (Vespa bicolor)ntnt+ntntntnt+nt
Pupa+++nt++nt++
AdultB3/Beijing (Vespa bicolor)ntntntntntntntntnt
Pupantntntntntntntntnt
AdultB4/Beijing (Polistes rothneyi)+ntntntntntntntnt
Pupa+ntntntntntnt+nt
AdultB5/Beijing (Polistes rothneyi)ntntntntntntntntnt
Pupa++ntntntntntntnt
AdultQ1/Qinghai (Vespa koreansis)ntntntntntntntntnt
AdultGD1/Guangdong (Vespa affinia)ntnt+ntntntntntnt
Pupant++ntntntntntnt
AdultGD2/Guangdong (Vespa velutina)ntntntntntntntntnt
Pupant+ntntntntntntnt
AdultGD3/Guangdong (Vespa velutina)ntntntntntntntntnt
Pupant+ntntntntntntnt
AdultGD4/Guangdong (Vespa sp.)ntntntntntntntntnt
Pupant++++ntnt+nt
AdultGD5/Guangdong (Vespa bicolor)nt++ntntntntntnt
Pupant++ntntntntntnt
AdultJX/Jiangxi (Vespa bicolor)ntntntntntntntntnt
Pupantntntntntntntntnt
FranceAdultFondettes (Vespa crabro)ntnt+ntntntntntnt
AdultSaché (Vespa crabro)ntnt+ntntntntntnt
AdultSaché (Vespa crabro)ntnt+ntntntntntnt
AdultLussault-sur-Loire (Vespa velutina)ntnt+ntntntntntnt
AdultLussault-sur-Loire (Vespa velutina)ntntntntntnt+ntnt
AdultFondettes (Vespa crabro)ntntntntntntnt+nt
AdultTours (Vespa crabro)ntntntntntntnt+nt
nt: not tested; only the positive samples are included in the table.
Table
Table 2. Co-occurrence of virus species in vespid colonies in China.
Table 2. Co-occurrence of virus species in vespid colonies in China.
Colony IDVespid SpeciesDevelopmental StageViruses Present
B2Vespa bicolorAdultDWV+BQCV
GD5Vespa bicolorAdultIAPV+DWV
B4Polistes rothneyiPupaAmFV+BQCV
B5Polistes rothneyiPupaAmFV+IAPV
GD1Vespa affiniaPupaDWV+IAPV
B1Vespa bicolorPupaAmFV+IAPV+DWV
GD4Vespa sp.PupaIAPV+DWV+BQCV
B2Vespa bicolorPupaAmFV+IAPV+DWV+CSBV+BQCV+LSV

Share and Cite

MDPI and ACS Style

Yang, S.; Gayral, P.; Zhao, H.; Wu, Y.; Jiang, X.; Wu, Y.; Bigot, D.; Wang, X.; Yang, D.; Herniou, E.A.; et al. Occurrence and Molecular Phylogeny of Honey Bee Viruses in Vespids. Viruses 2020, 12, 6. https://doi.org/10.3390/v12010006

AMA Style

Yang S, Gayral P, Zhao H, Wu Y, Jiang X, Wu Y, Bigot D, Wang X, Yang D, Herniou EA, et al. Occurrence and Molecular Phylogeny of Honey Bee Viruses in Vespids. Viruses. 2020; 12(1):6. https://doi.org/10.3390/v12010006

Chicago/Turabian Style

Yang, Sa, Philippe Gayral, Hongxia Zhao, Yaojun Wu, Xuejian Jiang, Yanyan Wu, Diane Bigot, Xinling Wang, Dahe Yang, Elisabeth A. Herniou, and et al. 2020. "Occurrence and Molecular Phylogeny of Honey Bee Viruses in Vespids" Viruses 12, no. 1: 6. https://doi.org/10.3390/v12010006

APA Style

Yang, S., Gayral, P., Zhao, H., Wu, Y., Jiang, X., Wu, Y., Bigot, D., Wang, X., Yang, D., Herniou, E. A., Deng, S., Li, F., Diao, Q., Darrouzet, E., & Hou, C. (2020). Occurrence and Molecular Phylogeny of Honey Bee Viruses in Vespids. Viruses, 12(1), 6. https://doi.org/10.3390/v12010006

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