# Viruses in the Invasive Hornet Vespa velutina

Anne Dalmon, Philippe Gayral, Damien Decante, Christophe C. Kropp, Diane Bigot, Maxime Thomasson, Elisabeth A. Herniou, Cédric Alaux, Yves Le Conte

► To cite this version:

Anne Dalmon, Philippe Gayral, Damien Decante, Christophe C. Kropp, Diane Bigot, et al. Viruses in the Invasive Hornet Vespa velutina. Viruses, 2019, 11 (11), pp.1041. 10.3390/v1111041. hal-02388796

HAL Id: hal-02388796

https://hal.science/hal-02388796

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Viruses in the Invasive Hornet *Vespa velutina*

Anne Dalmon<sup>1,2,*</sup>, Philippe Gayral<sup>3</sup>, Damien Decante<sup>2,4</sup>, Christophe Klopp<sup>5 ID</sup>, Diane Bigot<sup>3,†</sup>, Maxime Thomasson<sup>1</sup>, Elisabeth A Herniou<sup>3 ID</sup>, Cédric Alaux<sup>1,2</sup> and Yves Le Conte<sup>1,2 ID</sup>

Received: 31 July 2019; Accepted: 25 October 2019; Published: 8 November 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** The Asian yellow-legged hornet *Vespa velutina nigrithorax*, a major predator of honeybees, is spreading in Europe in part due to a lack of efficient control methods. In this study, as a first step to identify biological control agents, we characterized viral RNA sequences present in asymptomatic or symptomatic hornets. Among 19 detected viruses, the honey bee virus *Deformed wing virus-B* was predominant in all the samples, particularly in muscles from the symptomatic hornet, suggesting a putative cause of the deformed wing symptom. Interestingly, two new viruses closely related to *Acyrtosiphon pisum virus* and *Himetobi P virus* and viruses typically associated with honey bees, *Acute bee paralysis virus* and *Black queen cell virus*, were detected in the brain and muscles, and may correspond to the circulation and possible replication forms of these viruses in the hornet. *Aphid lethal paralysis virus*, *Bee Macula-like virus*, and *Moku virus*, which are known to infect honey bees, were also identified in the gut virus metagenome of hornets. Therefore, our study underlined the urgent need to study the host range of these newly discovered viruses in hornets to determine whether they represent a new threat for honey bees or a hope for the biocontrol of *V. velutina*.

**Keywords:** *Vespidae*; invasive species; honey bee viruses; new viruses; DWV; ABPV; BQCV; KBV; ALPV; BeeMLV

## 1. Introduction

The original distribution area of the Asian yellow-legged hornet *Vespa velutina nigrithorax*, native to Southeast Asia, ranges from Nepal to South Eastern China. However, in 2004, it was first observed near Agen, France, [1], and one year later, phylogenetically identified [2,3]. Its accidental introduction, through very few or single multi-mated gynes, most probably occurred before 2005 via boat transport of bonsai pots from the Zhejiang or Jiangsu provinces of Eastern China [4,5].

Since a single individual can found a colony, this outbreak had a high invasive potential, even though, like other social *Vespidae* under temperate climates, the colonies last only one season.

<footnote n="1">1 INRA—Unité Abeilles et Environnement, Site Agroparc, Domaine St Paul, 228, Route de l'aérodrome CS40509, 84914 Avignon Cedex 9, France; maxime.thomasson@hotmail.fr (M.T.); cedric.alaux@inra.fr (C.A.); yves.le-conte@inra.fr (Y.L.C.)</footnote>

<footnote n="2">2 UMT PRADE, Unité Mixte Technologique de la PRotation de l'Abeille Dans l'Environnement, Site Agroparc, Domaine St Paul, 228, Route de l'aérodrome CS40509, 84914 Avignon Cedex 9, France; d.decante.adapi@free.fr</footnote>

<footnote n="3">3 Institut de Recherche sur la Biologie de l'Insecte, UMR 7261, CNRS - Université de Tours, 37200 Tours, France; philippe.gayral@univ-tours.fr (P.G.); bigot.diane@gmail.Com (D.B.); Elisabeth.Herniou@Univ-Tours.Fr (E.A.H.)</footnote>

<footnote n="4">4 ITSAP, Site Agroparc, Domaine St Paul, 228, Route de l'aérodrome CS40509, 84914 Avignon Cedex 9, France</footnote>

<footnote n="5">5 INRA—Genotoul Bioinfo, Unité MIAT Mathématiques et Informatique Appliquées, 24 chemin de Borde-Rouge, 31326 Auzeville, France; christophe.klopp@inra.fr</footnote>

<footnote>* Correspondence: anne.dalmon@inra.fr</footnote>

<footnote>† Present address: General Zoology, Institute for Biology, Martin Luther University Halle-Wittenberg, Hoher Weg 8, 06120 Halle (Saale), Germany.</footnote>
Every year, the surviving yellow-legged hornet-mated queens disperse in order to initiate new nests after hibernation. Each colony, initiated by this single individual, develops until the following autumn, first producing several thousand workers, and then hundreds of males and new founders, which mate in autumn [6]. Despite the genetic bottleneck and inbreeding depression associated with this single introduction event [4,5,7], no population extinction has been observed yet. Instead, this first successful introduction of an exotic *Vespidae* in Europe resulted in a rapid colonization of France and neighbouring countries (Spain, Portugal, Italy, Belgium, Germany, England, and the Netherlands; INPN, 2019), with substantial population levels and impacts in most of the colonized regions [5,8–10]. Moreover, the niche models developed in order to determine its possible distribution area [1,11–17] converge to predict a much wider distribution in the future.

During nest development, *V. velutina* predates on a wide diversity of insects such as hoverflies and houseflies, but with a preference for aggregative species, including social Hymenoptera species such as honey bees and common wasps [18]. Although the mechanisms resulting in colony collapse remain poorly described, predictive models suggest that high *V. velutina* pressure strongly impacts the colonies [19]. To date, empirical methods, including spring foundress trapping and autumnal worker mass trapping, have been used to impair *V. velutina* establishment and dissemination [20,21]. However, the efficacy of these methods is insufficient to control *V. velutina* populations [18,22–26]. Nest destruction seems to be efficient, but despite recent technological advances for nest detection [27–29], these are cumbersome and expensive [30], thus limiting large-scale implementation.

As observed with other wasp introduction events [31], the lack of efficient control methods stems from a lack of knowledge on the biology, behaviour and pathology of *V. velutina* [6,9,21,23]. Regarding pathogens that could affect the Asian hornet, only *Deformed wing virus* [32], *Moku virus* [33] and *Israeli acute paralysis virus* (IAPV) [34] have been described in some Asian specimens. However, 11 viruses initially detected in managed Western honey bees (*Apis mellifera*) have also been identified in wild Hymenoptera species [35] including *Vespula* spp. [36–38], which belongs to the same sub-family as the *Vespa* genus, *Vespinae*. Because *V. velutina* is an important predator of *A. mellifera*, we can expect pathogen transfers from target prey to predator. Among the seven viruses commonly found in honey bee colonies, *Deformed wing virus* (DWV), *Black queen cell virus* (BQCV) and *Sacbrood virus* (SBV) have the highest prevalence, as they are detected in up to 90% of western honey bee colonies [39–41]. DWV presence has been positively correlated with colony losses [41–45]. Different genotypic variants have been described: the original first described DWV named DWV-A [46], *Varroa destructor virus 1* (VDV-1) recently renamed DWV-B, and DWV-C [47]. Acute paralysis viruses from the “AKI” complex (*Acute bee paralysis virus* (ABPV), *Kashmir bee virus* (KBV), *Israeli acute bee paralysis* (IAPV)) have a worldwide distribution. While ABPV is more prevalent in Europe and South America, KBV prevails in North America and South America and IAPV in the Middle East and Australia [48]. *Chronic bee paralysis virus* (CBPV) is widespread, although not highly prevalent in western honey bees [49], but has been detected in other social insects such as ants [50]. Such widespread viruses initially discovered in bees may therefore be transmitted to the predator *Vespa velutina*, which may also be affected by some specific pathogens.

In this study, we investigated viral agents of *V. velutina* in order to identify specific viruses or strains that could be used as potential biocontrol agents of *V. velutina* populations. A number of viral microbial control agents have been successfully used against invasive species [51] or pathogens [52]. Therefore, we searched for any viral sequences using RNA sequencing in order to detect either undescribed viruses or viruses already described in other host species. To further investigate virus infection and whether some are tissue-specific, RNA-seq analysis was performed on different body parts. Lastly, we analysed the genetic diversity of the honey bee viruses from both hornet and bee samples to infer possible host-specific strains.
## 2. Materials and Methods

### 2.1. *Asian Hornet and Honey Bee Sampling*

Hornets were sampled on five occasions (one to four hornets per sampling) in front of honey bee colonies in two apiaries from the south-east of France during autumn 2014 (Table S1). Live hornets were kept at room temperature for 3 to 4 days and provided with water and sugar syrup *ad libitum* to foster elimination of previous feeding before dissections and then frozen at 20 °C. One hornet sampled by a beekeeper (in southeast France) exhibited deformed wings (159+, Figure 1) and was similarly frozen at -20 °C. Pools of bees from the apiaries where hornets were preying on were simultaneously collected (Table S1) and immediately frozen at -20°C in the laboratory. All samples were then stored at -80 °C.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Dorsal view of deformed wings in hornet (sample 159+ naturally curled up by freezing).</figcaption></figure>

Dissections of two asymptomatic hornets (individuals 140-1 and 140-2) and the hornet 159+ exhibiting deformed wings were performed on ice and resulted in brain, muscle, gut and mandible samples. The remaining external parts (cuticle and legs) were pooled together. For the deformed wings hornet, only muscle, intestine and external parts could be sampled (Table S1). For the dissection of other body parts, we used RNAlater®-ICE (Ambion, ref # AM7030). Head, thorax and abdomen of non-symptomatic individuals from the same apiary (140-3 and 144-1) were cut (avoiding thawing) and separately incubated in 1.7 mL RNAlater®-ICE at -80 °C for 18 h, then for 20 h at -20 °C. We then dissected the gut, the abdomen cuticle and fat bodies, and the six legs, which were pooled with the wings and thorax cuticle, the thorax muscles, the mandibles and the brain. Leftover samples of thorax and head were also kept for later analysis. The detailed list of the hornet and honey bee samples can be found in Table S1.

### 2.2. RNA Extraction

All the hornet samples (whole, or dissected tissues) were ground in 500 to 800 µL TRIzol® Reagent (Invitrogen, ref #15596-026, Carlsbad, California, USA), depending on sample size) with a 0.8 cm-diameter. bead in a 2 mL tube and a TissueLyser (Qiagen) (four times for 30 s each at 30 Hz and 30 s intervals). The tube was then centrifuged for 2 min at 12,000× g at 4 °C, and the supernatant was collected into a new tube to be processed for RNA extraction. RNA was extracted using an RNeasy mini kit (Qiagen, Courtaboeuf, France) and DNase according to the manufacturer's instructions. The final suspension volume was 100 µL for whole hornet and 50 µL for the dissected parts. RNA was quantified with a Nanodrop2000® (Thermo Scientific, Waltham, Massachusetts, USA) and RNA
quality was assessed with RNA Chips (Agilent technologies, Santa Clara, California, USA) using an
Agilent 2100 Bioanalyzer.

## 2.3. Next Generation Sequencing (NGS)

Next generation RNA-sequencing was performed by MGX-Montpellier GenomiX (Montpellier, France) on four samples showing the best RNA yield and quality: muscle from the symptomatic specimen (159M), and brain (140C), muscle (140M) and gut (144I) from non-symptomatic specimens. To enrich the samples for polyA and mRNA, ribosomal RNA was first depleted twice with a Ribo-zero kit (Magnetic Gold Epidemiology, Illumina Netherlands, Eindhoven, The Netherlands). The library was prepared with a TruSeq Stranded mRNA Sample Preparation kit (Illumina). Samples were paired-end sequenced (2 × 125 nt) with a Illumina HiSeq 2500 sequencer.

## 2.4. General Analysis of Viral Content

Short reads were quality checked (fastQC and contaminant search including E. coli, yeast and phiX) and stored in ng6.toulouse.inra.fr [53]. Reads were adapter-cleaned using cutadapt [54] and each of the five libraries was independently assembled with Oases [55] (version 0.2.06 using 25, 31, 37, 43, 49, 55, 61, 65, 69 kmers). The assembled contigs were merged and the overlapping ones removed using CD-HIT-EST [56] (Version v4.6.1-2012-08-27 with parameter -c 1). The Sequence Read Archive in Genbank is accessible with SRA accession: PRJNA556788.

For a first general analysis of viral content, all contigs were renamed using their sample name as the prefix and gathered in a unique file. This file was reduced by CD-HIT-EST in the same way. The resulting contigs were annotated with kraken [57] (version 0.10.5-beta with standard parameters). Kraken viral annotations were manually checked and only contigs with at least 20 kmers-long specific blocks for the set of viruses presented in Figure S1 were kept for further analyses. The number of contigs and their identifier are presented in Table 1. The reads were remapped to the global contig set using bwa mem [58] (version 0.7.12-r1039 with standard parameters), sorted and compressed using samtools [59] (view, sort and index version 1.1 with standard parameters) and the coverage of the selected contigs was visually inspected in IGV.

<table><caption>**Table 1.** Sequencing and alignment statistics. RNA was extracted from brain (140C), muscle (159M (symptomatic hornet), 140M), or gut (144I).</caption><thead><tr><th>Sample</th><th>Number of Read Pairs</th><th>Number of Contigs</th><th>Total Size of Contigs (in Nucleotides)</th><th>Read Alignment Rate (*)</th><th>Pair Alignment Rate (*)</th></tr></thead><tbody><tr><td>159M</td><td>41,650,349</td><td>85,968</td><td>44,138,196</td><td>99.35%</td><td>71.34%</td></tr><tr><td>144I</td><td>50,537,821</td><td>468,927</td><td>276,819,205</td><td>98.44%</td><td>41.19%</td></tr><tr><td>140M</td><td>49,181,345</td><td>420,974</td><td>246,505,091</td><td>98.81%</td><td>51.54%</td></tr><tr><td>140C</td><td>49,526,919</td><td>820,766</td><td>374,354,293</td><td>99.08%</td><td>62.55%</td></tr></tbody></table>

## 2.5. Analysis of Viral Genomes

A two-step mapping method was used as in previous virus discovery studies in mosquitoes [60] and wild bees [61] as this improves the recovery of full-length genomes and ensures that final contig assembly stems from reads originating from the hornet viral metagenome only. In short, in order to reconstruct virus genomes from the four assembled RNA sequencing datasets, full and partial ORFs were first detected in hornet contigs using Prodigal V2_60 software for metagenomic data [62]. Contigs with relevant virus homologies (e-value < 10⁻⁵) were then detected using BLASTx program [63] on the non-redundant (nr) protein database from the NCBI (National Center for Biotechnology Information).

<footnote>(*) from samtools flagstat.</footnote>
For each of the four libraries, assembled contigs matching expected virus genome size, or the presence of multiple viral hits belonging to the same viral family were an indication of putative viral genomes to be reconstructed.

Virus genome reconstruction was performed in three steps. First, for each potential virus found, the nucleotide genomic sequence of the virus species detected by BLAST as the first hit was retrieved from GenBank and used as a viral reference genome for further reconstruction. Paired-end reads from each library were mapped using Bowtie2 V2.3.3 [64] against all viral reference genomes, and only the best matches were reported. When necessary, the consensus sequences of the mapped reads were re-assembled with the initial contigs to fill the gaps. Secondly, the reads were mapped again on the consensus sequences of the targeted virus obtained at the first step. Thirdly, ORFs were predicted using the ORF finder program of the Geneious 9.1.8 software [65] and annotated using InterProScan version 5, [66], CD-search web using the CCD 3030 database and e-value < 0.05, and BLASTp (e-value < 0.001) against the nucleotide database. The sequences have been submitted to GenBank (accession numbers MN565031 to MN565059).

A more stringent mapping was also performed to calculate the mean coverage and expression level FPKM of the reconstructed genome and to estimate the quantity of polymorphism. For this mapping, paired-reads ends were trimmed (five bases at 5' ends and 9 bases at 3' ends) using the fastx-trimmer program from the standalone FASTX-Toolkit 0.0.13 when the median quality showed the lowest values (Q < 38) and when the quality variance was high. The number of bases trimmed was chosen according to Illumina Quality plots generated with fastx_quality_stats and Fastq质量_boxplot_graph.sh programs implemented in the FASTX-Toolkit (Figures S2 and S3). Trimmed paired-end reads were mapped against the reconstructed genomes using Bowtie2 as described above. Polymorphism Pi statistics was estimated using a Popoolation software with the disable correction parameter set to 1 [67]. To analyse the specific single nucleotide polymorphism (SNP) distribution of the reconstructed virus genomes across the samples, SNP detection was performed using Geneious Prime 2019 software, using a minimum coverage of 10X, the minimum variant frequency set at 0.3 and the other parameters set by default.

A Spearman correlation test between polymorphism and expression level was performed using the PAST v3 program [68], using 0.05 as the significance threshold.

## 2.6. Phylogenetic Analysis

Sequences were aligned with the Muscle program [69]. Site selection was performed using GBlocks [70] available in Seaview 4.4.4 suite [71] when necessary. Maximum likelihood phylogenies were built using the PhyML 3.0 program [72], except for DWV complete genomes, for which we used SplitsTree4 [73] to represent putative recombination events. The best substitution models for nucleotide and protein substitution were first assessed using the SMS [74] program and AIC selection criterion. Branch supports were calculated using a LRT SH-like likelihood-based method. Depending on the amount of published sequence data of closely-related sequences, trees were built either at the species level from nucleotide whole genome sequences alignments, or at a broader taxonomic level (typically genera or families) from RNA-dependent RNA polymerase (RdRp) ORF or domain amino acid alignments.

## 2.7. RT-PCR

Reverse transcription of either 1 µg RNA from the whole hornet samples or 144 ng to 1 µg RNA from the dissected parts was performed in 20 reaction volume using random priming and according to the manufacturer's protocol (High capacity RNA to cDNA, Life technologies).

Multiplex PCR was performed as previously described [75,76] to detect DWV, SBV, BQCV, CBPV, ABPV, KBV and IAPV. Briefly, an 8-µM master mix of the 11 primers was prepared and stored at -20 °C. Half of the cDNA (10 µL) was added to a mix of 0.5 µM primers, 2.5 mM MgCl2, 0.2 mM dNTPs and 2.5 U of Go Taq® G2 Flexi DNA polymerase (Promega) in 25 µl total volume. PCR reactions were
performed using an Eppendorf thermocycler. The following program was used: one cycle at 94 °C for 5 min, 35 cycles at 94 °C for 30 s, 56 °C for 30 s, and 72 °C for 45 s, and a final elongation cycle at 72 °C for 10 min. The reactions were then held at 10 °C. Five microliters of PCR products were run in 4% agarose gel, or one microliter was analysed with DNA Chips using an Agilent DNA 1000 reagent (Agilent technologies). The presence of the virus was determined according to the expected sizes of the fragments (IAPV: 158 bp; DWV: 269 bp; SBV: 342 bp; ABPV: 460 bp; BQCV: 536 bp; KBV: 641 bp; CBPV: 774 bp).

For SANGER sequencing, PCR was carried out with specific virus primers in the same conditions, except for DWV sequencing, which was performed from the (conserved) helicase coding region for partial sequencing as described in [77]. PCR products quality (purity and band size) was controlled in 1% agarose gel prior to sequencing. The PCR products were sent to GenoScreen (France) for purification and Sanger sequencing.

## 3. Results

### 3.1. Virus Discovery from RNAseq Experiments

RNAseq data were used for identifying viruses from hornet tissues sampled from one symptomatic and two asymptomatic hornets (brain, muscle or gut, Table S1). Initially, 1.7 M contigs were assembled, and after redundancy removal, 509,572 contigs could be annotated with Kraken: 113,336 sequences were classified (22.24%), and 396,236 sequences resulted unclassified (77.76%) (Table 1).

Overall, 19 viral species, including seven complete genomes and 12 partial sequences, could be identified from the four libraries (Table 2). Apart from honey bee-associated viruses, a high number of contigs was initially assigned as dsDNA viruses in brain, muscle, or intestine samples (Figure S1) in spite of DNase treatment during RNA extraction. However, further BLAST analyses suggested that these sequences corresponded to insect genes comprising conserved functional domains shared between eukaryotes and large dsDNA viruses. Indeed, no assembled contigs from this category allowed the reconstruction of full-length or even partial virus genomes.

Of the 509,572 virus contigs identified by Kraken, DWV-B (=VDV-1) was clearly predominant both in muscle sample 159M (40.3% of the classified contigs) and gut sample 144I (93.2% of the classified contigs) (Figure S1).

### 3.2. DWV-B is Highly Prevalent in Symptomatic and Asymptomatic Hornets

Firstly, complete DWV-B genomes could be reconstructed from the RNA libraries of individual organs (muscles 159M, gut 144I, brain 140C) (Table 2, Figure S4a-d). In the symptomatic sample 159M, the DWV-B virus was the only virus detected and no other viral contig could be found as compared to the other libraries.

Phylogenetic analysis from the DWV-B assembled contigs showed that the brain (140C), muscle (140M, 159M), and gut (144I) samples clustered together with the DWV-B reference sequence (NC006494-DWV-B) and not with DWV recombinant variants (HM067437, HM067437, KJ437447 and KX373900) (Figure 2). In the library, 144I various contigs were found to correspond to the DWV-C strain but in this case, the genome was segmented (Table 2).
<table><caption>**Table 2.** Taxonomy and genomic characteristics of viruses detected in hornet samples.</caption><thead><tr><th>Library</th><th>Virus Acronym</th><th>Virus Name</th><th>Family</th><th>Genus</th><th>Genome</th><th>Number of. Contig (s) Assembled</th><th>Contig(s) Size</th><th>Expected Full-Length Genome Size (kb)</th><th>% Missing Genome</th><th>Coverage of Assembled Contigs *</th><th>Coverage of a Full-Length Expected Genome</th></tr></thead><tbody><tr><td>159M</td><td>DWV-B</td><td>Deformed wing virus type B</td><td>Iflaviridae</td><td>Iflavirus</td><td>Nearly full-length</td><td>1</td><td>9434</td><td>10.1</td><td>6.7</td><td>11.1</td><td>10.4</td></tr><tr><td rowspan="3">140C</td><td>na</td><td>Vespa velutina associated acypi-like virus</td><td>unclassified</td><td>unclassified</td><td>Full-length</td><td>1</td><td>9872</td><td>10</td><td>0</td><td>535.1</td><td>535.1</td></tr><tr><td>DWV-B</td><td>Deformed wing virus type B</td><td>Iflaviridae</td><td>Iflavirus</td><td>Nearly full-length</td><td>4</td><td>9665</td><td>10.1</td><td>4.1</td><td>8.5</td><td>8.2</td></tr><tr><td>KBV</td><td>Kashmir bee virus</td><td>Dicistroviridae</td><td>Aparavirus</td><td>Partial</td><td>16</td><td>3536</td><td>9.5</td><td>62.9</td><td>1.9</td><td>0.7</td></tr><tr><td rowspan="3">140M</td><td>na</td><td>Vespa velutina associated acypi-like virus</td><td>unclassified</td><td>unclassified</td><td>Full-length</td><td>1</td><td>9962</td><td>10</td><td>0</td><td>1186</td><td>1186</td></tr><tr><td>DWV-B</td><td>Deformed wing virus type B</td><td>Iflaviridae</td><td>Iflavirus</td><td>Nearly full-length</td><td>2</td><td>10000</td><td>10.1</td><td>0.8</td><td>13.2</td><td>13.2</td></tr><tr><td>na</td><td>Vespa velutina associated triato-like virus</td><td>Dicistroviridae</td><td>Triatovirus</td><td>Partial</td><td>4</td><td>2380</td><td>9.3</td><td>74.3</td><td>4.5</td><td>1.2</td></tr><tr><td rowspan="14">144I</td><td>DWV-B</td><td>Deformed wing virus type B</td><td>Iflaviridae</td><td>Iflavirus</td><td>Full-length</td><td>1</td><td>10113</td><td>10.1</td><td>0</td><td>110,316.3</td><td>110,316.3</td></tr><tr><td>na</td><td>Vespa velutina associated acypi-like virus</td><td>unclassified</td><td>unclassified</td><td>Full-length</td><td>1</td><td>9748</td><td>10</td><td>0</td><td>122</td><td>122</td></tr><tr><td>ABPV</td><td>Acute bee paralysis virus</td><td>Dicistroviridae</td><td>Aparavirus</td><td>Full-length</td><td>1</td><td>9474</td><td>9.5</td><td>0</td><td>1794.6</td><td>1794.6</td></tr><tr><td>BQCV</td><td>Black queen cell virus</td><td>Dicistroviridae</td><td>Triatovirus</td><td>Full-length</td><td>1</td><td>8464</td><td>8.6</td><td>0</td><td>268.4</td><td>268.4</td></tr><tr><td>na</td><td>Vespa velutina associated Menton virus</td><td>unclassified</td><td>unclassified</td><td>Full-length segment 1</td><td>1</td><td>3031</td><td>2.8 (seg.1)</td><td>0 (Seg 2 missing: 1.5kb)</td><td>316.7</td><td>316.7</td></tr><tr><td>na</td><td>Vespa velutina associated ifla-like virus</td><td>unclassified</td><td>unclassified</td><td>Full-length</td><td>1</td><td>10023</td><td>9.9</td><td>0</td><td>19.5</td><td>19.5</td></tr><tr><td>na</td><td>Vespa velutina associated permutotetra-like virus 1</td><td>unclassified</td><td>unclassified</td><td>Full-length</td><td>1</td><td>4851</td><td>4.8</td><td>0</td><td>34.7</td><td>34.7</td></tr><tr><td>na</td><td>Vespa velutina associated permutotetra-like virus 2</td><td>unclassified</td><td>unclassified</td><td>Full-length</td><td>1</td><td>4742</td><td>4.8</td><td>0</td><td>11.7</td><td>11.7</td></tr><tr><td>na</td><td>Vespa velutina associated partiti-like virus 1</td><td>Partitiviridae</td><td>unclassified</td><td>Full-length segment 1</td><td>1</td><td>1779</td><td>1.7</td><td>0 (Seg 2 missing: 1.7 kb)</td><td>1082.5</td><td>1082.5</td></tr><tr><td>na</td><td>Vespa velutina associated partiti-like virus 2</td><td>Partitiviridae</td><td>unclassified</td><td>Full-length segment 1</td><td>1</td><td>1415</td><td>1.4</td><td>0 (Seg 2 missing: 1.7 kb)</td><td>44.6</td><td>44.6</td></tr><tr><td>na</td><td>Vespa velutina associated partiti-like virus 3</td><td>Partitiviridae</td><td>unclassified</td><td>Full-length segment 1</td><td>1</td><td>1417</td><td>1.4</td><td>0 (Seg 2 missing: 1.7 kb)</td><td>100.3</td><td>100.3</td></tr><tr><td>na</td><td>Vespa velutina associated partiti-like virus 4</td><td>Partitiviridae</td><td>unclassified</td><td>Partial segment 1</td><td>1</td><td>467</td><td>1.5</td><td>Seg1: 63.2 (Seg 2 missing: 1.7 kb)</td><td>6.2</td><td>1.9</td></tr><tr><td>na</td><td>Vespa velutina associated nora-like virus</td><td>unclassified</td><td>unclassified</td><td>Partial</td><td>1</td><td>1406</td><td>8</td><td>17.5</td><td>6.6</td><td>1.2</td></tr></tbody></table>
<table><caption>Table 2. *Cont.*</caption><thead><tr><th>Library</th><th>Virus Acronym</th><th>Virus Name</th><th>Family</th><th>Genus</th><th>Genome</th><th>Number of. Contig (s) Assembled</th><th>Contig(s) size</th><th>Expected Full-Length Genome Size (kb)</th><th>% Missing Genome</th><th>Coverage of Assembled Contigs *</th><th>Coverage of a Full-Length Expected Genome</th></tr></thead><tbody><tr><td></td><td>na</td><td>Vespa velutina associated triato-like virus</td><td>Dicistroviridae</td><td>Triatovirus</td><td>Partial</td><td>3</td><td>3149</td><td>9.3</td><td>66</td><td>4.9</td><td>1.6</td></tr><tr><td>ALPV</td><td>na</td><td>Aphid lethal paralysis virus</td><td>Dicistroviridae</td><td>Cripavirus</td><td>Partial</td><td>20</td><td>4186</td><td>9.8</td><td>57.4</td><td>2</td><td>0.9</td></tr><tr><td>BMLV</td><td>na</td><td>Moku virus</td><td>Iflaviridae</td><td>Iflavirus</td><td>Partial</td><td>23</td><td>6232</td><td>10.1</td><td>38</td><td>2</td><td>1.3</td></tr><tr><td>DWV-C</td><td>na</td><td>Bee Macula-like virus</td><td>Tymoviridae</td><td>na</td><td>Partial</td><td>13</td><td>3102</td><td>6.3</td><td>50.3</td><td>2.4</td><td>1.2</td></tr><tr><td></td><td>na</td><td>Deformed wing virus type C</td><td>Iflaviridae</td><td>Iflavirus</td><td>Partial</td><td>15</td><td>3292</td><td>10.1</td><td>64.9</td><td>9,634.2</td><td>3,377</td></tr></tbody></table>

<footnote>* Coverage calculated from the corresponding full-length expected.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Split tree of *Deformed wing virus* from 23 nucleotide sequences. All ambiguous positions were removed for each sequence pair (pairwise deletion option). There was a total of 7560 positions from the 3' end of the DWV genome in the final dataset. GenBank accession numbers and samples 140C (brain), 140M and 159M (muscle) and 144I (intestine) are associated with the DWV strain. Contigs from the intestine sample assigned as DWV-C were too short to be included.</figcaption></figure>

Secondly, multiplex PCR was used to study the presence of DWV in both hornets and honey bees sampled at the same location. DWV was detected from all the hornets and bee samples except for the dissected parts of sample 140-1, where only gut and muscle were positive (Table 2, Figure S5).

In order to determine whether the DWV strains in hornets and bees sampled together were the same, 11 helicase coding regions were sequenced. This confirmed that all the sequences from this study belonged to DWV species (Genbank Submission 2248150), but no strongly supported clusters could be identified in the phylogenetic tree because helicase is too conserved for fine phylogenetic resolution.

### 3.3. Six Other Bee-Associated Viruses Detected in Hornets by RNAseq

Aphid lethal paralysis virus (ALPV, *Dicistroviridae*, *Cripavirus*), Moku virus (*Iflaviridae*, *Iflavirus*) and a Bee macula-like virus (*BeeMLV*, unclassified *Tymoviridae*) were found in the hornet gut (144I). Complete genome could not be assembled for these viruses, but 4.2, 6.3 and 3.1 Kb contigs of the corresponding genomes could be assembled (Figure S6a-c). Phylogenies from the RdRp domain clearly showed that these viruses belonged to the known ALPV, Moku virus and BeeMLV species (Figure 3a,b and Figure S7).

Kashmir bee virus (KBV, *Dicistroviridae*, *Aparavirus*) was found in the brain (sample 140C), but not in the muscle of the same hornet individual. A partial 3.5-kb genome could be assembled (Figure S6d), but the RdRp domain was too short to produce a qualitative alignment and no phylogenies could be built. However, homology search strongly suggested that this sequence belonged to KBV species (best BLASTx hit with, for example, one High-Scoring Segment Pair displaying e = 100, 94% identity over 560 amino acids).

A 9.5-Kb complete genome of *Acute bee paralysis virus* (ABPV, *Dicistroviridae*, *Aparavirus*) and an 8.5-Kb complete genome of Black queen cell virus (BQCV, *Dicistroviridae*, *Triatovirus*) were found in the gut of hornet (144I) (Figure S6e-f). The taxonomic assignment of both ABPV and BQCV was again clearly confirmed by RdRp phylogeny (Figure S8).
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Maximum likelihood phylogeny of the replicase genes of (a) ALPV and *Vespa velutina* associated triato-like virus (372 aa, LG+I+G+F model), and (b) Moku virus and *Vespa velutina* associated ifla-like virus (338 aa, LG+I+G+F model) detected in hornets. The complete virus names are shown in Table S2.</figcaption></figure>
### 3.4. Eleven Other Insect Viruses Detected in Hornets from RNAseq

Three complete 9.7–9.9 Kb genomes of *Vespa velutina* associated acypi-like virus were detected in both the brain (140C) and muscle (140M) from the same hornet and in the gut of the 144I individual (Table 2, Figure 4a–c). The phylogenetic tree indicated that all three isolated *Vespa velutina* associated acypi-like virus clustered in a single clade (aLRT node statistics = 0.99) and showed similarities to the *Acyrthosiphon pisum* virus and the Rosy apple aphid virus, both proposed as Acypiviridae (node aLRT statictic = 0.95) (Figure 5a). The *Vespa velutina* associated acypi-like virus genome contained two ORFs. ORF1 product contained sequence motifs characteristic of Helicase (superfamily 3, PR000605) and RNA-dependent RNA polymerase (PF00680), a genomic organization similar to that of APV [78] (Figure 4a–c). Thus, *Vespa velutina* associated acypi-like virus is most certainly a new hornet virus species from the *Picornavirales* order and belongs to the newly proposed Acypiviridae family.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Schematic representation, coverage and annotation of the complete or partial genomes of insect RNA viruses found in hornets. (a–c) *Vespa velutina* associated acypi-like virus; (d–e) *Vespa velutina* associated triato-like virus and (f–g) *Vespa velutina* associated permutotetra-like virus 1 and 2.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Maximum likelihood amino acid replicase phylogeny of (a) *Vespa velutina* associated acypi-like virus (277 aa, RtREV +G+I+Fmodel), and (b) *Vespa velutina* associated permutotetra-like virus 1 and 2 (219 aa, LG+G model). GenBank accessions are shown in Table S2.</figcaption></figure>

*Vespa velutina* associated triato-like virus was identified from hornet muscle (140M) and gut (144I) (Table 2, Figure 4d–e). Coverage for *Himetobi P virus* (HPV)-like sequences did not allow the reconstruction of the whole virus genome, but only four annotated contigs from both muscle (140M) and gut (144I) samples, covering 2.4 kb. RNA-dependent RNA polymerase domains (PF00680) were identified in contig 1 and 3 from samples 140M and 144I, respectively, and phylogeny clearly showed that both *Vespa velutina* associated triato-like viruses belonged to genus *Triatovirus*, and are closely related to the *Bat dicistrovirus* species (node aLRT = 1) (Figure 3a). Along with *Bat dicistrovirus*, *Vespa velutina* associated triato-like virus is a sister group of the bee virus *Himetobi P virus* (HPV, *Dicistroviridae, Triatovirus*) (node aLRT = 0.97) and shares only 50% identity with HPV.

The remaining nine insect viruses were only detected in the gut transcriptome of 144I, and not in brain or muscle tissues. These viruses likely originate from infected preys ingested by hornets, and not from infected hornets themselves. *Vespa velutina* associated permutotetra-like virus 1 and 2 are two full-length genomes of 4.7–4.8 Kb (Table 2). The genomic organization of *Vespa velutina* associated permutotetra-like virus 1 and 2 was similar to that of the *Drosophila* A virus (unassigned) and the other
permutotetra-like viruses [79] (Figure 4f-g). RdRP phylogeny showed that *Vespa velutina* associated permutotetra-like virus 1 and 2 are two distinct species and are the sister group of the *Drosophila* A virus (Figure 5b).

*Vespa velutina* associated partiti-like viruses 1–3 are three full-length genomic segment 1 of 1.4–1.8 Kb, (Table 2, Figure 6a-d). *Vespa velutina* associated partiti-like virus 4 is fragmented, but its 0.5 Kb contig contained RdRp domain, which enabled us to build a phylogeny. Figure 7 shows that these four assembled genomes belonged to four new viruses of distinct genera and are distributed across the *Partitiviridae* tree.

<figure><img src="image_10.png" /><figcaption>b Vespa velutina associated partiti-like virus 2
seg1 144I</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>c Vespa velutina associated partiti-like virus 3
seg1 144I</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>d Vespa velutina associated partiti-like virus 4
seg1 144I</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>e Vespa velutina associated ifla-like virus 144I</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>f Vespa velutina associated nora-like virus 144I</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>g Vespa velutina associated menton virus 144I</figcaption></figure>

<figure><img src="image_16.png" /></figure>

**Figure 6.** Schematic representation, coverage and annotation of complete or partial genomes of insect RNA viruses found in hornets. (a–d): *Vespa velutina* associated partiti-like virus 1-4, (e): *Vespa velutina* associated ifla-like virus, (f): *Vespa velutina* associated nora-like virus, (g): *Vespa velutina* associated Menton virus.

*Vespa velutina* associated ifla-like virus is a 10-Kb full-length genome (Table 2, Figure 6e). Replicase phylogeny indicated that this new virus is a member of the *Iflaviridae* family (Order *Picornavirales*) (Figure 3b). *Vespa velutina* associated nora-like virus is a 1.4-Kb genomic fragment containing parvo coat and dsRNA-binding domains (Table 2, Figure 6f). Maximum likelihood phylogeny showed that this virus is closely related to the Nora virus discovered in the *Drosophila* species [80], a member of a newly described picorna-like virus family infecting arthropods [81] (Figure 8a).
<figure><img src="image_17.png" /><figcaption>**Figure 7.** Maximum likelihood amino acid replicase phylogeny of *Vespa velutina* associated partiti-like viruses 1–4 (384 aa, LG+I+G+F model). GenBank accessions are in Table S2.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>**Figure 8.** Maximum likelihood amino acid replicase phylogeny of (a) *Vespa velutina* associated nora-like virus (54 aa, LG+G model), (b) *Vespa velutina* associated Menton virus (162 aa, LG+G model). GenBank accessions are shown in Table S2.</figcaption></figure>

*Vespa velutina* associated Menton virus is a 3-Kb full-length genomic segment 1 containing serine protease and RdRp domains (Table 2, Figure 6g). Phylogeny indicated that this virus was a close relative of Motts Mill virus, Teise virus and Prestney Burn virus, detected in several *Drosophila* species [82–84] (Figure 8b). These viruses are members of an yet unclassified new family related to *Luteoviridae*, *Sobemoviridae*, *Barnaviridae* and *Alvernaviridae*.

## 3.5. *Viral Polymorphism Study Confirmed Genuine Hornet Infections*

Viral nucleotide diversity was estimated using Pi polymorphism in the three virus species found in two or more samples, i.e., *Vespa velutina* associated acypi-like virus, DWV-B, and *Vespa velutina* associated triato-like virus. For a given virus species, polymorphism levels were clearly distinct between samples, and a two-fold change was observed between the brain and intestine samples in *Vespa velutina* associated acypi-like virus (Table S3). Furthermore, Pi estimates were on average higher in intestine samples compared to the other organs. Finally, we confirmed that Pi estimates were not biased and influenced by a difference in viral genomes coverage and showed that Pi was not correlated with Fragment per Kilobase of viral sequences per million reads (FPKM) expression levels (Spearman's rho = 0.1, *p* = 0.78). These results show that viral sequences found in brain and muscles have different
genetic diversity signatures compared to the ones found in intestine, and thus that they did not result
from simple cross-contaminations.

To further investigate this point, specific patterns of single nucleotide polymorphisms (SNPs) were analyzed in *Vespa velutina* associated acypi-like virus. DWV-B and *Vespa velutina* associated triato-like virus were not included in this analysis since conditions on coverage and/or minimum variant frequency thresholds were not met. Table S4 shows the eight loci distributed across the viral genome were the SNP patterns highlighted in bold were clearly different between samples. In six loci, both brain and muscle samples shared the same SNP patterns, whereas intestine displayed distinct SNPs. Moreover, the SNP patterns were different in the two other loci (SNPs numbers 6.089 and 8.597) between the brain and muscle from the same hornet individual.

3.6. *Distribution of Bee-Associated Viruses in Hornets Detected by Multiplex RT-PCR*

Eight hornet samples were positive for SBV, as determined by RT-PCR (Table S1). A similar pattern of SBV presence was found in three honey bee pools sampled at the same locations than SBV-positive hornets (bee-pool-142, -145, -149). Only the gut, and no other dissected parts, was positive for SBV in 140-2 hornet (Table S1).

From the multiplex PCR method, eleven hornets out of twelve tested positive for ABPV, (Table S1). However, ABPV was detected by PCR in all the dissected parts of hornet 140-2 and both in the muscle and gut of hornet 140-1. Partial sequences from dissected parts were completely similar, even though they were sampled from two hornets and showed 99% identity with *Acute bee paralysis virus* isolate Hungary 1 (GenBank AF486072.2) isolated from bees (Genbank Submission 2248150). However, they were different from other ABPV sequences obtained from bees in the same experiment (29 informative positions out of 431).

Seven hornet samples were positive for BQCV, as determined by RT-PCR, and this virus was detected not only in the gut, as suggested from RNA-sequencing data in the sample 144I, but also in muscles (hornet 140-1) (Figure S6e, Table S1).

For both BQCV and ABPV, similar patterns of BQCV or ABPV presence or absence were found in
bee samples at the same locations than hornets (Table S1).

**4. Discussion**

Polymorphism and SNP analyses of viral signatures between samples indicated that the detection of the *Vespa velutina* associated acypi-like virus, DWV-B, and *Vespa velutina* associated triato-like virus sequences in the brain and muscle samples were not a result of simple cross-contamination that could have occurred during sample processing or library preparation. Instead, this study indicates that viruses typically associated with honey bees were actively infecting several hornet tissues. To date, oral inoculation has been the only known virus transmission route in hornets. The presence of viruses in the brain and muscle tissues combined with intra-species levels of polymorphism (SNP variation) indicate virus multiplication in hornets. The putative DWV circulation in the insect was further emphasized by the identification of DWV in other organs like mandibles or thorax cuticle and legs. Furthermore, DWV-B was shown to be predominant in all the RNA-sequencing samples. DWV-B read counts in the muscle sample were much higher in the hornet exhibiting symptoms of deformed wings (159M) than in the sample that did not show any symptoms (140M). The high DWV-B expression level may be related to the symptoms, when a low level of expression of the virus may correspond to a latent infection, as defined by de Miranda and Genersch [85], especially since DWV-B has been characterized, in vitro, as a more virulent honey bee virus strain than DWV-A [86]. However, one needs to be cautious regarding this assumption since only one hornet with a deformed wing symptom was analyzed. DWV-C was detected in a low number and only in hornet gut, which suggests food contamination or latent infection. Furthermore, since only partial sequences of DWV-C could be assembled, we could assume that they resulted from recombination with other DWV strains [77]. To our knowledge, this is the first description of DWV-C in France since its description in UK [46]. Finally, when comparing
DWV-B sequences from hornets to previously described DWV sequences from honeybees, they seem to be very close, which raises the question about the historical host range of DWV. The detection of DWV combined with genotypic diversity studies in the native area of *V. velutina* would clarify whether *V. velutina* is an usual host for DWV or if it only became infected from *A. mellifera* honeybees after its introduction in Europe.

ABPV and, to a lesser extent, BQCV were both identified in muscles and/or cuticle, legs, and mandibles, as well as in the guts. Similarly, food contamination alone seems unlikely for the presence of these viruses in hornets. BQCV is known to infect many *Apis* species (up to 90% prevalence in *A. mellifera*, other honeybees, bumblebees, and solitary bees), but also social wasps, with prevalence ranging from 2.2% to 67% prevalence [39]. Conversely, the ABPV host range is very limited. To date, ABPV has only been described in *A. mellifera*, *A. ceranae* [87], and bumblebees [88]. Our data show that *Vespidae* may be another host of such viruses. As a follow up, experimental infections will be required to assess the pathogenicity of these honey bee viruses (BQCV and ABPV) in hornets.

If we assume that circulation of viruses in insect body should be related to the infection cycle, then *Vespa velutina* associated acypi-like virus sequences, identified in both brain and muscle from 140C and 140M and whose whole genome could be reconstructed, may be an infective virus. In aphids, APV is known to reduce growth, to increase developmental time and possibly related to an increased mortality rate [89]. Further experiments are needed to characterize its virulence in *Vespidae* and to bees, such as *A. mellifera*.

Additionally, even if the whole *Vespa velutina* associated triato-like virus genome could not be reconstructed, the identification of sequences of this virus both in muscles and guts suggests that it may infect *V. velutina*. However, larger samplings are required to confirm its prevalence. Furthermore, the phylogenetic tree showed that these sequences are more similar to the *Bat dicistrovirus*, which has, to date, only been described twice in the feces of the bat species *Pipistrellus pipistrellus* [90,91]. Of note, this bat species has an insectivorous diet, which can explain the presence of this potentially insect-specific virus in their feces.

In contrast, the presence of viruses exclusively in the hornet gut may result from food contamination (infected preys or contaminated pollen) [38,53]. The absence of feces in cages where hornets were reared before freezing suggests that they could not eliminate the last food they ingested before sampling and sugar feeding. This could potentially lead to the detection of infected preys in the hornet gut content. Since DWV-C was only detected in the gut, we cannot state if they infect hornets or only their preys. The same conclusion can be drawn for all the other viruses described in this analysis, like Partiti-like virus 1, the *Motts Mill virus*, *Drosophila A virus* (DAV-like), the *Aphid lethal paralysis virus* (ALPV), the Bee Macula-like virus (BeeMLV), and Moku virus. The presence of Moku virus has already been described in an Asian yellow-legged hornet (*Vespa velutina nigrithorax*) collected in Belgium [33], but the viral metagenomics analysis was performed from a pool of 10 hornets, from which food contamination or real hornet infection could not be discriminated.

Our results from the French samples showed the presence of a large diversity of viruses in the invasive hornet *Vespa velutina*. DWV-B was predominant and may cause symptoms. Because DWV variants exhibited various pathogenicity levels, it is of great importance to assess their virulence for this new host species and to study their evolution over time, as recombination was already shown among variants [77]. Other viruses initially discovered in bees could infect the Asian hornet, since they were detected in different parts of the hornet body. The transmission routes from or to *V. velutina* have to be investigated to better assess the impact they may have on European honeybee and hornet populations. Moreover, new viruses were described, and the *Vespa velutina* associated acypi-like and *Vespa velutina* associated triato-like viruses were, for the first time identified in *Vespidae* species. In sum, there is an urgent need to study their host range and evaluate whether they represent a new threat to honeybees or a potential solution for the biocontrol of *V. velutina* in Europe, where they are dramatically spreading.
**Supplementary Materials:** The following are available online at http://www.mdpi.com/1999-4915/11/11/1041/s1.
Table S1: Hornet and honey bee sample list and corresponding viral PCR detections. Table S2: Accessions number and virus names of sequences used in virus phylogenies. Table S3: Polymorphism and expression levels of viral sequences between samples. Table S4: Read counts (and read %) supporting SNP from *Vespa velutina*-associated acypi-like virus across hornet samples. Figure S1: Kraken counts from reduced fasta contigs for brain (140C), muscle (140M, 159M) or gut (144I) samples. Only the most significant values and honey bee-associated viruses are shown. Figure S2: Quality plots of RNAseq reads of *Vespa velutina* samples 140C and 140M before FastQC trimming. Figure S3: Quality plots of RNAseq reads of *Vespa velutina* samples 144I and 159M before FastQC trimming. Figure S4: Schematic representation, coverage and annotation of complete *Deformed wing virus* genome found in hornets. (a–d): *Deformed wing virus* strain B found in 140C, 140M, 144I, and 159M, respectively. Figure S5: Multiplex PCR virus detection for six hornet and four bee samples originating from each apiary (Table S1). Ladder 100 bp is from Promega, the most intense band corresponds to a 500-bp sized fragment. Figure S6: Schematic representation, coverage and annotation of complete or partial genomes of RNA honey bee-associated viruses found in hornets. (a) *Aphid lethal paralysis virus*-like, (b) Bee macula-like virus, (c) Moku virus, (d) *Kashmir bee virus*, (e) *Acute bee paralysis virus* and (f) *Black queen cell virus*. Figure S7: Maximum likelihood amino acid phylogeny of Bee macula-like virus replicase from hornets (243 aa, LG+I+G+F model). GenBank accessions are shown in Table S2. Figure S8: Maximum likelihood amino acid replicase phylogeny of (a) *Acute bee paralysis virus* (139 amino acids, JTT+G model) and (b) *Black queen cell virus* (8,683 aa, GTR+G+I model). GenBank accessions are shown in Table S2.

**Author Contributions:** Conceptualization, A.D., Y.L.C., and C.A.; formal analysis, A.D., M.T., C.K., P.G., and D.B.; investigation, A.D., P.G., and E.A.H.; data curation, C.K., P.G., D.B., and A.D.; writing—original draft preparation, A.D., P.G., and D.D.; writing—review and editing, A.D., P.G., D.B., D.D., C.K., C.A., E.A.H., and Y.L.C.; funding acquisition, A.D.

**Funding:** This research was funded by the scientific department "plant health" of INRA.

**Acknowledgments:** The authors thank Audrey Jouvance for her technical support and Jean-Louis Lautard, beekeeper, for providing the symptomatic hornet.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. 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. [CrossRef]

2. Haxaire, J.; Bouguet, J.-P.; Jean-Philippe, T. *Vespa velutina* Lepeletier, 1836, une redoutable nouveauté pour la faune de France (Hym., Vespidae). *Bull. Soc. Entomol. Fr.* **2006**, *111*, 194.

3. Villemant, C.; Haxaire, J.; Streito, J.-C. La découverte du frelon asiatique *Vespa velutina*, en France. *Insectes 2006*, *143*, 3–7.

4. Arca, M.; Capdevielle-Dulac, C.; Villemant, C.; Mougel, F.; Arnold, G.; Silvain, J.F. Development of microsatellite markers for the yellow-legged Asian hornet, *Vespa velutina*, a major threat for European bees. *Conserv. Genet. Resour.* **2012**, *4*, 283–286. [CrossRef]

5. Arca, M.; Mougel, F.; Guillemaud, T.; Dupas, S.; Rome, Q.; Perrard, A.; Muller, F.; Fossoud, A.; Capdevielle-Dulac, C.; Torres-Leguizamon, M.; et al. Reconstructing the invasion and the demographic history of the yellow-legged hornet, *Vespa velutina*, in Europe. *Biol. Invasions* **2015**, *17*, 2357–2371. [CrossRef]

6. Rome, Q.; Muller, F.J.; Touret-Alby, A.; Darrouzet, E.; Perrard, A.; Villemant, C. Caste differentiation and seasonal changes in *Vespa velutina* (Hym.: Vespidae) colonies in its introduced range. *J. Appl. Entomol.* **2015**, *139*, 771–782. [CrossRef]

7. Darrouzet, E.; Gevar, J.; Guignard, Q.; Aron, S. Production of early diploid males by European colonies of the invasive hornet *Vespa velutina* nigrithorax. *PLoS ONE* **2015**, *10*, e0136680. [CrossRef]

8. Rome, Q.; Dambrine, L.; Onate, C.; Muller, F.; Villemant, C.; Garcia Pérez, A.L.; Maia, M.; Carvalho Esteves, P.; Bruneau, E. Spread of the invasive hornet *Vespa velutina* Lepeletier, 1836, in Europe. *Bull. Soc. Entomol. Fr.* **2013**, *118*, 21–22.

9. Monceau, K.; Bonnard, O.; Moreau, J.; Thiery, D. Spatial distribution of *Vespa velutina* individuals hunting at domestic honeybee hives: Heterogeneity at a local scale. *Insect Sci.* **2014**, *21*, 765–774. [CrossRef]

10. Bertolino, S.; Lioy, S.; Laurino, D.; Manino, A.; Porporato, M. Spread of the invasive yellow-legged hornet *Vespa velutina* (Hymenoptera: Vespidae) in Italy. *Appl. Entomol. Zool.* **2016**, *51*, 589–597. [CrossRef]

11. Barbet-Massin, M.; Rome, Q.; Muller, F.; Perrard, A.; Villemant, C.; Jiguet, F. Climate change increases the risk of invasion by the yellow-legged hornet. *Biol. Conserv.* **2013**, *157*, 4–10. [CrossRef]
12. Goldarazena, A.; de Heredia, I.P.; Romon, P.; Iturrondobeitia, J.C.; Gonzalez, M.; Lopez, S. Spread of the yellow-legged hornet *Vespa velutina* nigrithorax du Buysson (Hymenoptera: Vespidae) across Northern Spain. *EPPO Bull.* **2015**, *45*, 133–138. [CrossRef]

13. Bessa, A.S.; Carvalho, J.; Gomes, A.; Santarem, F. Climate and land-use drivers of invasion: Predicting the expansion of *Vespa velutina* nigrithorax into the Iberian Peninsula. *Insect Conserv. Diver.* **2016**, *9*, 27–37. [CrossRef]

14. Keeling, M.J.; Franklin, D.N.; Datta, S.; Brown, M.A.; Budge, G.E. Predicting the spread of the Asian hornet (*Vespa velutina*) following its incursion into Great Britain. *Sci. Rep.* **2017**, *7*, 6240. [CrossRef]

15. Lioy, S.; Manino, A.; Porporato, M.; Laurino, D.; Romano, A.; Capello, M.; Bertolino, S. Establishing surveillance areas for tackling the invasion of *Vespa velutina* in outbreaks and over the border of its expanding range. *Neobiota* **2019**, *46*, 51–69. [CrossRef]

16. Robinet, C.; Suppo, C.; Darrouzet, E. Rapid spread of the invasive yellow-legged hornet in France: The role of human-mediated dispersal and the effects of control measures. *J. Appl. Ecol.* **2017**, *54*, 205–215. [CrossRef]

17. Fournier, A.; Barbet-Massin, M.; Rome, Q.; Courchamp, F. Predicting species distribution combining multi-scale drivers. *Glob. Ecol. Conserv.* **2017**, *12*, 215–226. [CrossRef]

18. Rome, Q.; Perrard, A.; Muller, F.; Villemant, C. Aliens: The Invasive Species Bulletin: Newsletter of the IUCN/SSC Invasive Species Specialist Group. Available online: http://www.issg.org/pdf/aliens_newsletters/A31.pdf (accessed on 7 November 2019).

19. 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. [CrossRef]

20. Monceau, K.; Moreau, J.; Poidatz, J.; Bonnard, O.; Thiéry, D. Behavioral syndrome in a native and an invasive hymenoptera species. *Insect Sci.* **2015**, *22*, 541–548. [CrossRef]

21. Turchi, L.; Derijard, B. Options for the biological and physical control of *Vespa velutina* nigrithorax (Hym.: Vespidae) in Europe: A review. *J. Appl. Entomol.* **2018**, *142*, 553–562. [CrossRef]

22. Monceau, K.; Bonnard, O.; Thiery, D. Chasing the queens of the alien predator of honeybees: A water drop in the invasiveness ocean. *Open J. Ecol.* **2012**, *2*, 183–191. [CrossRef]

23. Monceau, K.; Thiéry, D. *Vespa velutina* nest distribution at a local scale: An 8-year survey of the invasive honeybee predator. *Insect Sci.* **2017**, *24*, 663–674. [CrossRef] [PubMed]

24. Haxaire, J.; Villemant, C. Impact sur l'entomofaune des "pièges à frelon asiatique". *Insectes* **2010**, *159*, 6.

25. Rojas-Nossa, S.V.; Novoa, N.; Serrano, A.; Calvino-Cancela, M. Performance of baited traps used as control tools for the invasive hornet *Vespa velutina* and their impact on non-target insects. *Apidologie* **2018**, *49*, 872–885. [CrossRef]

26. Decante, D.; ITSAP. Lutte Contre le Frelon Asiatique *Vespa velutina*. Available online: http://itsap.asso.fr/wp-content/uploads/2016/03/cr_evaluation_piegeage_vvelutina_2014.pdf (accessed on 18 April 2015).

27. Milanesio, D.; Maurice, S.; Maggiora, R.; Laurino, D.; Porporato, M. Recent upgrades of the harmonic radar for the tracking of the Asian yellow-legged hornet. *Ecol. Evol.* **2017**, *7*, 4599–4606. [CrossRef]

28. Milanesio, D.; Saccani, M.; Maggiora, R.; Laurino, D.; Porporato, M. Design of an harmonic radar for the tracking of the Asian yellow-legged hornet. *Ecol. Evol.* **2016**, *6*, 2170–2178. [CrossRef]

29. Kennedy, P.J.; Ford, S.M.; Poidatz, J.; Thiery, D.; Osborne, J.L. Searching for nests of the invasive Asian hornet (*Vespa velutina*) using radio-telemetry. *Commun. Biol.* **2018**, *1*, 88. [CrossRef]

30. Blot, J.; ITSAP. Localisation et Destruction des Nids de Frelons Asiatiques. Available online: http://www.itsap.asso.fr/downloads/publications/frelondestuction_light.pdf (accessed on 9 September 2015).

31. Beggs, J.R.; Brockerhoff, E.G.; Corley, J.C.; Kenis, M.; Masciocchi, M.; Muller, F.; Rome, Q.; Villemant, C. Ecological effects and management of invasive alien Vespidae. *Biocontrol* **2011**, *56*, 505–526. [CrossRef]

32. 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**, *71*, 211–216.

33. Mutien, G.; Bernard, T.; Noëmie El, A.; Daniel, C.; Gautier, G.; Marie, H.; Daniel, D.; Georges, D.; Annick, L.; Frédéric, F.; et al. Moku virus in invasive Asian hornets, Belgium, 2016. *Emeg. Infect. Dis.* *J.* **2017**, *23*, 2109.

34. Yanez, 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**, *43*, 587–589. [CrossRef]
35. Tehel, A.; Brown, M.J.F.; Paxton, R.J. Impact of managed honey bee viruses on wild bees. *Curr. Opin. Virol.* **2016**, *19*, 16–22. [CrossRef] [PubMed]

36. Evison, S.E.; Roberts, K.E.; Laurenson, L.; Pietravalle, S.; Hui, J.; Biesmeijer, J.C.; Smith, J.E.; Budge, G.; Hughes, W.O. Pervasiveness of parasites in pollinators. *PLoS ONE* **2012**, *7*, e30641. [CrossRef] [PubMed]

37. 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. [CrossRef]

38. 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. [CrossRef]

39. 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. [CrossRef]

40. Tentcheva, D.; Gauthier, L.; Zappulla, N.; Dainat, B.; Cousserans, F.; Colin, M.E.; Bergoin, M. Prevalence and seasonal variations of six bee viruses in *Apis mellifera* L. and *Varroa destructor* mite populations in France. *Appl. Environ. Microbiol.* **2004**, *70*, 7185–7191. [CrossRef]

41. Traynor, K.S.; Rennich, K.; Forsgren, E.; Rose, R.; Pettis, J.; Kunkel, G.; Madella, S.; Evans, J.; Lopez, D.; van Engelsdorp, D. Multiyear survey targeting disease incidence in US honey bees. *Apidologie* **2016**, *47*, 325–347. [CrossRef]

42. Dainat, B.; Evans, J.D.; Chen, Y.P.; Gauthier, L.; Neumann, P. Dead or alive: Deformed wing virus and *Varroa destructor* reduce the life span of winter honeybees. *Appl. Environ. Microbiol.* **2012**, *78*, 981–987. [CrossRef]

43. Francis, R.M.; Nielsen, S.L.; Kryger, P. Varroa-virus interaction in collapsing honey bee colonies. *PLoS ONE* **2013**, *8*, e57540. [CrossRef]

44. Genersch, E.; Aubert, M. Emerging and re-emerging viruses of the honey bee (*Apis mellifera* L.). *Vet. Res.* **2010**, *41*, 41:54. [CrossRef] [PubMed]

45. Kielmanowicz, M.G.; Inberg, A.; Lerner, I.M.; Golani, Y.; Brown, N.; Turner, C.L.; Hayes, G.J.R.; Ballam, J.M. Prospective large-scale field study generates predictive model identifying major contributors to colony losses. *PLoS Pathog.* **2015**, *11*, e1004816. [CrossRef] [PubMed]

46. Mordecai, G.J.; Wilfert, L.; Martin, S.J.; Jones, I.M.; Schroeder, D.C. Diversity in a honey bee pathogen: First report of a third master variant of the deformed wing virus quasispecies. *ISME J.* **2015**, *10*, 1264–1273. [CrossRef]

47. Mordecai, G.J.; Brettell, L.E.; Martin, S.J.; Dixon, D.; Jones, I.M.; Schroeder, D.C. Superinfection exclusion and the long-term survival of honey bees in Varroa-infested colonies. *ISME J.* **2015**. [CrossRef]

48. De Miranda, J.R.; Cordoni, G.; Budge, G. The Acute bee paralysis virus—Kashmir bee virus—Israeli acute paralysis virus complex. *J. Invertbr. Pathol.* **2010**, *103*, S30–S47. [CrossRef]

49. Ribière, M.; Olivier, V.; Blanchard, P. Chronic bee paralysis: A disease and a virus like no other? *J. Invertbr. Pathol.* **2010**, *103*, S120–S131. [CrossRef]

50. Celle, O.; Blanchard, P.; Olivier, V.; Schuff, F.; Cougoule, N.; Faucon, J.P.; Ribiere, M. Detection of chronic bee paralysis virus (CBPV) genome and its replicative RNA form in various hosts and possible ways of spread. *Virus Res.* **2008**, *133*, 280–284. [CrossRef]

51. Pedler, R.D.; Brandle, R.; Read, J.L.; Southgate, R.; Bird, P.; Moseby, K.E. Rabbit biocontrol and landscape-scale recovery of threatened desert mammals. *Conserv. Biol.* **2016**, *30*, 774–782. [CrossRef]

52. Bardin, M.; Ajouz, S.; Comby, M.; Lopez-Ferber, M.; Graillot, B.; Siegwart, M.; Nicot, P.C. Is the efficacy of biological control against plant diseases likely to be more durable than that of chemical pesticides? *Front. Plant Sci.* **2015**, *6*. [CrossRef]

53. Murray, E.A.; Burand, J.; Trikoz, N.; Schnabel, J.; Grab, H.; Danforth, B.N. Viral transmission in honey bees and native bees, supported by a global black queen cell virus phylogeny. *Environ. Microbiol.* **2019**, *21*, 972–983. [CrossRef]

54. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. *EMBnet J.* **2011**, *17*, 3. [CrossRef]

55. Schulz, M.H.; Zerbino, D.R.; Vingron, M.; Birney, E. Oases: Robust de novo RNA-seq assembly across the dynamic range of expression levels. *Bioinformatics* **2012**, *28*, 1086–1092. [CrossRef] [PubMed]

56. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. *Bioinformatics* **2012**, *28*, 3150–3152. [CrossRef] [PubMed]

57. Wood, D.E.; Salzberg, S.L. Kraken: Ultrafast metagenomic sequence classification using exact alignments. *Genome Biol.* **2014**, *15*, R46. [CrossRef] [PubMed]
58. Li, H.; Durbin, R. Fast and accurate short read alignment with burrows-wheeler transform. *Bioinformatics* **2009**, *25*, 1754–1760. [CrossRef] [PubMed]

59. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The sequence alignment/map format and SAMtools. *Bioinformatics* **2009**, *25*, 2078–2079. [CrossRef] [PubMed]

60. Bigot, D.; Atyame, C.M.; Weill, M.; Justy, F.; Herniou, E.A.; Gayral, P. Discovery of *Culex pipiens* associated Tunisia virus: A new ssRNA(+) virus representing a new insect associated virus family. *Virus Evol.* **2018**, *4*. [CrossRef]

61. Bigot, D.; Dalmon, A.; Roy, B.; Hou, C.; Germain, M.; Romary, M.; Deng, S.; Diao, Q.; Weinert, L.A.; Cook, J.M.; et al. The discovery of HalictiVirus Resolves the Sinaivirus phylogeny. *J. General Virol.* **2017**, *98*, 2864–2875. [CrossRef]

62. Hyatt, D.; Chen, G.L.; Locascio, P.F.; Land, M.L.; Larimer, F.W.; Hauser, L.J. Prodigal: Prokaryotic gene recognition and translation initiation site identification. *BMC Bioinf.* **2010**, *11*, 119. [CrossRef]

63. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. *BMC Bioinf.* **2009**, *10*, 421. [CrossRef]

64. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. *Nat. Meth.* **2012**, *9*, 357. [CrossRef] [PubMed]

65. 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. [CrossRef] [PubMed]

66. Jones, P.; Binns, D.; Chang, H.-Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G.; et al. InterProScan 5: Genome-scale protein function classification. *Bioinformatics* **2014**, *30*, 1236–1240. [CrossRef] [PubMed]

67. Kofler, R.; Orozco-terWengel, P.; De Maio, N.; Pandey, R.V.; Nolte, V.; Futschik, A.; Kosiol, C.; Schlotterer, C. PoPoolation: A toolbox for population genetic analysis of next generation sequencing data from pooled individuals. *PLoS ONE* **2011**, *6*, e15925. [CrossRef]

68. Hammer, Ø.; Harper, D.A.T.; Ryan, P.D. PAST: Paleontological statistics software package for education and data analysis. *Palaeontol. Electron.* **2001**, *4*, 9.

69. Edgar, R.C. Muscle: Multiple sequence alignment with high accuracy and high throughput. *Nucl. Acids Res.* **2004**, *32*, 1792–1797. [CrossRef]

70. Castresana, J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. *Mol. Biol. Evol.* **2000**, *17*, 540–552. [CrossRef]

71. Gouy, M.; Guindon, S.; Gascuel, O. SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. *Mol. Biol. Evol.* **2010**, *27*, 221–224. [CrossRef]

72. 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**, *59*, 307–321. [CrossRef]

73. Huson, D.H.; Bryant, D. Application of phylogenetic networks in evolutionary studies. *Mol. Biol. Evol.* **2006**, *23*, 254–267. [CrossRef]

74. Lefort, V.; Longueville, J.E.; Gascuel, O. SMS: Smart model selection in PhyML. *Mol. Biol. Evol.* **2017**, *34*, 2422–2424. [CrossRef] [PubMed]

75. Reynaldi, F.J.; Sguazza, G.H.; Albicoro, F.J.; Pecoraro, M.R.; Galosi, C.M. First molecular detection of co-infection of honey bee viruses in asymptomatic *Bombus atratus* in South America. *Braz. J. Biol.* **2014**, *73*, 797–800. [CrossRef] [PubMed]

76. Sguazza, G.H.; Reynaldi, F.J.; Galosi, C.M.; Pecoraro, M.R. Simultaneous detection of bee viruses by multiplex PCR. *J. Virol. Meth.* **2013**, *194*, 102–106. [CrossRef] [PubMed]

77. Dalmon, A. Evidence for positive selection and recombination hotspots in deformed wing virus (DWV). *Sci. Rep.* **2017**, *7*, 41045. [CrossRef]

78. Van der Wilk, F.; Dullemans, A.M.; Verbeek, M.; van den Heuvel, J.F. Nucleotide sequence and genomic organization of *Acyrthosiphon pisum* virus. *Virology* **1997**, *238*, 353–362. [CrossRef]

79. 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. [CrossRef]

80. Ekstrom, J.O.; Habayeb, M.S.; Srivastava, V.; Kieselbach, T.; Wingsle, G.; Hultmark, D. Drosophila Nora virus capsid proteins differ from those of other picorna-like viruses. *Virus Res.* **2011**, *160*, 51–58. [CrossRef]
81. Habayeb, M.S.; Ekengren, S.K.; Hultmark, D. Nora virus, a persistent virus in Drosophila, defines a new picorna-like virus family. *J. General Virol.* **2006**, *87*, 3045–3051. [CrossRef]

82. Medd, N.C.; Fellous, S.; Waldron, F.M.; Xuereb, A.; Nakai, M.; Cross, J.V.; Obbard, D.J. The virome of *Drosophila suzukii*, an invasive pest of soft fruit. *Virus Evol.* **2018**, *4*, vey009. [CrossRef]

83. Webster, C.L.; Longdon, B.; Lewis, S.H.; Obbard, D.J. Twenty-Five new viruses associated with the Drosophilidae (Diptera). *Evol. Bioinf. Online* **2016**, *12*, 13–25. [CrossRef]

84. Webster, C.L.; Waldron, F.M.; Robertson, S.; Crowson, D.; Ferrari, G.; Quintana, J.F.; Brouqui, J.M.; Bayne, E.H.; Longdon, B.; Buck, A.H.; et al. The discovery, distribution, and evolution of viruses associated with *Drosophila melanogaster*. *PLoS Biol.* **2015**, *13*, e1002210. [CrossRef] [PubMed]

85. De Miranda, J.R.; Genersch, E. Deformed wing virus. *J. Invertbr. Pathol.* **2010**, *103*, S48–S61. [CrossRef] [PubMed]

86. McMahon, D.P.; Natsopoulou, M.E.; Doublet, V.; Fürst, M.; Weging, S.; Brown, M.J.F.; Gogol-Döring, A.; Paxton, R.J. Elevated virulence of an emerging viral genotype as a driver of honeybee loss. *Proc. R. Soc. B Biol. Sci.* **2016**, *283*. [CrossRef] [PubMed]

87. Forsgren, E.; Shi, W.; Ding, G.L.; Liu, Z.G.; Tran, T.V.; Tang, P.T.; Truong, T.A.; Dinh, T.Q.; Fries, I. Preliminary observations on possible pathogen spill-over from *Apis mellifera* to *Apis cerana*. *Apidologie* **2015**, *46*, 265–275. [CrossRef]

88. 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. [CrossRef]

89. Van den Heuvel, J.F.; Hummelen, H.; Verbeek, M.; Dullemans, A.M.; van der Wilk, F. Characteristics of acyrthosiphon pisum virus, a newly identified virus infecting the pea aphid. *J. Invertebr. Pathol.* **1997**, *70*, 169–176. [CrossRef]

90. Dacheux, L.; Cervantes-Gonzalez, M.; Guigon, G.; Thiberge, J.M.; Vandenbogaert, M.; Maufrais, C.; Caro, V.; Bourhy, H. A preliminary study of viral metagenomics of french bat species in contact with humans: Identification of new mammalian viruses. *PLoS ONE* **2014**, *9*, e87194. [CrossRef]

91. Reuter, G.; Pankovics, P.; Gyongyi, Z.; Delwart, E.; Boros, A. Novel dicistrovirus from bat guano. *Arch. Virol.* **2014**, *159*, 3453–3456. [CrossRef]

<figure><img src="image_19.png" /></figure>

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).