Dynamic evolution in the key honey bee pathogen deformed wing virus: Novel insights into virulence and competition using reverse genetics

Genotype B of deformed wing virus and related recombinant viruses become dominant in European honey bee colonies

The Varroa destructor mite's transmission of deformed wing virus (DWV) to honey bees is responsible for most winter mortalities of colonies worldwide. Four DWV genotypes (A, B, C and D) and numerous recombinants have been described. The most recent studies have reported the greater prevalence of DWV-B over DWV-A in several countries, including European ones, while C and D genotypes appear rare or extinct. However, no global evaluation of DWV-A and DWV-B distribution was available at the European level to date. In this study, we quantified both DWV genotypes by real-time PCR from pools or individual honey bees and from V. destructor mites sampled in 15 European countries between 2010 and 2017. These data and the sequencing of the viral RNA provide a first insight into DWV diversity, with a clear dominance of DWV-B and recombinants (A/B) in Europe. Chimeric sequencing reads were used to locate the recombinant junctions along the DWV genome. These were not randomly distributed, but mainly clustered in three genomic areas: the 5'UTR, leader peptide and helicase coding sequences. In our study, the DWV recombinant genomes shared at least the VP1-VP3 coding sequences with the DWV-B. Further studies are needed to explore the apicultural context explaining these differences in DWV genotype dominance.

Avoiding the tragedies of parasite tolerance in Darwinian beekeeping

Bee declines have been partly attributed to the impacts of invasive or emerging parasite outbreaks. For western honeybees, Apis mellifera, major losses are associated with the virus-vectoring mite, Varroa destructor. In response, beekeepers have focused breeding efforts aimed at conferring resistance to this key parasite. One method of many is survival-based beekeeping where colonies that survive despite significant Varroa infestations produce subsequent colonies. We argue that this 'hands-off' approach will not always lead to Varroa resistance evolving but rather tolerance. Tolerance minimizes host fitness costs of parasitism without reducing parasite abundance, whereas resistance either prevents parasitism outright or keeps parasitism intensity low. With clear epidemiological distinctions, and as honeybee disease dynamics impact other wild bees owing to shared pathogens, we discuss why tolerance outcomes in honeybee breeding have important implications for wider pollinator health. Crucially, we argue that unintentional selection for tolerance will not only lead to more spillover from honeybees but may also select for pathogens that are more virulent in wild bees leading to 'tragedies of tolerance'. These tragedies can be avoided through successful breeding regimes that specifically select for low Varroa. We emphasize how insights from evolutionary ecology can be applied in ecologically responsible honeybee management.

Deformed wing virus coopts the host arginine kinase to enhance its fitness in honey bees (Apis mellifera)

BACKGROUND

Deformed wing virus (DWV) is a major honey bee pathogen that is actively transmitted by the parasitic mite Varroa destructor and plays a primary role in Apis mellifera winter colony losses. Despite intense investigation on this pollinator, which has a unique environmental and economic importance, the mechanisms underlying the molecular interactions between DWV and honey bees are still poorly understood. Here, we report on a group of honey bee proteins, identified by mass spectrometry, that specifically co-immunoprecipitate with DWV virus particles.

RESULTS

Most of the proteins identified are involved in fundamental metabolic pathways. Among the co-immunoprecipitated proteins, one of the most interesting was arginine kinase (ArgK), a conserved protein playing multiple roles both in physiological and pathological processes and stress response in general. Here, we investigated in more detail the relationship between DWV and this protein. We found that argK RNA level positively correlates with DWV load in field-collected honey bee larvae and adults and significantly increases in adults upon DWV injection in controlled laboratory conditions, indicating that the argK gene was upregulated by DWV infection. Silencing argK gene expression in vitro, using RNAi, resulted in reduced DWV viral load, thus confirming that argK upregulation facilitates DWV infection, likely through interfering with the delicate balance between metabolism and immunity.

CONCLUSIONS

In summary, these data indicate that DWV modulates the host ArgK through transcriptional regulation and cooptation to enhance its fitness in honey bees. Our findings open novel perspectives on possible new therapies for DWV control by targeting specific host proteins.

Regional patterns and climatic predictors of viruses in honey bee (Apis mellifera) colonies over time

Honey bee viruses are serious pathogens that can cause poor colony health and productivity. We analyzed a multi-year longitudinal dataset of abundances of nine honey bee viruses (deformed wing virus A, deformed wing virus B, black queen cell virus, sacbrood virus, Lake Sinai virus, Kashmir bee virus, acute bee paralysis virus, chronic bee paralysis virus, and Israeli acute paralysis virus) in colonies located across Canada to describe broad trends in virus intensity and occurrence among regions and years. We also tested climatic variables (temperature, wind speed, and precipitation) as predictors in an effort to understand possible drivers underlying seasonal patterns in viral prevalence. Temperature was a significant positive predictor of the total number of viruses per sample, which was highest in British Columbia (mean = 5.0). Lake Sinai virus (LSV) was the most prevalent overall (at 89%) and had the highest infection intensity, at an average of 3.9 × 108 copies per bee. Acute bee paralysis virus was the least prevalent virus (at 4.7%) and had the lowest infection intensity (1.9 × 105 copies per bee). Surprisingly, including Varroa abundance as a covariate did not significantly improve model fit for any virus. All viruses, except Kashmir bee virus, varied by region, and one or more climatic variables were significant predictors for six of the nine viruses. Although climatic effects were often inconsistent among individual viruses, we show that climatic variables can be better predictors of virus intensity and occurrence than Varroa mite abundance, at least when infestation rates are low.

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.

New insights into honey bee viral and bacterial seasonal infection patterns using third-generation nanopore sequencing on honey bee haemolymph

Honey bees are rapidly declining, which poses a significant threat to our environment and agriculture industry. These vital insects face a disease complex believed to be caused by a combination of parasites, viruses, pesticides, and nutritional deficiencies. However, the real aetiology is still enigmatic. Due to the conventional analysis methods, we still lack complete insights into the honey bee virome and the presence of pathogenic bacteria. To fill this knowledge gap, we employed third-generation nanopore metagenomic sequencing on honey bee haemolymph to monitor the presence of pathogens over almost a year. This study provides valuable insights into the changes in bacterial and viral loads within honey bee colonies. We identified different pathogens in the honey bee haemolymph, which are not included in honey bee screenings. These pathogens comprise the Apis mellifera filamentous virus, Apis rhabdoviruses, and various bacteria such as Frischella sp. and Arsenophonus sp. Furthermore, a sharp contrast was observed between young and old bees. Our research proposes that transgenerational immune priming may play a role in shaping infection patterns in honey bees. We observed a significant increase in pathogen loads in the spring, followed by a notable decrease in pathogen presence during the summer and autumn months. However, certain pathogens seem to be able to evade this priming effect, making them particularly intriguing as potential factors contributing to mortality. In the future, we aim to expand our research on honey bee transgenerational immune priming and investigate its potential in natural settings. This knowledge will ultimately enhance honey bee health and decrease colony mortality.

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 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.

Ontogeny of immunity and natural viral infection in Apis mellifera drones and workers

The most common viral diseases affecting honey bees (Apis mellifera) in Israel include deformed wing viruses (DWV-A and DWV-B) and acute paralysis viruses (ABPV and IAPV). These viruses are transmitted within and between colonies, both horizontally and vertically. All members of the colony contribute to this transmission, on the other hand individual and social immunity, particularly hygienic behaviour, may affect the outcome of the process. In this study, we evaluated the ontogeny of natural infections of DWV-A, DWV-B, ABPV and IAPV, their prevalence and loads, in workers and drones from high (H) and low (L) hygienic colonies. In parallel, we evaluated the expression of two immune genes: peptidoglycan recognition protein S2(PGRP-S2) and hymenoptaecin. The prevalence of DWV-B and IAPV increased with age and was higher in workers than in drones. ABPV was not detected in drones. The expression of both immune genes was significantly affected by age and sex. Drones from H colonies had higher expression of these genes. The increased expression of immune genes with drones' age, particularly in hygienic colonies, suggest additional value of honey bee breeding for hygienic behaviour for sustainable beekeeping.

Data mining reveals tissue-specific expression and host lineage-associated forms of Apis mellifera filamentous virus

Background

Apis mellifera filamentous virus (AmFV) is a large double-stranded DNA virus of uncertain phylogenetic position that infects honey bees (Apis mellifera). Little is known about AmFV evolution or molecular aspects of infection. Accurate annotation of open-reading frames (ORFs) is challenged by weak homology to other known viruses. This study was undertaken to evaluate ORFs (including coding-frame conservation, codon bias, and purifying selection), quantify genetic variation within AmFV, identify host characteristics that covary with infection rate, and examine viral expression patterns in different tissues.

Methods

Short-read data were accessed from the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI). Sequence reads were downloaded from accessions meeting search criteria and scanned for kmers representative of AmFV genomic sequence. Samples with kmer counts above specified thresholds were downloaded in full for mapping to reference sequences and de novo assembly.

Results

At least three distinct evolutionary lineages of AmFV exist. Clade 1 predominates in Europe but in the Americas and Africa it is replaced by the other clades as infection level increases in hosts. Only clade 3 was found at high relative abundance in hosts with African ancestry, whereas all clades achieved high relative abundance in bees of non-African ancestry. In Europe and Africa, clade 2 was generally detected only in low-level infections but was locally dominant in some North American samples. The geographic distribution of clade 3 was consistent with an introduction to the Americas with 'Africanized' honey bees in the 1950s. Localized genomic regions of very high nucleotide divergence in individual isolates suggest recombination with additional, as-yet unidentified AmFV lineages. A set of 155 high-confidence ORFs was annotated based on evolutionary conservation in six AmFV genome sequences representative of the three clades. Pairwise protein-level identity averaged 94.6% across ORFs (range 77.1-100%), which generally exhibited low evolutionary rates and moderate to strong codon bias. However, no robust example of positive diversifying selection on coding sequence was found in these alignments. Most of the genome was detected in RNA short-read alignments. Transcriptome assembly often yielded contigs in excess of 50 kb and containing ORFs in both orientations, and the termini of long transcripts were associated with tandem repeats. Lower levels of AmFV RNA were detected in brain tissue compared to abdominal tissue, and a distinct set of ORFs had minimal to no detectable expression in brain tissue. A scan of DNA accessions from the parasitic mite Varroa destructor was inconclusive with respect to replication in that species.

Discussion

Collectively, these results expand our understanding of this enigmatic virus, revealing transcriptional complexity and co-evolutionary associations with host lineage.

Deformed wing virus of honey bees is inactivated by cold plasma ionized hydrogen peroxide

Deformed wing virus (DWV) is a widespread pathogen of Apis mellifera honey bees, and is considered a major causative factor for the collapse of infected honey bee colonies. DWV can be horizontally transmitted among bees through various oral routes, including via food sharing and by interactions of bees with viral-contaminated solid hive substrates. Cold plasma ionized hydrogen peroxide (iHP) is used extensively by the food production, processing and medical industries to clean surfaces of microbial contaminants. In this study, we investigated the use of iHP to inactivate DWV particles in situ on a solid substrate. iHP-treated DWV sources were ~105-fold less infectious when injected into naïve honey bee pupae compared to DWV receiving no iHP treatment, matching injected controls containing no DWV. iHP treatment also greatly reduced the incidence of overt DWV infections (i.e., pupae having >109 copies of DWV). The level of DWV inactivation achieved with iHP treatment was higher than other means of viral inactivation such as gamma irradiation, and iHP treatment is likely simpler and safer. Treatment of DWV contaminated hive substrates with iHP, even with honey bees present, may be an effective way to decrease the impacts of DWV infection on honey bees.

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.

Comparative analysis of 3 pollen sterilization methods for feeding bumble bees

Pollen is an essential component of bee diets, and rearing bumble bees (Bombus spp.) for commercial use necessitates feeding pollen in mass quantities. This pollen is collected from honey bee (Apis mellifera L.) colonies because neither an artificial diet nor an economical, large-scale pollen collection process from flowers is available. The provenance of honey bee-collected pollen is often unknown, and in some cases has crossed international borders. Both deformed wing virus (DWV) and the fungal pathogen Ascosphaera apis (Claussen) Olive & Spiltoir (cause of chalkbrood disease); occur in honey bee-collected pollen, and infections have been observed in bumble bees. We used these pathogens as general surrogates for viruses and spore-forming fungal diseases to test the efficacy of 3 sterilization methods, and assessed whether treatment altered pollen quality for the bumble bee. Using honey bee-collected pollen spiked with known doses of DWV and A. apis, we compared gamma irradiation (GI), ozone fumigation (OZ), and ethylene oxide fumigation (EO) against an untreated positive control and a negative control. Following sterilization treatments, we tested A. apis spore viability, detected viral presence with PCR, and tested palatability to the bumble bee Bombus impatiens Cresson. We also measured bacterial growth from pollens treated with EO and GI. GI and EO outperformed OZ treatment in pathogen suppression. EO had the highest sterilizing properties under commercial conditions and retained palatability and supported bee development better than other treatments. These results suggest that EO sterilization reduces pathogen risks while retaining pollen quality as a food source for rearing bumble bees.

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.

Molecular Detection and Phylogenetic Analysis of Deformed Wing Virus and Sacbrood Virus Isolated from Pollen

Among many pathogens and pests, honey bee viruses are known as one of the most common cause of diseases in honey bee colonies. In this study, we demonstrate that pollen grains and bee bread are potential sources of viral DNA. We extracted DNA from 3 types of pollen samples: directly provided by beekeepers (n = 12), purchased from trade markets (n = 5), and obtained from honeycombs (bee bread, n = 10). The extracted DNA was used for molecular detection (RT-PCR analysis) of six of the most widely distributed honey bee viruses: deformed wing virus, sacbrood virus, acute bee paralysis virus, black queen cell virus, Kashmir bee virus, Israeli acute paralysis virus, and chronic bee paralysis virus. We successfully managed to establish only the deformed wing virus (DWV) and the sacbrood virus (SBV), with different distribution frequencies depending on the territory of the country. The phylogenetic analyses of Bulgarian isolates were performed with the most similar sequences available in molecular databases from other countries. Phylogenies of Bulgarian viral strains demonstrated genetically heterogeneous populations of DWV and relatively homogenous populations of SBV. In conclusion, the results obtained from the current study have shown that pollen is a valuable source for molecular detection of honey bee pathogens. This allows epidemiological monitoring of honey bee diseases at a regional and a national level.

Promiscuous feeding on multiple adult honey bee hosts amplifies the vectorial capacity of Varroa destructor

Varroa destructor is a cosmopolitan pest and leading cause of colony loss of the European honey bee. Historically described as a competent vector of honey bee viruses, this arthropod vector is the cause of a global pandemic of Deformed wing virus, now endemic in honeybee populations in all Varroa-infested regions. Our work shows that viral spread is driven by Varroa actively switching from one adult bee to another as they feed. Assays using fluorescent microspheres were used to indicate the movement of fluids in both directions between host and vector when Varroa feed. Therefore, Varroa could be in either an infectious or naïve state dependent upon the disease status of their host. We tested this and confirmed that the relative risk of a Varroa feeding depended on their previous host's infectiousness. Varroa exhibit remarkable heterogeneity in their host-switching behavior, with some Varroa infrequently switching while others switch at least daily. As a result, relatively few of the most active Varroa parasitize the majority of bees. This multiple-feeding behavior has analogs in vectorial capacity models of other systems, where promiscuous feeding by individual vectors is a leading driver of vectorial capacity. We propose that the honeybee-Varroa relationship offers a unique opportunity to apply principles of vectorial capacity to a social organism, as virus transmission is both vectored and occurs through multiple host-to-host routes common to a crowded society.

Metagenomic analysis of viromes in honey bee colonies (Apis mellifera; Hymenoptera: Apidae) after mass disappearance in Korea

After the nationwide, massive winter losses of honey bees in Korea during the winter of 2021, samplings were conducted from live honey bees in colonies and dead honey bees nearby colonies in the same bee-farms in six regions in Korea. Each sample was subjected to virome analysis using high-throughput sequencing technology. The number of viral reads was the lowest in the live honey bee group sample with 370,503 reads and the highest in the dead honey bee group sample with 42,659,622 reads. Viral contigs were matched with the viral genomes of the black queen cell virus, deformed wing virus, Israeli acute paralysis virus, and sacbrood virus, all of which have been previously reported in Korea. However, Apis rhabdovirus 5, bee macula-like virus, Varroa orthomyxovirus-1, Hubei partiti-like virus 34, Lake Sinai virus 2, 3, and 4, and the Ditton virus, were also discovered in this study, which are the first records in Korea. Plant viral sequences resembling those of Arabidopsis latent virus 1, and a novel viral sequence was also discovered. In the present study 55 complete viral genome sequences were identified. This study is the first virome analysis of domestic honey bees and provides the latest information on the diversity of honey bee viruses in Korea.

Virus Prevalence in Egg Samples Collected from Naturally Selected and Traditionally Managed Honey Bee Colonies across Europe

Monitoring virus infections can be an important selection tool in honey bee breeding. A recent study pointed towards an association between the virus-free status of eggs and an increased virus resistance to deformed wing virus (DWV) at the colony level. In this study, eggs from both naturally surviving and traditionally managed colonies from across Europe were screened for the prevalence of different viruses. Screenings were performed using the phenotyping protocol of the 'suppressed in ovo virus infection' trait but with qPCR instead of end-point PCR and a primer set that covers all DWV genotypes. Of the 213 screened samples, 109 were infected with DWV, 54 were infected with black queen cell virus (BQCV), 3 were infected with the sacbrood virus, and 2 were infected with the acute bee paralyses virus. It was demonstrated that incidences of the vertical transmission of DWV were more frequent in naturally surviving than in traditionally managed colonies, although the virus loads in the eggs remained the same. When comparing virus infections with queen age, older queens showed significantly lower infection loads of DWV in both traditionally managed and naturally surviving colonies, as well as reduced DWV infection frequencies in traditionally managed colonies. We determined that the detection frequencies of DWV and BQCV in honey bee eggs were lower in samples obtained in the spring than in those collected in the summer, indicating that vertical transmission may be lower in spring. Together, these patterns in vertical transmission show that honey bee queens have the potential to reduce the degree of vertical transmission over time.

Procedures and potential pitfalls for constructing a bee-infecting RNA virus clone

Viruses are factors that can fluctuate insect populations, including honey bees. Most honey bee infecting viruses are single positive-stranded RNA viruses that may not specifically infect honey bees and can be hazardous to other pollinator insects. In addition, these viruses could synergize with other stressors to worsen the honey bee population decline. To identify the underlying detailed mechanisms, reversed genetic studies with infectious cDNA clones of the viruses are necessary. Moreover, an infectious cDNA clone can be applied to studies as an ideal virus isolate that consists of a single virus species with a uniform genotype. However, only a few infectious cDNA clones have been reported in honey bee studies since the first infectious cDNA clone was published four decades ago. This article discusses steps, rationales, and potential issues in bee-infecting RNA virus cloning. In addition, failed experiences of cloning a Deformed wing virus isolate that was phylogenetically identical to Kakugo virus were addressed. We hope the information provided in this article can facilitate further developments of reverse-genetic studies of bee-infecting viruses to clarify the roles of virus diseases in the current pollinator declines.

Transmission of deformed wing virus between Varroa destructor foundresses, mite offspring and infested honey bees

Background

Varroa destructor is the major ectoparasite of the western honey bee (Apis mellifera). Through both its parasitic life-cycle and its role as a vector of viral pathogens, it can cause major damage to honey bee colonies. The deformed wing virus (DWV) is the most common virus transmitted by this ectoparasite, and the mite is correlated to increased viral prevalence and viral loads in infested colonies. DWV variants A and B (DWV-A and DWV-B, respectively) are the two major DWV variants, and they differ both in their virulence and transmission dynamics.

Methods

We studied the transmission of DWV between bees, parasitic mites and their offspring by quantifying DWV loads in bees and mites collected in in vitro and in situ environments. In vitro, we artificially transmitted DWV-A to mites and quantified both DWV-A and DWV-B in mites and bees. In situ, we measured the natural presence of DWV-B in bees, mites and mites’ offspring.

Results

Bee and mite viral loads were correlated, and mites carrying both variants were associated with higher mortality of the infected host. Mite infestation increased the DWV-B loads and decreased the DWV-A loads in our laboratory conditions. In situ, viral quantification in the mite offspring showed that, after an initially non-infected egg stage, the DWV-B loads were more closely correlated with the foundress (mother) mites than with the bee hosts.

Conclusions

The association between mites and DWV-B was highlighted in this study. The parasitic history of a mite directly impacts its DWV infection potential during the rest of its life-cycle (in terms of variant and viral loads). Regarding the mite’s progeny, we hypothesize that the route of contamination is likely through the feeding site rather than by vertical transmission, although further studies are needed to confirm this hypothesis.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05463-9.

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.

Epidemiology of the Microsporidium Nosema ceranae in Four Mediterranean Countries

Nosema ceranae is a highly prevalent intracellular parasite of honey bees' midgut worldwide. This Microsporidium was monitored during a long-term study to evaluate the infection at apiary and intra-colony levels in six apiaries in four Mediterranean countries (France, Israel, Portugal, and Spain). Parameters on colony strength, honey production, beekeeping management, and climate were also recorded. Except for São Miguel (Azores, Portugal), all apiaries were positive for N. ceranae, with the lowest prevalence in mainland France and the highest intra-colony infection in Israel. A negative correlation between intra-colony infection and colony strength was observed in Spain and mainland Portugal. In these two apiaries, the queen replacement also influenced the infection levels. The highest colony losses occurred in mainland France and Spain, although they did not correlate with the Nosema infection levels, as parasitism was low in France and high in Spain. These results suggest that both the effects and the level of N. ceranae infection depends on location and beekeeping conditions. Further studies on host-parasite coevolution, and perhaps the interactions with other pathogens and the role of honey bee genetics, could assist in understanding the difference between nosemosis disease and infection, to develop appropriate strategies for its control.

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

The western honey bee (Apis mellifera) is of major economic and ecological importance, with elevated rates of colony losses in temperate regions over the last two decades thought to be largely caused by the exotic ectoparasitic mite Varroa destructor and deformed wing virus (DWV), which the mite transmits. DWV currently exists as two main genotypes: the formerly widespread DWV-A and the more recently described and rapidly expanding DWV-B. It is an excellent system to understand viral evolution and the replacement of one viral variant by another. Here we synthesise published results on the distribution and prevalence of DWV-A and -B over the period 2008-2021 and present novel data for Germany, Italy and the UK to suggest that (i) DWV-B has rapidly expanded worldwide since its first description in 2004 and (ii) that it is potentially replacing DWV-A. Both genotypes are also found in wild bee species. Based on a simple mathematical model, we suggest that interference between viral genotypes when co-infecting the same host is key to understanding their epidemiology. We finally discuss the consequences of genotype replacement for beekeeping and for wild pollinator species.

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.

Two Novel Iflaviruses Discovered in Bat Samples in Washington State

Arthropods are integral to ecosystem equilibrium, serving as both a food source for insectivores and supporting plant reproduction. Members of the Iflaviridae family in the order Picornavirales are frequently found in RNA sequenced from arthropods, who serve as their hosts. Here we implement a metagenomic deep sequencing approach followed by rapid amplification of cDNA ends (RACE) on viral RNA isolated from wild and captured bat guano in Washington State at two separate time points. From these samples we report the complete genomes of two novel viruses in the family Iflaviridae. The first virus, which we call King virus, is 46% identical by nucleotide to the lethal honeybee virus, deformed wing virus, while the second virus which we call Rolda virus, shares 39% nucleotide identity to deformed wing virus. King and Rolda virus genomes are 10,183 and 8934 nucleotides in length, respectively. Given these iflaviruses were detected in guano from captive bats whose sole food source was the Tenebrio spp. mealworm, we anticipate this invertebrate may be a likely host. Using the NCBI Sequence Read Archive, we found that these two viruses are located in six continents and have been isolated from a variety of arthropod and mammalian specimens.

Molecular Detection and Differentiation of Arthropod, Fungal, Protozoan, Bacterial and Viral Pathogens of Honeybees

The honeybee Apis mellifera is highly appreciated worldwide because of its products, but also as it is a pollinator of crops and wild plants. The beehive is vulnerable to infections due to arthropods, fungi, protozoa, bacteria and/or viruses that manage to by-pass the individual and social immune mechanisms of bees. Due to the close proximity of bees in the beehive and their foraging habits, infections easily spread within and between beehives. Moreover, international trade of bees has caused the global spread of infections, several of which result in significant losses for apiculture. Only in a few cases can infections be diagnosed with the naked eye, by direct observation of the pathogen in the case of some arthropods, or by pathogen-associated distinctive traits. Development of molecular methods based on the amplification and analysis of one or more genes or genomic segments has brought significant progress to the study of bee pathogens, allowing for: (i) the precise and sensitive identification of the infectious agent; (ii) the analysis of co-infections; (iii) the description of novel species; (iv) associations between geno- and pheno-types and (v) population structure studies. Sequencing of bee pathogen genomes has allowed for the identification of new molecular targets and the development of specific genotypification strategies.

A derived honey bee stock confers resistance to Varroa destructor and associated viral transmission

The ectoparasite Varroa destructor is the greatest threat to managed honey bee (Apis mellifera) colonies globally. Despite significant efforts, novel treatments to control the mite and its vectored pathogens have shown limited efficacy, as the host remains naïve. A prospective solution lies in the development of Varroa-resistant honey bee stocks, but a paucity of rigorous selection data restricts widespread adoption. Here, we characterise the parasite and viral dynamics of a Varroa-resistant honey bee stock, designated ‘Pol-line’, using a large-scale longitudinal study. Results demonstrate markedly reduced Varroa levels in this stock, diminished titres of three major viruses (DWV-A, DWV-B, and CBPV), and a two-fold increase in survival. Levels of a fourth virus that is not associated with Varroa—BQCV—do not differ between stocks, supporting a disruption of the transmission pathway. Further, we show that when decoupled from the influence of Varroa levels, viral titres do not constitute strong independent predictors of colony mortality risk. These findings highlight the need for a reassessment of Varroa etiology, and suggest that derived stocks represent a tractable solution to the Varroa pandemic.

An Insight Into the microRNA Profile of the Ectoparasitic Mite Varroa destructor (Acari: Varroidae), the Primary Vector of Honey Bee Deformed Wing Virus

The remarkably adaptive mite Varroa destructor is the most important honey bee ectoparasite. Varroa mites are competent vectors of deformed wing virus (DWV), and the Varroa-virus complex is a major determinant of annual honey bee colony mortality and collapse. MicroRNAs (miRNAs) are 22-24 nucleotide non-coding RNAs produced by all plants and animals and some viruses that influence biological processes through post-transcriptional regulation of gene expression. Knowledge of miRNAs and their function in mite biology remains limited. Here we constructed small RNA libraries from male and female V. destructor using Illumina's small RNA-Seq platform. A total of 101,913,208 and 91,904,732 small RNA reads (>18 nucleotides) from male and female mites were analyzed using the miRDeep2 algorithm. A conservative approach predicted 306 miRNAs, 18 of which were upregulated and 13 downregulated in female V. destructor compared with males. Quantitative real-time PCR validated the expression of selected differentially-expressed female Varroa miRNAs. This dataset provides a list of potential miRNA targets involved in regulating vital Varroa biological processes and paves the way for developing strategies to target Varroa and their viruses.

Critical View on the Importance of Host Defense Strategies on Virus Distribution of Bee Viruses: What Can We Learn from SARS-CoV-2 Variants?

Bees, both wild and domesticated ones, are hosts to a plethora of viruses, with most of them infecting a wide range of bee species and genera. Although viral discovery and research on bee viruses date back over 50 years, the last decade is marked by a surge of new studies, new virus discoveries, and reports on viral transmission in and between bee species. This steep increase in research on bee viruses was mainly initiated by the global reports on honeybee colony losses and the worldwide wild bee decline, where viruses are regarded as one of the main drivers. While the knowledge gained on bee viruses has significantly progressed in a short amount of time, we believe that integration of host defense strategies and their effect on viral dynamics in the multi-host viral landscape are important aspects that are currently still missing. With the large epidemiological dataset generated over the last two years on the SARS-CoV-2 pandemic, the role of these defense mechanisms in shaping viral dynamics has become eminent. Integration of these dynamics in a multi-host system would not only greatly aid the understanding of viral dynamics as a driver of wild bee decline, but we believe bee pollinators and their viruses provide an ideal system to study the multi-host viruses and their epidemiology.

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.

Haplotype Analysis of Varroa destructor and Deformed Wing Virus Using Long Reads

As a phoretic parasite and virus vector, the mite Varroa destructor and the associated Deformed wing virus (DWV) form a lethal combination to the honey bee, Apis mellifera. Routine acaricide treatment has been reported to reduce the diversity of mites and select for tolerance against these treatments. Further, different DWV strains face selective pressures when transmitted via mites. In this study, the haplotypes of Varroa mites and associated DWV variants were quantified using long reads. A single haplotype dominated the mite mitochondrial gene cytochrome oxidase subunit I, reflecting an ancient bottleneck. However, highly polymorphic genes were present across the mite genome, suggesting the diversity of mites could be actively maintained at a regional level. DWV detected in both mites and honey bees show a dominant variant with only a few low-frequency alternate haplotypes. The relative abundances of DWV haplotypes isolated from honey bees and mites were highly consistent, suggesting that some variants are favored by ongoing selection.

First come, first served: superinfection exclusion in Deformed wing virus is dependent upon sequence identity and not the order of virus acquisition

Deformed wing virus (DWV) is the most important globally distributed pathogen of honey bees and, when vectored by the ectoparasite Varroa destructor, is associated with high levels of colony losses. Divergent DWV types may differ in their pathogenicity and are reported to exhibit superinfection exclusion upon sequential infections, an inevitability in a Varroa-infested colony. We used a reverse genetic approach to investigate competition and interactions between genetically distinct or related virus strains, analysing viral load over time, tissue distribution with reporter gene-expressing viruses and recombination between virus variants. Transient competition occurred irrespective of the order of virus acquisition, indicating no directionality or dominance. Over longer periods, the ability to compete with a pre-existing infection correlated with the genetic divergence of the inoculae. Genetic recombination was observed throughout the DWV genome with recombinants accounting for ~2% of the population as determined by deep sequencing. We propose that superinfection exclusion, if it occurs at all, is a consequence of a cross-reactive RNAi response to the viruses involved, explaining the lack of dominance of one virus type over another. A better understanding of the consequences of dual- and superinfection will inform development of cross-protective honey bee vaccines and landscape-scale DWV transmission and evolution.

Global similarity, and some key differences, in the metagenomes of Swedish varroa-surviving and varroa-susceptible honeybees

There is increasing evidence that honeybees (Apis mellifera L.) can adapt naturally to survive Varroa destructor, the primary cause of colony mortality world-wide. Most of the adaptive traits of naturally varroa-surviving honeybees concern varroa reproduction. Here we investigate whether factors in the honeybee metagenome also contribute to this survival. The quantitative and qualitative composition of the bacterial and viral metagenome fluctuated greatly during the active season, but with little overall difference between varroa-surviving and varroa-susceptible colonies. The main exceptions were Bartonella apis and sacbrood virus, particularly during early spring and autumn. Bombella apis was also strongly associated with early and late season, though equally for all colonies. All three affect colony protein management and metabolism. Lake Sinai virus was more abundant in varroa-surviving colonies during the summer. Lake Sinai virus and deformed wing virus also showed a tendency towards seasonal genetic change, but without any distinction between varroa-surviving and varroa-susceptible colonies. Whether the changes in these taxa contribute to survival or reflect demographic differences between the colonies (or both) remains unclear.

Simulated vector transmission differentially influences dynamics of two viral variants of deformed wing virus in honey bees (Apis mellifera)

Understanding how vectors alter the interactions between viruses and their hosts is a fundamental question in virology and disease ecology. In honey bees, transmission of deformed wing virus (DWV) by parasitic Varroa mites has been associated with elevated disease and host mortality, and Varroa transmission has been hypothesized to lead to increased viral titres or select for more virulent variants. Here, we mimicked Varroa transmission by serially passaging a mixed population of two DWV variants, A and B, by injection through in vitro reared honey bee pupae and tracking these viral populations through five passages. The DWV-A and DWV-B variant proportions shifted dynamically through passaging, with DWV-B outcompeting DWV-A after one passage, but levels of both variants becoming equivalent by Passage 5. Sequencing analysis revealed a dominant, recombinant DWV-B strain (DWV-A derived 5' IRES region with the rest of the genome DWV-B), with low nucleotide diversity that decreased through passaging. DWV-A populations had higher nucleotide diversity compared to DWV-B, but this also decreased through passaging. Selection signatures were found across functional regions of the DWV-A and DWV-B genomes, including amino acid mutations in the putative capsid protein region. Simulated vector transmission differentially impacted two closely related viral variants which could influence viral interactions with the host, demonstrating surprising plasticity in vector-host-viral dynamics.

Beeporter: Tools for high-throughput analyses of pollinator-virus infections

Pollinators are in decline thanks to the combined stresses of disease, pesticides, habitat loss, and climate. Honey bees face numerous pests and pathogens but arguably none are as devastating as Deformed wing virus (DWV). Understanding host-pathogen interactions and virulence of DWV in honey bees is slowed by the lack of cost-effective high-throughput screening methods for viral infection. Currently, analysis of virus infection in bees and their colonies is tedious, requiring a well-equipped molecular biology laboratory and the use of hazardous chemicals. Here we describe virus clones tagged with green fluorescent protein (GFP) or nanoluciferase (nLuc) that provide high-throughput detection and quantification of virus infections. GFP fluorescence is measured noninvasively in living bees via commonly available long-wave UV light sources and a smartphone camera, or a standard ultraviolet transilluminator gel imaging system. Nonlethal monitoring with GFP allows continuous screening of virus growth and serves as a direct breeding tool for identifying honey bee parents with increased antiviral resistance. Expression using the nLuc reporter strongly correlates with virus infection levels and is especially sensitive. Using multiple reporters, it is also possible to visualize competition, differential virulence, and host tissue targeting by co-occuring pathogens. Finally, it is possible to directly assess the risk of cross-species "spillover" from honey bees to other pollinators and vice versa.

Variant A of the Deformed Wings Virus Alters the Olfactory Sensitivity and the Expression of Odorant Binding Proteins on Antennas of Apis mellifera

Simple Summary

Honey bees, Apis melllifera, are the most commonly managed bee in the world for pollination services. However, worldwide continuous colony losses have been reported for almost a decade. One factor of these losses is associated to pathogens being the virus one of the most important problems in honey bee health. One of the known viruses that affect the honey bee population is deformed wing virus (DWV). DWV causes physical malformation and behavioral disturbances, but also, this virus can be found in the antenna affecting the anatomical integrity of infected areas, which could compromise normal antennal functioning associated to aroma perception. Thus, we evaluate olfactory sensitivity and the expression of antenna-specific odorant-binding proteins (OBP) genes in honey bees inoculated with variant A of the DWV. We performed olfactory sensitivity analysis using the essential oils Eucalyptus globulus and Mentha piperita, but also, and molecular analysis of gene expression of nine OBPs. We found that the high level of replication of DWV-A in the antennae decreased the olfactory sensitivity and led to a down-regulation of some OBPs in middle- and forager-age worker bees. Thus, DWV-A infection in adults of honey bees could compromise volatile compound recognition inside the hive and outside the hive.

Abstract

Insects have a highly sensitive sense of smell, allowing them to perform complex behaviors, such as foraging and peer recognition. Their sense of smell is based on the recognition of ligands and is mainly coordinated by odorant-binding proteins (OBPs). In Apis mellifera, behavior can be affected by different pathogens, including deformed wing virus (DWV) and its variants. In particular, it has been shown that variant A of DWV (DWV-A) is capable of altering the ultra-cellular structure associated with olfactory activity. In this study was evaluated olfactory sensitivity and the expression of OBP genes in honey bees inoculated with DWV-A. Electroantennographic analyses (EAG) were carried out to determine the olfactory sensitivity to the essential oils Eucalyptus globulus and Mentha piperita. The expression of nine antenna-specific OBP genes and DWV-A load in inoculated bees was also quantified by qPCR. We observed an inverse relationship between viral load and olfactory sensitivity and the expression of some OBP proteins. Thus, high viral loads reduced olfactory sensitivity to essential oils and the gene expression of the OBP2, OBP5, OBP11, and OBP12 proteins on the antennas of middle- and forager-age bees. These results suggest that DWV-A could have negative effects on the processes of aroma perception by worker bees, affecting their performance in tasks carried out in and outside the colony.

A Survey of Common Pathogens of Apis spp. in Wild Non-Cave Honeybees in Southwest China

Honeybees are crucial pollinators with significant ecologic value. The decline of wild honeybee populations has been recognized and documented during recent decades. However, the health status of wild non-cave Apis spp., including giant and dwarf honeybees, remains generally unknown. We investigated eight common viruses and five bacterial or fungal pathogens in four wild non-cave honeybee species at 11 locations in Southwest China. As a result, Melissococcus plutonius, the pathogenic agent of European foulbrood, was detected in all the species, and the sequences were identical to the pathogen in managed cave honeybees. Only one virus, black queen cell virus (BQCV), was positive in one dwarf species, Apis florea, in our study. The positive BQCV infected three A. florea colonies in Guangxi Province, with distinct sequences from this virus reported in cave honeybees or in the same host in the nearby Yunnan Province. Although our results indicated a low pathogenic level of common diseases in the wild non-cave Apis spp. in Southwest China, the conservation of these wild pollinators is of importance in light of the noticeable decline in populations and the irreplaceable position of pollination.

Breeding for Virus Resistance and Its Effects on Deformed Wing Virus Infection Patterns in Honey Bee Queens

Viruses, and in particular the deformed wing virus (DWV), are considered as one of the main antagonists of honey bee health. The 'suppressed in ovo virus infection' trait (SOV) described for the first time that control of a virus infection can be achieved from genetically inherited traits and that the virus state of the eggs is indicative for this. This research aims to explore the effect of the SOV trait on DWV infections in queens descending from both SOV-positive (QDS+) and SOV-negative (QDS-) queens. Twenty QDS+ and QDS- were reared from each time four queens in the same starter-finisher colony. From each queen the head, thorax, ovaries, spermatheca, guts and eviscerated abdomen were dissected and screened for the presence of the DWV-A and DWV-B genotype using qRT-PCR. Queens descending from SOV-positive queens showed significant lower infection loads for DWV-A and DWV-B as well as a lower number of infected tissues for DWV-A. Surprisingly, differences were less expressed in the reproductive tissues, the ovaries and spermatheca. These results confirm that selection on the SOV trait is associated with increased virus resistance across viral genotypes and that this selection drives DWV towards an increased tissue specificity for the reproductive tissues. Further research is needed to explore the mechanisms underlying the interaction between the antiviral response and DWV.

Pupal cannibalism by worker honey bees contributes to the spread of deformed wing virus

Transmission routes impact pathogen virulence and genetics, therefore comprehensive knowledge of these routes and their contribution to pathogen circulation is essential for understanding host-pathogen interactions and designing control strategies. Deformed wing virus (DWV), a principal viral pathogen of honey bees associated with increased honey bee mortality and colony losses, became highly virulent with the spread of its vector, the ectoparasitic mite Varroa destructor. Reproduction of Varroa mites occurs in capped brood cells and mite-infested pupae from these cells usually have high levels of DWV. The removal of mite-infested pupae by worker bees, Varroa Sensitive Hygiene (VSH), leads to cannibalization of pupae with high DWV loads, thereby offering an alternative route for virus transmission. We used genetically tagged DWV to investigate virus transmission to and between worker bees following pupal cannibalisation under experimental conditions. We demonstrated that cannibalization of DWV-infected pupae resulted in high levels of this virus in worker bees and that the acquired virus was then transmitted between bees via trophallaxis, allowing circulation of Varroa-vectored DWV variants without the mites. Despite the known benefits of hygienic behaviour, it is possible that higher levels of VSH activity may result in increased transmission of DWV via cannibalism and trophallaxis.

Exposure to Inactivated Deformed Wing Virus Leads to Trans-Generational Costs but Not Immune Priming in Honeybees (Apis mellifera)

Pathogens are identified as one of the major drivers behind the honeybee colony losses, as well as one of the reasons for the reported declines in terrestrial insect abundances in recent decades. To fight infections, animals rely on their immune system. The immune system of many invertebrates can be primed by exposure to a pathogen, so that upon further exposure the animal is better protected. The protective priming effect can even extend to the next generation, but the species capable of priming the immune system of their offspring are still being investigated. Here we studied whether honeybees could prime their offspring against a viral pathogen, by challenging honeybee queens orally with an inactivated deformed wing virus (DWV), one of the most devastating honeybee viruses. The offspring were then infected by viral injection. The effects of immune priming were assayed by measuring viral loads and two typical symptoms of the virus, pupal mortality, and abnormal wing phenotype. We saw a low amount of wing deformities and low pupal mortality. While no clear priming effect against the virus was seen, we found that the maternal immune challenge, when combined with the stress caused by an injection during development, manifested in costs in the offspring, leading to an increased number of deformed wings.

Virus Prospecting in Crickets-Discovery and Strain Divergence of a Novel Iflavirus in Wild and Cultivated Acheta domesticus

Orthopteran insects have high reproductive rates leading to boom-bust population dynamics with high local densities that are ideal for short, episodic disease epidemics. Viruses are particularly well suited for such host population dynamics, due to their supreme ability to adapt to changing transmission criteria. However, very little is known about the viruses of Orthopteran insects. Since Orthopterans are increasingly reared commercially, for animal feed and human consumption, there is a risk that viruses naturally associated with these insects can adapt to commercial rearing conditions, and cause disease. We therefore explored the virome of the house cricket Acheta domesticus, which is both part of the natural Swedish landscape and reared commercially for the pet feed market. Only 1% of the faecal RNA and DNA from wild-caught A. domesticus consisted of viruses. These included both known and novel viruses associated with crickets/insects, their bacterial-fungal microbiome, or their plant food. Relatively abundant among these viral Operational Taxonomic Units (OTUs) was a novel Iflavirus, tentatively named Acheta domesticus Iflavirus (AdIV). Quantitative analyses showed that AdIV was also abundant in frass and insect samples from commercially reared crickets. Interestingly, the wild and commercial AdIV strains had short, extremely divergent variation hotspots throughout the genome, which may indicate specific adaptation to their hosts’ distinct rearing environments.

Direct Evidence for Infection of Varroa destructor Mites with the Bee-Pathogenic Deformed Wing Virus Variant B - but Not Variant A - via Fluorescence-in situ-Hybridization Analysis

Deformed wing virus (DWV) is a bee pathogenic, single- and positive-stranded RNA virus that has been involved in severe honey bee colony losses worldwide. DWV, when transmitted horizontally or vertically from bee to bee, causes mainly covert infections not associated with any visible symptoms or damage. Overt infections occur after vectorial transmission of DWV to the developing bee pupae through the ectoparasitic mite Varroa destructor Symptoms of overt infections are pupal death, bees emerging with deformed wings and shortened abdomens, or cognitive impairment due to brain infection. So far, three variants of DWV, DWV-A, DWV-B, and DWV-C, have been described. While it is widely accepted that V. destructor acts as vector of DWV, the question of whether the mite only functions as a mechanical vector or whether DWV can infect the mite thus using it as a biological vector is hotly debated, because in the literature data can be found that support both hypotheses. In order to settle this scientific dispute, we analyzed putatively DWV-infected mites with a newly established protocol for fluorescence-in situ-hybridization of mites and demonstrated DWV-specific signals inside mite cells. We provide compelling and direct evidence that DWV-B infects the intestinal epithelium and the salivary glands of V. destructor In contrast, no evidence for DWV-A infecting mite cells was found. Our data are key to understanding the pathobiology of DWV, the mite's role as a biological DWV vector and the quasispecies dynamics of this RNA virus when switching between insect and arachnid host species.IMPORTANCE Deformed wing virus (DWV) is a bee pathogenic, originally rather benign, single- and positive-stranded RNA virus. Only the vectorial transmission of this virus to honey bees by the ectoparasitic mite Varroa destructor leads to fatal or symptomatic infections of individuals, usually followed by collapse of the entire colony. Studies on whether the mite only acts as a mechanical virus vector or whether DWV can infect the mite and thus use it as a biological vector have led to disparate results. In our study using fluorescence-in situ-hybridization we provide compelling and direct evidence that at least the DWV-B variant infects the gut epithelium and the salivary glands of V. destructor Hence, the host range of DWV includes both, bees (Insecta) and mites (Arachnida). Our data contribute to a better understanding of the triangular relationship between honey bees, V. destructor and DWV and the evolution of virulence in this viral bee pathogen.

Deformed wing virus type a and b in managed honeybee colonies of Argentina

Apis mellifera is infected by more than 24 virus species worldwide, mainly positive-sense, single-stranded RNA viruses of the Dicistroviridae and Iflaviridae families. Among the viruses that infect honey bees, Deformed wing virus is the most prevalent and is present as three master variants DWV-A, B, and C. Given that the ectoparasitic mite Varroa destructor vectors these virus variants, recombination events between them are expected, and variants and their recombinants can co-exist in mites and honeybees at the same time. In this study, we detect, through RT-qPCR, the presence of DWV-A and B in the same samples of adult bees from colonies of Argentina. Total RNA was extracted from pools of ten adult bees from 45 apiaries distributed across the main beekeeping Provinces of Argentina (Buenos Aires, Santa Fe, Córdoba, Santiago del Estero, Río Negro, and Mendoza); then RT-qPCR reactions were performed to detect DWV-A and B, with specific primer pairs. After the amplifications, PCR products (204 and 660 bp amplicons for DWV-B, and ~250 bp for DWV-A) were purified and sequenced to verify that they corresponded to reported sequences, analyzing them using the Blast software. Of the 45 samples analyzed by RT-qPCR, over 90% were infected with DWV-A and 47% were also positive for DWV-B, where it was found in high prevalence specifically in colonies of A. mellifera of the Buenos Aires Province. Future studies will determine the impact of this type of the virus and its ability to recombine with the other DWV types in the apiaries of our country.

Possible Spillover of Pathogens between Bee Communities Foraging on the Same Floral Resource

Simple Summary

Floral resource availability is one of the keys to preserving the health of bee communities. However, flowers also present a risk of pathogen transmission, as infected pollinators could deposit pathogens while foraging, exposing other pollinators to infection via the consumption of contaminated nectar or pollen. Here, we studied, over time, the prevalence of seven viruses in bee communities that share the same small surface of floral resource in order to assess the risk of virus spillover. In total, 2057 bee specimens from 30 species were caught, identified and checked for the presence of viruses. Specimens from the Halictidae family were the dominant wild bees. The prevalence of viruses was quite high: at least one virus was detected in 78% of the samples, and co-infections were frequent. The genetic diversity of the viruses was also investigated to look for the possible association of geographic origin or host with shared ancestry.

Abstract

Viruses are known to contribute to bee population decline. Possible spillover is suspected from the co-occurrence of viruses in wild bees and honey bees. In order to study the risk of virus transmission between wild and managed bee species sharing the same floral resource, we tried to maximize the possible cross-infections using Phacelia tanacetifolia, which is highly attractive to honey bees and a broad range of wild bee species. Virus prevalence was compared over two years in Southern France. A total of 1137 wild bees from 29 wild bee species (based on COI barcoding) and 920 honey bees (Apis mellifera) were checked for the seven most common honey bee RNA viruses. Halictid bees were the most abundant. Co-infections were frequent, and Sacbrood virus (SBV), Black queen cell virus (BQCV), Acute bee paralysis virus (ABPV) and Israeli acute paralysis virus (IAPV) were widespread in the hymenopteran pollinator community. Conversely, Deformed wing virus (DWV) was detected at low levels in wild bees, whereas it was highly prevalent in honey bees (78.3% of the samples). Both wild bee and honey bee virus isolates were sequenced to look for possible host-specificity or geographical structuring. ABPV phylogeny suggested a specific cluster for Eucera bees, while isolates of DWV from bumble bees (Bombus spp.) clustered together with honey bee isolates, suggesting a possible spillover.

The Role of Pathogen Dynamics and Immune Gene Expression in the Survival of Feral Honey Bees

Studies of the ecoimmunology of feral organisms can provide valuable insight into how host–pathogen dynamics change as organisms transition from human-managed conditions back into the wild. Honey bees (Apis mellifera Linnaeus) offer an ideal system to investigate these questions as colonies of these social insects often escape management and establish in the wild. While managed honey bee colonies have low probability of survival in the absence of disease treatments, feral colonies commonly survive in the wild, where pathogen pressures are expected to be higher due to the absence of disease treatments. Here, we investigate the role of pathogen infections [Deformed wing virus (DWV), Black queen cell virus (BQCV), and Nosema ceranae] and immune gene expression (defensin-1, hymenoptaecin, pgrp-lc, pgrp-s2, argonaute-2, vago) in the survival of feral and managed honey bee colonies. We surveyed a total of 25 pairs of feral and managed colonies over a 2-year period (2017–2018), recorded overwintering survival, and measured pathogen levels and immune gene expression using quantitative polymerase chain reaction (qPCR). Our results showed that feral colonies had higher levels of DWV but it was variable over time compared to managed colonies. Higher pathogen levels were associated with increased immune gene expression, with feral colonies showing higher expression in five out of the six examined immune genes for at least one sampling period. Further analysis revealed that differential expression of the genes hymenoptaecin and vago increased the odds of overwintering survival in managed and feral colonies. Our results revealed that feral colonies express immune genes at higher levels in response to high pathogen burdens, providing evidence for the role of feralization in altering pathogen landscapes and host immune responses.

The Interplay Between Viruses and RNAi Pathways in Insects

As an overarching immune mechanism, RNA interference (RNAi) displays pathogen specificity and memory via different pathways. The small interfering RNA (siRNA) pathway is the primary antiviral defense mechanism against RNA viruses of insects and plays a lesser role in defense against DNA viruses. Reflecting the pivotal role of the siRNA pathway in virus selection, different virus families have independently evolved unique strategies to counter this host response, including protein-mediated, decoy RNA-based, and microRNA-based strategies. In this review, we outline the interplay between insect viruses and the different pathways of the RNAi antiviral response; describe practical application of these interactions for improved expression systems and for pest and disease management; and highlight research avenues for advancement of the field.

Deformed wing virus: using reverse genetics to tackle unanswered questions about the most important viral pathogen of honey bees

Deformed wing virus (DWV) is the most important viral pathogen of honey bees. It usually causes asymptomatic infections but, when vectored by the ectoparasitic mite Varroa destructor, it is responsible for the majority of overwintering colony losses globally. Although DWV was discovered four decades ago, research has been hampered by the absence of an in vitro cell culture system or the ability to culture pure stocks of the virus. The recent developments of reverse genetic systems for DWV go some way to addressing these limitations. They will allow the investigation of specific questions about strain variation, host tropism and pathogenesis to be answered, and are already being exploited to study tissue tropism and replication in Varroa and non-Apis pollinators. Three areas neatly illustrate the advances possible with reverse genetic approaches: (i) strain variation and recombination, in which reverse genetics has highlighted similarities rather than differences between virus strains; (ii) analysis of replication kinetics in both honey bees and Varroa, in studies that likely explain the near clonality of virus populations often reported; and (iii) pathogen spillover to non-Apis pollinators, using genetically tagged viruses to accurately monitor replication and infection.

OneHealth implications of infectious diseases of wild and managed bees

The OneHealth approach aims to further our understanding of the drivers of human, animal and environmental health, and, ultimately, to improve them by combining approaches and knowledge from medicine, biology and fields beyond. Wild and managed bees are essential pollinators of crops and wild flowers. Their health thus directly impacts on human and environmental health. At the same time, these bee species represent highly amenable and relevant model organisms for a OneHealth approach that aims to study fundamental epidemiological questions. In this review, we focus on how infectious diseases of wild and managed bees can be used as a OneHealth model system, informing fundamental questions on ecological immunology and disease transmission, while addressing how this knowledge can be used to tackle the issues facing pollinator health.

Varroa destructor mites vector and transmit pathogenic honey bee viruses acquired from an artificial diet

The ectoparasitic mite Varroa destructor is one of the most destructive pests of the honey bee (Apis mellifera) and the primary biotic cause of colony collapse in many regions of the world. These mites inflict physical injury on their honey bee hosts from feeding on host hemolymph and fat body cells/cellular components, and serve as the vector for deadly honey bee viruses, including Deformed wing virus (DWV) and the related Varroa destructor virus-1 (VDV-1) (i.e., DWV-like viruses). Studies focused on elucidating the dynamics of Varroa-mediated vectoring and transmission of DWV-like viruses may be confounded by viruses present in ingested host tissues or the mites themselves. Here we describe a system that includes an artificial diet free of insect tissue-derived components for maintaining Varroa mites for in vitro experimentation. Using this system, together with the novel engineered cDNA clone-derived genetically tagged VDV-1 and wild-type DWV, we demonstrated for the first time that Varroa mites provided an artificial diet supplemented with engineered viruses for 36 hours could acquire and transmit sufficient numbers of virus particles to establish an infection in virus-naïve hosts. While the in vitro system described herein provides for only up to five days of mite survival, precluding study of the long-term impacts of viruses on mite health, the system allows for extensive insights into the dynamics of Varroa-mediated vectoring and transmission of honey bee viruses.

RNAseq of Deformed Wing Virus and Other Honey Bee-Associated Viruses in Eight Insect Taxa with or without Varroa Infestation

The global spread of a parasitic mite (Varroa destructor) has resulted in Deformed wing virus (DWV), a previously rare pathogen, now dominating the viromes in honey bees and contributing to large-scale honey bee colony losses. DWV can be found in diverse insect taxa and has been implicated in spilling over from honey bees into associated ("apiary") and other ("non-apiary") insects. Here we generated next generation sequence data from 127 insect samples belonging to diverse taxa collected from Hawaiian islands with and without Varroa to identify whether the mite has indirectly affected the viral landscapes of key insect taxa across bees, wasps, flies and ants. Our data showed that, while Varroa was associated with a dramatic increase in abundance of (predominantly recombinant) DWV in honey bees (and no other honey bee-associated RNA virus), this change was not seen in any other taxa sampled. Honey bees share their environment with other insect populations and exist as a homogenous group, frequently sharing common viruses, albeit at low levels. Our data suggest that the threat of Varroa to increase viral load in an apiary does not automatically translate to an increase in virus load in other insects living in the wider community.