Abstract
Plant-pollinator interactions are essential for plant productivity but face growing threats from climate change, including vegetation loss and mismatches in flowering. Yet, the consequences for bee food resources remain poorly understood at continental scales. Here, we analyse 2 500 samples collected by honey bees (Apis mellifera) between May and August 2023 from 310 locations across Europe using ITS2 metabarcoding. We derive climatic response curves of floral resources and assess exceedance risks of interaction loss under projected climate scenarios. Our findings reveal that rising temperatures and reduced precipitation decrease the diversity of foraging resources across Europe, pushing many plants beyond critical limits. When both warming and drying coincide, the potential for resilience through temporal or spatial buffering is strongly constrained. These declines pose serious risks to bee nutrition, ecosystem functioning, and food security. Our study underscores the urgency of mitigating climate change to preserve vital plant-pollinator systems and the services they sustain.
Similar content being viewed by others
Honey bee colony loss linked to parasites, pesticides and extreme weather across the United States
A Geospatial atlas of honey bee forage plants and their distribution patterns in Africa and beyond
Poor air quality raises mortality in honey bees, a concern for all pollinators
Data availability
The raw sequencing data generated in this study have been deposited in the NCBI SRA database under accession code PRJNA1198597. The processed sequencing data (for process see code availability), climate data, and metadata of samples are available at Zenodo 10.5281/zenodo.17578272 [https://doi.org/10.5281/zenodo.17578272]77 and GitHub chiras/HoneyBee-ResistanceResilience [https://github.com/chiras/HoneyBee-ResistanceResilience].
Code availability
The pipeline for processing raw sequencing data for metabarcoding is publicly available at GitHub chiras/metabarcoding_pipeline [https://github.com/chiras/metabarcoding_pipeline]. Code for all downstream analyses is publicly available at Zenodo 10.5281/zenodo.17578272 [https://doi.org/10.5281/zenodo.17578272]77 and GitHub chiras/HoneyBee-ResistanceResilience [https://github.com/chiras/HoneyBee-ResistanceResilience]. All display items presented in the main manuscript and supplementary information can be reproduced from this public data and code.
References
Ollerton, J., Winfree, R. & Tarrant, S. How many flowering plants are pollinated by animals? Oikos 120, 321–326 (2011).
Potts, S. G. et al. IPBES. The Assessment Report of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services on Pollinators, Pollination and Food Production; http://www.ipbes.net/node/44781, (2016).
Klein, A.-M. et al. Importance of pollinators in changing landscapes for world crops. Proc. R. Soc. B. 274, 303–313 (2007).
Maracchi, G., Sirotenko, O. & Bindi, M. Impacts of present and future climate variability on agriculture and forestry in the temperate regions: Europe. Clim. Change 70, 117–135 (2005).
Olesen, J. E. & Bindi, M. Consequences of climate change for European agricultural productivity, land use and policy. Eur. J. Agron. 16, 239–262 (2002).
Tun, W., Yoon, J., Jeon, J.-S. & An, G. Influence of climate change on flowering time. J. Plant Biol. 64, 193–203 (2021).
Chan, W.-P. et al. Climate velocities and species tracking in global mountain regions. Nature 629, 114–120 (2024).
Nagy, L., Kreyling, J., Gellesch, E., Beierkuhnlein, C. & Jentsch, A. Recurring weather extremes alter the flowering phenology of two common temperate shrubs. Int. J. Biometeorol. 57, 579–588 (2013).
Scheper, J. et al. Museum specimens reveal loss of pollen host plants as key factor driving wild bee decline in the Netherlands. Proc. Natl. Acad. Sci. 111, 17552–17557 (2014).
Dominique, Z., Sabine, S., Herbert, Z., Christa, H.-R. & Sophie, K. Changes in the wild bee community (Hymenoptera: Apoidea) over 100 years in relation to land use: a case study in a protected steppe habitat in Eastern Austria. J. Insect Conserv. 27, 625–641 (2023).
Klein, A.-M., Boreux, V., Fornoff, F., Mupepele, A.-C. & Pufal, G. Relevance of wild and managed bees for human well-being. Curr. Opin. Insect Sci. 26, 82–88 (2018).
Burkle, L. A., Marlin, J. C. & Knight, T. M. Plant-pollinator interactions over 120 years: loss of species, co-occurrence, and function. Science 339, 1611–1615 (2013).
Thomas, C. D. et al. Extinction risk from climate change. Nature 427, 145–148 (2004).
Thuiller, W., Lavorel, S., Araújo, M. B., Sykes, M. T. & Prentice, I. C. Climate change threats to plant diversity in Europe. Proc. Natl. Acad. Sci. 102, 8245–8250 (2005).
Bakkenes, M., Eickhout, B. & Alkemade, R. Impacts of different climate stabilisation scenarios on plant species in Europe. Glob. Environ. Change 16, 19–28 (2006).
Copernicus Climate Change Service (C3S), European state of the climate 2023, summary, https://doi.org/10.24381/bs9v-8c66 (2024).
Feehan, J., Harley, M. & van Minnen, J. Climate change in Europe. 1. Impact on terrestrial ecosystems and biodiversity. A review. Agron. Sustain. Dev. 29, 409–421 (2009).
Cunze, S., Heydel, F. & Tackenberg, O. Are plant species able to keep pace with the rapidly changing climate? PLoS One 8, e67909 (2013).
Parreño, M. A. et al. Critical links between biodiversity and health in wild bee conservation. Trends Ecol. Evol. 37, 309–321 (2022).
Danner, N., Keller, A., Härtel, S. & Steffan-Dewenter, I. Honey bee foraging ecology: Season but not landscape diversity shapes the amount and diversity of collected pollen. PLoS One 12, e0183716 (2017).
Osterman, J. et al. Global trends in the number and diversity of managed pollinator species. Agric. Ecosyst. Environ. 322, 107653 (2021).
Sickel, W. et al. Increased efficiency in identifying mixed pollen samples by meta-barcoding with a dual-indexing approach. BMC Ecol. 15, 20 (2015).
Rafferty, N. E. Effects of global change on insect pollinators: multiple drivers lead to novel communities. Curr. Opin. Insect Sci. 23, 22–27 (2017).
Lee, H. & Romero, J. eds. Climate Change 2023: Synthesis Report. Contribution of Working Groups I, II and III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. IPCC, (2023).
Seneviratne, S. I., Donat, M. G., Pitman, A. J., Knutti, R. & Wilby, R. L. Allowable CO2 emissions based on regional and impact-related climate targets. Nature 529, 477–483 (2016).
Nieto, A. et al. European Red List of Bees. Luxembourg: Publication Office of the European Union, Belgium, (2014).
Menzel, A. & Fabian, P. Growing season extended in Europe. Nature 397, 659–659 (1999).
Vautard, R. et al. The European climate under a 2 °C global warming. Environ. Res. Lett. 9, 034006 (2014).
Ferraro, A. J. & Griffiths, H. G. Quantifying the temperature-independent effect of stratospheric aerosol geoengineering on global-mean precipitation in a multi-model ensemble. Environ. Res. Lett. 11, 034012 (2016).
Settele, J., Bishop, J. & Potts, S. G. Climate change impacts on pollination. Nat. Plants 2, 16092 (2016).
Rasmont, P. et al. Climatic risk and distribution atlas of European Bumblebees Biorisk 10 (Pensoft Publishers, (2015).
Poland, T. M. et al. Invasive Species in Forests and Rangelands of the United States: A Comprehensive Science Synthesis for the United States Forest Sector. Springer Nature (2021).
Uzunov, A. et al. Apis florea in Europe: first report of the dwarf honey bee in Malta. J. Apic. Res. 63, 1122–1125 (2024).
Menzel, A. et al. European phenological response to climate change matches the warming pattern. Glob. Chang. Biol. 12, 1969–1976 (2006).
Eckhardt, M., Haider, M., Dorn, S. & Müller, A. Pollen mixing in pollen generalist solitary bees: a possible strategy to complement or mitigate unfavourable pollen properties? J. Anim. Ecol. 83, 588–597 (2014).
Bell, K. L. et al. Plants, pollinators and their interactions under global ecological change: The role of pollen DNA metabarcoding. Mol. Ecol. 32, 6345–6362 (2023).
Goulson, D., Nicholls, E., Botías, C. & Rotheray, E. L. Bee declines driven by combined stress from parasites, pesticides, and lack of flowers. Science 347, 1255957 (2015).
Horn, J., Becher, M. A., Kennedy, P. J., Osborne, J. L. & Grimm, V. Multiple stressors: using the honeybee model BEEHAVE to explore how spatial and temporal forage stress affects colony resilience. Oikos 125, 1001–1016 (2016).
Timberlake, T. P., Vaughan, I. P. & Memmott, J. Phenology of farmland floral resources reveals seasonal gaps in nectar availability for bumblebees. J. Appl. Ecol. 56, 1585–1596 (2019).
Couvillon, M. J., Schürch, R. & Ratnieks, F. L. W. Waggle dance distances as integrative indicators of seasonal foraging challenges. PLoS One 9, e93495 (2014).
Requier, F. et al. Honey bee diet in intensive farmland habitats reveals an unexpectedly high flower richness and a major role of weeds. Ecol. Appl. 25, 881–890 (2015).
Elliott, B. et al. Pollen diets and niche overlap of honey bees and native bees in protected areas. Basic. Appl. Ecol. 50, 169–180 (2021).
Kaluza, B. F. et al. Social bees are fitter in more biodiverse environments. Sci. Rep. 8, 12353 (2018).
Rhodes, C. J. Pollinator decline – an ecological calamity in the making? Sci. Prog. 101, 121–160 (2018).
Menzel, F. & Feldmeyer, B. How does climate change affect social insects? Curr. Opin. Insect Sci. 46, 10–15 (2021).
Quaresma, A. et al. Preservation methods of honey bee-collected pollen are not a source of bias in ITS2 metabarcoding. Environ. Monit. Assess. 193, 785 (2021).
Chen, S. et al. Validation of the ITS2 Region as a Novel DNA Barcode for Identifying Medicinal Plant Species. PLoS One 5, e8613 (2010).
White, J. S. et al. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In PCR Protocols: A Guide to Methods Applications (eds Innis, M. A., Gelfand, D. H., Shinsky, J. J. & White, T. J.) 315–322 Academic Press, (1990).
Leonhardt, S. D., Peters, B. & Keller, A. Do amino and fatty acid profiles of pollen provisions correlate with bacterial microbiomes in the mason bee Osmia bicornis? Philos. Trans. R. Soc. B 377, 20210171 (2022).
Rognes, T., Flouri, T., Nichols, B., Quince, C. & Mahé, F. VSEARCH: a versatile open source tool for metagenomics. PeerJ 4, e2584 (2016).
Edgar, R. C. SINTAX: a simple non-Bayesian taxonomy classifier for 16S and ITS sequences. bioRxiv, 074161 (2016).
Quaresma, A. et al. Semi-automated sequence curation for reliable reference datasets in ITS2 vascular plant DNA (meta-)barcoding. Sci. Data 11, 129 (2024).
University of East Anglia Climatic Research Unit; Harris, I.C.; Jones, P.D.; Osborn, T. (2024): CRU TS4.08: Climatic Research Unit (CRU) Time-Series (TS) version 4.08 of high-resolution gridded data of month-by-month variation in climate (Jan. 1901- Dec. 2023). NERC EDS Centre for Environmental Data Analysis, date of citation. https://catalogue.ceda.ac.uk/uuid/715abce1604a42f396f81db83aeb2a4b.
R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria (2024).
RStudio Team RStudio: integrated development environment for R. RStudio, PBC, Boston, MA. http://www.rstudio.com/ (2024).
McMurdie, P. J. & Holmes, S. phyloseq: An R Package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8, e61217 (2013).
Wickham, H., Vaughan, D. & Girlich, M. tidyr: tidy messy data. R package version 1.3.1 https://CRAN.R-project.org/package=tidyr (2024).
McLaren, M. speedyseq: faster implementations of phyloseq functions. R package version 0.5.3.9021 https://github.com/mikemc/speedyseq, (2024).
Wickham, H., François, R., Henry, L., Müller, K. & Vaughan, D. dplyr: a grammar of data manipulation. R package version 1.1. 2. Computer software (2023).
Garnier, S. et al. viridis(Lite) – colorblind-friendly color maps for R. viridis package version 0.6.5. (2024).
Dormann, C. F., Fründ, J., Blüthgen, N. & Gruber, B. Indices, graphs and null models: analyzing bipartite ecological networks. Open Ecol. J. 2, 7–24 (2009).
Wickham, H. in ggplot2: elegant graphics for data analysis (ed Hadley Wickham) 189–201 Springer International Publishing, 2016).
Xiao, N. ggsci: scientific journal and sci-fi themed color palettes for ‘ggplot2’. R package version 30, https://CRAN.R-project.org/package=ggsci (2024).
Yutani, H. gghighlight: highlight lines and points in ‘ggplot2’. R package version 0.4.1, https://CRAN.R-project.org/package=gghighlight (2023).
Massicotte, P., South, A. rnatural earth: world map data from natural earth. R package version 11, https://CRAN.R-project.org/package=rnaturalearth (2023).
Pebesma, E. J. Simple features for R: standardized support for spatial vector data. R. J. 10, 439 (2018).
Hijmans, R. Geosphere: spherical trigonometry. R package version 1.5-18, https://CRAN.R-project.org/package=geosphere (2022).
Fitzpatrick, M., Mokany, K., Manion, G., Nieto-Lugilde, D., Ferrier, S. gdm: generalized dissimilarity modeling. R package version 1.5.0-9.1, https://CRAN.R-project.org/package=gdm (2022).
Wickham, H., Pedersen, T., Seidel, D. Scales: scale functions for visualization. R package version 10, (2023). https://CRAN.R-project.org/package=scales.
Microsoft Corporation & Weston, S. foreach: provides for each looping construct. R. package version 1, 2 (2022).
Microsoft Corporation & Weston, S. doParallel: foreach parallel adaptor for the ‘parallel’ Package. R. package version 1, 17 (2022).
Microsoft Corporation & Weston, S. doSNOW: foreach parallel adaptor for the ‘snow’ Package. R. package version 1, 20 (2022).
Csárdi, G. & FitzJohn, R. Progress: terminal progress bars. R. Package version 1, 3 (2023).
Pedersen T. patchwork: the composer of plots. R package version 1.3.0.9000, https://github.com/thomasp85/patchwork (2024).
Pimm, S. L., Donohue, I., Montoya, J. M. & Loreau, M. Measuring resilience is essential to understand it. Nat. Sustain. 2, 895–897 (2019).
Pimm, S. L. The complexity and stability of ecosystems. Nature 307, 321–326 (1984).
Keller, A., Quaresma, A. & Pinto, A. M. Data and code for “Honey bee food resources under threat from climate change”. Zenodo, https://doi.org/10.5281/zenodo.17578272 (2025).
Acknowledgements
We thank all Citizen Scientists and National Coordinators from the 27 EU countries for their invaluable contributions to pollen sampling. We also acknowledge Maíra Costa, Maria Celenza, and Maria João Caldeira for their assistance with pollen homogenisation and DNA extractions. The Portuguese Foundation for Science and Technology (FCT) supported A.Q.‘s PhD scholarship (2020.05155.BD; DOI: 10.54499/2020.05155.BD). FCT provided financial support through national funds to CIMO and SusTEC via FCT/MCTES (PIDDAC): CIMO under UIDB/00690/2020 (DOI: 10.54499/UIDB/00690/2020) and UIDP/00690/2020 (DOI: 10.54499/UIDP/00690/2020); and SusTEC, under LA/P/0007/2020 (DOI: 10.54499/LA/P/0007/2020). This work was conducted within the framework of INSIGNIA-EU: Preparatory action for monitoring of environmental pollution using honey bees (European Union service contract 09.200200/2021/864096/SER/ ENV.D.2).
Funding
Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Contributions
A.Q., A.K., and M.A.P. conceived the ideas and designed the methodology. A.K. and A.Q. conducted the analysis. A.Q., A.K., and M.A.P. drafted the manuscript. Particularly A.Q., M.A.P., A.K., and N.C., but all authors contributed to improving the manuscript. J.B., W.B.B., and I.R. modelled the apiaries’ location. R.B., K.G., F.H., O.K., M.A.P., I.R., F.V., N.C. and JvdS designed the pollen sampling methodology, prepared all the materials and manuals for pollen collection by Citizen Scientists, and obtained the INSIGNIA-EU funding. All the authors critically reviewed the manuscript for important intellectual content.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Hajnalka Szentgyörgyi, Ellen Rotheray, Pritam Banerjee, and the other anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Supplementary Information
Peer Review file
Reporting Summary
Rights and permissions
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/.
Reprints and permissions
About this article
Cite this article
Quaresma, A., Baveco, J.M., Brodschneider, R. et al. Honey bee food resources under threat from climate change.
Nat Commun (2025). https://doi.org/10.1038/s41467-025-68085-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-025-68085-6
Source: Ecology - nature.com
