Diversity and Global Distribution of Viruses of the Western Honey Bee, Apis mellifera

Tetracycline-induced gut community dysbiosis and Israeli Acute Paralysis Virus infection synergistically negatively affect honeybees

Antibiotics are frequently employed to control bacterial diseases in honeybees, but their broad-spectrum action can disrupt the delicate balance of the gut microbiome, leading to dysbiosis. This imbalance in the gut microbiota of honeybees adversely affects their physiological health and weakens their resistance to pathogens, including viruses that significantly threaten honeybee health. In this study, we investigated whether tetracycline-induced gut microbiome dysbiosis promotes the replication of Israeli acute paralysis virus (IAPV), a key virus associated with colony losses and whether IAPV infection exacerbates gut microbiome dysbiosis. Our results demonstrated that tetracycline-induced gut microbiome dysbiosis increases the susceptibility of honeybees to IAPV infection. The viral titer in worker bees with antibiotic-induced gut microbiome dysbiosis prior to IAPV inoculation was significantly higher than in those merely inoculated with IAPV. Furthermore, we observed a synergistic effect between tetracycline and IAPV on the disruption of the honeybee gut microbiome balance. The progression of IAPV replication could, in turn, exacerbate antibiotic-induced gut microbiome dysbiosis in honeybees. Our research provides novel insights into the role of the gut microbiota in host-virus interactions, emphasizing the complex interplay between antibiotic use, gut microbiome health, and viral susceptibility in honeybees. We highlight the crucial role of a balanced gut microbiota in honey bees for their immune response against pathogens and emphasize the importance of careful, safe antibiotic use in beekeeping to protect these beneficial microbes.

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.

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.

Quantitative microbiome profiling of honey bee (Apis mellifera) guts is predictive of winter colony loss in northern Virginia (USA)

For the past 15 years, the proportion of honey bee hives that fail to survive winter has averaged ~ 30% in the United States. Winter hive loss has significant negative impacts on agriculture, the economy, and ecosystems. Compared to other factors, the role of honey bee gut microbial communities in driving winter hive loss has received little attention. We investigate the relationship between winter survival and honey bee gut microbiome composition of 168 honey bees from 23 hives, nine of which failed to survive through winter 2022. We found that there was a substantial difference in the abundance and community composition of honey bee gut microbiomes based on hive condition, i.e., winter survival or failure. The overall microbial abundance, as assessed using Quantitative Microbiome Profiling (QMP), was significantly greater in hives that survived winter 2022 than in those that failed, and the average overall abundance of each of ten bacterial genera was also greater in surviving hives. There were no significant differences in alpha diversity based on hive condition, but there was a highly significant difference in beta diversity. The bacterial genera Commensalibacter and Snodgrassella were positively associated with winter hive survival. Logistic regression and random forest machine learning models on pooled ASV counts for the genus data were highly predictive of winter outcome, although model performance decreased when samples from the location with no hive failures were excluded from analysis. As a whole, our results show that the abundance and community composition of honey bee gut microbiota is associated with winter hive loss, and can potentially be used as a diagnostic tool in evaluating hive health prior to the onset of winter. Future work on the functional characterization of the honey bee gut microbiome's role in winter survival is warranted.

Biotic and abiotic stresses on honeybee health

Honeybees are the most critical pollinators providing key ecosystem services that underpin crop production and sustainable agriculture. Amidst a backdrop of rapid global change, this eusocial insect encounters a succession of stressors during nesting, foraging, and pollination. Ectoparasitic mites, together with vectored viruses, have been recognized as central biotic threats to honeybee health, while the spread of invasive giant hornets and small hive beetles also increasingly threatens colonies worldwide. Cocktails of agrochemicals, including acaricides used for mite treatment, and other pollutants of the environment have been widely documented to affect bee health in various ways. Additionally, expanding urbanization, climate change, and agricultural intensification often result in the destruction or fragmentation of flower-rich bee habitats. The anthropogenic pressures exerted by beekeeping management practices affect the natural selection and evolution of honeybees, and colony translocations facilitate alien species invasion and disease transmission. In this review, the multiple biotic and abiotic threats and their interactions that potentially undermine bee colony health are discussed, while taking into consideration the sensitivity, large foraging area, dense network among related nestmates, and social behaviors of honeybees.

Current honey bee stressor investigations and mitigation methods in the United States and Canada

Honey bees are the most important managed insect pollinators in the US and Canadian crop systems. However, the annual mortality of colonies in the past 15 years has been consistently higher than historical records. Because they are eusocial generalist pollinators and amenable to management, honey bees provide a unique opportunity to investigate a wide range of questions at molecular, organismal, and ecological scales. Here, the American Association of Professional Apiculturists (AAPA) and the Canadian Association of Professional Apiculturists (CAPA) created 2 collections of articles featuring investigations on micro and macro aspects of honey bee health, sociobiology, and management showcasing new applied research from diverse groups studying honey bees (Apis mellifera) in the United States and Canada. Research presented in this special issue includes examinations of abiotic and biotic stressors of honey bees, and evaluations and introductions of various stress mitigation measures that may be valuable to both scientists and the beekeeping community. These investigations from throughout the United States and Canada showcase the wide breadth of current work done and point out areas that need further research.

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.

Discovery and characterization of novel DNA viruses in Apis mellifera: expanding the honey bee virome through metagenomic analysis

To date, many viruses have been discovered to infect honey bees. In this study, we used high-throughput sequencing to expand the known virome of the honey bee, Apis mellifera, by identifying several novel DNA viruses. While the majority of previously identified bee viruses are RNA, our study reveals nine new genomes from the Parvoviridae family, tentatively named Bee densoviruses 1 to 9. In addition, we characterized a large DNA virus, Apis mellifera filamentous-like virus (AmFLV), which shares limited protein identities with the known Apis mellifera filamentous virus. The complete sequence of AmFLV, obtained by a combination of laboratory techniques and bioinformatics, spans 152,678 bp. Linear dsDNA genome encodes for 112 proteins, of which 49 are annotated. Another large virus we discovered is Apis mellifera nudivirus, which belongs to a group of Alphanudivirus. The virus has a length of 129,467 bp and a circular dsDNA genome, and has 106 protein encoding genes. The virus contains most of the core genes of the family Nudiviridae. This research demonstrates the effectiveness of viral binning in identifying viruses in honey bee virology, showcasing its initial application in this field.IMPORTANCEHoney bees contribute significantly to food security by providing pollination services. Understanding the virome of honey bees is crucial for the health and conservation of bee populations and also for the stability of the ecosystems and economies for which they are indispensable. This study unveils previously unknown DNA viruses in the honey bee virome, expanding our knowledge of potential threats to bee health. The use of the viral binning approach we employed in this study offers a promising method to uncovering and understanding the vast viral diversity in these essential pollinators.

Prevalence and genome features of lake sinai virus isolated from Apis mellifera in the Republic of Korea

Lake Sinai Virus (LSV) is an emerging pathogen known to affect the honeybee (Apis mellifera). However, its prevalence and genomic characteristics in the Republic of Korea (ROK) remain unexplored. This study aimed to assess the prevalence of and analyze the LSVs by examining 266 honeybee samples from the ROK. Our findings revealed that LSV exhibited the highest infection rate among the pathogens observed in Korean apiaries, particularly during the reported period of severe winter loss (SWL) in A. mellifera apiaries in 2022. Three LSV genotypes- 2, 3, and 4 -were identified using RNA-dependent RNA polymerase gene analysis. Importantly, the infection rates of LSV2 (65.2%) and LSV3 (73.3%) were significantly higher in colonies experiencing SWL than in those experiencing normal winter loss (NWL) (p < 0.03). Furthermore, this study provides the first near-complete genome sequences of the Korean LSV2, LSV3, and LSV4 strains, comprising 5,759, 6,040, and 5,985 nt, respectively. Phylogenetic analysis based on these near-complete genome sequences demonstrated a close relationship between LSVs in the ROK and China. The high LSV infection rate in colonies experiencing a heightened mortality rate during winter suggests that this pathogen might contribute to SWL in ROK. Moreover, the genomic characteristic information on LSVs in this study holds immense potential for epidemiological information and the selection of specific genes suitable for preventing and treating LSV, including the promising utilization of RNA interference medicine in the future.

The effects of protein supplementation, fumagillin treatment, and colony management on the productivity and long-term survival of honey bee (Apis mellifera) colonies

In this study, we intensively measured the longitudinal productivity and survival of 362 commercially managed honey bee colonies in Canada, over a two-year period. A full factorial experimental design was used, whereby two treatments were repeated across apiaries situated in three distinct geographic regions: Northern Alberta, Southern Alberta and Prince Edward Island, each having unique bee management strategies. In the protein supplemented treatment, colonies were continuously provided a commercial protein supplement containing 25% w/w pollen, in addition to any feed normally provided by beekeepers in that region. In the fumagillin treatment, colonies were treated with the label dose of Fumagilin-B® each year during the fall. Neither treatment provided consistent benefits across all sites and dates. Fumagillin was associated with a large increase in honey production only at the Northern Alberta site, while protein supplementation produced an early season increase in brood production only at the Southern Alberta site. The protein supplement provided no long-lasting benefit at any site and was also associated with an increased risk of death and decreased colony size later in the study. Differences in colony survival and productivity among regions, and among colonies within beekeeping operations, were far larger than the effects of either treatment, suggesting that returns from extra feed supplements and fumagillin were highly contextually dependent. We conclude that use of fumagillin is safe and sometimes beneficial, but that beekeepers should only consider excess protein supplementation when natural forage is limiting.

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.

Ecological and social factors influence interspecific pathogens occurrence among bees

The interspecific transmission of pathogens can occur frequently in the environment. Among wild bees, the main spillover cases are caused by pathogens associated with Apis mellifera, whose colonies can act as reservoirs. Due to the limited availability of data in Italy, it is challenging to accurately assess the impact and implications of this phenomenon on the wild bee populations. In this study, a total of 3372 bees were sampled from 11 Italian regions within the BeeNet project, evaluating the prevalence and the abundance of the major honey bee pathogens (DWV, BQCV, ABPV, CBPV, KBV, Nosema ceranae, Ascosphaera apis, Crithidia mellificae, Lotmaria passim, Crithidia bombi). The 68.4% of samples were positive for at least one pathogen. DWV, BQCV, N. ceranae and CBPV showed the highest prevalence and abundance values, confirming them as the most prevalent pathogens spread in the environment. For these pathogens, Andrena, Bombus, Eucera and Seladonia showed the highest mean prevalence and abundance values. Generally, time trends showed a prevalence and abundance decrease from April to July. In order to predict the risk of infection among wild bees, statistical models were developed. A low influence of apiary density on pathogen occurrence was observed, while meteorological conditions and agricultural management showed a greater impact on pathogen persistence in the environment. Social and biological traits of wild bees also contributed to defining a higher risk of infection for bivoltine, communal, mining and oligolectic bees. Out of all the samples tested, 40.5% were co-infected with two or more pathogens. In some cases, individuals were simultaneously infected with up to five different pathogens. It is essential to increase knowledge about the transmission of pathogens among wild bees to understand dynamics, impact and effects on pollinator populations. Implementing concrete plans for the conservation of wild bee species is important to ensure the health of wild and human-managed bees within a One-Health perspective.

Distribution of infectious and parasitic agents among three sentinel bee species across European agricultural landscapes

Infectious and parasitic agents (IPAs) and their associated diseases are major environmental stressors that jeopardize bee health, both alone and in interaction with other stressors. Their impact on pollinator communities can be assessed by studying multiple sentinel bee species. Here, we analysed the field exposure of three sentinel managed bee species (Apis mellifera, Bombus terrestris and Osmia bicornis) to 11 IPAs (six RNA viruses, two bacteria, three microsporidia). The sentinel bees were deployed at 128 sites in eight European countries adjacent to either oilseed rape fields or apple orchards during crop bloom. Adult bees of each species were sampled before their placement and after crop bloom. The IPAs were detected and quantified using a harmonised, high-throughput and semi-automatized qPCR workflow. We describe differences among bee species in IPA profiles (richness, diversity, detection frequencies, loads and their change upon field exposure, and exposure risk), with no clear patterns related to the country or focal crop. Our results suggest that the most frequent IPAs in adult bees are more appropriate for assessing the bees' IPA exposure risk. We also report positive correlations of IPA loads supporting the potential IPA transmission among sentinels, suggesting careful consideration should be taken when introducing managed pollinators in ecologically sensitive environments.

Seasonal trends of the ABPV, KBV, and IAPV complex in Italian managed honey bee (Apis mellifera L.) colonies

Acute bee paralysis virus (ABPV), Kashmir bee virus (KBV), and Israeli acute paralysis virus (IAPV) usually persist as covert infections in honey bee colonies. They can cause rapid bee mortality in cases of severe infection, often associated with high Varroa destructor infestation, by which they are transmitted. In various countries, these viruses have been associated with colony collapse. Despite their potential danger, these viruses are often disregarded, and little information is available on their occurrence in many countries, including Italy. In 2021, 370 apiaries representing all of the Italian regions were investigated in four different months (June, September, November, and March) for the presence of ABPV, KBV, and IAPV. IAPV was not found in any of the apiaries investigated, whereas 16.45% and 0.67% of the samples tested positive for ABPV and KBV, respectively. Most ABPV cases occurred in late summer-autumn in both northern and southern regions. We observed a scattered pattern of KBV-positive colonies that did not allow any seasonal or regional trends to be discerned. Differences observed among regions and months were potentially related to the dynamics of varroa infestation, viral genetic variations, and different climatic conditions resulting in variations in bee behaviour. This study improves our understanding of the circulation of bee viruses and will contribute to better disease prevention and preservation of bee health.

WIV, a protein domain found in a wide number of arthropod viruses, which probably facilitates infection

The most powerful approach to detect distant homologues of a protein is based on structure prediction and comparison. Yet this approach is still inapplicable to many viral proteins. Therefore, we applied a powerful sequence-based procedure to identify distant homologues of viral proteins. It relies on three principles: (1) traces of sequence similarity can persist beyond the significance cutoff of homology detection programmes; (2) candidate homologues can be identified among proteins with weak sequence similarity to the query by using 'contextual' information, e.g. taxonomy or type of host infected; (3) these candidate homologues can be validated using highly sensitive profile-profile comparison. As a test case, this approach was applied to a protein without known homologues, encoded by ORF4 of Lake Sinai viruses (which infect bees). We discovered that the ORF4 protein contains a domain that has homologues in proteins from >20 taxa of viruses infecting arthropods. We called this domain 'widespread, intriguing, versatile' (WIV), because it is found in proteins with a wide variety of functions and within varied domain contexts. For example, WIV is found in the NSs protein of tospoviruses, a global threat to food security, which infect plants as well as their arthropod vectors; in the RNA2 ORF1-encoded protein of chronic bee paralysis virus, a widespread virus of bees; and in various proteins of cypoviruses, which infect the silkworm Bombyx mori. Structural modelling with AlphaFold indicated that the WIV domain has a previously unknown fold, and bibliographical evidence suggests that it facilitates infection of arthropods.

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.

Quantitative trait loci mapping for survival of virus infection and virus levels in honey bees

Israeli acute paralysis virus (IAPV) is a highly virulent, Varroa-vectored virus that is of global concern for honey bee health. Little is known about the genetic basis of honey bees to withstand infection with IAPV or other viruses. We set up and analyzed a backcross between preselected honey bee colonies of low and high IAPV susceptibility to identify quantitative trait loci (QTL) associated with IAPV susceptibility. Experimentally inoculated adult worker bees were surveyed for survival and selectively sampled for QTL analysis based on SNPs identified by whole-genome resequencing and composite interval mapping. Additionally, natural titers of other viruses were quantified in the abdomen of these workers via qPCR and also used for QTL mapping. In addition to the full dataset, we analyzed distinct subpopulations of susceptible and non-susceptible workers separately. These subpopulations are distinguished by a single, suggestive QTL on chromosome 6, but we identified numerous other QTL for different abdominal virus titers, particularly in the subpopulation that was not susceptible to IAPV. The pronounced QTL differences between the susceptible and non-susceptible subpopulations indicate either an interaction between IAPV infection and the bees' interaction with other viruses or heterogeneity among workers of a single cohort that manifests itself as IAPV susceptibility and results in distinct subgroups that differ in their interaction with other viruses. Furthermore, our results indicate that low susceptibility of honey bees to viruses can be caused by both, virus tolerance and virus resistance. QTL were partially overlapping among different viruses, indicating a mixture of shared and specific processes that control viruses. Some functional candidate genes are located in the QTL intervals, but their genomic co-localization with numerous genes of unknown function delegates any definite characterization of the underlying molecular mechanisms to future studies.

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.

A preliminary survey reveals that common viruses are found at low titers in a wild population of honey bees (Apis mellifera)

A major threat to honey bee (Apis mellifera Linnaeus, Hymenoptera: Apidae) health continues to be parasitism by the mite Varroa destructor, which has been linked to high colony losses worldwide. Besides feeding on developing and adult bees, Varroa is also a prolific vector of honey bee-associated viruses. Because they live in unmanaged conditions, wild honey bee colonies are not treated against Varroa, which has enabled the natural selection of more mite-tolerant bees. To date, few studies have explored the prevalence of viruses in unmanaged colonies. The Welder Wildlife Refuge (WWR) in Texas is a unique site to study the viral landscape of unmanaged honey bees in the United States. The goals of this study were to identify and quantify viruses in wild colonies at the WWR, to examine changes in the prevalence of viruses in these colonies over time, and to compare the presence and titers of viruses between wild colonies at the WWR and those from the nearest managed apiary. We collected bees from colonies at the WWR in 2013, 2016, and 2021, and analyzed selected viruses for their presence and titers via quantitative polymerase chain reaction. In 2021, we also sampled bees from the nearest managed apiary for comparison. We found low average virus titers in all wild colonies sampled, and no difference in virus titers between colonies at the WWR and those from the managed apiary. Our study indicates that virus titers in wild colonies at the WWR are similar to those found in nearby colonies, and that these titers fluctuate over time.

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.

Honey bee retinue workers respond similarly to queens despite seasonal differences in Queen Mandibular Pheromone (QMP) signaling

Honey bee colonies maintain viable queens in part through communication with Queen Mandibular Pheromone (QMP), a mixture that signals the queen's presence and reproductive quality to workers. In turn, workers are thought to provide retinue queen care or replace queens partially based on QMP profiles. We examined the effects of seasonal dearth (overwintering in a warm subtropical location) on queen-worker interactions. Retinue worker responses to continuously ovipositing queens were considered in view of QMP signaling and queen reproductive quality. QMP signaling was estimated from QMP residues recovered from nest worker bodies, which is the primary mode of QMP transfer from the queen to the colony at large. QMP residues varied seasonally but not at all with queen reproductive quality (spermatheca sperm storage, ovary protein and lipid contents). 9-HDA and 9-ODA were lower in January than other months. HOB decreased from July to January, while HVA, a component associated with mated queens, increased sharply in January. Despite these seasonal signaling differences, retinue workers attended queens at similar levels through the months. In terms of reproductive quality, queens did not differ over the months in matedness (spermatheca sperm storage) or physiological age (protein carbonyl content), but varied in nutrient allocation to reproductive and non-reproductive tissues. Queen ovaries contained more protein in September than in November, and more lipid in July and September than in November and January. Queen fat bodies had more protein in July than September or November, but less lipid in July and September than November or January. Retinue worker responses did not vary with seasonal QMP changes, but reflected overall continuous brood rearing efforts and queen matedness throughout the year. The absence of seasonal differences in worker responses to QMP should be considered in the broader context of continuous reproductive efforts in warm subtropical colonies.

Effects of dialkoxybenzenes against Varroa destructor and identification of 1-allyloxy-4-propoxybenzene as a promising acaricide candidate

The honey bee is responsible for pollination of a large proportion of crop plants, but the health of honey bee populations has been challenged by the parasitic mite Varroa destructor. Mite infestation is the main cause of colony losses during the winter months, which causes significant economic challenges in apiculture. Treatments have been developed to control the spread of varroa. However, many of these treatments are no longer effective due to acaricide resistance. In a search of varroa-active compounds, we tested the effect of dialkoxybenzenes on the mite. A structure-activity relationship revealed that 1-allyloxy-4-propoxybenzene is most active of a series of dialkoxybenzenes tested. We found that three compounds (1-allyloxy-4-propoxybenzene, 1,4-diallyloxybenzene and 1,4-dipropoxybenzene) cause paralysis and death of adult varroa mites, whereas the previously discovered compound, 1,3-diethoxybenzene, which alters host choice of adult mites in certain conditions, did not cause paralysis. Since paralysis can be caused by inhibition of acetylcholinesterase (AChE), a ubiquitous enzyme in the nervous system of animals, we tested dialkoxybenzenes on human, honey bee and varroa AChE. These tests revealed that 1-allyloxy-4-propoxybenzene had no effects on AChE, which leads us to conclude that 1-allyloxy-4-propoxybenzene does not exert its paralytic effect on mites through AChE. In addition to paralysis, the most active compounds affected the ability of the mites to find and remain at the abdomen of host bees provided during assays. A test of 1-allyloxy-4-propoxybenzene in the field, during the autumn of 2019 in two locations, showed that this compound has promise in the treatment of varroa infestations.

A systematic review of honey bee (Apis mellifera, Linnaeus, 1758) infections and available treatment options

BACKGROUND

Honey bees and honeycomb bees are very valuable for wild flowering plants and economically important crops due to their role as pollinators. However, these insects confront many disease threats (viruses, parasites, bacteria and fungi) and large pesticide concentrations in the environment. Varroa destructor is the most prevalent disease that has had the most negative effects on the fitness and survival of different honey bees (Apis mellifera and A. cerana). Moreover, honey bees are social insects and this ectoparasite can be easily transmitted within and across bee colonies.

OBJECTIVE

This review aims to provide a survey of the diversity and distribution of important bee infections and possible management and treatment options, so that honey bee colony health can be maintained.

METHODS

We used PRISMA guidelines throughout article selection, published between January 1960 and December 2020. PubMed, Google Scholar, Scopus, Cochrane Library, Web of Science and Ovid databases were searched.

RESULTS

We have collected 132 articles and retained 106 articles for this study. The data obtained revealed that V. destructor and Nosema spp. were found to be the major pathogens of honey bees worldwide. The impact of these infections can result in the incapacity of forager bees to fly, disorientation, paralysis, and death of many individuals in the colony. We find that both hygienic and chemical pest management strategies must be implemented to prevent, reduce the parasite loads and transmission of pathogens. The use of an effective miticide (fluvalinate-tau, coumaphos and amitraz) now seems to be an essential and common practice required to minimise the impact of Varroa mites and other pathogens on bee colonies. New, alternative biofriendly control methods, are on the rise, and could be critical for maintaining honey bee hive health and improving honey productivity.

CONCLUSIONS

We suggest that critical health control methods be adopted globally and that an international monitoring system be implemented to determine honey bee colony safety, regularly identify parasite prevalence, as well as potential risk factors, so that the impact of pathogens on bee health can be recognised and quantified on a global scale.

Transcriptomic Responses Underlying the High Virulence of Black Queen Cell Virus and Sacbrood Virus following a Change in Their Mode of Transmission in Honey Bees (Apis mellifera)

BACKGROUND

Over the last two decades, honey bees (Apis mellifera) have suffered high rates of colony losses that have been attributed to a variety of factors, chief among which are viral pathogens, such as deformed wing virus (DWV), whose virulence has increased because of vector-based transmission by the invasive, ectoparasitic varroa mite (Varroa destructor). A shift in the experimental mode of transmission of the black queen cell virus (BQCV) and sacbrood virus (SBV) from fecal/food-oral (direct horizontal) to vector-mediated (indirect horizontal) transmission also results in high virulence and viral titers in pupal and adult honey bees. Agricultural pesticides represent another factor that acts independently or in interaction with pathogens, and they are also thought to cause colony loss. Understanding the molecular mechanisms underlying the higher virulence following a vector-based mode of transmission provides deeper insight into honey bee colony losses, as does determining whether or not host-pathogen interactions are modulated by exposure to pesticides.

METHODS

Through an experimental design with controlled laboratory, we investigated the effects of the modes of transmission of BQCV and SBV (feeding vs. vector-mediated via injection) alone or in combination with chronic exposure to sublethal and field-realistic concentrations of flupyradifurone (FPF), a novel agricultural insecticide, on honey bee survival and transcription responses by using high-throughput RNA sequencing (RNA-seq) analysis.

RESULTS

Co-exposure to viruses via feeding (VF) or injection (VI) and FPF insecticide had no statistically significant interactive effect on their survival compared to, respectively, VF or VI treatments alone. Transcriptomic analysis revealed a distinct difference in the gene expression profiles of bees inoculated with viruses via injection (VI) and exposed to FPF insecticide (VI+FPF). The number of differentially expressed genes (DEGs) at log2 (fold-change) > 2.0 in VI bees (136 genes) or/and VI+FPF insecticide (282 genes) was very high compared to that of VF bees (8 genes) or the VF+FPF insecticide treatment (15 genes). Of these DEGs, the expression in VI and VI+FPF bees of some immune-related genes, such as those for antimicrobial peptides, Ago2, and Dicer, was induced. In short, several genes encoding odorant binding proteins, chemosensory proteins, odor receptors, honey bee venom peptides, and vitellogenin were downregulated in VI and VI+FPF bees.

CONCLUSIONS

Given the importance of these suppressed genes in honey bees' innate immunity, eicosanoid biosynthesis, and olfactory associative function, their inhibition because of the change in the mode of infection with BQCV and SBV to vector-mediated transmission (injection into haemocoel) could explain the high virulence observed in these viruses when they were experimentally injected into hosts. These changes may help explain why other viruses, such as DWV, represent such a threat to colony survival when transmitted by varroa mites.

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.

Area Wide Monitoring of Plant and Honey Bee (Apis mellifera) Viruses in Blueberry (Vaccinium corymbosum) Agroecosystems Facilitated by Honey Bee Pollination

Healthy agroecosystems are dependent on a complex web of factors and inter-species interactions. Flowers are hubs for pathogen transmission, including the horizontal or vertical transmission of plant-viruses and the horizontal transmission of bee-viruses. Pollination by the European honey bee (Apis mellifera) is critical for industrial fruit production, but bees can also vector viruses and other pathogens between individuals. Here, we utilized commercial honey bee pollination services in blueberry (Vaccinium corymbosum) farms for a metagenomics-based bee and plant virus monitoring system. Following RNA sequencing, viruses were identified by mapping reads to a reference sequence database through the bioinformatics portal Virtool. In total, 29 unique plant viral species were found at two blueberry farms in British Columbia (BC). Nine viruses were identified at one site in Ontario (ON), five of which were not identified in BC. Ilarviruses blueberry shock virus (BlShV) and prune dwarf virus (PDV) were the most frequently detected viruses in BC but absent in ON, while nepoviruses tomato ringspot virus and tobacco ringspot virus were common in ON but absent in BC. BlShV coat protein (CP) nucleotide sequences were nearly identical in all samples, while PDV CP sequences were more diverse, suggesting multiple strains of PDV circulating at this site. Ten bee-infecting viruses were identified, with black queen cell virus frequently detected in ON and BC. Area-wide bee-mediated pathogen monitoring can provide new insights into the diversity of viruses present in, and the health of, bee-pollination ecosystems. This approach can be limited by a short sampling season, biased towards pollen-transmitted viruses, and the plant material collected by bees can be very diverse. This can obscure the origin of some viruses, but bee-mediated virus monitoring can be an effective preliminary monitoring approach.

Expression of Molecular Markers of Resilience against Varroa destructor and Bee Viruses in Ethiopian Honey Bees (Apis mellifera simensis) Focussing on Olfactory Sensing and the RNA Interference Machinery

Varroa destructor mites and the viruses it vectors are two major factors leading to high losses of honey bees (Apis mellifera) colonies worldwide. However, honey bees in some African countries show resilience to varroa infestation and/or virus infections, although little is known about the mechanisms underlying this resilience. In this study, we investigated the expression profiles of some key molecular markers involved in olfactory sensing and RNA interference, as these processes may contribute to the bees' resilience to varroa infestation and virus infection, respectively. We found significantly higher gene expression of the odorant binding protein, OBP14, in the antennae of Ethiopian bees compared to Belgian bees. This result suggests the potential of OBP14 as a molecular marker of resilience to mite infestation. Scanning electron microscopy showed no significant differences in the antennal sensilla occurrence and distribution, suggesting that resilience arises from molecular processes rather than morphological adaptations. In addition, seven RNAi genes were upregulated in the Ethiopian honey bees and three of them-Dicer-Drosha, Argonaute 2, and TRBP2-were positively correlated with the viral load. We can conclude that the antiviral immune response was triggered when bees were experiencing severe viral infection and that this might contribute to the bees' resilience to viruses.

The final frontier: ecological and evolutionary dynamics of a global parasite invasion

Studying rapid biological changes accompanying the introduction of alien organisms into native ecosystems can provide insights into fundamental ecological and evolutionary theory. While powerful, this quasi-experimental approach is difficult to implement because the timing of invasions and their consequences are hard to predict, meaning that baseline pre-invasion data are often missing. Exceptionally, the eventual arrival of Varroa destructor (hereafter Varroa) in Australia has been predicted for decades. Varroa is a major driver of honeybee declines worldwide, particularly as vectors of diverse RNA viruses. The detection of Varroa in 2022 at over a hundred sites poses a risk of further spread across the continent. At the same time, careful study of Varroa's spread, if it does become established, can provide a wealth of information that can fill knowledge gaps about its effects worldwide. This includes how Varroa affects honeybee populations and pollination. Even more generally, Varroa invasion can serve as a model for evolution, virology and ecological interactions between the parasite, the host and other organisms.

First molecular detection of the presence of honey bee viruses in insects, Varroa destructor mites, and pollinated plants in an isolated region of Armenia

Background and Aim

Recently, viral diseases of honey bees (Apis mellifera) have presented an increasing threat to beekeeping. This study aimed to examine the presence of honey bee viruses in Apis and non-Apis bee species, the mite Varroa destructor, and pollinated plants in Armenia.

Materials and Methods

Sampling was performed in Tavush Province, in the northeast of the Republic of Armenia, from August to November 2019. Overall, 200 A. mellifera bees, 50 V. destructor mites, and 20 wasps were collected (corresponding to three bees, five mites, and 2-11 wasps in each investigated sample) and homogenized for RNA isolation and detection of viruses. Ten pollinated plants were taken from each plant, and 2 g of each sample was used for homogenization. In each investigated case Apis mellifera, Varroa destructor, Vespula germanica and plants received percentages of the virus presence.

Results

Six important honey bee viruses (acute bee paralysis virus [ABPV], deformed wing virus [DWV], A. mellifera norovirus [ANV], Lake Sinai virus-2 [LSV-2], Big Sioux River virus [BSRV], and A. mellifera filamentous virus [AmFV]) were detected in samples by polymerase chain reaction. Our results showed that DWV, ANV, and ABPV were the most common viruses in honey bees. All viruses were detected in wasps, but LSV-2 and ANV were present in almost all samples.

Conclusion

Our results showed that almost all viruses were present in V. destructor. Although ANV is very common in honey bees, it did not appear in any mite samples. Our study indicates that viruses typically associated with honey bees were also actively infecting wasps. Our data suggest that the survival of viruses in plants can be an important source of seasonal transmission of viruses to bees. In addition, pollinated plants can potentially serve as reservoirs for honey bee viruses.

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

A Preliminary Study on “Personalised Treatment” against Varroa destructor Infestations in Honey Bee (Apis mellifera) Colonies

The requirement for the protection of bee colonies against Varroa destructor invasions has been noted by many breeders and is included as an aspect of the development of beekeeping. This research aimed to check the effect of the development of a colony exposed to laying eggs (brood surface) by queen bees with similar chemical potential (sisters) on the effect of a preparation combating V. destructor depending on the number of mites found in a given colony. We chose this as a standard model of conduct that treats each bee colony as one organism subjected to individual parasite control. For this purpose, we created a bee colony with a mother-of-one breeding line and fertilised drones from one colony. Infection with V. destructor occurred naturally and uncontrollably. Without interfering with the colony’s development, the frame insulator helped each colony’s brood (mothers’ reproductive potential) and the initial and final individuals from the mites themselves. The study was carried out in four species (two control species and two species with up to 20 and over 21 mites, respectively). Treatments with amitraz to combat damage were divided into four treatment subgroups: two treatments every four days or four treatments every two days. We observed the number of individuals that were protected in all subgroups in the average brood area. The reproductive potential of the sisters’ mothers did not change after the treatments with amitraz, which indicated that amitraz did not affect the delegation of egg laying. The invasion rate was also tracked relative to the control group, which allowed us to conclude that a two-time treatment with amitraz reduced the frequency of mites and a four-time treatment checked the effectiveness. Tailoring the control of V. destructor in bee colonies may be an effective measure in the fight against this parasite.

Virus transmission via honey bee prey and potential impact on cocoon-building in labyrinth spiders (Agelena labyrinthica)

Interspecies transmission of RNA viruses is a major concern for human and animal health. However, host-range, transmission routes and especially the possible impact of these viruses on alternative hosts are often poorly understood. Here, we investigated the role of the labyrinth spider, Agelena labyrinthica, as a potential alternative host of viruses commonly known from western honey bees, Apis mellifera. Field-collected spiders were screened for Acute bee paralysis virus (ABPV), Black queen cell virus, Chronic bee paralysis virus, Deformed wing virus type A and B (DWV-B), Israeli acute paralysis virus, Lake Sinai virus and Sacbrood virus. In a laboratory experiment, labyrinth spiders were fed with ABPV and DWV-B infected honey bees or virus free control food. Our results show that natural infections of A. labyrinthica with these viruses are common in the field, as 62.5% of the samples were positive for at least one virus, supporting their wide host range. For DWV-B, the laboratory data indicate that foodborne transmission occurs and that high virus titres may reduce cocoon building, which would be the first report of clinical symptoms of DWV in Araneae. Since cocoons are tokens of fitness, virus transmission from honey bees might affect spider populations, which would constitute a concern for nature conservation.

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.

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.

Nationwide genomic surveillance reveals the prevalence and evolution of honeybee viruses in China

BACKGROUND

The economic and environmental value of honeybees has been severely challenged in recent years by the collapse of their colonies worldwide, often caused by outbreaks of infectious diseases. However, our understanding of the diversity, prevalence, and transmission of honeybee viruses is largely obscure due to a lack of large-scale and longitudinal genomic surveillance on a global scale.

RESULTS

We report the meta-transcriptomic sequencing of nearly 2000 samples of the two most important economic and widely maintained honeybee species, as well as an associated ectoparasite mite, collected across China during 2016-2019. We document the natural diversity and evolution of honeybee viruses in China, providing evidence that multiple viruses commonly co-circulate within individual bee colonies. We also expanded the genomic data for 12 important honeybee viruses and revealed novel genetic variants and lineages associated with China. We identified more than 23 novel viruses from the honeybee and mite viromes, with some exhibiting ongoing replication in their respective hosts. Together, these data provide additional support to the idea that mites are an important reservoir and spill-over host for honeybee viruses.

CONCLUSIONS

Our data show that honeybee viruses are more widespread, prevalent, and genetically diverse than previously realized. The information provided is important in mitigating viral infectious diseases in honeybees, in turn helping to maintain sustainable productive agriculture on a global scale. Video Abstract.

Trans-generational immune priming against American Foulbrood does not affect the performance of honeybee colonies

Honeybees are major pollinators for our food crops, but at the same time they face many stressors all over the world. One of the major threats to honeybee health are bacterial diseases, the most severe of which is the American Foulbrood (AFB). Recently a trans-generational vaccination approach against AFB has been proposed, showing strong potential in protecting the colonies from AFB outbreaks. Yet, what remains unstudied is whether the priming of the colony has any undesired side-effects. It is widely accepted that immune function is often a trade-off against other life-history traits, hence immune priming could have an effect on the colony performance. In this experiment we set up 48 hives, half of them with primed queens and half of them as controls. The hives were placed in six apiaries, located as pair of apiaries in three regions. Through a 2-year study we monitored the hives and measured their health and performance. We measured hive weight and frame contents such as brood amount, worker numbers, and honey yield. We studied the prevalence of the most common honeybee pathogens in the hives and expression of relevant immune genes in the offspring at larval stage. No effect of trans-generational immune priming on any of the hive parameters was found. Instead, we did find other factors contributing on various hive performance parameters. Interestingly not only time but also the region, although only 10 km apart from each other, had an effect on the performance and health of the colonies, suggesting that the local environment plays an important role in hive performance. Our results suggest that exploiting the trans-generational priming could serve as a safe tool in fighting the AFB in apiaries.

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.

A Reciprocal Transplant Experiment Confirmed Mite-Resistance in a Honey Bee Population from Uruguay

Simple Summary

In Uruguay, as in many countries around the world, the Varroa destructor mite is the main biotic threat to honey bees (Apis mellifera). Most beekeepers regularly apply acaricides to their colonies to have good honey harvests and avoid large losses, with the exception of beekeepers in the east of the country where bees coexist with V. destructor without suffering significant damage. To unravel the different A. mellifera–V. destructor relationships found in the country, a reciprocal transplant experiment was performed between the mite-resistant bee colonies and the mite-susceptible bee colonies from the east and the west of the country, respectively. The differences between the two groups of bees in the control of V. destructor were maintained in the two environments. No mite-susceptible colonies survived the winter. The behavioral resistance of bees (hygienic behavior) and reproductive aspects of V. destructor (phoretic mites/reproductive mites and mites in drone cells/mites in worker cells ratio) could explain the results obtained.

Abstract

In the past few years there has been an increasing interest for the study of honey bee populations that are naturally resistant to the ectoparasitic mite Varroa destructor, aiming to identify the mechanisms that allow the bees to limit the reproduction of the mite. In eastern Uruguay there are still bees resistant to mites that survive without acaricides. In order to determine if the differential resistance to V. destructor was maintained in other environments, a reciprocal transplant experiment was performed between the mite-resistant bee colonies and the mite-susceptible bee colonies from the east and the west of the country, respectively, infesting bees with local mites. In both regions, the mite-resistant colonies expressed a higher hygienic behavior and presented a higher phoretic mites/reproductive mites and mites in drone cells/mites in worker cells ratio than the mite-susceptible colonies. All the mite-susceptible colonies died during fall–winter, while a considerable number of mite-resistant colonies survived until spring, especially in the east of the country. This study shows that the bees in the east of the country maintain in good measure the resistance to V. destructor in other regions and leaves open the possibility that the mites of the two populations have biases in the reproductive behavior.

Acute and chronic viruses mediated by an ectoparasite targeting different developmental stages of honeybee (Apis mellifera and Apis cerana) brood

The health of the western honeybee, Apis mellifera, the most crucial pollinator, has been challenged globally over the past decades. An ectoparasitic mite, Varroa destructor, together with the viruses it vectored, is generally regarded as the vital pathogenic agent. Although the poor health status of A. mellifera compared to its eastern counterpart, Apis cerana, has been broadly identified, the underlying mechanism remains poorly understood and comparison between susceptible and resistant hosts will potentially ameliorate this predicament. Here, we investigated the impacts of two widespread viruses—deformed wing virus type A (DWV-A) and Israeli acute paralysis virus (IAPV), mediated by V. destructor mite, on the capped developing honeybee brood, in the absence of adult workers, of A. mellifera and A. cerana, with positive and negative controls. Our results demonstrated that the endogenous viruses imposed limited damage on the hosts even if the brood was wounded. In contrast, the exogenous viruses introduced by ectoparasites triggered variable mortality of the infested brood between host species. Intriguingly, death causes of both honeybee species presented a similar trend: the acute IAPV generally causes morbidity and mortality of late larvae, while the chronic DWV-A typically leads to brood mortality during and after pupation. Notably, the susceptible immature A. cerana individuals, supported by higher observed mortality and a lower virus tolerance, serve the interests of the colony and foster the overall survival of a resistant honeybee superorganism. These results improve our understanding of the interactions between viruses carried by ectoparasites and their developing hosts, and the novel insight of weak individuals fostering strong colonies may promote breeding efforts to mitigate the indefensible colony losses globally.

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.

Sacbrood Virus: A Growing Threat to Honeybees and Wild Pollinators

Sacbrood virus (SBV) is one of the many viruses that infect both the Western honeybee (Apis mellifera) and the Eastern honeybee (Apis cerana). Recently, the interspecies transmission of SBV has been discovered, especially among wild pollinators. This newly discovered evolutionary occurrence regarding SBV indicates a much wider host range than previously believed, causing further concern about the future sustainability of agriculture and the resilience of ecosystems. Over the past few decades, vast numbers of studies have been undertaken concerning SBV infection in honeybees, and remarkable progress has been made in our understanding of the epidemiology, pathogenesis, transmission, and manifestations of SBV infection in honeybees and other pollinators. Meanwhile, some methods, including Chinese medicine, have been established to control and prevent sacbrood disease in A. cerana in Asian countries. In this review, we summarize the existing knowledge of SBV and address the gaps in the knowledge within the existing literature in the hope of providing future directions for the research and development of management strategies for controlling the spread of this deadly disease.

The Epidemiological Situation of the Managed Honey Bee (Apis mellifera) Colonies in the Italian Region Emilia-Romagna

The recent decades witnessed the collapse of honey bee colonies at a global level. The major drivers of this collapse include both individual and synergic pathogen actions, threatening the colonies’ survival. The need to define the epidemiological pattern of the pathogens that are involved has led to the establishment of monitoring programs in many countries, Italy included. In this framework, the health status of managed honey bees in the Emilia–Romagna region (northern Italy) was assessed, throughout the year 2021, on workers from 31 apiaries to investigate the presence of major known and emerging honey bee pathogens. The prevalence and abundance of DWV, KBV, ABPV, CBPV, Nosema ceranae, and trypanosomatids (Lotmaria passim, Crithidia mellificae, Crithidia bombi) were assessed by molecular methods. The most prevalent pathogen was DWV, followed by CBPV and N. ceranae. Trypanosomatids were not found in any of the samples. Pathogens had different peaks in abundance over the months, showing seasonal trends that were related to the dynamics of both bee colonies and Varroa destructor infestation. For some of the pathogens, a weak but significant correlation was observed between abundance and geographical longitude. The information obtained in this study increases our understanding of the epidemiological situation of bee colonies in Emilia–Romagna and helps us to implement better disease prevention and improved territorial management of honey bee health.

Total RNA sequencing of Phlebotomus chinensis, a neglected vector in China, simultaneously revealed viral, bacterial, and eukaryotic microbes that are potentially pathogenic to humans

Phlebotomus chinensis sandfly is a neglected insect vector in China that is well-known for carrying Leishmania. Recent studies have expanded its pathogen repertoire with two novel arthropod-borne phleboviruses capable of infecting humans and animals. Despite these discoveries, our knowledge of the general pathogen diversity and overall microbiome composition of this vector species is still very limited. Here we carried out a meta-transcriptomics analysis that revealed the actively replicating/transcribing RNA viruses, DNA viruses, bacteria, and eukaryotic microbes, namely, the “total microbiome”, of several sandfly populations in China. Strikingly, “microbiome” made up 1.8% of total non-ribosomal RNA and comprised more than 87 species, among which 70 were novel, including divergent members of the genera Flavivirus and of the family Trypanosomatidae. Importantly, among these microbes we were able to reveal four distinguished types of human and/or mammalian pathogens, including two phleboviruses (hedi and wuxiang viruses), one novel Spotted fever group rickettsia, as well as a member of Leishmania donovani complex, among which hedi virus and Leishmania each had > 50% pool prevalence rate and relatively high abundance levels. Our study also showed the ubiquitous presence of an endosymbiont, namely Wolbachia, although no anti-viral or anti-pathogen effects were detected based on our data. In summary, our results uncovered the much un-explored diversity of microbes harboured by sandflies in China and demonstrated that high pathogen diversity and abundance are currently present in multiple populations, implying disease potential for exposed local human population or domestic animals.

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.

Recovery and genetic characterization of black queen cell virus

Black queen cell virus (BQCV) is a severe threat to the honeybee (Apis mellifera) worldwide. Although several BQCV strains have been reported in China, the molecular basis for BQCV pathogenicity has not been well understood. Thus, a reverse genetic system of BQCV is required for studying viral replication and its pathogenic mechanism. Here, the complete genome sequence of BQCV was obtained from honeybees using reverse transcription PCR (RT-PCR), namely a BQCV China-GS1 strain (KY741959). Then, a phylogenetic tree was built to analyse the genetic relationships among BQCV strains from different regions. Our results showed that the BQCV China-GS1 contained two ORFs, consistent with the known reference strains, except for the BQCV China-JL1 strain (KP119603). Furthermore, the infectious clone of BQCV was constructed based on BQCV China-GS1 using a low copy vector pACYC177 and gene recombination. Due to the lack of culture cells for bee viruses, we infected the healthy bees with infectious clone of BQCV, and the rescued BQCV resulted in the recovery of recombinant virus, which induced higher mortality than those of the control group. Immune response after inoculated with BQCV further confirmed that the infectious clone of BQCV caused the cellular and humoral immune response of honeybee (A. mellifera). In conclusion, the full nucleotide sequence of BQCV China-GS1 strain was determined, and the infectious clone of BQCV was constructed in this study. These data will improve the understanding of pathogenesis and the host immune responses to viral infection.