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

Development of a multiplex real-time quantitative reverse-transcription polymerase chain reaction for the detection of four bee viruses

Viruses in the families Dicistroviridae and Iflaviridae are among the main threats to western honey bees (Apis mellifera) and native bee species. Polymerase chain reaction (PCR) is the gold standard for pathogen detection in bees. However, high throughput screening for bee virus infections in singleplex PCR reactions is cumbersome and limited by the high quantities of sample RNA required. Thus, the development of a sensitive and specific multiplex PCR detection method for screening for multiple viruses simultaneously is necessary. Here, we report the development of a one-step multiplex reverse-transcription quantitative polymerase chain reaction (RT-qPCR) assay to detect four viruses commonly encountered in pollinator species. The optimized multiplex RT-qPCR protocol described in this study allows simultaneous detection of two dicistroviruses (Israeli acute paralysis virus and Black queen cell virus) and two iflaviruses (Sacbrood virus and Deformed wing virus) with high efficiency and specificity comparable to singleplex detection assays. This assay provides a broad range of detection and quantification, and the results of virus quantification in this study are similar to those performed in other studies using singleplex detection assays. This method will be particularly useful for data generation from small-bodied insect species that yield low amounts of RNA.

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.

Interaction of acetamiprid, Varroa destructor, and Nosema ceranae in honey bees

Health of honey bees is threatened by a variety of stressors, including pesticides and parasites. Here, we investigated effects of acetamiprid, Varroa destructor, and Nosema ceranae, which act either alone or in combination. Our results suggested that interaction between the three factors was additive, with survival risk increasing as the number of stressors increased. Although exposure to 150 μg/L acetamiprid alone did not negatively impact honey bee survival, it caused severe damage to midgut tissue. Among the three stressors, V. destructor posed the greatest threat to honey bee survival, and N. ceranae exacerbated intestinal damage and increased thickness of the midgut wall. Transcriptomic analysis indicated that different combinations of stressors elicited specific gene expression responses in honey bees, and genes involved in energy metabolism, immunity, and detoxification were altered in response to multiple stressor combinations. Additionally, genes associated with Toll and Imd signalling, tyrosine metabolism, and phototransduction pathway were significantly suppressed in response to different combinations of multiple stressors. This study enhances our understanding of the adaptation mechanisms to multiple stressors and aids in development of suitable protective measures for honey bees. ENVIRONMENTAL IMPLICATION: We believe our study is environmentally relevant for the following reasons: This study investigates combined effects of pesticide, Varroa destructor, and Nosema ceranae. These stressors are known to pose a threat to long-term survival of honey bees (Apis mellifera) and stability of the ecosystems. The research provides valuable insights into the adaptive mechanisms of honey bees in response to multiple stressors and developing effective conservation strategies. Further research can identify traits that promote honey bee survival in the face of future challenges from multiple stressors to maintain the overall stability of environment.

Parasite and virus dynamics in the honeybee Apis mellifera unicolor on a tropical island recently invaded by Varroa destructor

In La Réunion, the established honeybee subspecies Apis mellifera unicolor, an endemic subspecies of African lineage, is facing considerable challenges. Since the introduction of the Varroa destructor mite in 2017 high colony losses have been recorded. We investigated the dynamics of V. destructor and two viruses, the Deformed Wing Virus (DWV), known to be transmitted by the mite, and the Chronic Bee Paralysis Virus (CBPV), in A. m. unicolor. Colonies from two apiaries located at 300 and 900 m a.s.l were monitored twice for one year without any acaricide treatment. The brood area, V. destructor infestation rates, DWV and CBPV prevalence and load were recorded monthly. A. m. unicolor maintained brood rearing throughout the year. Varroa destructor infestation resulted in high colony mortality (up to 85 %) and high phoretic mite rates (up to 52 mites per hundred bees). The establishment of DWV in colonies occurred after that of V. destructor and the mite infestation rate had a significant effect on the virus prevalence and load. CBPV appeared only transiently throughout the surveys. The data showed that, in tropical colonies with permanent brood rearing, V. destructor and DWV can reach high levels, but are still subject to seasonal variations that appear to be influenced by environmental conditions. This suggests that beekeeping practices could be adapted by favouring sites and periods for transhumance or acaricide treatment.

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

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

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

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

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

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

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.

Toll receptor ligand Spätzle 4 responses to the highly pathogenic Enterococcus faecalis from Varroa mites in honeybees

Honeybees play a major role in crop pollination, which supports the agricultural economy and international food supply. The colony health of honeybees is threatened by the parasitic mite Varroa destructor, which inflicts physical injury on the hosts and serves as the vector for variable viruses. Recently, it shows that V. destructor may also transmit bacteria through the feeding wound, yet it remains unclear whether the invading bacteria can exhibit pathogenicity to the honeybees. Here, we incidentally isolate Enterococcus faecalis, one of the most abundant bacteria in Varroa mites, from dead bees during our routine generation of microbiota-free bees in the lab. In vivo tests show that E. faecalis is only pathogenic in Apis mellifera but not in Apis cerana. The expression of antimicrobial peptide genes is elevated following infection in A. cerana. The gene-based molecular evolution analysis identifies positive selection of genes encoding Späetzle 4 (Spz4) in A. cerana, a signaling protein in the Toll pathway. The amino acid sites under positive selection are related to structural changes in Spz4 protein, suggesting improvement of immunity in A. cerana. The knock-down of Spz4 in A. cerana significantly reduces the survival rates under E. faecalis challenge and the expression of antimicrobial peptide genes. Our results indicate that bacteria associated with Varroa mites are pathogenic to adult bees, and the positively selected gene Spz4 in A. cerana is crucial in response to this mite-related pathogen.

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.

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.

Establishment and Application of CRISPR-Cas12a-Based Recombinase Polymerase Amplification and a Lateral Flow Dipstick and Fluorescence for the Detection and Distinction of Deformed Wing Virus Types A and B

Deformed wing virus (DWV) is one of the important pathogens of the honey bee (Apis mellifera), which consists of three master variants: types A, B, and C. Among them, DWV types A (DWV-A) and B (DWV-B) are the most prevalent variants in honey bee colonies and have been linked to colony decline. DWV-A and DWV-B have different virulence, but it is difficult to distinguish them via traditional methods. In this study, we established a visual detection assay for DWV-A and DWV-B using recombinase polymerase amplification (RPA) and a lateral flow dipstick (LFD) coupled with the clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein (Cas) 12a fluorescence system (RPA-CRISPR-Cas12a-LFD). The limit of detection of this system was ~6.5 × 100 and 6.2 × 101 copies/μL for DWV-A and DWV-B, respectively. The assays were specific and non-cross-reactive against other bee viruses, and the results could be visualized within 1 h. The assays were validated by extracting cDNA from 36 clinical samples of bees that were suspected to be infected with DWV. The findings were consistent with those of traditional reverse transcription-quantitative polymerase chain reaction, and the RPA-CRISPR-Cas12a assay showed the specific, sensitive, simple, and appropriate detection of DWV-A and DWV-B. This method can facilitate the visual and qualitative detection of DWV-A and DWV-B as well as the monitoring of different subtypes, thereby providing potentially better control and preventing current and future DWV outbreaks.

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.

Consequences of microsporidian prior exposure for virus infection outcomes and bumble bee host health

Host-parasite interactions do not occur in a vacuum, but in connected multi-parasite networks that can result in co-exposures and coinfections of individual hosts. These can affect host health and disease ecology, including disease outbreaks. However, many host-parasite studies examine pairwise interactions, meaning we still lack a general understanding of the influence of co-exposures and coinfections. Using the bumble bee Bombus impatiens, we study the effects of larval exposure to a microsporidian Nosema bombi, implicated in bumble bee declines, and adult exposure to Israeli Acute Paralysis Virus (IAPV), an emerging infectious disease from honey bee parasite spillover. We hypothesize that infection outcomes will be modified by co-exposure or coinfection. Nosema bombi is a potentially severe, larval-infecting parasite, and we predict that prior exposure will result in decreased host resistance to adult IAPV infection. We predict double parasite exposure will also reduce host tolerance of infection, as measured by host survival. Although our larval Nosema exposure mostly did not result in viable infections, it partially reduced resistance to adult IAPV infection. Nosema exposure also negatively affected survival, potentially due to a cost of immunity in resisting the exposure. There was a significant negative effect of IAPV exposure on survivorship, but prior Nosema exposure did not alter this survival outcome, suggesting increased tolerance given the higher IAPV infections in the bees previously exposed to Nosema. These results again demonstrate that infection outcomes can be non-independent when multiple parasites are present, even when exposure to one parasite does not result in a substantial infection.

Gut microbiota assembly of Gotland varroa-surviving honey bees excludes major viral pathogens

The spread of the parasite Varroa destructor and associated viruses has resulted in massive honey bee colony losses with considerable economic and ecological impact. The gut microbiota has a major role in shaping honey bees tolerance and resistance to parasite infestation and viral infection, but the contribution of viruses to the assembly of the host microbiota in the context of varroa resistance and susceptibility remains unclear. Here, we used a network approach including viral and bacterial nodes to characterize the impact of five viruses, Apis Rhabdovirus-1 (ARV-1), Black Queen Cell virus (BQCV), Lake Sinai virus (LSV), Sacbrood virus (SBV) and Deformed wing virus (DWV) on the gut microbiota assembly of varroa-susceptible and Gotland varroa-surviving honey bees. We found that microbiota assembly was different in varroa-surviving and varroa-susceptible honey bees with the network of the latter having a whole module not present in the network of the former. Four viruses, ARV-1, BQCV, LSV, and SBV, were tightly associated with bacterial nodes of the core microbiota of varroa-susceptible honey bees, while only two viruses BQCV and LSV, appeared correlated with bacterial nodes in varroa-surviving honey bees. In silico removal of viral nodes caused major re-arrangement of microbial networks with changes in nodes centrality and significant reduction of the networks' robustness in varroa-susceptible, but not in varroa-surviving honey bees. Comparison of predicted functional pathways in bacterial communities using PICRUSt2 showed the superpathway for heme b biosynthesis from uroporphyrinogen-III and a pathway for arginine, proline, and ornithine interconversion as significantly increased in varroa-surviving honey bees. Notably, heme and its reduction products biliverdin and bilirubin have been reported as antiviral agents. These findings show that viral pathogens are differentially nested in the bacterial communities of varroa-surviving and varroa-susceptible honey bees. These results suggest that Gotland honey bees are associated with minimally-assembled and reduced bacterial communities that exclude viral pathogens and are resilient to viral nodes removal, which, together with the production of antiviral compounds, may explain the resiliency of Gotland honey bees to viral infections. In contrast, the intertwined virus-bacterium interactions in varroa-susceptible networks suggest that the complex assembly of microbial communities in this honey bee strain favor viral infections, which may explain viral persistence in this honey bee strain. Further understanding of protective mechanisms mediated by the microbiota could help developing novel ways to control devastating viral infections affecting honey bees worldwide.

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.

Local habitat type influences bumble bee pathogen loads and bee species distributions

Bumble bees (Hymenoptera: Apidae, Bombus Latreille) perform important ecological services in both managed and natural ecosystems. Anthropogenically induced change has altered floral resources, climate, and insecticide exposure, factors that impact health and disease levels in these bees. Habitat management presents a solution for improving bee health and biodiversity, but this requires better understanding of how different pathogens and bee species respond to habitat conditions. We take advantage of the washboard of repeated ridges (forested) and valleys (mostly developed) in central Pennsylvania to examine whether local variation in habitat type and other landscape factors influence bumble bee community composition and levels of 4 leading pathogens in the common eastern bumble bee, Bombus impatiens Cresson. Loads of viruses (DWV and BQCV) were found to be lowest in forest habitats, whereas loads of a gut parasite, Crithidia bombi, were highest in forests. Ridgetop forests hosted the most diverse bumble bee communities, including several habitat specialists. B. impatiens was most abundant in valleys, and showed higher incidence in areas of greater disturbance, including more developed, unforested, and lower floral resource sites, a pattern which mirrors its success in the face of anthropogenic change. Additionally, DNA barcoding revealed that B. sandersoni is much more common than is apparent from databases. Our results provide evidence that habitat type can play a large role in pathogen load dynamics, but in ways that differ by pathogen type, and point to a need for consideration of habitat at both macro-ecological and local spatial scales.

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.

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.

Sunflower-Associated Reductions in Varroa Mite Infestation of Honey Bee Colonies

Landscapes can affect parasite epidemiology in wild and agricultural animals. Honey bees are threatened by loss of floral resources and by parasites, principally the mite Varroa destructor and the viruses it vectors. Existing mite control relies heavily on chemical treatments that can adversely affect bees. Alternative, pesticide-free control methods are needed to mitigate infestation with these ectoparasites. Many flowering plants provide nectar and pollen that confer resistance to parasites. Enrichment of landscapes with antiparasitic floral resources could therefore provide a sustainable means of parasite control in pollinators. Floral rewards of Asteraceae plants can reduce parasitic infection in diverse bee species, including honey and bumble bees. Here, we tested the effects of sunflower (Helianthus annuus) cropland and pollen supplementation on honey bee resistance to macro- and microparasites. Although sunflower had nonsignificant effects on microparasites, We found that increased sunflower pollen availability correlated with reduced Varroa mite infestation in landscapes and pollen-supplemented colonies. At the landscape level, each doubling of sunflower crop area was associated with a 28% reduction in mite infestation. In field trials, late-summer supplementation of colonies with sunflower pollen reduced mite infestation by 2.75-fold relative to artificial pollen. United States sunflower crop acreage has declined by 2% per year since 1980, however, suggesting reduced availability of this floral resource. Although further research is needed to determine whether the observed effects represent direct inhibition of mite fecundity or mite-limiting reductions in honey bee brood-rearing, our findings suggest the potential for sunflower plantings or pollen supplements to counteract a major driver of honey bee losses worldwide.

A Novel Variant of Deformed Wing Virus (DWV) from the Invasive Honeybee Apis florea (Apidae, Hymenoptera) and Its Ectoparasite Euvarroa sinhai (Acarina, Mesostigmata) in Taiwan

The invasion of Apis florea in Taiwan was first recorded in 2017. The deformed wing virus (DWV) has been identified as a common bee virus in apiculture around the world. Ectoparasitic mites are the main DWV vector for horizontal transmission. However, there are few studies about the ectoparasitic mite of Euvarroa sinhai, which has been found in A. florea. In this study, the prevalence of DWV among four hosts, including A. florea, Apis mellifera, E. sinhai, and Varroa destructor, was determined. The results showed that a high DWV-A prevalence rate in A. florea, ranging from 69.2% to 94.4%, was detected. Additionally, the genome of DWV isolates was sequenced and subjected to phylogenetic analysis based on the complete polyprotein sequence. Furthermore, isolates from A. florea and E. sinhai both formed a monophyletic group for the DWV-A lineage, and the sequence identity was 88% between the isolates and DWV-A reference strains. As noted above, two isolates could be the novel DWV strain. It cannot be excluded that novel DWV strains could pose an indirect threat to sympatric species, such as A. mellifera and Apis cerana.

Current Knowledge on Bee Innate Immunity Based on Genomics and Transcriptomics

As important pollinators, bees play a critical role in maintaining the balance of the ecosystem and improving the yield and quality of crops. However, in recent years, the bee population has significantly declined due to various pathogens and environmental stressors including viruses, bacteria, parasites, and increased pesticide application. The above threats trigger or suppress the innate immunity of bees, their only immune defense system, which is essential to maintaining individual health and that of the colony. In addition, bees can be divided into solitary and eusocial bees based on their life traits, and eusocial bees possess special social immunities, such as grooming behavior, which cooperate with innate immunity to maintain the health of the colony. The omics approach gives us an opportunity to recognize the distinctive innate immunity of bees. In this regard, we summarize innate bee immunity from a genomic and transcriptomic perspective. The genetic characteristics of innate immunity were revealed by the multiple genomes of bees with different kinds of sociality, including honeybees, bumblebees, wasps, leaf-cutter bees, and so on. Further substantial transcriptomic data of different tissues from diverse bees directly present the activation or suppression of immune genes under the infestation of pathogens or toxicity of pesticides.

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.

Honey Bee Pathogen Prevalence and Interactions within the Marmara Region of Turkey

Beekeeping has yet to reach its full potential in terms of productivity in Turkey where it has a relatively large role in the economy. Poor colony health is suspected to be the reason for this, but comprehensive disease monitoring programs are lacking to support this notion. We sampled a total of 115 colonies across five different apiaries throughout the Marmara region of Turkey and screened for all of the major bee pathogens using PCR and RNA-seq methods. We found that Varroa mites are more prevalent in comparison to Nosema infections. The pathogens ABPV, DWV, KV, and VDV1 are near 100% prevalent and are the most abundant across all locations, which are known to be vectored by the Varroa mite. We therefore suspect that controlling Varroa mites will be key for improving bee health in Turkey moving forward. We also documented significant interactions between DWV, KV, and VDV1, which may explain how the more virulent strain of the virus becomes abundant. ABPV had a positive interaction with VDV1, thereby possibly facilitating this more virulent viral strain, but a negative interaction with Nosema ceranae. Therefore, these complex pathogen interactions should be taken into consideration in the future to improve bee health.

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.

Complete genome analysis of a novel iflavirus from the spotted lanternfly Lycorma delicatula

The spotted lanternfly (Lycorma delicatula) is an invasive pest that causes serious economic losses in fruit and wood production. Here, we identified a novel iflavirus named "Lycorma delicatula iflavirus 1" (LDIV1), in a spotted lanternfly. The full genome sequence of LDIV1 is 10,222 nt in length and encodes a polyprotein containing a picornavirus capsid-protein-domain-like domain, a cricket paralysis virus capsid superfamily domain, an RNA helicase domain, a peptidase C3 superfamily domain, and an RNA-dependent RNA polymerase (RdRp) domain. LDIV1 replicates in the host insect and activates small interfering RNA (siRNA)-based host antiviral immunity. Phylogenetic analysis demonstrated that LDIV1 is most closely related to an unspecified member of the order Picornavirales, with 61.7% sequence identity in the RdRp region and 57.6% sequence identity in the coat protein region, and thus meets the demarcation criteria for new species in the genus Iflavirus. To the best of our knowledge, LDIV1 is the first virus discovered in L. delicatula.

Honey Bee Genetic Stock Determines Deformed Wing Virus Symptom Severity but not Viral Load or Dissemination Following Pupal Exposure

Honey bees exposed to Varroa mites incur substantial physical damage in addition to potential exposure to vectored viruses such as Deformed wing virus (DWV) that exists as three master variants (DWV-A, DWV-B, and DWV-C) and recombinants. Although mite-resistant bees have been primarily bred to mitigate the impacts of Varroa mites, mite resistance may be associated with increased tolerance or resistance to the vectored viruses. The goal of our study is to determine if five honey bee stocks (Carniolan, Italian, Pol-Line, Russian, and Saskatraz) differ in their resistance or tolerance to DWV based on prior breeding for mite resistance. We injected white-eyed pupae with a sublethal dose (10⁵) of DWV or exposed them to mites and then evaluated DWV levels and dissemination and morphological symptoms upon adult emergence. While we found no evidence of DWV resistance across stocks (i.e., similar rates of viral replication and dissemination), we observed that some stocks exhibited reduced symptom severity suggestive of differential tolerance. However, DWV tolerance was not consistent across mite-resistant stocks as Russian bees were most tolerant, while Pol-Line exhibited the most severe symptoms. DWV variants A and B exhibited differential dissemination patterns that interacted significantly with the treatment group but not bee stock. Furthermore, elevated DWV-B levels reduced adult emergence time, while both DWV variants were associated with symptom likelihood and severity. These data indicate that the genetic differences underlying bee resistance to Varroa mites are not necessarily correlated with DWV tolerance and may interact differentially with DWV variants, highlighting the need for further work on mechanisms of tolerance and bee stock–specific physiological interactions with pathogen variants.

Colony-Level Viral Load Influences Collective Foraging in Honey Bees

Nutrition is an important component of social insect colony health especially in the face of stressors such as parasitism and viral infections. Honey bees are known to preferentially select nectar and pollen based on macronutrient and phytochemical contents and in response to pathogen loads. However, given that honey bees live in colonies, collective foraging decisions may be impacted directly by forager infection status but also by colony health. This field experiment was conducted to determine if honey bee viral infections are correlated with pollen and nectar foraging and if these associations are impacted more by colony or forager infection. By comparing regressions with and without forager and colony variables and through structural equation models, we were able to determine the relative contributions of colony and forager virus loads on forager decisions. We found that foragers had higher numbers and levels of BQCV and CBPV but lower levels of DWV viruses than their respective colonies. Overall, individuals appeared to forage based a combination of their own and colony health but with greater weight given to colony metrics. Colony parasitism by Varroa mites, positively correlated with both forager and colony DWV-B levels, was negatively associated with nectar weight. Further, colony DWV-B levels were negatively associated with individually foraged pollen protein: lipid ratios but positively correlated with nectar weight and sugar content. This study shows that both colony and forager health can simultaneously mediate individual foraging decisions and that the importance of viral infections and parasite levels varies with foraging metrics. Overall, this work highlights the continued need to explore the interactions of disease, nutrition, and genetics in social interactions and structures.

The Virome of Healthy Honey Bee Colonies: Ubiquitous Occurrence of Known and New Viruses in Bee Populations

The Western honey bee, Apis mellifera , is a vital part of our ecosystem as well as cultural heritage. Annual colony losses endanger beekeeping.

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.

The antiviral role of NF-κB-mediated immune responses and their antagonism by viruses in insects

The antiviral role of innate immune responses mediated by the NF-κB family of transcription factors is well established in vertebrates but was for a long time less clear in insects. Insects encode two canonical NF-κB pathways, the Toll and Imd ('immunodeficiency') pathways, which are best characterised for their role in antibacterial and antifungal defence. An increasing body of evidence has also implicated NF-κB-mediated innate immunity in antiviral responses against some, but not all, viruses. Specific pattern recognition receptors (PRRs) and molecular events leading to NF-κB activation by viral pathogen-associated molecular patterns (PAMPs) have been elucidated for a number of viruses and insect species. Particularly interesting are recent findings indicating that the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway detects viral RNA to activate NF-κB-regulated gene expression. We summarise the literature on virus-NF-κB pathway interactions across the class Insecta, with a focus on the dipterans Drosophila melanogaster and Aedes aegypti. We discuss potential reasons for differences observed between different virus-host combinations, and highlight similarities and differences between cGAS-STING signalling in insects versus vertebrates. Finally, we summarise the increasing number of known molecular mechanisms by which viruses antagonise NF-κB responses, which suggest that NF-κB-mediated immunity exerts strong evolutionary pressures on viruses. These developments in our understanding of insect antiviral immunity have relevance to the large number of insect species that impact on humans through their transmission of human, livestock and plant diseases, exploitation as biotechnology platforms, and role as parasites, pollinators, livestock and pests.

The Multiplicity of Infection of Recombinant Vaccinia Virus Expressing the T7 RNA Polymerase Determines the Rescue Efficiency of Vesicular Stomatitis Virus

Vesicular stomatitis virus (VSV) has a wide range of cell tropism, making it a prototype of studying the negative-strand RNA virus (NSRV), including virus–host interactions and vaccine development. Although VSV rescue systems have been progressively optimized throughout time, the T7-based expression system is the most commonly utilized to rescue VSV. However, it remains a significant barrier for many labs. In our study, we found that rescue VSV’s efficiency is associated with the various multiplicities of infection (MOIs) of recombinant vaccinia virus expressing the T7 RNA polymerase (vTF-7.3). It works at maximum efficiency while the MOI of vTF-7.3 is 5, which is analyzed by quantitative PCR, Western blot, and flow cytometry, compared to the lowest rescue level with MOI of 1. Meanwhile, our data also suggest that purification of vTF-7.3 prior to transfection is a prerequisite for VSV rescue. Overall, our study reveals for the first time a precise correlation between vTF-7.3 and rescue efficiency, which may aid in resolving the uncertainties in the quest to build the VSV reverse genetic system.

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.

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.

Virome Analysis Reveals Diverse and Divergent RNA Viruses in Wild Insect Pollinators in Beijing, China

Insect pollinators provide major pollination services for wild plants and crops. Honeybee viruses can cause serious damage to honeybee colonies. However, viruses of other wild pollinating insects have yet to be fully explored. In the present study, we used RNA sequencing to investigate the viral diversity of 50 species of wild pollinating insects. A total of 3 pathogenic honeybee viruses, 8 previously reported viruses, and 26 novel viruses were identified in sequenced samples. Among these, 7 novel viruses were shown to be closely related to honeybee pathogenic viruses, and 4 were determined to have potential pathogenicity for their hosts. The viruses detected in wild insect pollinators were mainly from the order Picornavirales and the families Orthomyxoviridae, Sinhaliviridae, Rhabdoviridae, and Flaviviridae. Our study expanded the species range of known insect pollinator viruses, contributing to future efforts to protect economic honeybees and wild pollinating insects.

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.

Context-Dependent Viral Transgenerational Immune Priming in Honey Bees (Hymenoptera: Apidae)

Transgenerational immune priming is the process of increased resistance to infection in offspring due to parental pathogen exposure. Honey bees (Apis mellifera L. (Hymenoptera: Apidae)) are hosts to multiple pathogens, and this complex immune function could help protect against overwhelming infection. Honey bees have demonstrated transgenerational immune priming for the bacterial pathogen Paenibacillus larvae; however, evidence for viral transgenerational immune priming is lacking across insects in general. Here we test for the presence of transgenerational immune priming in honey bees with Deformed wing virus (DWV) by injecting pupae from DWV-exposed queens and measuring virus titer and immune gene expression. Our data suggest that there is evidence for viral transgenerational immune priming in honey bees, but it is highly context-dependent based on route of maternal exposure and potentially host genetics or epigenetic factors.

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.

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.

Chemical Stimulants and Stressors Impact the Outcome of Virus Infection and Immune Gene Expression in Honey Bees (Apis mellifera)

Western honey bees (Apis mellifera) are ecologically, agriculturally, and economically important plant pollinators. High average annual losses of honey bee colonies in the US have been partially attributed to agrochemical exposure and virus infections. To examine the potential negative synergistic impacts of agrochemical exposure and virus infection, as well as the potential promise of phytochemicals to ameliorate the impact of pathogenic infections on honey bees, we infected bees with a panel of viruses (i.e., Flock House virus, deformed wing virus, or Sindbis virus) and exposed to one of three chemical compounds. Specifically, honey bees were fed sucrose syrup containing: (1) thyme oil, a phytochemical and putative immune stimulant, (2) fumagillin, a beekeeper applied fungicide, or (3) clothianidin, a grower-applied insecticide. We determined that virus abundance was lower in honey bees fed 0.16 ppb thyme oil augmented sucrose syrup, compared to bees fed sucrose syrup alone. Parallel analysis of honey bee gene expression revealed that honey bees fed thyme oil augmented sucrose syrup had higher expression of key RNAi genes (argonaute-2 and dicer-like), antimicrobial peptide expressing genes (abaecin and hymenoptaecin), and vitellogenin, a putative honey bee health and age indicator, compared to bees fed only sucrose syrup. Virus abundance was higher in bees fed fumagillin (25 ppm or 75 ppm) or 1 ppb clothianidin containing sucrose syrup relative to levels in bees fed only sucrose syrup. Whereas, honey bees fed 10 ppb clothianidin had lower virus levels, likely because consuming a near lethal dose of insecticide made them poor hosts for virus infection. The negative impact of fumagillin and clothianidin on honey bee health was indicated by the lower expression of argonaute-2, dicer-like, abaecin, and hymenoptaecin, and vitellogenin. Together, these results indicate that chemical stimulants and stressors impact the outcome of virus infection and immune gene expression in honey bees.

Processing Technologies for Bee Products: An Overview of Recent Developments and Perspectives

Increased demand for a more balanced, healthy, and safe diet has accelerated studies on natural bee products (including honey, bee bread, bee collected pollen royal jelly, propolis, beeswax, and bee venom) over the past decade. Advanced food processing techniques, such as ultrasonication and microwave and infrared (IR) irradiation, either has gained popularity as alternatives or combined with conventional processing techniques for diverse applications in apiculture products at laboratory or industrial scale. The processing techniques used for each bee products have comprehensively summarized in this review, including drying (traditional drying, infrared drying, microwave-assisted traditional drying or vacuum drying, and low temperature high velocity-assisted fluidized bed drying), storage, extraction, isolation, and identification; the assessment methods related to the quality control of bee products are also fully mentioned. The different processing techniques applied in bee products aim to provide more healthy active ingredients largely and effectively. Furthermore, improved the product quality with a shorter processing time and reduced operational cost are achieved using conventional or emerging processing techniques. This review will increase the positive ratings of the combined new processing techniques according to the needs of the bee products. The importance of the models for process optimization on a large scale is also emphasized in the future.

Occurrence and Phylogenetic Analysis of DWV in Stingless Bee (Apidae sp.) in China: A Case Report

Honey bees play a vital role in providing pollination services for agricultural crops and wild flowering plants. However, the spillover risk of their pathogens to other pollinators or wild insects is becoming a cause for concern. There is some evidence that stingless bees can carry honey bee viruses, but little is known about the presence of honey bee viruses in stingless bees in China. Here, we investigate the occurrence of major honey bee pathogens including bacteria, fungi, and viruses in stingless bees (Apidae: sp.). Our results show that the stingless bees (Apidae: sp.) were mainly infected with DWV-A, but no DWV-B and DWV-C. Phylogenetic analysis on fragments of lp, RdRp, and VP3 of DWV-A indicated that genetic variation in VP3 might an important indicator for host-specific viruses, but it requires further study. Our results indicated that DWV-A is not only the major strain of virus currently circulating in managed bee colonies in China and globally, but in stingless bee species as a whole.