The wild solitary bees Andrena vaga, Anthophora plumipes, Colletes cunicularius, and Osmia cornuta microbiota are host specific and dominated by endosymbionts and environmental microorganisms.

We characterized the microbial communities of the crop, midgut, hindgut, and ovaries of the wild solitary bees Andrena vaga, Anthophora plumipes, Colletes cunicularius, and Osmia cornuta through 16S rRNA gene and ITS2 amplicon sequencing and a large-scale isolation campaign. The bacterial communities of these bees were dominated by endosymbionts of the genera Wolbachia and Spiroplasma. Bacterial and yeast genera representing the remaining predominant taxa were linked to an environmental origin. While only a single sampling site was examined for Andrena vaga, Anthophora plumipes, and Colletes cunicularius, and two sampling sites for Osmia cornuta, the microbiota appeared to be host specific: bacterial, but not fungal, communities generally differed between the analyzed bee species, gut compartments and ovaries. This may suggest a selective process determined by floral and host traits. Many of the gut symbionts identified in the present study are characterized by metabolic versatility. Whether they exert similar functionalities within the bee gut and thus functional redundancy remains to be elucidated.

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.

Linking remote sensing data to the estimation of pollination services in agroecosystems.

Wild bees are key providers of pollination services in agroecosystems. The abundance of these pollinators and the ecosystem services they provide rely on supporting resources in the landscape. Spatially explicit models that quantify wild bee abundance and pollination services in food crops are built on the foundations of foraging and nesting resources. This dependence limits model implementation as land-cover maps and pollination experts capable of evaluating habitat resource quality are scarce. This study presents a novel approach to assessing crop pollination services using remote sensing data (RSD) as an alternative to the more conventional use of land-cover data and local expertise on spatially explicit models. We used landscape characteristics derived from remote sensors to qualify nesting resources in the landscape and to evaluate the delivery of pollination services by mining bees (Andrena spp.) in 30 fruit orchards located in the Flemish region of Belgium. For this study, we selected mining bees for their importance as local pollinators and underground nesting behavior. We compared the estimated pollination services derived from RSD with those derived from the conventional qualification of nesting resources. We did not observe significant differences (p = 0.68) in the variation in mining bee activity predicted by the two spatial models. Estimated pollination services derived from RSD and conventional characterizations explained 69% and 72% of the total variation, respectively. These results confirmed that RSD can deliver nesting suitability characterizations sufficient for estimating pollination services. This research also illustrates the importance of nesting resources and landscape characteristics when estimating pollination services delivered by insects like mining bees. Our results support the development of holistic agroenvironmental policies that rely on modern tools like remote sensors and promote pollinators by considering nesting resources.

Managed bumble bees acquire parasites from their foraging environment: A case study on parasite spillback.

The use of commercially reared bumble bees in agricultural environments has been recognized as a potential threat to wild pollinators due to competition, genetic contamination, and most notably, disease transmission. Higher parasite prevalence near greenhouses where managed bumble bees are used has been linked to parasite spillover from managed to wild bees. However, pathogen transmission is not unidirectional, and can also flow from wild to managed bees. These newly infected managed bees can subsequently re-infect (other) wild bees, in a process known as spillback, which is an alternative explanation for the increased parasite prevalence near greenhouses. Reducing parasite prevalence in managed bees is key to controlling host-parasite dynamics in cases of spillover; in spillback, producing managed bees that are resilient to infection is important. Here we establish that the managed bumble bee Bombus terrestris can acquire parasites from their foraging environment, which is the major infection route for Apicystis spp. and Crithidia spp., but not for Nosema spp.. Managed B. terrestris were found to have a higher prevalence of Crithdia and a higher load of Apicystis than local wild conspecifics, showing that for these parasites, spillback is a possible risk scenario.

Establishment of wildflower fields in poor quality landscapes enhances micro-parasite prevalence in wild bumble bees.

The current worldwide pollinator decline is caused by the interplay of different drivers. Several strategies have been undertaken to counteract or halt this decline, one of which is the implementation of wildflower fields. These supplementary flowers provide extra food resources and have proven their success in increasing pollinator biodiversity and abundance. Yet such landscape alterations could also alter the host-pathogen dynamics of pollinators, which could affect the populations. In this study, we investigated the influence of sown wildflower fields on the prevalence of micro-parasites and viruses in the wild bumble bee Bombus pascuorum, one of the most abundant bumble bee species in Europe and the Netherlands. We found that the effect of sown wildflower fields on micro-parasite prevalence is affected by the composition of the surrounding landscape and the size of the flower field. The prevalence of micro-parasites increases with increasing size of sown wildflower fields in landscapes with few semi-natural landscape elements. This effect was not observed in landscapes with a high amount of semi-natural landscape elements. We elaborate on two mechanisms which can support these findings: (1) "transmission hot spots" within the altered flower-networks, which could negatively impact hosts experiencing an increased exposure; (2) improved tolerance of the hosts, withstanding higher parasite populations.

Double-stranded RNA reduces growth rates of the gut parasite Crithidia mellificae.

Parasites of managed bees can disrupt the colony success of the host, but also influence local bee-parasite dynamics, which is regarded as a threat for wild bees. Therapeutic measures have been suggested to improve the health of managed bees, for instance, exploiting the bees' RNA interference (RNAi) pathway to treat against viral pathogens. Gut trypanosomes are an important group of bee parasites in at least two common managed bee species, i.e., managed Apis mellifera and reared Bombus terrestris. In several trypanosomes, RNAi activity is present, while in other associated genes of RNAi, such as Dicer-like (DCL) and Argonaute (AGO), it is lost. Up to date, the ability to exploit the RNAi of gut trypanosomes of bees has remained unexplored. Here, we screened parasite genomes of two honey bee protozoa (Crithidia mellificae and Lotmaria passim) and two bumble bee protozoa (Crithidia bombi and Crithidia expoeki) for the presence of DCL and AGO proteins. For C. mellificae, we constructed a double-stranded RNA (dsRNA) targeting kinetoplastid membrane protein-11 (KMP-11) to test the RNAi potential to kill this parasite. Transfection with KMP-11 dsRNA, but also adding it to the growth medium resulted in small growth reduction of the trypanosome C. mellificae, thereby showing the limited potential to apply dsRNA therapeutics to control trypanosome infection in managed honey bee species. Within bumble bees, there seems to be no application potentials against C. bombi, as we could only retrieve non-functional DCL- and AGO-related genes within the genome of this bumble bee parasite.

Large-scale cultivation of the bumblebee gut microbiota reveals an underestimated bacterial species diversity capable of pathogen inhibition.

A total of 1940 isolates from gut samples of 60 bumblebees representing Bombus pascuorum, Bombus terrestris, Bombus lucorum and Bombus lapidarius was collected and identified through state-of the-art taxonomic methods. The bacterial species diversity in these Bombus species exceeded that suggested by phylotype analysis through 16S rRNA amplicon sequencing, and revealed that B. pascuorum and B. terrestris had a unique microbiota composition, each. Representatives of most phylotypes reported earlier and detected in the present study were effectively isolated, and included several novel bacterial taxa and species reported for the first time in the bumblebee gut. Isolates were screened in pectin degradation assays and growth inhibition assays against the honeybee pathogens Paenibacillus larvae, Melissococcus plutonius and Ascosphaera apis and the bumblebee parasite Crithidia bombi. While inhibitory activity against each of these pathogens was observed, only one single culture was able to degrade pectin and polygalacturonic acid in vitro. The availability of accurately identified microbial isolates will facilitate future evaluation of the functional potential of the bumblebee gut microbiota.

Infection with the multi-host micro-parasite Apicystis bombi (Apicomplexa: Neogregarinorida) decreases survival of the solitary bee Osmia bicornis.

The micro-parasites present in the Apidae pollinator community are mostly multi-host pathogens. To study the impact of these multi-host pathogens on the Apidae pollinator community, as a driver of wild bee decline, pathological studies are needed for different hosts. Yet data on the pathogenesis of these pathogens for different genera and species is scarce or lacking, especially for solitary bee species. In this study, we monitored the effect of the multi-host micro-parasite Apicystis bombi on the lifespan of Osmia bicornis, a common solitary bee. Our results show that A. bombi infection significantly reduces O. bicornis survival.

Systemic Israeli acute paralysis virus (IAPV) infection in bumblebees (Bombus terrestris) through feeding and injection.

Israeli acute paralysis virus (IAPV) can cause a systemic infection, resulting in mortality in both Apis and Bombus spp. bees. However, little is known about the virus infection dynamics within bee tissues. Here, we established systemic IAPV infections in reared bumblebee Bombus terrestris workers through feeding and injection and investigated the mortality, tissue tropism and viral localization. Injection of approximately 500 IAPV (IAPVinj stock) particles resulted in acute infection, viral loads within tissues that were relatively stable from bee to bee, and a distinctive tissue tropism, making this method suitable for studying systemic IAPV infection in bumblebees. Feeding with approximately 1 × 106 particles of the same virus stock did not result in systemic infection. A high-concentration stock of IAPV (IAPVfed stock) allowed us to feed bumblebees with approximately 1 × 109 viral particles, which induced both chronic and acute infection. We also observed a higher variability in viral titers within tissues and less clear tissue tropism during systemic infection, making feeding with IAPVfed stock less optimal for studying IAPV systemic infection. Strikingly, both infection methods and stocks with different viral loads gave a similar viral localization pattern in the brain and midgut of bumblebees with an acute infection. The implications of these findings in the study of the local immunity in bees and barriers to oral transmission are discussed. Our data provide useful information on the establishment of a systemic viral infection in bees.

Israeli acute paralysis virus associated paralysis symptoms, viral tissue distribution and Dicer-2 induction in bumblebee workers (Bombus terrestris).

Although it is known that Israeli acute paralysis virus (IAPV) can cause bee mortality, the symptoms of paralysis and the distribution of the virus in different body tissues and their potential to respond with an increase of the siRNA antiviral immune system have not been studied. In this project we worked with Bombus terrestris, which is one of the most numerous bumblebee species in Europe and an important pollinator for wild flowers and many crops in agriculture. Besides the classic symptoms of paralysis and trembling prior to death, we report a new IAPV-related symptom, crippled/immobilized forelegs. Reverse-transcriptase quantitative PCR showed that IAPV accumulates in different body tissues (midgut, fat body, brain and ovary). The highest levels of IAPV were observed in the fat body. With fluorescence in situ hybridization (FISH) we detected IAPV in the Kenyon cells of mushroom bodies and neuropils from both antennal and optic lobes of the brain in IAPV-infected workers. Finally, we observed an induction of Dicer-2, a core gene of the RNAi antiviral immune response, in the IAPV-infected tissues of B. terrestris workers. According to our results, tissue tropism and the induction strength of Dicer-2 could not be correlated with virus-related paralysis symptoms.

Apibacter mensalis sp. nov.: a rare member of the bumblebee gut microbiota.

Isolates LMG 28357T (=R-53146T) and LMG 28623 were obtained from gut samples of Bombus lapidarius bumblebees caught in Ghent, Belgium. They had identical 16S rRNA gene sequences which were 95.7 % identical to that of Apibacter adventoris wkB301T, a member of the family Flavobacteriaceae. Both isolates had highly similar matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) MS and randomly amplified polymorphic DNA (RAPD) profiles. A draft genome sequence was obtained for strain LMG 28357T (Gold ID Gp0108260); its DNA G+C content was 30.4%, which is within the range reported for members of the family Flavobacteriaceae (27 to 56 mol%) and which is similar to that of the type strain of A. adventoris (29.0 mol%). Whole-cell fatty acid methyl ester analysis of strain LMG 28357T revealed many branched-chain fatty acids, a typical characteristic of bacteria of the family Flavobacteriaceae and a profile that was similar to that reported for A. adventoris wkB301T. MK6 was the major respiratory quinone, again conforming to bacteria of the family Flavobacteriaceae. The isolates LMG 28357T and LMG 28623 could be distinguished from A. adventoris strains through their oxidase activity. On the basis of phylogenetic, genotypic and phenotypic data, we propose to classify both isolates as representatives of a novel species of the genus Apibacter, Apibacter mensalis sp. nov., with LMG 28357T (=DSM 100903T=R-53146T) as the type strain.

The effects of single and mixed infections of Apicystis bombi and deformed wing virus in Bombus terrestris.

Many pollinators are currently suffering from declines, diminishing their gene pool and increasing their vulnerability to parasites. Recently, an increasing diversity of parasites has been recorded in bumblebees, yet for many, knowledge of their virulence and hence the risk their presence poses, is lacking. The deformed wing virus (DWV), known to be ubiquitous in honey bees, has now been detected in bumblebees. In addition, the neogregarine Apicystis bombi has been discovered to be more prevalent than previously thought. Here, we assess for the first time the lethal and sublethal effects of these parasites during single and mixed infections of worker bumblebees (Bombus terrestris). Fifteen days after experimental exposure, 22% of bees exposed to A. bombi, 50% of bees exposed to DWV and 86% of bees exposed to both parasites had died. Bumblebees that had ingested A. bombi had increased sucrose sensitivity (SS) and a lower lipid:body size ratio than control bees. While dual infected bumblebees showed no increase in SS. Overall, we find that A. bombi exhibits both lethal and sublethal effects. DWV causes lethal effect and may reduce the sub lethal effects imposed by A. bombi. The results show that both parasites have significant, negative effects on bumblebee health, making them potentially of conservation concern.

Varroa destructor Macula-like virus, Lake Sinai virus and other new RNA viruses in wild bumblebee hosts (Bombus pascuorum, Bombus lapidarius and Bombus pratorum).

Pollinators such as bumblebees (Bombus spp.) are in decline worldwide which poses a threat not only for ecosystem biodiversity but also to human crop production services. One main cause of pollinator decline may be the infection and transmission of diseases including RNA viruses. Recently, new viruses have been discovered in honeybees, but information on the presence of these in wild bumblebees is largely not available. In this study, we investigated the prevalence of new RNA viruses in Bombus species, and can report for the first time Varroa destructor Macula-like virus (VdMLV) and Lake Sinai virus (LSV) infection in multiple wild bumblebee hosts of Bombus pascuorum, Bombus lapidarius and Bombus pratorum. We sampled in 4 locations in Flanders, Belgium. Besides, we confirmed Slow bee paralysis virus (SBPV) in wild bumblebees, but no positive samples were obtained for Big Sioux river virus (BSRV). Secondly, we screened for the influence of apiaries on the prevalence of these viruses. Our results indicated a location effect for the prevalence of VdMLV in Bombus species, with a higher prevalence in the proximity of honeybee apiaries mainly observed in one location. For LSV, the prevalence was not different in the proximity or at a 1.5 km-distance of apiaries, but we reported a different isolate with similarities to LSV-2 and "LSV-clade A" as described by Ravoet et al. (2015), which was detected both in Apis mellifera and Bombus species. In general, our results indicate the existence of a disease pool of new viruses that seems to be associated to a broad range of Apoidae hosts, including multiple Bombus species.

Large pathogen screening reveals first report of Megaselia scalaris (Diptera: Phoridae) parasitizing Apis mellifera intermissa (Hymenoptera: Apidae).

As it is most likely that global warming will also lead to a shift in pollinator-habitats northwards, the study of southern species becomes more and more important. Pathogen screenings in subspecies of Apis mellifera capable of withstanding higher temperatures, provide an insight into future pathogen host interactions. Screenings in different climate regions also provide a global perspective on the prevalence of certain pathogens. In this project, we performed a pathogen screening in Apis mellifera intermissa, a native subspecies of Algeria in northern Africa. Colonies were sampled from different areas in the region of Annaba over a period of two years. Several pathogens were detected, among them Apicystis bombi, Crithidia mellificae, Nosema ceranae, Paenibacillus larvae, Lake Sinai Virus, Sacbrood Virus and Deformed Wing Virus (DWV). Our screening also revealed a phoroid fly, Megaselia scalaris, parasitizing honey bee colonies, which we report here for the first time. In addition, we found DWV to be present in the adult flies and replicating virus in the larval stages of the fly, which could indicate that M. scalaris acts as a vector of DWV.

Bombella intestini gen. nov., sp. nov., an acetic acid bacterium isolated from bumble bee crop.

In the frame of a bumble bee gut microbiota study, acetic acid bacteria (AAB) were isolated using a combination of direct isolation methods and enrichment procedures. MALDI-TOF MS profiling of the isolates and a comparison of these profiles with profiles of established AAB species identified most isolates as Asaia astilbis or as 'Commensalibacter intestini', except for two isolates (R-52486 and LMG 28161(T)) that showed an identical profile. A nearly complete 16S rRNA gene sequence of strain LMG 28161(T) was determined and showed the highest pairwise similarity to Saccharibacter floricola S-877(T) (96.5%), which corresponded with genus level divergence in the family Acetobacteraceae. Isolate LMG 28161(T) was subjected to whole-genome shotgun sequencing; a 16S-23S rRNA internal transcribed spacer (ITS) sequence as well as partial sequences of the housekeeping genes dnaK, groEL and rpoB were extracted for phylogenetic analyses. The obtained data confirmed that this isolate is best classified into a new genus in the family Acetobacteraceae. The DNA G+C content of strain LMG 28161(T) was 54.9 mol%. The fatty acid compositions of isolates R-52486 and LMG 28161(T) were similar to those of established AAB species [with C18:1ω7c (43.1%) as the major component], but the amounts of fatty acids such as C19:0 cyclo ω8c, C14:0 and C14:0 2-OH enabled to differentiate them. The major ubiquinone was Q-10. Both isolates could also be differentiated from the known genera of AAB by means of biochemical characteristics, such as their inability to oxidize ethanol to acetic acid, negligible acid production from melibiose, and notable acid production from d-fructose, sucrose and d-mannitol. In addition, they produced 2-keto-d-gluconate, but not 5-keto-d-gluconate from d-glucose. Therefore, the name Bombella intestini gen nov., sp. nov. is proposed for this new taxon, with LMG 28161(T) ( =DSM 28636(T) =R-52487(T)) as the type strain of the type species.

Novel lactic acid bacteria isolated from the bumble bee gut: Convivina intestini gen. nov., sp. nov., Lactobacillus bombicola sp. nov., and Weissella bombi sp. nov.

Twelve isolates of lactic acid bacteria (LAB) were obtained in the course of a bumble bee gut microbiota study and grouped into four matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry clusters. Comparative 16S rRNA gene sequence analysis revealed that cluster 1 isolates, represented by strain LMG 28288(T), are most closely related to Lactobacillus apis (97.0% sequence similarity to that of L. apis LMG 26964(T)). Cluster 2 isolates represented by strain LMG 28290(T) are most closely related to Weissella hellenica (99.6% sequence similarity to that of W. hellenica LMG 15125(T)). The single cluster 3 and 4 isolates had identical 16S rRNA gene sequences which were 94.8% similar to that of Leuconostoc mesenteroides subsp. mesenteroides LMG 6893(T), their nearest phylogenetic neighbour. A polyphasic taxonomic study additionally including comparative pheS sequence analysis, DNA-DNA hybridization experiments, DNA G+C content analysis, (GTG)5-PCR fingerprinting and a biochemical characterization, demonstrated that cluster 1 isolates represent a novel Lactobacillus species for which we propose the name Lactobacillus bombicola sp. nov. with LMG 28288(T) (= DSM 28793(T)) as the type strain; and that cluster 2 isolates represent a novel Weissella species for which we propose the name Weissella bombi sp. nov. with LMG 28290(T) (= DSM 28794(T)) as the type strain. Cluster 3 and 4 isolates, in contrast, represented a very distinct, novel taxon that could be distinguished from members of the genera Leuconostoc and Fructobacillus, its nearest phylogenetic neighbours, by its cellular morphology, non-fructophilic metabolism and DNA G+C content. We therefore classify both isolates into a novel species representing a novel LAB genus for which the name Convivina intestini gen. nov., sp. nov. is proposed with LMG 28291(T) (= DSM 28795(T)) as the type strain.

Bifidobacterium commune sp. nov. isolated from the bumble bee gut.

Bifidobacteria were isolated from the gut of Bombus lapidarius, Bombus terrestris and Bombus hypnorum bumble bees by direct isolation on modified trypticase phytone yeast extract agar. The MALDI-TOF MS profiles of four isolates (LMG 28292(T), R-53560, R-53124, LMG 28626) were found to be identical and did not cluster with the profiles of established Bifidobacterium species. Analysis of the 16S rRNA gene sequence of strain LMG 28292(T) revealed that LMG 28292(T) is most closely related to the Bifidobacterium bohemicum type strain (96.8%), which was also isolated from bumble bee gut specimens. The hsp60 gene of strain LMG 28292(T) shows 85.8% sequence similarity to that of the B. bohemicum type strain. The (GTG)5-PCR profiles and the hsp60 sequences of all four isolates were indistinguishable; however, three different phenotypes were observed among the four isolates by means of the API 50CHL microtest system. Based on the phylogenetic, genotypic and phenotypic data, we propose to classify the four isolates within the novel species Bifidobacterium commune sp. nov., with LMG 28292(T) (= DSM 28792(T)) as the type strain.

Bee pathogens found in Bombus atratus from Colombia: A case study.

Bombus atratus bumblebees from Colombia that were caught in the wild and from breeding programs were screened for a broad set of bee pathogens. We discovered for the first time Lake Sinai Virus and confirmed the infection by other common viruses. The prevalence of Apicystis bombi, Crithidia bombi and Nosema ceranae was remarkably high. According to other studies the former two could have been co-introduced in South America with exotic bumble bees as Bombus terrestris or Bombus ruderatus. Given the fact that none of these species occur in Colombia, our data puts a new light on the spread of these pathogens over the South American continent.

Analysis of reference gene stability after Israeli acute paralysis virus infection in bumblebees Bombus terrestris.

To date, there are no validated internal reference genes for the normalization of RT-qPCR data from virus infection experiments with pollinating insects. In this study we evaluated the stability of five candidate internal reference genes: elongation factor-1-alpha (ELF1α), peptidylprolyl isomerase A (PPIA), 60S ribosomal protein L23 (RPL23), TATA-binding protein (TBP) and polyubiquitin (UBI), in relation to Israeli acute paralysis virus (IAPV) infection of Bombus terrestris. We investigated the stability of these genes: in whole bodies and individual body parts, as well as in whole bodies collected at different time intervals after infection with IAPV. Our data identified PPIA as the single, most-optimal internal reference gene and the combination of PPAI-RPL23-UBI as a fully-sufficient multiple internal reference genes set for IAPV infection experiments in B. terrestris.

Effect of oral infection with Kashmir bee virus and Israeli acute paralysis virus on bumblebee (Bombus terrestris) reproductive success.

Israeli acute paralysis virus (IAPV) together with Acute bee paralysis virus (ABPV) and Kashmir bee virus (KBV) constitute a complex of closely related dicistroviruses. They are infamous for their high mortality after injection in honeybees. These viruses have also been reported in non-Apis hymenopteran pollinators such as bumblebees, which got infected with IAPV when placed in the same greenhouse with IAPV infected honeybee hives. Here we orally infected Bombus terrestris workers with different doses of either IAPV or KBV viral particles. The success of the infection was established by analysis of the bumblebees after the impact studies: 50days after infection. Doses of 0.5×10(7) and 1×10(7) virus particles per bee were infectious over this period, for IAPV and KBV respectively, while a dose of 0.5×10(6) IAPV particles per bee was not infectious. The impact of virus infection was studied in micro-colonies consisting of 5 bumblebees, one of which becomes a pseudo-queen which proceeds to lay unfertilized (drone) eggs. The impact parameters studied were: the establishment of a laying pseudo-queen, the timing of egg-laying, the number of drones produced, the weight of these drones and worker mortality. In this setup KBV infection resulted in a significant slower colony startup and offspring production, while only the latter can be reported for IAPV. Neither virus increased worker mortality, at the oral doses used. We recommend further studies on how these viruses transmit between different pollinator species. It is also vital to understand how viral prevalence can affect wild bee populations because disturbance of the natural host-virus association may deteriorate the already critically endangered status of many bumblebee species.