Pesticide use negatively affects bumble bees across European landscapes.

Sustainable agriculture requires balancing crop yields with the effects of pesticides on non-target organisms, such as bees and other crop pollinators. Field studies demonstrated that agricultural use of neonicotinoid insecticides can negatively affect wild bee species1,2, leading to restrictions on these compounds3. However, besides neonicotinoids, field-based evidence of the effects of landscape pesticide exposure on wild bees is lacking. Bees encounter many pesticides in agricultural landscapes4-9 and the effects of this landscape exposure on colony growth and development of any bee species remains unknown. Here we show that the many pesticides found in bumble bee-collected pollen are associated with reduced colony performance during crop bloom, especially in simplified landscapes with intensive agricultural practices. Our results from 316 Bombus terrestris colonies at 106 agricultural sites across eight European countries confirm that the regulatory system fails to sufficiently prevent pesticide-related impacts on non-target organisms, even for a eusocial pollinator species in which colony size may buffer against such impacts10,11. These findings support the need for postapproval monitoring of both pesticide exposure and effects to confirm that the regulatory process is sufficiently protective in limiting the collateral environmental damage of agricultural pesticide use.

Varroa destructor shapes the unique viral landscapes of the honey bee populations of the Azores archipelago.

The worldwide dispersal of the ectoparasitic mite Varroa destructor from its Asian origins has fundamentally transformed the relationship of the honey bee (Apis mellifera) with several of its viruses, via changes in transmission and/or host immunosuppression. The extent to which honey bee-virus relationships change after Varroa invasion is poorly understood for most viruses, in part because there are few places in the world with several geographically close but completely isolated honey bee populations that either have, or have not, been exposed long-term to Varroa, allowing for separate ecological, epidemiological, and adaptive relationships to develop between honey bees and their viruses, in relation to the mite's presence or absence. The Azores is one such place, as it contains islands with and without the mite. Here, we combined qPCR with meta-amplicon deep sequencing to uncover the relationship between Varroa presence, and the prevalence, load, diversity, and phylogeographic structure of eight honey bee viruses screened across the archipelago. Four viruses were not detected on any island (ABPV-Acute bee paralysis virus, KBV-Kashmir bee virus, IAPV-Israeli acute bee paralysis virus, BeeMLV-Bee macula-like virus); one (SBV-Sacbrood virus) was detected only on mite-infested islands; one (CBPV-Chronic bee paralysis virus) occurred on some islands, and two (BQCV-Black queen cell virus, LSV-Lake Sinai virus,) were present on every single island. This multi-virus screening builds upon a parallel survey of Deformed wing virus (DWV) strains that uncovered a remarkably heterogeneous viral landscape featuring Varroa-infested islands dominated by DWV-A and -B, Varroa-free islands naïve to DWV, and a refuge of the rare DWV-C dominating the easternmost Varroa-free islands. While all four detected viruses investigated here were affected by Varroa for one or two parameters (usually prevalence and/or the Richness component of ASV diversity), the strongest effect was observed for the multi-strain LSV. Varroa unambiguously led to elevated prevalence, load, and diversity (Richness and Shannon Index) of LSV, with these results largely shaped by LSV-2, a major LSV strain. Unprecedented insights into the mite-virus relationship were further gained from implementing a phylogeographic approach. In addition to enabling the identification of a novel LSV strain that dominated the unique viral landscape of the easternmost islands, this approach, in combination with the recovered diversity patterns, strongly suggests that Varroa is driving the evolutionary change of LSV in the Azores. This study greatly advances the current understanding of the effect of Varroa on the epidemiology and adaptive evolution of these less-studied viruses, whose relationship with Varroa has thus far been poorly defined.

Cold case: The disappearance of Egypt bee virus, a fourth distinct master strain of deformed wing virus linked to honeybee mortality in 1970's Egypt.

In 1977, a sample of diseased adult honeybees (Apis mellifera) from Egypt was found to contain large amounts of a previously unknown virus, Egypt bee virus, which was subsequently shown to be serologically related to deformed wing virus (DWV). By sequencing the original isolate, we demonstrate that Egypt bee virus is in fact a fourth unique, major variant of DWV (DWV-D): more closely related to DWV-C than to either DWV-A or DWV-B. DWV-A and DWV-B are the most common DWV variants worldwide due to their close relationship and transmission by Varroa destructor. However, we could not find any trace of DWV-D in several hundred RNA sequencing libraries from a worldwide selection of honeybee, varroa and bumblebee samples. This means that DWV-D has either become extinct, been replaced by other DWV variants better adapted to varroa-mediated transmission, or persists only in a narrow geographic or host range, isolated from common bee and beekeeping trade routes.

Managed bumble bees acquire parasites from their foraging environment: A case study on parasite spillback.

The use of commercially reared bumble bees in agricultural environments has been recognized as a potential threat to wild pollinators due to competition, genetic contamination, and most notably, disease transmission. Higher parasite prevalence near greenhouses where managed bumble bees are used has been linked to parasite spillover from managed to wild bees. However, pathogen transmission is not unidirectional, and can also flow from wild to managed bees. These newly infected managed bees can subsequently re-infect (other) wild bees, in a process known as spillback, which is an alternative explanation for the increased parasite prevalence near greenhouses. Reducing parasite prevalence in managed bees is key to controlling host-parasite dynamics in cases of spillover; in spillback, producing managed bees that are resilient to infection is important. Here we establish that the managed bumble bee Bombus terrestris can acquire parasites from their foraging environment, which is the major infection route for Apicystis spp. and Crithidia spp., but not for Nosema spp.. Managed B. terrestris were found to have a higher prevalence of Crithdia and a higher load of Apicystis than local wild conspecifics, showing that for these parasites, spillback is a possible risk scenario.

Feeding Honeybee Colonies with Honeybee-Specific Lactic Acid Bacteria (Hbs-LAB) Does Not Affect Colony-Level Hbs-LAB Composition or Paenibacillus larvae Spore Levels, Although American Foulbrood Affected Colonies Harbor a More Diverse Hbs-LAB Community.

The main current methods for controlling American Foulbrood (AFB) in honeybees, caused by the bacterial pathogen Paenibacillus larvae, are enforced incineration or prophylactic antibiotic treatment, neither of which is fully satisfactory. This has led to an increased interest in the natural relationships between the pathogenic and mutualistic microorganisms of the honeybee microbiome, in particular, the antagonistic effects of Honeybee-Specific Lactic Acid Bacteria (hbs-LAB) against P. larvae. We investigated whether supplemental administration of these bacteria affected P. larvae infection at colony level over an entire flowering season. Over the season, the supplements affected neither colony-level hbs-LAB composition nor naturally subclinical or clinical P. larvae spore levels. The composition of hbs-LAB in colonies was, however, more diverse in apiaries with a history of clinical AFB, although this was also unrelated to P. larvae spore levels. During the experiments, we also showed that qPCR could detect a wider range of hbs-LAB, with higher specificity and sensitivity than mass spectrometry. Honeybee colonies are complex super-organisms where social immune defenses, natural homeostatic mechanisms, and microbiome diversity and function play a major role in disease resistance. This means that observations made at the individual bee level cannot be simply extrapolated to infer similar effects at colony level. Although individual laboratory larval assays have clearly demonstrated the antagonistic effects of hbs-LAB on P. larvae infection, the results from the experiments presented here indicate that direct conversion of such practice to colony-level administration of live hbs-LAB is not effective.

Correction to: American foulbrood in a honeybee colony: spore-symptom relationship and feedbacks between disease and colony development.

Unfortunately, the original version of the article [1] contained an error. The author has brought to our attention that the article title is truncated in the published version. The correct title is American foulbrood in a honeybee colony: spore-symptom relationship and feedbacks between disease and colony development. Instead, it was published inadvertently as American foulbrood in a honeybee colony: spore symptom relationship and feedbacks due to an error occurred during the production process.

American foulbrood in a honeybee colony: spore-symptom relationship and feedbacks.

BACKGROUND: The most severe bacterial disease of honeybees is American foulbrood (AFB). The epidemiology of AFB is driven by the extreme spore resilience, the difficulty of bees to remove these spores, and the considerable incidence of undetected spore-producing colonies. The honeybee collective defence mechanisms and their feedback on colony development, which involves a division of labour at multiple levels of colony organization, are difficult to model. To better predict disease outbreaks we need to understand the feedback between colony development and disease progression within the colony. We therefore developed Bayesian models with data from forty AFB-diseased colonies monitored over an entire foraging season to (i) investigate the relationship between spore production and symptoms, (ii) disentangle the feedback loops between AFB epidemiology and natural colony development, and (iii) discuss whether larger insect societies promote or limit within-colony disease transmission. RESULTS: Rather than identifying a fixed spore count threshold for clinical symptoms, we estimated the probabilities around the relationship between spore counts and symptoms, taking into account modulators such as brood amount/number of bees and time post infection. We identified a decrease over time in the bees-to-brood ratio related to disease development, which should ultimately induce colony collapse. Lastly, two contrasting theories predict that larger colonies could promote either higher (classical epidemiological SIR-model) or lower (increasing spatial nest segregation and more effective pathogen removal) disease prevalence. CONCLUSIONS: AFB followed the predictions of the SIR-model, partly because disease prevalence and brood removal are decoupled, with worker bees acting more as disease vectors, infecting new brood, than as agents of social immunity, by removing infected brood. We therefore established a direct link between disease prevalence and social group size for a eusocial insect. We furthermore provide a probabilistic description of the relationship between AFB spore counts and symptoms, and how disease development and colony strength over a season modulate this relationship. These results help to better understand disease development within honeybee colonies, provide important estimates for further epidemiological modelling, and gained important insights into the optimal sampling strategy for practical beekeeping and honeybee research.

ICTV Virus Taxonomy Profile: Solinviviridae.

Solinviviridae is a family of picorna/calici-like viruses with non-segmented, linear, positive-sense RNA genomes of approximately 10-11 kb. Unusually, their capsid proteins are encoded towards the 3'-end of the genome where they can be expressed both from a subgenomic RNA and as an extension of the replication (picorna-like helicase-protease-polymerase) polyprotein. Members of two species within the family infect ants, but related unclassified virus sequences derive from a large variety of insects and other arthropods. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the Solinviviridae, which is available at www.ictv.global/report/solinviviridae.

ICTV Virus Taxonomy Profile: Polycipiviridae.

Polycipiviridae is a family of picorna-like viruses with non-segmented, linear, positive-sense RNA genomes of approximately 10-12 kb. Unusually for viruses within the order Picornavirales, their genomes are polycistronic, with four (or more) consecutive 5'-proximal open reading frames (ORFs) encoding structural (and possibly other) proteins and a long 3' ORF encoding the replication polyprotein. Members of species within the family have all been detected in ants or via arthropod transcriptomic datasets. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the Polycipiviridae, which is available at www.ictv.global/report/polycipiviridae.

Honeybee-Specific Lactic Acid Bacterium Supplements Have No Effect on American Foulbrood-Infected Honeybee Colonies.

Paenibacillus larvae, the causative agent of American foulbrood (AFB), is the primary bacterial pathogen affecting honeybees and beekeeping. The main methods for controlling AFB are incineration of diseased colonies or prophylactic antibiotic treatment (e.g., with tylosin), neither of which is fully satisfactory. The search for superior means for controlling AFB has led to an increased interest in the natural relationships between the honeybee-pathogenic and mutualistic microorganisms and, in particular, the antagonistic effects of honeybee-specific lactic acid bacteria (hbs-LAB) against P. larvae These effects have been demonstrated only on individual larvae in controlled laboratory bioassays. Here we investigated whether supplemental administration of hbs-LAB had a similar beneficial effect on P. larvae infection at colony level. We compared experimentally AFB-infected colonies treated with hbs-LAB supplements to untreated and tylosin-treated colonies and recorded AFB symptoms, bacterial spore levels, and two measures of colony health. To account for the complexity of a bee colony, we focused on (Bayesian) probabilities and magnitudes of effect sizes. Tylosin reduced AFB disease symptoms but also had a negative effect on colony strength. The tylosin treatment did not, however, affect P. larvae spore levels and might therefore "mask" the potential for disease. hbs-LAB tended to reduce brood size in the short term but was unlikely to affect AFB symptoms or spores. These results do not contradict demonstrated antagonistic effects of hbs-LAB against P. larvae at the individual bee level but rather suggest that supplementary administration of hbs-LAB may not be the most effective way to harness these beneficial effects at the colony level.IMPORTANCE The previously demonstrated antagonistic effects of honeybee-derived bacterial microbiota on the infectivity and pathogenicity of P. larvae in laboratory bioassays have identified a possible new approach to AFB control. However, honeybee colonies are complex superorganisms where social immune defenses play a major role in resistance against disease at the colony level. Few studies have investigated the effect of beneficial microorganisms on bee diseases at the colony level. Effects observed at the individual bee level do not necessarily translate into similar effects at the colony level. This study partially fills this gap by showing that, unlike at the individual level, hbs-LAB supplements did not affect AFB symptoms at the colony level. The inference is that the mechanisms regulating the honeybee microbial dynamics within a colony are too strong to manipulate positively through supplemental feeding of live hbs-LAB and that new potential remedies identified through laboratory research have to be tested thoroughly in situ, in colonies.

Viral infections alter antennal epithelium ultrastructure in honey bees.

Varroa destructor and its associated viruses, in particular deformed wing virus (DWV), have been identified as probable causes of honey bee (Apis mellif era L.) colony losses. Evidence suggests that elevated DWV titres in bees could compromise sensory and communication abilities resulting in negative consequences for hygienic behaviour. As antennae play a central role in this behaviour, we compared antennal ultrastructure in DWV-symptomatic and asymptomatic bees. The results show that virus capsids accumulate in the basal regions of the antennal epithelium, close to the haemolymph. No virus particles were detected at the level of sensory sensilla, such as pore plates, nor within the sensory cell dendrites associated with these sensilla. However, membranous structures appeared to be more prevalent in supporting cells surrounding the dendrites of DWV-symptomatic bees. Para-crystalline arrays containing large numbers of virus particles were detected in the antennae of DWV-symptomatic bees but not in asymptomatic bees.

Virion Structure of Black Queen Cell Virus, a Common Honeybee Pathogen.

Viral diseases are a major threat to honeybee (Apis mellifera) populations worldwide and therefore an important factor in reliable crop pollination and food security. Black queen cell virus (BQCV) is the etiological agent of a fatal disease of honeybee queen larvae and pupae. The virus belongs to the genus Triatovirus from the family Dicistroviridae, which is part of the order Picornavirales Here we present a crystal structure of BQCV determined to a resolution of 3.4 Å. The virion is formed by 60 copies of each of the major capsid proteins VP1, VP2, and VP3; however, there is no density corresponding to a 75-residue-long minor capsid protein VP4 encoded by the BQCV genome. We show that the VP4 subunits are present in the crystallized virions that are infectious. This aspect of the BQCV virion is similar to that of the previously characterized triatoma virus and supports the recent establishment of the separate genus Triatovirus within the family Dicistroviridae The C terminus of VP1 and CD loops of capsid proteins VP1 and VP3 of BQCV form 34-Å-tall finger-like protrusions at the virion surface. The protrusions are larger than those of related dicistroviruses.IMPORTANCE The western honeybee is the most important pollinator of all, and it is required to sustain the agricultural production and biodiversity of wild flowering plants. However, honeybee populations worldwide are suffering from virus infections that cause colony losses. One of the most common, and least known, honeybee pathogens is black queen cell virus (BQCV), which at high titers causes queen larvae and pupae to turn black and die. Here we present the three-dimensional virion structure of BQCV, determined by X-ray crystallography. The structure of BQCV reveals large protrusions on the virion surface. Capsid protein VP1 of BQCV does not contain a hydrophobic pocket. Therefore, the BQCV virion structure provides evidence that capsid-binding antiviral compounds that can prevent the replication of vertebrate picornaviruses may be ineffective against honeybee virus infections.

Virion Structure of Iflavirus Slow Bee Paralysis Virus at 2.6-Angstrom Resolution.

UNLABELLED: The western honeybee (Apis mellifera) is the most important commercial insect pollinator. However, bees are under pressure from habitat loss, environmental stress, and pathogens, including viruses that can cause lethal epidemics. Slow bee paralysis virus (SBPV) belongs to the Iflaviridae family of nonenveloped single-stranded RNA viruses. Here we present the structure of the SBPV virion determined from two crystal forms to resolutions of 3.4 Å and 2.6 Å. The overall structure of the virion resembles that of picornaviruses, with the three major capsid proteins VP1 to 3 organized into a pseudo-T3 icosahedral capsid. However, the SBPV capsid protein VP3 contains a C-terminal globular domain that has not been observed in other viruses from the order Picornavirales The protruding (P) domains form "crowns" on the virion surface around each 5-fold axis in one of the crystal forms. However, the P domains are shifted 36 Å toward the 3-fold axis in the other crystal form. Furthermore, the P domain contains the Ser-His-Asp triad within a surface patch of eight conserved residues that constitutes a putative catalytic or receptor-binding site. The movements of the domain might be required for efficient substrate cleavage or receptor binding during virus cell entry. In addition, capsid protein VP2 contains an RGD sequence that is exposed on the virion surface, indicating that integrins might be cellular receptors of SBPV. IMPORTANCE: Pollination by honeybees is needed to sustain agricultural productivity as well as the biodiversity of wild flora. However, honeybee populations in Europe and North America have been declining since the 1950s. Honeybee viruses from the Iflaviridae family are among the major causes of honeybee colony mortality. We determined the virion structure of an Iflavirus, slow bee paralysis virus (SBPV). SBPV exhibits unique structural features not observed in other picorna-like viruses. The SBPV capsid protein VP3 has a large C-terminal domain, five of which form highly prominent protruding "crowns" on the virion surface. However, the domains can change their positions depending on the conditions of the environment. The domain includes a putative catalytic or receptor binding site that might be important for SBPV cell entry.

Virion Structure of Israeli Acute Bee Paralysis Virus.

UNLABELLED: The pollination services provided by the western honeybee (Apis mellifera) are critical for agricultural production and the diversity of wild flowering plants. However, honeybees suffer from environmental pollution, habitat loss, and pathogens, including viruses that can cause fatal diseases. Israeli acute bee paralysis virus (IAPV), from the family Dicistroviridae, has been shown to cause colony collapse disorder in the United States. Here, we present the IAPV virion structure determined to a resolution of 4.0 Å and the structure of a pentamer of capsid protein protomers at a resolution of 2.7 Å. IAPV has major capsid proteins VP1 and VP3 with noncanonical jellyroll ß-barrel folds composed of only seven instead of eight ß-strands, as is the rule for proteins of other viruses with the same fold. The maturation of dicistroviruses is connected to the cleavage of precursor capsid protein VP0 into subunits VP3 and VP4. We show that a putative catalytic site formed by the residues Asp-Asp-Phe of VP1 is optimally positioned to perform the cleavage. Furthermore, unlike many picornaviruses, IAPV does not contain a hydrophobic pocket in capsid protein VP1 that could be targeted by capsid-binding antiviral compounds. IMPORTANCE: Honeybee pollination is required for agricultural production and to sustain the biodiversity of wild flora. However, honeybee populations in Europe and North America are under pressure from pathogens, including viruses that cause colony losses. Viruses from the family Dicistroviridae can cause honeybee infections that are lethal, not only to individual honeybees, but to whole colonies. Here, we present the virion structure of an Aparavirus, Israeli acute bee paralysis virus (IAPV), a member of a complex of closely related viruses that are distributed worldwide. IAPV exhibits unique structural features not observed in other picorna-like viruses. Capsid protein VP1 of IAPV does not contain a hydrophobic pocket, implying that capsid-binding antiviral compounds that can prevent the replication of vertebrate picornaviruses may be ineffective against honeybee virus infections.

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.

Bees under stress: sublethal doses of a neonicotinoid pesticide and pathogens interact to elevate honey bee mortality across the life cycle.

Microbial pathogens are thought to have a profound impact on insect populations. Honey bees are suffering from elevated colony losses in the northern hemisphere possibly because of a variety of emergent microbial pathogens, with which pesticides may interact to exacerbate their impacts. To reveal such potential interactions, we administered at sublethal and field realistic doses one neonicotinoid pesticide (thiacloprid) and two common microbial pathogens, the invasive microsporidian Nosema ceranae and black queen cell virus (BQCV), individually to larval and adult honey bees in the laboratory. Through fully crossed experiments in which treatments were administered singly or in combination, we found an additive interaction between BQCV and thiacloprid on host larval survival likely because the pesticide significantly elevated viral loads. In adult bees, two synergistic interactions increased individual mortality: between N. ceranae and BQCV, and between N. ceranae and thiacloprid. The combination of two pathogens had a more profound effect on elevating adult mortality than N. ceranae plus thiacloprid. Common microbial pathogens appear to be major threats to honey bees, while sublethal doses of pesticide may enhance their deleterious effects on honey bee larvae and adults. It remains an open question as to whether these interactions can affect colony survival.

Analysis of reference gene stability after Israeli acute paralysis virus infection in bumblebees Bombus terrestris.

To date, there are no validated internal reference genes for the normalization of RT-qPCR data from virus infection experiments with pollinating insects. In this study we evaluated the stability of five candidate internal reference genes: elongation factor-1-alpha (ELF1α), peptidylprolyl isomerase A (PPIA), 60S ribosomal protein L23 (RPL23), TATA-binding protein (TBP) and polyubiquitin (UBI), in relation to Israeli acute paralysis virus (IAPV) infection of Bombus terrestris. We investigated the stability of these genes: in whole bodies and individual body parts, as well as in whole bodies collected at different time intervals after infection with IAPV. Our data identified PPIA as the single, most-optimal internal reference gene and the combination of PPAI-RPL23-UBI as a fully-sufficient multiple internal reference genes set for IAPV infection experiments in B. terrestris.

Effect of oral infection with Kashmir bee virus and Israeli acute paralysis virus on bumblebee (Bombus terrestris) reproductive success.

Israeli acute paralysis virus (IAPV) together with Acute bee paralysis virus (ABPV) and Kashmir bee virus (KBV) constitute a complex of closely related dicistroviruses. They are infamous for their high mortality after injection in honeybees. These viruses have also been reported in non-Apis hymenopteran pollinators such as bumblebees, which got infected with IAPV when placed in the same greenhouse with IAPV infected honeybee hives. Here we orally infected Bombus terrestris workers with different doses of either IAPV or KBV viral particles. The success of the infection was established by analysis of the bumblebees after the impact studies: 50days after infection. Doses of 0.5×10(7) and 1×10(7) virus particles per bee were infectious over this period, for IAPV and KBV respectively, while a dose of 0.5×10(6) IAPV particles per bee was not infectious. The impact of virus infection was studied in micro-colonies consisting of 5 bumblebees, one of which becomes a pseudo-queen which proceeds to lay unfertilized (drone) eggs. The impact parameters studied were: the establishment of a laying pseudo-queen, the timing of egg-laying, the number of drones produced, the weight of these drones and worker mortality. In this setup KBV infection resulted in a significant slower colony startup and offspring production, while only the latter can be reported for IAPV. Neither virus increased worker mortality, at the oral doses used. We recommend further studies on how these viruses transmit between different pollinator species. It is also vital to understand how viral prevalence can affect wild bee populations because disturbance of the natural host-virus association may deteriorate the already critically endangered status of many bumblebee species.

Origins, diversity, and adaptive evolution of DWV in the honey bees of the Azores: the impact of the invasive mite Varroa destructor.

Deformed wing virus (DWV) is a honey bee virus, whose emergence from relative obscurity is driven by the recent host-switch, adaptation, and global dispersal of the ectoparasitic mite Varroa destructor (a highly efficient vector of DWV) to reproduction on honey bees (Apis mellifera). Our study examines how varroa affects the continuing evolution of DWV, using the Azores archipelago, where varroa is present on only three out of the eight Islands, as a natural experimental system for comparing different evolutionary conditions and trajectories. We combined qPCR of 494 honey bee colonies sampled across the archipelago with amplicon deep sequencing to reveal how the DWV genetic landscape is altered by varroa. Two of the varroa-free Islands were also free of DWV, while a further two Islands were intriguingly dominated by the rare DWV-C major variant. The other four Islands, including the three varroa-infested Islands, were dominated by the common DWV-A and DWV-B variants. The varroa-infested Islands had, as expected, an elevated DWV prevalence relative to the uninfested Islands, but not elevated DWV loads, due the relatively high prevalence and loads of DWV-C on the varroa-free Islands. This establishes the Azores as a stable refuge for DWV-C and provides the most convincing evidence to date that at least some major strains of DWV may be capable of not just surviving, but actually thriving in honey bees in the absence of varroa-mediated transmission. We did not detect any change in DWV genetic diversity associated with island varroa status but did find a positive association of DWV diversity with virus load, irrespective of island varroa status.