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.

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.

An invasive ant increases deformed wing virus loads in honey bees.

The majority of invasive species are best known for their effects as predators. However, many introduced predators may also be substantial reservoirs for pathogens. Honey bee-associated viruses are found in various arthropod species including invasive ants. We examined how the globally invasive Argentine ant (Linepithema humile), which can reach high densities and infest beehives, is associated with pathogen dynamics in honey bees. Viral loads of deformed wing virus (DWV), which has been linked to millions of beehive deaths around the globe, and black queen cell virus significantly increased in bees when invasive ants were present. Microsporidian and trypanosomatid infections, which are more bee-specific, were not affected by ant invasion. The bee virome in autumn revealed that DWV was the predominant virus with the highest infection levels and that no ant-associated viruses were infecting bees. Viral spillback from ants could increase infections in bees. In addition, ant attacks could pose a significant stressor to bee colonies that may affect virus susceptibility. These viral dynamics are a hidden effect of ant pests, which could have a significant impact on disease emergence in this economically important pollinator. Our study highlights a perhaps overlooked effect of species invasions: changes in pathogen dynamics.

Effect of Chronic Exposure to Sublethal Doses of Imidacloprid and Nosema ceranae on Immunity, Gut Microbiota, and Survival of Africanized Honey Bees.

Large-scale honey bee colony losses reported around the world have been associated with intoxication with pesticides, as with the presence of pests and pathogens. Among pesticides, neonicotinoid insecticides are the biggest threat. Due to their extensive use, they can be found in all agricultural environments, including soil, water, and air, are persistent in the environment, and are highly toxic for honey bees. In addition, infection by different pests and pathogens can act synergistically, weakening bees. In this study, we investigated the effects of chronic exposure to sublethal doses of imidacloprid alone or combined with the microsporidia Nosema ceranae on the immune response, deformed wing virus infection (DWV), gut microbiota, and survival of Africanized honey bees. We found that imidacloprid affected the expression of some genes associated with immunity generating an altered physiological state, although it did not favor DWV or N. ceranae infection. The pesticide alone did not affect honey bee gut microbiota, as previously suggested, but when administered to N. ceranae infected bees, it generated significant changes. Finally, both stress factors caused high mortality rates. Those results illustrate the negative impact of imidacloprid alone or combined with N. ceranae on Africanized honey bees and are useful to understand colony losses in Latin America.

No impact of hygienic behavior and viral coinfection on the development of European foulbrood in honey bee (Apis mellifera) colonies during blueberry pollination in Michigan.

European foulbrood (EFB) is a severe disease of honey bee (Apis mellifera) larvae caused by the bacterium Linnaeus [Hymenoptera: Apidae]) Melissococcus plutonius (ex White) Bailey and Collins (Lactobacillales: Enterococcaceae). Many beekeepers in North America report severe EFB following blueberry pollination, but it is not clear what factors during pollination are related to clinical disease. Additionally, the impact that other factors such as viral load and hygienic behavior have on EFB has not been studied. In Spring of 2020 we enrolled 60 commercial honey bee colonies in a prospective cohort study. Colonies were inspected 3 times over the season with hive metrics and samples taken for viral testing. Each colony was tested for hygienic behavior twice and the score was averaged. Viral loads were determined by qPCR for deformed wing virus (DWV) A and B. We found no statistical difference in the EFB prevalence or severity between the 2 yards at any timepoint; 50% (n = 16) of the colonies in the holding yard and 63% (n = 17) in blueberry developed moderate to severe EFB over the study period. When colonies from both yards were pooled, we found no relationship between viral load or hygienic behavior and development of EFB. These results suggest that other factors may be responsible for driving EFB virulence and hygienic behavior is not likely helpful in managing this disease.

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.

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.

Multiple benefits of breeding honey bees for hygienic behavior.

Honey bee colonies are prone to invasion by pests and pathogens. The combination of the parasitic mite Varroa destructor (Varroa) and the multiple viruses it vectors, is a major driver of colony losses. Breeding for hygienic behavior to reduce Varroa populations is considered a sustainable way to reduce the impact of Varroa on honey bee health. However, hygienic behavior may have a cost to the health of individual bees, both in terms of viral infection risk and immune function. To determine whether selection for hygienic behavior at the colony level is associated with trade-offs in honey bee viral infection and immune function, we compared Varroa populations, viral loads, and individual immune function between honey bee colonies that were bred for high and low hygienic behavior. Specifically, we measured Varroa infestation, Deformed wing virus DWV-A, DWV-B, Acute bee paralysis virus (ABPV), and Israeli acute paralysis virus IAPV viral genome levels in bee samples from artificially inseminated queens in our bi-directional selection program for hygienic behavior in Israel. In addition, we evaluated the expression of 12 genes from the Jak-STAT, Toll, IMD and RNAi immune pathways. We found significantly lower Varroa infestation and DWV loads in highly hygienic colonies than in colonies exhibiting low hygienic behavior. In addition, workers of the hygienic colonies had significantly higher expression of the immune genes PGRP-S2 and hymenoptaecin compared to workers from low hygienic colonies. These results indicate no trade-offs in breeding for hygienic behavior. Hygienic honey bees were associated with reduced Varroa populations and reduced DWV prevalence or load at the colony level. Individual immunity of hygienic bees was increased, which could also contribute to lower virus levels, although lower Varroa levels due to social immunity presumably contributed as well. In sum, we demonstrate multiple health benefits of breeding for honey bee hygiene.

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.

Impact of Honey Bee Migratory Management on Pathogen Loads and Immune Gene Expression is Affected by Complex Interactions With Environment, Worker Life History, and Season.

The effects of honey bee management, such as intensive migratory beekeeping, are part of the ongoing debate concerning causes of colony health problems. Even though comparisons of disease and pathogen loads among differently managed colonies indicate some effects, the direct impact of migratory practices on honey bee pathogens is poorly understood. To test long- and short-term impacts of managed migration on pathogen loads and immunity, experimental honey bee colonies were maintained with or without migratory movement. Individuals that experienced migration as juveniles (e.g., larval and pupal development), as adults, or both were compared to control colonies that remained stationary and therefore did not experience migratory relocation. Samples at different ages and life-history stages (hive bees or foragers), taken at the beginning and end of the active season, were analyzed for pathogen loads and physiological markers of health. Bees exposed to migratory management during adulthood had increased levels of the AKI virus complex (Acute bee paralysis, Kashmir bee, and Israeli acute bee paralysis viruses) and decreased levels of antiviral gene expression (dicer-like). However, those in stationary management as adults had elevated gut parasites (i.e. trypanosomes). Effects of environment during juvenile development were more complex and interacted with life-history stage and season. Age at collection, life-history stage, and season all influenced numerous factors from viral load to immune gene expression. Although the factors that we examined are not independent, the results illuminate potential factors in both migratory and nonmigratory beekeeping that are likely to contribute to colony stress, and also indicate potential mitigation measures.

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.

The Effect of Covert and Overt Infections on Disease Dynamics in Honey-Bee Colonies.

Viral diseases of honey bees are important economically and ecologically and have been widely modelled. The models reflect the fact that, in contrast to the typical case for vertebrates, invertebrates cannot acquire immunity to a viral disease, so they are of SIS or (more often) SI type. Very often, these diseases may be transmitted vertically as well as horizontally, by vectors as well as directly, and through the environment, although models do not generally reflect all these transmission mechanisms. Here, we shall consider an important additional complication the consequences of which have yet to be fully explored in a model, namely that both infected honey bees and their vectors may best be described using more than one infection class. For honey bees, we consider three infection classes. Covert infections occur when bees have the virus under control, such that they do not display symptoms of the disease, and are minimally or not at all affected by it. Acutely overtly infected bees often exhibit severe symptoms and have a greatly curtailed lifespan. Chronically overtly infected bees typically have milder symptoms and a moderately shortened lifespan. For the vector, we consider just two infection classes which are covert infected and overt infected as has been observed in deformed-wing virus (DWV) vectored by varroa mites. Using this structure, we explore the impact of spontaneous transition of both mites and bees from a covertly to an overtly infected state, which is also a novel element in modelling viral diseases of honey bees made possible by including the different infected classes. The dynamics of these diseases are unsurprisingly rather different from the dynamics of a standard SI or SIS disease. In this paper, we highlight how our compartmental structure for infection in honey bees and their vectors impact the disease dynamics observed, concentrating in particular on DWV vectored by varroa mites. If there is no spontaneous transition, then a basic reproduction number [Formula: see text] exists. We derive a condition for [Formula: see text] that reflects the complexities of the system, with components for vertical and for direct and vector-mediated horizontal transmission, using the directed graph of the next-generation matrix of the system. Such a condition has never previously been derived for a honey-bee-mite-virus system. When spontaneous transitions do occur, then [Formula: see text] no longer exists, but we introduce a modification of the analysis that allows us to determine whether (i) the disease remains largely covert or (ii) a substantial outbreak of overt disease occurs.

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.

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

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

Tropilaelaps species identification and viral load evaluation of Tropilaelaps and Varroa mites and their Apis mellifera hosts in Palawan, Philippines.

Apis mellifera pupae and their parasites Tropilaelaps and Varroa destructor were collected from honey bee hives in Palawan, Philippines for species identification of the Tropilaelaps and viral analyses. Genetic analysis identified Tropilaelaps mercedesae infesting A. mellifera on the island. Viral analyses showed that all pupae and their infesting Tropilaelaps or Varroa shared the same Deformed Wing Virus (DWV) variant infections with DWV-B being more prevalent than DWV-A. Pupae infested with either Varroa or Tropilaelaps had higher levels of both DWV variants than uninfested pupae. Vigilance is needed to prevent the spread of Tropilaelaps clareae into Palawan and T. mercedesae and DWV variants from Palawan to other provinces.

Deformed wings of small hive beetle independent of virus infections and mites.

Insect wing deformities can be caused by viruses, mites and other environmental stressors during development. Here we conducted differential diagnostics of deformed wings in small hive beetles, Aethina tumida (SHB). Adult SHB with and without deformed wings from individual and mass reared scenarios were evaluated for the mite Tyrophagus putrescentiae and for deformed wing virus. Viral load and mite number were similar for SHB with deformed wings and unaffected beetles. Because deformed wings were only observed in individually pupating SHB, a humidity challenge most likely explains the observed clinical symptoms. Our observations support the importance of differential diagnostics.

Knock-on community impacts of a novel vector: spillover of emerging DWV-B from Varroa-infested honeybees to wild bumblebees.

Novel transmission routes can directly impact the evolutionary ecology of infectious diseases, with potentially dramatic effect on host populations and knock-on effects on the wider host community. The invasion of Varroa destructor, an ectoparasitic viral vector in Western honeybees, provides a unique opportunity to examine how a novel vector affects disease epidemiology in a host community. This specialist honeybee mite vectors deformed wing virus (DWV), an important re-emerging honeybee pathogen that also infects wild bumblebees. Comparing island honeybee and wild bumblebee populations with and without V. destructor, we show that V. destructor drives DWV prevalence and titre in honeybees and sympatric bumblebees. Viral genotypes are shared across hosts, with the potentially more virulent DWV-B overtaking DWV-A in prevalence in a current epidemic. This demonstrates disease emergence across a host community driven by the acquisition of a specialist novel transmission route in one host, with dramatic community level knock-on effects.

Pathogen shifts in a honeybee predator following the arrival of the Varroa mite.

Emerging infectious diseases (EIDs) are a global threat to honeybees, and spillover from managed bees threaten wider insect populations. Deformed wing virus (DWV), a widespread virus that has become emergent in conjunction with the spread of the mite Varroa destructor, is thought to be partly responsible for global colony losses. The arrival of Varroa in honeybee populations causes a dramatic loss of viral genotypic diversity, favouring a few virulent strains. Here, we investigate DWV spillover in an invasive Hawaiian population of the wasp, Vespula pensylvanica, a honeybee predator and honey-raider. We show that Vespula underwent a parallel loss in DWV variant diversity upon the arrival of Varroa, despite the mite being a honeybee specialist. The observed shift in Vespula DWV and the variant-sharing between Vespula and Apis suggest that these wasps can acquire DWV directly or indirectly from honeybees. Apis prey items collected from Vespula foragers were positive for DWV, indicating predation is a possible route of transmission. We also sought cascading effects of DWV shifts in a broader Vespula pathogen community. We identified concurrent changes in a suite of additional pathogens, as well as shifts in the associations between these pathogens in Vespula. These findings reveal how hidden effects of the Varroa mite can, via spillover, transform the composition of pathogens in interacting species, with potential knock-on effects for entire pathogen communities.

Haemolymph removal by Varroa mite destabilizes the dynamical interaction between immune effectors and virus in bees, as predicted by Volterra's model.

The association between the deformed wing virus and the parasitic mite Varroa destructor has been identified as a major cause of worldwide honeybee colony losses. The mite acts as a vector of the viral pathogen and can trigger its replication in infected bees. However, the mechanistic details underlying this tripartite interaction are still poorly defined, and, particularly, the causes of viral proliferation in mite-infested bees. Here, we develop and test a novel hypothesis that mite feeding destabilizes viral immune control through the removal of both virus and immune effectors, triggering uncontrolled viral replication. Our hypothesis is grounded on the predator-prey theory developed by Volterra, which predicts prey proliferation when both predators and preys are constantly removed from the system. Consistent with this hypothesis, we show that the experimental removal of increasing volumes of haemolymph from individual bees results in increasing viral densities. By contrast, we do not find consistent support for alternative proposed mechanisms of viral expansion via mite immune suppression or within-host viral evolution. Our results suggest that haemolymph removal plays an important role in the enhanced pathogen virulence observed in the presence of feeding Varroa mites. Overall, these results provide a new model for the mechanisms driving pathogen-parasite interactions in bees, which ultimately underpin honeybee health decline and colony losses.

Pollen reverses decreased lifespan, altered nutritional metabolism and suppressed immunity in honey bees (Apis mellifera) treated with antibiotics.

Nutrition is involved in regulating multiple aspects of honey bee biology such as caste, immunity, lifespan, growth and behavioral development. Deformed wing virus (DWV) is a major pathogenic factor which threatens honey bee populations, and its replication is regulated by the nutrition status and immune response of honey bees. The alimentary canal of the honey bee is home to a diverse microbial community that provides essential nutrients and serves to bolster immune responses. However, to what extent gut bacteria affect honey bee nutrition metabolism and immunity with respect to DWV has not been investigated fully. In this study, newly emerged worker bees were subjected to four diets that contained (1) pollen, (2) pollen and antibiotics, (3) neither pollen nor antibiotics or (4) antibiotics alone. The expression level of two nutrition genes target of rapamycin (tor) and insulin like peptide (ilp1), one nutritional marker gene vitellogenin (vg), five major royal jellyprotein genes (mrjp1-5), one antimicrobial peptide regulating gene relish (rel), and DWV virus titer and its replication intermediate, negative RNA strand, were determined by qRT-PCR from the honey bees at 7 days post-antibiotic treatment. Additionally, honey bee head mass and survival rate were measured. We observed that antibiotics decreased the expression of tor and rel, and increased DWV titer and its replication activity. Expression of ilp1, mrjp1-5 and vg, and honey bee head mass were also reduced compared with bees on a pollen diet. Antibiotics also caused a significant drop in survivorship, which could be rescued by addition of pollen to the diet. Of importance, pollen could partially rescue the loss of vg and mrjp2 while also increasing the head mass of antibiotic-treated bees. Our results illuminate the roles of bacteria in honey bee nutrition, metabolism and immunity, which confer the ability to inhibit virus replication, extend honey bee lifespan and improve overall health.