Epidemiology of a major honey bee pathogen, deformed wing virus: potential worldwide replacement of genotype A by genotype B

Dominance of recombinant DWV genomes with changing viral landscapes as revealed in national US honey bee and varroa mite survey

Honey bees are essential pollinators for global agriculture. The viromes of US commercial apiaries and their ectoparasitic mites are poorly characterized at a strain level and there is a need to integrate genomics into pathogen surveillance. We sequenced RNA viromes from 383 adult bees and 173 mites pooled samples from 11 major US beekeeping hubs in 2021, assembling 45 complete and 1702 partial genomes. Protein sequence similarity networks and recombinant genome identification revealed a new viral landscape. Sinaivirus (n = 312), Iflavirus sacbroodi (n = 280), and Iflavirus aladeformis (DWV, n = 135) genomes were common. Recombinant DWV genomes with high nucleotide identity were widespread, and DWV type A master variants were rare, with an indication that RT-PCR surveillance may over-detect type A due to the prevalence of recombinant DWV genomes. Future work should use genomic strategies to avoid misidentification of common honey bee virus genomes and their impact on colony health.

Honey bee viruses in the yellow-legged hornet Vespa velutina (Lepelieter 1836): Prevalence, loads, and detection of replicative DWV and LSV forms

Apiaries in Galicia, northwestern Spain, are currently facing the invasive alien species Vespa velutina, which is well established in the region. The pressure on honey bee colonies is high, resulting in both economic and ecological losses. Honey bee colonies also face the challenge of viruses, which are becoming increasingly diverse. In recent years, honey bee viruses have been spreading across taxonomic groups beyond Apoidea, infecting the Vespoidea superfamily. This cross-species spillover has raised concerns in the scientific community due to the potential risk of viruses spreading in ecosystems. Currently, there is a lack of knowledge on this topic, and further research is needed to address this issue. This study employed qPCR and sequencing to investigate the prevalence, loads, and presence of replicative forms of important honey bee viruses in V. velutina individuals collected from 11 apiaries in Galicia. All V. velutina individuals tested positive for DWV, BQCV, AKI complex (ABPV, KBV, and IAPV), or LSV but not for CBPV. DWV showed the highest prevalence (97.0 %) and loads, with both DWV-A (67.4 %) and DWV-B (32.6 %) being detected. The AKI complex (46.3 %) and LSV (43.3 %) were also common, whereas BQCV (11.9 %) was rarer. LSV is detected for the first time in V. velutina. LSV-2 was the dominant strain (82.1 %), and two less frequent (17.9 %) unknown strains were also detected. All 44 screened V. velutina samples carried the replicative form of DWV, and six of these also carried the replicative form of LSV, raising for the first time the possibility of co-infection in the hornet. The detection of honey bee viruses in V. velutina, and the ability of these viruses to spread to other species, may indicate a potential risk of spillover in the apiaries.

Molecular Characterization and Phylogenetic Analysis of Honeybee (Apis mellifera) Mite-Borne Pathogen DWV-A and DWV-B Isolated from Lithuania

Deformed wing virus (DWV) is known as one of the main viruses that affect honeybees' health all around the world. The virus has two widespread genotypes, DWV-A and DWV-B (VDV-1), transmitted mainly by V. destructor mites. In this study, we collected honeycombs with covered broods from 73 apiaries in eight Lithuanian regions and initially investigated the prevalence of V. destructor mites. Mites were collected from May to the end of July in 2021 from 124 hives. The prevalence of V. destructor infestations in beehives reached 30% and 63% in investigated apiaries. The presence of DWV-A and DWV-B pathogens in mites and broods was examined by RT-qPCR targeting the CRPV-capsid region. The molecular characterization of the virus in mite samples was based on sequence analysis of the RNA-dependent RNA polymerase (RdRp) region. In addition, leader polypeptide (LP), structural protein (Vp3), Helicase, and RdRp genes were used for phylogenetic characterization of dual infection. The prevalences of DWV-B in mites and broods were 56.5% and 31.5%, respectively, while DWV-A was detected in 12.9% of mite samples and 24.7% of brood samples. Some of the examined mite samples harboured dual virus infections. Our findings showed that bee colonies from the same apiary were not always infected by the same viruses. Some bee colonies were virus-free, while others were highly infected. Phylogenetic analysis of 21 sequences demonstrated the presence of highly variable DWV-B and DWV-A genotypes in Lithuania and possible recombinant variants of the virus. This study represents the first molecular characterization of mite-borne pathogens hosted by honeybees (Apis mellifera) in Lithuania.

Common viral infections inhibit egg laying in honey bee queens and are linked to premature supersedure

With their long lives and extreme reproductive output, social insect queens have escaped the classic trade-off between fecundity and lifespan, but evidence for a trade-off between fecundity and immunity has been inconclusive. This is in part because pathogenic effects are seldom decoupled from effects of immune induction. We conducted parallel, blind virus infection experiments in the laboratory and in the field to interrogate the idea of a reproductive immunity trade-off in honey bee (Apis mellifera) queens and to better understand how these ubiquitous stressors affect honey bee queen health. We found that queens injected with infectious virus had smaller ovaries and were less likely to recommence egg-laying than controls, while queens injected with UV-inactivated virus displayed an intermediate phenotype. In the field, heavily infected queens had smaller ovaries and infection was a meaningful predictor of whether supersedure cells were observed in the colony. Immune responses in queens receiving live virus were similar to queens receiving inactivated virus, and several of the same immune proteins were negatively associated with ovary mass in the field. This work supports the hypothesized relationship between virus infection and symptoms associated with queen failure and suggests that a reproductive-immunity trade-off is partially, but not wholly responsible for these effects.

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.

Varroa destructor and deformed wing virus interaction increases incidence of winter mortality in honey bee colonies

Winter mortality of honey bee colonies represents a major source of economic loss for the beekeeping industry. The objectives of this prospective study were to estimate the incidence risk of winter colony mortality in southwestern Quebec, Canada and to evaluate and quantify the impact of the associated risk factors. A total of 242 colonies from 31 apiaries was selected for sampling in August 2017. The presence of Varroa destructor, Vairimorpha (Nosema) spp., Melissococcus plutonius, deformed wing virus (DWV), and viruses of the acute-Kashmir-Israeli complex (AKI complex) was investigated in each colony. Management practices of the various colonies were obtained from a questionnaire. The incidence risk of colony mortality during the winter of 2017-2018 was estimated to be 26.5% [95% confidence interval (CI): 15.4 to 40.3]. In logistic regression modeling of winter mortality in colonies, an interaction was discovered between V. destructor and DWV; the detection of ≥ 1 V. destructor mites per 100 bees was associated with higher odds of mortality (3.46, 95% CI: 1.35 to 8.90) compared to colonies with < 1 mite per 100 bees, but only in DWV-positive colonies. There were more colony losses in apiaries from beekeepers owning 1 to 5 colonies than in apiaries from beekeepers owning over 100 colonies, which suggests that beekeeper experience and/or type of management are important contributors to winter colony mortality. Assuming a causal relationship, the results of this study suggest that up to 9% of all colony mortalities in the population could have been prevented by reducing the level of V. destructor to < 1 mite per 100 bees in all colonies.

Characterization of a Molecular Clone of Deformed Wing Virus B

Honey bees (Apis mellifera) play a crucial role in agriculture through their pollination activities. However, they have faced significant health challenges over the past decades that can limit colony performance and even lead to collapse. A primary culprit is the parasitic mite Varroa destructor, known for transmitting harmful bee viruses. Among these viruses is deformed wing virus (DWV), which impacts bee pupae during their development, resulting in either pupal demise or in the emergence of crippled adult bees. In this study, we focused on DWV master variant B. DWV-B prevalence has risen sharply in recent decades and appears to be outcompeting variant A of DWV. We generated a molecular clone of a typical DWV-B strain to compare it with our established DWV-A clone, examining RNA replication, protein expression, and virulence. Initially, we analyzed the genome using RACE-PCR and RT-PCR techniques. Subsequently, we conducted full-genome RT-PCR and inserted the complete viral cDNA into a bacterial plasmid backbone. Phylogenetic comparisons with available full-length sequences were performed, followed by functional analyses using a live bee pupae model. Upon the transfection of in vitro-transcribed RNA, bee pupae exhibited symptoms of DWV infection, with detectable viral protein expression and stable RNA replication observed in subsequent virus passages. The DWV-B clone displayed a lower virulence compared to the DWV-A clone after the transfection of synthetic RNA, as evidenced by a reduced pupal mortality rate of only 20% compared to 80% in the case of DWV-A and a lack of malformations in 50% of the emerging bees. Comparable results were observed in experiments with low infection doses of the passaged virus clones. In these tests, 90% of bees infected with DWV-B showed no clinical symptoms, while 100% of pupae infected with DWV-A died. However, at high infection doses, both DWV-A and DWV-B caused mortality rates exceeding 90%. Taken together, we have generated an authentic virus clone of DWV-B and characterized it in animal experiments.

Covert RNA viruses in medflies differ in their mode of transmission and tissue tropism

Numerous studies have demonstrated the presence of covert viral infections in insects. These infections can be transmitted in insect populations via two main routes: vertical from parents to offspring, or horizontal between nonrelated individuals. Thirteen covert RNA viruses have been described in the Mediterranean fruit fly (medfly). Some of these viruses are established in different laboratory-reared and wild medfly populations, although variations in the viral repertoire and viral levels have been observed at different time points. To better understand these viral dynamics, we characterized the prevalence and levels of covert RNA viruses in two medfly strains, assessed the route of transmission of these viruses, and explored their distribution in medfly adult tissues. Altogether, our results indicated that the different RNA viruses found in medflies vary in their preferred route of transmission. Two iflaviruses and a narnavirus are predominantly transmitted through vertical transmission via the female, while a nodavirus and a nora virus exhibited a preference for horizontal transmission. Overall, our results give valuable insights into the viral tropism and transmission of RNA viruses in the medfly, contributing to the understanding of viral dynamics in insect populations.

IMPORTANCE

The presence of RNA viruses in insects has been extensively covered. However, the study of host-virus interaction has focused on viruses that cause detrimental effects to the host. In this manuscript, we uncovered which tissues are infected with covert RNA viruses in the agricultural pest Ceratitis capitata, and which is the preferred transmission route of these viruses. Our results showed that vertical and horizontal transmission can occur simultaneously, although each virus is transmitted more efficiently following one of these routes. Additionally, our results indicated an association between the tropism of the RNA virus and the preferred route of transmission. Overall, these results set the basis for understanding how viruses are established and maintained in medfly populations.

Selection and Multiplexing of Reverse Transcription-Quantitative PCR Tests Targeting Relevant Honeybee Viral Pathogens

Verifying the inclusivity of molecular detection methods gives indications about the reliability of viral infection diagnosis because of the tendency of viral pathogens to undergo sequence variation. This study was aimed at selecting inclusive probes based on reverse transcription-quantitative PCR (RT-qPCR) assays for the diagnosis of the most widespread and detrimental viruses infecting honeybees, namely the acute bee paralysis virus (ABPV), the black queen cell virus (BQCV), the chronic paralysis bee virus (CBPV), the deformed wing virus variants A (DWVA) and B (DWVB), and the sacbrood virus (SBV). Therefore, previously described detection methods were re-evaluated in silico for their specificity and inclusivity. Based on this evaluation, selected methods were modified, or new ones were designed and tested in duplex RT-qPCR reactions. The limits of detection (LODs), effect of multiplexing on sensitivity and the viral RNA quantification potential in bees and hive debris were assessed. This study made available diagnostic assays able to detect an increased number of virus variants compared with previously described tests and two viral pathogens in a single PCR reaction.

Caught in the act: the invasion of a viral vector changes viral prevalence and titre in native honeybees and bumblebees

Novel transmission routes change pathogen landscapes and may facilitate disease emergence. The varroa mite is a virus vector that switched to western honeybees at the beginning of the last century, leading to hive mortality, particularly in combination with RNA viruses. A recent invasion of varroa on the French island of Ushant introduced vector-mediated transmission to one of the last varroa-naive native honeybee populations and caused rapid changes in the honeybee viral community. These changes were characterized by a drastic increase in deformed wing virus type B prevalence and titre in honeybees, as well as knock-on effects in bumblebees, particularly in the year following the invasion. Slow bee paralysis virus also appeared in honeybees and bumblebees, with a 1 year delay, while black queen cell virus declined in honeybees. This study highlights the rapid and far-reaching effects of vector-borne transmission that can extend beyond the directly affected host species, and that the direction of the effect depends on the pathogen's virulence.

Varroa destructor parasitism and Deformed wing virus infection in honey bees are linked to peroxisome-induced pathways

The ectoparasitic mite Varroa destructor transmits and triggers viral infections that have deleterious effects on honey bee colonies worldwide. We performed a manipulative experiment in which worker bees collected at emergence were exposed to Varroa for 72 h, and their proteomes were compared with those of untreated control bees. Label-free quantitative proteomics identified 77 differentially expressed A. mellifera proteins (DEPs). In addition, viral proteins were identified by orthogonal analysis, and most importantly, Deformed wing virus (DWV) was found at high levels/intensity in Varroa-exposed bees. Pathway enrichment analysis suggested that the main pathways affected included peroxisomal metabolism, cyto-/exoskeleton reorganization, and cuticular proteins. Detailed examination of individual DEPs revealed that additional changes in DEPs were associated with peroxisomal function. In addition, the proteome data support the importance of TGF-β signaling in Varroa-DWV interaction and the involvement of the mTORC1 and Hippo pathways. These results suggest that the effect of DWV on bees associated with Varroa feeding results in aberrant autophagy. In particular, autophagy is selectively modulated by peroxisomes, to which the observed proteome changes strongly corresponded. This study complements previous research with different study designs and suggests the importance of the peroxisome, which plays a key role in viral infections.

Eggs sampling as an effective tool for identifying the incidence of viruses in honey bees involved in artificial queen rearing

The Carniolan honey bee (Apis mellifera carnica) plays an essential role in crop pollination, environment diversity, and the production of honey bee products. However, the health of individual honey bees and their colonies is under pressure due to multiple stressors, including viruses as a significant threat to bees. Monitoring various virus infections could be a crucial selection tool during queen rearing. In the present study, samples from all developmental stages (eggs, larvae, pupae, and queens) were screened for the incidence of seven viruses during queen rearing in Slovenia. The screening of a total of 108 samples from five queen breeders was performed by the RT-qPCR assays. The results showed that the highest incidence was observed for black queen cell virus (BQCV), Lake Sinai virus 3 (LSV3), deformed wing virus B (DWV-B), and sacbrood virus (SBV). The highest viral load was detected in queens (6.07 log10 copies/queen) and larvae (5.50 log10 copies/larva) for BQCV, followed by SBV in larvae (5.47 log10 copies/larva). When comparing all the honey bee developmental stages, the eggs exhibited general screening for virus incidence and load in queen mother colonies. The results suggest that analyzing eggs is a good indicator of resilience to virus infection during queen development.

Host Barriers Limit Viral Spread in a Spillover Host: A Study of Deformed Wing Virus in the Bumblebee Bombus terrestris

The transmission of pathogens from reservoir to recipient host species, termed pathogen spillover, can profoundly impact plant, animal, and public health. However, why some pathogens lead to disease emergence in a novel species while others fail to establish or do not elicit disease is often poorly understood. There is strong evidence that deformed wing virus (DWV), an (+)ssRNA virus, spills over from its reservoir host, the honeybee Apis mellifera, into the bumblebee Bombus terrestris. However, the low impact of DWV on B. terrestris in laboratory experiments suggests host barriers to virus spread in this recipient host. To investigate potential host barriers, we followed the spread of DWV genotype B (DWV-B) through a host's body using RT-PCR after experimental transmission to bumblebees in comparison to honeybees. Inoculation was per os, mimicking food-borne transmission, or by injection into the bee's haemocoel, mimicking vector-based transmission. In honeybees, DWV-B was present in both honeybee faeces and haemolymph within 3 days of inoculation per os or by injection. In contrast, DWV-B was not detected in B. terrestris haemolymph after inoculation per os, suggesting a gut barrier that hinders DWV-B's spread through the body of a B. terrestris. DWV-B was, however, detected in B. terrestris faeces after injection and feeding, albeit at a lower abundance than that observed for A. mellifera, suggesting that B. terrestris sheds less DWV-B than A. mellifera in faeces when infected. Barriers to viral spread in B. terrestris following oral infection may limit DWV's impact on this spillover host and reduce its contribution to the community epidemiology of DWV.

Prevalence of pathogens in honey bee colonies and association with clinical signs in southwestern Quebec, Canada

Honey bees can be affected by a variety of pathogens, which impacts their vital role as pollinators in agriculture. A cross-sectional study was conducted in southwestern Quebec to: i) estimate the prevalence of 11 bee pathogens; ii) assess the agreement between beekeeper suspicion of a disease and laboratory detection of the causative pathogen; and iii) explore the association between observed clinical signs and pathogen detection in a colony. A total of 242 colonies in 31 apiaries owned by 15 beekeepers was sampled in August 2017. The prevalence of Varroa destructor detection was estimated as 48% for colonies and 93% for apiaries. The apparent prevalence of colonies infected by Nosema spp. and Melissococcus plutonius was estimated as 40% and 21%, respectively. At least 180 colonies were tested by polymerase chain reaction (PCR) for deformed wing virus (DWV), acute-Kashmir-Israeli complex (AKI complex), and black queen cell virus (BQCV), which were detected in 33%, 9%, and 95% of colonies, respectively. Acarapis woodi, Paenibacillus larvae, and Aethina tumida were not detected. Varroasis was suspected by beekeepers in 14 of the 15 beekeeping operations in which the mite was detected. However, no correlation was found between suspected European foulbrood and detection of M. plutonius or between suspected nosemosis and detection of Nosema spp. Colony weakness was associated with Nosema spore counts of at least 0.5 × 106 per bee. Melissococcus plutonius was more frequently detected in colonies showing scattered brood.

Seasonal Variability in the Prevalence of DWV Strains in Individual Colonies of European Honeybees in Hawaii

The most prevalent viral pathogen of honeybees is Deformed Wing Virus (DWV) and its two most widely studied and common master-variants are DWV-A and DWV-B. The prevalence of DWV variants in the UK and in the US is changing, with the prevalence of the DWV-A strain declining and DWV-B increasing over time. In 2012, only DWV-A was detected on the Hawaiian Islands of Oahu. In this study we focused on a colony-level survey of DWV strains in a single apiary and examined the prevalence of DWV variants over the course of two years. In 2018 and 2019, a total of 16 colonies underwent viral testing in January, May, and September. Of those 16 colonies, four were monitored in both 2018 and 2019. Individual colonies showed variability of DWV master variants throughout the sampling period. DWV-A was consistently detected; however, the detection of DWV-B was variable across time in individual colonies. Ultimately, this study demonstrated a seasonal variation in both viral prevalence and load for DWV-B, providing a perspective on the dynamic nature of DWV master variants emerging in Hawaii.

A semi-automated and high-throughput approach for the detection of honey bee viruses in bee samples

Deformed wing virus (DWV) was first detected in dead honey bees in 1982 but has been in honey bees for at least 300 years. Due to its high prevalence and virulence, they have been linked with the ongoing decline in honey bee populations worldwide. A rapid, simple, semi-automated, high-throughput, and cost-effective method of screening colonies for viruses would benefit bee research and the beekeeping industry. Here we describe a semi-automated approach that combines an RNA-grade liquid homogenizer followed by magnetic bead capture for total virus nucleic acid extraction. We compare it to the more commonly applied nucleic acid column-based purification method and use qPCR plus Oxford Nanopore Technologies sequencing to evaluate the accuracy of analytical results for both methods. Our results showed high reproducibility and accuracy for both approaches. The semi-automated method described here allows for faster screening of viral loads in units of 96 samples at a time. We developed this method to monitor viral loads in honey bee colonies, but it could be easily applied for any PCR or genomic-based screening assays.

Integrated Pest Management Strategies to Control Varroa Mites and Their Effect on Viral Loads in Honey Bee Colonies

Honey bee viruses in combination with varroa mite are very damaging for honey bee colonies worldwide. There are no effective methods to control the viral load in honey bee colonies except regular and effective control of mites. Integrated Pest Management strategies are required to effectively control mites with veterinary medicines based on organic compounds. We evaluated the effect of two brood interruption techniques, queen caging (QC) and trapping comb (TC), followed by an oxalic acid treatment, on the mite fall, colony strength, and viral load of Deformed Wing Virus (DWV) and Acute Bee Paralysis Virus (ABPV). In this paper, we report the data obtained in two experimental sites, in Slovenia and Italy, in terms of the varroacide efficacy, colony strength, and viral load. The number of adult bees after the adoption of the two techniques showed similar decreasing trends in both locations. The viral load of Acute Bee Paralysis Virus did not show any significant reduction after 25 days, reported as the number of Real-Time PCR cycles needed to detect the virus. The viral load of DWV also did not show a significant reduction after 25 days. The acaricidal efficacy of the applied protocols was high in both experimental groups and in both apiaries. Both the queen caging and trapping comb techniques, followed by an oxalic acid treatment, can be considered effective varroa treatment strategies, but further studies should be carried out to evaluate the long-term effects on viral loads to plan the Integrated Pest Management strategy with the right timing before wintering.

Shift in virus composition in honeybees (Apis mellifera) following worldwide invasion by the parasitic mite and virus vector Varroa destructor

Invasive vectors can induce dramatic changes in disease epidemiology. While viral emergence following geographical range expansion of a vector is well known, the influence a vector can have at the level of the host's pathobiome is less well understood. Taking advantage of the formerly heterogeneous spatial distribution of the ectoparasitic mite Varroa destructor that acts as potent virus vector among honeybees Apis mellifera, we investigated the impact of its recent global spread on the viral community of honeybees in a retrospective study of historical samples. We hypothesized that the vector has had an effect on the epidemiology of several bee viruses, potentially altering their transmissibility and/or virulence, and consequently their prevalence, abundance, or both. To test this, we quantified the prevalence and loads of 14 viruses from honeybee samples collected in mite-free and mite-infested populations in four independent geographical regions. The presence of the mite dramatically increased the prevalence and load of deformed wing virus, a cause of unsustainably high colony losses. In addition, several other viruses became more prevalent or were found at higher load in mite-infested areas, including viruses not known to be actively varroa-transmitted, but which may increase opportunistically in varroa-parasitized bees.

Virus replication in the honey bee parasite, Varroa destructor

IMPORTANCE

The parasitic mite Varroa destructor is a significant driver of worldwide colony losses of our most important commercial pollinator, the Western honey bee Apis mellifera. Declines in honey bee health are frequently attributed to the viruses that mites vector to honey bees, yet whether mites passively transmit viruses as a mechanical vector or actively participate in viral amplification and facilitate replication of honey bee viruses is debated. Our work investigating the antiviral RNA interference response in V. destructor demonstrates that key viruses associated with honey bee declines actively replicate in mites, indicating that they are biological vectors, and the host range of bee-associated viruses extends to their parasites, which could impact virus evolution, pathogenicity, and spread.

Signatures of adaptive decreased virulence of deformed wing virus in an isolated population of wild honeybees (Apis mellifera)

Understanding the ecological and evolutionary processes that drive host-pathogen interactions is critical for combating epidemics and conserving species. The Varroa destructor mite and deformed wing virus (DWV) are two synergistic threats to Western honeybee (Apis mellifera) populations across the globe. Distinct honeybee populations have been found to self-sustain despite Varroa infestations, including colonies within the Arnot Forest outside Ithaca, NY, USA. We hypothesized that in these bee populations, DWV has been selected to produce an avirulent infection phenotype, allowing for the persistence of both host and disease-causing agents. To investigate this, we assessed the titre of viruses in bees from the Arnot Forest and managed apiaries, and assessed genomic variation and virulence differences between DWV isolates. Across groups, we found viral abundance was similar, but DWV genotypes were distinct. We also found that infections with isolates from the Arnot Forest resulted in higher survival and lower rates of symptomatic deformed wings, compared to analogous isolates from managed colonies, providing preliminary evidence to support the hypothesis of adaptive decreased viral virulence. Overall, this multi-level investigation of virus genotype and phenotype indicates that host ecological context can be a significant driver of viral evolution and host-pathogen interactions in honeybees.

The presence of identical deformed wing virus sequence variants in co-occurring Apis species in Northern Thailand may represent a potential epidemiological threat to native honey bees of Southeast Asia

Widespread native honey bee species in South and East Asia (Apis cerana, Apis dorsata and Apis florea) and the imported western honey bee (Apis mellifera) share habitats and potentially also share pathogens. Chief among the threats facing A. mellifera in Europe and North America is deformed wing virus (DWV), including its two principal genotypes: A and B (DWV-A and DWV-B respectively). Though DWV-A has been recorded in Asia's native Apis species, it is not known if DWV-B, or both DWV-A and DWV-B, are currently widespread in Asia and, if so, whether viral transmission is primarily intraspecific or interspecific. This study aims to fill these knowledge gaps by (i) determining the DWV genotype in four co-occurring Apis host species using qPCR and (ii) inferring viral transmission between them using nucleotide sequences of DWV from Apis host species collected at three independent localities in Northern Thailand. We found DWV-A and -B in all four Apis species, the exotic A. mellifera and the native A. cerana, A. dorsata and A. florea. That DWV-A sequences were identical across Apis species at the same locality, with a similar pattern for DWV-B sequences, suggests that DWV's epidemiology is largely driven by ongoing interspecific transmission (spillover) of DWV across co-occurring native and exotic Apis species. Both genotypes of DWV represent a serious threat to Asia's exotic and native honey bee species.

Age-related response to mite parasitization and viral infection in the honey bee suggests a trade-off between growth and immunity

Host age at parasites' exposure is often neglected in studies on host-parasite interactions despite the important implications for epidemiology. Here we compared the impact of the parasitic mite Varroa destructor, and the associated pathogenic virus DWV on different life stages of their host, the western honey bee Apis mellifera. The pre-imaginal stages of the honey bee proved to be more susceptible to mite parasitization and viral infection than adults. The higher viral load in mite-infested bees and DWV genotype do not appear to be the drivers of the observed difference which, instead, seems to be related to the immune-competence of the host. These results support the existence of a trade-off between immunity and growth, making the pupa, which is involved in the highly energy-demanding process of metamorphosis, more susceptible to parasites and pathogens. This may have important implications for the evolution of the parasite's virulence and in turn for honey bee health. Our results highlight the important role of host's age and life stage at exposure in epidemiological modelling. Furthermore, our study could unravel new aspects of the complex honey bee-Varroa relationship to be addressed for a sustainable management of this parasite.

Doomed drones? Using passage experiments and mathematical modelling to determine Deformed wing virus population dynamics in male honeybees

Varroa destructor is an ectoparasitic mite of honeybees which vectors a range of pathogenic viruses, the most notable being Deformed wing virus (DWV). Mites parasitise bees during pupal development and male honeybees, drones, have a longer development cycle than female workers (24 versus 21 days), allow for more progeny mites to develop per foundress (1.6-2.5 compared to 0.7-1.45). How this longer exposure time influences evolution of the transmitted virus population is unknown. Using uniquely tagged viruses recovered from cDNA we investigated the replication, competition and morbidity of DWV genotypes in drones. Assays examining virus replication and morbidity revealed drones are highly susceptible to both predominant genotypes of DWV. In virus passage studies using an equimolar inocula of major DNA genotypes and their recombinants, the recombinant form dominated but did not reach 100% of the virus population within 10 passages. Using an in-silico model of the virus-mite-bee system we examined bottlenecks during virus acquisition by the mite and subsequent injection of viruses into the host, which may play a significant role in shaping virus diversity. This study furthers our understanding of the variables influencing DWV diversity changes and provides insight into areas of future research in the mite-virus-bee system.

Viral Co-Infections and Antiviral Immunity in Honey Bees

Over the past few decades, honey bees have been facing an increasing number of stressors. Beyond individual stress factors, the synergies between them have been identified as a key factor in the observed increase in colony mortality. However, these interactions are numerous and complex and call for further research. Here, in line with our need for a systemic understanding of the threats that they pose to bee health, we review the interactions between honey bee viruses. As viruses are obligate parasites, the interactions between them not only depend on the viruses themselves but also on the immune responses of honey bees. Thus, we first summarise our current knowledge of the antiviral immunity of honey bees. We then review the interactions between specific pathogenic viruses and their interactions with their host. Finally, we draw hypotheses from the current literature and suggest directions for future research.

All together now: Geographically coordinated miticide treatment benefits honey bee health

Deformed wing virus (DWV) is a pathogenic virus of honey bees transmitted by the ectoparasitic mite Varroa destructor. Annual overwintering colony losses, accounting for ~25% of all colonies, are associated with high levels of Varroa-DWV infestation. Effective miticide treatments are available to control Varroa. However, the absence of coordinated treatment means environmental transmission of mites continues unchecked. We aimed to determine whether rational, coordinated treatment is beneficial, and characterized the DWV population as an indicator of colony health.This study uses coordinated treatment of Varroa in a geographically isolated environment (Isle of Arran, Scotland) over 3 years. The study area contained 50-84 colonies managed by ~20 amateur beekeepers. Sampling and virus analysis to assess strain diversity and viral loads were conducted before and after treatments, and changes in population diversity were quantified by sequence analysis.Over the 3 years analysis of the virus population revealed that the dominant DWV variant shifted from Type A to Type B in all apiaries, regardless of mite levels or proximity to other colonies. During this period the number of managed colonies increased by 47% (57-84 colonies), but despite this, we estimate total mite numbers decreased by 58%. Synthesis and applications. In this study, the beekeepers in Arran significantly improved the number of colonies they managed, without importing any bees onto the island, indicating that an improved focus on management techniques, through the combination of a coordinated miticide programme and an improved understanding of bee diseases, could yield positive results for bee health and sustainability.

Conservation measures or hotspots of disease transmission? Agri-environment schemes can reduce disease prevalence in pollinator communities

Insects are under pressure from agricultural intensification. To protect pollinators, conservation measures such as the EU agri-environment schemes (AES) promote planting wildflowers along fields. However, this can potentially alter disease ecology by serving as transmission hubs or by diluting infections. We tested this by measuring plant-pollinator interactions and virus infections (DWV-A, DWV-B and ABPV) across pollinator communities in agricultural landscapes over a year. AES had a direct effect on DWV-B, reducing prevalence and load in honeybees, with a tentative general dilution effect on load in early summer. DWV-A prevalence was reduced both under AES and with increasing niche overlap between competent hosts, likely via a dilution effect. By contrast, AES had no impact on ABPV, its prevalence driven by the proportion of bumblebees in the community. Epidemiological differences were also reflected in the virus phylogenies, with DWV-B showing recent rapid expansion, while DWV-A and ABPV showed slower growth rates and geographical population structure. Phylogenies indicate that all three viruses freely circulate across their host populations. Our study illustrates how complex interactions between environmental, ecological and evolutionary factors may influence wildlife disease dynamics. Supporting pollinator nutrition can mitigate the transmission of important bee diseases, providing an unexpected boost to pollinator conservation. This article is part of the theme issue 'Infectious disease ecology and evolution in a changing world'.

Nationwide Screening for Bee Viruses in Apis mellifera Colonies in Egypt

Honey bees are essential for crop and wild plant pollination. However, many countries have reported high annual colony losses caused by multiple possible stressors. Diseases, particularly those caused by viruses, are a major cause of colony losses. However, little is known about the prevalence of honey bee pathogens, particularly virus prevalence, in Egyptian honey bees. To address this shortfall, we determined the prevalence of widespread bee viruses in honey bee colonies in Egypt-whether it is affected by geography, the season, or infestation with Varroa destructor (varroa) mites. Honey bee worker samples were collected from 18 geographical regions across Egypt during two seasons: winter and summer of 2021. Three apiaries were chosen in each region, and a pooled sample of 150 worker bees was collected from five colonies in each apiary then screened by qPCR for 10 viral targets: acute bee paralysis virus (ABPV), black queen cell virus (BQCV), chronic bee paralysis virus (CBPV), deformed wing virus (DWV) genotypes A (DWV-A), B (DWV-B) and D (Egyptian bee virus), Israeli acute paralysis virus (IAPV), Kashmir bee virus (KBV), sacbrood virus (SBV), and slow bee paralysis virus (SBPV). Our results revealed that DWV-A was the most prevalent virus, followed by BQCV and ABPV; the DWV genotype now spreading across the world, DWV-B, was not detected. There was no difference in varroa infestation rates as well as virus prevalence between winter and summer. However, colonies infected with BQCV had a significantly higher varroa count (adjusted p < 0.05) in the winter season, indicating that there is a seasonal association between the intensity of infestation by varroa and the presence of this virus. We provide data on the current virus prevalence in Egypt, which could assist in the protection of Egypt's beekeeping industry. Moreover, our study aids in the systematic assessment of the global honey bee virome by filling a knowledge gap about the prevalence of honey bee viruses in Egypt.

The vectoring competence of the mite Varroa destructor for deformed wing virus of honey bees is dynamic and affects survival of the mite

The ectoparasitic mite, Varroa destructor and the viruses it vectors, including types A and B of Deformed wing virus (DWV), pose a major threat to honey bees, Apis mellifera. Analysis of 256 mites collected from the same set of field colonies on five occasions from May to October 2021 showed that less than a half of them, 39.8% (95% confidence interval (CI): 34.0 - 46.0%), were able to induce a high (overt) level DWV infection with more than 109 viral genomes per bee in the pupa after 6 days of feeding, with both DWV-A and DWV-B being vectored at similar rates. To investigate the effect of the phoretic (or dispersal) stage on adult bees on the mites' ability to vector DWV, the mites from two collection events were divided into two groups, one of which was tested immediately for their infectiveness, and the other was kept with adult worker bees in cages for 12 days prior to testing their infectiveness. We found that while 39.2% (95% CI: 30.0 - 49.1%) of the immediately tested mites induced overt-level infections, 12-day passage on adult bees significantly increased the infectiousness to 89.8% (95% CI: 79.2 - 95.6%). It is likely that Varroa mites that survive brood interruptions in field colonies are increasingly infectious. The mite lifespan was affected by the DWV type it transmitted to pupae. The mites, which induced high DWV-B but not DWV-A infection had an average lifespan of 15.5 days (95% CI: 11.8 - 19.2 days), which was significantly shorter than those of the mites which induced high DWV-A but not DWV-B infection, with an average lifespan of 24.3 days (95% CI: 20.2 - 28.5), or the mites which did not induce high levels of DWV-A or DWV-B, with an average survival of 21.2 days (95% CI: 19.0 - 23.5 days). The mites which transmitted high levels of both DWV-A and DWV-B had an intermediate average survival of 20.5 days (95% CI: 15.1 - 25.9 days). The negative impact of DWV-B on mite survival could be a consequence of the ability of DWV-B, but not DWV-A to replicate in Varroa.

Quantitative and Qualitative Changes in the Deformed Wing Virus Population in Honey Bees Associated with the Introduction or Removal of Varroa destructor

Varroa destructor is an ectoparasitic mite associated with significant losses of honeybee colonies globally. The mite vectors a range of pathogenic viruses, the most important of which is the Deformed wing virus (DWV). In the absence of Varroa, DWV exists as a low-level, highly diverse virus population. However, when transmitted by Varroa, certain variants become highly elevated, and may become near-clonal and cause symptomatic infections. Mite transmission between colonies can occur when parasitised workers drift from or rob adjacent hives. These activities can result in elevated mite levels, but the resulting change in the DWV population, the primary determinant of winter colony losses, has not been determined. In reciprocal studies, we investigated the influence of the removal of mites, or their acquisition, on the DWV population. When mites were removed from heavily infested colonies, there was a striking and rapid reduction in virus load. Conversely, siting Varroa-naïve colonies in a mite-infested apiary resulted in the acquisition of mites and concomitant changes in the virus population. We observed both near-clonal and highly divergent virus populations regardless of titre, suggesting changes were stochastic and colony-specific. Our findings have implications for the outcome of strategies in areas with total or patchy implementation of Varroa control plans.