# Honey bees and climate explain viral prevalence in wild bee communities on a continental scale

Niels Piot, Oliver Schweiger, Ivan Meeus, Orlando Yañez, Lars Straub, Laura Villamar-Bouza, Pilar de La Rúa, Laura Jara, Carlos Ruiz, Martin Malmström, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Niels Piot, Oliver Schweiger, Ivan Meeus, Orlando Yañez, Lars Straub, et al. Honey bees and climate explain viral prevalence in wild bee communities on a continental scale. Scientific Reports, 2022, 12 (1), pp.1904. 10.1038/s41598-022-05603-2. hal-03586244

HAL Id: hal-03586244

https://hal.inrae.fr/hal-03586244

Submitted on 23 Feb 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Honey bees and climate explain viral prevalence in wild bee communities on a continental scale

Niels Piot<fnref n="1" />✉, Oliver Schweiger²,³, Ivan Meeus<fnref n="1" />, Orlando Yañez⁴, Lars Straub⁴, Laura Villamar-Bouza⁴,¹⁸, Pilar De la Rúa⁵, Laura Jara<fnref n="1" />,⁵, Carlos Ruiz⁵,¹⁹, Martin Malmstrøm⁶, Sandra Mustafa¹⁶,²⁰, Anders Nielsen⁶,⁷, Marika Mänd⁸, Reet Karise⁸, Ivana Tlak-Gajger⁹, Erkay Özgör¹⁰,¹⁸, Nevin Keskin¹¹, Virginie Diévart¹², Anne Dalmon¹², Anna Gajda¹³, Peter Neumann⁴, Guy Smagghe<fnref n="1" />, Peter Graystock¹⁴, Rita Radzevičiūtė¹⁵, Robert J. Paxton¹⁶✉ & Joachim R. de Miranda¹⁷✉

Viruses are omnipresent, yet the knowledge on drivers of viral prevalence in wild host populations is often limited. Biotic factors, such as sympatric managed host species, as well as abiotic factors, such as climatic variables, are likely to impact viral prevalence. Managed and wild bees, which harbor several multi-host viruses with a mostly fecal-oral between-species transmission route, provide an excellent system with which to test for the impact of biotic and abiotic factors on viral prevalence in wild host populations. Here we show on a continental scale that the prevalence of three broad host viruses: the AKI-complex (Acute bee paralysis virus, Kashmir bee virus and Israeli acute paralysis virus), Deformed wing virus, and Slow bee paralysis virus in wild bee populations (bumble bees and solitary bees) is positively related to viral prevalence of sympatric honey bees as well as being impacted by climatic variables. The former highlights the need for good beekeeping practices, including *Varroa destructor* management to reduce honey bee viral infection and hive placement. Furthermore, we found that viral prevalence in wild bees is at its lowest at the extreme ends of both temperature and precipitation ranges. Under predicted climate change, the frequency of extremes in precipitation and temperature will continue to increase and may hence impact viral prevalence in wild bee communities.

Even though the presence of viral pathogens often goes unnoticed, they form an indispensable facet of ecosystems<fnref n="1" />⁻³. However, when the natural dynamic interactions between hosts and their viral pathogens are disturbed they can have devastating effects on their hosts, often modulated by host shifts, as recently demonstrated by the SARS-CoV-2 pandemic⁴. The interactions between viruses and their hosts are a complex and can

<footnote n="1">¹Department of Plants and Crops, University of Gent, Coupure Links 653, 9000 Gent, Belgium. ²Department Community Ecology, UFZ-Helmholtz Centre for Environmental Research, Theodor Lieser Str. 4, Halle, Germany. ³iDiv, German Centre for Integrative Biodiversity Research, Halle-Jena-Leipzig, Leipzig, Germany. ⁴Institute of Bee Health, Vetsuisse Faculty, University of Bern, Bern, Switzerland. ⁵Departamento de Zoología y Antropología Física, Facultad de Veterinaria, Universidad de Murcia, Campus de Espinardo, Murcia, Spain. ⁶Department of Biosciences, Center for Ecological and Evolutionary Synthesis (CEES), University of Oslo, Oslo, Norway. ⁷Department of Landscape and Biodiversity, Norwegian Institute of Bioeconomy Research (NIBIO), Ås, Norway. ⁸Estonian University of Life Sciences, Kreutzwaldi 1, Tartu, Estonia. ⁹Department for Biology and Pathology of Fish and Bees, Faculty of Veterinary Medicine University of Zagreb, Zagreb, Croatia. ¹⁰Department of Molecular Biology and Genetics, Cyprus International University, Mersin 10, Nicosia, Turkey. ¹¹Department of Biology, Hacettepe University, Beytepe, Ankara, Turkey. ¹²INRAE, Unité Abeilles et Environnement, Avignon, France. ¹³Laboratory of Bee Diseases, Institute of Veterinary Medicine, Warsaw University of Life Sciences, Warsaw, Poland. ¹⁴Department of Life Sciences, Imperial College London, Silwood Park Campus, Ascot, UK. ¹⁵Molecular Evolution and Animal Systematics, Institute of Biology, University of Leipzig, Talstraße 33, Leipzig, Germany. ¹⁶Institute for Biology, Martin Luther University Halle-Wittenberg, 06120 Halle (Saale), Germany. ¹⁷Department of Ecology, Swedish University of Agricultural Sciences, 750 07 Uppsala, Sweden. ¹⁸European Food Safety Authority (EFSA), Parma, Italy. ¹⁹Departamento de Biología Animal, Edafología y Geología, Facultad de Ciencias, Universidad de La Laguna, La Laguna, Tenerife, Spain. ²⁰Department of Zoology and Entomology, University of Pretoria, Pretoria 0002, South Africa. ²email: niels.piot@ugent.be; robert.paxton@zoologie.uni-halle.de; joachim.de.miranda@slu.se</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Map with the geographical distribution of the 12 sample sites (indicated with black circles) across Europe with color representing the climatic gradients (mean temperature of the warmest month [°C], as an example). Map created using R (version 4.0.4)⁵¹.</figcaption></figure>

be affected by biotic factors, such as contact between wild and domesticated animals or humans⁵,⁶, as well as abiotic factors, such as temperature or precipitation⁷⁻¹³. Hence, to fully understand what drives viral prevalence and what triggers its negative effects on host populations, a thorough understanding of the role of biotic and abiotic factors is crucial for both wildlife as well as human welfare.

The virosphere of bees is very diverse, yet most knowledge on bee viruses and pathogens in general originates from studies of the managed Western honey bee, *Apis mellifera*¹⁴. Data for managed honey bees show that the currently documented viruses appear to have a global distribution¹⁴,¹⁵. In this managed species the viral landscape is highly impacted by the presence of the ectoparasitic mite, *Varroa destructor*, an effective vector of several viruses found in honey bees¹⁶,¹⁷. In the last decade, sparked by reported global declines of wild bees, research focus has shifted to viruses in wild bees¹⁸⁻²⁰. As is true for the majority of pathogens²¹,²², most viruses in bees can be detected in multiple hosts¹⁶,²³,²⁴. Furthermore, co-infections with multiple pathogens are common, where different pathogens can impact one another inside the host²⁰,²⁵⁻²⁸. Moreover, the bee microbiome may modulate within host pathogen dynamics²⁹⁻³¹. Various studies have also highlighted the impact of agrochemical exposure of bees on their interaction with pathogens (reviewed in ³²).

One omnipresent factor impacting the whole bee community and its interactions with pathogens is climate. Albeit, several studies having shown that climatic variables can impact both the honey bee as a host as well as the interactions with its (viral) pathogens³³⁻³⁸, these studies are restricted to local scales and no study has investigated these impacts at large geographical scales. To understand the impact of viral pathogens on bees, we need to understand that a multitude of pathogens interact with different bee species within the pollinator community and how these interactions are influenced by climatic conditions across large spatial scales. The influence of climatic conditions can either be direct, affecting a host's lifestyle, or indirect, by affecting pathogen transmission. Since most bee pathogens are transmitted via the fecal-oral route, enabling inter- and intra-species transmission via shared flowers¹⁶,³⁹⁻⁴², climatic variables, such as UV-exposure, temperature and precipitation could potentially influence pathogen survival on flowers and thereby their transmission⁴⁰. Furthermore, climatic variables affect vegetation phenology⁴³, flower attractiveness⁴⁴ and diversity⁴⁵, which may alter the transmission network via flowers²³, as well as quality, quantity of floral resources⁴⁵⁻⁴⁷, which may impact host immunity⁴⁸,⁴⁹ and consequently pathogen susceptibility and transmission.

Here, we performed a pan-European assessment of the virosphere of 12 bee communities, each consisting of three bee groups: bumble bees, solitary bee species and sympatric managed honey bees, on a continental scale and across different climatic zones (Fig. 1). Using an AICc-based multi-model inference approach⁵⁰, we related viral prevalence, i.e. fraction of specimens with positive virus detection per bee group, in wild bees to that of managed honey bees, climatic conditions, i.e. temperature and precipitation, and vegetation phenology and tested for their interactions, particularly addressing the following two hypotheses:

*H1* Viral prevalence in wild bees is entirely determined by abiotic climatic conditions.

*H2* Viral prevalence in wild bees is related to both biotic factors, namely viral prevalence in sympatric managed honey bees, as well as abiotic climatic conditions.
## Results

In total, we screened 1227 bee specimens (495 honey bees, 476 bumble bees, 256 solitary bees) for the following three virus species: *Deformed wing virus* (DWV), *Slow bee paralysis virus* (SBPV), AKI-virus complex (i.e. *Acute bee paralysis virus* (ABPV), *Kashmir bee virus* (KBV) and *Israeli acute paralysis virus* (IAPV)) and calculated viral prevalence for each of the three bee groups. Across Europe (see Fig. 1) and all three bee groups, DWV was the most prevalent virus, while all viruses showed considerable geographic variation in their prevalence (Supplementary Fig. 1). An exclusive impact of climatic conditions or length of the vegetation period (i.e., approximation for temporal extent of flower availability) on the observed geographic variation in virus prevalence (H1) was not supported. While virus prevalence in honey bees was not related to environmental conditions (see Supplementary Table 1), there was a strong positive relation with viral prevalence in wild bees (Table 1; column 'Apis' present each model, with a positive coefficient), providing support for H2, i.e. a strong relationship between viral prevalence in honey bees and wild bees. The relationship with viral prevalence in honey bees was similar for both bumble bees and solitary bees (there was no interaction effect of species group). However, the relation was different for the three viruses (see Fig. 2). While the slope of this relationship was steep for the AKI-virus complex and SBPV, it was shallower for DWV, i.e., only higher prevalence of DWV in honey bees was reflected by an increase in DWV prevalence in wild bees (Fig. 2). In addition to the strong relation with viral prevalence in honey bees, wild bee viral prevalence depended also on climatic and vegetation conditions (Table 1), further supporting H2. The viral prevalence in wild bees generally declined across the year (Table 1; column 'Phen' indicated by a negative coefficient in all models). We also found virus species-specific responses to the temperature of the warmest quarter. This impact was negligible for DWV and SBPV, while the AKI-virus complex prevalence decreased with the temperature of the warmest quarter (Supplementary Fig. 3). The impact of precipitation of the warmest and driest quarter (Fig. 3A,B) and temperature of the driest quarter (Fig. 3C) was hump-shaped. Viral prevalence in wild bees was particularly low at both high and low levels of precipitation and temperature but without differences among the viruses.

The impacts of vegetation period or environmental conditions on the relation between viral prevalence in managed honey bees and wild bees were only weakly supported ($2 < \Delta AICc < 4$; Table 1: column 'Apis: Env'). Of all environmental conditions, only temperature of the warmest quarter had an impact on the relation in viral prevalence between honey bees and wild bees (Supplementary Fig. 2a). Further, we found that the relation between honey bee viral prevalence and viral prevalence in wild bees was less steep in areas with shorter vegetation periods, while it increased with the length of the vegetation period (Supplementary Fig. 2b).

## Discussion

Within this study we specifically address viral prevalence in wild bees with respect to viral prevalence of sympatric honey bees as a biotic factor, as well as temperature and precipitation as climatic abiotic factors. We find that wild bee viral prevalence is affected by both the biotic factor as well as the abiotic factors, supporting our second hypothesis. We find only a weakly supported interaction effect between honey bee viral prevalence and the climatic factors on wild bee viral prevalence (see Table 1 column 'Apis: Env'; the models with an interaction term all have $\Delta AICc > 2$).

The viral prevalence of multi-host viruses in the bee community is, amongst others, determined by the presence and density of different host species. As bee hosts are all somehow connected through their shared floral resources and viral transmission can be multidirectional for multi-host viruses in the bee community, we cannot demonstrate the directionality of viral transmission in our study. Nevertheless, we find a strong link between the viral prevalence in wild bees and that of sympatric honey bees. This result confirms earlier studies that demonstrate that viral prevalence in wild bees is related to viral prevalence of sympatric honey bees¹⁹,²⁰,⁵³. While these former studies were restricted to local⁵³ and national scales¹⁹,²⁰, we expand these relationships to a continental scale.

Interestingly, the effects of managed honey bee viral prevalence on wild bee viral prevalence did not differ between bumble bees and solitary bees, despite their different lifestyles and social organization, thereby underlining the importance of interspecific pathogen transmission, e.g. via shared floral resources¹⁶,³⁹,⁴²,⁵⁴,⁵⁵. Even though some solitary bee species are dietary specialists, visiting flowers of only a few plant species, the presence of generalist species such as honey bees connect many species in one network module, increasing the likelihood of viral transmission to all bee species present in the same environment⁵⁶⁻⁶⁰. However, we found that the relation between viral prevalence in honey bees and wild bees differs between viruses. There is a tight relationship for both AKI-virus complex and SBPV, with increasing wild bee prevalence when honey bee prevalence increases. Yet, although a similar relation was found for DWV, this relation only becomes apparent when viral prevalence in honey bees was high (see Fig. 2).

In contrast to AKI-virus complex and SBPV, DWV is a virus which is predominantly found in honey bees²⁰,⁵⁹,⁶¹. This discrepancy in viral prevalence is also reflected in our results, where DWV is always more prevalent in honey bees at our different sampling sites compared to the other viruses, namely the AKI-virus complex and SBPV. Experimental infection assays have shown that all three viruses can infect both wild bees as well as honey bees¹⁹,⁶²⁻⁶⁵. However, in contrast to wild bees, honey bees face the additional pressure of the parasitic mite *Varroa destructor*, which besides its virus-vectoring capacity can also weaken colonies, increasing the impact of viral infections⁶⁶,⁶⁷. A recent study highlighted the role of *V. destructor* in the transmission dynamics of DWV⁵⁵. Incorporation of within-honey bee colony transmission, mediated by *V. destructor* vectoring and bee-to-bee transfer into a transmission model could explain the higher DWV prevalence in managed honey bees compared to wild bumble bees observed in the field⁵⁵. These results may also explain why we only observed DWV in wild bees when its prevalence in honey bees was high and the environmental viral load is likely to have been high.
<table><thead><tr><th>Apis</th><th>Phen</th><th>Env</th><th>Env²</th><th>Virus</th><th>Apis:Env</th><th>Apis:Virus</th><th>Env:Virus</th><th>AICc</th><th>Delta</th><th>Weight</th><th>R²m</th><th>R²c</th></tr></thead><tbody><tr><td colspan="13" style="text-align: left; font-weight: bold;">Precipitation warmest quarters</td></tr><tr><td>1.29</td><td>-0.78</td><td></td><td>-1.04</td><td></td><td></td><td></td><td></td><td>228.94</td><td>0.00</td><td>0.65</td><td>0.33</td><td>0.67</td></tr><tr><td>1.30</td><td>-0.84</td><td>0.15</td><td>-0.98</td><td></td><td></td><td></td><td></td><td>231.60</td><td>2.66</td><td>0.17</td><td>0.35</td><td>0.67</td></tr><tr><td>1.45</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>232.84</td><td>3.90</td><td>0.09</td><td>0.27</td><td>0.58</td></tr><tr><td>1.53</td><td>-0.41</td><td></td><td></td><td></td><td></td><td></td><td></td><td>232.85</td><td>3.92</td><td>0.09</td><td>0.31</td><td>0.59</td></tr><tr><td colspan="13" style="text-align: left; font-weight: bold;">Temperature warmest quarter</td></tr><tr><td>2.52</td><td>-0.55</td><td>-1.40</td><td></td><td>+</td><td></td><td></td><td>+</td><td>230.89</td><td>0.00</td><td>0.18</td><td>0.51</td><td>0.64</td></tr><tr><td>1.34</td><td>-0.75</td><td>-0.92</td><td></td><td></td><td></td><td></td><td></td><td>231.83</td><td>0.94</td><td>0.11</td><td>0.26</td><td>0.66</td></tr><tr><td>2.68</td><td></td><td>-1.16</td><td></td><td>+</td><td></td><td></td><td>+</td><td>232.15</td><td>1.27</td><td>0.10</td><td>0.53</td><td>0.66</td></tr><tr><td>2.78</td><td>-0.56</td><td>-1.22</td><td>-0.34</td><td>+</td><td></td><td></td><td>+</td><td>232.64</td><td>1.76</td><td>0.07</td><td>0.56</td><td>0.66</td></tr><tr><td>1.45</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>232.84</td><td>1.95</td><td>0.07</td><td>0.27</td><td>0.58</td></tr><tr><td>1.53</td><td>-0.41</td><td></td><td></td><td></td><td></td><td></td><td></td><td>232.85</td><td>1.97</td><td>0.07</td><td>0.31</td><td>0.59</td></tr><tr><td>2.38</td><td>-0.58</td><td>-1.57</td><td></td><td>+</td><td>0.30</td><td></td><td>+</td><td>233.00</td><td>2.11</td><td>0.06</td><td>0.52</td><td>0.64</td></tr><tr><td>3.00</td><td>-0.48</td><td>-1.55</td><td></td><td>+</td><td></td><td>+</td><td>+</td><td>233.17</td><td>2.28</td><td>0.06</td><td>0.56</td><td>0.66</td></tr><tr><td>3.10</td><td></td><td>-1.52</td><td></td><td>+</td><td></td><td>+</td><td>+</td><td>233.32</td><td>2.43</td><td>0.05</td><td>0.55</td><td>0.66</td></tr><tr><td>1.36</td><td>-0.73</td><td>-0.91</td><td></td><td></td><td>0.17</td><td></td><td></td><td>233.59</td><td>2.71</td><td>0.05</td><td>0.28</td><td>0.65</td></tr><tr><td>2.68</td><td>-0.57</td><td>-1.40</td><td>-0.52</td><td>+</td><td>0.43</td><td></td><td>+</td><td>233.83</td><td>2.95</td><td>0.04</td><td>0.58</td><td>0.67</td></tr><tr><td>2.71</td><td></td><td>-1.06</td><td>-0.25</td><td>+</td><td></td><td></td><td>+</td><td>234.32</td><td>3.44</td><td>0.03</td><td>0.54</td><td>0.66</td></tr><tr><td>1.33</td><td>-0.79</td><td>-0.91</td><td>-0.16</td><td></td><td></td><td></td><td></td><td>234.39</td><td>3.50</td><td>0.03</td><td>0.26</td><td>0.67</td></tr><tr><td>2.44</td><td></td><td>-1.25</td><td></td><td>+</td><td>0.25</td><td></td><td>+</td><td>234.55</td><td>3.66</td><td>0.03</td><td>0.51</td><td>0.64</td></tr><tr><td>2.16</td><td>-0.48</td><td></td><td></td><td>+</td><td></td><td></td><td></td><td>234.67</td><td>3.79</td><td>0.03</td><td>0.42</td><td>0.63</td></tr><tr><td>1.37</td><td></td><td>-0.27</td><td></td><td></td><td></td><td></td><td></td><td>234.73</td><td>3.84</td><td>0.03</td><td>0.23</td><td>0.59</td></tr><tr><td colspan="13" style="text-align: left; font-weight: bold;">Precipitation driest quarter</td></tr><tr><td>1.36</td><td>-0.54</td><td></td><td>-1.02</td><td></td><td></td><td></td><td></td><td>231.02</td><td>0.00</td><td>0.30</td><td>0.33</td><td>0.63</td></tr><tr><td>1.34</td><td></td><td>-0.71</td><td></td><td></td><td></td><td></td><td></td><td>232.72</td><td>1.71</td><td>0.13</td><td>0.28</td><td>0.59</td></tr><tr><td>1.45</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>232.84</td><td>1.82</td><td>0.12</td><td>0.27</td><td>0.58</td></tr><tr><td>1.53</td><td>-0.41</td><td></td><td></td><td></td><td></td><td></td><td></td><td>232.85</td><td>1.84</td><td>0.12</td><td>0.31</td><td>0.59</td></tr><tr><td>1.40</td><td>-0.68</td><td>0.24</td><td>-0.90</td><td></td><td></td><td></td><td></td><td>233.46</td><td>2.45</td><td>0.09</td><td>0.36</td><td>0.62</td></tr><tr><td>2.32</td><td>-0.89</td><td>0.60</td><td></td><td>+</td><td></td><td></td><td></td><td>233.65</td><td>2.63</td><td>0.08</td><td>0.50</td><td>0.66</td></tr><tr><td>1.60</td><td>-0.73</td><td>0.47</td><td></td><td></td><td></td><td></td><td></td><td>233.89</td><td>2.87</td><td>0.07</td><td>0.37</td><td>0.61</td></tr><tr><td>1.29</td><td></td><td>-0.32</td><td>-0.98</td><td></td><td></td><td></td><td></td><td>234.35</td><td>3.33</td><td>0.06</td><td>0.27</td><td>0.62</td></tr><tr><td>2.16</td><td>-0.48</td><td></td><td></td><td></td><td>+</td><td></td><td></td><td>234.67</td><td>3.66</td><td>0.05</td><td>0.42</td><td>0.63</td></tr><tr><td colspan="13" style="text-align: left; font-weight: bold;">Temperature driest quarter</td></tr><tr><td>1.41</td><td>-0.75</td><td></td><td>-0.98</td><td></td><td></td><td></td><td></td><td>231.19</td><td>0.00</td><td>0.32</td><td>0.30</td><td>0.62</td></tr><tr><td>1.45</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>232.84</td><td>1.65</td><td>0.14</td><td>0.27</td><td>0.58</td></tr><tr><td>1.53</td><td>-0.41</td><td></td><td></td><td></td><td></td><td></td><td></td><td>232.85</td><td>1.67</td><td>0.14</td><td>0.31</td><td>0.59</td></tr><tr><td>1.43</td><td>-0.81</td><td>0.41</td><td>-1.29</td><td></td><td></td><td></td><td></td><td>233.54</td><td>2.36</td><td>0.10</td><td>0.32</td><td>0.61</td></tr><tr><td>1.45</td><td>-0.45</td><td>-0.50</td><td></td><td></td><td></td><td></td><td></td><td>234.16</td><td>2.97</td><td>0.07</td><td>0.28</td><td>0.61</td></tr><tr><td>1.39</td><td></td><td></td><td>-0.34</td><td></td><td></td><td></td><td></td><td>234.48</td><td>3.30</td><td>0.06</td><td>0.25</td><td>0.58</td></tr><tr><td>1.40</td><td></td><td>-0.35</td><td></td><td></td><td></td><td></td><td></td><td>234.59</td><td>3.41</td><td>0.06</td><td>0.24</td><td>0.59</td></tr><tr><td>1.83</td><td>-0.72</td><td></td><td>-0.79</td><td>+</td><td></td><td></td><td></td><td>234.59</td><td>3.41</td><td>0.06</td><td>0.37</td><td>0.61</td></tr><tr><td>2.16</td><td>-0.48</td><td></td><td></td><td>+</td><td></td><td></td><td></td><td>234.67</td><td>3.49</td><td>0.06</td><td>0.42</td><td>0.63</td></tr><tr><td colspan="13" style="text-align: left; font-weight: bold;">Length of vegetation period</td></tr><tr><td>1.44</td><td>-1.24</td><td>1.12</td><td></td><td></td><td></td><td></td><td></td><td>231.69</td><td>0.00</td><td>0.26</td><td>0.40</td><td>0.67</td></tr><tr><td>1.45</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>232.84</td><td>1.15</td><td>0.15</td><td>0.27</td><td>0.58</td></tr><tr><td>1.53</td><td>-0.41</td><td></td><td></td><td></td><td></td><td></td><td></td><td>232.85</td><td>1.17</td><td>0.15</td><td>0.31</td><td>0.59</td></tr><tr><td>1.38</td><td>-1.39</td><td>0.65</td><td>0.49</td><td></td><td></td><td></td><td></td><td>233.76</td><td>2.07</td><td>0.09</td><td>0.33</td><td>0.66</td></tr><tr><td>1.46</td><td>-1.25</td><td>1.15</td><td></td><td></td><td>-0.04</td><td></td><td></td><td>234.37</td><td>2.68</td><td>0.07</td><td>0.41</td><td>0.67</td></tr><tr><td>2.16</td><td>-0.48</td><td></td><td>+</td><td></td><td></td><td></td><td></td><td>234.67</td><td>2.98</td><td>0.06</td><td>0.42</td><td>0.63</td></tr><tr><td>1.87</td><td>-1.07</td><td>0.86</td><td></td><td>+</td><td></td><td></td><td></td><td>234.88</td><td>3.20</td><td>0.05</td><td>0.44</td><td>0.65</td></tr><tr><td>1.49</td><td></td><td></td><td>-0.11</td><td></td><td></td><td></td><td></td><td>235.09</td><td>3.41</td><td>0.05</td><td>0.29</td><td>0.59</td></tr><tr><td>1.94</td><td>-0.55</td><td></td><td></td><td>+</td><td></td><td>+</td><td></td><td>235.28</td><td>3.59</td><td>0.04</td><td>0.49</td><td>0.64</td></tr><tr><td>1.94</td><td></td><td></td><td></td><td>+</td><td></td><td></td><td></td><td>235.31</td><td>3.62</td><td>0.04</td><td>0.36</td><td>0.60</td></tr><tr><td>1.44</td><td></td><td>0.05</td><td></td><td></td><td></td><td></td><td></td><td>235.34</td><td>3.65</td><td>0.04</td><td>0.26</td><td>0.58</td></tr></tbody></table>
**Table 1.** Multi-model inference results for separate analyses of climatic and vegetation phenological conditions ordered with increasing AICc. *Apis*, viral prevalence in *Apis mellifera*; Phen, vegetation phenology during sampling; Env, linear term of respective environmental variable; Env², quadratic term of respective environmental variable; Virus, virus species; colon indicates interaction terms. Numbers are coefficient estimates of continuous variables; '+' indicates relevance of categorical variables. AICc, Akaike information criterion corrected for small sample sizes; Delta; delta AICc; Weight, Akaike weight; R²m, marginal pseudo-R² (only fixed effects⁵²); R²c, conditional pseudo-R² (fixed and random effects). Models with delta AICc < 2 indicate strong support and are marked with grey background. White background: models with delta AICc between 2 and 4, indicating weaker support.

Although both the AKI-virus complex and SBPV can also be vectored and amplified to lethal levels in honey bee colonies by *V. destructor*, these viruses are more commonly found in wild bees²⁰,⁵⁹,⁶¹,⁶⁸,⁶⁹, suggesting key differences in host susceptibility and competence between these three viruses. The strong relations found in our study between viral prevalence in honey bees and wild bees together with the results from previous studies¹⁹,⁵⁵ highlight the need for good beekeeping practices, which include proactive varroa monitoring and control by beekeepers as well as appropriate hive placement, where the beekeeper takes into account the presence of endangered wild bees, to prevent managed honey bee colonies becoming sources of virus amplification and dissemination to wild bees¹⁹,⁵⁵,⁷⁰.

In spite of the high variation in viral prevalence in managed honey bees across Europe (Supplementary Fig. 1), we found no relation between viral prevalence in managed honey bees and the investigated environmental conditions, temperature and precipitation (Supplementary Table 1). In contrast to most wild bee species, honey bees are highly social insects, who actively control the climate within the hive⁷¹, and can maintain large food reserves to overcome poor foraging conditions; these traits may make their viral dynamics independent from environmental conditions. Moreover, the active management of honey bees may contribute to the absence of an environmental effect on viral prevalence. Beekeeping practices, such as hive placement and water provisioning as well as varroa mite control and nutritional support⁷², could have a large impact on the health status and viral prevalence in colonies, hence they may obscure potentially smaller effects of environmental conditions in modulating viral prevalence.

Contrary to honey bees, we found that viral prevalence in wild bees was related to environmental conditions, whereby intermediate temperatures and intermediate precipitation of the warmest and driest quarters led to the highest viral prevalence in wild bees.

The impact of climatic conditions, in particular those related to heat and drought stress (precipitation in the driest and in the warmest quarter, temperature in the driest quarter), was consistent across all three viruses. For these environmental conditions, viral prevalence in wild bees was highest in the most moderate (intermediate) conditions and lowest at the extreme ends. Here one might argue that these extremes of both precipitation in the warmest and driest quarter as well as temperature in the driest quarter are likely suboptimal for most pollinator species as well as their floral resources. This is either due to a lack of rain-free days for foraging or drought or a combination of drought and heat, which may have a negative impact on the availability and attractiveness of floral resources⁴⁴,⁴⁷. These suboptimal conditions can induce physiological and nutritional stress, which can either result in reduced host populations⁷³,⁷⁴ and hence lower viral transmission and prevalence or result in an increased virulence of viral infections due to stressed hosts⁶³. Malnutrition has been shown to negatively affect a host's defense system⁷⁵, which can increase virus-induced mortality⁷⁶, as sown for Slow bee paralysis virus and DWV infections in bumble bees⁶³,⁷⁷, resulting in lower viral prevalence⁷⁸,⁷⁹. In the light of climate change, which has been shown to impact host-pathogen interaction⁸⁰,⁸¹, we anticipate that these extremes in both temperature and precipitation will increase⁸²,⁸³. This, together with the ongoing loss of suitable naturally occurring floral resources⁸⁴, can further increase the nutritional stress on wild bees, which may amplify pathogen stress⁸⁵, and hence warrants further research. Besides their direct effect on hosts, these environmental conditions may also affect the viruses themselves. Heat, UV-exposure and drought may impact virus viability on the flowers, and increased rainfall may increase wash-off or dilute the virus below its minimal infection dose, leading to lower viral transmission at flowers.

We found only weak support for environmental modulation of the relationship between viral prevalence in honey bees and wild bees. Here the link between viral prevalence in honey bees and wild bees increased with the length of the vegetation period (Supplementary Fig. 2b), which may be attributed to prolonged foraging activity, facilitating transmission across a longer period of time together with the increase in viral prevalence in honey bee colonies as the foraging season progresses¹⁴,⁸⁶.

Overall, our study shows at a continental scale that wild bee viral prevalence is affected by both biotic factors, namely viral prevalence in sympatric honey bees, and abiotic factors, specifically climatic conditions. The role of managed honey bees as an impact on wild bee viral prevalence has been identified before on a smaller scale, yet the link with climate has to our knowledge not been addressed at a continental scale. Temperature and precipitation extremes will likely continue to increase in the coming years due to climate change, and pathogens can exert negative effects when their interaction dynamics change⁸⁵. Understanding pathogen prevalence in an environmental context, especially that of threatened species, can improve conservation strategies and hence deserves more attention.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Relationship between viral prevalence in managed honey bees and wild bees split by virus. Analysis indicated a strong relation between the viral prevalence in honey bees and wild bees, the slope of the relation differs, however, between the viruses. For both SBPV (top panel) and the AKI-complex (bottom left panel) we see a clear increase in viral prevalence in wild bees as the viral prevalence in managed honey bees increases. For DWV (bottom right panel) we see a similar relation, yet only when the viral prevalence in honey bees increases above 50%. Shaded blue indicates the 95% confidential interval.</figcaption></figure>

## Methods

The bee sampling, identification, processing and virus detection protocols were performed as described in Miranda et al.⁸⁷, and described in short here below. Sample collection, identification, processing, and cDNA synthesis was all done by the local partner. The virus assays were conducted by the partners in Belgium, France, Germany, Sweden, and Switzerland.

**Bee collection.** Foraging adults were collected at 12 sites, in 11 countries across a climatic gradient in Europe (Fig. 1). The field sampling was performed according to a standardized sampling protocol⁸⁷, where each partner sampled in their respective country. At each site 30 honey bees and 30 wild bees were collected in the order they were encountered as well an additional 15 specimens of the most common wild bee. The presence of wild honey bees is close to non-existing in Europe, with the exception of a few records. It is assumed that these wild colonies disappeared due to the introduced parasitic mite *Varroa destructor* and its viral vectoring capacity.⁸⁸,⁸⁹ It is therefore reasonable to presume that nearly all the caught honey bees in this study originate from managed colonies and not from wild ones. At each site, all specimens were collected on the same day and in the same flower-rich area (ca. 100 m²) and then stored individually on ice after capture. Upon return to the laboratory, all samples were stored immediately at -80 °C until further analysis. Sampling date varied between sites as bees were collected between April and September on a day with high foraging activity.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Relationship between viral prevalence (all viruses together) in wild bees (y-axis) and mean precipitation of the warmest quarter (A), mean precipitation of the driest quarter (B) and mean temperature of the driest quarter (C). Analysis indicated that the viral prevalence in wild bees was significantly affected by these three climatic variables. Precipitation and temperature means were obtained over a period of more than two decades at a resolution of 1–3 km. Shaded blue indicates the 95% confidential interval.</figcaption></figure>

**RNA extraction from individual bees and cDNA production.** For RNA extraction of the bees, we used the protocol described by de Miranda et al.⁸⁷. For each bee we only used the abdomen, which was dissected sterile and crushed in a TBS-buffer (50 mM TRIS.HCl pH 7.4, 150 mM NaCl). The volume of TBS buffer used differed depending on the size of the abdomen, i.e., 800 µl for bumble bee sized abdomens, 500 µl for honey bee sized abdomens and 200 µl for small sweat bee sized abdomens. The homogenate was centrifuged to spin down the exoskeleton then 100 µl of supernatant was mixed with 350 µl RLT buffer supplemented with 1% beta-mercaptoethanol. RNA extraction was done using the Qiagen Plant RNeasy Plant kit (Qiagen) according to manufacturer's protocol and was eluted with 50 µl (for large bees) or 30 µl (for small bees) nuclease free water. RNA concentrations were measured using NanoDrop and adjusted to 100 ng/µl using nuclease free water. cDNA was prepared using random hexamer primers and 1 µg RNA and cDNA kit (#K1612, ThermoFisher Scientific, Waltham, Massachusetts, United States of America) containing M-MLV reverse transcriptase and RNAse inhibitor in a 20 µl volume according to the manufacturer's protocol.

**Virus screening.** All bees were screened for three common bee virus targets, namely the AKI-virus complex (acute bee paralysis virus, Kashmir bee virus and Israeli acute paralysis virus), deformed wing virus and slow bee paralysis virus.
These viruses were chosen based on their wide distribution, relative high prevalence, and broad host range where they are detected in both wild bees and honey bees, which was a perquisite to test our hypothesis on the impact of honey bees as a biotic factor on wild bee viral prevalence.

Each PCR reaction contained 2 µl of diluted (1/10) cDNA and 18 µl qPCR mixture, containing 0.2 µl of both forward and reverse primers. For primer sequences, see Supplementary Table 2. Although the used primers are broad-range, capable of detecting several viral strains and species within a complex, we are aware that they are not flawless and may have missed some viral strains/variants due to the high mutation and recombination rates of these RNA-viruses.

**Statistical analysis.** Viral prevalence was calculated for honey bees and wild bees, aggregated into two taxonomic groups, bumble bees and solitary bees, as the fraction of specimens scored positive in virus detection. We used a Cq threshold of 35 to define positive viral detection (Cq<35) per species group. As we report the detection of viruses and classify a bee as positive if Cq values are below the set threshold, we are aware that we cannot state that positive bees are truly infected. This would require the detection of replicating virus (e.g. by detecting the negative strand of these positive strand RNA viruses), which we have not done. All analysis were performed using R (version 4.0.4)⁵¹ and the following packages (rgdal, raster, sp, dplyr, Hmisc, glmmTMB, MuMln, ggplot2, effects)⁹⁰⁻⁹⁸.

**Explanatory variables.** We selected an initial set of the following nine bioclimatic variables covering relevant temperature and precipitation conditions: annual mean temperature (BIO1), temperature seasonality (BIO4), mean temperature of driest quarter (BIO9), mean temperature of warmest quarter (BIO10), mean temperature of coldest quarter (BIO11), annual precipitation (BIO12), precipitation seasonality (BIO15), precipitation of driest quarter (BIO17), and precipitation of warmest quarter (BIO18).

Bioclimatic variables were obtained from CHELSA (v1.2) at a resolution of 30 arc seconds (about 1 km) averaged across the years 1979–2013⁹⁹,¹⁰⁰. To analyze the impact of activity period and time of sampling in the respective sampling year, we extracted data on the start and length of the vegetation period for each site as a proxy for the time of flower availability. Calculation of the vegetation period was based on the Normalized Difference Vegetation Index (NDVI) and provided by the Vegetation Index and Phenology Lab at a 3 arc min (about 5.6 km) resolution (https://vip.arizona.edu/). The actual sampling date was related to the start of the vegetation period per site (number of days since start of the vegetation period) and is referred to as ‘sampling phenology’.

All variables were tested for collinearity with a hierarchical cluster analysis based on pairwise Spearman rank correlation coefficients using UPGMA (unweighted pair-group method with arithmetic averages) agglomeration. Only non-collinear variables were retained (Spearman r < 0.7). Selection from clusters of collinear variables was based on ecological relevance (e.g. mean temperature of the warmest quarter was selected instead of annual mean temperature). After selection, the following six variables remained: mean temperature of driest quarter, mean temperature of warmest quarter, precipitation of driest quarter, precipitation of warmest quarter, length of vegetation period, and sampling phenology.

**Model development.** To test the impact of climate and phenology on honey bee viral prevalence, we used a generalized linear mixed effects model (GLMM) with a binomial error structure (prevalence between 0 and 1) and a logit link, weighted by sample size. Climatic variables were included with their linear and quadratic terms. All three viruses were considered in one model to test for differential responses among the viruses, i.e., interaction between virus type and environmental variable. As crossed random effects, we considered virus type to avoid pseudo-replication, assay laboratory to account for potential systematic differences among the analyzing laboratories and an observer term to address overdispersion where necessary.

To test the relationship between viral prevalence in wild bees and honey bees and the, potentially modulating, impact of environmental variables, we used the same approach as described above but added main effects, two-way and three-way interactions of viral prevalence in honey bees, environment, and virus type. Bumble bees and solitary bees were considered in one model. Group-specific responses to viral prevalence in honey bees were initially tested with interaction effects but without environmental variables to avoid model overfitting. No interaction was evident and thus it was excluded from subsequent analyses. We also added bee taxon (bumble bees, solitary bees) as an additional random effect.

To avoid model overfitting, we limited the number of fixed effects and each environmental variable was tested separately while keeping sampling phenology as a covariate to control for potential differences in viral prevalence due to differences in the time of sampling across the year. We used a multi-model inference approach based on the Akaike Information Criterion corrected for small sample sizes (AICc)¹⁰¹ for model simplification. We considered models with a delta AICc lower than 2, towards the best model, having strong support and models with delta AICc between 2 and 4 as having weaker support. We disregarded a set of models for a particular environmental variable if the intercept-only model was within the respective subset (set of models with delta AICc<2 and set of models with delta 2<AICc<4).

**Data availability**

The dataset is included as a supplementary file.

Received: 25 October 2021; Accepted: 23 December 2021
Published online: 03 February 2022
## References

1. Lafferty, K. D., Dobson, A. P. & Kuris, A. M. Parasites dominate food web links. *Proc. Natl. Acad. Sci. USA* **103**, 11211–11216 (2006).

2. Wood, C. L. & Johnson, P. T. A world without parasites: Exploring the hidden ecology of infection. *Front. Ecol. Environ.* **13**, 425–434 (2015).

3. Preston, D. L., Mischler, J. A., Townsend, A. R. & Johnson, P. T. J. Disease ecology meets ecosystem science. *Ecosystems* **19**, 737–748 (2016).

4. Petersen, E. et al. Comparing SARS-CoV-2 with SARS-CoV and influenza pandemics. *Lancet Infect. Dis.* **20**, e238–e244 (2020).

5. Daszak, P. Emerging infectious diseases of wildlife—Threats to biodiversity and human health. *Science* (80-). **287**, 443–449 (2000).

6. Mack, R. The great African cattle plague epidemic of the 1890s. *Trop. Anim. Health Prod.* **2**, 210–219 (1970).

7. Van Hemert, C., Pearce, J. M. & Handel, C. M. Wildlife health in a rapidly changing North: Focus on avian disease. *Front. Ecol. Environ.* **12**, 548–556 (2014).

8. Hueffer, K., O'Hara, T. M. & Follmann, E. H. Adaptation of mammalian host-pathogen interactions in a changing arctic environment. *Acta Vet. Scand.* **53**, 17 (2011).

9. Musgrave, K., Bartlow, A. W. & Fair, J. M. Long-term variation in environmental conditions influences host-parasite fitness. *Ecol. Evol.* **9**, 7688–7703 (2019).

10. Poisot, T., Guéveneux-Julien, C., Fortin, M. J., Gravel, D. & Legendre, P. Hosts, parasites and their interactions respond to different climatic variables. *Glob. Ecol. Biogeogr.* **26**, 942–951 (2017).

11. Elliot, S. L., Blanford, S. & Thomas, M. B. Host-pathogen interactions in a varying environment: Temperature, behavioural fever and fitness. *Proc. R. Soc. Lond. Ser. B Biol. Sci.* **269**, 1599–1607 (2002).

12. Lambertini, C. et al. Biotic and abiotic determinants of *Batrachochytrium dendrobatidis* infections in amphibians of the Brazilian Atlantic Forest. *Fungal Ecol.* **49**, 100995 (2021).

13. Hoarau, A. O. G., Mavingui, P. & Lebarbenchon, C. Coinfections in wildlife: Focus on a neglected aspect of infectious disease epidemiology. *PLOS Pathog.* **16**, e1008790 (2020).

14. Beaurepaire, A. et al. Diversity and global distribution of viruses of the western honey bee, *Apis mellifera*. *Insects* **11**, 239 (2020).

15. Boncristiani, H. et al. World honey bee health: The global distribution of western honey bee (*Apis mellifera* L.) pests and pathogens. *Bee World.* **98**, 1–5. https://doi.org/10.1080/0005772X.2020.1800330 (2021).

16. Yañez, O. et al. Bee viruses: Routes of infection in hymenoptera. *Front. Microbiol.* **11**, 1–22 (2020).

17. Martin, S. J. et al. Global honey bee viral landscape altered by a parasitic mite. *Science* (80-). **336**, 1304–1306 (2012).

18. Ravoet, J. et al. Widespread occurrence of honey bee pathogens in solitary bees. *J. Invertebr. Pathol.* **122**, 55–58 (2014).

19. Fürst, M. A., McMahon, D. P., Osborne, J. L., Paxton, R. J. & Brown, M. J. F. Disease associations between honeybees and bumblebees as a threat to wild pollinators. *Nature* **506**, 364–366 (2014).

20. McMahon, D. P. et al. A sting in the spit: Widespread cross-infection of multiple RNA viruses across wild and managed bees. *J. Anim. Ecol.* **84**, 615–624 (2015).

21. Rigaud, T., Perrot-Minnot, M.-J. & Brown, M. J. F. Parasite and host assemblages: Embracing the reality will improve our knowledge of parasite transmission and virulence. *Proc. R. Soc. B Biol. Sci.* **277**, 3693–3702 (2010).

22. Rodrigues, R. A. L. et al. An anthropocentric view of the virosphere-host relationship. *Front. Microbiol.* **8**, 1–11 (2017).

23. Daughenbaugh, K. F. et al. Metatranscriptome analysis of sympatric bee species identifies bee virus variants and a new virus, Andrena-associated bee virus-1. *Viruses* **13**, 291 (2021).

24. Olgun, T., Everhart, S. E., Anderson, T. & Wu-Smart, J. Comparative analysis of viruses in four bee species collected from agricultural, urban, and natural landscapes. *PLoS ONE* **15**, e0234431 (2020).

25. Carrillo-Tripp, J. et al. In vivo and in vitro infection dynamics of honey bee viruses. *Sci. Rep.* **6**, 22265 (2016).

26. Remnant, E. J., Mather, N., Gillard, T. L., Yagound, B. & Beekman, M. Direct transmission by injection affects competition among RNA viruses in honeybees. *Proc. Biol. Sci.* **286**, 20182452 (2019).

27. Jones, C. M. & Brown, M. J. F. Parasites and genetic diversity in an invasive bumblebee. *J. Anim. Ecol.* https://doi.org/10.1111/1365-2656.12235 (2014).

28. Faurot-Daniels, C. et al. Longitudinal monitoring of honey bee colonies reveals dynamic nature of virus abundance and indicates a negative impact of Lake Sinai virus 2 on colony health. *PLoS ONE* **15**, e0237544 (2020).

29. Engel, P. et al. The bee microbiome: Impact on bee health and model for evolution and ecology of host-microbe interactions. *MBio* **7**, e02164-e2215 (2016).

30. Maebe, K. et al. The holobiont as a key to the adaptation and conservation of wild bees in the anthropocene. *Front. Ecol. Evol.* **9**, 1–5 (2021).

31. McMenamin, A. J., Brutscher, L. M., Glenny, W. & Flenniken, M. L. Abiotic and biotic factors affecting the replication and pathogenicity of bee viruses. *Curr. Opin. Insect Sci.* **16**, 14–21 (2016).

32. Harwood, G. P. & Dolezal, A. G. Pesticide-virus interactions in honey bees: Challenges and opportunities for understanding drivers of bee declines. *Viruses* **12**, 566 (2020).

33. Dalmon, A., Peruzzi, M., Le Conte, Y., Alaux, C. & Pioz, M. Temperature-driven changes in viral loads in the honey bee *Apis mellifera*. *J. Invertebr. Pathol.* **160**, 87–94 (2019).

34. Le Conte, Y. & Navajas, M. Climate change: Impact on honey bee populations and diseases. *Rev. Sci. Tech.* **27**(485–497), 499–510 (2008).

35. Switanek, M., Crailsheim, K., Truhetz, H. & Brodschneider, R. Modelling seasonal effects of temperature and precipitation on honey bee winter mortality in a temperate climate. *Sci. Total Environ.* **579**, 1581–1587 (2017).

36. Flores, J. M. et al. Effect of the climate change on honey bee colonies in a temperate Mediterranean zone assessed through remote hive weight monitoring system in conjunction with exhaustive colonies assessment. *Sci. Total Environ.* **653**, 1111–1119 (2019).

37. Martin-Hernandez, R. et al. Effect of temperature on the biotic potential of honeybee microsporidia. *Appl. Environ. Microbiol.* **75**, 2554–2557 (2009).

38. Appler, R., Frank, S. & Tarpy, D. Within-colony variation in the immunocompetency of managed and feral honey bees (*Apis mellifera* L.) in different urban landscapes. *Insects* **6**, 912–925 (2015).

39. Alger, S. A., Burnham, P. A. & Brody, A. K. Flowers as viral hot spots: Honey bees (*Apis mellifera*) unevenly deposit viruses across plant species. *PLoS ONE* **14**, e0221800 (2019).

40. Figueroa, L. L. et al. Bee pathogen transmission dynamics: Deposition, persistence and acquisition on flowers. *Proc. R. Soc. B Biol. Sci.* **286**, 20190603 (2019).

41. McArt, S. H., Koch, H., Irwin, R. E. & Adler, L. S. Arranging the bouquet of disease: Floral traits and the transmission of plant and animal pathogens. *Ecol. Lett.* **17**, 624–636 (2014).

42. Durrer, S. & Schmid-Hempel, P. Shared use of flowers leads to horizontal pathogen transmission. *Proc. R. Soc. B Biol. Sci.* **258**, 299–302 (1994).

43. Zhang, X., Tarpley, D. & Sullivan, J. T. Diverse responses of vegetation phenology to a warming climate. *Geophys. Res. Lett.* **34**, L19405 (2007).
44. Descamps, C., Jambrek, A., Quinet, M. & Jacquemart, A.-L. Warm temperatures reduce flower attractiveness and bumblebee foraging. *Insects* **12**, 493 (2021).

45. Suggitt, A. J., Lister, D. G. & Thomas, C. D. Widespread effects of climate change on local plant diversity. *Curr. Biol.* **29**, 2905-2911.e2 (2019).

46. Descamps, C., Quinet, M. & Jacquemart, A.-L. Climate change-induced stress reduce quantity and alter composition of nectar and pollen from a bee-pollinated species (*Borago officinalis*, Boraginaceae). *Front. Plant Sci.* **12**, 2264 (2021).

47. Rering, C. C., Franco, J. G., Yeater, K. M. & Mallinger, R. E. Drought stress alters floral volatiles and reduces floral rewards, pollinator activity, and seed set in a global plant. *Ecosphere* **11**, e03254 (2020).

48. Alaux, C., Ducloz, F., Crauser, D. & Le Conte, Y. Diet effects on honeybee immunocompetence. *Biol. Lett.* **6**, 562-565 (2010).

49. Roger, N., Michez, D., Wattiez, R., Sheridan, C. & Vanderplanck, M. Diet effects on bumblebee health. *J. Insect Physiol.* **96**, 128-133 (2017).

50. Burnham, K. P. & Anderson, D. R. *Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach* (Springer, 2002).

51. R Core Team. *R A Lang. Environ. Stat. Comput. R Found. Stat. Comput. Vienna, Austria*. http://www.R-project.org (2021).

52. Nakagawa, S., Johnson, P. C. D. & Schielzeth, H. The coefficient of determination R 2 and intra-class correlation coefficient from generalized linear mixed-effects models revisited and expanded. *J. R. Soc. Interface* **14**, 20170213 (2017).

53. Alger, S. A., Burnham, P. A., Boncristiani, H. F. & Brody, A. K. RNA virus spillover from managed honeybees (*Apis mellifera*) to wild bumblebees (*Bombus* spp.). *PLoS ONE* **14**, e0217822 (2019).

54. Singh, R. et al. RNA viruses in hymenopteran pollinators: evidence of inter-taxa virus transmission via pollen and potential impact on non-apis hymenopteran species. *PLoS ONE* **5**, e14357 (2010).

55. Burnham, P. A. et al. Flowers as dirty doorknobs: Deformed wing virus transmitted between *Apis mellifera* and Bombus impatiens through shared flowers. *J. Appl. Ecol.* https://doi.org/10.1111/1365-2664.13962 (2021).

56. Figueroa, L. L. et al. Landscape simplification shapes pathogen prevalence in plant-pollinator networks. *Ecol. Lett.* **23**, 1212-1222 (2020).

57. Piot, N., Smagghe, G. & Meeus, I. Network centrality as an indicator for pollinator parasite transmission via flowers. *Insects* **11**, 872 (2020).

58. Truitt, L. L., McArt, S. H., Vaughn, A. H. & Ellner, S. P. Trait-based modeling of multihost pathogen transmission: Plant-pollinator networks. *Am. Nat.* **193**, E149-E167 (2019).

59. Dalmon, A. et al. Possible spillover of pathogens between bee communities foraging on the same floral resource. *Insects* **12**, 122 (2021).

60. Rasmussen, C. et al. Evaluating competition for forage plants between honey bees and wild bees in Denmark. *PLoS ONE* **16**, e0250056 (2021).

61. Manley, R. et al. Knock-on community impacts of a novel vector: Spillover of emerging DWV-B from Varroa-infested honeybees to wild bumblebees. *Ecol. Lett.* https://doi.org/10.1111/ele.13323 (2019).

62. Gusachenko, O. N., Woodford, L., Balbirnie-Cumming, K., Ryabov, E. V. & Evans, D. J. Evidence for and against deformed wing virus spillover from honey bees to bumble bees: A reverse genetic analysis. *Sci. Rep.* **10**, 16847 (2020).

63. Manley, R., Boots, M. & Wilfert, L. Condition-dependent virulence of slow bee paralysis virus in *Bombus terrestris*: Are the impacts of honeybee viruses in wild pollinators underestimated? *Oecologia* **184**, 305-315 (2017).

64. Meeus, I., de Miranda, J. R., de Graaf, D. C., Wäckers, F. & Smagghe, G. Effect of oral infection with Kashmir bee virus and Israeli acute paralysis virus on bumblebee (*Bombus terrestris*) reproductive success. *J. Invertebr. Pathol.* **121**, 64-69 (2014).

65. Wang, H., Meeus, I., Piot, N. & Smagghe, G. Systemic Israeli acute paralysis virus (IAPV) infection in bumblebees (*Bombus terrestris*) through feeding and injection. *J. Invertebr. Pathol.* **151**, 158-164 (2018).

66. Genersch, E. & Aubert, M. Emerging and re-emerging viruses of the honey bee (*Apis mellifera* L.). *Veterinary Res.* **41**, 54 (2010).

67. Martin, S. J. The role of Varroa and viral pathogens in the collapse of honeybee colonies: A modelling approach. *J. Appl. Ecol.* **38**, 1082-1093 (2001).

68. Carreck, N. L., Ball, B. V. & Martin, S. J. Honey bee colony collapse and changes in viral prevalence associated with Varroa destructor. *J. Apic. Res.* **49**, 93-94 (2010).

69. de Miranda, J., Cordoni, G. & Budge, G. The Acute bee paralysis virus–Kashmir bee virus–Israeli acute paralysis virus complex. *J. Invertebr. Pathol.* **103**, S30-S47 (2010).

70. de Pereira, K. et al. Managed bumble bees acquire parasites from their foraging environment: A case study on parasite spillback. *J. Invertebr. Pathol.* **182**, 107583 (2021).

71. Jones, J. C. Honey bee nest thermoregulation: Diversity promotes stability. *Science* (80-)., 305, 402–404 (2004).

72. Döke, M. A., Frazier, M. & Grozinger, C. M. Overwintering honey bees: Biology and management. *Curr. Opin. Insect Sci.* **10**, 185-193 (2015).

73. Rasmont, P. & Iserbyt, S. The Bumblebees Scarcity Syndrome: Are heat waves leading to local extinctions of bumblebees (Hymenoptera: Apidae: Bombus)? *Ann. la Société Entomol. Fr.* **48**, 275-280 (2012).

74. Soroye, P., Newbold, T. & Kerr, J. Climate change contributes to widespread declines among bumble bees across continents. *Science* (80-)., 367, 685–688 (2020).

75. Brunner, F. S., Schmid-Hempel, P. & Barribeau, S. M. Protein-poor diet reduces host-specific immune gene expression in *Bombus terrestris*. *Proc. R. Soc. B Biol. Sci.* **281**, 20140128 (2014).

76. Dolezal, A. G. et al. Interacting stressors matter: Diet quality and virus infection in honeybee health. *R. Soc. Open Sci.* **6**, 181803 (2019).

77. Tehel, A., Streicher, T., Tragust, S. & Paxton, R. J. Experimental infection of bumblebees with honeybee-associated viruses: No direct fitness costs but potential future threats to novel wild bee hosts. *R. Soc. Open Sci.* **7**, 200480 (2020).

78. De Castro, F. & Bolker, B. Mechanisms of disease-induced extinction. *Ecol. Lett.* **8**, 117-126 (2004).

79. Swinton, J., Harwood, J., Grenfell, B. T. & Gilligan, C. A. Persistence thresholds for phocine distemper virus infection in harbour seal *Phoca vitulina* metapopulations. *J. Anim. Ecol.* **67**, 54-68 (1998).

80. Gallana, M., Ryser-Degiorgis, M.-P., Wahli, T. & Segner, H. Climate change and infectious diseases of wildlife: Altered interactions between pathogens, vectors and hosts. *Curr. Zool.* **59**, 427-437 (2013).

81. Zhan, J., Ericson, L. & Burdon, J. J. Climate change accelerates local disease extinction rates in a long-term wild host-pathogen association. *Glob. Chang. Biol.* **24**, 3526-3536 (2018).

82. Trenberth, K. E. Changes in precipitation with climate change. *Clim. Res.* **47**, 123-138 (2011).

83. Fischer, E. M. & Knutti, R. Anthropogenic contribution to global occurrence of heavy-precipitation and high-temperature extremes. *Nat. Clim. Change.* **5**, 560–564 (2015).

84. Goulson, D. et al. Bee declines driven by combined stress from parasites, pesticides, and lack of flowers. *Science* **347**, 1255957 (2015).

85. Meeus, I., Pisman, M., Smagghe, G. & Piot, N. Interaction effects of different drivers of wild bee decline and their influence on host-pathogen dynamics. *Curr. Opin. Insect Sci.* **26**, 136-141 (2018).

86. Tentcheva, D. et al. Prevalence and seasonal variations of six bee viruses in *Apis mellifera* L. and Varroa destructor mite populations in France. *Appl. Environ. Microbiol.* **70**, 7185-7191 (2004).
87. de Miranda, J. R. et al. Protocols for assessing the distribution of pathogens in individual Hymenopteran pollinators. *Protoc. Exch.* https://doi.org/10.21203/rs.3.pex-1453/v1 (2021).

88. Potts, S. G. et al. Declines of managed honey bees and beekeepers in Europe. *J. Apic. Res.* **49**(15), 22 (2010).

89. Moritz, R. F. A., Kraus, F. B., Kryger, P. & Crewe, R. M. The size of wild honeybee populations (*Apis mellifera*) and its implications for the conservation of honeybees. *J. Insect Conserv.* **11**, 391–397 (2007).

90. Bivand, R., Keitt, T. & Rowlingson, B. rgdal: Bindings for the 'geospatial' data abstraction library. (2019).

91. Hijmans, R. J. raster: Geographic data analysis and modeling. (2019).

92. Pebesma, E. J. & Bivand, R. S. Classes and methods for spatial data in R. (2005).

93. Wickham, H., François, R., Henry, L. & Müller, K. dplyr: A grammar of data manipulation. (2020).

94. Harrell, F. E. Hmisc: Harrell miscellaneous. (2020).

95. Brooks, M. E. et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. *R J.* **9**, 378–400 (2017).

96. Barton, K. MuMln: Multi-model inference. (2020).

97. Wickham, H. ggplot2: *Elegant Graphics for Data Analysis* (Springer Verlag, 2016).

98. Fox, J. & Weisberg, S. *An R Companion to Applied Regression* (Sage publishing, 2019).

99. Karger, D. N. et al. Climatologies at high resolution for the earth's land surface areas. *Sci. Data* https://doi.org/10.1038/sdata.2017.122 (2017).

100. Karger, D. N. et al. Data from: Climatologies at high resolution for the earth's land surface areas. Dryad Digit. Repos. https://doi.org/10.5061/dryad.kd1d4 (2018).

101. Burnham, K. P., Anderson, D. R. & Huyvaert, K. P. AIC model selection and multimodel inference in behavioral ecology: Some background, observations, and comparisons. *Behav. Ecol. Sociobiol.* **65**, 23–35 (2011).

## Acknowledgements

We thank the referees for valuable and insightful comments which helped improve the clarity of the manuscript. The authors acknowledge the support of the EU COST-Action FA1307 Sustainable pollination in Europe-joint research on bees and other pollinators (Super-B) for providing the platform to initiate this project, as well as internal funders from each country and institute for financing the sample collection and processing costs. We also would like to thank the technicians and students at the collecting-processing labs and the assaying labs for their invaluable help and competence. RJP thanks Kerstin Gößel and Anja Manigk for technical help. JRM would like to thank the excellent technical support provided by MSc student Sofie Tinggren, BMA student Nellie Svedin and BMA Emilia Semberg during the course of this project. GS, NP and IM would like to thank student Tian Tian for her practical assistance and the Research Foundation-Flanders (FWO-Vlaanderen) for funding. PG would like to thank Crystal Frost for technical assistance. AD would like to thank Maxime Thomasson and Elodie Boeuf for their technical support, Bernard E. Vaissière and Laurent Guilbaud for their support in taxonomy for wild bees, and the scientific department SPE of INRAE for funding. RJP thanks the DFG (Pa632/10) for funding.

## Author contributions

All authors participated in the conception of the project. IM and RJP designed the sample collection strategy. J.R.M., I.M., G.S., P.N. and O.Y. designed the initial sample processing methods. J.R.M., I.M., O.Y., G.S., N.P. and P.N. designed the initial pathogen assay methods. R.J.P. and J.R.M. designed the barcoding identification method. All authors contributed to the optimization and validation of the protocols. All authors contributed to the sample collection and processing. J.R.M., I.M., O.Y., S.M., P.N., R.R., A.D., L.J. and V.D. performed the pathogen assays. JRM curated the data. O.S. performed the statistical analyses. N.P. wrote the first draft of the manuscript. All authors contributed to editing the manuscript. All authors read and approved the final manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-05603-2.

**Correspondence** and requests for materials should be addressed to N.P., R.J.P. or J.R.d.M.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022