RT-PCR analysis of Deformed wing virus in honeybees (Apis mellifera) and mites (Varroa destructor)

Short reads from honey bee (Apis sp.) sequencing projects reflect microbial associate diversity

High throughput (or ‘next generation’) sequencing has transformed most areas of biological research and is now a standard method that underpins empirical study of organismal biology, and (through comparison of genomes), reveals patterns of evolution. For projects focused on animals, these sequencing methods do not discriminate between the primary target of sequencing (the animal genome) and ‘contaminating’ material, such as associated microbes. A common first step is to filter out these contaminants to allow better assembly of the animal genome or transcriptome. Here, we aimed to assess if these ‘contaminations’ provide information with regard to biologically important microorganisms associated with the individual. To achieve this, we examined whether the short read data from Apis retrieved elements of its well established microbiome. To this end, we screened almost 1,000 short read libraries of honey bee (Apis sp.) DNA sequencing project for the presence of microbial sequences, and find sequences from known honey bee microbial associates in at least 11% of them. Further to this, we screened ∼500 Apis RNA sequencing libraries for evidence of viral infections, which were found to be present in about half of them. We then used the data to reconstruct draft genomes of three Apis associated bacteria, as well as several viral strains de novo. We conclude that ‘contamination’ in short read sequencing libraries can provide useful genomic information on microbial taxa known to be associated with the target organisms, and may even lead to the discovery of novel associations. Finally, we demonstrate that RNAseq samples from experiments commonly carry uneven viral loads across libraries. We note variation in viral presence and load may be a confounding feature of differential gene expression analyses, and as such it should be incorporated as a random factor in analyses.

Persistence of subclinical deformed wing virus infections in honeybees following Varroa mite removal and a bee population turnover

Deformed wing virus (DWV) is a lethal virus of honeybees (Apis mellifera) implicated in elevated colony mortality rates worldwide and facilitated through vector transmission by the ectoparasitic mite Varroa destructor. Clinical, symptomatic DWV infections are almost exclusively associated with high virus titres during pupal development, usually acquired through feeding by Varroa mites when reproducing on bee pupae. Control of the mite population, generally through acaricide treatment, is essential for breaking the DWV epidemic and minimizing colony losses. In this study, we evaluated the effectiveness of remedial mite control on clearing DWV from a colony. DWV titres in adult bees and pupae were monitored at 2 week intervals through summer and autumn in acaricide-treated and untreated colonies. The DWV titres in Apistan treated colonies was reduced 1000-fold relative to untreated colonies, which coincided with both the removal of mites and also a turnover of the bee population in the colony. This adult bee population turnover is probably more critical than previously realized for effective clearing of DWV infections. After this initial reduction, subclinical DWV titres persisted and even increased again gradually during autumn, demonstrating that alternative non-Varroa transmission routes can maintain the DWV titres at significant subclinical levels even after mite removal. The implications of these results for practical recommendations to mitigate deleterious subclinical DWV infections and improving honeybee health management are discussed.

Single-stranded RNA viruses infecting the invasive Argentine ant, Linepithema humile

Social insects host a diversity of viruses. We examined New Zealand populations of the globally widely distributed invasive Argentine ant (Linepithema humile) for RNA viruses. We used metatranscriptomic analysis, which identified six potential novel viruses in the Dicistroviridae family. Of these, three contigs were confirmed by Sanger sequencing as Linepithema humile virus-1 (LHUV-1), a novel strain of Kashmir bee virus (KBV) and Black queen cell virus (BQCV), while the others were chimeric or misassembled sequences. We extended the known sequence of LHUV-1 to confirm its placement in the Dicistroviridae and categorised its relationship to closest relatives, which were all viruses infecting Hymenoptera. We examined further for known viruses by mapping our metatranscriptomic sequences to all viral genomes, and confirmed KBV, BQCV, LHUV-1 and Deformed wing virus (DWV) presence using qRT-PCR. Viral replication was confirmed for DWV, KBV and LHUV-1. Viral titers in ants were higher in the presence of honey bee hives. Argentine ants appear to host a range of' honey bee' pathogens in addition to a virus currently described only from this invasive ant. The role of these viruses in the population dynamics of the ant remain to be determined, but offer potential targets for biocontrol approaches.

A Comparison of Wolbachia Infection Frequencies in Varroa With Prevalence of Deformed Wing Virus

Wolbachia are widely distributed bacterial endosymbionts of arthropods and filarial nematodes. These bacteria can affect host fitness in a variety of ways, such as protecting hosts against viruses and other pathogens. Here, we investigate the possible role of Wolbachia in the prevalence of the deformed wing virus (DWV), a highly virulent pathogen of honey bees (Apis mellifera) that is transmitted by parasitic Varroa mites (Varroa destructor). About 180 Varroa mites from 18 beehives were tested for infection with Wolbachia and DWV. We first screened for Wolbachia using two standard primers (wsp and 16S rDNA), and found 26% of the mites to be positive for Wolbachia using the wsp primer and 64% of the mites to be positive using the 16S rDNA primer. Using these intermediate Wolbachia frequencies, we then tested for statistical correlations with virus infection frequencies. The analysis revealed a significant positive correlation between DWV and Wolbachia using the wsp primer, but no significant association between DWV and Wolbachia using the 16S rDNA primer. In conclusion, there is no evidence for an anti-pathogenic effect of Wolbachia in V. destructor, but weak evidence for a pro-pathogenic effect. These results encourage further examination of Wolbachia-virus interactions in Varroa mites since an increased vector competence of the mites may significantly impact disease outbreaks in honey bees.

Prevalence of honey bee (Apis mellifera) viruses in temperate and subtropical regions from Argentina

In Argentina, bee virus studies are still incipient, and there are no studies regarding the climatic effect. The aim of this study was to assess and compare the presence of honeybee viruses in different climatic regions from Argentina. A total of 385 colonies distributed in five Argentinean eco-regions were examined to evaluate the percentage of infestation with Varroa destructor and the presence of seven virus species (Deformed wing virus, DWV; Acute bee paralysis virus, ABPV; Chronic bee paralysis virus, CBPV; Black queen cell virus, BQCV; Kashmer bee virus, KBV; Israeli acute bee paralysis virus, IAPV; and Sacbrood bee virus, SBV) after honey yield. Two viruses, KBV and IAPV, were not detected. The other five viruses were found in different prevalences: DWV (35%), ABPV (21.5%), BQCV (8.0%), CBPV (2.2%), and SBV (1.1%). We found double and triple viral associations in approximately 25% of the sampled colonies. The mean V. destructor infestation in the colonies prior to the acaricide treatment was 7.12%±8.7%. The knowledge of the prevalence of these viruses in the region and their relation with the mite and other possible influencing factors is important for preventing colony losses. Further studies are necessary to identify the risk factors associated with virus presence and its relationship with other pathogens such as V. destructor.

Replication of honey bee-associated RNA viruses across multiple bee species in apple orchards of Georgia, Germany and Kyrgyzstan

The essential ecosystem service of pollination is provided largely by insects, which are considered threatened by diverse biotic and abiotic global change pressures. RNA viruses are one such pressure, and have risen in prominence as a major threat for honey bees (Apis mellifera) and global apiculture, as well as a risk factor for other bee species through pathogen spill-over between managed honey bees and sympatric wild pollinator communities. Yet despite their potential role in global bee decline, the prevalence of honey bee-associated RNA viruses in wild bees is poorly known from both geographic and taxonomic perspectives. We screened members of pollinator communities (honey bees, bumble bees and other wild bees belonging to four families) collected from apple orchards in Georgia, Germany and Kyrgyzstan for six common honey bee-associated RNA virus complexes encompassing nine virus targets. The Deformed wing virus complex (DWV genotypes A and B) had the highest prevalence across all localities and host species and was the only virus complex found in wild bee species belonging to all four studied families. Based on amplification of negative-strand viral RNA, we found evidence for viral replication in wild bee species of DWV-A/DWV-B (hosts: Andrena haemorrhoa and several Bombus spp.) and Black queen cell virus (hosts: Anthophora plumipes, several Bombus spp., Osmia bicornis and Xylocopa spp.). Viral amplicon sequences revealed that DWV-A and DWV-B are regionally distinct but identical in two or more bee species at any one site, suggesting virus is shared amongst sympatric bee taxa. This study demonstrates that honey bee associated RNA viruses are geographically and taxonomically widespread, likely infective in wild bee species, and shared across bee taxa.

A Comparison of Deformed Wing Virus in Deformed and Asymptomatic Honey Bees

Deformed wing virus (DWV) in association with Varroa destructor is currently attributed to being responsible for colony collapse in the western honey bee (Apis mellifera). The appearance of deformed individuals within an infested colony has long been associated with colony losses. However, it is unknown why only a fraction of DWV positive bees develop deformed wings. This study concerns two small studies comparing deformed and non-deformed bees. In Brazil, asymptomatic bees (no wing deformity) that had been parasitised by Varroa as pupae had higher DWV loads than non-parasitised bees. However, we found no greater bilateral asymmetry in wing morphology due to DWV titres or parasitisation. As expected, using RT-qPCR, deformed bees were found to contain the highest viral loads. In a separate study, next generation sequencing (NGS) was applied to compare the entire DWV genomes from paired symptomatic and asymptomatic bees from three colonies on two different Hawaiian islands. This revealed no consistent differences between DWV genomes from deformed or asymptomatic bees, with the greatest variation seen between locations, not phenotypes. All samples, except one, were dominated by DWV type A. This small-scale study suggests that there is no unique genetic variant associated with wing deformity; but that many DWV variants have the potential to cause deformity.

Covert deformed wing virus infections have long-term deleterious effects on honeybee foraging and survival

Several studies have suggested that covert stressors can contribute to bee colony declines. Here we provide a novel case study and show using radiofrequency identification tracking technology that covert deformed wing virus (DWV) infections in adult honeybee workers seriously impact long-term foraging and survival under natural foraging conditions. In particular, our experiments show that adult workers injected with low doses of DWV experienced increased mortality rates, that DWV caused workers to start foraging at a premature age, and that the virus reduced the workers' total activity span as foragers. Altogether, these results demonstrate that covert DWV infections have strongly deleterious effects on honeybee foraging and survival. These results are consistent with previous studies that suggested DWV to be an important contributor to the ongoing bee declines in Europe and the USA. Overall, our study underlines the strong impact that covert pathogen infections can have on individual and group-level performance in bees.

Characterization of a novel member of genus Iflavirus in Helicoverpa armigera

The cotton bollworm, Helicoverpa armigera, is one of the most important agricultural pests of many economic crops worldwide. Herein, we found a novel single-strand RNA virus by RNA-Seq and Polymerase Chain Reaction (PCR) method in H. armigera named Helicoverpa armigera iflavirus (HaIV), which possessed a genome with 10,017 nucleotides in length and contained a single large open reading frame (ORF) encoding a putative polyprotein of 3021 amino acids with a predicted molecular mass of 344.16kDa and a theoretical isoelectric point (pI) of 6.45. The deduced amino acid sequence showed highest similarity (61.0%) with the protein of Lymantria dispar Iflavirus 1. Phylogenetic analysis with putative RdRp amino acid sequences indicated that the virus clustered with members of the genus Iflavirus. The virus was mainly distributed in the fat body of its host and was found to be capable of both horizontal and vertical transmission. The efficiency of perorally horizontal transmission was dose dependent (100% infection rate with a viral dose of 10⁸copies/μl) while vertical transmission efficiency was found to be relatively low (<28.57%). These results suggest that we have found a novel member of genus Iflavirus in H. armigera.

Distribution and variability of deformed wing virus of honeybees (Apis mellifera) in the Middle East and North Africa

Three hundred and eleven honeybee samples from 12 countries in the Middle East and North Africa (MENA) (Jordan, Lebanon, Syria, Iraq, Egypt, Libya, Tunisia, Algeria, Morocco, Yemen, Palestine, and Sudan) were analyzed for the presence of deformed wing virus (DWV). The prevalence of DWV throughout the MENA region was pervasive, but variable. The highest prevalence was found in Lebanon and Syria, with prevalence dropping in Palestine, Jordan, and Egypt before increasing slightly moving westwards to Algeria and Morocco Phylogenetic analysis of a 194 nucleotide section of the DWV Lp gene did not identify any significant phylogenetic resolution among the samples, although the sequences did show consistent regional clustering, including an interesting geographic gradient from Morocco through North Africa to Jordan and Syria. The sequences revealed several clear variability hotspots in the deduced amino acid sequence, which furthermore showed some patterns of regional identity. Furthermore, the sequence variants from the Middle East and North Africa appear more numerous and diverse than those from Europe.

Importance of Ecological Factors and Colony Handling for Optimizing Health Status of Apiaries in Mediterranean Ecosystems

We analyzed six apiaries in several natural environments with a Mediterranean ecosystem in Madrid, central Spain, in order to understand how landscape and management characteristics may influence apiary health and bee production in the long term. We focused on five criteria (habitat quality, landscape heterogeneity, climate, management and health), as well as 30 subcriteria, and we used the analytic hierarchy process (AHP) to rank them according to relevance. Habitat quality proved to have the highest relevance, followed by beehive management. Within habitat quality, the following subcriteria proved to be most relevant: orographic diversity, elevation range and important plant species located 1.5 km from the apiary. The most important subcriteria under beehive management were honey production, movement of the apiary to a location with a higher altitude and wax renewal. Temperature was the most important subcriterion under climate, while pathogen and Varroa loads were the most significant under health. Two of the six apiaries showed the best values in the AHP analysis and showed annual honey production of 70 and 28 kg/colony. This high productivity was due primarily to high elevation range and high orographic diversity, which favored high habitat quality. In addition, one of these apiaries showed the best value for beehive management, while the other showed the best value for health, reflected in the low pathogen load and low average number of viruses. These results highlight the importance of environmental factors and good sanitary practices to maximize apiary health and honey productivity.

Flight performance of actively foraging honey bees is reduced by a common pathogen

Sudden and severe declines in honey bee (Apis mellifera) colony health in the US and Europe have been attributed, in part, to emergent microbial pathogens, however, the mechanisms behind the impact are unclear. Using roundabout flight mills, we measured the flight distance and duration of actively foraging, healthy-looking honey bees sampled from standard colonies, before quantifying the level of infection by Nosema ceranae and Deformed Wing Virus complex (DWV) for each bee. Neither the presence nor the quantity of N. ceranae were at low, natural levels of infection had any effect on flight distance or duration, but presence of DWV reduced flight distance by two thirds and duration by one half. Quantity of DWV was shown to have a significant, but weakly positive relation with flight distance and duration, however, the low amount of variation that was accounted for suggests further investigation by dose-response assays is required. We conclude that widespread, naturally occurring levels of infection by DWV weaken the flight ability of honey bees and high levels of within-colony prevalence are likely to reduce efficiency and increase the cost of resource acquisition. Predictions of implications of pathogens on colony health and function should take account of sublethal effects on flight performance.

Deformed wing virus can be transmitted during natural mating in honey bees and infect the queens

Deformed wing virus is an important contributor to honey bee colony losses. Frequently queen failure is reported as a cause for colony loss. Here we examine whether sexual transmission during multiple matings of queens is a possible way of virus infection in queens. In an environment with high prevalence of deformed wing virus, queens (n = 30) were trapped upon their return from natural mating flights. The last drone's endophallus (n = 29), if present, was removed from the mated queens for deformed wing virus quantification, leading to the detection of high-level infection in 3 endophalli. After oviposition, viral quantification revealed that seven of the 30 queens had high-level deformed wing virus infections, in all tissues, including the semen stored in the spermathecae. Two groups of either unmated queens (n = 8) with induced egg laying, or queens (n = 12) mated in isolation with drones showing comparatively low deformed wing virus infections served as control. None of the control queens exhibited high-level viral infections. Our results demonstrate that deformed wing virus infected drones are competitive to mate and able to transmit the virus along with semen, which occasionally leads to queen infections. Virus transmission to queens during mating may be common and can contribute noticeably to queen failure.

Pathogenesis of varroosis at the level of the honey bee (Apis mellifera) colony

The parasitic mite Varroa destructor, in interaction with different viruses, is the main cause of honey bee colony mortality in most parts of the world. Here we studied how effects of individual-level parasitization are reflected by the bee colony as a whole. We measured disease progression in an apiary of 24 hives with differing degree of mite infestation, and investigated its relationship to 28 biometrical, physiological and biochemical indicators. In early summer, when the most heavily infested colonies already showed reduced growth, an elevated ratio of brood to bees, as well as a strong presence of phenoloxidase/prophenoloxidase in hive bees were found to be predictors of the time of colony collapse. One month later, the learning performance of worker bees as well as the activity of glucose oxidase measured from head extracts were significantly linked to the timing of colony collapse. Colonies at the brink of collapse were characterized by reduced weight of winter bees and a strong increase in their relative body water content. Our data confirm the importance of the immune system, known from studies of individually-infested bees, for the pathogenesis of varroosis at colony level. However, they also show that single-bee effects cannot always be extrapolated to the colony as a whole. This fact, together with the prominent role of colony-level factors like the ratio between brood and bees for disease progression, stress the importance of the superorganismal dimension of Varroa research.

Elevated virulence of an emerging viral genotype as a driver of honeybee loss

Emerging infectious diseases (EIDs) have contributed significantly to the current biodiversity crisis, leading to widespread epidemics and population loss. Owing to genetic variation in pathogen virulence, a complete understanding of species decline requires the accurate identification and characterization of EIDs. We explore this issue in the Western honeybee, where increasing mortality of populations in the Northern Hemisphere has caused major concern. Specifically, we investigate the importance of genetic identity of the main suspect in mortality, deformed wing virus (DWV), in driving honeybee loss. Using laboratory experiments and a systematic field survey, we demonstrate that an emerging DWV genotype (DWV-B) is more virulent than the established DWV genotype (DWV-A) and is widespread in the landscape. Furthermore, we show in a simple model that colonies infected with DWV-B collapse sooner than colonies infected with DWV-A. We also identify potential for rapid DWV evolution by revealing extensive genome-wide recombination in vivo The emergence of DWV-B in naive honeybee populations, including via recombination with DWV-A, could be of significant ecological and economic importance. Our findings emphasize that knowledge of pathogen genetic identity and diversity is critical to understanding drivers of species decline.

Varroa destructor (Mesostigmata: Varroidae) Parasitism and Climate Differentially Influence the Prevalence, Levels, and Overt Infections of Deformed Wing Virus in Honey Bees (Hymenoptera: Apidae)

The prevalence and loads of deformed wing virus (DWV) between honey bee (Apis mellifera L.) colonies from a tropical and a temperate environment were compared. The interaction between these environments and the mite Varroa destructor in relation to DWV prevalence, levels, and overt infections, was also analyzed. V. destructor rates were determined, and samples of mites, adult bees, brood parasitized with varroa mites and brood not infested by mites were analyzed. DWV was detected in 100% of the mites and its prevalence and loads in honey bees were significantly higher in colonies from the temperate climate than in colonies from the tropical climate. Significant interactions were found between climate and type of sample, with the highest levels of DWV found in varroa-parasitized brood from temperate climate colonies. Additionally, overt infections were observed only in the temperate climate. Varroa parasitism and DWV loads in bees from colonies with overt infections were significantly higher than in bees from colonies with covert infections. These results suggest that interactions between climate, V. destructor, and possibly other factors, may play a significant role in the prevalence and levels of DWV in honey bee colonies, as well as in the development of overt infections. Several hypotheses are discussed to explain these results.

Bay laurel (Laurus nobilis) as potential antiviral treatment in naturally BQCV infected honeybees

Viral diseases are one of the multiple factors associated with honeybee colony losses. Apart from their innate immune system, including the RNAi machinery, honeybees can use secondary plant metabolites to reduce or fully cure pathogen infections. Here, we tested the antiviral potential of Laurus nobilis leaf ethanolic extracts on forager honeybees naturally infected with BQCV (Black queen cell virus). Total viral loads were reduced even at the lowest concentration tested (1mg/ml). Higher extract concentrations (≥5mg/ml) significantly reduced virus replication. Measuring vitellogenin gene expression as an indicator for transcript homeostasis revealed constant RNA levels before and after treatment, suggesting that its expression was not impacted by the L. nobilis treatment. In conclusion, plant secondary metabolites can reduce virus loads and virus replication in naturally infected honeybees.

A Toolbox for Quantitative Gene Expression in Varroa destructor: RNA Degradation in Field Samples and Systematic Analysis of Reference Gene Stability

Varroa destructor is the major pest of Apis mellifera and contributes to the global honey bee health crisis threatening food security. Developing new control strategies to combat Varroa will require the application of molecular biology, including gene expression studies by quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR). Both high quality RNA samples and suitable stable internal reference genes are required for accurate gene expression studies. In this study, ten candidate genes (succinate dehydrogenase (SDHA), NADH dehydrogenase (NADH), large ribsosmal subunit, TATA-binding protein, glyceraldehyde-3-phosphate dehydrogenase, 18S rRNA (18S), heat-shock protein 90 (HSP90), cyclophilin, α-tubulin, actin), were evaluated for their suitability as normalization genes using the geNorm, Normfinder, BestKeeper, and comparative ΔCq algorithims. Our study proposes the use of no more than two of the four most stable reference genes (NADH, 18S, SDHA and HSP90) in Varroa gene expression studies. These four genes remain stable in phoretic and reproductive stage Varroa and are unaffected by Deformed wing virus load. When used for determining changes in vitellogenin gene expression, the signal-to-noise ratio (SNR) for the relatively unstable genes actin and α-tubulin was much lower than for the stable gene combinations (NADH + HSP90 +18S; NADH + HSP90; or NADH). Using both electropherograms and RT-qPCR for short and long amplicons as quality controls, we demonstrate that high quality RNA can be recovered from Varroa up to 10 days later stored at ambient temperature if collected into RNAlater and provided the body is pierced. This protocol allows the exchange of Varroa samples between international collaborators and field sample collectors without requiring frozen collection or shipping. Our results make important contributions to gene expression studies in Varroa by proposing a validated sampling protocol to obtain high quality Varroa RNA and the validation of suitable reference genes for expression studies in this globally important pest.

The ectoparasitic mite Tropilaelaps mercedesae reduces western honey bee, Apismellifera, longevity and emergence weight, and promotes Deformed wing virus infections

Historically an ectoparasite of the native Giant honey bee Apis dorsata, the mite Tropilaelaps mercedesae has switched hosts to the introduced western honey bee Apis mellifera throughout much of Asia. Few data regarding lethal and sub-lethal effects of T. mercedesae on A. mellifera exist, despite its similarity to the devastating mite Varroa destructor. Here we artificially infested worker brood of A. mellifera with T. mercedesae to investigate lethal (longevity) and sub-lethal (emergence weight, Deformed wing virus (DWV) levels and clinical symptoms of DWV) effects of the mite on its new host. The data show that T. mercedesae infestation significantly reduced host longevity and emergence weight, and promoted both DWV levels and associated clinical symptoms. Our results suggest that T. mercedesae is a potentially important parasite to the economically important A. mellifera honey bee.

Impact of the Phoretic Phase on Reproduction and Damage Caused by Varroa destructor (Anderson and Trueman) to Its Host, the European Honey Bee (Apis mellifera L.)

Varroa destructor is a parasitic mite of the honeybee that causes thousands of colony losses worldwide. The parasite cycle is composed of a phoretic and a reproductive phase. During the former, mites stay on adult bees, mostly on nurses, to feed on hemolymph. During the latter, the parasites enter brood cells and reproduce. We investigated if the type of bees on which Varroa stays during the phoretic phase and if the duration of this stay influenced the reproductive success of the parasite and the damage caused to bees. For that purpose, we used an in vitro rearing method developed in our laboratory to assess egg laying rate and the presence and number of fully molted daughters. The expression level of two Varroa vitellogenin genes (VdVg1 and VdVg2), known to vary throughout reproduction, was also quantified. Results showed that the status of the bees or time spent during the phoretic phase impacts neither reproduction parameters nor the Varroa vitellogenin genes levels of expression. However, we correlated these parameters to the gene expression and demonstrated that daughters expressed the vitellogenin genes at lower levels than their mother. Regarding the damage to bees, the data indicated that a longer stay on adult bees during the phoretic phase resulted in more frequent physical deformity in newborn bees. We showed that those mites carry more viral loads of the Deformed Wing Virus and hence trigger more frequently overt infections. This study provides new perspectives towards a better understanding of the Varroa-honeybee interactions.

Hidden Infections and Changing Environments

It is increasingly evident that cryptic stages of many parasites cause asymptomatic infections in a diversity of hosts. This review examines what may cause these infectious agents to persist as asymptomatic infections in invertebrates and how environmental change is linked with the subsequent development of overt infection and disease. In many systems, disease dynamics are closely associated with host condition which, in turn, is linked with environmental change. Symbionts (commensals and mutualists) display similar dynamics when environmental change causes them to exert negative effects on their hosts. Although such asymptomatic infections are demonstrated in a range of invertebrate hosts they are greatly undersampled because most invertebrate diseases are uninvestigated, infections are difficult to detect, and many parasite groups are poorly characterized. A better understanding of the diversity and distribution of parasites that cause asymptomatic infections and of their complex relationships with invertebrate hosts will enable a fuller appreciation of context-dependent host-parasite interactions and will address the biased focus on diseases of invertebrates of practical importance. The existence of such infections could underlie novel disease outbreaks that might otherwise be attributed to invasives while altered disease dynamics may provide an additional and complementary indicator of ecosystem change.

Evidence of Apis cerana Sacbrood virus Infection in Apis mellifera

Sacbrood virus(SBV) is one of the most destructive viruses in the Asian honeybee Apis cerana but is much less destructive in Apis mellifera In previous studies, SBV isolates infecting A. cerana(AcSBV) and SBV isolates infecting A. mellifera(AmSBV) were identified as different serotypes, suggesting a species barrier in SBV infection. In order to investigate this species isolation, we examined the presence of SBV infection in 318A. mellifera colonies and 64A. cerana colonies, and we identified the genotypes of SBV isolates. We also performed artificial infection experiments under both laboratory and field conditions. The results showed that 38A. mellifera colonies and 37A. cerana colonies were positive for SBV infection. Phylogenetic analysis based on RNA-dependent RNA polymerase (RdRp) gene sequences indicated that A. cerana isolates and most A. mellifera isolates formed two distinct clades but two strains isolated fromA. mellifera were clustered with theA. cerana isolates. In the artificial-infection experiments, AcSBV negative-strand RNA could be detected in both adult bees and larvae ofA. mellifera, although there were no obvious signs of the disease, demonstrating the replication of AcSBV inA. mellifera Our results suggest that AcSBV is able to infectA. melliferacolonies with low prevalence (0.63% in this study) and pathogenicity. This work will help explain the different susceptibilities ofA. cerana and A. melliferato sacbrood disease and is potentially useful for guiding beekeeping practices.

Transcriptome analysis of the synganglion from the honey bee mite, Varroa destructor and RNAi knockdown of neural peptide targets

Varroa mites (Varroa destructor) and the viruses that they transmit are one of the major contributing factors to the global honey bee crisis. Gene products within the nervous system are the targets of all the insecticides currently used to control Varroa but there is a paucity of transcriptomic data available for Varroa neural tissues. A cDNA library from the synganglia ("brains") of adult female Varroa was constructed and 600 ESTs sequenced and analysed revealing several current and potential druggable targets. Contigs coding for the deformed wing virus (DWV) variants V. destructor virus-1 (VDV-1) and the recombinant (VDV-1DVD) were present in the synganglion library. Negative-sense RNA-specific PCR indicated that VDV-1 replicates in the Varroa synganglion and all other tissues tested, but we could not detect DWV replicating in any Varroa tissue. Two neuropeptides were identified in the synganlion EST library: a B-type allatostatin and a member of the crustacean hyperglycaemic hormone (CHH) superfamily. Knockdown of the allatostatin or the CHH-like gene by double-stranded RNA-interference (dsRNAi) resulted in 85% and 55% mortality, respectively, of Varroa. Here, we present the first transcriptomic survey in Varroa and demonstrate that neural genes can be targeted by dsRNAi either for genetic validation of putative targets during drug discovery programmes or as a potential control measure in itself.

The Iflaviruses Sacbrood virus and Deformed wing virus evoke different transcriptional responses in the honeybee which may facilitate their horizontal or vertical transmission

Sacbrood virus (SBV) and Deformed wing virus (DWV) are evolutionarily related positive-strand RNA viruses, members of the Iflavirus group. They both infect the honeybee Apis mellifera but have strikingly different levels of virulence when transmitted orally. Honeybee larvae orally infected with SBV usually accumulate high levels of the virus, which halts larval development and causes insect death. In contrast, oral DWV infection at the larval stage usually causes asymptomatic infection with low levels of the virus, although high doses of ingested DWV could lead to DWV replicating to high levels. We investigated effects of DWV and SBV infection on the transcriptome of honeybee larvae and pupae using global RNA-Seq and real-time PCR analysis. This showed that high levels of SBV replication resulted in down-regulation of the genes involved in cuticle and muscle development, together with changes in expression of putative immune-related genes. In particular, honeybee larvae with high levels of SBV replication, with and without high levels of DWV replication, showed concerted up-regulated expression of antimicrobial peptides (AMPs), and down-regulated expression of the prophenoloxidase activating enzyme (PPAE) together with up-regulation of the expression of a putative serpin, which could lead to the suppression of the melanisation pathway. The effects of high SBV levels on expression of these immune genes were unlikely to be a consequence of SBV-induced developmental changes, because similar effects were observed in honeybee pupae infected by injection. In the orally infected larvae with high levels of DWV replication alone we observed no changes of AMPs or of gene expression in the melanisation pathway. In the injected pupae, high levels of DWV alone did not alter expression of the tested melanisation pathway genes, but resulted in up-regulation of the AMPs, which could be attributed to the effect of DWV on the regulation of AMP expression in response to wounding. We propose that the difference in expression of the honeybee immune genes induced by SBV and DWV may be an evolutionary adaptation to the different predominant transmission routes used by these viruses.

Superinfection exclusion and the long-term survival of honey bees in Varroa-infested colonies

Over the past 50 years, many millions of European honey bee (Apis mellifera) colonies have died as the ectoparasitic mite, Varroa destructor, has spread around the world. Subsequent studies have indicated that the mite's association with a group of RNA viral pathogens (Deformed Wing Virus, DWV) correlates with colony death. Here, we propose a phenomenon known as superinfection exclusion that provides an explanation of how certain A. mellifera populations have survived, despite Varroa infestation and high DWV loads. Next-generation sequencing has shown that a non-lethal DWV variant 'type B' has become established in these colonies and that the lethal 'type A' DWV variant fails to persist in the bee population. We propose that this novel stable host-pathogen relationship prevents the accumulation of lethal variants, suggesting that this interaction could be exploited for the development of an effective treatment that minimises colony losses in the future.

Four Categories of Viral Infection Describe the Health Status of Honey Bee Colonies

Honey bee virus prevalence data are an essential prerequisite for managing epidemic events in a population. A survey study was carried out for seven viruses in colonies representing a healthy Danish honey bee population. In addition, colonies from apiaries with high level Varroa infestation or high level of winter mortality were also surveyed. Results from RT-qPCR showed a considerable difference of virus levels between healthy and sick colonies. In the group of healthy colonies, no virus was detected in 36% of cases, while at least one virus was found in each of the sick colonies. Virus titers varied among the samples, and multiple virus infections were common in both groups with a high prevalence of Sacbrood virus (SBV), Black queen cell virus (BQCV) and Deformed wing virus (DWV). Based on the distribution of virus titers, we established four categories of infection: samples free of virus (C = 0), samples with low virus titer (estimated number of virus copies 0 < C < 103), samples with medium virus titer (103 ≤ C < 107) and samples with high virus titer (C ≥ 107). This allowed us to statistically compare virus levels in healthy and sick colonies. Using categories to communicate virus diagnosis results to beekeepers may help them to reach an informed decision on management strategies to prevent further spread of viruses among colonies.

In-depth proteomic analysis of Varroa destructor: Detection of DWV-complex, ABPV, VdMLV and honeybee proteins in the mite

We investigated pathogens in the parasitic honeybee mite Varroa destructor using nanoLC-MS/MS (TripleTOF) and 2D-E-MS/MS proteomics approaches supplemented with affinity-chromatography to concentrate trace target proteins. Peptides were detected from the currently uncharacterized Varroa destructor Macula-like virus (VdMLV), the deformed wing virus (DWV)-complex and the acute bee paralysis virus (ABPV). Peptide alignments revealed detection of complete structural DWV-complex block VP2-VP1-VP3, VDV-1 helicase and single-amino-acid substitution A/K/Q in VP1, the ABPV structural block VP1-VP4-VP2-VP3 including uncleaved VP4/VP2, and VdMLV coat protein. Isoforms of viral structural proteins of highest abundance were localized via 2D-E. The presence of all types of capsid/coat proteins of a particular virus suggested the presence of virions in Varroa. Also, matches between the MWs of viral structural proteins on 2D-E and their theoretical MWs indicated that viruses were not digested. The absence/scarce detection of non-structural proteins compared with high-abundance structural proteins suggest that the viruses did not replicate in the mite; hence, virions accumulate in the Varroa gut via hemolymph feeding. Hemolymph feeding also resulted in the detection of a variety of honeybee proteins. The advantages of MS-based proteomics for pathogen detection, false-positive pathogen detection, virus replication, posttranslational modifications, and the presence of honeybee proteins in Varroa are discussed.

Consequences of Nosema apis infection for male honey bees and their fertility

The queens of eusocial bees, ants and wasps mate only during a very short period early in life and males therefore produce ejaculates consisting of large numbers of high quality sperm. Such extreme selection for high fecundity resulted in males investing minimally into their somatic survival, including their immune system. However, if susceptible males are unable to protect their reproductive tissue from infections, they compromise queen fitness if they transfer pathogens during mating. We used the honey bee Apis mellifera and investigated the course of infection of the sexually transmitted pathogen Nosema apis. We predicted that honey bee males are susceptible but protect their reproductive tissues from infections. We investigated the effects of N. apis infections on the midgut, the accessory glands and the accessory testes and quantified the consequences of infection on male survival and fecundity. We found that N. apis is able to infect males, and as infections progressed, it significantly impacted fertility and survival in older males. Even though we confirm males to be able to minimize N. apis infections of their reproductive tissues, the parasite is present in ejaculates of older males. Consequently N. apis evolved alternative routes to successfully infect ejaculates and get sexually transmitted.

Virus Infections of Honeybees Apis Mellifera

The health and vigour of honeybee colonies are threatened by numerous parasites (such as Varroa destructor and Nosema spp.) and pathogens, including viruses, bacteria, protozoa. Among honeybee pathogens, viruses are one of the major threats to the health and well-being of honeybees and cause serious concern for researchers and beekeepers. To tone down the threats posed by these invasive organisms, a better understanding of bee viral infections will be of crucial importance in developing effective and environmentally benign disease control strategies. Here we summarize recent progress in the understanding of the morphology, genome organization, transmission, epidemiology and pathogenesis of eight honeybee viruses: Deformed wing virus (DWV) and Kakugo virus (KV); Sacbrood virus (SBV); Black Queen cell virus (BQCV); Acute bee paralysis virus (ABPV); Kashmir bee virus (KBV); Israeli Acute Paralysis Virus (IAPV); Chronic bee paralysis virus (CBPV). The review has been designed to provide researchers in the field with updated information about honeybee viruses and to serve as a starting point for future research.

Strand-Specific Quantitative Reverse Transcription-Polymerase Chain Reaction Assay for Measurement of Arenavirus Genomic and Antigenomic RNAs

Arenaviruses are bi-segmented, single-stranded RNA viruses that cause significant human disease. The manner in which they regulate the replication of their genome is not well-understood. This is partly due to the absence of a highly sensitive assay to measure individual species of arenavirus replicative RNAs. To overcome this obstacle, we designed a quantitative reverse transcription (RT)-PCR assay for selective quantitation of each of the lymphocytic choriomeningitis virus (LCMV) genomic or antigenomic RNAs. During the course of assay design, we identified a nonspecific priming phenomenon whereby, in the absence of an RT primer, cDNAs complementary to each of the LCMV replicative RNA species are generated during RT. We successfully circumvented this nonspecific priming event through the use of biotinylated primers in the RT reaction, which permitted affinity purification of primer-specific cDNAs using streptavidin-coated magnetic beads. As proof of principle, we used the assay to map the dynamics of LCMV replication at acute and persistent time points and to determine the quantities of genomic and antigenomic RNAs that are incorporated into LCMV particles. This assay can be adapted to measure total S or L segment-derived viral RNAs and therefore represents a highly sensitive diagnostic platform to screen for LCMV infection in rodent and human tissue samples and can also be used to quantify virus-cell attachment.

Honey bee colony losses and associated viruses

Recent large-scale colony losses among managed Western honey bees (Apis mellifera) have alarmed researchers and apiculturists alike. Here, the existing correlative evidence provided by monitoring studies is reviewed which (i) identified members of the deformed wing virus and acute bee paralysis virus clades as lethal pathogens for entire colonies, and (ii) identified novel viruses whose impact on honey bee health remains elusive. Also discussed in this review is related evidence obtained via controlled experimental infection assays and RNAi approaches underscoring the damage inflicted by some of these viruses on individuals and colonies. The relevance of the ectoparasitic mite Varroa destructor acting as mechanical and biological virus vector for the enhanced virulence of certain viruses or mite selected virus strains is carefully considered.

Prospective large-scale field study generates predictive model identifying major contributors to colony losses

Over the last decade, unusually high losses of colonies have been reported by beekeepers across the USA. Multiple factors such as Varroa destructor, bee viruses, Nosema ceranae, weather, beekeeping practices, nutrition, and pesticides have been shown to contribute to colony losses. Here we describe a large-scale controlled trial, in which different bee pathogens, bee population, and weather conditions across winter were monitored at three locations across the USA. In order to minimize influence of various known contributing factors and their interaction, the hives in the study were not treated with antibiotics or miticides. Additionally, the hives were kept at one location and were not exposed to potential stress factors associated with migration. Our results show that a linear association between load of viruses (DWV or IAPV) in Varroa and bees is present at high Varroa infestation levels (>3 mites per 100 bees). The collection of comprehensive data allowed us to draw a predictive model of colony losses and to show that Varroa destructor, along with bee viruses, mainly DWV replication, contributes to approximately 70% of colony losses. This correlation further supports the claim that insufficient control of the virus-vectoring Varroa mite would result in increased hive loss. The predictive model also indicates that a single factor may not be sufficient to trigger colony losses, whereas a combination of stressors appears to impact hive health.

Lower virus infections in Varroa destructor-infested and uninfested brood and adult honey bees (Apis mellifera) of a low mite population growth colony compared to a high mite population growth colony

A comparison was made of the prevalence and relative quantification of deformed wing virus (DWV), Israeli acute paralysis virus (IAPV), black queen cell virus (BQCV), Kashmir bee virus (KBV), acute bee paralysis virus (ABPV) and sac brood virus (SBV) in brood and adult honey bees (Apis mellifera) from colonies selected for high (HMP) and low (LMP) Varroa destructor mite population growth. Two viruses, ABPV and SBV, were never detected. For adults without mite infestation, DWV, IAPV, BQCV and KBV were detected in the HMP colony; however, only BQCV was detected in the LMP colony but at similar levels as in the HMP colony. With mite infestation, the four viruses were detected in adults of the HMP colony but all at higher amounts than in the LMP colony. For brood without mite infestation, DWV and IAPV were detected in the HMP colony, but no viruses were detected in the LMP colony. With mite infestation of brood, the four viruses were detected in the HMP colony, but only DWV and IAPV were detected and at lower amounts in the LMP colony. An epidemiological explanation for these results is that pre-experiment differences in virus presence and levels existed between the HMP and LMP colonies. It is also possible that low V. destructor population growth in the LMP colony resulted in the bees being less exposed to the mite and thus less likely to have virus infections. LMP and HMP bees may have also differed in susceptibility to virus infection.

Genome sequence heterogeneity of Lake Sinai Virus found in honey bees and Orf1/RdRP-based polymorphisms in a single host

Honey bees (Apis mellifera) are susceptible to a wide range of pathogens, including a broad set of viruses. Recently, next-generation sequencing has expanded the list of viruses with, for instance, two strains of Lake Sinai Virus. Soon after its discovery in the USA, LSV was also discovered in other countries and in other hosts. In the present study, we assemble four almost complete LSV genomes, and show that there is remarkable sequence heterogeneity based on the Orf1, RNA-dependent RNA polymerase and capsid protein sequences in comparison to the previously identified LSV 1 and 2 strains. Phylogenetic analyses of LSV sequences obtained from single honey bee specimens further revealed that up to three distinctive clades could be present in a single bee. Such superinfections have not previously been identified for other honey bee viruses. In a search for the putative routes of LSV transmission, we were able to demonstrate the presence of LSV in pollen pellets and in Varroa destructor mites. However, negative-strand analyses demonstrated that the virus only actively replicates in honey bees and mason bees (Osmia cornuta) and not in Varroa mites.

Differential viral levels and immune gene expression in three stocks of Apis mellifera induced by different numbers of Varroa destructor

The viral levels and immune responses of Italian honey bees (IHB), Russian honey bees (RHB) and an outcross of Varroa Sensitive Hygienic bees (POL) deliberately infested with one or two foundress Varroa were compared. We found that the Deformed wing virus (DWV) level in IHB inoculated with one or two foundress Varroa increased to about 10(3) or 10(5) fold the levels of their uninfested brood. In contrast, POL (10(2) or 10(4) fold) and RHB (10(2) or l0(4) fold) supported a lower increase in DWV levels. The feeding of different stages of Varroa nymphs did not increase DWV levels of their pupal hosts. Analyses of their corresponding Varroa mites showed the same trends: two foundress Varroa yielded higher DWV levels than one foundress, and the addition of nymphs did not increase viral levels. Using the same pupae examined for the presence of viruses, 16 out of 24 genes evaluated showed significant differential mRNA expression levels among the three honey bee stocks. However, only four genes (Defensin, Dscam, PPOact and spaetzle), which were expressed at similar levels in uninfested pupae, were altered by the number of feeding foundress Varroa and levels of DWV regardless of stocks. This research provides the first evidence that immune response profiles of different honey bee stocks are induced by Varroa parasitism.

The Toll-dorsal pathway is required for resistance to viral oral infection in Drosophila

Pathogen entry route can have a strong impact on the result of microbial infections in different hosts, including insects. Drosophila melanogaster has been a successful model system to study the immune response to systemic viral infection. Here we investigate the role of the Toll pathway in resistance to oral viral infection in D. melanogaster. We show that several Toll pathway components, including Spätzle, Toll, Pelle and the NF-kB-like transcription factor Dorsal, are required to resist oral infection with Drosophila C virus. Furthermore, in the fat body Dorsal is translocated from the cytoplasm to the nucleus and a Toll pathway target gene reporter is upregulated in response to Drosophila C Virus infection. This pathway also mediates resistance to several other RNA viruses (Cricket paralysis virus, Flock House virus, and Nora virus). Compared with control, viral titres are highly increased in Toll pathway mutants. The role of the Toll pathway in resistance to viruses in D. melanogaster is restricted to oral infection since we do not observe a phenotype associated with systemic infection. We also show that Wolbachia and other Drosophila-associated microbiota do not interact with the Toll pathway-mediated resistance to oral infection. We therefore identify the Toll pathway as a new general inducible pathway that mediates strong resistance to viruses with a route-specific role. These results contribute to a better understanding of viral oral infection resistance in insects, which is particularly relevant in the context of transmission of arboviruses by insect vectors.

Pathogenic microbes can enter their hosts through different routes. This can have a strong impact on which host defensive mechanisms are elicited and in disease outcome. We used the model organism Drosophila melanogaster to understand how resistance to viruses differs between infection by direct virus entry into the body cavity and infection through feeding on food with the virus. We show that the Toll pathway is required to resist oral infection with different RNA viruses. On the other hand this pathway does not influence the outcome of viral infection performed by injection. Together our results show that the Toll pathway has a route-specific general antiviral effect. Our work expands the role of this classical innate immunity pathway and contributes to a better understanding of viral oral infection resistance in insects. This is particularly relevant because insect vectors of emerging human viral diseases, like dengue, are infected through feeding on contaminated hosts.

Infectivity of DWV associated to flower pollen: experimental evidence of a horizontal transmission route

Deformed wing virus (DWV) is a honeybee pathogen whose presence is generally associated with infestation of the colony by the mite Varroa destructor, leading to the onset of infections responsible for the collapse of the bee colony. DWV contaminates bee products such as royal jelly, bee-bread and honey stored within the infected hive. Outside the hive, DWV has been found in pollen loads collected directly from infected as well as uninfected forager bees. It has been shown that the introduction of virus-contaminated pollen into a DWV-free hive results in the production of virus-contaminated food, whose role in the development of infected bees from virus-free eggs has been experimentally demonstrated. The aim of this study was twofold: (i) to ascertain the presence of DWV on pollen collected directly from flowers visited by honeybees and then quantify the viral load and (ii) determine whether the virus associated with pollen is infective. The results of our investigation provide evidence that DWV is present on pollen sampled directly from visited flowers and that, following injection in individuals belonging to the pollinator species Apis mellifera, it is able to establish an active infection, as indicated by the presence of replicating virus in the head of the injected bees. We also provide the first indication that the pollinator species Osmia cornuta is susceptible to DWV infection.

On the front line: quantitative virus dynamics in honeybee (Apis mellifera L.) colonies along a new expansion front of the parasite Varroa destructor

Over the past fifty years, annual honeybee (Apis mellifera) colony losses have been steadily increasing worldwide. These losses have occurred in parallel with the global spread of the honeybee parasite Varroa destructor. Indeed, Varroa mite infestations are considered to be a key explanatory factor for the widespread increase in annual honeybee colony mortality. The host-parasite relationship between honeybees and Varroa is complicated by the mite's close association with a range of honeybee viral pathogens. The 10-year history of the expanding front of Varroa infestation in New Zealand offered a rare opportunity to assess the dynamic quantitative and qualitative changes in honeybee viral landscapes in response to the arrival, spread and level of Varroa infestation. We studied the impact of de novo infestation of bee colonies by Varroa on the prevalence and titres of seven well-characterised honeybee viruses in both bees and mites, using a large-scale molecular ecology approach. We also examined the effect of the number of years since Varroa arrival on honeybee and mite viral titres. The dynamic shifts in the viral titres of black queen cell virus and Kashmir bee virus mirrored the patterns of change in Varroa infestation rates along the Varroa expansion front. The deformed wing virus (DWV) titres in bees continued to increase with Varroa infestation history, despite dropping infestation rates, which could be linked to increasing DWV titres in the mites. This suggests that the DWV titres in mites, perhaps boosted by virus replication, may be a major factor in maintaining the DWV epidemic after initial establishment. Both positive and negative associations were identified for several pairs of viruses, in response to the arrival of Varroa. These findings provide important new insights into the role of the parasitic mite Varroa destructor in influencing the viral landscape that affects honeybee colonies.

Honeybees currently face a dramatic decline worldwide. The main honeybee parasite - Varroa destructor - plays a key role in these mortalities, since uncontrolled infestation inevitably results in the death of the colony. The pathological effects of Varroa infestations are partly attributed to the association of the mite with several honeybee viruses, primarily deformed wing virus (DWV). However the exact role that Varroa plays in the spread of honeybee viruses is still unknown. The recent arrival of Varroa in New Zealand provided a timely opportunity to gain insights into the complex relationship between bees, Varroa and viruses. Our data reveal that the different viruses have unique quantitative dynamics in relation to Varroa infestation, resulting in a shifting succession of virus infections that ultimately leaves DWV as the predominant infection. Assumption-free analysis shows consistent clustering of the data according to Varroa-infestation history, confirming a progressive change in the overall virus landscape co-incident with Varroa infestation. We also highlight possible interactions between several viruses. Our findings may have implications for the beekeeping industry, by highlighting the dynamic changes in the virus infections due to the arrival of Varroa, and how these are maintained.

Israeli acute paralysis virus: epidemiology, pathogenesis and implications for honey bee health

Israeli acute paralysis virus (IAPV) is a widespread RNA virus of honey bees that has been linked with colony losses. Here we describe the transmission, prevalence, and genetic traits of this virus, along with host transcriptional responses to infections. Further, we present RNAi-based strategies for limiting an important mechanism used by IAPV to subvert host defenses. Our study shows that IAPV is established as a persistent infection in honey bee populations, likely enabled by both horizontal and vertical transmission pathways. The phenotypic differences in pathology among different strains of IAPV found globally may be due to high levels of standing genetic variation. Microarray profiles of host responses to IAPV infection revealed that mitochondrial function is the most significantly affected biological process, suggesting that viral infection causes significant disturbance in energy-related host processes. The expression of genes involved in immune pathways in adult bees indicates that IAPV infection triggers active immune responses. The evidence that silencing an IAPV-encoded putative suppressor of RNAi reduces IAPV replication suggests a functional assignment for a particular genomic region of IAPV and closely related viruses from the Family Dicistroviridae, and indicates a novel therapeutic strategy for limiting multiple honey bee viruses simultaneously and reducing colony losses due to viral diseases. We believe that the knowledge and insights gained from this study will provide a new platform for continuing studies of the IAPV–host interactions and have positive implications for disease management that will lead to mitigation of escalating honey bee colony losses worldwide.

The mysterious outbreak of honey bee Colony Collapse Disorder (CCD) in the US in 2006–2007 has attracted massive media attention and created great concerns over the effects of various risk factors on bee health. Understanding the factors that are linked to the honey bee colony declines may provide insights for managing similar incidents in the future. We conducted this study to elucidate traits of a key honey bee virus, Israeli acute paralysis virus. We then developed an innovative strategy to control virus levels. The knowledge and insights gained from this study will have positive implications for bee disease management, helping to mitigate worldwide colony losses.

A virulent strain of deformed wing virus (DWV) of honeybees (Apis mellifera) prevails after Varroa destructor-mediated, or in vitro, transmission

The globally distributed ectoparasite Varroa destructor is a vector for viral pathogens of the Western honeybee (Apis mellifera), in particular the Iflavirus Deformed Wing Virus (DWV). In the absence of Varroa low levels DWV occur, generally causing asymptomatic infections. Conversely, Varroa-infested colonies show markedly elevated virus levels, increased overwintering colony losses, with impairment of pupal development and symptomatic workers. To determine whether changes in the virus population were due Varroa amplifying and introducing virulent virus strains and/or suppressing the host immune responses, we exposed Varroa-naïve larvae to oral and Varroa-transmitted DWV. We monitored virus levels and diversity in developing pupae and associated Varroa, the resulting RNAi response and transcriptome changes in the host. Exposed pupae were stratified by Varroa association (presence/absence) and virus levels (low/high) into three groups. Varroa-free pupae all exhibited low levels of a highly diverse DWV population, with those exposed per os (group NV) exhibiting changes in the population composition. Varroa-associated pupae exhibited either low levels of a diverse DWV population (group VL) or high levels of a near-clonal virulent variant of DWV (group VH). These groups and unexposed controls (C) could be also discriminated by principal component analysis of the transcriptome changes observed, which included several genes involved in development and the immune response. All Varroa tested contained a diverse replicating DWV population implying the virulent variant present in group VH, and predominating in RNA-seq analysis of temporally and geographically separate Varroa-infested colonies, was selected upon transmission from Varroa, a conclusion supported by direct injection of pupae in vitro with mixed virus populations. Identification of a virulent variant of DWV, the role of Varroa in its transmission and the resulting host transcriptome changes furthers our understanding of this important viral pathogen of honeybees.

Honeybees are the most important managed pollinating insect, contributing billions of dollars to annual global agricultural production. Over the last century a parasitic mite, Varroa, has spread worldwide, with significant impacts on honeybee colony health as a consequence of its transmission of a cocktail of viruses while feeding on honeybee ‘blood’. The most important virus for colony health is deformed wing virus (DWV), high levels of which cause developmental deformities and premature ageing resulting in high overwintering colony losses. In experiments on individual Varroa-exposed pupae we demonstrate that a single type of virulent DWV is amplified 1,000–10,000 times in the recipient pupae, despite the mite containing a high diversity of replicating DWV strains. We could recapitulate this by direct injection of pupae with mixed virus populations, showing the virulent strain is advantaged by the route of transmission. In parallel, we detected changes in the immune response and developmental gene expression of the honeybee and propose that these contribute to the characteristic pathogenesis of DWV. Identification of a virulent strain of DWV has implications for therapeutic or prophylactic interventions to improve honeybee colony health, as well as contributing to our understanding of the biology of this important honeybee viral pathogen.

Honey bee Apis mellifera parasites in the absence of Nosema ceranae fungi and Varroa destructor mites

Few areas of the world have western honey bee (Apis mellifera) colonies that are free of invasive parasites Nosema ceranae (fungi) and Varroa destructor (mites). Particularly detrimental is V. destructor; in addition to feeding on host haemolymph, these mites are important vectors of several viruses that are further implicated as contributors to honey bee mortality around the world. Thus, the biogeography and attendant consequences of viral communities in the absence of V. destructor are of significant interest. The island of Newfoundland, Province of Newfoundland and Labrador, Canada, is free of V. destructor; the absence of N. ceranae has not been confirmed. Of 55 Newfoundland colonies inspected visually for their strength and six signs of disease, only K-wing had prevalence above 5% (40/55 colonies = 72.7%). Similar to an earlier study, screenings again confirmed the absence of V. destructor, small hive beetles Aethina tumida (Murray), tracheal mites Acarapis woodi (Rennie), and Tropilaelaps spp. ectoparasitic mites. Of a subset of 23 colonies screened molecularly for viruses, none had Israeli acute paralysis virus, Kashmir bee virus, or sacbrood virus. Sixteen of 23 colonies (70.0%) were positive for black queen cell virus, and 21 (91.3%) had some evidence for deformed wing virus. No N. ceranae was detected in molecular screens of 55 colonies, although it is possible extremely low intensity infections exist; the more familiar N. apis was found in 53 colonies (96.4%). Under these conditions, K-wing was associated (positively) with colony strength; however, viruses and N. apis were not. Furthermore, black queen cell virus was positively and negatively associated with K-wing and deformed wing virus, respectively. Newfoundland honey bee colonies are thus free of several invasive parasites that plague operations in other parts of the world, and they provide a unique research arena to study independent pathology of the parasites that are present.

Increased tolerance and resistance to virus infections: a possible factor in the survival of Varroa destructor-resistant honey bees (Apis mellifera)

The honey bee ectoparasitic mite, Varroa destructor, has a world-wide distribution and inflicts more damage than all other known apicultural diseases. However, Varroa-induced colony mortality is more accurately a result of secondary virus infections vectored by the mite. This means that honey bee resistance to Varroa may include resistance or tolerance to virus infections. The aim of this study was to see if this is the case for a unique population of mite-resistant (MR) European honey bees on the island of Gotland, Sweden. This population has survived uncontrolled mite infestation for over a decade, developing specific mite-related resistance traits to do so. Using RT-qPCR techniques, we monitored late season virus infections, Varroa mite infestation and honey bee colony population dynamics in the Gotland MR population and compared this to mite-susceptible (MS) colonies in a close by apiary. From summer to autumn the deformed wing virus (DWV) titres increased similarly between the MR and MS populations, while the black queen cell virus (BQCV) and sacbrood virus (SBV) titres decreased substantially in the MR population compared to the MS population by several orders of magnitude. The MR colonies all survived the following winter with high mite infestation, high DWV infection, small colony size and low proportions of autumn brood, while the MS colonies all perished. Possible explanations for these changes in virus titres and their relevance to Varroa resistance and colony winter survival are discussed.

The ability to cause infection in a pathogenic fungus uncovers a new biological feature of honey bee viruses

We demonstrated that honey bee viruses including Deformed wing virus (DWV), Black queen cell virus (BQCV) and Israeli acute paralysis virus (IAPV) could infect and replicate in the fungal pathogen Ascosphaera apis that causes honey bee chalkbrood disease, revealing a novel biological feature of honey bee viruses. The phylogenetic analysis show that viruses of fungal and honey bee origins form two clusters in the phylogenetic trees distinctly and that host range of honey bee viruses is dynamic. Further studies are warranted to investigate the impact of the viruses on the fitness of their fungal host and phenotypic effects the virus-fungus combination has on honey bee hosts.

Dynamics of the presence of israeli acute paralysis virus in honey bee colonies with colony collapse disorder

The determinants of Colony Collapse Disorder (CCD), a particular case of collapse of honey bee colonies, are still unresolved. Viruses including the Israeli acute paralysis virus (IAPV) were associated with CCD. We found an apiary with colonies showing typical CCD characteristics that bore high loads of IAPV, recovered some colonies from collapse and tested the hypothesis if IAPV was actively replicating in them and infectious to healthy bees. We found that IAPV was the dominant pathogen and it replicated actively in the colonies: viral titers decreased from April to September and increased from September to December. IAPV extracted from infected bees was highly infectious to healthy pupae: they showed several-fold amplification of the viral genome and synthesis of the virion protein VP3. The health of recovered colonies was seriously compromised. Interestingly, a rise of IAPV genomic copies in two colonies coincided with their subsequent collapse. Our results do not imply IAPV as the cause of CCD but indicate that once acquired and induced to replication it acts as an infectious factor that affects the health of the colonies and may determine their survival. This is the first follow up outside the US of CCD-colonies bearing IAPV under natural conditions.

Chronic bee paralysis virus in honeybee queens: evaluating susceptibility and infection routes

Chronic bee paralysis virus (CBPV) is known as a disease of worker honey bees. To investigate pathogenesis of the CBPV on the queen, the sole reproductive individual in a colony, we conducted experiments regarding the susceptibility of queens to CBPV. Results from susceptibility experiment showed a similar disease progress in the queens compared to worker bees after infection. Infected queens exhibit symptoms by Day 6 post infection and virus levels reach 10¹¹ copies per head. In a transmission experiment we showed that social interactions may affect the disease progression. Queens with forced contact to symptomatic worker bees acquired an overt infection with up to 10¹¹ virus copies per head in six days. In contrast, queens in contact with symptomatic worker bees, but with a chance to receive food from healthy bees outside the cage appeared healthy. The virus loads did not exceed 10⁷ in the majority of these queens after nine days. Symptomatic worker bees may transmit sufficient active CBPV particles to the queen through trophallaxis, to cause an overt infection.