Ecology

Clutton-Brock, T. & Sheldon, B. C. Individuals and populations: The role of long-term, individual-based studies of animals in ecology and evolutionary biology. Trends Ecol. Evol. 25, 562–573 (2010).PubMed 
Article 

Google Scholar 
Keuling, O., Lauterbach, K., Stier, N. & Roth, M. Hunter feedback of individually marked wild boar Sus scrofa L.: Dispersal and efficiency of hunting in northeastern Germany. Eur. J. Wildl. Res. 56, 159–167 (2010).Article 

Google Scholar 
Trondrud, L. M. et al. Fat storage influences fasting endurance more than body size in an ungulate. Funct. Ecol. 35, 1470–1480 (2021).CAS 
Article 

Google Scholar 
Wilmers, C. C. et al. The golden age of bio-logging: How animal-borne sensors are advancing the frontiers of ecology. Ecology 96, 1741–1753 (2015).PubMed 
Article 

Google Scholar 
Kukalová, M., Gazárková, A. & Adamík, P. Should i stay or should i go? The influence of handling by researchers on den use in an arboreal nocturnal rodent. Ethology 119, 848–859 (2013).Article 

Google Scholar 
Holt, R. D. et al. Estimating duration of short-term acute effects of capture handling and radiomarking. J. Wildl. Manag. 73, 989–995 (2009).Article 

Google Scholar 
Marco, I., Viñas, L., Velarde, R., Pastor, J. & Lavin, S. Effects of capture and transport on blood parameters in free-ranging mouflon (Ovis ammon). J. Zoo Wildl. Med. 28, 428–433 (1997).CAS 
PubMed 

Google Scholar 
Cattet, M., Boulanger, J., Stenhouse, G., Powell, R. A. & Reynolds-Hogland, M. J. An evaluation of long-term capture effects in ursids: Implications for wildlife welfare and research. J. Mammal. 89, 973–990 (2008).Article 

Google Scholar 
Mortensen, R. M. & Rosell, F. Long-term capture and handling effects on body condition, reproduction and survival in a semi-aquatic mammal. Sci. Rep. 10, 1–16 (2020).Article 

Google Scholar 
Soulsbury, C. D. et al. The welfare and ethics of research involving wild animals: A primer. Methods Ecol. Evol. 11, 1164–1181 (2020).Article 

Google Scholar 
Herman, J. P. et al. Regulation of the hypothalamic-pituitary- adrenocortical stress response. Compr. Physiol. 6, 603–621 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Sapolsky, R. M., Romero, L. M. & Munck, A. U. How do glucocorticoids influence stress responses? Integrating permissive, suppressive, stimulatory, and preparative actions. Endocr. Rev. 21, 55–89 (2000).CAS 
PubMed 

Google Scholar 
Sjaastad, V. Ø., Hove, K. & Sand, O. Physiology of Domestic Animals (Scandinavian Veterinary Press, 2016).
Google Scholar 
Omsjø, E. H. et al. Evaluating capture stress and its effects on reproductive success in Svalbard reindeer. Can. J. Zool. 87, 73–85 (2009).Article 

Google Scholar 
Marco, I., Viñas, L., Velarde, R., Pastor, J. & Lavin, S. The stress response to repeated capture in mouflon (Ovis ammon): Physiological, haematological and biochemical parameters. J. Vet. Med. Ser. A Physiol. Pathol. Clin. Med. 45, 243–253 (1998).CAS 
Article 

Google Scholar 
Hattingh, J., Pitts, N. I. & Ganhao, M. F. Immediate response to repeated capture and handling of wild impala. J. Exp. Zool. 248, 109–112 (1988).CAS 
PubMed 
Article 

Google Scholar 
Ortega, A. C. et al. Effectiveness of partial sedation to reduce stress in captured mule deer. J. Wildl. Manag. 84, 1445–1456 (2020).Article 

Google Scholar 
Arnemo, J. M. & Caulkett, N. Stress. In Zoo Animal and Wildlife Anesthesia and Immobilization (eds West, G. et al.) 103–109 (Blackwell Publications, 2007).
Google Scholar 
Sinclair, M. D. A review of the physiological effects of α2-agonists related to the clinical use of medetomidine in small animal practice. Can. Vet. J. 44, 885–897 (2003).CAS 
PubMed 
PubMed Central 

Google Scholar 
Ranheim, B. et al. The effects of medetomidine and its reversal with atipamezole on plasma glucose, cortisol and noradrenaline in cattle and sheep. J. Vet. Pharmacol. Ther. 23, 379–387 (2000).CAS 
PubMed 
Article 

Google Scholar 
Carroll, G. L. et al. Effect of medetomidine and its antagonism with atipamezole on stress-related hormones, metabolites, physiologic responses, sedation, and mechanical threshold in goats. Vet. Anaesth. Analg. 32, 147–157 (2005).CAS 
PubMed 
Article 

Google Scholar 
Rode, K. D. et al. Effects of capturing and collaring on polar bears: finDings from long-term research on the southern Beaufort Sea population. Wildl. Res. 41, 311–322 (2014).Article 

Google Scholar 
Sakamoto, H., Misumi, K., Nakama, M. & Aoki, Y. The effects of xylazine on intrauterine pressure, uterine blood flow, maternal and fetal cardiovascular and pulmonary function in pregnant goats. J. Vet. Med. Sci. 58, 211–217 (1996).CAS 
PubMed 
Article 

Google Scholar 
Katila, T. & Oijala, M. The effect of detomidine (Domosedan) on the maintenance of equine pregnancy and foetal development: ten cases. Equine Vet. J. 20, 323–326 (1988).CAS 
PubMed 
Article 

Google Scholar 
Larsen, D. G. & Gauthier, D. A. Effects of capturing pregnant moose and calves on calf survivorship. J. Wildl. Manag. 53, 564 (1989).Article 

Google Scholar 
Côté, S. D., Festa-Bianchet, M. & Fournier, F. Life-history effects of chemical immobilization and radiocollars on mountain goats. J. Wildl. Manage. 62, 745–752 (1998).Article 

Google Scholar 
DelGiudice, G. D., Mech, L. D., Paul, W. J. & Karns, P. D. Effects on fawn survival of multiple immobilizations of captive pregnant white-tailed deer. J. Wildl. Dis. 22, 245–248 (1986).CAS 
PubMed 
Article 

Google Scholar 
Brivio, F., Grignolio, S., Sica, N., Cerise, S. & Bassano, B. Assessing the impact of capture on wild animals: The case study of chemical immobilisation on alpine ibex. PLoS ONE 10, e0130957 (2015).PubMed 
PubMed Central 
Article 

Google Scholar 
Wingfield, J. C. et al. Ecological bases of hormone-behavior interactions: The ‘emergency life history stage’. Am. Zool. 38, 191–206 (1998).CAS 
Article 

Google Scholar 
Huber, S., Palme, R. & Arnold, W. Effects of season, sex, and sample collection on concentrations of fecal cortisol metabolites in red deer (Cervus elaphus). Gen. Comp. Endocrinol. 130, 48–54 (2003).CAS 
PubMed 
Article 

Google Scholar 
Morellet, N. et al. The effect of capture on ranging behaviour and activity of the European roe deer Capreolus capreolus. Wildlife Biol. 15, 278–287 (2009).Article 

Google Scholar 
Tarlow, E. M. & Blumstein, D. T. Evaluating methods to quantify anthropogenic stressors on wild animals. Appl. Anim. Behav. Sci. 102, 429–451 (2007).Article 

Google Scholar 
Hik, D. S. Does risk of predation influence the cyclic decline of snowshoe hares. Wildl. Res. 22, 115–129 (1995).Article 

Google Scholar 
Ordiz, A. et al. Lasting behavioural responses of brown bears to experimental encounters with humans. J. Appl. Ecol. 50, 306–314 (2013).Article 

Google Scholar 
Dechen Quinn, A. C., Williams, D. M. & Porter, W. F. Postcapture movement rates can inform data-censoring protocols for GPS-collared animals. J. Mammal. 93, 456–463 (2012).Article 

Google Scholar 
Cattet, M. R. L. Falling through the cracks: Shortcomings in the collaboration between biologists and veterinarians and their consequences for wildlife. ILAR J. 54, 33–40 (2013).CAS 
PubMed 
Article 

Google Scholar 
Albon, S. D. et al. Contrasting effects of summer and winter warming on body mass explain population dynamics in a food-limited Arctic herbivore. Glob. Change Biol. 23, 1374–1389 (2017).ADS 
Article 

Google Scholar 
Ovejero, R. et al. Do cortisol and corticosterone play the same role in coping with stressors? Measuring glucocorticoid serum in free-ranging guanacos (Lama guanicoe). J. Exp. Zool. Part A Ecol. Genet. Physiol. 319, 539–547 (2013).CAS 
Article 

Google Scholar 
Bonacic, C., Feber, R. E. & Macdonald, D. W. Capture of the vicuña (Vicugna vicugna) for sustainable use: Animal welfare implications. Biol. Conserv. 129, 543–550 (2006).Article 

Google Scholar 
Romero, L. M. & Beattie, U. K. Common myths of glucocorticoid function in ecology and conservation. J. Exp. Zool. Part A Ecol. Integr. Physiol. 337, 7–14 (2022).CAS 
Article 

Google Scholar 
Sire, J. E. et al. The effect of blood sampling on plasma cortisol in female reindeer (Rangifer tarandus tarandus L). Acta Vet. Scand. 36, 583–587 (1995).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Harlow, H. J., Thorne, E. T., Williams, E. S., Belden, E. L. & Gern, W. A. Adrenal responsiveness in domestic sheep ( Ovis aries ) to acute and chronic stressors as predicted by remote monitoring of cardiac frequency. Can. J. Zool. 65, 2021–2027 (1987).Article 

Google Scholar 
Pottinger, T. G. & Moran, T. A. Differences in plasma cortisol and cortisone dynamics during stress in two strains of rainbow trout (Oncorhynchus mykiss). J. Fish Biol. 43, 121–130 (1993).CAS 
Article 

Google Scholar 
Arnemo, J. M. & Ranheim, B. Effects of medetomidine and atipamezole on serum glucose and cortisol levels in captive reindeer (Rangifer tarandus tarandus). Rangifer 19, 85–89 (1999).Article 

Google Scholar 
Mentaberre, G. et al. Effects of azaperone and haloperidol on the stress response of drive-net captured Iberian ibexes (Capra pyrenaica). Eur. J. Wildl. Res. 56, 757–764 (2010).Article 

Google Scholar 
Northrup, J. M., Anderson, C. R. & Wittemyer, G. Effects of helicopter capture and handling on movement behavior of mule deer. J. Wildl. Manag. 78, 731–738 (2014).Article 

Google Scholar 
Jung, T. S. et al. Short-term effect of helicopter-based capture on movements of a social ungulate. J. Wildl. Manag. 83, 830–837 (2019).Article 

Google Scholar 
Nurmi, H., Laaksonen, S., Raekallio, M. & Hänninen, L. Wintertime pharmacokinetics of intravenously and orally administered meloxicam in semi-domesticated reindeer (Rangifer tarandus tarandus). Vet. Anaesth. Analg. 49, 423–428 (2022).CAS 
PubMed 
Article 

Google Scholar 
Chapple, R. S., English, A. W., Mulley, R. C. & Lepherd, E. E. Haematology and serum biochemistry of captive unsedated chital deer (Axis axis) in Australia. J. Wildl. Dis. 27, 396–406 (1991).CAS 
PubMed 
Article 

Google Scholar 
Brosh, A. Heart rate measurements as an index of energy expenditure and energy balance in ruminants: A review1. J. Anim. Sci. 85, 1213–1227 (2007).CAS 
PubMed 
Article 

Google Scholar 
Suazo, A. A., Delong, A. T., Bard, A. A. & Oddy, D. M. Repeated capture of beach mice (Peromyscus polionotus phasma and P. P. niveiventris) reduces body mass. J. Mammal. 86, 520–523 (2005).Article 

Google Scholar 
Hoyle, S. D., Horsup, A. B., Johnson, C. N., Crossman, D. G. & McCallum, H. Live-trapping of the northern hairy-nosed wombat (Lasiorhinus krefftii): Population-size estimates and effects on individuals. Wildl. Res. 22, 741–755 (1995).Article 

Google Scholar 
Estruelas, N. F. Short- and long-term physiological effects of capture and handling on free-ranging brown bears (Ursus arctos). PhD Thesis. (Inland Norway University of Applied Sciences, 2017).Veiberg, V. et al. Maternal winter body mass and not spring phenology determine annual calf production in an Arctic herbivore. Oikos 126, 980–987 (2017).Article 

Google Scholar 
Loe, L. E. et al. The neglected season: Warmer autumns counteract harsher winters and promote population growth in Arctic reindeer. Glob. Change Biol. 27, 993–1002 (2021).ADS 
CAS 
Article 

Google Scholar 
Larsen, T. S., Nilsson, N. & Blix, A. S. Seasonal changes in lipogenesis and lipolysis in isolated adipocytes from Svalbard and Norwegian reindeer. Acta Physiol. Scand. 123, 97–104 (1985).CAS 
PubMed 
Article 

Google Scholar 
Colman, J. E., Jacobsen, B. W. & Reimers, E. Summer response distances of Svalbard reindeer (Rangifer tarandus platyrhynchus) to provocation by humans on foot. Wildlife Biol. 7, 275–283 (2001).Article 

Google Scholar 
Trondrud, L. M. et al. Determinants of heart rate in Svalbard reindeer reveal mechanisms of seasonal energy management. Philos. Trans. R. Soc. B Biol. Sci. 376, 20200215 (2021).Article 

Google Scholar 
Pigeon, G. et al. Context-dependent fitness costs of reproduction despite stable body mass costs in an Arctic herbivore. J. Anim. Ecol. 91, 61–73 (2022).PubMed 
Article 

Google Scholar 
Peeters, B., Pedersen, Å., Veiberg, V. & Hansen, B. Hunting quotas, selectivity and stochastic population dynamics challenge the management of wild reindeer. Clim. Res. https://doi.org/10.3354/cr01668 (2021).Article 

Google Scholar 
Loe, L. E. et al. Activity pattern of arctic reindeer in a predator-free environment: No need to keep a daily rhythm. Oecologia 152, 617–624 (2007).ADS 
PubMed 
Article 

Google Scholar 
Dahl, S. R. et al. Assay of steroids by liquid chromatography–tandem mass spectrometry in monitoring 21-hydroxylase deficiency. Endocr. Connect. 7, 1542–1550 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Loe, L. E. et al. Testing five hypotheses of sexual segregation in an arctic ungulate. J. Anim. Ecol. 75, 485–496 (2006).PubMed 
Article 

Google Scholar 
Reimers, E., Lund, S. & Ergon, T. Vigilance and fright behaviour in the insular Svalbard reindeer (Rangifer tarandus platyrhynchus). Can. J. Zool. 89, 753–764 (2011).Article 

Google Scholar 
The R Core Team. R: A language and environment for statistical computing (2021).Burnham, K. P. & Anderson, D. R. in Model selection and multimodel inference. A Practical Information-Theoretic Approach. Ecological Modelling (Springer, 2002).Blanchet, F. G., Tikhonov, G. & Norberg, A. HMSC: Hierarchical modelling of species community. R package version 2.2-0 (2019).Ovaskainen, O. et al. How to make more out of community data? A conceptual framework and its implementation as models and software. Ecol. Lett. 20, 561–576 (2017).PubMed 
Article 

Google Scholar 
Legendre, P. & Legendre, L. Numerical Ecology (Elsevier Science BV, 2012).MATH 

Google Scholar 
Diggle, P. J., Heagerty, P., Liang, K.-Y. & Zeger, S. L. Analysis of Longitudinal Data (Oxford University Press, 2013).MATH 

Google Scholar  More

Longhurst, A., Sathyendranath, S., Platt, T. & Caverhill, C. An estimate of global primary production in the ocean from satellite radiometer data. J. Plankton Res. 17, 1245–1271 (1995).
Google Scholar 
Karl, D. M. et al. Seasonal and interannual variability in primary production and particle flux at station ALOHA. Deep Res. Part II Top. Stud. Oceanogr. 43, 539–568 (1996).CAS 

Google Scholar 
Yang, B., Emerson, S. R. & Quay, P. D. The subtropical ocean’s biological carbon pump determined from O2 and DIC/DI13C tracers. Geophys. Res. Lett. 46, 5361–5368 (2019).
Google Scholar 
Nowicki, M., DeVries, T. & Siegel, D. A. Quantifying the carbon export and sequestration pathways of the ocean’s biological carbon pump. Glob. Biogeochem. Cycles 36, 1–22 (2022).
Google Scholar 
Chávez, F. P., Messié, M. & Pennington, J. T. Marine primary production in relation to climate variability and change. Annu. Rev. Mar. Sci. 3, 227–260 (2011).
Google Scholar 
Polovina, J. J., Howell, E. A. & Abecassis, M. Ocean’s least productive waters are expanding. Geophys. Res. Lett. 35, 2–6 (2008).
Google Scholar 
Irwin, A. J. & Oliver, M. J. Are ocean deserts getting larger? Geophys. Res. Lett. 36, 1–5 (2009).
Google Scholar 
Signorini, S. R., Franz, B. A. & McClain, C. R. Chlorophyll variability in the oligotrophic gyres: Mechanisms, seasonality and trends. Front. Mar. Sci. 2, 1–11 (2015).
Google Scholar 
Sarmiento, J. L., Hughes, T. M. C., Stouffer, R. J. & Manabe, S. Simulated response of the ocean carbon cycle to anthropogenic climate warming. Nature 393, 245–249 (1998).CAS 

Google Scholar 
Bopp, L. et al. Potential impact of climate change on marine export production. Glob. Biogeochem. Cycles 15, 81–99 (2001).CAS 

Google Scholar 
Taucher, J. & Oschlies, A. Can we predict the direction of marine primary production change under global warming? Geophys. Res. Lett. 38, 1–6 (2011).
Google Scholar 
Flombaum, P., Wang, W. L., Primeau, F. W. & Martiny, A. C. Global picophytoplankton niche partitioning predicts overall positive response to ocean warming. Nat. Geosci. 13, 116–120 (2020).CAS 

Google Scholar 
Behrenfeld, M. Uncertain future for ocean algae. Nat. Clim. Chang. 1, 33–34 (2011).CAS 

Google Scholar 
Flombaum, P. & Martiny, A. C. Diverse but uncertain responses of picophytoplankton lineages to future climate change. Limnol. Oceanogr. 66, 4171–4181 (2021).
Google Scholar 
Eppley, R. W. Temperature and phytoplankton growth in the sea. Fish. Bull. 10, 1063–1085 (1972).
Google Scholar 
Falkowski, P. G. & Oliver, M. J. Mix and max: how climate selects phytoplankton. Nature. Rev. Microbiol. 5, 813–819 (2007).CAS 

Google Scholar 
van de Waal, D. B. & Litchman, E. Multiple global change stressor effects on phytoplankton nutrient acquisition in a future ocean. Philos. Trans. R. Soc. B Biol. Sci. 375, 1–8 (2020).
Google Scholar 
Kremer, C. T., Thomas, M. K. & Litchman, E. Temperature- and size-scaling of phytoplankton population growth rates: reconciling the Eppley curve and the metabolic theory of ecology. Limnol. Oceanogr. 62, 1658–1670 (2017).
Google Scholar 
Cross, W. F., Hood, J. M., Benstead, J. P., Huryn, A. D. & Nelson, D. Interactions between temperature and nutrients across levels of ecological organization. Glob. Chang. Biol. 21, 1025–1040 (2015).PubMed 

Google Scholar 
Marañón, E., Lorenzo, M. P., Cermeño, P. & Mouriño-Carballido, B. Nutrient limitation suppresses the temperature dependence of phytoplankton metabolic rates. ISME J. 12, 1836–1845 (2018).PubMed 
PubMed Central 

Google Scholar 
Skau, L. F., Andersen, T., Thrane, J.-E. & Hessen, D. O. Growth, stoichiometry and cell size; temperature and nutrient responses in haptophytes. PeerJ 5, e3743 (2017).PubMed 
PubMed Central 

Google Scholar 
Fernández‐González, C. et al. Effects of temperature and nutrient supply on resource allocation, photosynthetic strategy and metabolic rates of Synechococcus sp. J. Phycol. 56, 818–829 (2020).PubMed 

Google Scholar 
O’Connor, M. I., Piehler, M. F., Leech, D. M., Anton, A. & Bruno, J. F. Warming and resource availability shift food web structure and metabolism. PLoS Biol. 7, e1000178 (2009).PubMed 
PubMed Central 

Google Scholar 
Liu, K., Suzuki, K., Chen, B. & Liu, H. Are temperature sensitivities of Prochlorococcus and Synechococcus impacted by nutrient availability in the subtropical northwest Pacific? Limnol. Oceanogr. 66, 639–651 (2020).
Google Scholar 
Hayashida, H., Matear, R. J. & Strutton, P. G. Background nutrient concentration determines phytoplankton bloom response to marine heatwaves. Glob. Chang. Biol. 26, 4800–4811 (2020).PubMed 

Google Scholar 
Davey, M. et al. Nutrient limitation of picophytoplankton photosynthesis and growth in the tropical North Atlantic. Limnol. Oceanogr. 53, 1722–1733 (2008).CAS 

Google Scholar 
Moore, C. M. et al. Processes and patterns of oceanic nutrient limitation. Nat. Geosci. 6, 701–710 (2013).CAS 

Google Scholar 
Browning, T. J. et al. Nutrient co-limitation at the boundary of an oceanic gyre. Nature 551, 242–246 (2017).CAS 
PubMed 

Google Scholar 
Ustick, L. J. et al. Metagenomic analysis reveals global-scale patterns of ocean nutrient limitation. Science 372, 287–291 (2021).CAS 
PubMed 

Google Scholar 
Zubkov, M. V., Sleigh, M. A., Tarran, G. A., Burkill, P. H. & Leakey, R. J. G. Picoplanktonic community structure on an Atlantic transect from 50°N to 50°S. Deep Res. Part I Oceanogr. Res. Pap. 45, 1339–1355 (1998).
Google Scholar 
Marañón, E., Behrenfeld, M. J., González, N., Mouriño, B. & Zubkov, M. V. High variability of primary production in oligotrophic waters of the Atlantic Ocean: Uncoupling from phytoplankton biomass and size structure. Mar. Ecol. Prog. Ser. 257, 1–11 (2003).
Google Scholar 
Marañón, E. Cell size as a key determinant of phytoplankton metabolism and community structure. Annu. Rev. Mar. Sci. 7, 241–264 (2015).
Google Scholar 
Worden, A. Z., Nolan, J. K. & Palenik, B. Assessing the dynamics and ecology of marine picophytoplankton: the importance of the eukaryotic component. Limnol. Oceanogr. 49, 168–179 (2004).CAS 

Google Scholar 
Visintini, N., Martiny, A. C. & Flombaum, P. Prochlorococcus, Synechococcus, and picoeukaryotic phytoplankton abundances in the global ocean. Limnol. Oceanogr. Lett. 6, 207–215 (2021).
Google Scholar 
Chen, B., Liu, H., Huang, B. & Wang, J. Temperature effects on the growth rate of marine picoplankton. Mar. Ecol. Prog. Ser. 505, 37–47 (2014).
Google Scholar 
Stawiarski, B., Buitenhuis, E. T. & Le Quéré, C. The physiological response of picophytoplankton to temperature and its model representation. Front. Mar. Sci. 3, 1–13 (2016).
Google Scholar 
Marañón, E. et al. Unimodal size scaling of phytoplankton growth and the size dependence of nutrient uptake and use. Ecol. Lett. 16, 371–379 (2013).PubMed 

Google Scholar 
Duhamel, S., Kim, E., Sprung, B. & Anderson, O. R. Small pigmented eukaryotes play a major role in carbon cycling in the P-depleted western subtropical North Atlantic, which may be supported by mixotrophy. Limnol. Oceanogr. 64, 2424–2440 (2019).CAS 

Google Scholar 
Berthelot, H. et al. NanoSIMS single cell analyses reveal the contrasting nitrogen sources for small phytoplankton. ISME J. 13, 651–662 (2019).CAS 
PubMed 

Google Scholar 
Berthelot, H., Duhamel, S., L’Helguen, S., Maguer, J. F. & Cassar, N. Inorganic and organic carbon and nitrogen uptake strategies of picoplankton groups in the northwestern Atlantic Ocean. Limnol. Oceanogr. 66, 3682–3696 (2021).CAS 

Google Scholar 
Marañón, E. et al. Degree of oligotrophy controls the response of microbial plankton to Saharan dust. Limnol. Oceanogr. 55, 2339–2352 (2010).
Google Scholar 
Mouriño-Carballido, B. et al. Nutrient supply controls picoplankton community structure during three contrasting seasons in the northwestern Mediterranean Sea. Mar. Ecol. Prog. Ser. 543, 1–19 (2016).
Google Scholar 
Thomas, M. K., Kremer, C. T., Klausmeier, C. A. & Litchman, E. A global pattern of thermal adaptation in marine phytoplankton. Science 338, 1085–1088 (2012).CAS 
PubMed 

Google Scholar 
Doney, S. C. et al. Climate change impacts on marine ecosystems. Annu. Rev. Mar. Sci. 4, 11–37 (2012).
Google Scholar 
Frölicher, T. L., Fischer, E. M. & Gruber, N. Marine heatwaves under global warming. Nature 560, 360–364 (2018).PubMed 

Google Scholar 
Gruber, N., Boyd, P. W., Frölicher, T. L. & Vogt, M. Biogeochemical extremes and compound events in the ocean. Nature 600, 395–407 (2021).CAS 
PubMed 

Google Scholar 
Babin, S. M., Carton, J. A., Dickey, T. D. & Wiggert, J. D. Satellite evidence of hurricane-induced phytoplankton blooms in an oceanic desert. J. Geophys. Res. Oceans 109, 1–21 (2004).
Google Scholar 
Walker, N. D., Leben, R. R. & Balasubramanian, S. Hurricane-forced upwelling and chlorophyll a enhancement within cold-core cyclones in the Gulf of Mexico. Geophys. Res. Lett. 32, 1–5 (2005).
Google Scholar 
Boyd, P. W. et al. Experimental strategies to assess the biological ramifications of multiple drivers of global ocean change—a review. Glob. Chang. Biol. 24, 2239–2261 (2018).PubMed 

Google Scholar 
Mills, M. M., Ridame, C., Davey, M., La Roche, J. & Geider, R. J. Iron and phosphorus co-limit nitrogen fixation in the eastern tropical North Atlantic. Nature 429, 292–294 (2004).CAS 
PubMed 

Google Scholar 
Marañón, E. Phytoplankton growth rates in the Atlantic subtropical gyres. Limnol. Oceanogr. 50, 299–310 (2005).
Google Scholar 
Halsey, K. H. & Jones, B. M. Phytoplankton strategies for photosynthetic energy allocation. Annu. Rev. Mar. Sci. 7, 265–297 (2015).
Google Scholar 
Quevedo, M. & Anadón, R. Protist control of phytoplankton growth in the subtropical north-east Atlantic. Mar. Ecol. Prog. Ser. 221, 29–38 (2001).
Google Scholar 
Schmoker, C., Hernández-León, S. & Calbet, A. Microzooplankton grazing in the oceans: Impacts, data variability, knowledge gaps and future directions. J. Plankton Res. 35, 691–706 (2013).
Google Scholar 
Landry, M. R. & Hassett, R. P. Estimating the grazing impact of marine micro-zooplankton. Mar. Biol. 67, 283–288 (1982).
Google Scholar 
Kiørboe, T. Turbulence, phytoplankton cell size, and the structure of pelagic food webs. Adv. Mar. Biol. 29, 1–72 (1993).
Google Scholar 
Cermeño, P. et al. Marine primary productivity is driven by a selection effect. Front. Mar. Sci. 3, 1–10 (2016).Browning, T. J. et al. Nutrient co-limitation in the subtropical Northwest Pacific. Limnol. Oceanogr. Lett. 7, 52–61 (2022).
Google Scholar 
Klausmeier, C. A., Litchman, E. & Levin, S. A. Phytoplankton growth and stoichiometry under multiple nutrient limitation. Limnol. Oceanogr. 49, 1463–1470 (2004).
Google Scholar 
Behrenfeld, M. J. & Milligan, A. J. Photophysiological expressions of iron stress in phytoplankton. Annu. Rev. Mar. Sci. 5, 217–246 (2013).
Google Scholar 
Geider, R. J. Light and temperature dependence of the carnon to chlorophyll a ratio in microalgae and cyanobacteria: implications for physiology and growth of phytoplankton. N. Phytol. 106, 1–34 (1987).CAS 

Google Scholar 
Maxwell, D. P., Laudenbach, D. E. & Huner, N. P. Redox regulation of light-harvesting complex II and cab mRNA abundance in Dunaliella salina. Plant Physiol. 109, 787–795 (1995).CAS 
PubMed 
PubMed Central 

Google Scholar 
Ye, H. J., Sui, Y., Tang, D. L. & Afanasyev, Y. D. A subsurface chlorophyll a bloom induced by typhoon in the South China Sea. J. Mar. Syst. 128, 138–145 (2013).
Google Scholar 
Zhang, H., He, H., Zhang, W. Z. & Tian, D. Upper ocean response to tropical cyclones: a review. Geosci. Lett. 8, 1–12 (2021).
Google Scholar 
Lin, I. et al. New evidence for enhanced ocean primary production triggered by tropical cyclone. Geophys. Res. Lett. 30, 1–4 (2003).Chai, F. et al. A limited effect of sub-tropical typhoons on phytoplankton dynamics. Biogeosciences 18, 849–859 (2021).
Google Scholar 
Gillooly, J. F., Brown, J. H., West, G. B., Savage, V. M. & Charnov, E. L. Effects of size and temperature on metabolic rate. Science 293, 2248–2251 (2001).CAS 
PubMed 

Google Scholar 
Brown, J. H., Gillooly, J. F., Allen, A. P., Savage, V. M. & West, G. B. Toward a metabolic theory of ecology. Ecology 85, 1771–1789 (2004).
Google Scholar 
Davidson, E. A. & Janssens, I. A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 440, 165–173 (2006).CAS 
PubMed 

Google Scholar 
Somero, G. N. Adaptation of enzymes to temperature: Searching for basic ‘strategies’. Comp. Biochem. Physiol.—B Biochem. Mol. Biol. 139, 321–333 (2004).PubMed 

Google Scholar 
Rose, J. M. & Caron, D. A. Does low temperature constrain the growth rates of heterotrophic protists? Evidence and implications for algal blooms in cold waters. Limnol. Oceanogr. 52, 886–895 (2007).
Google Scholar 
Harvey, B. P., Marshall, K. E., Harley, C. D. G. & Russell, B. D. Predicting responses to marine heatwaves using functional traits. Trends Ecol. Evol. 37, 20–29 (2022).PubMed 

Google Scholar 
Staehr, P. A. & Birkeland, M. J. Temperature acclimation of growth, photosynthesis and respiration in two mesophilic phytoplankton species. Phycologia 45, 648–656 (2006).
Google Scholar 
Morán, X. A. G., Calvo-Díaz, A., Arandia-Gorostidi, N. & Huete-Stauffer, T. M. Temperature sensitivities of microbial plankton net growth rates are seasonally coherent and linked to nutrient availability. Environ. Microbiol. 20, 3798–3810 (2018).PubMed 

Google Scholar 
Courboulès, J. et al. Effects of experimental warming on small phytoplankton, bacteria and viruses in autumn in the Mediterranean coastal Thau Lagoon. Aquat. Ecol. 55, 647–666 (2021).
Google Scholar 
López-Sandoval, D. C., Duarte, C. M. & Agustí, S. Nutrient and temperature constraints on primary production and net phytoplankton growth in a tropical ecosystem. Limnol. Oceanogr. 66, 2923–2935 (2021).
Google Scholar 
Landry, M. R., Selph, K. E., Hood, R. R., Davies, C. H. & Beckley, L. E. Low temperature sensitivity of picophytoplankton P:B ratios and growth rates across a natural 10 °C temperature gradient in the oligotrophic Indian Ocean. Limnol. Oceanogr. Lett. https://doi.org/10.1002/lol2.10224 (2021)Martiny, A. C. et al. Strong latitudinal patterns in the elemental ratios of marine plankton and organic matter. Nat. Geosci. 6, 279–283 (2013).CAS 

Google Scholar 
Fernández-González, C. & Marañón, E. Effect of temperature on the unimodal size scaling of phytoplankton growth. Sci. Rep. 11, 1–9 (2021).
Google Scholar 
Marañón, E. et al. Patterns of phytoplankton size structure and productivity in contrasting open-ocean environments. Mar. Ecol. Prog. Ser. 216, 43–56 (2001).
Google Scholar 
Tarran, G. A., Heywood, J. L. & Zubkov, M. V. Latitudinal changes in the standing stocks of nano- and picoeukaryotic phytoplankton in the Atlantic Ocean. Deep Res. Part II Top. Stud. Oceanogr. 53, 1516–1529 (2006).
Google Scholar 
Hillebrand, H. et al. Cell size as driver and sentinel of phytoplankton community structure and functioning. Funct. Ecol. 1–18 https://doi.org/10.1111/1365-2435.13986 (2021).Partensky, F. & Garczarek, L. Prochlorococcus: advantages and limits of minimalism. Annu. Rev. Mar. Sci. 2, 305–331 (2010).
Google Scholar 
Landry, M. R. et al. Biological response to iron fertilization in the eastern equatorial Pacific (IronEx II). I. Microplankton community abundances and biomass. Mar. Ecol. Prog. Ser. 201, 27–42 (2000).CAS 

Google Scholar 
Morel, A. et al. Examining the consistency of products derived from various ocean color sensors in open ocean (Case 1) waters in the perspective of a multi-sensor approach. Remote Sens. Environ. 111, 69–88 (2007).
Google Scholar 
Fofonoff, N. P. & Millard, R. C. Algorithms for computation of fundamental properties of seawater. UNESCO Tech. Pap. Mar. Sci. 44, 1–53 (1983).
Google Scholar 
Becker, S. et al. GO-SHIP repeat hydrography nutrient manual: the precise and accurate determination of dissolved inorganic nutrients in seawater, using continuous flow analysis methods. Front. Mar. Sci. https://doi.org/10.3389/fmars.2020.581790 (2020).Marañón, E. et al. Resource supply overrides temperature as a controlling factor of marine phytoplankton growth. PLoS ONE 9, 20–23 (2014).
Google Scholar 
Schuback, N. et al. Single-turnover variable chlorophyll fluorescence as a tool for assessing phytoplankton photosynthesis and primary productivity: opportunities, caveats and recommendations. Front. Mar. Sci. 8, 1–24 (2021).Piggott, J. J., Townsend, C. R. & Matthaei, C. D. Reconceptualizing synergism and antagonism among multiple stressors. Ecol. Evol. 5, 1538–1547 (2015).PubMed 
PubMed Central 

Google Scholar  More

Fauna, culture and chronology datasetsA geo-referenced dataset of chronometric dates covering the late MIS 3 (55–30 kyr cal bp) was compiled from the literature (dataset 1). The dataset included 363 radiocarbon, thermoluminescence, optically stimulated luminescence and uranium series dates obtained from 62 archaeological sites and seven palaeontological sites. These chronological determinations were obtained from ten palaeontological levels and 138 archaeological levels. The archaeological levels were culturally attributed to the Mousterian (n = 75), Châtelperronian (n = 6) and Aurignacian (n = 57) technocomplexes. A number of issues can potentially hamper the chronological assessment of Palaeolithic technocomplexes from radiocarbon dates, such as pretreatment protocols that do not remove sufficient contaminants or the quality of the bone collagen extracted. Moreover, discrepancies in cultural attributions or stratigraphic inconsistencies are commonly detected in Palaeolithic archaeology. Information regarding the quality of date determinations and cultural attribution or stratigraphic issues is provided in the Supplementary Information.Our dataset also included the presence of herbivore species recovered from each archaeo-palaeontological site (hereafter referred to as local faunal assemblages (LFAs)), their body masses and their chronology. The mean body mass of both sexes, for each species, was obtained from the PHYLACINE database53 and used in the macroecological modelling approach described below (see ‘Carrying capacity of herbivores’). For visual representation purposes, the herbivore species were grouped into four weight categories: small (500 kg). The chronology of the occurrence of each herbivore species was assumed to be the same as the dated archaeo-palaeontological layer where the species remains were recovered. Thus, to estimate the chronological range of each species in each region, all radiocarbon determinations were calibrated with the IntCal20 calibration curve54 and OxCAL4.2 software55. The BAMs were run to compute the upper and lower chronological boundaries at a CI of 95.4% of each LFA (see ‘Chronological assessment’ for more details). One of the purposes of the current study was to estimate the potential fluctuations in herbivore biomass during the stadial and interstadial periods of the late MIS 3. Accordingly, the time spans of the LFAs were classified into the discrete GS and GI phases provided by Rasmussen et al.51.Geographic settingsThe Iberian Peninsula locates at the southwestern edge of Europe (Fig. 1). It constitutes a large geographic area that exhibits a remarkable diversity of ecosystems, climates and landscapes. Both now and in the past, altitudinal, latitudinal and oceanic gradients affected the conformation of two biogeographical macroregions with different flora and fauna species pools: the Eurosiberian and Mediterranean regions13,46. In the north, along the Pyrenees and Cantabrian strip, the Eurosiberian region is characterized by oceanic influence and mild temperatures in the present day, whereas the Mediterranean region features drier summers and milder winters (Fig. 1). Between the Eurosiberian and Mediterranean regions, there is a transitional area termed Submediterranean or Supramediterranean. Lastly, the Mediterranean region is divided into two distinctive bioclimatic belts: (1) the Thermomediterranean region, located at lower latitudes, with high evapotranspiration rates and affected by its proximity to the coast; and (2) the Mesomediterranean region, with lower temperatures and wetter conditions (Fig. 1).Previous studies have shown that zoocoenosis and phytocenosis differed between these macroregions in the Pleistocene13,46. However, flora and fauna distributions changed during the stadial–interstadial cycles in the Iberian Peninsula, which suggests potential alterations in the boundaries of these biogeographical regions. The modelling approach used in this study to estimate the biomass of primary consumers is dependent on the reconstructed NPP and the herbivore guild structure in each biogeographical region. To test the suitability of the present-day biogeographical demarcations of the Iberian Peninsula during MIS 3, we assessed whether the temporal trends of NPP and the composition of each herbivore palaeocommunity differed between these biogeographical regions during the MUPT.Chouakria and Nagabhusan56 proposed a dissimilarity index to compare time series data by taking into consideration the proximity of values and the temporal correlation of the time series:$${rm{CORT}}(S_1,S_2) = frac{{mathop {sum}nolimits_{i = 1}^{p – 1} {left( {u_{left( {i + 1} right)} – u_i} right)} (v_{(i + 1)} – v_i)}}{{sqrt {mathop {sum}nolimits_{i = 1}^{p – 1} {(u_{(i + 1)} – u_i)^2} } sqrt {mathop {sum}nolimits_{i = 1}^{p – 1} {(v_{(i + 1)} – v)^2} } }}$$
(1)
where S1 and S2 are the time series of data, u and v represent the values of S1 and S2, respectively, and p is the length of values of each time series. CORT(S1, S2) belongs to the interval (−1,1). The value CORT(S1, S2) = 1 indicates that in any observed period (ti, ti+1), the values of the sequence S1 and those of S2 increase or decrease at the same rate, whereas CORT = −1 indicates that when S1 increases, S2 decreases or vice versa. Lastly, CORT(S1, S2) = 0 indicates that the observed trends in S1 are independent of those observed in S2. To complement this approach by considering not only the temporal correlation between each pair of time series but also the proximity between the raw values, these authors proposed an adaptive tuning function defined as follows:$$d{rm{CORT}}left( {S_1,S_2} right) = fleft({{rm{CORT}}left( {S_1,S_2} right)} right)times dleft( {S_1,S_2} right)$$
(2)
where$$fleft( x right) = frac{2}{{1 + exp left( {k,x} right)}},k ge 0$$
(3)
In this study, k was 2, meaning that the behaviour contribution was 76% and the contribution of the proximity between values was 24%57. Hence, f(x) modulates a conventional pairwise raw data distance (d(S1,S2)) according to the observed temporal correlation56. Consequently, dCORT adjusts the degree of similarity between each pair of observations according to the temporal correlation and the proximity between values. This function was used to compare the reconstructed NPP between biogeographical regions during MIS 3 in the Iberian Peninsula. However, two different biogeographical regions could have experienced similar evolutionary trends in their NPP, even though their biota composition was different. Therefore, this analysis was complemented with a JSI to assess whether the reconstructed herbivore species composition in each palaeocommunity differed among biogeographical regions during the late MIS 3. The JSI was based on presence–absence data and was calculated as follows:$${rm{JSI}} = frac{c}{{(a + b + c)}}$$
(4)
where c is the number of shared species in both regions and a and b are the numbers of species that were only present in one of the biogeographical regions. Therefore, the higher the value the more similar the palaeocommunities of both regions were.Chronological assessmentPivotal to any hypothesis of Neanderthal replacement patterns by AMHs is the chronology of that population turnover. To this end, we used three different approaches to provide greater confidence in the results: BAMs, the OLE model and SPD of archaeological assemblages. As detailed below, each of these approaches provides complementary information about the MUPT.First, we built a set of BAMs for the Mousterian, Châtelperronian and Aurignacian technocomplexes in each region during the MIS 3. As stated above, we compiled the available radiocarbon dates for Iberia between 55 and 30 kyr cal bp. However, not all dates or levels were included in the Bayesian chronology models. Radiocarbon determinations obtained from shell remains were incorporated in the dataset (dataset 1); however, the local variation of the reservoir age was unknown from 55 to 30 kyr bp. Because of uncertainties related to marine reservoir offsets, all BAMs that incorporated dates from marine shells were run twice: including and excluding these dates. All of the archaeological levels with cultural attribution issues or stratigraphic inconsistencies were excluded. The Supplementary Note provides a detailed description of the sites, levels and dates excluded and their justification. All BAMs were built for each technocomplex using the OxCAL4.2 software55 and IntCal20 calibration curve54.Bayesian chronology models were built for each archaeological and palaeontological level. Then, the dates associated with each technocomplex were grouped within a single phase to determine each culture’s regional appearance or disappearance. Our interest was not focused on the chronological duration of the Mousterian, Châtelperronian and Aurignacian cultures, but on the probability distribution function of the temporal boundaries of these cultures in each region. Thus, this chronological assessment aims to provide an updated chronological frame for Neanderthal replacement by AMHs in Iberia. For this reason, we did not differentiate between proto- and early Aurignacian cultures, since both are attributed to AMHs.In each BAM, we inserted into the same sequence the radiocarbon dates associated with a given technocomplex within a start and end boundary to bracket each culture, which allowed us to determine the probability distribution function for the beginning and end moment of each cultural phase6. The resolution of all models was set at 20 years. We used a t-type outlier model with an initial 5% probability for each determination, but when more than one radiocarbon date was obtained from the same bone remain, we used an s-type outlier model and the combine function. The thermoluminescence dating likelihoods were included in the models, together with their associated 1σ uncertainty ranges. When dates with low agreement ( More

The impacts of climate change on host–parasite dynamics are particularly complex to predict, as they involve an interplay of both physiological and behavioural factors, from both host and parasite. For example, while warming may increase parasite developmental rates and thus increase transmission, excessive heat may instead exceed thermal limits, leading to higher parasite mortality. Transmission also relates to both the distribution and abundance of host species, which may also shift under changing climates. More

Berenbaum, M. R. & Calla, B. Honey as a functional food for Apis mellifera. Annu. Rev. Entomol. 66, 185–208. https://doi.org/10.1146/annurev-ento-040320-074933 (2021).CAS 
Article 
PubMed 

Google Scholar 
Crane, E. Honey: A Comprehensive Survey (Heinemann, 1975).
Google Scholar 
Park, O. W. The storing and ripening of honey by honeybees. J. Econ. Entomol. 18, 405–410 (1925).Article 

Google Scholar 
Reinhardt, J. F. Ventilating the bee colony to facilitate the honey ripening process. J. Econ. Entomol. 32, 654–660. https://doi.org/10.1093/jee/32.5.654 (1939).Article 

Google Scholar 
Eyer, M., Neumann, P. & Dietemann, V. A look into the cell: Honey storage in honey bees, Apis mellifera. PLoS ONE 11(8), e0161059 (2016).Article 

Google Scholar 
Oertel, E., Fieger, E. A., Williams, V. R. & Andrews, E. A. Inversion of cane sugar in the honey stomach of the bee. J. Econ. Entomol. 44, 487–492 (1951).CAS 
Article 

Google Scholar 
Park, O. W. Studies on the changes in nectar concentration produced by the honeybee, Apis mellifera. Part I. Changes which occur between the flower and the hive. Res. Bull. Iowa Agric. Exp. Station 151, 211–243 (1932).
Google Scholar 
Nicolson, S. W. & Human, H. Bees get a head start on honey production. Biol. Let. 4, 299–301. https://doi.org/10.1098/rsbl.2008.0034 (2008).Article 

Google Scholar 
Nicolson, S. W. & Louw, G. N. Simultaneous measurement of evaporative water loss, oxygen consumption, and thoracic temperature during flight in a carpenter bee. J. Exp. Zool. 222, 287–296 (1982).Article 

Google Scholar 
Schmid-Hempel, P., Kacelnik, A. & Houston, A. I. Honeybees maximize efficiency by not filling their crop. Behav. Ecol. Sociobiol. 17, 61–66 (1985).Article 

Google Scholar 
Kacelnik, A., Houston, A. I. & Schmid-Hempel, P. Central-place foraging in honey bees: The effect of travel time and nectar flow on crop filling. Behav. Ecol. Sociobiol. 19, 19–24. https://doi.org/10.1007/BF00303838 (1986).Article 

Google Scholar 
Wolf, T. J., Schmid-Hempel, P., Ellington, C. P. & Stevenson, R. D. Physiological correlates of foraging efforts in honey-bees: Oxygen consumption and nectar load. Funct. Ecol. 3, 417–424 (1989).Article 

Google Scholar 
Mitchell, D. Thermal efficiency extends distance and variety for honeybee foragers: Analysis of the energetics of nectar collection and desiccation by Apis mellifera. J. R. Soc. Interface 16, 20180879. https://doi.org/10.1098/rsif.2018.0879 (2019).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Corbet, S. A. et al. Native or exotic? Double or single? Evaluating plants for pollinator-friendly gardens. Ann. Bot. 87, 219–232 (2001).Article 

Google Scholar 
Harano, K. & Nakamura, J. Nectar loads as fuel for collecting nectar and pollen in honeybees: Adjustment by sugar concentration. J. Comp. Physiol. A. https://doi.org/10.1007/s00359-016-1088-x (2016).Article 

Google Scholar 
Nicolson, S. W. & van Wyk, B.-E. Nectar sugars in Proteaceae: Patterns and processes. Aust. J. Bot. 46, 489–504 (1998).Article 

Google Scholar 
Corbet, S. A. Nectar sugar content: Estimating standing crop and secretion rate in the field. Apidologie 34, 1–10. https://doi.org/10.1051/apido:2002049 (2003).CAS 
Article 

Google Scholar 
Southwick, E. E. & Pimentel, D. Energy efficiency of honey production by bees. Bioscience 31, 730–732. https://doi.org/10.2307/1308779 (1981).Article 

Google Scholar 
Mitchell, D. Nectar, humidity, honey bees (Apis mellifera) and varroa in summer: A theoretical thermofluid analysis of the fate of water vapour from honey ripening and its implications on the control of Varroa destructor. J. R. Soc. Interface 16, 20190048. https://doi.org/10.1098/rsif.2019.0048 (2019).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Human, H., Nicolson, S. W. & Dietemann, V. Do honeybees, Apis mellifera scutellata, regulate humidity in their nest?. Naturwissenschaften 93, 397–401 (2006).ADS 
CAS 
Article 

Google Scholar 
Ellis, M. B. Homeostasis: Humidity and water relations in honeybee colonies, MSc thesis, University of Pretoria (2008).Ellis, M., Nicolson, S., Crewe, R. & Dietemann, V. Hygropreference and brood care in the honeybee (Apis mellifera). J. Insect Physiol. 54, 1516–1521. https://doi.org/10.1016/j.jinsphys.2008.08.011 (2008).CAS 
Article 
PubMed 

Google Scholar 
Portman, Z. M., Ascher, J. S. & Cariveau, D. P. Nectar concentrating behavior by bees (Hymenoptera: Anthophila). Apidologie 52, 1169–1194. https://doi.org/10.1007/s13592-021-00895-1 (2021).Article 

Google Scholar 
Nicolson, S. W. Water homeostasis in bees, with the emphasis on sociality. J. Exp. Biol. 212, 429–434. https://doi.org/10.1242/jeb.022343 (2009).Article 
PubMed 

Google Scholar 
Pokorny, T., Lunau, K. & Eltz, T. Raising the sugar content – orchid bees overcome the constraints of suction feeding through manipulation of nectar and pollen provisions. PLoS ONE 9(11), e113823. https://doi.org/10.1371/journal.pone.0113823 (2014).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Lindauer, M. The water economy and temperature regulation of the honeybee colony. Bee World 36, 81–92 (1955).Article 

Google Scholar 
Heinrich, B. Mechanisms of body-temperature regulation in honeybees, Apis mellifera. I. Regulation of head temperature. J. Exp. Biol. 85, 61–72 (1980).Article 

Google Scholar 
Cooper, P. D., Schaffer, W. M. & Buchmann, S. L. Temperature regulation of honeybees (Apis mellifera) foraging in the Sonoran desert. J. Exp. Biol. 114, 1–15 (1985).Article 

Google Scholar 
Louw, G. N. & Hadley, N. F. Water economy of the honeybee: A stoichiometric accounting. J. Exp. Zool. 235, 147–150 (1985).Article 

Google Scholar 
Rodney, S. & Purdy, J. Dietary requirements of individual nectar foragers, and colony-level pollen and nectar consumption: A review to support pesticide exposure assessment for honey bees. Apidologie 51, 163–179. https://doi.org/10.1007/s13592-019-00694-9 (2020).Article 

Google Scholar 
Drezner-Levy, T., Smith, B. & Shafir, S. The effect of foraging specialization on various learning tasks in the honey bee (Apis mellifera). Behav. Ecol. Sociobiol. 64, 135–148. https://doi.org/10.1007/s00265-009-0829-z (2009).Article 

Google Scholar 
Afik, O. & Shafir, S. Effect of ambient temperature on crop loading in the honey bee, Apis mellifera (Hymenoptera: Apidae). Entomologia Generalis 29, 135–148 (2007).Article 

Google Scholar 
Seeley, T. D. Honey bee foragers as sensory units of their colonies. Behav. Ecol. Sociobiol. 34, 51–62 (1994).Article 

Google Scholar 
Waller, G. D. Evaluating responses of honeybees to sugar solutions using an artificial-flower feeder. Ann. Entomol. Soc. Am. 65, 857–862 (1972).CAS 
Article 

Google Scholar 
Nicolson, S. W., de Veer, L., Köhler, A. & Pirk, C. W. W. Honeybees prefer warmer nectar and less viscous nectar, regardless of sugar concentration. Proc. R. Soc. B: Biol. Sci. 280, 20131597. https://doi.org/10.1098/rspb.2013.1597 (2013).Article 

Google Scholar 
Neff, J. L. & Simpson, B. B. The roles of phenology and reward structure in the pollination biology of wild sunflower (Helianthus annuus L., Asteraceae). Israel J. Bot. 39, 197–216 (1990).
Google Scholar 
Waller, G. D., Carpenter, E. W. & Ziehl, O. A. Potassium in onion nectar and its probable effect on attractiveness of onion flowers to honey bees. J. Am. Soc. Hortic. Sci. 97, 535–539 (1972).CAS 
Article 

Google Scholar 
Roubik, D. W., Yanega, D., Aluja, M., Buchmann, S. L. & Inouye, D. W. On optimal nectar foraging by some tropical bees (Hymenoptera: Apidae). Apidologie 26, 197–211 (1995).Article 

Google Scholar 
Power, E. F., Stabler, D., Borland, A. M., Barnes, J. & Wright, G. A. Analysis of nectar from low-volume flowers: A comparison of collection methods for free amino acids. Methods Ecol. Evol. 9, 734–743. https://doi.org/10.1111/2041-210X.12928 (2018).Article 
PubMed 

Google Scholar 
Pattrick, J. G., Symington, H. A., Federle, W. & Glover, B. J. The mechanics of nectar offloading in the bumblebee Bombus terrestris and implications for optimal concentrations during nectar foraging. J. R. Soc. Interface 17, 20190632. https://doi.org/10.1098/rsif.2019.0632 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Strauss, U., Dietemann, V., Human, H., Crewe, R. M. & Pirk, C. W. W. Resistance rather than tolerance explains survival of savannah honeybees (Apis mellifera scutellata) to infestation by the parasitic mite Varroa destructor. Parasitology 143, 374–387. https://doi.org/10.1017/s0031182015001754 (2016).Article 
PubMed 

Google Scholar 
Dyer, F. C. & Seeley, T. D. Interspecific comparisons of endothermy in honey-bees (Apis): Deviations from the expected size-related patterns. J. Exp. Biol. 127, 1–26. https://doi.org/10.1242/jeb.127.1.1 (1987).Article 

Google Scholar  More

Barrientos, R. et al. A review of searcher efficiency and carcass persistence in infrastructure-driven mortality assessment studies. Biol. Conserv. 222, 146–153 (2018).
Google Scholar 
Stevens, B. S., Reese, K. P. & Connelly, J. W. Survival and detectability bias of avian fence collision surveys in sagebrush steppe. J. Wildl. Manag. 75, 437–449 (2011).
Google Scholar 
Hunting, K. A Roadmap for PIER Research on Avian Collisions with Power Lines in California. (2002).Barrientos, R. et al. Wire marking results in a small but significant reduction in avian mortality at power lines: A baci designed study. PLoS ONE 7, e32569 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Costantini, D., Gustin, M., Ferrarini, A. & Dell’Omo, G. Estimates of avian collision with power lines and carcass disappearance across differing environments. Anim. Conserv. 20, 173–181 (2017).
Google Scholar 
Jenkins, A. R. et al. Estimating the impacts of power line collisions on Ludwig’s Bustards Neotis ludwigii. Bird Conserv. Int. 21, 303–310 (2011).
Google Scholar 
Shaw, J. M., Reid, T. A., Schutgens, M., Jenkins, A. R. & Ryan, P. G. High power line collision mortality of threatened bustards at a regional scale in the Karoo, South Africa. Ibis (Lond. 1859) 1859(160), 431–446 (2018).
Google Scholar 
Gómez-Catasús, J. et al. Factors affecting differential underestimates of bird collision fatalities at electric lines: a case study in the Canary Islands. Ardeola 68, 71–94 (2021).
Google Scholar 
Ponce, C., Alonso, J. C., Argandoña, G., García Fernández, A. & Carrasco, M. Carcass removal by scavengers and search accuracy affect bird mortality estimates at power lines. Anim. Conserv. 13, 603–612 (2010).
Google Scholar 
Bernardino, J. et al. Bird collisions with power lines: State of the art and priority areas for research. Biol. Conserv. 222, 1–13 (2018).
Google Scholar 
Brooks, J. W. & Sutton, L. in Veterinary Forensic Pathology (ed. Brooks, J. W.) 43–63 (2018). https://doi.org/10.1007/978-3-319-67172-7_4Brooks, J. W. Postmortem changes in animal carcasses and estimation of the postmortem interval. Vet. Pathol. 53, 929–940 (2016).CAS 
PubMed 

Google Scholar 
Ascensão, F. et al. Beware that the lack of wildlife mortality records can mask a serious impact of linear infrastructures. Glob. Ecol. Conserv. 19, e00661 (2019).
Google Scholar 
Hau, T. C., Hamzah, N. H., Lian, H. H. & Amir Hamzah, S. P. A. Decomposition process and post mortem changes: Review. Sains Malaysiana 43, 1873–1882 (2014).
Google Scholar 
Cooper, J. E. in Wildlife Forensic Investigation: Principles and Practice (eds. Cooper, J. & Cooper, M.) 237–324 (CRC Press, 2013). https://doi.org/10.1201/b14553Sutherland, A., Myburgh, J., Steyn, M. & Becker, P. J. The effect of body size on the rate of decomposition in a temperate region of South Africa. Forensic Sci. Int. 231, 257–262 (2013).CAS 
PubMed 

Google Scholar 
Valverde, I., Espín, S., María-Mojica, P. & García-Fernández, A. J. Protocol to classify the stages of carcass decomposition and estimate the time of death in small-size raptors. Eur. J. Wildl. Res. 66, 1–13 (2020).
Google Scholar 
Goff, M. L. in Current Concepts in Forensic Entomology (eds. Amendt, J., Goff, M., Campobasso, C. & Grassberger, M.) 1–24 (Springer, 2010). https://doi.org/10.1007/978-1-4020-9684-6_1Pittner, S. et al. A field study to evaluate PMI estimation methods for advanced decomposition stages. Int. J. Legal Med. 134, 1361–1373 (2020).PubMed 
PubMed Central 

Google Scholar 
Probst, C. et al. Estimating the postmortem interval of wild boar carcasses. Vet. Sci. 7, 6 (2020).PubMed Central 

Google Scholar 
Cambra-Moo, Ó., Delgado-Buscalioni, Á. & Delgado-Buscalioni, R. An approach to the study of variations in early stages of Gallus gallus decomposition. J. Taphon. 6, 21–40 (2008).
Google Scholar 
Oates, D., Coggin, J., Hartman, F. & Hoilien, G. Guide to Time of Death in Selected Wildlife Species. (Nebraska Technical Series No. 14. Lincoln, N.E., Nebraska Game and Parks Commission, 1984).Hewadikaram, K. A. & Goff, M. L. Effect of carcass size on rate of decomposition and arthropod succession patterns. Am. J. Forensic Med. Pathol. 12, 240–265 (1991).
Google Scholar 
Zhou, C. & Byard, R. W. Factors and processes causing accelerated decomposition in human cadavers—An overview. J. Forensic Leg. Med. 18, 6–9 (2011).PubMed 

Google Scholar 
Cockle, D. L. & Bell, L. S. Human decomposition and the reliability of a ‘Universal’ model for post mortem interval estimations. Forensic Sci. Int. 253(136), e1-136.e9 (2015).
Google Scholar 
Azevedo, R. R. & Krüger, R. F. The influence of temperature and humidity on abundance and richness of Calliphoridae (Diptera). Iheringia. Série Zool. 103, 145–152 (2013).
Google Scholar 
Barnes, K. M. in Wildlife Forensic Investigation: Principles and Practice (eds. Cooper, J. & Cooper, M.) 149–160 (CRC Press, 2013).Mann, R. W., Bass, W. M. & Meadows, L. Time since death and decomposition of the human body: Variables and observations in case and experimental field studies. J. Forensic Sci. 35, 103–111 (1990).CAS 
PubMed 

Google Scholar 
Gliksman, D. et al. Biotic degradation at night, abiotic degradation at day: Positive feedbacks on litter decomposition in drylands. Glob. Change Biol. 23, 1564–1574 (2017).ADS 

Google Scholar 
Araujo, P. I., Grasso, A. A., González-Arzac, A., Méndez, M. S. & Austin, A. T. Sunlight and soil biota accelerate decomposition of crop residues in the Argentine Pampas. Agric. Ecosyst. Environ. 330, 107908 (2022).
Google Scholar 
Fernández-Palacios, J. M. & Martín-Esquivel, J. L. Naturaleza de las Islas Canarias: Ecología y Conservación. (Turquesa, 2001).Kenward, M. G. & Roger, J. H. An improved approximation to the precision of fixed effects from restricted maximum likelihood. Comput. Stat. Data Anal. 53, 2583–2595 (2009).MathSciNet 
MATH 

Google Scholar 
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. http://www.r-project.org (2020).Bates, D., Mächler, M., Bolker, B. M. & Walker, S. C. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).
Google Scholar 
Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. B. lmerTest package: Tests in linear mixed effects models. J. Stat. Softw. 82, 1–26 (2017).
Google Scholar 
Zeileis, A. & Hothorn, T. Diagnostic checking in regression relationships. R News 2, 7–10 (2002).
Google Scholar 
Halekoh, U. & Højsgaard, S. A Kenward–Roger approximation and parametric bootstrap methods for tests in linear mixed models-the R package pbkrtest. J. Stat. Softw. 59, 1–30 (2014).
Google Scholar 
Fox, J. & Weisberg, S. An {R} Companion to Applied Regression, Second Edition. (Sage, 2011).Bartoń, K. MuMIn: Multi-Model Inference. (R Package Version 1.43.6, 2019).De Rosario-Martinez, H., Fox, J. & R Core Team. Package ‘phia’ Title Post-Hoc Interaction Analysis. (R Package Version 0.2–1, 2015).Nakagawa, S. & Schielzeth, H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol. Evol. 4, 133–142 (2013).
Google Scholar 
Vass, A. Beyond the grave—Understanding human decomposition. Microbiol. Today 28, 190–192 (2001).
Google Scholar 
Gill-King, H. in Forensic Taphonomy: The Postmortem Fate of Human Remains (eds. Haglund, W. D. & Sorg, M. H.) 93–104 (CRC Press, 1996). https://doi.org/10.1201/9781439821923.sec2Campobasso, C. P., Di Vella, G. & Introna, F. Factors affecting decomposition and Diptera colonization. Forensic Sci. Int. 12, 18–27 (2001).
Google Scholar 
Austin, A. T., Araujo, P. I. & Leva, P. E. Interaction of position, litter type, and water pulses on decomposition of grasses from the semiarid Patagonian steppe. Ecology 90, 2642–2647 (2009).PubMed 

Google Scholar 
Brandt, L. A., Bonnet, C. & King, J. Y. Photochemically induced carbon dioxide production as a mechanism for carbon loss from plant litter in arid ecosystems. J. Geophys. Res. Biogeosci. 114, G02004 (2009).ADS 

Google Scholar 
Lee, H., Rahn, T. & Throop, H. An accounting of C-based trace gas release during abiotic plant litter degradation. Glob. Chang. Biol. 18, 1185–1195 (2012).ADS 

Google Scholar 
Zepp, R. G., Erickson, D. J., Paul, N. D. & Sulzberger, B. Interactive effects of solar UV radiation and climate change on biogeochemical cycling. Photochem. Photobiol. Sci. 6, 286–300 (2007).CAS 
PubMed 

Google Scholar 
Archer, M. S. Rainfall and temperature effects on the decomposition rate of exposed neonatal remains. Sci. Justice J. Forensic Sci. Soc. 44, 35–41 (2004).Simmons, T., Adlam, R. E. & Moffatt, C. Debugging decomposition data—Comparative taphonomic studies and the influence of insects and carcass size on decomposition rate. J. Forensic Sci. 55, 8–13 (2010).PubMed 

Google Scholar 
Spicka, A., Johnson, R., Bushing, J., Higley, L. G. & Carter, D. O. Carcass mass can influence rate of decomposition and release of ninhydrin-reactive nitrogen into gravesoil. Forensic Sci. Int. 209, 80–85 (2011).CAS 
PubMed 

Google Scholar 
Tracqui. in Encyclopaedia of Forensic Sciences (eds. Siegel, J. A., Saukko, P. J. & Max, M. H.) 1357–1363 (Academic Press, 2000).Riding, C. S. & Loss, S. R. Factors influencing experimental estimation of scavenger removal and observer detection in bird–window collision surveys. Ecol. Appl. 28, 2119–2129 (2018).PubMed 

Google Scholar  More

Dataset 1: Analysis of truffle growers archivesWe selected eleven T. melanosporum orchards located across the South-West France, from Montpellier (43°44′01.4″N 3°42′13.2″E) to Jonzac (45°27′17.7″N, 0°25′26.9″W; Fig. 2). These sites were selected for (1) the quality of the records of fruitbody production and practices by truffle growers (Table S1), including the detail of inoculations since plantation (amount and frequency of added crushed sporocarps), (2) the use of truffle traps by the owners and the quality of the record from these devices, and (3) the presence of oaks (Quercus ilex, Q. pubescens and Q. suber) as the only hosts tree species. Based on the archives of truffle growers, including a systematic recording of truffle production within and outside traps, we reported at each study site the contribution of truffle traps to the annual fruitbody production of the entire truffle grounds, by using number and/or weight of collected fruitbodies within (Pin) and outside (Pout) truffle traps.Dataset 2: In situ experiment tracing the inoculation effectThree orchards located near Angoulème (45°74′35.5″N, − 0°63′78.4″W), Jonzac (45°44′09.8″N, 0°43′96.7″W), and Arles-sur-Tech (42°45′44.9″N, 2°62′89.4″W), hereafter referred to Site 1 to 3 (Fig. 2) were selected for testing both disturbance effect and inoculum effect on fruitbody production in truffle traps. These sites presented a high fruitbody production and a high Pin/Pout ratio, thus optimum conditions to test mechanisms underlying how truffle traps influence fruitbody production. Host trees were between 5 and 18 years old at the beginning of the experiment (Fig. 2). At each site, we selected three non-adjacent trees (four on Site 3) that displayed a continuous fruitbody production over the three previous years. Under each selected tree, we excavated, at two-thirds of the distance between the tree trunk and the limit of brûlé (a vegetation-poor zone that shows the extension mycelia in the soil40, eight equidistant truffle traps [20 × 20 cm large × 20 cm deep] as shown in Fig. 3a. Under each tree, two traps were filled with only a mixture of peat and vermiculite (hereafter referred as non-inoculated controls) to test for disturbance effect. The used mixture was identical to that which is currently applied in commercial orchards. In three other traps, 5 g of crushed material from a single black truffle fruitbody (including its gleba and spores) were added to the previous mixture (hereafter referred as one mating-type inoculum). In the three last traps, 5 g of crushed material from two ascocarps with gleba of opposite mating types (hereafter referred as two mating-type inoculum) were added to the previous mixture. We added the two mating-type condition to accurately test a potential contribution of the gleba (haploid and thus with a single mating type) on future production. As quoted in Introduction, maternal individuals with opposite mating types tend to exclude each other locally (spatial segregation of clusters of individuals of same mating types26. Thus, the two mating-type inoculum allows us to detect in each trap a maternal contribution by the introduced gleba, despite potential exclusion by pre-installed individuals of the locally dominant mating type in the surrounding. Moreover, it allows us to detect a paternal contribution by the introduced gleba of the mating type opposite to the locally dominant. The eight truffle traps were randomly arranged, so that two repetitions of same modality were always separated by a repetition of another modality (Fig. 3a).In March 2013, six freshly collected truffles (weighting  > 60 g) were molecularly analyzed for the mating type of their gleba as in18. On Site 1 and Site 2, the inoculum was made of fruitbodies collected at Site 1. On Site 3, fruitbodies used as inoculum originated from truffle grounds in Sarrion (Spain). In April 2013, truffles traps were installed as explained above (in all, 8 traps × 3 (or 4) trees × 3 sites) and monitored for two years by truffle growers. Harvesting was performed by trained dogs (one different dog per site) checking truffle traps and the surrounding brûlés at each visit of the orchard by truffle growers. When dogs detected truffles, a small hole was excavated to collect ascocarps without disturbing the trap further. At the end of January, 2015, all truffle traps were completely excavated, remnant truffles overlooked by dogs were systematically collected (Fig. 3b). Three soil aliquots were collected within all traps and pooled. All truffles and soil aliquots were frozen for subsequent DNA analysis.Molecular and genetic analysesDNA extractions, mating typing and genotyping were done as in18. Briefly, DNA was extracted from the gleba and from spores of each fruitbody to get access to the maternal and zygotic DNA, respectively. Simple sequence repeat (SSRs) genotyping was performed using 12 polymorphic markers and the mating-type locus as in18. Gleba extracts displaying apparent heterozygous genotypes, likely due to contamination by spore DNA were systematically discarded from further analyses. For each fruitbody, the haploid paternal genotype was then deduced by subtracting the haploid maternal genotype from the zygotic diploid genotype. This data set was used for relatedness estimations. We discarded from all further analysis the marker me11, which displayed more than 39% missing data, as well as all samples with missing data for at any locus.Multilocus genotypes comparisonsBased on the 11 remaining SSRs and the mating-type (Table S5 and Figure S2), MLGs were identified on all maternal and paternal haploid genomes using GenClone v.2.041, and the probability that MLGs represented more than once resulted from independent events of sexual reproduction was calculated (PSex41,42). On each site, clonal diversity was measured as R = (G − 1)/(N − 1) according to43, where N is the number of fruitbodies and G the number of MLGs. For testing whether the gleba of the inoculated fruitbody contributed, either paternally (H1) or maternally (H2) to the harvested fruitbodies (Fig. 1c), the inoculated maternal MLG was compared to the paternal and maternal MLG of the harvested fruitbodies.Relatedness estimationFor testing whether the spores of the inoculum, which carry many distinct haploid MLGs due to meiosis, had paternal or maternal contribution(s) to the harvested fruitbodies (H3; Fig. 1c), we used relatedness estimation.For testing whether spores of the inoculum had a paternal contribution, an individual relatedness estimate to the spore inoculum was computed for each paternal genome detected in truffle traps. Relatedness r here describes the expected frequency E[p_offpat] of each allele in a given genome, E[p_offpat] = p_pop + r * (p_inoc − p_pop), where p_pop is the allele frequency in the local population (here estimated from the glebas of other truffles collected under the focal tree), and p_inoc is the frequency of the allele in the inoculum. Thus, p_offpat takes values 0 or 1, and p_inoc takes values 0, 0.5 or 1, except when two fruitbodies were used as inoculum (two gleba mating types traps). Thus r = (p_offpat − p_pop)/(p_inoc − p_pop). An individual relatedness estimate for each genome is then obtained by summing over alleles and loci the observed values of the numerator and denominator in this expression. A population-level estimate is further obtained by summing numerators and denominators over the paternity events in each population.To test whether such estimates are compatible with the hypothesis that the paternal individuals are not from the inocula, we obtained the distribution of population-level relatedness estimates by simulating samples under this hypothesis: paternal genotypes were randomly simulated according to alleles frequencies in the local population. For each population, 10,000 samples were simulated, and p-values were estimated as the proportion of simulations with higher population-level relatedness with inocula than the observed one. Confidence intervals for these p-values were computed from the binomial distribution for 10,000 draws, and Bonferroni-corrected over the three populations.For testing whether spores of the inoculum had a maternal contribution (H4, Fig. 1c), we estimated the relatedness of the locally used spore inoculum to each maternal genome detected in truffle traps (deduced from the gleba), and we confronted it to simulated samples as previously but with one modification: if the focal fruitbody was harvested in a trap inoculated with the inoculum A1, all genomes of truffles from traps inoculated with the same inoculum (A1 or A1 + A2 + A3, see Fig. 3c.) were discarded from the estimation of p_pop.Assessment of T. melanosporum mycelium concentration in truffle trapsOn Sites 1, 2 and 3, soil samples were collected in all traps and in the surrounding brûlés at harvesting date (January, 2015). In collected soils, total DNA was extracted and quantified as in19. Briefly, after sieving and homogenizing soil collected in each trap and from out of the brûlés, aliquots (10 g) were analyzed as follows. After extraction with the kit Power Soil (MoBio Laboratories, Carlsbad, CA, USA), the extra-radical mycelium of T. melanosporum was quantified using quantitative Taqman™ PCR (qPCR) with the primers and probe described in44. Triplicate real-time PCR were performed on each sample using the same concentration of primer and the same thermocycling program as in19. Standards were prepared using fresh immature T. melanosporum ascocarp, and a standard curve was generated for each site by plotting serial tenfold dilutions against corresponding initial amount of ascocarp. Absolute quantification of mycelium biomass of T. melanosporum was expressed in mg of mycelium per g of soil.Statistical analysesStatistics were done using R version 4.0.445.Effect of truffle traps on fruitbody production—The contribution of truffle traps to the overall production of orchards was assessed by (1) data mining of truffle growers’ archives (Dataset 1) and (2) comparing the density of truffles harvested in traps (expressed in number of truffles per m2 per orchard; for each sampled tree, traps correspond to an investigated soil surface of s = 8 × 0.2 x 0.2 = 0.32 m2) with the density measured within surrounding brûlés (Dataset 1). On Dataset2, at each site, the area occupied by brûlés was evaluated by measuring in the field the surface of soil devoid of vegetation consecutively to spontaneous T. melanosporum brûlé.Fruitbody production under different conditions (i.e. non-inoculated controls versus one gleba mating type traps versus two gleba mating type traps) were compared using generalized linear mixed models with negative binominal family and log link (R, spam package46). The full model included the logarithm of the sampled area as offset to account for variations in this sampled area, interactions of trap-modality effects with site effect. Formal likelihood ratio tests are based on one-step deletions from this full model, applied to subsets of the data relevant for each hypothesis tested. Additional bootstrap tests (1000 iterations) were run to correct any bias in small sample likelihood ratio tests.Concentrations of T. melanosporum mycelium in soil—Similarly as above, the inoculum effect on mycelium concentrations was compared using generalized linear mixed models with Gamma log family.Plant materialThe use of plants in the present study complies with international, national and/or institutional guidelines. All permissions to collect T. melanosporum fruitbodies in truffle orchards were obtained. The formal identification of biological material used in the study (T. melanosporum fruitbodies) was undertaken by F. Richard and E. Taschen. Voucher specimens of all collected fruitbodies have been deposited in the Centre d’Ecologie Fonctionnelle et Evolutive herbarium in Montpellier (France).Ethical approvalAll co-authors approve the ethical statement regarding the submitted manuscript.Consent to participateAll co-authors consent to participate to the research and agree with the content of the submitted manuscript. All authors reviewed and submitted manuscript. More

Changes in temperature and extreme environmental conditions can affect the dynamics of vector-borne pathogens. These include leishmaniasis, transmitted by phlebotomine sandflies, as well as mosquitoes that spread arboviruses like dengue, encephalitis, yellow fever, West Nile fever, and lymphatic filariasis19,20,21.The CCA analysis showed that maximum temperature significantly influenced the abundance of mosquito populations in the study area. In addition, the NMDS showed two different groupings that consisted of samples collected during the rainy and dry seasons. Accordingly, Refs.22,23 report that changes in temperature and relative humidity determine the abundance of mosquitoes, which can disappear entirely during the dry season. Moreover, Refs.22,24,25 note that certain species of mosquitoes increase proportionally with the regional rainfall regime. This is consistent with Ref.10, who find alternating patterns in tropical and temperate climates in some Brazilian regions.As shown by the geometric regression, there is a positive correlation between cumulative rainfall in the days before collection and the number of species found in the study period. Likewise, Ref.26 reported that under the conditions observed in the Serra do Mar State Park, climate variables directly influenced the abundance of Cq. chrysonotum and Cq. venezuelensis, favoring the occurrence of culicids during the more warm, wet, and rainy months.The current climate scenario and future projections about climate, environmental, demographic, and meteorological factors directly influence the distribution and abundance of mosquito vectors and/or diseases27,28,29,30. Environmental temperature alters mosquito population dynamics, thereby affecting the development of immature stages as well as reproduction31. While temperature has an important effect on population dynamics, rainfall and drought also affect the density and dispersal of mosquitoes in temperate and tropical regions32.To be sure, environmental changes other than climate can modify the behavior of vector insects and, subsequently, the mechanism of transmission of parasites20. Specifically, human impacts on the environment can result in drastically different disease transmission cycles in and around inhabited areas33.A previous study34 reported that changes in land use influence the mosquito communities with potential implications for the emergence of arboviruses. Another study35 noted that environmental changes negatively affect natural ecosystems with accelerated biodiversity loss. This is due to the modification and loss of natural habitat and unsustainable land use, which leads to the spread of pathogens and disease vectors.Hence, understanding the relationship between humans and the environment becomes increasingly critical, given the way in which climate changes can lead to alterations in the epidemiology of diseases such as dengue in areas considered free of the disease, as well as in endemic areas36.We found that the abundance and diversity of Mansoniini were directly influenced by the effect of the rainy season and other climatic factors. The rainfall regime has been shown to affect the development of immature forms12,37; explaining the greater frequency of these specimens in the warmer and wetter months38,39,40. According to Ref.41, stable ecosystems such as forests contain great species diversity. On the other hand, diversity tends to be reduced in biotic communities suffering from stress.Studies of insect populations in natural areas are important because they allow a direct analysis of how environmental factors influence phenomena such as the choice of breeding sites by females for oviposition, hematophagous behavior, and the distribution of species along a vegetation gradient12,26,42,43.Throughout the experimental period of the present study, we observed that Shannon light traps are an effective method for catching mosquitoes from the Mansoniini tribe. Interestingly, Ref.44 reported a species richness pattern strongly influenced by Coquillettidia fasciolata (Lynch Arribálzaga, 1891) on mosquito samples from different capture points by using CDC and Shannon light traps as sampling methods. In contrast to the results of Ref.44, where the highest population density of mosquitoes was captured with CDC traps, we observed that these traps were not effective at capturing specimens of Mansoniini in spite of being used in large numbers in the present study. Moreover, Ref.45 conducted another study on faunal diversity in an Atlantic Forest remnant of the state of Rio de Janeiro and observed the highest abundance of Cq. chrysonotum (Peryassú, 1922) and Cq. venezuelensis by using Shannon light traps, while the numbers of captures of Ma. titillans were very similar using CDC and Shannon traps.The results of this study indicate that the makeup of culicid fauna remains quite similar throughout the year, despite seasonal variations in abundance, though there was a lower variability of fauna in the dry season. Therefore, although the seasonality did not affect the temporal variation of the faunal composition in a generalized way, it was possible to detect a partial effect of the seasonality on fauna abundance.
Reference46 report that the incidence peaks of mosquitoes in the warmer and wetter months, as well as mosquito populations remaining between tolerance limits for most of the year, indicate the sensitivity of some species to the local climate.The elevated abundance and diversity of species of Mansoniini in the study area were influenced by the favorable maintenance of breeding sites, including specific water accumulations with emerging vegetation that remain present throughout the year and the well-defined rainy season in the region. In addition, the representatives of Mansoniini, which prefer breeding sites containing macrophytes, made up nearly all of the species collected7.Besides providing a greater awareness of mosquito populations’ ecological and biological aspects, research carried out in wild areas also provides information on the relationship between species diversity and the area in which they are found. Considering that wild insects may become potential vectors of diseases, research in wild areas also provides helpful information for understanding relevant epidemiological aspects. These studies facilitate the identification, monitoring, and control of mosquito populations following environmental changes caused by direct human action, which can lead to major epidemics26.We observed considerable heterogeneity among Mansoniini fauna, and the months with the highest rainfall directly influence the structure of the communities and contribute to the increase in mosquito diversity and abundance, possibly due to variations in the availability of habitat for their immature forms. More