Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Minireviews
    • JVI Classic Spotlights
    • Archive
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About JVI
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
Journal of Virology
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Minireviews
    • JVI Classic Spotlights
    • Archive
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About JVI
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
Genetic Diversity and Evolution

A Diverse Range of Novel RNA Viruses in Geographically Distinct Honey Bee Populations

Emily J. Remnant, Mang Shi, Gabriele Buchmann, Tjeerd Blacquière, Edward C. Holmes, Madeleine Beekman, Alyson Ashe
Susan R. Ross, Editor
Emily J. Remnant
aBehaviour and Genetics of Social Insects Laboratory, School of Life and Environmental Sciences, The University of Sydney, Sydney, Australia
bMarie Bashir Institute for Infectious Diseases and Biosecurity, Sydney Medical School, The University of Sydney, Sydney, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mang Shi
cCharles Perkins Centre, The University of Sydney, Sydney, Australia
dSchool of Life and Environmental Sciences, The University of Sydney, Sydney, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gabriele Buchmann
aBehaviour and Genetics of Social Insects Laboratory, School of Life and Environmental Sciences, The University of Sydney, Sydney, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tjeerd Blacquière
eWageningen University & Research, Wageningen, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Edward C. Holmes
bMarie Bashir Institute for Infectious Diseases and Biosecurity, Sydney Medical School, The University of Sydney, Sydney, Australia
cCharles Perkins Centre, The University of Sydney, Sydney, Australia
dSchool of Life and Environmental Sciences, The University of Sydney, Sydney, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Edward C. Holmes
Madeleine Beekman
aBehaviour and Genetics of Social Insects Laboratory, School of Life and Environmental Sciences, The University of Sydney, Sydney, Australia
bMarie Bashir Institute for Infectious Diseases and Biosecurity, Sydney Medical School, The University of Sydney, Sydney, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alyson Ashe
dSchool of Life and Environmental Sciences, The University of Sydney, Sydney, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Susan R. Ross
University of Illinois at Chicago
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/JVI.00158-17
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Understanding the diversity and consequences of viruses present in honey bees is critical for maintaining pollinator health and managing the spread of disease. The viral landscape of honey bees (Apis mellifera) has changed dramatically since the emergence of the parasitic mite Varroa destructor, which increased the spread of virulent variants of viruses such as deformed wing virus. Previous genomic studies have focused on colonies suffering from infections by Varroa and virulent viruses, which could mask other viral species present in honey bees, resulting in a distorted view of viral diversity. To capture the viral diversity within colonies that are exposed to mites but do not suffer the ultimate consequences of the infestation, we examined populations of honey bees that have evolved naturally or have been selected for resistance to Varroa. This analysis revealed seven novel viruses isolated from honey bees sampled globally, including the first identification of negative-sense RNA viruses in honey bees. Notably, two rhabdoviruses were present in three geographically diverse locations and were also present in Varroa mites parasitizing the bees. To characterize the antiviral response, we performed deep sequencing of small RNA populations in honey bees and mites. This provided evidence of a Dicer-mediated immune response in honey bees, while the viral small RNA profile in Varroa mites was novel and distinct from the response observed in bees. Overall, we show that viral diversity in honey bee colonies is greater than previously thought, which encourages additional studies of the bee virome on a global scale and which may ultimately improve disease management.

IMPORTANCE Honey bee populations have become increasingly susceptible to colony losses due to pathogenic viruses spread by parasitic Varroa mites. To date, 24 viruses have been described in honey bees, with most belonging to the order Picornavirales. Collapsing Varroa-infected colonies are often overwhelmed with high levels of picornaviruses. To examine the underlying viral diversity in honey bees, we employed viral metatranscriptomics analyses on three geographically diverse Varroa-resistant populations from Europe, Africa, and the Pacific. We describe seven novel viruses from a range of diverse viral families, including two viruses that are present in all three locations. In honey bees, small RNA sequences indicate that these viruses are processed by Dicer and the RNA interference pathway, whereas Varroa mites produce strikingly novel small RNA patterns. This work increases the number and diversity of known honey bee viruses and will ultimately contribute to improved disease management in our most important agricultural pollinator.

INTRODUCTION

Viruses are some of the most common pathogens contributing to declining honey bee health and colony losses worldwide (1–4), and at least 24 viruses have been described in the western honey bee, Apis mellifera (1, 5–12). Historically, viruses were identified based on pathological symptoms present in diseased bees by using serological methods (6, 7, 13, 14). Other than two DNA viruses (10, 11, 15), all currently characterized honey bee virus genomes comprise positive-sense, single-stranded RNA molecules (4, 16–19). Indeed, most viruses from honey bees fall into two families within the order Picornavirales (12, 16–21): the iflaviruses, including Sacbrood virus (SBV) and Deformed wing virus (DWV) (22, 23), and the dicistroviruses, including Black queen cell virus (BQCV) and Acute bee paralysis virus (ABPV) (24). More recently, genomic methods have identified additional viruses (5, 9, 10), which led to the discovery of a new genetic variant of DWV (23) and the Lake Sinai virus (LSV) group (5, 25). LSVs are common and widespread in honey bees, although whether they are associated with overt disease is unknown (26).

The prevalence, distribution, and virulence of honey bee viruses seem closely associated with the ectoparasitic mite Varroa destructor (2, 3, 22–24, 27). Varroa mites can act as an important virus vector, causing a dramatic change in both the viral landscape and virulence (3, 28, 29). As Varroa mites spread globally due to the human-mediated translocation of honey bees during the early to mid-20th century, so did virulent viruses, leading to the widespread loss of managed and wild honey bee colonies (2, 3, 27, 29, 30). Viruses from the Picornavirales appear to have a particularly close association with Varroa mites (8).

While Picornavirales are commonplace, it is striking that negative-sense RNA viruses apparently seem to be absent from honey bee colonies. This is puzzling, as negative-sense viruses, such as members of the Rhabdoviridae and Bunyaviridae, are widespread in other arthropods (31–33). Similarly, other categories of positive-sense RNA viruses, such as the Flaviviridae, are also common in insects but seemingly absent from studies of honey bees performed to date (34, 35). The Flaviviridae are notable, as they include mosquito and tick vector-borne viruses responsible for a number of important human diseases, including dengue, West Nile, and Zika, as well as other insect-specific viruses (35, 36).

Honey bees exhibit multiple antiviral defense mechanisms, including one of the key innate immune responses in insects, the RNA interference (RNAi) pathway (19, 37, 38). During virus replication, double-stranded RNA (dsRNA) intermediates are formed, which are recognized and cleaved by the endonuclease enzyme Dicer into 21- to 22-nucleotide (nt) fragments called small inhibitory RNAs (siRNAs) (39). These siRNAs are then bound by Argonaute proteins, which guide the RNA-induced silencing complex (RISC) and degrade complementary RNA molecules such as viral genomes (38). Honey bees produce siRNAs that match the predominant viruses in collapsing colonies (40) and can also produce a small RNA response when experimentally infected with double-stranded RNA (41, 42).

The use of RNA sequencing in a broad range of arthropods has greatly enriched our understanding of virus biodiversity (31, 33, 35, 43–46). It therefore seems opportune to revisit the viral diversity of honey bee colonies using similarly high-powered techniques. Accordingly, we screened honey bee populations for the presence of RNA viruses using total RNA transcriptome sequencing, so-called “metatranscriptomics” (33). We specifically focused on honey bee populations that have Varroa parasites but do not appear to suffer any negative consequences (47) and populations without Varroa mites. These populations were selected to avoid the possibility that novel viruses are outcompeted by highly virulent viral strains associated with Varroa mites. Notably, we also sampled bees from geographically diverse locations in Europe, Africa, and the Pacific to determine how viral diversity varies on a spatial scale.

RESULTS

Viral diversity in Varroa-resistant and Varroa-free honey bee populations.We examined the viral diversity of three Apis mellifera populations from Europe, Africa, and the Pacific by sequencing total ribosome-depleted RNA extracted from worker honey bees. We sampled honey bee colonies at the Amsterdam Water Dunes (The Netherlands), Stellenbosch and Robben Island (South Africa), and Vava'u and Tongatapu islands (Kingdom of Tonga) (Table 1). The colonies from The Netherlands are part of a selection program that started in 2008. These colonies are not treated for Varroa infestation but carry low numbers of mites (48). Colonies of Apis mellifera capensis in South Africa are similarly not treated for Varroa infestation and are naturally resistant to mites (49). The honey bee population on Robben Island became infested with Varroa mites 2 years prior to sampling, and numbers of mites per colony remain low (49). In the Pacific islands of Tonga, honey bees were most likely introduced during the 19th century, and large numbers of feral colonies are found in multiple island groups. Varroa mites were introduced to the island of Vava'u in 2006. Due to the lack of commercial beekeeping, colonies are not treated to remove mites, and, like the South African population, honey bees appear to be naturally tolerant to Varroa mites. The honey bees on Tongatapu island have never been exposed to Varroa mites but are derived from the same original source population as those on the island of Vava'u (50). We synthesized libraries from pooled RNA extracted from five individuals per colony, and 100-bp paired-end sequencing yielded a range of 4 to 9 Gb of data per library (Table 1). We assembled reads into contigs de novo using Trinity (51) and compared the resulting contigs to available viral protein sequences from GenBank with BLASTx.

View this table:
  • View inline
  • View popup
TABLE 1

Data generated in this study and summary of the viruses identified

We first examined the assembled contigs that matched previously characterized honey bee pathogens (Table 1). We found similarities in the presence and absence of known viruses in all three Varroa-resistant populations. Contigs for DWV, BQCV, and sacbrood virus (SBV) were present in all three locations. ABPV was present in one colony from Robben Island, and LSV was present in one colony each from The Netherlands and Tonga and in two colonies from South Africa (Table 1). The LSV genomes from The Netherlands, Tonga, and South Africa show 4 predicted open reading frames (ORFs), similarly to previously characterized LSV-1 and -2 genomes (Fig. 1). However, the nucleotide sequences from each location exhibit significant divergence, with 69 to 91% identity to previously characterized variants. We also examined our samples for contigs that matched common fungal, bacterial, and protozoan parasites of honey bees, such as the fungi Nosema apis and Nosema ceranae, the bacterial agents of European and American foulbrood, and the trypanosomes Crithidia mellificae and Lotmaria passim (4, 52). We observed contigs for L. passim in three colonies from Robben Island, South Africa, and one colony from The Netherlands (Table 1). The two colonies from mainland South Africa and one colony from The Netherlands contained contigs for Nosema apis, and one colony from Tonga contained contigs for Nosema ceranae, along with a single contig with similarity to Leishmania sp. (Table 1).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Lake Sinai virus variants identified in this study. (A) Genome structures of LSV strains identified in The Netherlands (GenBank accession number KY354242 ), Tonga (accession number KY354241 ), and South Africa (accession numbers KY354243 to KY354244 ), compared to previously characterized genomes of LSV-1 and LSV-2. Open reading frames are blue, and conserved functional domains are indicated (NCBI protein sequence accession numbers ARO50053 to ARO50067 ). (B) Maximum likelihood phylogenetic tree of nucleotide alignment of LSV strains from The Netherlands, Tonga, and South Africa with LSV-1 and -2 and other strains described previously (26).

Next, we examined contigs that showed similarity to previously characterized virus sequences from positive-sense, negative-sense, and double-stranded RNA genomes available in GenBank. We found genomic evidence of seven previously undescribed viruses in the honey bee transcriptomes, including four negative-sense and three positive-sense RNA viruses (Tables 1 and 2 and Fig. 2 and 3). Of the negative-sense RNA viruses, two belong to the family Rhabdoviridae of the order Mononegavirales (Table 2 and Fig. 2A), and two belong to the family Bunyaviridae (Table 2 and Fig. 2B). The three positive-sense RNA viruses include one virus belonging to the family Flaviviridae, one related to Nora viruses (picorna-like) found in Drosophila (53), and one belonging to the family Dicistroviridae of the order Picornavirales, with homology to Drosophila C virus (DCV) (Table 2 and Fig. 2C to E). For the new viruses described here, individual abbreviations are derived from the host species (e.g., Apis mellifera), followed by the name or category of virus and number if more than 1 (e.g., rhabdovirus 1).

View this table:
  • View inline
  • View popup
TABLE 2

Classifications, genome characteristics, and geographic locations of the novel viruses discovered in this study

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Genome structures of novel viruses. (A) Genome structures of rhabdo-like viruses, showing the genome size (nucleotides) and rhabdovirus open reading frames (N, P, M, G, and L/RdRp proteins) of ARV-1 (GenBank accession numbers KY354230 to KY354233 ) and ARV-2 (accession number KY354234 ) relative to the structure of the previously characterized Farmington virus. (B) Genome structure of bunya-like viruses, showing the identified L segment sizes (nucleotides) and the ORFs of ABV-1 (GenBank accession number KY354236 ) relative to LepmorLBV1 and of ABV-2 (accession number KY354237 ) relative to Wuhan mosquito virus 1. (C) Genome structure of a flavi-like virus, AFV (GenBank accession number KY354238 ), showing the genome size (nucleotides) and ORF relative to GKaV. (D) Genome structure of a dicistro-like virus, ADV (GenBank accession number KY354239 ), showing the genome size (nucleotides) and two ORFs of ADV relative to Drosophila C virus. (E) Genome structure of a Nora-like virus, ANV (GenBank accession number KY354240 ), showing the putative 5′-truncated genome size (nucleotides) and four ORFs relative to Drosophila Nora virus. (The NCBI protein sequence accession numbers are ARO50020 to ARO50052 .)

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Evolutionary relationships of novel viruses. Shown are maximum likelihood phylogenies of the novel rhabdoviruses ARV-1 and ARV-2 (A), the novel bunyaviruses ABV-1 and ABV-2 (B), the novel flavivirus AFV (C), the novel dicistrovirus ADV (D), and the novel Nora virus ANV (E). (See Fig. S1 to S3 in the supplemental material for detailed trees for panels A to C.)

Novel negative-sense RNA viruses. (i) Rhabdo-like viruses.Apis mellifera rhabdovirus 1 (ARV-1) has a 14,613-nt genome with a prototypic rhabdovirus structure (54) corresponding to the conserved gene order with five ORFs (Fig. 2A). The most conserved ORF encodes a 2,143-amino-acid (aa) protein containing the RNA-dependent RNA polymerase (RdRp) domain. ARV-1 RdRp was most similar to that of Farmington virus (FARV), a virus originally isolated from birds (55), with 30% amino acid identity (Table 2). The predicted G protein of ARV-1 was also most similar to that of FARV (18%); however, BLAST searches of the other ORFs showed no significant similarity to any known sequences. ARV-2 consists of a 14,029-nt genome with five predicted ORFs (Fig. 2A). The RdRp of ARV-2 was also related to that of FARV but with much lower sequence similarity (23% amino acid identity) (Table 2) than for ARV-1. Unlike the ARV-1 G protein, the ARV-2 G protein showed structural and sequence similarities to the hemagglutinin protein of Quaranfil virus of the Orthomyxoviridae family (22% identity), which suggests that this glycoprotein gene may be acquired through an inter-virus-family horizontal transfer event (33). Interestingly, both ARV-1 and ARV-2 show evidence of being widespread. ARV-1 was found in six colonies and in all three geographically diverse locations (The Netherlands, South Africa, and Tonga) (Tables 2 and 3), with a high abundance ranging from 50 to 500 transcripts per million (TPM) (Table 3). ARV-2 was moderately abundant in South Africa and Tonga (5 to 17 TPM) and was also detected in The Netherlands (1.8 TPM) (Tables 2 and 3). The ARV-1 and ARV-2 genomes from each location exhibit 98 to 99% nucleotide identity to each other (Fig. 4). The RdRp protein sequences of ARV-1 and ARV-2 formed a monophyletic group with that of FARV, which were distantly related to other members of the order Mononegavirales (Fig. 3A; see also Fig. S1 in the supplemental material).

View this table:
  • View inline
  • View popup
TABLE 3

Abundances and prevalences of novel viruses in sampled colonies

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Apis mellifera rhabdovirus 1 variants. Shown is a maximum likelihood phylogenetic tree of the nucleotide alignment of ARV-1 variant genomes isolated from The Netherlands (GenBank accession number KY354230 ), South Africa (accession number KY354231 ), and Tonga (accession number KY354232 ).

(ii) Bunya-like viruses.We identified two bunya-like virus sequences, both of which were obtained from colonies from Robben Island (Fig. 2B). Members of the Bunyaviridae are negative-sense RNA viruses with three genome segments (L, M, and S), each containing a separate ORF (31). Due to a lack of sequence similarity to known sequences, for both viruses, we identified only the larger L segments encoding the RdRp domains. Apis mellifera bunyavirus 1 (ABV-1) was present in three colonies (Table 3). ABV-1 was most similar to a recently characterized class of bunyaviruses, the leishbunyaviruses, identified in the insect trypanosomatid parasite Leptomonas moramango (56). The RdRp protein of ABV-1 exhibits 56% amino acid identity to Leptomonas moramango leishbunyavirus 1 (LepmorLBV1) (Table 2). ABV-2 was present in one colony from Robben Island. The RdRp protein has 42% identity to Wuhan mosquito virus 1 (Table 2) (31). Phylogenetic analysis involving representative members from the Bunyaviridae placed ABV-1 in a clade with the leishbunyaviruses as a clade basal to other bunyaviruses found within invertebrates (31, 56) and placed ABV-2 in the cluster of Phasmavirus-like bunyaviruses (Fig. 3B and Fig. S2).

Novel positive-sense RNA viruses. (i) Flavi-like virus.Apis mellifera flavivirus (AFV) was identified in one colony from Robben Island. The 20,414-nt positive-sense RNA genome contains a single ORF of 6,615 aa (Fig. 2C). BLAST searches indicated that AFV has 20% amino acid identity to gentian Kobu-sho-associated virus (GKaV), a recently identified flavi-like virus originally thought to be a dsRNA virus (Table 2) (35, 57, 58). Similarly to GKaV and other newly identified flaviviruses, the 20.4-kb AFV-1 genome is longer than the typical length of previously characterized members of the Flaviviridae (35). The phylogeny based on the RdRp/NS5 domain of AFV and other members of the Flaviviridae placed AFV in a clade of other recently discovered flavi-like viruses with large genomes (35) (Fig. 3C; see also Fig. S3 in the supplemental material) and flavi-like virus segments identified in Drosophila species (Takaungu virus and Hermitage virus) (34).

(ii) Dicistro-like virus.We identified a novel dicistrovirus from one colony at the Amsterdam Water Dunes (Apis mellifera dicistrovirus [ADV]). The 9,126-nt genome contains two ORFs encoding the replication enzyme polypeptide and the capsid proteins, respectively, which is typical of dicistroviruses (Fig. 2D). The polypeptide containing the RdRp exhibited the highest genetic identity to Drosophila C virus (57% amino acid identity) (Table 2). Phylogenetic analysis placed ADV in the same clade as Cricket paralysis virus (CrPV), Drosophila C virus, and Anopheles C virus (Fig. 3D).

(iii) Nora-like virus.Twelve separate contigs were assembled from one colony from Robben Island, each of which showed similarity to Drosophila Nora virus after BLASTx analysis. To assemble a full-length genome, these contigs were ordered according to their most closely related virus, Drosophila pseudoobscura Nora virus, and gaps were filled by using reverse transcription-PCR (RT-PCR) and Sanger sequencing using primers spanning the neighboring contigs. The resulting Apis mellifera Nora virus (ANV) partial genome sequence is 10,091 nt long and covers the entire replicase, although it is missing the first ORF at the 5′ end of a typical Nora virus genome (Fig. 2E) (53). In the phylogenetic tree, the ANVs were closely related (52% to ∼54% amino acid identity) to Nora viruses isolated from different Drosophila species (Fig. 3E).

Small RNA profiles of ARV-1 and ARV-2 in honey bees.One way of confirming that a putative virus genuinely infects the host from which it is sampled is the presence of an antiviral immune response. In insects, likely candidates are the small RNA pathways that are utilized in viral defense (37, 40, 59). We therefore determined the presence of antiviral small RNAs in bees infected with our novel viruses. We focused on ARV-1 and ARV-2 because they are the first negative-sense RNA viruses described in honey bees and were present in colonies sampled from all three of our geographically diverse locations (Table 1).

We generated small RNA libraries from the abdomens of four A. mellifera samples: two from the Amsterdam Water Dunes and two from Robben Island. For each of these geographic locations, we used PCR to screen for individuals that were positive for ARV-1 and used one individual that tested positive for ARV-1 (Amsterdam Water Dunes positive [AWD+] and Robben Island positive [RI+]) and one individual that tested negative (Amsterdam Water Dunes negative [AWD−] and Robben Island negative [RI−]). The RI+ sample also tested positive for ARV-2. The four libraries were subjected to 50-bp single-end sequencing, resulting in between 10 million and 18 million reads per sample.

We first mapped the resulting small RNA reads to the Apis mellifera genome and then aligned the unmapped reads to the ARV-1 and ARV-2 genomes (Table 4). From this, we found highly abundant small RNAs mapping to ARV-1 and ARV-2 (Table 5). Such small RNAs could either be random degradation products of viral RNA or result from the antiviral immune response of honey bees. Random degradation products of negative-sense RNA viruses would show a mixed size distribution of predominantly negative-sense fragments spanning the entire viral genome. In contrast, our small RNA reads have a size distribution of 21 to 22 nt, occur equally in sense and antisense orientations, and map predominantly to the 5′ and 3′ ends of the ARV-1 and -2 genomes (Fig. 5A to D). These features are typical signatures of Dicer-produced antiviral RNAs, which occur when Dicer binds to a double-stranded RNA intermediate and cleaves the double-stranded RNA into viral siRNA (39).

View this table:
  • View inline
  • View popup
TABLE 4

Small RNA samples and data generated

View this table:
  • View inline
  • View popup
TABLE 5

Number of small RNA reads mapped to ARV-1 and ARV-2

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Small RNA analysis of ARV-1 and ARV-2 in honey bees. (A and B) Size distributions (15 to 37 nt) and 5′-nucleotide compositions of small RNAs arising from ARV-1 (A) and ARV-2 (B). The sample from which the small RNA library was produced is shown above each graph. Bars plotted above the x axis represent reads that map to the positive strand, and those plotted below represent those that map to the negative strand. Bars are colored according to the proportions of reads starting with A, C, G, and T. (C and D) Mapping of 20- to 23-nt-long viral RNAs (viRNAs) to the genomes of ARV-1 (C) and ARV-2 (D). The cartoon shows the domains of the viral genomes as shown in Fig. 2. (E and F) Phasing score over 8 phasing cycles for each position within a 22-nt phasing window. The top and bottom graphs show the phasing scores for the sense and antisense reads, respectively. A high score indicates that many small RNAs fall into that phase position (indicated with arrowheads). This analysis was performed by using the 21- to 23-nt-long reads from panels A and B. (G and H) Observed 5′ nucleotide (Obs) compared with that expected (Exp) from the base compositions of the viral genomes for ARV-1 (G) and ARV-2 (H). Sense (S) and antisense (AS) reads were compared by using a chi-squared test. **, P value of <0.01.

The 5′ and 3′ genome biases suggest that replication intermediates at the ends of the ARV-1 and ARV-2 genomes provide dsRNA termini for Dicer to bind. Dicer-produced antiviral RNAs occur at regularly spaced 21- to 23-nt intervals (phases) starting at the dsRNA termini and fading with increasing distance from the termini, with a characteristic 2- to 3-nt overhang (37). We looked for evidence of small RNAs occurring at regularly spaced intervals from the 5′ end of the ARV-1 and ARV-2 genomes using phasing analysis (60). We detected a strong phasing signature for ARV-1, 7 nt from the 5′ end, with a 2-nt overhang between the sense and antisense strands (Fig. 5E). We also detected a phasing signature for ARV-2, 2 nt from the 5′ end, with an offset of 3 nt between the sense and antisense strands (Fig. 5F). These data strongly indicate that Dicer is responsible for producing the 22-nt small RNAs.

Finally, the antiviral immune response is also mediated by RNA-binding proteins such as Argonaute proteins, which bind to small RNAs and induce the degradation of RNA sequences complementary to the small RNA (61). Argonaute proteins often show a 5′-nucleotide preference in small RNA molecules (61), so we looked for nucleotide bias at the 5′ end of our small RNA compared to the base composition of the viral genome. Both the sense and antisense small RNAs against ARV-1 and ARV-2 display a highly significant reduction in 5′-G and a highly significant increase in 5′-U as the 5′ nucleotides (P < 0.01 for both by a chi-squared test) (Fig. 5G and H).

Taken together, our data suggest that the ARV-1 and -2 small RNAs have been generated by Dicer acting on a double-stranded RNA replication intermediate and that the small RNAs are subsequently bound by Argonaute proteins, indicating that the bees have an active antiviral immune response against ARV-1 and ARV-2.

Small RNA profiles of ARV-1 and ARV-2 in mites.We next wanted to determine if ARV-1 and -2 are also found in mites feeding on infected bees. To this end, we generated small RNA libraries from two V. destructor mites collected from A. mellifera individuals from Robben Island and performed Illumina 50-bp single-end sequencing, resulting in 18 million to 20 million reads per sample (Table 4). We found small RNA reads mapping to ARV-1 and ARV-2 in both mite samples (Table 5 and Fig. 6).

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Small RNA analysis of ARV-1 and ARV-2 in Varroa. Left panels show ARV-1, and right panels show ARV-2. (A and B) Size distributions (15 to 37 nt) and 5′-nucleotide compositions of small RNAs in mite 1 (M1) arising from ARV-1 (A) and ARV-2 (B). Bars plotted above the x axis represent reads that map to the positive strand, and those plotted below represent those that map to the negative strand. Bars are colored according to the proportions of reads starting with A, C, G, and T. (C and D) Mapping of 23- to 25-nt-long viRNAs to the genomes of ARV-1 (C) and ARV-2 (D). The cartoon shows the domains of the viral genomes as shown in Fig. 2. (E and F) Size distributions and 5′-nucleotide compositions of the sense small RNAs from panels A and B, respectively, normalized to the number of sense reads present. Mapping of the 19- to 24-nt sense reads to the viral genomes is also shown. (G and H) Phasing scores over 8 phasing cycles for each position within a 24-nt phasing window for ARV-1 (G) and ARV-2 (H). The top and bottom graphs show the phasing scores for the sense and antisense reads, respectively. This analysis was performed by using the 24-nt-long reads only. (I and J) Observed 5′ nucleotides (Obs) compared with those expected (Exp) from the base compositions of the viral genomes of ARV-1 (I) and ARV-2 (J). Sense (S) and antisense (AS) reads were compared by using a chi-squared test. *, P value of <0.05; ns, not significant. (K and L) Distances between the 5′ ends of overlapping reads on opposite strands (left) and the base compositions of each nucleotide position (right) for the 23- to 25-nt-long reads of ARV-1 (K) and ARV-2 (L).

The small RNA reads present in the mites show markedly different characteristics compared to the small RNA patterns in the honey bee samples. In mites, the antisense reads vastly outnumber the sense reads, span the length of the genome, and have a distinct length distribution centered at 24 nt. Furthermore, the antisense reads correspond to the first four ORFs encoding the ARV-1 and ARV-2 proteins (Fig. 6A to D). The sense reads, while much less abundant, show a broader size range than the antisense reads, with peaks at 20 and 23 nt, and do not localize specifically to the ORFs (Fig. 6E and F and data not shown). Phasing analysis from the 5′ end of the viral RNA did not show any evidence of phased RNAs of any size for either sense or antisense reads (Fig. 6G and H). Furthermore, there was only a weak 5′-nucleotide bias against 5′-G and toward 5′-U compared to the base composition of the viral genome in sense reads for ARV-1 (P < 0.05 by a chi-squared test), no 5′-nucleotide bias for antisense reads for ARV-1 (not significant [NS] by a chi-squared test), and no 5′-nucleotide bias for sense or antisense reads mapping to ARV-2 (NS by a chi-squared test) (Fig. 6I and J). The lack of phasing, weak 5′-nucleotide bias, and small quantity of sense reads suggest that the Dicer and Argonaute antiviral pathways do not act on replicating ARV-1 and ARV-2 dsRNAs in mites.

We also tested the 23- to 25-nt-long reads for signatures of the piwi-interacting RNA (piRNA) ping-pong amplification pathway, as this pathway has been implicated in viral defense in Aedes aegypti (62, 63). These signatures include a bias for uridine at the 5′ end and an adenine at position 10 on the complementary piRNA as well as a 10-nucleotide distance between the 5′ ends of overlapping sense and antisense reads (64). We generated a heat map showing the relative enrichment of each nucleotide at each position in the 23- to 25-nt-long small RNAs. Although we saw some evidence of U enrichment at position 1, there was no evidence to suggest an enrichment of A at position 10. We also plotted the distance between the 5′ ends of overlapping 23- to 25-nt-long RNAs that map to ARV-1 and ARV-2 but found no evidence for a peak at 10 nt. Thus, we could not detect any evidence of a “ping-pong” signature in the reads originating from either ARV-1 or ARV-2 in either mite sample (Fig. 6K and L and data not shown). Taken together, these data could mean that the small RNAs from ARV-1 and ARV-2 detected in mites are simply the products of random degradation, which is most likely the case for the sense small RNAs. However, in the antisense reads, the strong bias for 24-nt RNAs argues against random degradation.

The small RNA profiles in mites do not mimic those in honey bees, indicating that the mite reads do not come from ingested honey bee material. Thus, the high abundance of reads in the mites (Table 4) and the presence of sense reads (showing that a sense genome has been produced) suggest that both mites and honey bees are bona fide hosts of ARV-1 and ARV-2.

DISCUSSION

We describe a diverse set of new viruses in Varroa-resistant or Varroa-free honey bee populations from three locations in Europe, Africa, and the Pacific. We present genomic evidence of seven new RNA viruses, including three novel positive-sense and four novel negative-sense viruses. Our study therefore increases the number of known honey bee viruses from 24 to 31. ARV-1 and -2 were found in three geographically distinct populations. Using small RNA sequencing, we show that honey bees exhibit classic Dicer-mediated siRNA profiles, suggesting an active bee immune response. We also report the first analysis of small RNAs in Varroa destructor mites and show that ARV-1 and -2 are present in these mites, although the small RNA profile is distinctly different from that for honey bees.

To our knowledge, this is the first identification of negative-sense viruses in honey bees. Three of the four novel negative-sense viruses (ARV-1, ARV-2, and ABV-2) are related to viruses known to be present in insects (65). Our findings are thus consistent with data from recent studies describing the wide distribution of negative-sense viruses in arthropod hosts (31, 32, 66). Indeed, two negative-sense viruses were also recently found in the wild solitary bee Osmia cornuta, including one virus from the order Mononegavirales and one from the family Bunyaviridae, indicating that viral diversity in other hymenopteran species extends to negative-sense viruses (67). One of the novel negative-sense viruses from the Bunyaviridae identified in our South African population, ABV-1, was closely related to a recently isolated protist-infecting virus, the leishbunyavirus LepmorLBV1 (56). LepmorLBV1 was isolated from the insect trypanosomatid parasite Leptomonas moramango, a parasite of Brachycera flies (68). Interestingly, the three colonies from Robben Island that contained ABV-1 contigs were also positive for the honey bee trypanosome Lotmaria passim (52). We therefore cannot exclude that ABV-1 is a virus of protists that infect bees.

The small RNA patterns of the rhabdoviruses in honey bees show classical Dicer-mediated degradation profiles, providing strong evidence that ARV-1 and -2 enter the cells of bees and begin to replicate. During viral replication, a double-stranded RNA replication intermediate is formed, which can be recognized by the RNAi machinery and chopped by Dicer. The resulting small RNAs become part of an antiviral immune response (59). This strongly suggests that these novel negative-sense viruses are bona fide viruses capable of replication in honey bees.

In contrast, in Varroa, the strong bias for 24-nt antisense small RNAs and the lack of phasing suggest that replicating dsRNA is not a template for Dicer. Internal secondary structures within a negative-sense virus genome can provide dsRNA templates for Dicer (69, 70). An alternative Dicer-mediated antiviral response seems to be the most parsimonious explanation for the small RNA profile found for the mites. Interestingly, however, the predominant size of 24-nt antisense reads is larger than expected for canonical Dicer products. It is unclear whether Dicer in Varroa destructor produces longer-than-usual RNA fragments or if these 24-nt RNAs are generated by a different viral degradation pathway. Other mites, such as Tetranychus and Metaseiulus mites, contain the components of the RNAi machinery, with considerable variation in gene copy numbers for Dicer and Argonaute proteins (71). An important step toward understanding the Varroa RNA interference pathway will be determining if the Varroa genome contains similar variations in key RNAi-mediating genes. It is intriguing that the 24-nt antisense reads correspond to the ORFs of ARV-1 and -2. Similar “hot spots” have been observed previously, but their functional relevance is unclear (72).

Our study substantially expands the taxonomic diversity of honey bee viruses. Until now, most characterized honey bee viruses were restricted to the order Picornavirales (18, 19). The invertebrate-specific Dicistroviridae and Iflaviridae classes are evidently well adapted to parasitizing insects, and many viruses in these groups show extremely broad host ranges (73–76), which may facilitate spread, allow viruses to replicate more readily in multiple hosts, and thus allow viruses to become more prevalent and easily detected. Positive-sense RNA viruses are also more abundant in eukaryotes generally (77), which likely contributes to the frequency with which they are detected in honey bees. Indeed, we found three novel positive-sense RNA viruses. Two of these viruses, ANV and ADV, show relatively close evolutionary relationships with similar viruses from Drosophila. Interestingly, early serological characterizations of DCV and the related CrPV included honey bee samples for cross-reactivity to CrPV and DCV sera and identified a honey bee variant (73, 78, 79). The third positive-sense virus identified here, AFV, follows from the recent identification of larger flavivirus genomes (35), suggesting that arthropods could harbor a variety of viruses that will further illuminate the evolutionary origins of common viral categories (31, 35).

Until now, viral surveys in honey bees predominantly focused on PCR-based approaches, which were heavily dependent on existing virus diversity (2, 3, 16, 25, 27) or on infectivity tests in honey bees, where viruses were identified based on their ability to multiply after injection into adult bees or pupae (18). This approach would necessarily exclude viruses that require different preparation methods or that are not amenable to crude extraction. Importantly, our use of metagenomic techniques was crucial in revealing a more complete bee virome, as many highly divergent viruses can be detected only at the sequence level (43, 46). Our study was also aided by the recent explosion in novel virus genomes (31–33, 35, 45), which provides a more comprehensive database for BLAST searches. As this database continues to expand, it is likely that more new viruses in a range of host species will be isolated. Finally, our detection of novel viruses may in part reflect our focus on bee colonies that are resistant to, or free of, Varroa, as the rapid spread of some viruses in the context of Varroa, most notably virulent variants of DWV, may have resulted in a general reduction in virus diversity. Clearly, the discovery of the new viruses here suggests that the bee virome will continue to expand following more extensive metagenomic surveys in diverse geographic regions.

MATERIALS AND METHODS

Sample collection.Honey bee colonies and mites from Europe, Africa, and the Pacific were sampled in 2013 to 2015. Seven A. mellifera capensis colonies from Robben Island and five colonies from mainland South Africa were sampled in March 2013. Varroa mites were first identified in South Africa in 1997, and after initial deleterious effects, A. mellifera capensis colonies exhibited natural mite resistance after 3 to 5 years (49). On Robben Island, Varroa mites were first detected 2 years prior to sampling (M. Allsopp, personal communication). In July 2014, 10 colonies of A. mellifera bees from an apiary at the Amsterdam Water Dunes, The Netherlands, where Varroa has been present since the 1980s and where natural selection for Varroa-resistant colonies has been ongoing since 2008, were sampled (48). In October 2015, feral and managed A. mellifera colonies from islands in the Kingdom of Tonga were sampled. Nine colonies from Vava'u island, where Varroa destructor was introduced in 2006, were sampled. Honey bees on Vava'u exhibit a natural tolerance that has enabled the survival of colonies in the decade following the introduction of the mites. Four colonies from Tongatapu island, where Varroa is not yet present, were sampled. In Tonga, 10 adult worker honey bee thoraxes per colony were crushed individually in 500 μl of RNAlater (Qiagen) and transported at room temperature prior to storage at −80°C until processing. In Africa and The Netherlands, a minimum of 10 adult worker honey bees were sampled per colony, frozen immediately on dry ice, and stored at −80°C until processing.

Sample processing.Thorax and abdomen tissues of five adult bees from each colony (thorax only from Tonga) were homogenized in TRIzol reagent (Thermo Fisher), and total RNA was extracted according to the manufacturer's protocol. RNA from each individual was diluted to 200 ng/μl, and 2.65 ng RNA from each of the five individuals from each colony was pooled prior to DNase treatment (Ambion), followed by column purification (RNeasy minikit; Qiagen). Total RNA was transported to the Australian Genome Research Facility (AGRF) on dry ice. Sample RNA integrity was confirmed by using a Bioanalyzer (Agilent). RNA was subjected to a ribosome depletion step (Ribo-Zero-Gold Human/Mouse/Rat) prior to the preparation of Illumina TruSeq Stranded Total RNA paired-end libraries according to the manufacturer's instructions. Libraries were run on an Illumina HiSeq2000 100-bp paired-end sequencing system for a total data yield of 4 to 9 Gb per sample.

Sequence read assembly and virus discovery.Sequencing reads were assembled de novo by using Trinity (80). The resulting contigs were compared to reference protein sequences of all previously characterized viruses downloaded from GenBank by using BLASTx (with an E value of 1E−5 to maximize sensitivity while minimizing false-positive results [35]). The resulting virus-like contigs were then compared to a nonredundant database using BLAST to remove nonviral hits, such as host contigs with similarity to viral sequences. We also removed any contigs with high similarity to plant viruses, which were more likely to be derived from food sources (although there is the possibility that these viruses could be replicating in bees [see reference 81]).

Virus sequences were aligned to sequences present in the current NCBI databases of homologous viral proteins by using MAFFT (82). Alignments were viewed, manually trimmed to remove large gaps and nonconserved regions, and further trimmed with TrimAL to remove ambiguously aligned regions (83). Maximum likelihood phylogenetic trees for each data set were inferred by using PhyML (84) in parallel mode, using Message Passaging Interface (MPI) with 12 threads and 12 random starting trees. We used a best-fit model of amino acid substitution determined by using ProtTest (85), a Subtree Pruning and Regrafting (SPR) branch-swapping algorithm, and an approximate likelihood ratio test (aLRT) with the Shimodaira-Hasegawa-like procedure to assess branch support.

Validation of novel viruses.We used PCR and sequencing to confirm the presence of each novel virus, designing primers based on the contigs assembled from next-generation sequencing. One virus, ANV, produced fragmented contigs spanning an incomplete genome. Each contig was ordered based on translated homology to the most closely related virus, Drosophila pseudoobscura Nora virus, and primers spanning contig breaks were designed to confirm the correct genome order and to sequence any unassembled regions (data not shown). We also confirmed the arrangement of the ARV-2 G and L proteins and the putative inter-virus-family horizontal transfer event, using PCR for amplification across the two ORFs and Sanger sequencing for confirmation.

We used Bowtie2 (86) to map reads to each of the novel virus genomes, Samtools to determine the sequencing depth and coverage (87), and the RSEM program implemented in Trinity to estimate the abundance of virus transcripts per million (TPM) (Table 3) (51). ORFs (Fig. 1) were annotated based on predicted amino acid sequences that were more than 200 nucleotides long as well as conserved positions in the genome compared to the most closely related viruses. Conserved domains were identified by using NCBI CDD BLAST searches (88).

Small RNA library preparation and sequencing.One microgram of total RNA (prepared as described above) from individual bees and Varroa mites was used to generate a small RNA library by using a TrueSeq small RNA kit (Illumina) according to the manufacturer's protocol. Samples were barcoded appropriately for pooling on an Illumina HiSeq2000 instrument, with 50-bp single-end sequencing.

Small RNA analysis.Small RNA reads were quality checked, trimmed to remove the TruSeq adapter, and then mapped by using CLC Genomics Workbench (Qiagen). The reads were first mapped to the A. mellifera genome, allowing for up to two mismatches (length fraction of 1 and similarity fraction of 0.9) due to the divergence between honey bee strains and the reference genome. The unmapped reads were subsequently mapped to ARV-1 and ARV-2 with the same stringency settings. The Amsterdam Water Dunes and Robben Island samples were mapped to the consensus viral genome generated by RNA sequencing (see above) from their geographic location. For further small RNA analysis (nucleotide size, genome position, base composition, and 5′ read distance), the mapped reads were exported as BAM files, indexed by using Samtools (87), and then analyzed in RStudio using viRome (89) and custom scripts.

Phasing analysis of small RNAs.To determine whether the small RNAs produced from the viral genomes showed evidence of Dicer phasing, we adapted an algorithm designed for detecting phased RNAs in plants (60). The logic behind this analysis is that if an RNA molecule is being cut every 22 nucleotides, there should be a relative accumulation of small RNA reads every 22 nucleotides. Thus, we summed the number of reads every x + 22 nucleotides along the sequence (Pi) (where the phase cycle position [x] equals 1 to 22) using custom R scripts and then divided this value by the number of remaining out-of-phase reads (U). Phasing scores were calculated for each phase cycle position in the 22-nt window over 8 cycles, as the signal is expected to degrade over time due to imprecision in the cut length. The relative accumulation of reads at a particular phase cycle position will give a high phasing score. For the honey bee data, we specifically looked for a phase cycle of 22, since this was the predominant size of the small RNAs. For the Varroa data, we analyzed four different phase cycle lengths (20, 22, 23, or 24 nucleotides) using the following equation: phasing score=ln[(1+10Σi=18Pi1 + ΣU)n − 2](1) where i is the phasing cycle, Pi is the number of small RNA reads at a given phase cycle position, U is the number of small RNA reads within the phase cycle not at the phase cycle position (out of phase), and n is the number of phase cycle positions occupied by at least one small RNA within the 8 cycles.

Accession number(s).All virus genome sequences generated in this study have been deposited in the GenBank database under accession numbers KY354230 to KY354234 and KY354236 to KY354244 . Annotated NCBI protein sequence accession numbers are ARO50020 to ARO50067 . Raw sequence data have been deposited in the Sequence Read Archive (SRA) (accession number SRP095071 ) (BioProject accession number PRJNA357165 and BioSample accession numbers SAMN06140203 to SAMN06140219 ).

ACKNOWLEDGMENTS

We thank Bruce White, Lamorna Osborne, Mike Allsopp, and Benjamin Oldroyd for assistance with sample collection. We thank John Roberts for useful discussions of the results. We acknowledge the Sydney Informatics Hub and The University of Sydney's high performance computing cluster Artemis for providing the high-performance computing resources that have contributed to the research results reported in this paper.

This work was funded by the Australian Research Council (grants DP170100844 to M.B. and DE140100199 to A.A.), the Australian & Pacific Science Foundation (grants APSF15-02 to E.J.R. and M.B.), and a Marie Bashir Institute for Infectious Diseases and Biosecurity seed funding grant (to E.J.R., E.C.H., and M.B.). E.C.H. is funded by an NHMRC Australia fellowship (GNT1037231). Natural selection for Varroa resistance in the Netherlands was funded by the Dutch Ministry of Economic Affairs and the European Commission as part of the “National Honey Programs” running from 2004 to date (grants to T.B.). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

FOOTNOTES

    • Received 27 January 2017.
    • Accepted 26 April 2017.
    • Accepted manuscript posted online 17 May 2017.
  • Supplemental material for this article may be found at https://doi.org/10.1128/JVI.00158-17 .

  • Copyright © 2017 Remnant et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license .

REFERENCES

  1. 1.↵
    1. Ellis JD,
    2. Munn PA
    . 2005. The worldwide health status of honey bees. Bee World86:88–101. doi:10.1080/0005772X.2005.11417323.
    OpenUrlCrossRef
  2. 2.↵
    1. Mondet F,
    2. de Miranda JR,
    3. Kretzschmar A,
    4. Le Conte Y,
    5. Mercer AR
    . 2014. On the front line: quantitative virus dynamics in honeybee (Apis mellifera L.) colonies along a new expansion front of the parasite Varroa destructor. PLoS Pathog10:e1004323. doi:10.1371/journal.ppat.1004323.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Martin SJ,
    2. Highfield AC,
    3. Brettell L,
    4. Villalobos EM,
    5. Budge GE,
    6. Powell M,
    7. Nikaido S,
    8. Schroeder D
    . 2012. Global honey bee viral landscape altered by a parasitic mite. Science336:1304–1306. doi:10.1126/science.1220941.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Evans JD,
    2. Schwarz RS
    . 2011. Bees brought to their knees: microbes affecting honey bee health. Trends Microbiol19:614–620. doi:10.1016/j.tim.2011.09.003.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    1. Runckel C,
    2. Flenniken ML,
    3. Engel JC,
    4. Ruby JG,
    5. Ganem D,
    6. Andino R,
    7. DeRisi JL
    . 2011. Temporal analysis of the honey bee microbiome reveals four novel viruses and seasonal prevalence of known viruses, Nosema, and Crithidia. PLoS One6:e20656. doi:10.1371/journal.pone.0020656.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Bailey L,
    2. Ball BV
    . 1991. Viruses, honey bee pathology, 2nd ed. Academic Press, London, United Kingdom.
  7. 7.↵
    1. Allen M,
    2. Ball B
    . 1996. The incidence and world distribution of honey bee viruses. Bee World77:141–162. doi:10.1080/0005772X.1996.11099306.
    OpenUrlCrossRefWeb of Science
  8. 8.↵
    1. Chen YP,
    2. Siede R
    . 2007. Honey bee viruses. Adv Virus Res70:33–80. doi:10.1016/S0065-3527(07)70002-7.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. de Miranda J,
    2. Cornman R,
    3. Evans J,
    4. Semberg E,
    5. Haddad N,
    6. Neumann P,
    7. Gauthier L
    . 2015. Genome characterization, prevalence and distribution of a macula-like virus from Apis mellifera and Varroa destructor. Viruses7:3586–3602. doi:10.3390/v7072789.
    OpenUrlCrossRef
  10. 10.↵
    1. Gauthier L,
    2. Cornman RS,
    3. Hartmann U,
    4. Cousserans F,
    5. Evans JD,
    6. de Miranda JR,
    7. Neumann P
    . 2015. The Apis mellifera filamentous virus genome. Viruses7:3798–3815. doi:10.3390/v7072798.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. de Miranda JR,
    2. Gauthier L,
    3. Ribiére M,
    4. Chen YP
    . 2012. Honey bee viruses and their effect on bee and colony health, p 71–102. InSammataro D, Yoder JA (ed), Honey bee colony health: challenges and sustainable solutions. CRC Press, Boca Raton, FL.
  12. 12.↵
    1. McMenamin AJ,
    2. Genersch E
    . 2015. Honey bee colony losses and associated viruses. Curr Opin Insect Sci8:121–129. doi:10.1016/j.cois.2015.01.015.
    OpenUrlCrossRef
  13. 13.↵
    1. Bailey L,
    2. Carpenter JM,
    3. Woods RD
    . 1979. Egypt bee virus and Australian isolates of Kashmir bee virus. J Gen Virol43:641–647. doi:10.1099/0022-1317-43-3-641.
    OpenUrlCrossRef
  14. 14.↵
    1. Bailey L,
    2. Woods RD
    . 1977. Two more small RNA viruses from honey bees and further observations on Sacbrood and acute bee-paralysis viruses. J Gen Virol37:175–182. doi:10.1099/0022-1317-37-1-175.
    OpenUrlCrossRef
  15. 15.↵
    1. Bromenshenk JJ,
    2. Henderson CB,
    3. Wick CH,
    4. Stanford MF,
    5. Zulich AW,
    6. Jabbour RE,
    7. Deshpande SV,
    8. McCubbin PE,
    9. Seccomb RA,
    10. Welch PM,
    11. Williams T,
    12. Firth DR,
    13. Skowronski E,
    14. Lehmann MM,
    15. Bilimoria SL,
    16. Gress J,
    17. Wanner KW,
    18. Cramer RA Jr
    . 2010. Iridovirus and microsporidian linked to honey bee colony decline. PLoS One5:e13181. doi:10.1371/journal.pone.0013181.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Baker A,
    2. Schroeder D
    . 2008. Occurrence and genetic analysis of picorna-like viruses infecting worker bees of Apis mellifera L. populations in Devon, South West England. J Invertebr Pathol98:239–242. doi:10.1016/j.jip.2008.02.010.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Baker A,
    2. Schroeder DC
    . 2008. The use of RNA-dependent RNA polymerase for the taxonomic assignment of picorna-like viruses (order Picornavirales) infecting Apis mellifera L. populations. Virol J5:10. doi:10.1186/1743-422X-5-10.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Ribière M,
    2. Ball B,
    3. Aubert MFA
    . 2008. Natural history and geographical distribution of honey bee viruses, p 19–25. InAubert M, Ball B, Fries I, Moritz R, Milani N, Bernardinelli I (ed), Virology and the honey bee. EEC Publications, Luxembourg City, Luxembourg.
  19. 19.↵
    1. Brutscher LM,
    2. McMenamin AJ,
    3. Flenniken ML
    . 2016. The buzz about honey bee viruses. PLoS Pathog12:e1005757. doi:10.1371/journal.ppat.1005757.
    OpenUrlCrossRef
  20. 20.↵
    1. Aubert M,
    2. Ball BV,
    3. Fries I,
    4. Moritz RFA,
    5. Milani N,
    6. Bernardinelli I
    (ed). 2008. Virology and the honey bee. EEC Publications, Luxembourg City, Luxembourg.
  21. 21.↵
    1. Le Gall O,
    2. Christian P,
    3. Fauquet CM,
    4. King AMQ,
    5. Knowles NJ,
    6. Nakashima N,
    7. Stanway G,
    8. Gorbalenya AE
    . 2008. Picornavirales, a proposed order of positive-sense single-stranded RNA viruses with a pseudo-T=3 virion architecture. Arch Virol153:715–727. doi:10.1007/s00705-008-0041-x.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    1. de Miranda JR,
    2. Genersch E
    . 2010. Deformed wing virus. J Invertebr Pathol103:S48–S61. doi:10.1016/j.jip.2009.06.012.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Mordecai GJ,
    2. Wilfert L,
    3. Martin SJ,
    4. Jones IM,
    5. Schroeder DC
    . 15November2015. Diversity in a honey bee pathogen: first report of a third master variant of the deformed wing virus quasispecies. ISME J doi:10.1038/ismej.2015.178.
    OpenUrlCrossRef
  24. 24.↵
    1. de Miranda JR,
    2. Cordoni G,
    3. Budge GE
    . 2010. The acute bee paralysis virus-Kashmir bee virus-Israeli acute paralysis virus complex. J Invertebr Pathol103:S30–S47. doi:10.1016/j.jip.2009.06.014.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Cornman RS,
    2. Tarpy DR,
    3. Chen Y,
    4. Jeffreys L,
    5. Lopez D,
    6. Pettis JS,
    7. vanEngelsdorp D,
    8. Evans JD
    . 2012. Pathogen webs in collapsing honey bee colonies. PLoS One7:e43562. doi:10.1371/journal.pone.0043562.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Daughenbaugh KF,
    2. Martin M,
    3. Brutscher LM,
    4. Cavigli I,
    5. Garcia E,
    6. Lavin M,
    7. Flenniken ML
    . 2015. Honey bee infecting Lake Sinai viruses. Viruses7:3285–3309. doi:10.3390/v7062772.
    OpenUrlCrossRef
  27. 27.↵
    1. Wilfert L,
    2. Long G,
    3. Leggett HC,
    4. Schmid-Hempel P,
    5. Butlin R,
    6. Martin SJM,
    7. Boots M
    . 2016. Deformed wing virus is a recent global epidemic in honeybees driven by Varroa mites. Science351:594–597. doi:10.1126/science.aac9976.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Ryabov EV,
    2. Wood GR,
    3. Fannon JM,
    4. Moore JD,
    5. Bull JC,
    6. Chandler D,
    7. Mead A,
    8. Burroughs N,
    9. Evans DJ
    . 2014. A virulent strain of deformed wing virus (DWV) of honeybees (Apis mellifera) prevails after Varroa destructor-mediated, or in vitro, transmission. PLoS Pathog10:e1004230. doi:10.1371/journal.ppat.1004230.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Genersch E,
    2. Aubert M
    . 2010. Emerging and re-emerging viruses of the honey bee (Apis mellifera L.). Vet Res41:54. doi:10.1051/vetres/2010027.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Rosenkranz P,
    2. Aumeier P,
    3. Ziegelmann B
    . 2010. Biology and control of Varroa destructor. J Invertebr Pathol103(Suppl):S96–S119. doi:10.1016/j.jip.2009.07.016.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Li C-X,
    2. Shi M,
    3. Tian J-H,
    4. Lin X-D,
    5. Kang Y-J,
    6. Chen L-J,
    7. Qin X-C,
    8. Xu J,
    9. Holmes EC,
    10. Zhang Y-Z
    . 2015. Unprecedented genomic diversity of RNA viruses in arthropods reveals the ancestry of negative-sense RNA viruses. eLife4:e05378. doi:10.7554/eLife.05378.
    OpenUrlCrossRef
  32. 32.↵
    1. Longdon B,
    2. Murray GGR,
    3. Palmer WJ,
    4. Day JP,
    5. Parker DJ,
    6. Welch JJ,
    7. Obbard DJ,
    8. Jiggins FM
    . 2015. The evolution, diversity, and host associations of rhabdoviruses. Virus Evol1:vev014. doi:10.1093/ve/vev014.
    OpenUrlCrossRef
  33. 33.↵
    1. Shi M,
    2. Lin X-D,
    3. Tian J-H,
    4. Chen L-J,
    5. Chen X,
    6. Li C-X,
    7. Qin X-C,
    8. Li J,
    9. Cao J-P,
    10. Eden J-S,
    11. Buchmann J,
    12. Wang W,
    13. Xu J,
    14. Holmes EC,
    15. Zhang Y-Z
    . 2016. Redefining the invertebrate RNA virosphere. Nature540:539–543. doi:10.1038/nature20167.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Webster CL,
    2. Longdon B,
    3. Lewis SH,
    4. Obbard DJ
    . 2016. Twenty-five new viruses associated with the Drosophilidae (Diptera). Evol Bioinform Online12:13–25. doi:10.4137/EBO.S39454.
    OpenUrlCrossRef
  35. 35.↵
    1. Shi M,
    2. Lin X-D,
    3. Vasilakis N,
    4. Tian J-H,
    5. Li C- X,
    6. Chen L-J,
    7. Eastwood G,
    8. Diao X-N,
    9. Chen M-H,
    10. Chen X,
    11. Qin X-C,
    12. Widen SG,
    13. Wood TG,
    14. Tesh RB,
    15. Xu J,
    16. Holmes EC,
    17. Zhang Y-Z
    . 2015. Divergent viruses discovered in arthropods and vertebrates revise the evolutionary history of the Flaviviridae and related viruses. J Virol90:659–669. doi:10.1128/JVI.02036-15.
    OpenUrlAbstract/FREE Full Text
  36. 36.↵
    1. Cook S,
    2. Moureau G,
    3. Kitchen A,
    4. Gould EA,
    5. de Lamballerie X,
    6. Holmes EC,
    7. Harbach RE
    . 2012. Molecular evolution of the insect-specific flaviviruses. J Gen Virol93:223–234. doi:10.1099/vir.0.036525-0.
    OpenUrlCrossRefPubMedWeb of Science
  37. 37.↵
    1. Ding S-W
    . 2010. RNA-based antiviral immunity. Nat Rev Immunol10:632–644. doi:10.1038/nri2824.
    OpenUrlCrossRefPubMedWeb of Science
  38. 38.↵
    1. Brutscher LM,
    2. Flenniken ML
    . 2015. RNAi and antiviral defense in the honey bee. J Immunol Res2015:941897. doi:10.1155/2015/941897.
    OpenUrlCrossRef
  39. 39.↵
    1. Li Y,
    2. Lu J,
    3. Han Y,
    4. Fan X,
    5. Ding S-W
    . 2013. RNA interference functions as an antiviral immunity mechanism in mammals. Science342:231–234. doi:10.1126/science.1241911.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    1. Chejanovsky N,
    2. Ophir R,
    3. Schwager MS,
    4. Slabezki Y,
    5. Grossman S,
    6. Cox-Foster D
    . 2014. Characterization of viral siRNA populations in honey bee colony collapse disorder. Virology454–455:176–183. doi:10.1016/j.virol.2014.02.012.
    OpenUrlCrossRef
  41. 41.↵
    1. Maori E,
    2. Paldi N,
    3. Shafir S,
    4. Kalev H,
    5. Tsur E,
    6. Glick E,
    7. Sela I
    . 2009. IAPV, a bee-affecting virus associated with colony collapse disorder can be silenced by dsRNA ingestion. Insect Mol Biol18:55–60. doi:10.1111/j.1365-2583.2009.00847.x.
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Liu X,
    2. Zhang Y,
    3. Yan X,
    4. Han R
    . 2010. Prevention of Chinese Sacbrood virus infection in Apis cerana using RNA interference. Curr Microbiol61:422–428. doi:10.1007/s00284-010-9633-2.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Liu S,
    2. Vijayendran D,
    3. Bonning BC
    . 2011. Next generation sequencing technologies for insect virus discovery. Viruses3:1849–1869. doi:10.3390/v3101849.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Tokarz R,
    2. Williams SH,
    3. Sameroff S,
    4. Sanchez Leon M,
    5. Jain K,
    6. Lipkin WI
    . 2014. Virome analysis of Amblyomma americanum, Dermacentor variabilis, and Ixodes scapularis ticks reveals novel highly divergent vertebrate and invertebrate viruses. J Virol88:11480–11492. doi:10.1128/JVI.01858-14.
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    1. Webster CL,
    2. Waldron FM,
    3. Robertson S,
    4. Crowson D,
    5. Ferrari G,
    6. Quintana JF,
    7. Brouqui J-M,
    8. Bayne EH,
    9. Longdon B,
    10. Buck AH,
    11. Lazzaro BP,
    12. Akorli J,
    13. Haddrill PR,
    14. Obbard DJ
    . 2015. The discovery, distribution, and evolution of viruses associated with Drosophila melanogaster. PLoS Biol13:e1002210. doi:10.1371/journal.pbio.1002210.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Mokili JL,
    2. Rohwer F,
    3. Dutilh BE
    . 2012. Metagenomics and future perspectives in virus discovery. Curr Opin Virol2:63–77. doi:10.1016/j.coviro.2011.12.004.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Locke B
    . 2015. Natural Varroa mite-surviving Apis mellifera honeybee populations. Apidologie47:467–482. doi:10.1007/s13592-015-0412-8.
    OpenUrlCrossRef
  48. 48.↵
    1. Blacquiere T
    . 2014. Natural selection for Varroa resistant honeybee colonies. Wageningen University & Research, Wageningen, The Netherlands. https://www.wur.nl/en/newsarticle/Natural-selection-for-varroa-resistant-honeybee-colonies.htm. Accessed 2 June 2014.
  49. 49.↵
    1. Allsopp M
    . 2006. Analysis of Varroa destructor infestation of southern African honeybee populations. MSc thesis. University of Pretoria, Pretoria, South Africa.
  50. 50.↵
    1. Taylor B,
    2. Roper T
    . 2006. Honey bee disease survey and beekeeper training in Tonga. AgriQuality Limited, Auckland, New Zealand.
  51. 51.↵
    1. Haas BJ,
    2. Papanicolaou A,
    3. Yassour M,
    4. Grabherr M,
    5. Blood PD,
    6. Bowden J,
    7. Couger MB,
    8. Eccles D,
    9. Li B,
    10. Lieber M,
    11. MacManes MD,
    12. Ott M,
    13. Orvis J,
    14. Pochet N,
    15. Strozzi F,
    16. Weeks N,
    17. Westerman R,
    18. William T,
    19. Dewey CN,
    20. Henschel R,
    21. LeDuc RD,
    22. Friedman N,
    23. Regev A
    . 2013. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc8:1494–1512. doi:10.1038/nprot.2013.084.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Schwarz RS,
    2. Bauchan GR,
    3. Murphy CA,
    4. Ravoet J,
    5. de Graaf DC,
    6. Evans JD
    . 2015. Characterization of two species of Trypanosomatidae from the honey bee Apis mellifera: Crithidia mellificae Langridge and McGhee, and Lotmaria passim n. gen., n. sp. J Eukaryot Microbiol62:567–583. doi:10.1111/jeu.12209.
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. Habayeb MS,
    2. Ekengren SK,
    3. Hultmark D
    . 2006. Nora virus, a persistent virus in Drosophila, defines a new picorna-like virus family. J Gen Virol87:3045–3051. doi:10.1099/vir.0.81997-0.
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. Walker PJ,
    2. Firth C,
    3. Widen SG,
    4. Blasdell KR,
    5. Guzman H,
    6. Wood TG,
    7. Paradkar PN,
    8. Holmes EC,
    9. Tesh RB,
    10. Vasilakis N
    . 2015. Evolution of genome size and complexity in the Rhabdoviridae. PLoS Pathog11:e1004664. doi:10.1371/journal.ppat.1004664.
    OpenUrlCrossRefPubMed
  55. 55.↵
    1. Palacios G,
    2. Forrester NL,
    3. Savij N,
    4. Travassos da Rosa APA,
    5. Guzman H,
    6. DeToy K,
    7. Popov VL,
    8. Walker PJ,
    9. Vasilakis N,
    10. Tesh RB
    . 2013. Characterization of Farmington virus, a novel virus from birds that is distantly related to members of the family Rhabdoviridae. Virol J10:219. doi:10.1186/1743-422X-10-219.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. Akopyants NS,
    2. Lye L-F,
    3. Dobson DE,
    4. Lukeš J,
    5. Beverley SM
    . 2016. A novel bunyavirus-like virus of trypanosomatid protist parasites. Genome Announc4:e00715-16. doi:10.1128/genomeA.00715-16.
    OpenUrlAbstract/FREE Full Text
  57. 57.↵
    1. Kobayashi K,
    2. Atsumi G,
    3. Iwadate Y,
    4. Tomita R,
    5. Chiba K,
    6. Akasaka S,
    7. Nishihara M,
    8. Takahashi H,
    9. Yamaoka N,
    10. Nishiguchi M,
    11. Sekine K-T
    . 2013. Gentian Kobu-sho-associated virus: a tentative, novel double-stranded RNA virus that is relevant to gentian Kobu-sho syndrome. J Gen Plant Pathol79:56–63. doi:10.1007/s10327-012-0423-5.
    OpenUrlCrossRef
  58. 58.↵
    1. Bekal S,
    2. Domier LL,
    3. Gonfa B,
    4. McCoppin NK,
    5. Lambert KN,
    6. Bhalerao K
    . 2014. A novel flavivirus in the soybean cyst nematode. J Gen Virol95:1272–1280. doi:10.1099/vir.0.060889-0.
    OpenUrlCrossRefPubMed
  59. 59.↵
    1. Gammon DB,
    2. Mello CC
    . 2015. RNA interference-mediated antiviral defense in insects. Curr Opin Insect Sci8:111–120. doi:10.1016/j.cois.2015.01.006.
    OpenUrlCrossRef
  60. 60.↵
    Anonymous. 2010. Genome sequencing and analysis of the model grass Brachypodium distachyon. Nature463:763–768. doi:10.1038/nature08747.
    OpenUrlCrossRefPubMedWeb of Science
  61. 61.↵
    1. Wilson RC,
    2. Doudna JA
    . 2013. Molecular mechanisms of RNA interference. Annu Rev Biophys42:217–239. doi:10.1146/annurev-biophys-083012-130404.
    OpenUrlCrossRefPubMed
  62. 62.↵
    1. Hess AM,
    2. Prasad AN,
    3. Ptitsyn A,
    4. Ebel GD,
    5. Olson KE,
    6. Barbacioru C,
    7. Monighetti C,
    8. Campbell CL
    . 2011. Small RNA profiling of Dengue virus-mosquito interactions implicates the PIWI RNA pathway in anti-viral defense. BMC Microbiol11:45. doi:10.1186/1471-2180-11-45.
    OpenUrlCrossRefPubMed
  63. 63.↵
    1. Morazzani EM,
    2. Wiley MR,
    3. Murreddu MG,
    4. Adelman ZN,
    5. Myles KM
    . 2012. Production of virus-derived ping-pong-dependent piRNA-like small RNAs in the mosquito soma. PLoS Pathog8:e1002470. doi:10.1371/journal.ppat.1002470.
    OpenUrlCrossRefPubMed
  64. 64.↵
    1. Gunawardane LS,
    2. Saito K,
    3. Nishida KM,
    4. Miyoshi K,
    5. Kawamura Y,
    6. Nagami T,
    7. Siomi H,
    8. Siomi MC
    . 2007. A Slicer-mediated mechanism for repeat-associated siRNA 5′ end formation in Drosophila. Science315:1587–1590. doi:10.1126/science.1140494.
    OpenUrlAbstract/FREE Full Text
  65. 65.↵
    1. King AMQ,
    2. Adams MJ,
    3. Carstens EB,
    4. Lefkowitz EJ
    (ed). 2011. Virus taxonomy. Classification and nomenclature of viruses. Ninth report of the International Committee on Taxonomy of Viruses. Elsevier Academic Press, San Diego, CA.
  66. 66.↵
    1. Marklewitz M,
    2. Zirkel F,
    3. Kurth A,
    4. Drosten C,
    5. Junglen S
    . 2015. Evolutionary and phenotypic analysis of live virus isolates suggests arthropod origin of a pathogenic RNA virus family. Proc Natl Acad Sci U S A112:7536–7541. doi:10.1073/pnas.1502036112.
    OpenUrlAbstract/FREE Full Text
  67. 67.↵
    1. Schoonvaere K,
    2. De Smet L,
    3. Smagghe G,
    4. Vierstraete A,
    5. Braeckman BP,
    6. de Graaf DC
    . 2016. Unbiased RNA shotgun metagenomics in social and solitary wild bees detects associations with eukaryote parasites and new viruses. PLoS One11:e0168456. doi:10.1371/journal.pone.0168456.
    OpenUrlCrossRef
  68. 68.↵
    1. Yurchenko V,
    2. Votypka J,
    3. Tesarova M,
    4. Klepetkova H,
    5. Kraeva N,
    6. Jirku M,
    7. Lukes J
    . 2014. Ultrastructure and molecular phylogeny of four new species of monoxenous trypanosomatids from flies (Diptera: Brachycera) with redefinition of the genus Wallaceina. Folia Parasitol (Praha)61:97–112.
    OpenUrl
  69. 69.↵
    1. Ashe A,
    2. Bélicard T,
    3. Le Pen J,
    4. Sarkies P,
    5. Frézal L,
    6. Lehrbach NJ,
    7. Félix M-A,
    8. Miska EA
    . 2013. A deletion polymorphism in the Caenorhabditis elegans RIG-I homolog disables viral RNA dicing and antiviral immunity. eLife2:e00994. doi:10.7554/eLife.00994.
    OpenUrlCrossRefPubMed
  70. 70.↵
    1. Sabin LR,
    2. Zheng Q,
    3. Thekkat P,
    4. Yang J,
    5. Hannon GJ,
    6. Gregory BD,
    7. Tudor M,
    8. Cherry S
    . 2013. Dicer-2 processes diverse viral RNA species. PLoS One8:e55458. doi:10.1371/journal.pone.0055458.
    OpenUrlCrossRefPubMed
  71. 71.↵
    1. Palmer WJ,
    2. Jiggins FM
    . 2015. Comparative genomics reveals the origins and diversity of arthropod immune systems. Mol Biol Evol32:2111–2129. doi:10.1093/molbev/msv093.
    OpenUrlCrossRefPubMed
  72. 72.↵
    1. Jayachandran B,
    2. Hussain M,
    3. Asgari S
    . 2012. RNA interference as a cellular defense mechanism against the DNA virus baculovirus. J Virol86:13729–13734. doi:10.1128/JVI.02041-12.
    OpenUrlAbstract/FREE Full Text
  73. 73.↵
    1. Christian PD,
    2. Scotti PD
    . 1998. Picornalike viruses of insects, p 301–336. InMiller LK, Ball LA (ed), The insect viruses. Plenum Publishing Corporation, New York, NY.
  74. 74.↵
    1. Zhang X,
    2. He SY,
    3. Evans JD,
    4. Pettis JS,
    5. Yin GF,
    6. Chen YP
    . 2012. New evidence that deformed wing virus and black queen cell virus are multi-host pathogens. J Invertebr Pathol109:156–159. doi:10.1016/j.jip.2011.09.010.
    OpenUrlCrossRefPubMed
  75. 75.↵
    1. Singh R,
    2. Levitt AL,
    3. Rajotte EG,
    4. Holmes EC,
    5. Ostiguy N,
    6. vanEngelsdorp D,
    7. Lipkin WI,
    8. dePamphillis CW,
    9. Toth AL,
    10. Cox-Foster DL
    . 2010. RNA viruses in hymenopteran pollinators: evidence of inter-taxa virus transmission via pollen and potential impact on non-Apis hymenopteran species. PLoS One5:e14357. doi:10.1371/journal.pone.0014357.
    OpenUrlCrossRefPubMed
  76. 76.↵
    1. McMahon DP,
    2. Fürst MA,
    3. Caspar J,
    4. Theodorou P,
    5. Brown MJF,
    6. Paxton RJ
    . 2015. A sting in the spit: widespread cross-infection of multiple RNA viruses across wild and managed bees. J Anim Ecol84:615–624. doi:10.1111/1365-2656.12345.
    OpenUrlCrossRefPubMed
  77. 77.↵
    1. Koonin EV,
    2. Dolja VV,
    3. Krupovic M
    . 2015. Origins and evolution of viruses of eukaryotes: the ultimate modularity. Virology479–480:2–25. doi:10.1016/j.virol.2015.02.039.
    OpenUrlCrossRefPubMed
  78. 78.↵
    1. Anderson DL,
    2. Gibbs AJ
    . 1988. Inapparent virus infections and their interactions in pupae of the honey bee (Apis mellifera Linnaeus) in Australia. J Gen Virol69:1617–1625. doi:10.1099/0022-1317-69-7-1617.
    OpenUrlCrossRefWeb of Science
  79. 79.↵
    1. Plus N,
    2. Scotti PD
    . 1984. The biological properties of eight different isolates of cricket paralysis virus. Ann Virol (Inst Pasteur)135E:257–268.
    OpenUrl
  80. 80.↵
    1. Grabherr MG,
    2. Haas BJ,
    3. Yassour M,
    4. Levin JZ,
    5. Thompson DA,
    6. Amit I,
    7. Adiconis X,
    8. Fan L,
    9. Raychowdhury R,
    10. Zeng Q,
    11. Chen Z,
    12. Mauceli E,
    13. Hacohen N,
    14. Gnirke A,
    15. Rhind N,
    16. di Palma F,
    17. Birren BW,
    18. Nusbaum C,
    19. Lindblad-Toh K,
    20. Friedman N,
    21. Regev A
    . 2011. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol29:644–652. doi:10.1038/nbt.1883.
    OpenUrlCrossRefPubMed
  81. 81.↵
    1. Li JL,
    2. Cornman RS,
    3. Evans JD,
    4. Pettis JS,
    5. Zhao Y,
    6. Murphy C,
    7. Peng WJ,
    8. Wu J,
    9. Hamilton M,
    10. Boncristiani HF,
    11. Zhou L,
    12. Hammond J,
    13. Chen YP
    . 2014. Systemic spread and propagation of a plant-pathogenic virus in European honeybees, Apis mellifera. mBio5:e00898-13. doi:10.1128/mBio.00898-13.
    OpenUrlAbstract/FREE Full Text
  82. 82.↵
    1. Katoh K,
    2. Standley DM
    . 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol30:772–780. doi:10.1093/molbev/mst010.
    OpenUrlCrossRefPubMedWeb of Science
  83. 83.↵
    1. Capella-Gutiérrez S,
    2. Silla-Martínez JM,
    3. Gabaldón T
    . 2009. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics25:1972–1973. doi:10.1093/bioinformatics/btp348.
    OpenUrlCrossRefPubMedWeb of Science
  84. 84.↵
    1. Guindon S,
    2. Dufayard JF,
    3. Lefort V,
    4. Anisimova M,
    5. Hordijk W,
    6. Gascuel O
    . 2010. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol59:307–321. doi:10.1093/sysbio/syq010.
    OpenUrlCrossRefPubMedWeb of Science
  85. 85.↵
    1. Darriba D,
    2. Taboada GL,
    3. Doallo R,
    4. Posada D
    . 17February2011. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics doi:10.1093/bioinformatics/btr088.
    OpenUrlCrossRefPubMedWeb of Science
  86. 86.↵
    1. Langmead B,
    2. Salzberg SL
    . 2012. Fast gapped-read alignment with Bowtie 2. Nat Methods9:357–359. doi:10.1038/nmeth.1923.
    OpenUrlCrossRefPubMedWeb of Science
  87. 87.↵
    1. Li H,
    2. Handsaker B,
    3. Wysoker A,
    4. Fennell T,
    5. Ruan J,
    6. Homer N,
    7. Marth G,
    8. Abecasis G,
    9. Durbin R
    . 2009. The Sequence Alignment/Map format and SAMtools. Bioinformatics25:2078–2079. doi:10.1093/bioinformatics/btp352.
    OpenUrlCrossRefPubMedWeb of Science
  88. 88.↵
    1. Marchler-Bauer A,
    2. Derbyshire MK,
    3. Gonzales NR,
    4. Lu S,
    5. Chitsaz F,
    6. Geer LY,
    7. Geer RC,
    8. He J,
    9. Gwadz M,
    10. Hurwitz DI,
    11. Lanczycki CJ,
    12. Lu F,
    13. Marchler GH,
    14. Song JS,
    15. Thanki N,
    16. Wang Z,
    17. Yamashita RA,
    18. Zhang D,
    19. Zheng C,
    20. Bryant SH
    . 2015. CDD: NCBI's conserved domain database. Nucleic Acids Res43:D222–D226. doi:10.1093/nar/gku1221.
    OpenUrlCrossRefPubMed
  89. 89.↵
    1. Watson M,
    2. Schnettler E,
    3. Kohl A
    . 2013. viRome: an R package for the visualization and analysis of viral small RNA sequence datasets. Bioinformatics29:1902–1903. doi:10.1093/bioinformatics/btt297.
    OpenUrlCrossRefPubMedWeb of Science
PreviousNext
Back to top
Download PDF
Citation Tools
A Diverse Range of Novel RNA Viruses in Geographically Distinct Honey Bee Populations
Emily J. Remnant, Mang Shi, Gabriele Buchmann, Tjeerd Blacquière, Edward C. Holmes, Madeleine Beekman, Alyson Ashe
Journal of Virology Jul 2017, 91 (16) e00158-17; DOI: 10.1128/JVI.00158-17

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this Journal of Virology article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
A Diverse Range of Novel RNA Viruses in Geographically Distinct Honey Bee Populations
(Your Name) has forwarded a page to you from Journal of Virology
(Your Name) thought you would be interested in this article in Journal of Virology.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
A Diverse Range of Novel RNA Viruses in Geographically Distinct Honey Bee Populations
Emily J. Remnant, Mang Shi, Gabriele Buchmann, Tjeerd Blacquière, Edward C. Holmes, Madeleine Beekman, Alyson Ashe
Journal of Virology Jul 2017, 91 (16) e00158-17; DOI: 10.1128/JVI.00158-17
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Bees
Biodiversity
genetic variation
RNA viruses
RNA interference
arthropod vectors
insect viruses
metagenomics
negative-strand RNA virus
phylogenetic analysis
plus-strand RNA virus

Related Articles

Cited By...

About

  • About JVI
  • Editor in Chief
  • Editorial Board
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Article Types
  • Ethics
  • Contact Us

Follow #Jvirology

@ASMicrobiology

       

 

JVI in collaboration with

American Society for Virology

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Print ISSN: 0022-538X; Online ISSN: 1098-5514