Adapted tolerance to virus infections in four geographically distinct Varroa destructor-resistant honeybee populations

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.

Processing of the 3C/D Region of the Deformed Wing Virus (DWV)

The deformed wing virus (DWV) belongs to the genus Iflavirus and the family Iflaviridae within the order Picornavirales. It is an important pathogen of the Western honey bee, Apis mellifera, causing major losses among honey bee colonies in association with the ectoparasitic mite Varroa destructor. Although DWV is one of the best-studied insect viruses, the mechanisms of viral replication and polyprotein processing have been poorly studied in the past. We investigated the processing of the protease-polymerase region at the C-terminus of the polyprotein in more detail using recombinant expression, novel serological reagents, and virus clone mutagenesis. Edman degradation of purified maturated polypeptides uncovered the C- and N-termini of the mature 3C-like (3CL) protease and RNA-dependent RNA polymerase (3DL, RdRp), respectively. Autocatalytic processing of the recombinant DWV 3CL protease occurred at P1 Q2118 and P1' G2119 (KPQ/GST) as well as P1 Q2393 and P1' S2394 (HAQ/SPS) cleavage sites. New monoclonal antibodies (Mab) detected the mature 3CL protease with an apparent molecular mass of 32 kDa, mature 3DL with an apparent molecular mass of 55 kDa as well as a dominant 3CDL precursor of 90 kDa in DWV infected honey bee pupae. The observed pattern corresponds well to data obtained via recombinant expression and N-terminal sequencing. Finally, we were able to show that 3CL protease activity and availability of the specific protease cleavage sites are essential for viral replication, protein synthesis, and establishment of infection using our molecular clone of DWV-A.

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.

Epidemiology, factors influencing prevalence and level of varroosis infestation (Varroa destructor) in honeybee (Apis mellifera) colonies in different agroecologies of Southwest Ethiopia

Little information is available on the epidemiology of varroosis caused by Varroa mite, Varroa destructor infestation in Ethiopia, although it is a devastating honeybee disease that results in significant economic losses in beekeeping. Therefore, between October 2021 and October 2022, a cross-sectional study was carried out in different agroecology zones in Southwest Ethiopia to determine the prevalence and associated risk factors for varroosis, as well as the effects of this disease on honeybee colonies and honey production. A multivariate logistic regression analysis was performed to identify possible risk factors for the prevalence of V. destructor. A total of 384 adult honeybee and worker or drone brood samples were collected from honeybee colonies and examined using standard diagnostic techniques in the laboratory. The result shows that the prevalence of V. destructor was found to be 39.3% (95% CI 34.44-44.21) and 43.2% (38.27-48.18) in adult honeybees and brood, respectively. The major risk factors for the prevalence of V. destructor in the study areas included agroecology (OR = 5.2, 95% CI 1.75-14.85), type of hive (OR = 2.9, 95% CI 1.17-17.03), management system (OR = 4.3, 95% CI 1.23-14.70), and colony management (OR = 3.5, 95% CI 1.31-9.14). The lower level of colony infestation in adult bees and brood was measured as 1.97 ± 0.14 and 3.19 ± 0.25, respectively. Season, colony status, colony management, and agroecology were among the determinant factors of the level of varroa mite infestation in adult bees and brood. The results of the study demonstrated that honey production losses are largely attributable to V. destructor infestation. Therefore, it is critical to inform the community about the effects of V. destructor on honey production and develop and implement effective management strategies for this disease. In addition, further research should be done to identify and isolate additional factors that contribute to varroosis in honeybees in different regions.

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.

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.

Using Surrogate Insects in Acid Bioassays for Development of New Controls for Varroa destructor (Arachnida: Varroidae)

Resistance to traditional synthetic compounds by Varroa destructor Anderson and Trueman and shortcomings of the organic acid class of acaracides commonly used in varroa management requires continual development of new controls. V. destructor, however, are difficult to obtain for use in control bioassays because they are obligate parasites that cannot be easily reared outside of a honey bee colony. We conducted bioassays using other, more easily obtainable species to find organisms that could be used as surrogates for V. destructor when testing new potential controls. We compared the toxicities of acetic acid, lactic acid, formic acid, and oxalic acid at 0.005%, 0.05%, 0.5%, 5%, and 50% (20% oxalic acid only) concentrations based on natural volatility (nonheated) for the control of two beetle species, Oryzaephilus surinamensis L. and Alphitobius diaperinus Panzer, greater wax moth larvae, Galleria mellonella L., and V. destructor. The assay results were consistent across all species with formic acid and acetic acid showing 100% mortality of all four test species at 50% concentration. The assays also provided insight into the method of application (vaporization or contact) needed to cause mortality. Our results show that other organisms can be used in place of V. destructor for initial testing of acids and possibly other chemicals for control of the ectoparasite.

Residues of Pesticides and Heavy Metals in Polish Varietal Honey

The levels of chemical pollutants were determined in 30 samples of varietal honey from southeastern Poland, including 223 pesticides (insecticides, herbicides, fungicides, acaricides, plant growth regulators, and veterinary drugs) and 5 heavy metals (Pb, Cd, Hg, Cu, and Zn). In 10% of the samples, no pesticide residues were found. The most frequently identified pesticides were thiacloprid (90% of the samples, max 0.337 mg/kg), acetamiprid (86.6%, max 0.061 mg/kg), carbendazim (60%, max 0.049 mg/kg), DMF (56.6%, max 0.038 mg/kg), total amitraz (53.3%, max 0.075 mg/kg), thiamethoxam (26.6%, max 0.004 mg/kg), thiacloprid-amide (13.3%, max 0.012 mg/kg), dimethoate (10%, max 0.003 mg/kg), azoxystrobin (10%, max 0.002 mg/kg), tebuconazole (6.66%, max 0.002 mg/kg), and boscalid (3.33%, max 0.001 mg/kg). The acceptable limits for the compounds were not exceeded in any sample. The Pb content ranged between 0.044 and 0.081 mg/kg. The concentration of Hg and Cd did not exceed 5.0 µg/kg and 0.02 mg/kg, respectively. The honey variety significantly (p < 0.01) influenced the content of Cu, which ranged from 0.504 (rapeseed honey) to 1.201 mg/kg (buckwheat). A similar tendency (p > 0.05) was observed for the Zn content, which ranged from 0.657 mg/kg (linden) to 2.694 mg/kg (buckwheat). Honey produced in southeastern Poland was shown to be safe for human consumption.

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.

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.

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

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

Honey bee viruses are highly prevalent but at low intensities in wild pollinators of cucurbit agroecosystems

Managed and wild bee populations are in decline around the globe due to several biotic and abiotic stressors. Pathogenic viruses associated with the Western honey bee (Apis mellifera) have been identified as key contributors to losses of managed honey bee colonies, and are known to be transmitted to wild bee populations through shared floral resources. However, little is known about the prevalence and intensity of these viruses in wild bee populations, or how bee visitation to flowers impacts viral transmission in agroecosystems. This study surveyed honey bee, bumble bee (Bombus impatiens) and wild squash bee (Eucera (Peponapis) pruinosa) populations in Cucurbita agroecosystems across Pennsylvania (USA) for the prevalence and intensity of five honey bee viruses: acute bee paralysis virus (ABPV), deformed wing virus (DWV), Israeli acute paralysis virus (IAPV), Kashmir bee virus (KBV), and slow bee paralysis virus (SBPV). We investigated the potential role of bee visitation rate to flowers on DWV intensity among species in the pollinator community, with the expectation that increased bee visitation to flowers would increase the opportunity for transmission events between host species. We found that honey bee viruses are highly prevalent but in lower titers in wild E. pruinosa and B. impatiens than in A. mellifera populations throughout Pennsylvania (USA). DWV was detected in 88% of B. impatiens, 48% of E. pruinosa, and 95% of A. mellifera. IAPV was detected in 5% of B. impatiens and 4% of E. pruinosa, compared to 9% in A. mellifera. KBV was detected in 1% of B. impatiens and 5% of E. pruinosa, compared to 32% in A. mellifera. Our results indicate that DWV titers are not correlated with bee visitation in Cucurbita fields. The potential fitness impacts of these low viral titers detected in E. pruinosa remain to be investigated.

Parallel evolution of Varroa resistance in honey bees: a common mechanism across continents?

The near-globally distributed ecto-parasitic mite of the Apis mellifera honeybee, Varroa destructor, has formed a lethal association with Deformed wing virus, a once rare and benign RNA virus. In concert, the two have killed millions of wild and managed colonies, particularly across the Northern Hemisphere, forcing the need for regular acaricide application to ensure colony survival. However, despite the short association (in evolutionary terms), a small but increasing number of A. mellifera populations across the globe have been surviving many years without any mite control methods. This long-term survival, or Varroa resistance, is consistently associated with the same suite of traits (recapping, brood removal and reduced mite reproduction) irrespective of location. Here we conduct an analysis of data extracted from 60 papers to illustrate how these traits connect together to explain decades of mite resistance data. We have potentially a unified understanding of natural Varroa resistance that will help the global industry achieve widespread miticide-free beekeeping and indicate how different honeybee populations across four continents have resolved a recent threat using the same suite of behaviours.