Philippa Kaur

1.Vittadini C. Monographia lycoperdineorum. Augustae Taurinorum, Torino, 1842.2.Frank B. On the nutrition of certain trees by underground fungi based on root symbiosis. Plant Biol. 1885;3:128–45.
Google Scholar 
3.Gadgil RL, Gadgil P. Mycorrhiza and litter decomposition. Nature 1971;233:133–133.CAS 
PubMed 

Google Scholar 
4.Harley J. Problems of mycotrophy. London: Academic Press; 1975.5.Clemmensen KE, Finlay RD, Dahlberg A, Stenlid J, Wardle DA, Lindahl BD. Carbon sequestration is related to mycorrhizal fungal community shifts during long‐term succession in boreal forests. N. Phytol. 2015;205:1525–36.CAS 

Google Scholar 
6.Crowther TW, Van den Hoogen J, Wan J, Mayes MA, Keiser A, Mo L, et al. The global soil community and its influence on biogeochemistry. Science 2019;365:eaav0550.CAS 
PubMed 

Google Scholar 
7.Bueno CG, Moora M, Gerz M, Davison J, Öpik M, Pärtel M, et al. Plant mycorrhizal status, but not type, shifts with latitude and elevation in Europe. Glob Ecol Biogeo. 2017;26:690–9.
Google Scholar 
8.Steidinger BS, Crowther TW, Liang J, Nuland MEV, Werner GDA, Reich PB, et al. Climatic controls of decomposition drive the global biogeography of forest-tree symbioses. Nature 2019;569:404–8.CAS 
PubMed 

Google Scholar 
9.Dixon RK, Garrett HE, Cox GS, Marx DH, Sander IL. Inoculation of three Quercus species with eleven isolates of ectomycorrhizal fungi. I. inoculation success and seedling growth relationships. Science. 1984;30:364–72.
Google Scholar 
10.Sim M-Y, Eom A-H. Effects of ectomycorrhizal fungi on growth of seedlings of Pinus densiflora. Mycobiology 2006;34:191–5.PubMed 
PubMed Central 

Google Scholar 
11.Dickie IA. Host preference, niches and fungal diversity. N. Phytol. 2007;174:230–3.
Google Scholar 
12.Alberton O, Kuyper TW, Gorissen A. Competition for nitrogen between Pinus sylvestris and ectomycorrhizal fungi generates potential for negative feedback under elevated CO2. Plant Soil. 2007;296:159–72.CAS 

Google Scholar 
13.Karst J, Marczak L, Jones MD, Turkington R. The mutualism–parasitism continuum in ectomycorrhizas: a quantitative assessment using meta-analysis. Ecology 2008;89:1032–42.PubMed 

Google Scholar 
14.Dalong M, Luhe W, Guoting Y, Liqiang M, Chun L. Growth response of Pinus densiflora seedlings inoculated with three indigenous ectomycorrhizal fungi in combination. Braz J Microbiol. 2011;42:1197–203.CAS 
PubMed 
PubMed Central 

Google Scholar 
15.Bennett JA, Maherali H, Reinhart KO, Lekberg Y, Hart MM, Klironomos J. Plant-soil feedbacks and mycorrhizal type influence temperate forest population dynamics. Science 2017;355:181–4.CAS 
PubMed 

Google Scholar 
16.Read DJ, Perez‐Moreno J. Mycorrhizas and nutrient cycling in ecosystems – a journey towards relevance? N. Phytol. 2003;157:475–92.CAS 

Google Scholar 
17.Buscot F, Weber G, Oberwinkler F. Interactions between Cylindrocarpon destructans and ectomycorrhizas of Picea abies with Laccaria laccata and Paxillus involutes. Trees. 1992;6:83–90.
Google Scholar 
18.Morin C, Samson J, Dessureault M. Protection of black spruce seedlings against Cylindrocladium root rot with ectomycorrhizal fungi. Can J Bot. 1999;77:169–74.
Google Scholar 
19.Abuzinadah RA, Read DJ. The role of proteins in the nitrogen nutrition of ectomycorrhizal plants. N. Phytol. 1989;112:55–60.CAS 

Google Scholar 
20.Jongbloed RH, Clement JMAM, Borst-Pauwels GWFH. Kinetics of NH4+ and K+ uptake by ectomycorrhizal fungi. effect of NH4+ on K+ uptake. Phys Plant 1991;83:427–32.CAS 

Google Scholar 
21.Selosse M, Bouchard D, Martin F, Tacon F. Effect of Laccaria bicolor strains inoculated on Douglas-fir (Pseudotsuga menziesii) several years after nursery inoculation. Can J Res. 2000;30:360–71.
Google Scholar 
22.Hoeksema JD, Chaudhary VB, Gehring CA, Johnson NC, Karst J, Koide RT, et al. A meta-analysis of context-dependency in plant response to inoculation with mycorrhizal fungi. Ecol Let. 2010;13:394–407.
Google Scholar 
23.Kipfer T, Wohlgemuth T, Heijden MGA, van der, Ghazoul J, Egli S. Growth response of drought-stressed Pinus sylvestris seedlings to single- and multi-species inoculation with ectomycorrhizal Fungi. PLoS ONE. 2012;7:e35275.CAS 
PubMed 
PubMed Central 

Google Scholar 
24.Pena R, Polle A. Attributing functions to ectomycorrhizal fungal identities in assemblages for nitrogen acquisition under stress. ISME J. 2014;8:321–30.CAS 
PubMed 

Google Scholar 
25.Mueller RC, Scudder CM, Whitham TG, Gehring CA. Legacy effects of tree mortality mediated by ectomycorrhizal fungal communities. N. Phytol. 2019;224:155–65.CAS 

Google Scholar 
26.Policelli N, Horton TR, Hudon AT, Patterson TR, Bhatnagar JM. Back to roots: the role of ectomycorrhizal fungi in boreal and temperate forest restoration. Front Glob Change. 2020;3:97.
Google Scholar 
27.Bever JD, Schultz PA, Pringle A, Morton JB. Arbuscular mycorrhizal fungi: more diverse than meets the eye, and the ecological tale of Why: the high diversity of ecologically distinct species of arbuscular mycorrhizal fungi within a single community has broad implications for plant ecology. BioScience 2001;51:923–31.
Google Scholar 
28.Delgado‐Baquerizo M, Giaramida L, Reich PB, Khachane AN, Hamonts K, Edwards C, et al. Lack of functional redundancy in the relationship between microbial diversity and ecosystem functioning. J Ecol. 2016;104:936–46.
Google Scholar 
29.Nelson MB, Martiny AC, Martiny JBH. Global biogeography of microbial nitrogen-cycling traits in soil. PNAS 2016;113:8033–40.CAS 
PubMed 
PubMed Central 

Google Scholar 
30.Louca S, Parfrey LW, Doebeli M. Decoupling function and taxonomy in the global ocean microbiome. Science 2016;353:1272–7.CAS 
PubMed 

Google Scholar 
31.Louca S, Jacques SMS, Pires APF, Leal JS, Srivastava DS, Parfrey LW, et al. High taxonomic variability despite stable functional structure across microbial communities. Nat Ecol Evol. 2016;1:1–12.
Google Scholar 
32.Louca S, Polz MF, Mazel F, Albright MBN, Huber JA, O’Connor MI, et al. Function and functional redundancy in microbial systems. Nat Ecol Evol. 2018;2:936–43.
Google Scholar 
33.Lindahl BD, Kyaschenko J, Varenius K, Clemmensen KE, Dahlberg A, Karltun E, et al. A group of ectomycorrhizal fungi restricts organic matter accumulation in boreal forest. Ecol Lett. 2021;24:1341–51.PubMed 

Google Scholar 
34.Rineau F, Courty P-E. Secreted enzymatic activities of ectomycorrhizal fungi as a case study of functional diversity and functional redundancy. Ann Sci. 2011;68:69–80.
Google Scholar 
35.Talbot JM, Bruns TD, Taylor JW, Smith DP, Branco S, Glassman SI, et al. Endemism and functional convergence across the North American soil mycobiome. PNAS 2014;111:6341–6.CAS 
PubMed 
PubMed Central 

Google Scholar 
36.Banerjee S, Kirkby CA, Schmutter D, Bissett A, Kirkegaard JA, Richardson AE. Network analysis reveals functional redundancy and keystone taxa amongst bacterial and fungal communities during organic matter decomposition in an arable soil. Soil Bio Biochem. 2016;97:188–98.CAS 

Google Scholar 
37.Hall EK, Bernhardt ES, Bier RL, Bradford MA, Boot CM, Cotner JB, et al. Understanding how microbiomes influence the systems they inhabit. Nat Microbiol. 2018;3:977–82.CAS 
PubMed 

Google Scholar 
38.Etzold S, Ferretti M, Reinds GJ, Solberg S, Gessler A, Waldner P, et al. Nitrogen deposition is the most important environmental driver of growth of pure, even-aged and managed European forests. Ecol Man. 2020;458:117762.
Google Scholar 
39.Van der Linde S, Suz LM, Orme CDL, Cox F, Andreae H, Asi E, et al. Environment and host as large-scale controls of ectomycorrhizal fungi. Nature 2018;558:243–8.PubMed 

Google Scholar 
40.Ferretti M, Fischer R Forest Monitoring: Methods for Terrestrial Investigations in Europe with an Overview of North America and Asia in Developments in Environmental Science. vol. 12. Elsevier, Amsterdam, 2013. pp 2-507.41.Dobbertin M, Neumann M Part V: Tree Growth. In: UNECE ICP Forests, Programme Co- ordinating Centre (ed.): Manual on methods and criteria for harmonized sampling, assessment, monitoring and analysis of the effects of air pollution on forests. Thünen Institute of Forest Ecosystems. Eberswalde. 2016. https://www.icp-forests.org/pdf/manual/2016/ICP_Manual_2016_01_part05.pdf.42.Averill C, Cates LL, Dietze MC, Bhatnagar JM. Spatial vs. temporal controls over soil fungal community similarity at continental and global scales. ISME J. 2019;13:2082–93.PubMed 
PubMed Central 

Google Scholar 
43.Pellitier PT, Ibáñez I, Zak DR, Argiroff WA, Acharya K. Ectomycorrhizal access to organic nitrogen mediates CO2 fertilization response in a dominant temperate tree. Nat Commun. 2021;12:5403.CAS 
PubMed 
PubMed Central 

Google Scholar 
44.Henry M, Bombelli A, Trotta C, Alessandrini A, Birigazzi L, Sola G, et al. GlobAllomeTree: international platform for tree allometric equations to support volume, biomass and carbon assessment. iFor – Biogeo. 2013;6:326–30.
Google Scholar 
45.Penman J, Gytarsky M, Hiraishi T, Krug T, Kruger D, Pipatti R, et al. Good practice guidance for land use, land-use change and forestry. Good practice guidance for land use, land-use change and forestry. Institute for Global Environmental Strategies (IGES) for the IPCC. 2003. https://www.ipcc-nggip.iges.or.jp/public/gpglulucf/gpglulucf_files/GPG_LULUCF_FULL.pdf.46.Waldner P. Detection of temporal trends in atmospheric deposition of inorganic nitrogen and sulphate to forests in Europe. Atm Env. 2014;95:363–74.CAS 

Google Scholar 
47.Nieminen T, De Vos B, Cools N, König N, Fischer R, Lost S, et al. Part XI: Soil Solution Collection and Analysis. In: UNECE ICP Forests Programme Co-ordinating Centre (eds): Manual on methods and criteria for harmonized sampling, assessment, monitoring and analysis of the effects of air pollution on forests. Thünen Institute of Forest Ecosystems. Eberswalde. 2016. https://www.icp-forests.org/pdf/manual/2016/ICP_Manual_2016_01_part11.pdf.48.Fick SE, Hijmans RJ. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Inter J Clim. 2017;37:4302–15.
Google Scholar 
49.Cox F, Barsoum N, Lilleskov EA, Bidartondo MI. Nitrogen availability is a primary determinant of conifer mycorrhizas across complex environmental gradients. Ecol Lett. 2010;13:1103–13.PubMed 

Google Scholar 
50.Okonechnikov K, Golosova O, Fursov M. Unipro GENE: a unified bioinformatics toolkit. Bioinf. 2012;28:1166–7.CAS 

Google Scholar 
51.Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinf. 2010;26:2460–1.CAS 

Google Scholar 
52.Abarenkov K, Nilsson RH, Larsson K-H, Alexander IJ, Eberhardt U, Erland S, et al. The UNITE database for molecular identification of fungi – recent updates and future perspectives. N. Phytol. 2010;186:281–5.
Google Scholar 
53.Grigoriev IV, Nikitin R, Haridas S, Kuo A, Ohm R, Otillar R, et al. MycoCosm portal: gearing up for 1000 fungal genomes. Nucleic Acids Res. 2014;42:D699–D704.CAS 
PubMed 

Google Scholar 
54.Douglas GM, Beiko RG, Langille MG Predicting the Functional Potential of the Microbiome from Marker Genes Using PICRUSt (eds). Microbiome Analysis. Methods in Molecular Biology. Vol 1849. Humana Press, New York, 2018. pp 169–77.55.Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat biot. 2013;31:814–21.CAS 

Google Scholar 
56.Treseder KK, Lennon JT. Fungal traits that drive ecosystem dynamics on land. Microbiol Mol Biol Rev. 2015;79:243–62.CAS 
PubMed 
PubMed Central 

Google Scholar 
57.Saifuddin M, Bhatnagar JM, Segrè D, Finzi AC. Microbial carbon use efficiency predicted from genome-scale metabolic models. Nat Commun. 2019;10:3568.PubMed 
PubMed Central 

Google Scholar 
58.Revell LJ, Revell MLJ Package ‘phytools’. 2020. https://github.com/liamrevell/phytools.59.Thompson LR, Sanders JG, McDonald D, Amir A, Ladau J, Locey KJ, et al. A communal catalogue reveals Earth’s multiscale microbial diversity. Nature 2017;551:457–63.CAS 
PubMed 
PubMed Central 

Google Scholar 
60.Gibbons SM, Lekberg Y, Mummey DL, Sangwan N, Ramsey PW, Gilbert JA. Invasive plants rapidly reshape soil properties in a grassland ecosystem. mSystems 2017;2:e00178–16.CAS 
PubMed 
PubMed Central 

Google Scholar 
61.Pold G, Domeignoz-Horta LA, Morrison EW, Frey SD, Sistla SA, DeAngelis KM. Carbon use efficiency and its temperature sensitivity covary in soil bacteria. MBio 2020;11:e02293–19.CAS 
PubMed 
PubMed Central 

Google Scholar 
62.Stewart JD, Shakya KM, Bilinski T, Wilson JW, Ravi S, Choi CS. Variation of near surface atmosphere microbial communities at an urban and a suburban site in Philadelphia, PA, USA. Sci Tot Env. 2020;724:138353.CAS 

Google Scholar 
63.Sun S, Jones RB, Fodor AA. Inference-based accuracy of metagenome prediction tools varies across sample types and functional categories. Microbiome 2020;8:46.PubMed 
PubMed Central 

Google Scholar 
64.Fierer N, Leff JW, Adams BJ, Nielsen UN, Bates ST, Lauber CL, et al. Cross-biome metagenomic analyses of soil microbial communities and their functional attributes. PNAS 2012;109:21390–5.CAS 
PubMed 
PubMed Central 

Google Scholar 
65.Moore JAM, Anthony MA, Pec GJ, Trocha LK, Trzebny A, Geyer KM, et al. Fungal community structure and function shifts with atmospheric nitrogen deposition. Glob Chan Bio. 2021;27:1349–64.
Google Scholar 
66.Team RC R: A language and environment for statistical computing. 2013.67.Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’hara R, et al. Package ‘vegan’. 2013. https://github.com/vegandevs/vegan.68.Morris MH, Pérez-Pérez MA, Smith ME, Bledsoe CS. Multiple species of ectomycorrhizal fungi are frequently detected on individual oak root tips in a tropical cloud forest. Mycorrhiza 2008;18:375–83.PubMed 

Google Scholar 
69.Lindner DL, Banik MT. Effects of cloning and root-tip size on observations of fungal ITS sequences from Picea glauca roots. Mycologia 2009;101:157–65.CAS 
PubMed 

Google Scholar 
70.Orlovich DA, Draffin SJ, Daly RA, Stephenson SL. Piracy in the high trees: ectomycorrhizal fungi from an aerial ‘canopy soil’ microhabitat. Mycologia 2013;105:52–60.PubMed 

Google Scholar 
71.Tedersoo L, Nilsson RH, Abarenkov K, Jairus T, Sadam A, Saar I, et al. 454 Pyrosequencing and Sanger sequencing of tropical mycorrhizal fungi provide similar results but reveal substantial methodological biases. N Phytol. 2010;188:291–301.CAS 

Google Scholar 
72.Morrison EW, Frey SD, Sadowsky JJ, van Diepen LT, Thomas WK, Pringle A. Chronic nitrogen additions fundamentally restructure the soil fungal community in a temperate forest. Fungal Ecol. 2016;23:48–57.
Google Scholar 
73.Paradis E, Schliep K. ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinf. 2019;35:526–8.CAS 

Google Scholar 
74.Madhulatha TS An overview on clustering methods. arXiv preprint. 2012;arXiv:1205.111775.Pallmann P, Hothorn LA. Analysis of means: a generalized approach using R. J Ap Stat 2016;43:1541–60.
Google Scholar 
76.De Caceres M, Jansen F, De Caceres MM Package ‘indicspecies’. 2016. https://vegmod.github.io/software/indicspecies.77.Wood S, Wood MS Package ‘mgcv’. 2015. https://cran.r-project.org/web/packages/mgcv/mgcv.pdf.78.Larsen WA, McCleary SJ. The use of partial residual plots in regression analysis. Technometrics 1972;14:781–90.
Google Scholar 
79.Gower ST, McMurtrie RE, Murty D. Aboveground net primary production decline with stand age: potential causes. Tr Eco Evol. 1996;11:378–82.CAS 

Google Scholar 
80.O’brien RM. A caution regarding rules of thumb for variance inflation factors. Qual Quant. 2007;41:673–90.
Google Scholar 
81.Koide RT, Fernandez CW. The continuing relevance of “older” mycorrhiza literature: insights from the work of John Laker Harley (1911–1990). Mycorrhiza 2018;28:577–86.CAS 
PubMed 

Google Scholar 
82.Anthony MA, Stinson KA, Moore JAM, Frey SD. Plant invasion impacts on fungal community structure and function depend on soil warming and nitrogen enrichment. Oecologia 2020;194:659–72.CAS 
PubMed 
PubMed Central 

Google Scholar 
83.Jonsson LM, Nilsson M-C, Wardle DA, Zackrisson O. Context dependent effects of ectomycorrhizal species richness on tree seedling productivity. Oikos 2001;93:353–64.
Google Scholar 
84.Hazard C, Kruitbos L, Davidson H, Taylor AFS, Johnson D. Contrasting effects of intra- and interspecific identity and richness of ectomycorrhizal fungi on host plants, nutrient retention and multifunctionality. N. Phytol. 2017;213:852–63.CAS 

Google Scholar 
85.Gehring CA, Sthultz CM, Flores-Rentería L, Whipple AV, Whitham TG. Tree genetics defines fungal partner communities that may confer drought tolerance. PNAS 2017;114:11169–74.CAS 
PubMed 
PubMed Central 

Google Scholar 
86.Liang J, Crowther TW, Picard N, Wiser S, Zhou M, Alberti G, et al. Positive biodiversity-productivity relationship predominant in global forests. Science. 2016;354.87.Baxter JW, Dighton J. Ectomycorrhizal diversity alters growth and nutrient acquisition of grey birch (Betula populifolia) seedlings in host–symbiont culture conditions. N. Phytol. 2001;152:139–49.
Google Scholar 
88.Dighton J, White JF. The fungal community: its organization and role in the ecosystem. 3rd ed. CRC Press, Boca Raton, 2005.89.Diagne N, Thioulouse J, Sanguin H, Prin Y, Krasova-Wade T, Sylla S, et al. Ectomycorrhizal diversity enhances growth and nitrogen fixation of Acacia mangium seedlings. Soil Bio Biochem. 2013;57:468–76.CAS 

Google Scholar 
90.Köhler J, Yang N, Pena R, Raghavan V, Polle A, Meier IC. Ectomycorrhizal fungal diversity increases phosphorus uptake efficiency of European beech. N. Phytol. 2018;220:1200–10.
Google Scholar 
91.Nygren CMR, Eberhardt U, Karlsson M, Parrent JL, Lindahl BD, Taylor AFS. Growth on nitrate and occurrence of nitrate reductase-encoding genes in a phylogenetically diverse range of ectomycorrhizal fungi. N. Phytol. 2008;180:875–89.CAS 

Google Scholar 
92.Wallenda T, Stober C, Högbom L, Schinkel H, George E, Högberg P, et al. Nitrogen Uptake Processes in Roots and Mycorrhizas (eds). Carbon and Nitrogen Cycling in European Forest Ecosystems. Springer, Berlin, 2000. pp 122–43.93.Jilling A, Keiluweit M, Contosta AR, Frey S, Schimel J, Schnecker J, et al. Minerals in the rhizosphere: overlooked mediators of soil nitrogen availability to plants and microbes. Biogeoch. 2018;139:103–22.CAS 

Google Scholar 
94.Marzluf GA. Regulation of nitrogen metabolism and gene expression in fungi. Microbi Rev. 1981;45:437–61.CAS 

Google Scholar 
95.Sinsabaugh RL, Moorhead DL. Resource allocation to extracellular enzyme production: a model for nitrogen and phosphorus control of litter decomposition. Soil Bio Biochem. 1994;26:1305–11.
Google Scholar 
96.Bödeker ITM, Clemmensen KE, de Boer W, Martin F, Olson Å, Lindahl BD. Ectomycorrhizal Cortinarius species participate in enzymatic oxidation of humus in northern forest ecosystems. N. Phytol. 2014;203:245–56.
Google Scholar 
97.Lilleskov E, Hobbie EA, Horton T. Conservation of ectomycorrhizal fungi: exploring the linkages between functional and taxonomic responses to anthropogenic N deposition. Fungal Ecol. 2011;4:174–83.
Google Scholar 
98.Franklin O, Näsholm T, Högberg P, Högberg MN. Forests trapped in nitrogen limitation – an ecological market perspective on ectomycorrhizal symbiosis. N. Phytol. 2014;203:657–66.CAS 

Google Scholar 
99.Rocca JD, Hall EK, Lennon JT, Evans SE, Waldrop MP, Cotner JB, et al. Relationships between protein-encoding gene abundance and corresponding process are commonly assumed yet rarely observed. ISME J. 2015;9:1693–9.PubMed 

Google Scholar 
100.Põlme S, Abarenkov K, Henrik Nilsson R, Lindahl BD, Clemmensen KE, Kauserud H, et al. FungalTraits: a user-friendly traits database of fungi and fungus-like stramenopiles. Fungal Div. 2020;105:1–16.
Google Scholar 
101.Ekblad A, Wallander H, Godbold DL, Cruz C, Johnson D, Baldrian P, et al. The production and turnover of extramatrical mycelium of ectomycorrhizal fungi in forest soils: role in carbon cycling. Plant Soil. 2013;366:1–27.CAS 

Google Scholar 
102.Agerer R. Exploration types of ectomycorrhizae. Mycorrhiza 2001;11:107–14.
Google Scholar 
103.Suz LM, Bidartondo MI, van der Linde S, Kuyper TW. Ectomycorrhizas and tipping points in forest ecosystems. N. Phytol 2021;231:1700–7.
Google Scholar 
104.Wasyliw J, Karst J. Shifts in ectomycorrhizal exploration types parallel leaf and fine root area with forest age. J Ecol. 2020;108:2270–82.CAS 

Google Scholar 
105.LeDuc SD, Lilleskov EA, Horton TR, Rothstein DE. Ectomycorrhizal fungal succession coincides with shifts in organic nitrogen availability and canopy closure in post-wildfire jack pine forests. Oecologia 2013;172:257–69.PubMed 

Google Scholar 
106.Struck C. Amino acid uptake in rust fungi. Front Plant Sci. 2015;6:40.PubMed 
PubMed Central 

Google Scholar 
107.Wen Z, Shi L, Tang Y, Shen Z, Xia Y, Chen Y. Effects of Pisolithus tinctorius and Cenococcum geophilum inoculation on pine in copper-contaminated soil to enhance phytoremediation. Int J Phytorem. 2017;19:387–94.CAS 

Google Scholar 
108.Garbaye J, Churin J-L. Effect of ectomycorrhizal inoculation at planting on growth and foliage quality of Tilia tomentosa. J Arbor 1996;22:29–34.
Google Scholar 
109.Fernandez CW, Koide RT. The function of melanin in the ectomycorrhizal fungus Cenococcum geophilum under water stress. Fungal Ecol. 2013;6:479–86.
Google Scholar 
110.Heinonsalo J, Sun H, Santalahti M, Bäcklund K, Hari P, Pumpanen J. Evidences on the ability of mycorrhizal genus Piloderma to use organic nitrogen and deliver it to Scots Pine. PLoS ONE. 2015;10:e0131561.PubMed 
PubMed Central 

Google Scholar 
111.Tedersoo L, Bahram M, Põlme S, Kõljalg U, Yorou NS, Wijesundera R, et al. Global diversity and geography of soil fungi. Science. 2014; 346.112.Polley H, Kroiher F, Riedel T Beech and spruce popular and in-demand. Thünen-Institut, Bundesforschungsinstitut für Ländliche Räume, Wald und Fischerei. 2015. https://literatur.thuenen.de/digbib_extern/dn055748.pdf.113.Brzostek ER, Fisher JB, Phillips RP. Modeling the carbon cost of plant nitrogen acquisition: mycorrhizal trade-offs and multipath resistance uptake improve predictions of retranslocation. J Geophy Res. 2014;119:1684–97.
Google Scholar 
114.Ley RE, Backhed F, Turnbaugh P, Lozupone CA, Knight RD, Gordon JI. Obesity alters gut microbial ecology. PNAS 2005;102:11070–5.CAS 
PubMed 
PubMed Central 

Google Scholar 
115.NIH Human Microbiome Portfolio Analysis Team. A review of 10 years of human microbiome research activities at the US National Institutes of Health, Fiscal Years 2007-2016. Microbiome 2019;7:31.
Google Scholar  More

Trends in cephalopods tradeSince 2000, trade in fresh octopus has been constantly dominated by the flow from China to Korea, followed by Vietnam to Japan, Portugal to Spain and Spain to Italy. However, there has been a marked decrease in the traded volume and monetary value over time, with a 50% reduction in the top 5 traders (Tables S3–S5).Over the last 20 years, fresh octopus exports have been strongly dominated by China, followed by Spain, Vietnam, Portugal, and France, and recently by Morocco and Thailand (Table S3). While Vietnam was the most important exporter in the first period (2000–2005), it was not within the top 5 traders in the last 5 years. Imports have been dominated by Korea, Italy, and Portugal, with no notable changes in the whole period (Table S4). Regarding trade of processed octopus, the largest transactions have been performed from Morocco to Spain, Morocco to Japan, Mauritania to Japan (and more recently also to Spain) and China to Korea (Table S5). Since 2000, exports of processed products have been dominated by Morocco, Mauritania, China, Spain, and Vietnam (Table S3), while imports have been led by Japan, Spain, Italy, Korea, and the United States (Table S4).Trade in volume of fresh cuttlefish and squid includes fewer clear relationships over time, such as transactions from Malaysia to Singapore (2000–2004 and 2005–2009); from Myanmar to Thailand (2005–2009 and 2015–2019); and from Yemen to Vietnam (2010–2015) (Tables S6-S8). For the first 5 years, exports of fresh commodities were dominated by Vietnam. However, since 2005, India, Spain and France have increased their exports in both monetary value and volume, displacing Vietnam from the top rank (Table S6). The main importing traders were Spain and Italy. Although China was important in the first decade, it was replaced by Vietnam in the last decade (Table S7). The trade of elaborated cuttlefish and squid products was dominated by monetary value flow from Thailand to Japan and from Malvinas/Falkland Islands to Spain in the first decade, while in the last decade the flow from China to other traders (Japan, the USA, and Thailand) gained relevance. However, the volume follows a different pattern, with flow from the Malvinas/Falkland Islands to Spain and from Korea to China in the first decade, while in the last decade, flows from Peru to China and from China to Thailand were important (Table S8). A disparity exists between the top five traders in terms of flow of monetary value and volume in the first 15 years; but in the last 5 years the top positions are constantly represented by, China, Peru, India, and Spain. Italy, Japan, China, and the USA are important importers in terms of monetary value and volume, although in the last decade Thailand has increased its importance, replacing the USA in the top 5 in the last 5 years (Table S7).The CGTN involved 220 traders (countries or territories) from around the world with exports greater than or equal to 500 kg between 2000 and 2019 (Fig. 3). The remaining 32 traders either did not report exports or their exports were below 500 kg. The most important cluster of traders was composed by 8 countries that dominate the cephalopod global markets in Asia (China, India, Republic of Korea, Thailand, Vietnam), Europe (the Netherlands, Spain) and the USA. The second and third most relevant clusters were composed of 8 and 12 traders, respectively. These two clusters involve 9 developed countries (Belgium, Canada, Denmark, France, Germany, Italy, Japan, United Kingdom and Portugal) and 11 developing countries (Morocco, Malaysia, UAE, Senegal, South Africa, Peru, Indonesia, Philippines, Argentina, Chile and New Zealand). Some of these traders have the most productive cephalopod fisheries in the world (e.g., Patagonian shortfin squid in the Southwest Atlantic Ocean and Patagonian squid in the Southeastern Pacific Ocean).Figure 3The Cephalopod Global Trade Network. The top 220 traders of the CGTN as nodes (circles) and their trade links as lines. The colour and the size of the nodes represent, respectively, the cluster membership and relative importance of the trader in the CGTN, estimated from the number of trade links with other traders (i.e., degree). The colour of the edges represents the origin, destination and the proportion of trade links for all years between each pair of traders. The clusters were made using Ward’s method. The figure was created with R12 (https://cran.r-project.org) packages: “ggraph” v.2.0.030 (https://ggraph.data-imaginist.com) and “ggtree” v3.0.232 (https://guangchuangyu.github.io/ggtree-book/chapter-ggtree.html).Full size imageOctopus trade networkLive, fresh, or chilled octopusThe normalised strength (Fig. 4) revealed the importance of China and Republic of Korea in the trade of fresh octopus in monetary value, with high importance of flows between these two traders over time (Supplementary Fig. S1). Other relevant traders over time were Spain, Portugal and Italy, in Europe; and Vietnam and Japan, in Asia from 2000 to 2004 (Supplementary Fig. S1). The network based in volume showed similar results.Figure 4Global trade network for octopus live, fresh or chilled between 1 January 2000, and 31 December 2019 in monetary value (USD). The numbers correspond to the normalised strength for the monetary value. Each node represents a trader, and each edge represents the export–import relationship between two traders. The size and colour of the node represent the relative importance of the trader in the network in terms of its strength. The width and colour of the edge represent the relative importance of the relationship between two traders in terms of their edge strength. The figure was created with R12 (https://cran.r-project.org) packages: “ggplot2” v.3.2.113 (https://ggplot2.tidyverse.org), “ggmap” v.3.0.029 (https://github.com/dkahle/ggmap) and “ggraph” v.2.0.030 (https://ggraph.data-imaginist.com).Full size imageThe Betweenness identified important actors facilitating flow through the network. For fresh octopus, the most relevant traders in the last two decades were Spain, France, and Italy, followed by Thailand, Portugal and the USA. Again, no major differences exist between the monetary value and volume networks. However, the ranking of traders changed over time, with Italy replacing Spain in the first place during the period 2005–2009, and then Spain consolidating again in the first place in the following periods. Also, in the last 20 years there were many changes in Asia, with Vietnam losing and Korea gaining prominence over time (Supplementary Fig. S2).In a global trade network, there are countries that are essential to the network structure because they are connected to other countries critical to the network and those critical countries, in turn, have no other significant connections. PageRank is a centrality measure that identifies these important countries, resulting from an iterative algorithm that assigns higher values to countries with a greater number of import connections with other countries that move large quantities of goods or money33. In the last 20 years, Italy, Germany France and Spain have occupied those central positions in the global trade network of live, fresh or chilled octopus. Their dominance has not changed over the four 5-year periods analysed (Fig. 5).Figure 5Global trade network for octopus live, fresh or chilled between 1 January 2000, and 31 December 2019 in monetary value (USD). The numbers correspond to the normalised PageRank for the monetary value. Each node represents a trader. The size and colour of the node represent the relative importance of the trader in the network in terms of its PageRank. The figure was created with R12 (https://cran.r-project.org) packages: “ggplot2” v.3.2.113 (https://ggplot2.tidyverse.org), “ggmap” v.3.0.029 (https://github.com/dkahle/ggmap) and “ggraph” v.2.0.030 (https://ggraph.data-imaginist.com).Full size imageElaborated octopusThe normalised strength revealed a diversified trade network for elaborated octopus products. This network has remained relatively stable over the last 20 years, with some exceptions. In the 2000s, there was an intense flow of exports from North Africa (i.e., Morocco and Mauritania) to Japan, although the most important flow was from Morocco to Spain. In the 2010s, Mauritania changed its preferential partner and exported large quantities of elaborated octopus to Spain, with the latter gaining dominance. In this last decade, flows from China and Vietnam to the Republic of Korea also became important. Other relevant actors were distributed globally (e.g., Italy, Portugal, Senegal, the USA). However, the most important routes showed a common pattern: the origin was in developing countries or territories (that emerged as producers) while developed countries showed a high and stable consumer demand (Supplementary Fig. S3). The network based on volume was highly similar to the monetary value network. However, Italy, China, Korea, Vietnam, and the USA reduced their importance compared to the top-ranked traders (i.e., Spain, Japan, and Morocco). The most important routes of the volume network were from China to Korea; Morocco to Spain; Morocco and Mauritania to Japan and Vietnam to Korea.The Betweenness measure highlighted the role of Spain as a facilitating actor in the trade network of elaborated octopus, followed by Italy, China, and the USA. These countries have maintained their importance as structurers of the world trade network in elaborated octopus over the past 20 years. Similarly, the routes from Italy to Spain, and from Spain to China and the USA emerged as relevant in the network structure, with a special mention to the route between Japan and China in the last 5-year period (Supplementary Fig. S4). There are no major differences between the most central traders in this network and the volume-based one.PageRank revealed the importance of Spain and Italy as leading traders in the elaborated octopus market. These two countries concentrated a large number of import relationships, which also concentrated a large monetary and volume flow. This importance has been maintained over time. Other reference actors were Greece, Japan, the USA and Portugal, in the first decade; and the USA, Portugal, Greece and Korea in the second decade. Note how in the second decade, Greece and Japan lost relevance, while Northern European countries and the Republic of Korea gained relevance (Supplementary Fig. S5).Squid and cuttlefish trade networkLive, fresh, or chilled squid and cuttlefishThe normalised strength revealed the importance of Spain, France, Italy, and India in the trade network of fresh squid and cuttlefish products, especially the route between east Asia and Spain (Fig. 6). The volume-based network is highly similar to the monetary value network. Over the four 5-year periods analysed, Vietnam and Japan have gradually lost relative importance in the network (Supplementary Fig. S6).Figure 6Global trade network for squid and cuttlefish live, fresh or chilled between 1 January 2000, and 31 December 2019 in monetary value (USD). The numbers correspond to the normalised strength for the monetary value. Each node represents a trader, and each edge represents the export–import relationship between two traders. The size and colour of the node represent the relative importance of the trader in the network in terms of its strength. The width and colour of the edge represent the relative importance of the relationship between two traders in terms of their edge strength. The figure was created with R12 (https://cran.r-project.org) packages: “ggplot2” v.3.2.113 (https://ggplot2.tidyverse.org), “ggmap” v.3.0.029 (https://github.com/dkahle/ggmap) and “ggraph” v.2.0.030 (https://ggraph.data-imaginist.com).Full size imageFor fresh squid and cuttlefish, Betweenness identified Spain as the most important structurer of the network, in both monetary value and volume, over the time. Europe emerged as a major region structuring the global trade network. While in Asia, the exchange of countries with higher betweenness over time evidences the strong struggle for control of trade in the region and the more fragile sub-networks. In the monetary network, in the 2000s the bridge between the United Kingdom and Korea stood out, while in the 2010s the Europe-Asia bridge was established between Spain and India. In the volume-based network, it is noteworthy that in the 2010s there were many critical routes for the stability of the network, even in the Europe–Asia connection. Note, for example how the link between Netherlands and Myanmar stand out (Supplementary Fig. S7).PageRank revealed the importance of Italy, Spain, Germany and France as leading traders in the fresh squid and cuttlefish market. In the last 20 years, Europe concentrated the largest number of import relationships, which also concentrated a large monetary and volume flow. Europe leadership has been maintained over time (Fig. 7).Figure 7Global trade network for squid and cuttlefish live, fresh or chilled between 1 January 2000, and 31 December 2019 in monetary value (USD). The numbers correspond to the normalised PageRank for the monetary value. Each node represents a trader. The size and colour of the node represent the relative importance of the trader in the network in terms of its PageRank. The figure was created with R12 (https://cran.r-project.org) packages: “ggplot2” v.3.2.113 (https://ggplot2.tidyverse.org), “ggmap” v.3.0.029 (https://github.com/dkahle/ggmap) and “ggraph” v.2.0.030 (https://ggraph.data-imaginist.com).Full size imageElaborated squid and cuttlefishThe trade networks based on monetary value and volume for elaborated squid and cuttlefish emerge as global and complex, where several far distant traders have relevant roles in the import/export network (Supplementary Fig. S8). Although the most important nodes in the volume-based network reflected important nodes in the monetary value network, the strengths of the links, i.e., the flow of value and volume, did not. For example, in the volume-based network, Peru exported the largest quantities of squid and cuttlefish to China (Supplementary Fig. S8b,d), but the flow of money for these transactions was less important (Supplementary Fig. S8a,c).The betweenness centrality metric (based both on monetary value and volume) showed the importance of China, the USA, and Spain (followed by Italy, Korea, and Thailand) as facilitators in the elaborated goods trade network (Supplementary Fig. S9). While the main important bridges in volume transactions were between Italy and Spain, Spain and China, and China and the USA (Supplementary Fig. S9b,d), the main monetary bridges were from the USA to China, followed by the routes from Spain to the USA and from Italy to Spain (Supplementary Fig. S9a,c). The key traders structuring the network were the same, but they follow different directions.Closeness centrality highlighted the main actors in a regional context (Fig. 8). In both the monetary value and volume-based networks, China, North and South Korea, India, Indonesia, Thailand, and Vietnam form a strong trade network for squid and cuttlefish elaborated in Asia. Key players include South America (Peru, Argentina, Chile, the Malvinas/Falkland Islands); the USA; the Mediterranean (Morocco, Spain); Africa (South Africa, Mauritania); and the West Pacific region (New Zealand, Japan). Note how the highest values of closeness were slightly different in the money-based network (Fig. 8a) and in the volume-based network (Fig. 8b), mainly for those countries that are historically large producers of elaborated squid and cuttlefish (e.g., Peru and Argentina).Figure 8Global trade network for squid and cuttlefish elaborated between 1 January 2000, and 31 December 2019 in monetary value (USD) above, and volume (kg) below. The numbers correspond to the normalised closeness for the monetary value (USD) and volume (kg) traded, respectively. Each node represents a trader. The size and colour of the node represent the relative importance of the trader in the network in terms of its closeness. The figure was created with R12 (https://cran.r-project.org) packages: “ggplot2” v.3.2.113 (https://ggplot2.tidyverse.org), “ggmap” v.3.0.029 (https://github.com/dkahle/ggmap) and “ggraph” v.2.0.030 (https://ggraph.data-imaginist.com).Full size imagePageRank revealed the importance of Greece, Spain, China, USA and Japan as leading traders in the global elaborated squid and cuttlefish market for the first 10-years (Supplementary Fig. S10a). Globally, in the last decade, China, the USA and Japan have declined in importance, while Europe has consolidated its importance. However, Spain’s relevance has declined over the last decade while Germany’s has increased (Supplementary Fig. S10b). More

1.Chin, A. Z. et al. Malaria elimination in Malaysia and the rising threat of Plasmodium knowlesi. J. Physiol. Anthropol. 39, 36 (2020).PubMed 
PubMed Central 

Google Scholar 
2.Vythilingam, I., Wong, M. L. & Wan-Yussof, W. S. Current status of Plasmodium knowlesi vectors: A public health concern?. Parasitol 145, 32–40 (2018).CAS 

Google Scholar 
3.World Health Organization. Global technical strategy for malaria 2016–2030. World Health Org (2015).4.Sallum, M. A. M. et al. Vector competence, vectorial capacity of Nyssorhynchus darlingi and the basic reproduction number of Plasmodium vivax in agricultural settlements in the Amazonian region of Brazil. Malar. J. 18, 117 (2019).PubMed 
PubMed Central 

Google Scholar 
5.Garrett-Jones, C. The human blood index of malaria vectors in relation to epidemiological assessment. Bull. World Health Org. 30, 241 (1964).CAS 
PubMed 
PubMed Central 

Google Scholar 
6.De Oliveira, C. D., Tadei, W. P., Abdalla, F. C., Paolucci Pimenta, P. F. & Marinotti, O. Multiple blood meals in Anopheles darlingi (Diptera: Culicidae). J. Vector Ecol. 37, 351–358 (2012).PubMed 

Google Scholar 
7.Subbarao, S. K. Anopheline Species Complexes in Southeast Asia (WHO Technical Publication, 1998).
Google Scholar 
8.Hansen, I. A., Attardo, G. M., Rodriguez, S. D. & Drake, L. L. Four-way regulation of mosquito yolk protein precursor genes by juvenile hormone-, ecdysone-, nutrient-, and insulin-like peptide signaling pathways. Front. Physiol. 5, 103 (2014).PubMed 
PubMed Central 

Google Scholar 
9.Bashar, K., Tuno, N., Ahmed, T. & Howlader, A. Blood-feeding patterns of Anopheles mosquitoes in a malaria-endemic area of Bangladesh. Parasites Vectors 5, 39 (2012).PubMed 
PubMed Central 

Google Scholar 
10.Escobar, D. et al. Blood meal sources of Anopheles spp. in malaria endemic areas of Honduras. Insects 11, 1–12 (2020).
Google Scholar 
11.Boreham, P. F. L. & Garrett-Jones, C. Prevalence of mixed blood meals and double feeding in a malaria vector (Anopheles sacharovi Favre). Bull. World Health Organ. 48, 605 (1973).CAS 
PubMed 
PubMed Central 

Google Scholar 
12.Kent, R. J. & Norris, D. E. Identification of mammalian blood meals in mosquitoes by a multiplexed polymerase chain reaction targeting cytochrome B. Am. J. Trop. Med. Hyg. 73, 336–342 (2005).CAS 
PubMed 

Google Scholar 
13.Brown, R., Hing, C. T., Fornace, K. & Ferguson, H. M. Evaluation of resting traps to examine the behaviour and ecology of mosquito vectors in an area of rapidly changing land use in Sabah, Malaysian Borneo. Parasites Vectors 11, 1–15 (2018).
Google Scholar 
14.Brant, H. L. et al. Vertical stratification of adult mosquitoes (Diptera: Culicidae) within a tropical rainforest in Sabah, Malaysia. Malar. J. 15, 370 (2016).PubMed 
PubMed Central 

Google Scholar 
15.Young, K. I. et al. Identification of mosquito bloodmeals collected in diverse habitats in Malaysian borneo using COI barcoding. Trop. Med. Infect. Dis. 5, 51 (2020).PubMed Central 

Google Scholar 
16.Rohani, A. et al. Comparative human landing catch and CDC light trap in mosquito sampling in knowlesi malaria endemic areas in Peninsula Malaysia. Adv. Entomol. 04, 1–10 (2016).
Google Scholar 
17.Jeyaprakasam, N. K. et al. Evaluation of Mosquito Magnet and other collection tools for Anopheles mosquito vectors of simian malaria. Parasites Vectors 14, 1–13 (2021).
Google Scholar 
18.Detinova, T. S. Age-grouping methods in Diptera of medical importance with special reference to some vectors of malaria. Monogr. Ser. World Health Organ. 47, 13–191 (1962).CAS 
PubMed 

Google Scholar 
19.Brugman, V. A. et al. Blood-feeding patterns of native mosquitoes and insights into their potential role as pathogen vectors in the Thames estuary region of the United Kingdom. Parasites Vectors 10, 1–12 (2017).
Google Scholar 
20.WHO. Training module on malaria control: Malaria entomology and vector control. Guide for participants. World Health Organisation (2013).21.Derek Charlwood, J. et al. ‘Nature or nurture’: Survival rate, oviposition interval, and possible gonotrophic discordance among South East Asian anophelines. Malar. J. 15, 356 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
22.Scott, T. W. et al. Detection of multiple blood feeding in Aedes aegypti (Diptera: Culicidae) during a single gonotrophic cycle using a histologic technique. J. Med. Entomol. 30, 94–99 (1993).CAS 
PubMed 

Google Scholar 
23.Ramasamy, M. S., Srikrishnaraj, K. A., Hadjirin, N., Perera, S. & Ramasamy, R. Physiological aspects of multiple blood feeding in the malaria vector Anopheles tessellatus. J. Insect Physiol. 46, 1051–1059 (2000).CAS 
PubMed 

Google Scholar 
24.Fernandes, L. & Briegel, H. Reproductive physiology of Anopheles gambiae and Anopheles atroparvus. J. Vector Ecol. 30, 11 (2005).PubMed 

Google Scholar 
25.Nirmala, X., Marinotti, O. & James, A. A. The accumulation of specific mRNAs following multiple blood meals in Anopheles gambiae. Insect Mol. Biol. 14, 95–103 (2005).CAS 
PubMed 

Google Scholar 
26.Briegel, H. & Horler, E. Multiple blood meals as a reproductive strategy in Anopheles (Diptera: Culcidae). J. Med. Entomol. 30, 975–985 (1993).CAS 
PubMed 

Google Scholar 
27.Klowden, M. J. & Lea, A. O. Effect of defensive host behavior on the blood meal size and feeding success of natural populations of mosquitoes (Diptera: Culicidae). J. Med. Entomol. 15, 514–517 (1979).CAS 
PubMed 

Google Scholar 
28.Briegel, H. Fecundity, metabolism, and body size in Anopheles (Diptera: Culicidae), vectors of malaria. J. Med. Entomol. 27, 839–850 (1990).CAS 
PubMed 

Google Scholar 
29.Service, M. W. Effects of wind on the behaviour and distribution of mosquitoes and blackflies. Int. J. Biometeorol. 24, 347–353 (1980).ADS 

Google Scholar 
30.Rahman, W. A., Che’Rus, A. & Ahmad, A. H. Malaria and Anopheles mosquitos in Malaysia. Southeast Asian J. Trop. Med. Public Health 28, 599–605 (1997).CAS 
PubMed 

Google Scholar 
31.Sandosham, A. A. & Thomas, V. Malariology: With Special Reference to Malaya (Singapore University Press, 1983).
Google Scholar 
32.Self, L.S. & De Datta, S.K. The impact of water management practices in rice production on mosquito vector propagation. In Vector-Borne Disease Control in Humans Through Rice Agroecosystem Management: Proceedings of the Workshop on Research and Training Needs in the Field of Integrated Vector-Borne Disease Control in Riceland Agroecosystems of Developing Countries. 67 (1988).33.Gooding, R. H. Digestive processes of haematophagous insects I. A literature review. Quaest. Entomol. 8, 5–60 (1972).
Google Scholar 
34.Hocking, K. S. & Macinnes, D. G. Notes on the bionomics of Anopheles gambiae and A funestus in East Africa. Bull. Entomol. Res. 39, 453–465 (1948).CAS 
PubMed 

Google Scholar 
35.Jiram, A. I. et al. Entomologic investigation of Plasmodium knowlesi vectors in Kuala Lipis, Pahang, Malaysia. Malar. J. 11, 1 (2012).
Google Scholar 
36.Suwonkerd, W. et al. Vector biology and malaria transmission in Southeast Asia. In: S. Manguin (Ed.), Anopheles Mosquitoes—New Insights into Malaria Vectors (IntechOpen, 2013).37.Reid, J. A. Anopheline mosquitoes of Malaya and Borneo. Studies from the Institute for Medical Research, Malaysia. Anopheline mosquitoes Malaya Borneo. Stud. Inst. Med. Res. Malaysia 31, 520 (1968).
Google Scholar 
38.Vythilingam, I. et al. Plasmodium knowlesi malaria an emerging public health problem in Hulu Selangor, Selangor, Malaysia (2009–2013): Epidemiologic and entomologic analysis. Parasites Vectors 7, 1–14 (2014).
Google Scholar 
39.Garrett-Jones, C., Boreham, P. F. L. & Pant, C. P. Feeding habits of Anophelines (Diptera: Culicidae) in 1971–78, with reference to the human blood index: A review. Bull. Entomol. Res. 70, 165–185 (2017).
Google Scholar 
40.Lyimo, I. N. & Ferguson, H. M. Ecological and evolutionary determinants of host species choice in mosquito vectors. Trends Parasitol. 25, 189–196 (2009).PubMed 

Google Scholar 
41.Molaei, G., Andreadis, T. G., Armstrong, P. M. & Diuk-Wasser, M. Host-feeding patterns of potential mosquito vectors in Connecticut, USA: molecular analysis of bloodmeals from 23 species of Aedes, Anopheles, Culex, Coquillettidia, Psorophora, and Uranotaenia. J. Med. Entomol. 45, 1143–1151 (2008).CAS 
PubMed 

Google Scholar 
42.Ang, J. X. et al. New vectors in northern Sarawak, Malaysian Borneo, for the zoonotic malaria parasite, Plasmodium knowlesi. Parasites Vectors 13, 472 (2020).PubMed 
PubMed Central 

Google Scholar 
43.Ang, J. X. et al. New vectors that are early feeders for Plasmodium knowlesi and other simian malaria parasites in Sarawak, Malaysian Borneo. Sci. Rep. 11, 7739 (2021).ADS 
PubMed 
PubMed Central 

Google Scholar 
44.Stark, D. J. et al. Long-tailed macaque response to deforestation in a Plasmodium knowlesi-endemic area. EcoHealth 16, 638–646 (2019).PubMed 
PubMed Central 

Google Scholar 
45.Jeyaprakasam, N. K., Liew, J. W. K., Low, V. L., Wan-Sulaiman, W.-Y.Y. & Vythilingam, I. Plasmodium knowlesi infecting humans in Southeast Asia: What’s next?. PLoS Negl. Trop. Dis. 14, 1–16 (2020).
Google Scholar 
46.Amir, A., Sum, J. S., Lau, Y. L., Vythilingam, I. & Fong, M. Y. Colonization of Anopheles cracens: A malaria vector of emerging importance. Parasites Vectors 6, 2–5 (2013).
Google Scholar 
47.Andolina, C. et al. The suitability of laboratory-bred Anopheles cracens for the production of Plasmodium vivax sporozoites. Malar. J. 14, 312 (2015).PubMed 
PubMed Central 

Google Scholar 
48.Rund, S. S. C., O’Donnell, A. J., Gentile, J. E. & Reece, S. E. Daily rhythms in mosquitoes and their consequences for malaria transmission. Insects. 7, 14 (2016).PubMed Central 

Google Scholar 
49.Aviles, E. I., Rotenberry, R. D., Collins, C. M., Dotson, E. M. & Benedict, M. Q. Fluorescent markers rhodamine B and uranine for Anopheles gambiae adults and matings. Malar. J. 19, 1–9 (2020).
Google Scholar 
50.Alcaide, M. et al. Disentangling vector-borne transmission networks: A universal DNA barcoding method to identify vertebrate hosts from arthropod bloodmeals. PLoS ONE 4, e7092 (2009).ADS 
PubMed 
PubMed Central 

Google Scholar 
51.Gunathilaka, N., Denipitiya, T., Hapugoda, M., Abeyewickreme, W. & Wickremasinghe, R. Determination of the foraging behaviour and blood meal source of malaria vector mosquitoes in Trincomalee district of Sri Lanka using a multiplex real time polymerase chain reaction assay. Malar. J. 15, 1–10 (2016).
Google Scholar 
52.Sallum, M. A. M., Peyton, E. L., Harrison, B. A. & Wilkerson, R. C. Revision of the Leucosphyrus group of Anopheles (Cellia) (Diptera, Culicidae). Rev. Bras. Entomol. 49, 1–152 (2005).
Google Scholar 
53.Sum, J. S. et al. Phylogenetic study of six species of Anopheles mosquitoes in Peninsular Malaysia based on inter-transcribed spacer region 2 (ITS2) of ribosomal DNA. Parasites Vectors 7, 309 (2014).PubMed 
PubMed Central 

Google Scholar 
54.Pappa, V., Reddy, M., Overgaard, H. J., Abaga, S. & Caccone, A. Short report: Estimation of the human blood index in malaria mosquito vectors in Equatorial Guinea after indoor antivector interventions. Am. J. Trop. Med. Hyg. 84, 298–301 (2011).PubMed 
PubMed Central 

Google Scholar  More

1.Scheu, S. Plants and generalist predators as links between the below-ground and above-ground system. Basic Appl. Ecol. 2, 3–13 (2001).
Google Scholar 
2.Wardle, D. A. et al. Ecological linkages between aboveground and belowground biota. Science 304, 1629–1633 (2004).CAS 
PubMed 
ADS 

Google Scholar 
3.Van Der Putten, W. H. et al. Empirical and theoretical challenges in aboveground-belowground ecology. Oecologia 161, 1–14 (2009).PubMed 
PubMed Central 
ADS 

Google Scholar 
4.Miyashita, T., Takada, M. & Shimazaki, A. Experimental evidence that aboveground predators are sustained by underground detritivores. Oikos 103, 31–36 (2003).
Google Scholar 
5.Moore, J. C. et al. Detritus, trophic dynamics and biodiversity. Ecol. Lett. 7, 584–600 (2004).ADS 

Google Scholar 
6.Haraguchi, T. F., Uchida, M., Shibata, Y. & Tayasu, I. Contributions of detrital subsidies to aboveground spiders during secondary succession, revealed by radiocarbon and stable isotope signatures. Oecologia 171, 935–944 (2013).PubMed 
ADS 

Google Scholar 
7.Halaj, J. & Wise, D. H. Impact of a detrital subsidy on trophic cascades in a terrestrial grazing food web. Ecology 83, 3141 (2002).
Google Scholar 
8.Zakharov, A. A. Ants at the border between the epigean and soil blocks of a forest coenosis. Rej 21, 219–222 (2012).
Google Scholar 
9.Iakovlev, I. K., Novgorodova, T. A., Tiunov, A. V. & Reznikova, Z. I. Trophic position and seasonal changes in the diet of the red wood ant Formica aquilonia as indicated by stable isotope analysis. Ecol. Entomol. 42, 263–272 (2017).
Google Scholar 
10.Chan, E. K. W., Zhang, Y. & Dudgeon, D. Arthropod ‘rain’ into tropical streams: The importance of intact riparian forest and influences on fish diets. Mar. Freshw. Res. 59, 653–660 (2008).
Google Scholar 
11.Pringle, R. M. & Fox-Dobbs, K. Coupling of canopy and understory food webs by ground-dwelling predators. Ecol. Lett. 11, 1328–1337 (2008).PubMed 

Google Scholar 
12.Rozanova, O. L., Tsurikov, S. M., Tiunov, A. V. & Semenina, E. E. Arthropod rain in a temperate forest: Intensity and composition. Pedobiologia 75, 52–56 (2019).
Google Scholar 
13.Goncharov, A. A., Tsurikov, S. M., Potapov, A. M. & Tiunov, A. V. Short-term incorporation of freshly fixed plant carbon into the soil animal food web: Field study in a spruce forest. Ecol. Res. 31, 923–933 (2016).CAS 

Google Scholar 
14.Potapov, A. M., Goncharov, A. A., Tsurikov, S. M., Tully, T. & Tiunov, A. V. Assimilation of plant-derived freshly fixed carbon by soil collembolans: Not only via roots?. Pedobiologia 59, 189–193 (2016).
Google Scholar 
15.Hyodo, F., Kohzu, A. & Tayasu, I. Linking aboveground and belowground food webs through carbon and nitrogen stable isotope analyses. Ecol. Res. 25, 745–756 (2010).CAS 

Google Scholar 
16.Potapov, A. M., Tiunov, A. V. & Scheu, S. Uncovering trophic positions and food resources of soil animals using bulk natural stable isotope composition. Biol. Rev. 94, 37–59 (2019).
Google Scholar 
17.Potapov, A. M., Semenina, E. E., Kurakov, A. V. & Tiunov, A. V. Large 13C/12C and small 15N/14N isotope fractionation in an experimental detrital foodweb (litter-fungi-collembolans). Ecol. Res. 28, 1069–1079 (2013).CAS 

Google Scholar 
18.McCutchan, J. H., Lewis, W. M., Kendall, C. & McGrath, C. C. Variation in trophic shift for stable isotope ratios of carbon, nitrogen, and sulfur. Oikos 102, 378–390 (2003).CAS 

Google Scholar 
19.Hyodo, F. Use of stable carbon and nitrogen isotopes in insect trophic ecology. Entomol. Sci. 18, 295–312 (2015).
Google Scholar 
20.Potapov, A. M., Korotkevich, A. Y. & Tiunov, A. V. Non-vascular plants as a food source for litter-dwelling Collembola: Field evidence. Pedobiologia 66, 11–17 (2018).
Google Scholar 
21.Tozer, W. C., Hackell, D., Miers, D. B. & Silvester, W. B. Extreme isotopic depletion of nitrogen in New Zealand lithophytes and epiphytes; the result of diffusive uptake of atmospheric ammonia?. Oecologia 144, 628–635 (2005).CAS 
PubMed 
ADS 

Google Scholar 
22.Delgado, V., Ederra, A. & Santamaría, J. M. Nitrogen and carbon contents and δ15N and δ13C signatures in six bryophyte species: Assessment of long-term deposition changes (1980–2010) in Spanish beech forests. Glob. Chang. Biol. 19, 2221–2228 (2013).PubMed 
ADS 

Google Scholar 
23.Eskov, A. K. et al. Dependence of epiphytic community on autochthonous and allochthonous sources of nitrogen in three forest habitats of southern Vietnam. Plant Soil 443, 565–574 (2019).CAS 

Google Scholar 
24.Potapov, A. A., Semenina, E. E., Korotkevich, A. Y., Kuznetsova, N. A. & Tiunov, A. V. Connecting taxonomy and ecology: Trophic niches of collembolans as related to taxonomic identity and life forms. Soil Biol. Biochem. 101, 20–31 (2016).CAS 

Google Scholar 
25.Wallwork, J. A. Ecology of Soil Animals (McGraw-Hill, 1970).
Google Scholar 
26.Behan-Pelletier, V. & Winchester, N. Arboreal oribatid mite diversity: Colonizing the canopy. Appl. Soil Ecol. 9, 45–51 (1998).
Google Scholar 
27.Yoshida, T. & Hijii, N. Microarthropod colonization of litter in arboreal and soil environments of a Japanese cedar (Cryptomeria japonica) plantation. J. For. Res. 16, 46–54 (2011).
Google Scholar 
28.Korobushkin, D. I., Gongalsky, K. B. & Tiunov, A. V. Isotopic niche (δ13C and δ15N values) of soil macrofauna in temperate forests. Rapid Commun. Mass Spectrom. 28, 1303–1311 (2014).CAS 
PubMed 
ADS 

Google Scholar 
29.Post, D. M. et al. Getting to the fat of the matter: Models, methods and assumptions for dealing with lipids in stable isotope analyses. Oecologia 152, 179–189 (2007).PubMed 
ADS 

Google Scholar 
30.Spence, K. O. & Rosenheim, J. A. Isotopic enrichment in herbivorous insects: A comparative field-based study of variation. Oecologia 146, 89–97 (2005).PubMed 
ADS 

Google Scholar 
31.Southwood, T. R. E., Wint, G. R. W., Kennedy, C. E. J. & Greenwood, S. R. The composition of the arthropod fauna of the canopies of some species of oak (Quercus). Eur. J. Entomol. 102, 65–72 (2005).
Google Scholar 
32.Chahartaghi, M., Langel, R., Scheu, S. & Ruess, L. Feeding guilds in Collembola based on nitrogen stable isotope ratios. Soil Biol. Biochem. 37, 1718–1725 (2005).CAS 

Google Scholar 
33.Winchester, N. N., Behan-Pelletier, V. M. & Ring, R. A. Arboreal specificity, diversity and abundance of canopy-dwelling oribatid mites (Acari: Oribatida). Pedobiologia 43, 391–400 (1999).
Google Scholar 
34.Grove, S. J. Saproxylic insect ecology and the sustainable management of forests. Annu. Rev. Ecol. Syst. 33, 1–23 (2002).
Google Scholar 
35.Erdmann, G., Otte, V., Langel, R., Scheu, S. & Maraun, M. The trophic structure of bark-living oribatid mite communities analysed with stable isotopes (15N, 13C) indicates strong niche differentiation. Exp. Appl. Acarol. 41, 1–10 (2007).PubMed 

Google Scholar 
36.Klarner, B., Maraun, M. & Scheu, S. Trophic diversity and niche partitioning in a species rich predator guild—Natural variations in stable isotope ratios (13C/12C, 15N/14N) of mesostigmatid mites (Acari, Mesostigmata) from Central European beech forests. Soil Biol. Biochem. 57, 327–333 (2013).CAS 

Google Scholar 
37.Brooks, J. R., Flanagan, L. B., Buchmann, N. & Ehleringer, J. R. Carbon isotope composition of boreal plants: Functional grouping of life forms. Oecologia 110, 301–311 (1997).CAS 
PubMed 
ADS 

Google Scholar 
38.Pollierer, M. M., Langel, R., Körner, C., Maraun, M. & Scheu, S. The underestimated importance of belowground carbon input for forest soil animal food webs. Ecol. Lett. 10, 729–736 (2007).PubMed 

Google Scholar 
39.Pollierer, M. M., Langel, R., Scheu, S. & Maraun, M. Compartmentalization of the soil animal food web as indicated by dual analysis of stable isotope ratios (15N/14N and 13C/12C). Soil Biol. Biochem. 41, 1221–1226 (2009).CAS 

Google Scholar 
40.Tsurikov, S. M., Goncharov, A. A. & Tiunov, A. V. Intra-body variation and ontogenetic changes in the isotopic composition (13C/12C and 15N/14N) of beetles (Coleoptera). Entomol. Rev. 95, 326–333 (2015).
Google Scholar 
41.Ventura, M. & Jeppesen, E. Effects of fixation on freshwater invertebrate carbon and nitrogen isotope composition and its arithmetic correction. Hydrobiologia 632, 297–308 (2009).CAS 

Google Scholar 
42.Krab, E. J., Van Logtestijn, R. S. P., Cornelissen, J. H. C. & Berg, M. P. Reservations about preservations: Storage methods affect δ13C signatures differently even in closely related soil fauna. Methods Ecol. Evol. 3, 138–144 (2012).
Google Scholar 
43.Striganova, B. R. Nutrition of Soil Saprophages (Nauka, 1980) (In Russian).
Google Scholar 
44.Jackson, A. L., Inger, R., Parnell, A. C. & Bearhop, S. Comparing isotopic niche widths among and within communities: SIBER—Stable Isotope Bayesian Ellipses in R. J. Anim. Ecol. 80, 595–602 (2011).
Google Scholar 
45.R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2014).
Google Scholar  More

1.Luyster, F. S., Strollo, P. J., Zee, P. C. & Walsh, J. K. Sleep: A health imperative. Sleep 35, 727–734 (2012).PubMed 
PubMed Central 

Google Scholar 
2.Siegel, J. M. Clues to the functions of mammalian sleep. Nature 437, 1264–1271 (2005).CAS 
PubMed 
ADS 

Google Scholar 
3.Cirelli, C. & Tononi, G. Is sleep essential?. PLoS Biol. 6, 1605–1611 (2008).CAS 

Google Scholar 
4.Lesku, J. A., Roth, T. C. II., Amlaner, C. J. & Lima, S. L. A phylogenetic analysis of sleep architecture in mammals: The integration of anatomy, physiology, and ecology. Am. Nat. 168, 441–453 (2006).PubMed 

Google Scholar 
5.Banks, S. & Dinges, D. F. Behavioral and physiological consequences of sleep restriction. J. Clin. Sleep Med. 3, 519–528 (2007).PubMed 
PubMed Central 

Google Scholar 
6.Tougeron, K. & Abram, P. K. An ecological perspective on sleep disruption. Am. Nat. 190, E55–E66 (2017).PubMed 

Google Scholar 
7.Vyazovskiy, V. V. & Tobler, I. The temporal structure of behaviour and sleep homeostasis. PLoS ONE 7, 1–11 (2012).
Google Scholar 
8.Tobler, I. Is sleep fundamentally different between mammalian species?. Behav. Brain Res. 69, 35–41 (1995).CAS 
PubMed 

Google Scholar 
9.Vyazovskiy, V. Sleep, recovery, and metaregulation: Explaining the benefits of sleep. Nat. Sci. Sleep 7, 171 (2015).PubMed 
PubMed Central 

Google Scholar 
10.Vyazovskiy, V. V. & Delogu, A. NREM and REM sleep. Neuroscience 20, 203–219 (2014).
Google Scholar 
11.Orzeł-Gryglewska, J. Consequences of sleep deprivation. Int. J. Occup. Med. Environ. Health 23, 95–114 (2010).PubMed 

Google Scholar 
12.Meerlo, P., Sgoifo, A. & Suchecki, D. Restricted and disrupted sleep: Effects on autonomic function, neuroendocrine stress systems and stress responsivity. Sleep Med. Rev. 12, 197–210 (2008).PubMed 

Google Scholar 
13.Touitou, Y., Reinberg, A. & Touitou, D. Association between light at night, melatonin secretion, sleep deprivation, and the internal clock: Health impacts and mechanisms of circadian disruption. Life Sci. 173, 94–106 (2017).CAS 
PubMed 

Google Scholar 
14.Pires, G. N., Bezerra, A. G., Tufik, S. & Andersen, M. L. Effects of acute sleep deprivation on state anxiety levels: A systematic review and meta-analysis. Sleep Med. 68, 575–589 (2016).
Google Scholar 
15.Hudson, A. N., Van Dongen, H. P. A. A. & Honn, K. A. Sleep deprivation, vigilant attention, and brain function: A review. Neuropsychopharmacology https://doi.org/10.1038/s41386-019-0432-6 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
16.Tobler, I. & Sigg, H. Long-term motor activity recording of dogs and the effect of sleep deprivation. Experientia 42, 987–991 (1986).CAS 
PubMed 

Google Scholar 
17.Hänninen, L. Sleep and Rest in Calves Relationship To Welfare, Housing and Hormonal Activity (University of Helsinki, Helsinki, 2007).
Google Scholar 
18.Hsieh, W.-H.H. et al. Simulated shift work in rats perturbs multiscale regulation of locomotor activity. J. R. Soc. Interface 11, 20140318 (2014).PubMed 
PubMed Central 

Google Scholar 
19.Storch, C., Höhne, A., Holsboer, F. & Ohl, F. Activity patterns as a correlate for sleep–wake behaviour in mice. J. Neurosci. Methods 133, 173–179 (2004).PubMed 

Google Scholar 
20.Hicks, R. A., Moore, J. D., Hayes, C., Phillips, N. & Hawkins, J. REM sleep deprivation increases aggressiveness in male rats. Physiol. Behav. 22, 1097–1100 (1979).CAS 
PubMed 

Google Scholar 
21.Pires, G. N., Tufik, S. & Andersen, M. L. Grooming analysis algorithm: Use in the relationship between sleep deprivation and anxiety-like behavior. Prog. Neuro-Psychopharmacol. Biol. Psychiatry 41, 6–10 (2013).
Google Scholar 
22.Vorster, A. P. & Born, J. Sleep and memory in mammals, birds and invertebrates. Neurosci. Biobehav. Rev. 50, 103–119 (2015).PubMed 

Google Scholar 
23.Greives, T. J. et al. Costs of sleeping in: Circadian rhythms influence cuckoldry risk in a songbird. Funct. Ecol. 29, 1300–1307 (2015).
Google Scholar 
24.Fuchs, T., Haney, A., Jechura, T. J., Moore, F. R. & Bingman, V. P. Daytime naps in night-migrating birds: Behavioural adaptation to seasonal sleep deprivation in the Swainson’s thrush, Catharus ustulatus. Anim. Behav. 72, 951–958 (2006).
Google Scholar 
25.Klein, B. A., Klein, A., Wray, M. K., Mueller, U. G. & Seeley, T. D. Sleep deprivation impairs precision of waggle dance signaling in honey bees. Proc. Natl. Acad. Sci. U. S. A. 107, 22705–22709 (2010).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
26.Owczarczak-Garstecka, S. C. & Burman, O. H. P. P. Can sleep and resting behaviours be used as indicators of welfare in shelter dogs (Canis lupus familiaris)?. PLoS ONE 11, e0163620 (2016).PubMed 
PubMed Central 

Google Scholar 
27.Langford, F. M. & Cockram, M. S. Is sleep in animals affected by prior waking experiences?. Anim. Welf. 19, 215–222 (2010).CAS 

Google Scholar 
28.Toth, L. A. & Bhargava, P. Animal models of sleep disorders. Comp. Med. 63, 91–104 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
29.Zanghi, B. M., Kerr, W., DeRivera, C., Araujo, J. & Milgram, B. Sleep and biorhythms as a function of age in the dog. J. Vet. Behav. 5, 159 (2010).
Google Scholar 
30.Zanghi, B. M. Circadian Biorhythms of Sleep/Wake Dogs, Activity/Rest Cycles In Adult and Aged dogs. in Nestle Purina Companion Aninmal Nutrition Summit vol. 31 114–117 (Purina Institute, 2010).31.Zanghi, B. M. et al. Characterizing behavioral sleep using actigraphy in adult dogs of various ages fed once or twice daily. J. Vet. Behav. 8, 195–203 (2013).
Google Scholar 
32.Bódizs, R., Kis, A., Gácsi, M. & Topál, J. Sleep in the dog: Comparative, behavioral and translational relevance. Curr. Opin. Behav. Sci. 33, 25–33 (2020).
Google Scholar 
33.Kis, A. et al. Development of a non-invasive polysomnography technique for dogs (Canis familiaris). Physiol. Behav. 130, 149–156 (2014).CAS 
PubMed 

Google Scholar 
34.Takeuchi, T. & Harada, E. Age-related changes in sleep-wake rhythm in dog. Behav. Brain Res. 136, 193–199 (2002).PubMed 

Google Scholar 
35.Adams, G. J. & Johnson, K. G. Sleep-wake cycles and other night-time behaviours of the domestic dog Canis familiaris. Appl. Anim. Behav. Sci. 36, 233–248 (1993).
Google Scholar 
36.Adams, G. J. & Johnson, K. G. Behavioural responses to barking and other auditory stimuli during night-time sleeping and waking in the domestic dog (Canis familiaris). Appl. Anim. Behav. Sci. 39, 151–162 (1994).
Google Scholar 
37.Bunford, N. et al. Differences in pre-sleep activity and sleep location are associated with variability in daytime/nighttime sleep electrophysiology in the domestic dog. Sci. Rep. 8, 7109 (2018).PubMed 
PubMed Central 
ADS 

Google Scholar 
38.Iotchev, I. B. et al. Age-related differences and sexual dimorphism in canine sleep spindles. Sci. Rep. 9, 10092 (2019).PubMed 
PubMed Central 
ADS 

Google Scholar 
39.Kis, A. et al. Sleep macrostructure is modulated by positive and negative social experience in adult pet dogs. Proc. R. Soc. B Biol. Sci. 284, 20171883 (2017).
Google Scholar 
40.Mong, J. A. & Cusmano, D. M. Sex differences in sleep: Impact of biological sex and sex steroids. Philos. Trans. R. Soc. B Biol. Sci. 371, 20150110 (2016).
Google Scholar 
41.Andersen, M. L. et al. Effects of sleep loss on sleep architecture in Wistar rats: Gender-specific rebound sleep. Prog. Neuro-Psychopharmacol. Biol. Psychiatry 32, 975–983 (2008).CAS 

Google Scholar 
42.McKillop, L. E. et al. Effects of aging on cortical neural dynamics and local sleep homeostasis in mice. J. Neurosci. 38, 3911–3928 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
43.Nakamura, T. J., Takasu, N. N. & Nakamura, W. The suprachiasmatic nucleus: Age-related decline in biological rhythms. J. Physiol. Sci. 66, 367–374 (2016).CAS 
PubMed 

Google Scholar 
44.Abou-Ismail, U. A., Burman, O. H. P., Nicol, C. J. & Mendl, M. Can sleep behaviour be used as an indicator of stress in group-housed rats (Rattus norvegicus)?. Anim. Welf. 16, 185–188 (2007).CAS 

Google Scholar 
45.Sadeh, A., Keinan, G. & Daon, K. Effects of stress on sleep: The moderating role of coping style. Heal. Psychol. 23, 542–545 (2004).
Google Scholar 
46.Okun, M. L. Biological consequences of disturbed sleep: Important mediators of health?. Jpn. Psychol. Res. 53, 163–176 (2011).PubMed 
PubMed Central 

Google Scholar 
47.Van Reeth, O. et al. Interactions between stress and sleep: From basic research to clinical situations. Sleep Med. Rev. 4, 201–219 (2000).
Google Scholar 
48.Novati, A. et al. Chronically restricted sleep leads to depression-like changes in neurotransmitter receptor sensitivity and neuroendocrine stress reactivity in rats. Sleep 31, 1579–1585 (2008).PubMed 
PubMed Central 

Google Scholar 
49.Taylor, K. D. & Mills, D. S. The effect of the kennel environment on canine welfare: A critical review of experimental studies. Anim. Welf. 16, 435–447 (2007).CAS 

Google Scholar 
50.Sales, G., Hubrecht, R., Peyvandi, A., Milligan, S. & Shield, B. Noise in dog kennelling: Is barking a welfare problem for dogs?. Appl. Anim. Behav. Sci. 52, 321–329 (1997).
Google Scholar 
51.Hewison, L. F., Wright, H. F., Zulch, H. E. & Ellis, S. L. H. Short term consequences of preventing visitor access to kennels on noise and the behaviour and physiology of dogs housed in a rescue shelter. Physiol. Behav. 133, 1–7 (2014).CAS 
PubMed 

Google Scholar 
52.Schwartz, J. R. L. & Roth, T. Neurophysiology of sleep and wakefulness: Basic science and clinical implications. Curr. Neuropharmacol. 6, 367–378 (2008).CAS 
PubMed 
PubMed Central 

Google Scholar 
53.Reicher, V. et al. Repeated afternoon sleep recordings indicate first-night-effect-like adaptation process in family dogs. J. Sleep Res. 29, 1–10 (2020).
Google Scholar 
54.Piccione, G. et al. Comparison of daily distribution of rest/activity in companion cats and dogs. Biol. Rhythm Res. 45, 615–623 (2014).
Google Scholar 
55.Kredlow, M. A., Capozzoli, M. C., Hearon, B. A., Calkins, A. W. & Otto, M. W. The effects of physical activity on sleep: A meta-analytic review. J. Behav. Med. 38, 427–449 (2015).PubMed 

Google Scholar 
56.Barber, N. Play and energy regulation in mammals. Q. Rev. Biol. 66, 129–147 (1991).CAS 
PubMed 

Google Scholar 
57.Ahloy-Dallaire, J., Espinosa, J. & Mason, G. Play and optimal welfare: Does play indicate the presence of positive affective states?. Behav. Processes 156, 3–15 (2018).PubMed 

Google Scholar 
58.Penev, P. D. Sleep deprivation and energy metabolism: To sleep, perchance to eat?. Curr. Opin. Endocrinol. Diabetes. Obes. 14, 374–381 (2007).PubMed 

Google Scholar 
59.Mavanji, V., Billington, C. J., Kotz, C. & Teske, J. A. Sleep and obesity: A focus on animal models. Neurosci. Biobehav. Rev. 36, 1015–1029 (2012).PubMed 
PubMed Central 

Google Scholar 
60.Zanghi, B. M., Kerr, W., de Rivera, C., Araujo, J. A. & Milgram, N. W. Effect of age and feeding schedule on diurnal rest/activity rhythms in dogs. J. Vet. Behav. 7, 339–347 (2012).
Google Scholar 
61.Knutson, K. L., Spiegel, K., Penev, P. & Van Cauter, E. The metabolic consequences of sleep deprivation. Sleep Med. Rev. 11, 163–178 (2007).PubMed 
PubMed Central 

Google Scholar 
62.Kiddie, J. & Collins, L. Identifying environmental and management factors that may be associated with the quality of life of kennelled dogs (Canis familiaris). Appl. Anim. Behav. Sci. 167, 43–55 (2015).
Google Scholar 
63.Taylor, K. D. & Mills, D. S. The effect of the kennel environment on canine welfare: A critical review of experimental studies (2007).
64.CONCEA. Resolução normativa no 12. Diretriz brasileira para o cuidado e a utilização de animais para fins científicos e didáticos. Diário Oficial da União No 186 (2013).65.Bateson, M. & Martin, P. Measuring Behaviour (Cambridge University Press, 2021) https://doi.org/10.1017/CBO9780511810893.Book 

Google Scholar 
66.Broom, D. M. & Fraser, A. F. Domestic Animal Behaviour and Welfare (CABI, 2015).
Google Scholar 
67.Luescher, U. A., McKeown, D. B. & Halip, J. Stereotypic or obsessive-compulsive disorders in dogs and cats. Vet. Clin. North Am. Small Anim. Pract. 21, 401–413 (1991).CAS 
PubMed 

Google Scholar 
68.Friard, O. & Gamba, M. BORIS: A free, versatile open-source event-logging software for video/audio coding and live observations. Methods Ecol. Evol. 7, 1325–1330 (2016).
Google Scholar 
69.Crawley, M. J. The R Book (Wiley, 2007). https://doi.org/10.1002/9780470515075.Book 
MATH 

Google Scholar 
70.Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting Linear Mixed-Effects Models Using lme4. J. Stat. Softw. 67, 1–48 (2015).
Google Scholar 
71.Becker, R. A., Chambers, J. M. & Wilks, A. R. The New S Language (Chapman and Hall/CRC, 1988).MATH 

Google Scholar  More

Concentrations and sources of organic compoundsMid- and long-chain n-alkanesThroughout most sections, the biomarker composition of the Hala Hu record is dominated by long-chain n-alkanes (Fig. 2a). In most samples, concentrations were highest for nC31 (ca 5–35 µg/g d.w.) and decreasing with chain-length (Supplementary Fig. S3–1). Generally, the concentrations were low in the glacial period, and increased between ca. 10 and 8 ka cal BP, while gradually decreasing after ca. 5 ka cal BP. These alkanes can be attributed to vascular higher plants from the lake catchment27, which is characterized by alpine meadows. Mid-chain nC23 and nC25-alkanes are frequently attributed to aquatic macrophytes28. However, n-alkane patterns of aquatic and terrestrial plants overlap. Moreover, we consider the contribution from macrophyte to the sedimentary n-alkane pool at the coring location to be minor, because of the specific n-alkane pattern of the samples, the overall low concentrations of mid-chain n-alkanes, and the deep water depth of the core site. An exception is a section in the core dated to ca. 15–14 ka cal BP, that has concentrations from ca. 20 up to 50 µg/g d.w. for individual mid-chain n-alkanes (Fig. 2a). This can be explained by enhanced contribution from submerged aquatic plants29 when the lake level was much lower than present9,23. Alternatively, a contribution from other mid-chain producers, such as Sphagnum species30,31, is possible during phases of lake regression and the potential formation of peatlands around the lake shore.Fig. 2: Summary of measured proxy parameters in cores H7 and H11 (overlap 7.4−9.1 kyrs BP).a δD values and concentrations of mid- and long-chain n-alkanes. Arrows (ASM) indicate the maximum strength of the summer monsoon over Asia. b Concentrations of aquatic biomarkers (C20-HBI, alkenones, n-alkenes) and of microbial derived PMI; alkenone index (Uk´3738) and C37-alkenone δD values. c Concentrations of firemarkers (ΣM, L, G). d Titanium, sulfur, strontium, and calcium contents from XRF-scanning22,23. e Lake level reconstructed from ostracod assemblages9. Dark dashed interval 8.4−8.0 kyrs BP indicates mass flow layer. Light and medium grey shaded areas mark episodes of late glacial and mid-Holocene regime shifts within the aquatic ecosystem.Full size imageAlgal biomarkersAside from n-alkanes, a range of compounds of mostly aquatic origin can be identified in the aliphatic and ketone fraction of the sediment extracts. Unsaturated mid-chain alkenes nC21:1, nC23:1, nC25:1, and nC27:1 were abundant in traces in large parts of the core, but exhibit very high concentrations ( >100 µg/g d.w.) in the glacial period in a single sample at 17.5 ka cal BP and between 15 and 14 ka cal BP. Relatively high concentrations up to 25 µg/d.w. were also observed in the mid-Holocene sequence from ca. 9 to 5 ka cal BP (Fig. 2b). The origin of those compounds is not fully resolved and so far n-alkenes have not been detected in samples from aquatic and terrestrial vascular plants on the TP. However, algae such as eustigmatophytes and chlorophytes have been suggested as possible precursors in the open freshwater Lake Lugu, southeastern TP32, and in Lake Challa, eastern Africa33,34. Therefore, we consider phytoplankton as the most likely source of those compounds in Hala Hu.It is notable that different n-alkene distribution patterns are visible throughout the core, with a predominance of nC27:1 and nC25:1 in the mid-Holocene section and nC23:1 and nC25:1 in other core sections (Supplementary Fig. S3–1). This indicates either different source organisms for at least nC23:1 compared to nC27:1 (supported by more depleted δD-values, see below and Supplementary Fig. S3–3) or a change towards the synthesis of longer chain-lengths by the same source organisms, triggered by changing environmental conditions.The C20-highly branched isoprenoid (HBI) compound is another observed phytoplankton biomarker, which is widely absent in the glacial sequences but shows traces throughout the Holocene with peak abundances (up to 30 µg/g d.w) in the mid-Holocene (8−5 ka cal BP; Fig. 2b). As often been found in cyanobacterial and algal mats (e.g.,35,36,37), it has recently been assigned as a trophic indicator derived from diatoms in lake systems38.Other observed biomarkers of algal origin are alkenones, which are primarily produced by haptophytes, even though a variety of species are possible precursors39,40 The summed concentrations of C37-, C38-, and C39-compounds were low (300 ng/g d.w.; ca. 7.8–6.3 ka cal BP) (Fig. 2b). Here, the di- and tri-unsaturated homologues of all chain lengths appear to elute as peaks undisturbed by contamination, while the tetra-unsaturated alkenones show co-elution with a fatty acid ethyl ester (FAEE) at least for the C37-compound.Pentamethylicosane (PMI) is a compound that was detected in most samples in minor concentrations (400 ng/g d.w.) during an episode in the glacial period (ca 16.6−14 ka cal BP) and during the late Holocene (4.5 ka cal BP to present). This compound has been assigned to microorganisms, i.e bacteria and archaea, often related to the methane cycle41,42,43,44.Fire markersAnother analyzed biomarker group, anhydrous sugars levoglucosan (L), galactosan (G), and mannosan (M) are generated by combustion and pyrolysis of cellulose and hemicellulose, thus, are often referenced as “pyromarkers” or “firemarkers”45,46,47. They occur in low concentrations ( More

Study areaThe North Sea is a continental sea connected to the Atlantic Ocean through the English Channel in the southwest and between northern Shetland along the 61° latitude parallel to Norway in the north (Fig. 1). It is bordered by Norway, Denmark, Germany, the Netherlands, Belgium, France and the UK, and has a surface of 570,000 km2. The North Sea has an average depth of 95 m, yet maximum depths of ca. 700 m are found in the Norwegian Trench. The maximum tidal amplitude of the North Sea can reach up to 8 m, average winter sea surface temperatures are ca. 6 °C and average summer temperatures reach ca. 17°C33. The English Channel encompasses the marine strait between the UK and France. It covers 75,000 km2, has an average depth of 63 m, a maximum depth of 174 m and can reach a maximum tidal amplitude up to 12 m. The average winter and summer sea surface temperatures in the English Channel are ca. 5 and 20 °C, respectively54.TaggingIn total, 320 silver eels were tagged with pop-off archival tags (Table 1; Supplementary Table S2). In Belgium, 238 eels were caught and tagged at a drainage system upstream of the Yser Estuary (hereafter referred to as the Belgian eels) in 2018–2020 via nets that were attached to gravitational discharge sluice gates (coordinates: 51.127 N, 2.761 E) in October, November and December (n2018 = 102, n2019 = 60 and n2020 = 76). In Germany, 82 eels were tagged in 2011 and 2012. In early December 2011, seven eels were caught at Lake Plön (coordinates: 54.137 N, 10.334 E) with fyke nets. During September, October and November 2012, eels were caught in the Rivers Eider (n = 30; coordinates: 54.190 N, 9.093 E) and Havel (n = 45; coordinates: 52.419 N, 12.571 E) with fyke and stow nets, respectively.Upon capture, the eels were anaesthetized with 0.3 ml/L clove oil (Belgium), 0.4 ml/L ethylene glycol monophenyl ether (Germany 2011) or 120 mg/L MS-222 (Germany 2012), and various morphometric characteristics were measured to identify the life stage55: total length (to the nearest mm), weight (to the nearest g), horizontal and vertical eye diameter (to the nearest 0.01 mm in Belgium and to the nearest 0.1 mm in Germany) and pectoral fin length (to the nearest 0.01 mm and 0.1 mm in Belgium and Germany, respectively). Given that their total body length was  > 450 mm, all eels were considered female55. According to the morphometrics, five Belgian eels could be considered in the premigratory stage (FIII); however, based on visual inspection, they were considered silver eels (i.e. silver-coloured abdomen, dark grey on the dorsal side, jaw hinge not proceeding beyond the eye, enlarged eyes and dark coloured pectoral fins). The other 315 eels identified as silver eels based on both morphometry and visual inspection (201 FIV stage and 114 FV stage).Eels weighing ≥ 550 g were externally fitted with a G5 PDST (CEFAS Technology Ltd, UK), which log temperature and pressure (providing information on depth). They were attached applying the three-point Westerberg attachment method56. Two tag types were used: one with a separate tag and pop-off mechanism (Germany) and one where both mechanisms were integrated (Belgium). The flotation collar of the PDSTs was painted bright red, contained contact information and a cash reward to stimulate retrieval by the general public (e.g. beach combers and fishermen). The seven eels caught in 2011 in Germany (minimum 1220 g) were fitted with PSATs (X-Tag, Microwave Telemetry Inc., USA), also using the Westerberg-method56. Like the PDSTs, the PSATs record temperature and pressure. After release, they drift to the surface and transmit the data to the user via the ARGOS satellite system (www.argos-system.org). For the specifications of the different tags, we refer to Supplementary Table S3.Upon recovery from the anaesthetic, eels tagged with PDSTs were released close to their capture locations in the rivers Eider (coordinates 2011: 54.381 N, 9.009 E; coordinates 2012: 54.379 N, 9.013 E), Elbe (coordinates 1: 53.793 N, 9.402 E; coordinates 2: 53.569 N, 9.700 E; coordinates 3: 53.396 N, 10.171 E) and Yser (coordinates: 51.135 N, 2.757 E) (Table 1). The seven eels captured for PSAT tagging in 2011 were held for several weeks in the Thünen Institute of Fisheries Ecology, then tagged and released the same day; others were tagged in the field.PreprocessingOnce downloaded, the temperature and pressure data obtained from the PDSTs was subsampled to 1-min (Belgian eels) or 2-min (German eels) intervals to reduce the datasets and improve geolocation calculation time; this discrepancy is due to the minimum logging rate of the tags (Supplementary Table S3). Linear regression was applied to correct for pressure sensor drift over time. Indeed, pressure values increased over time even if the tag was kept at atmospheric pressure level. The regression was applied between 15 min before release and the moment the tag popped off and reached the surface, since the tag was then considered at sea level and hence to be under zero pressure.The PSAT data were retrieved through the ARGOS satellite system as a subset with 15-min intervals and converted to values of pressure and temperature. Contemporaneous values of temperature and depth were not always transmitted due to the transmission method. As a consequence of the tag release programming, the transmission of the first position for one of the tags was only received five days after the tag reached the sea surface.GeolocationThe daily movements of each electronically tagged European eel were reconstructed using an adapted version of the tidal geolocation model of Pedersen et al.57. The geolocation model uses a novel Fokker–Planck based method to combine the tidal location method of Metcalfe and Arnold58 with a hidden Markov model (HMM), such that an individual’s daily location d is modelled conditionally on its previous location (d − 1), its inferred behavioural state ds, where behaviour is defined by a single diffusivity parameter (i.e. the maximum amount of movement permitted in a given day), and the observations made between d and d − 1. In this case, observations consisted of the recorded depth (m; D1, …, Dn) and temperature (°C; T1, …, Tn), where n is the number of measurements made per day (the HMM down-samples to 10-min intervals, hence 144 measurements per day), and any hydrostatic (tidal) data which are derived from the sinusoidal pressure cycle recorded in the depth data when a fish is at rest on the seafloor. In addition to bathymetry and tidal amplitude with phase, the model was developed to include sea surface temperature (SST), which can provide additional validation when fish are swimming at or near the surface (i.e. depth ≤ 20 m)59,60, and temperature at depth, which can provide additional validation when fish remain at depths well below the sea surface61,62.The model was run in three different configurations for each recovered dataset: (i) using the tidal location model only (as for Pedersen et al.57), hereafter termed TLM geolocation; (ii) using the TLM plus sea surface temperature (as for Wright et al.60), hereafter termed SST geolocation, and (iii) using temperature at the surface and sub-surface, hereafter termed 3D geolocation (Supplemental Fig. S3). The final trajectory output for the PDST Belgian eels and PSAT German eels was obtained via 3D geolocation, while SST geolocation was used for the PDST German eels. The reason for this discrepancy is that the German PDST eels stayed closer to the coast and in shallower water. Consequently, the 3D geolocation results were more prone to error due to coastal influences on water temperature. As a result, we used the SST geolocation method for these datasets to obtain more reliable results.Data for the model were derived from publicly available resources. Gridded global bathymetry data were obtained from the general bathymetric chart of the oceans (Gebco; British Oceanographic Data Centre, Liverpool, United Kingdom, 2009). Tidal constituents were obtained from the Oregon State University Tidal Prediction model, as described in Egbert and Erofeeva63. Sea surface temperature data were sourced from OSTIA64, while temperature at depth data were sourced from the operational Mercator global ocean analysis and forecast system65. These datasets were downloaded from the Copernicus Marine Environmental Model Service (CMEMS: documented here http://resources.marine.copernicus.eu/documents/PUM/CMEMS-GLO-PUM-001-024.pdf). Data were sourced so as to fit the spatial scale of the model (30°N to 80°N and from 110°W to 60°E) and coarsened to reduce model run-time by modifying the spatial grid to a 1/10th of a degree resolution. The output of the model is a nonparametric probability distribution of the geographical position from which a most probable location, for each day at liberty, and a most probable movement path can be estimated.Prior to running the model, a number of constraints and input parameters were defined to ensure that the model ran effectively. The recapture information was either set as (a) the latitude and longitude where the tag was recaptured, with a high confidence ( 200 m) and hence did not exhibit diel vertical migrations, the input estimates of longitude were based on a simple linear interpolation from release to estimated pop-up. However, for eels that did reach oceanic depths, the time of local noon was estimated (based on the timing of significant diel vertical migrations, as for Righton et al.16), and used to estimate longitude. Geolocation was conducted with MATLAB software66.Migration routesOnly datasets containing ≥ 100 km of net tracking distance were included for further analysis, leading to 54 datasets from the 96 retrieved tags and 320 tagged eels. The net tracking distance was identified as the distance along the reconstructed trajectory between the release of the tagged eel and the pop-off event. When an eel was ingested by a predator, leading to the tag tracking the predator rather than the eel, the data were excluded from the day the eel was predated. The 100 km cut-off point was arbitrarily chosen to select migration paths of sufficient length for further analysis (e.g. migration direction); tracks had a minimum deployment duration of 4 days.Migration speedTo exclude a size-effect, we first applied an independent two-sample t-test to confirm eel sizes (i.e. weight) did not differ between Belgian and German eels. The assumptions of normality (Shapiro–Wilk test), homogeneity of variances (F-test) and independence were met (weight measurements are individual-specific and therefore independent).Next, an independent two-sample t-test was conducted to test if the total migration speeds (i.e. the ground speed along the reconstructed trajectory between the release of the tagged eel and the pop-off or predation event) differed between Belgian and German eels. The assumptions were tested and met as described above.Finally, we tested if the daily migration speed (i.e. the ground speed along the reconstructed trajectory per day) differed according to the eel’s position (i.e. modelled latitude and longitude) via a linear mixed effects model. The tag IDs were implemented as a random effect to account for autocorrelation. Since the two-sample t-test showed a significant difference between Belgian and German total migration speeds, we performed a separate analysis on eels from both countries. Assumptions of normality, homogeneity of variances and independence were tested and met.The migration speed analyses were conducted in R (version 3.6.3)67. The packages ‘lme4’ and ‘nlme’ were used to conduct the linear mixed effects model.Ethical statementEels were tagged using approved protocols by trained and individually licensed scientists working under national project authority in accordance with institutional and national guides for the care and use of laboratory animals. These guidelines are consistent with Institutional Review Board/Institutional Animal Care and Use Committee guidelines. Tagging in Belgium was carried out in accordance with the Belgian national and regional regulations for animal welfare and treatment (Permit ID: EC INBO-011). Tagging in Germany followed German legislation concerning care and use of laboratory animals, and ethical permission for the experiments was given by the Ministry of Energy, Agriculture, the Environment, and Rural Areas of the federal state Schleswig–Holstein (reference numbers V312-72241.123-34 (90-8/11) and V311-7224.123.3 (93-6/12) for tagging in 2011 and 2012 respectively). More

Plant defence traitsWe compiled species level data for five plant traits: wood density (WD), leaf and stem spinescence, latex production, and leaf size, for tropical and extra-tropical South and Central American woody species (i.e., the Neotropical biogeographic realm). WD was obtained for 2577 species from ref. 44. We only used wood density data from Zanne et al.44, because this study used WD measured in stems, whereas most other studies with available data used WD measured in branches. Leaf size data were obtained for 2660 woody species from Wright et al.37. We did not include leaf size from herbaceous species because herbaceous and woody species are influence by different megafauna guilds, suggesting distinct mechanisms, and because this dataset37 only included data for 253 Neotropical herbaceous species. The presence or absence of stem (and/or branch) spines (mostly thorns, but also prickles) were obtained from Dantas and Pausas45 for Neotropical savanna and forest species (1004 species) and complemented with other literature sources for other ecoregions (listed in the supplementary materials) using the names of the species for which we had WD and Leaf Size data. Our final stem spines dataset included 2843 woody species. We also compiled data on the presence of latex in plant stems and leaves for all the species for which we had data on other traits (3160 species; references in the supplementary materials). Finally, we also compiled data on leaf spines. While we managed to find leaf spine data for a total of 2173 woody species, we found spinescence in leaves to be especially concentrated in the palm Family (Arecaceae; 198 out of 221 species with leaf spines). Moreover, out of the non-palm species, all but three species also presented stem spines, indicating that, for other taxa, leaf spines might be dependent on the presence of stem spines at the region (in palms, 51% have stem spines). Thus, we only used leaf spinescence data of palm species (694 species) from the global Palm Traits Database 1.046.For wood density and leaf size, we often had more than one trait value per species (1005 and 831 species with more than one trait value, respectively). Thus, we computed the species mean trait value. This rarely occurred for binary traits (spinescence and latex) and, when occurred, the maximum value was used (0 for absence and 1 for presence). This later decision was based in the assumption that omitting the presence of spines or latex is more likely than incorrectly reporting the presence when it is absent. Moreover, some of these traits can be plastic18.From species to ecoregionsWe searched for geographical distribution data (coordinates) from the Global Biodiversity Information Facility (GBIF) for all of the species in each species-trait dataset (Data available from GBIF using the following doi: WD: https://doi.org/10.15468/dl.3vua3x; Stem spines: https://doi.org/10.15468/dl.ar5ddj; Latex: https://doi.org/10.15468/dl.m8dzjd; Leaf spines: https://doi.org/10.15468/dl.vv8gw4; Leaf size: https://doi.org/10.15468/dl.k98nxc). For this search, we used tools provided by the “rgbif” package for R in which species names are updated to the most recent classification and the returned occurrences also include those associated with synonyms (i.e., the “backbone” method). We labelled the obtained geographical coordinates according to their ecoregion and biogeographical realm (following Dinerstein et al.47) and cropped out occurrences falling outside of the Neotropical realm. Since occurrence data was not available to all the species in our initial trait dataset, the number of species used to calculate ecoregion level means was reduced to 2110 species, for wood density, 2133, for leaf size, 2629, for stem spines, 2714, for latex, and 657, for leaf spines. A detailed evaluation of the representativity of this data in relation to ecoregion- and Neotropical- level patterns can be found in the Supplementary Methods. Based on the occurrence data and their ecoregion label, we built a species abundance (columns) by ecoregion (rows) matrix for each trait.We obtained ecoregion scale abundance-weighted means for continuous traits (WD and Leaf Size) by: (1) Multiplying species abundance in each grid cell of the ecoregion by the mean species value; (2) Summing up the row values; (3) dividing the resulting row sum by the total species abundance (row sum prior to trait multiplication), and (4) calculating the ecoregions’ means (across all of the grid cells). For Stem Spines and Latex (binary traits), we used a similar procedure, but the maximum (0 for absence and 1 for presence) value was used instead of the mean in step (1), and step (2) was directly used to calculate the number of presences (i.e., 1 s). Moreover, instead of the steps (3) and (4), we calculated the number of absences as the difference between the total abundance (row sums before trait multiplication) and the values obtained in step (2). This process resulted in weighted means for WD and stem spinescence for 173 ecoregions, and Leaf Size and Latex for 174, out of the 179 Neotropical ecoregions. For leaf spinescence, we used a similar approach, although, because of the fewer species, the abundance estimate from GBIF was less reliable. Thus, we transformed the ecoregion species abundance to presence/absence before multiplying the trait values (0/1 for absence/presence). We obtained leaf spinescence data for 159 out of the 179 Neotropical ecoregions. The species- and ecoregion- level data is provided in the Supplementary Data and in ref. 47.Historical megafauna distributionWe obtained data on historical distribution of megafauna species from the MegaPast2Future/PHYLACINE_1.2 dataset24, a dataset containing distribution maps (96.5 km of spatial resolution) and functional traits for mammal species of the last 130,000 years. From this dataset, we obtained the probable past distribution of extinct large mammal herbivore (hereafter, “megafauna”) species, if these species were still alive today (“Present Natural” scenario; see details below). The “Present Natural” distribution of extinct species in this database is based on the estimated historical distribution (i.e., preceding anthropogenic range modifications) of extant species that are known (from the fossil record) to have coexisted with the extinct species. In this approach, an extinct species is considered to have been present in a given grid cell if at least 50% of the extant species that were found coexisting with the extinct species in the fossil (and subfossil) record was predicted to have occurred in the same cell prior to anthropogenic range modifications24,48. This approach assumes that, since extant and extinct species coexisted in the same locations, they must have had similar ecological requirements. It also assumes that megafauna extinction had anthropogenic causes, instead of causes related to climate change49, which is largely accepted in the literature50.We extracted the “Present Natural” distribution of extinct mammal (coded “EP” for IUCN status; i.e., “extinct in prehistory”, meaning before 1500 CE) whose body mass was higher than 50 kg (megafauna), and for which at least 90% of their diet consisted of plants (i.e., strict herbivores). For each Ecoregion, we began by calculating two megafauna-related metrics: extinct megafauna species richness (Mrich) and their mean body mass (Mbm). For this, we cropped the distribution maps of the megafauna species (containing 1 for presence and 0 for absence of each species) to the Neotropical realm. To calculate Mrich, we (1) counted species presences within each of the grid cells in the global grid (i.e., calculated the cell’s megafauna richness); (2) assigned the corresponding ecoregion label to the resulting richness grid cells, subset the richness cell values corresponding to the Neotropical region; and (3) calculated the mean for each Neotropical ecoregion. For Mbm, we replaced the presences of the megafauna species in the initial raster object (grid cell map of each megafauna species) by their body masses and calculated the grid cell-level mean body mass, before calculating the ecoregion-level means. We also calculated megafauna density and secondary productivity based on allometric equations that relate these metrics to megafauna body mass. However, we did not used megafauna density and secondary productivity because they were strongly correlated to megafauna richness (Supplementary Fig. 3). More details on how these metrics were calculated can be found in the Supplementary Methods.We also obtained diet preference information from the literature for most megafauna species that occurred in the Neotropical region (details and references in the Supplementary Material). Based on these information, we calculated the richness of large browser (MBrich for megabrowser richness), grazer (MGrich for megagrazer richness), and mixed-feeder (MMfrich for mega mixed-feeder richness) species by sub setting the megafauna species by grid cell array before the richness calculation in order to select only species that were classified within the correspondent subgroup.Extant herbivore mammal distributionWe also compiled data on the distribution, body mass and diet of extant and recently extinct (i.e., extinct after 1500 CE) herbivore mammal species (for simplicity, called ‘extant’ species in this study). As with megafauna maps, the distributions used represented reconstructions for periods preceding anthropogenic reduction of extant herbivores ranges (“Present Natural” scenario), based on abiotic, biotic and geographic variables48, rather than the currently observed distribution. This scenario was used because modern anthropogenic range reductions are too recent to produce substantial geographic effects at this spatial scale. These data were obtained by sub setting the MegaPast2Future/PHYLACINE_1.2 dataset to exclude species that were coded “EP” for IUCN status and that were not strict herbivores (at least 90% of the diet constituting of plants). We subsequently associated diet information to these species using data from ref. 51 and excluded all species that did not feed mainly on aboveground vegetative plant tissues (i.e., species that fed mostly on fruits, seed, roots were excluded). This later filtering was because the number of herbivores that feed mostly on seed and fruit increase with decreasing size (and this dataset included small mammals). We subsequently calculated the same metrics as for the extinct megafauna species (except for the richness of mixed-feeders as our source for diets50 labelled species according to dominant feeding pattern). For this, we used the same approach described for extinct megafauna species. We did not use a size threshold for extant species because there were only 13 extant mammal herbivore species with over 50 kg in the Neotropical region, most of which were grazers (9 species; 4 species were mixed-feeders and none were browsers). Therefore, we relied on the mean body mass metric calculated for extant mammals to detect potential size-related effects.Climate, soil, fire, insularity, and hurricanesFor each Ecoregion, we obtained data on climate (mean annual precipitation and temperature, and rainfall seasonality) and soil (sand content, pH, and cation exchange capacity) variables. Climate data was obtained from WorldClim 2.1 (10 min spatial resolution) and was based on climate data from 1970 to 200052. Soil data were obtained from SoilGrids (5 km of spatial resolution)53, and consisted of mean values for two depths, 0.05 and 2 m. We calculated Ecoregion level means for all of the soil and climate variables after intersecting the climate and soil grid maps with the ecoregion map.We obtained the number (a proxy for frequency) and intensity of wildfires per ecoregion area using the MODIS active fire location product (MCD14ML)54. We only considered fires (i.e., hotspots) with detection confidence of 95% or higher occurring from November 2000 to December 2019 (both included). To ensure that only wildfires were considered, we associated each fire pixel with a land cover type (300 m of spatial resolution) from ref. 55 for a buffer area of 1000 m surrounding the fire pixel centroid. We excluded all of the fires occurring in areas in which more than 10% of the surrounding land cover pixels corresponded to agricultural, urban and water classes. We calculated the number of wildfires per ecoregion area by dividing the fire count of each Ecoregion by the ecoregion area, and multiplying the resulting value by the proportion of vegetated land cover pixels (same classes used to exclude fires in anthropogenic areas and water bodies above). Fire intensity was estimated as the average fire radiative power across all detected MODIS hotspots in the ecoregion. Ecoregions lacking large preserved vegetated areas (criteria above) were excluded from subsequent analyses.Using the ecoregion map, we also classified ecoregions into insular (1), when most of the ecoregion area was located in islands, vs. continental (0), otherwise. This was performed because island biogeography theory predicts that, in island, species richness should be low due to low colonization and high extinction rates. Insularity has also been shown to reduce megafauna body size (i.e., the island rule), even though the mechanisms are not fully understood56. We also compiled data on hurricane activity, as woody density was suggested to confer resistance against this disturbance57. We used data from 1990 to 2019 from the HURDAT2 dataset58, containing six-hourly information about the location of all of the known tropical and subtropical cyclones (0.1° latitude/longitude). We used the sum of hurricane occurrences per ecoregions divided by ecoregion area as an indicator of hurricane activity.Statistical analysesTo understand megafauna patterns, we began by fitting (multiple) regression models with habitat-related (fire, climate, soil) and geographical (insularity) variables as predictors. We expected that megafauna richness in general was higher under savanna conditions (arid nutrient-rich or mesic nutrient-poor environments with frequent fires)1,22. We also expected that megafauna richness and body mass were affected negatively by insularity (i.e., following the island biogeography theory and island rule). Before the analyses, we tested the correlations among all of the variables that would eventually be entered as predictors in the same model for both the megafauna and trait models (Supplementary Table 1), in order to avoid multicollinearity associated with highly correlated variables (here, r ≥ 0.60). Since mean annual precipitation and soil pH were strongly positively correlated (r = −0.78), for all of the analyses (including the analyses with functional traits, described below), model selection was performed separately for these two variables (i.e., two different model selection procedures, one containing each of the two variables among the initial set of predictors). We selected the best among the two resulting models as that with the lowest AIC (differences higher than two points in all of the cases). To make sure that no multicollinearity remained we also calculated the Variation Inflation Factor (VIF) for all of the predictor variables as 1/tolerance, where tolerance is calculated as 1 minus the R2 of all of the model regressing a predictive variable against all of the other predictors. In all of the models, VIF was 3.33 or smaller (i.e., a tolerance of 0.30 or higher), indicating absence of multicollinearity.Model simplification was carried interactively using stepwise (both forward and backward) searching for the model with the lowest AIC (using R’s “step” function) and subsequently retaining only the significant variables (p ≤ 0.05). We calculated the Pearson r statistics as a measure of effect size for the selected variables as well as the associated confidence intervals, using the packages “parameters” and “effectsize” for R. The average contribution of each predictor variable was also calculated, using the package “dominanceanalysis”, as the mean difference in R2 before and after removing the target variable from models containing all of the possible subset combinations of the selected predictor variables, including the full selected model.For testing whether the studied plant functional traits were related to our megafauna indicators, we fit linear models to WD and leaf size, and generalized linear models (GLM; binomial family) for spinescence and latescence, using ecoregion as the unit. For spinescence and latescence, we used the matrix containing the count of spiny/latex and non-spiny/non-latex plants (species abundance; for stem spines and latex) or number of species with or without spines (for leaf spines; see above) as response variables. The predictor variables included the animal indicators for extinct megafauna and extant herbivores, as well as climate, soil, and fire predictors (and, for WD, hurricane counts). Because total, as well as megagrazer, megabrowser, and mega mixed-feeder species richness were strongly positively correlated (Supplementary Table 1), we used the richness difference between grazers and browsers to evaluate the effect of diet (Supplementary Fig. 1). For consistency, we used the same diet variable for extant and extinct species. Since we did not identify strong correlations among extinct megafauna and extant herbivore indicators (Supplementary Table 1), these variables were all entered simultaneously in the same initial models. As with the analyses of the megafauna indices, we also used r as effect size and calculated the average predictor contribution in terms of R2 for these models. For the later, we used the MacFadden Pseudo-R2 in the GLM models as implemented in the “pscl” and “dominanceanalysis” packages for R, as this statistic is the most comparable with R2 from linear multiple regression (Maximum Likelihood and Cragg and Uhler’s Pseudo-R2 were also calculated for the logistic models), and adjusted R2 for continuous traits. Islands were not included in these models, as island plants were expected to respond differently due to the effects of insularity on animal species richness, precluding megafauna and extant mammal richness from being accurate proxies for consumer abundance. For stem spines, we always included a quadratic term to both megafauna and extant mammal herbivore body mass, as evidence suggest that medium-size herbivores (i.e., approximately 250 kg) are important selective drivers of this trait12. If a significant relationship with our herbivory indicators (both extant and extinct) were significant but not indicative of a selective effect by herbivores (for more defended plants), this relationship was discarded (along with related variables, such as diet); this happened only once, for leaf size, which increased with extant herbivore richness (Supplementary Table 8).For all of the general linear regression models, assumptions of normality, homoscesticity and lack of spatial autocorrelation in the residuals were checked using the Kolmogorov–Smirnov, Breusch–Pagan and Moran’s I tests, respectively. For the later, ecoregions were considered neighbours when they were adjacent and non-neighbour otherwise. In some cases, heteroscesticity was detected and, thus, the significance of the coefficients was tested using heteroskedasticity-consistent covariance matrix estimation. If one or more variable lost their significance they were stepwise removed from the final model, beginning by the least significant, until all remaining variables had a significant effect. Overdispersion in the generalized linear model was also detected and dealt with using overdispersed binomial logit models, as implemented in the “dispmod” package for R, in which weights are interactively calculated and used to maintain the residual deviance lower than the degrees of freedom. To confirm that the detected associations between megafauna indices and plant traits were robust, we also tested the coefficient significance using randomization of the plant species by ecoregion matrices (see Supplementary Methods for details).To test the prediction that Neotropical ecoregions could be broadly classified into the three hypothesised antiherbiomes, we used hierarchical clustering on principal component axes of the ecoregion by trait matrix (five plant traits, standardized to zero mean and unit variance). We selected the number of clusters associated with the highest loss of inertia (within group variability) when progressively increasing the number of clusters, using the R package “FactoMineR”. This procedure allowed the recognition of large regions characterised by specific patterns of defence strategies (‘antiherbiomes’). We subsequently tested for axes score, megafauna and environmental differences among the resulting antiherbiomes to verify whether and how trait, climate and soil patterns matched those described for African ecosystems, and to understand the megafaunal differences among the antiherbiomes. For these comparisons, we used Kruskal-Wallis and post-hoc pairwise Dunn tests, using the Benjamini & Hochberg59 (1995) correction of P-values for multiple comparisons in both cases, and exclusively included continental ecoregions. For spines, we used the proportion of spinescent plants/species (rather than the number of “yes” and “no” used on previous analyses) in the principal component analysis. Because palms were missing from 20 ecoregions, we completed the values for these ecoregions using predicted model probabilities. To better understand these associations between traits and the environmental and megafauna variables, we also regressed the PCA axes against the same predictors used for traits.We also developed a framework to identify forest ecoregions most likely to have experienced a biome shift after megafauna extinction using antiherbiome, biome and megafauna distribution data. Ecoregions likely to have experienced a savanna-to-forest shift since the Pleistocene are those that: (1) are currently forest-dominated; (2) are classified in antiherbiomes analogous to African arid nutrient-rich or mesic nutrient-poor savannas; and (3) were megafauna- and, especially, megagrazer- rich during the Pleistocene (richness equal or greater than the 0.75 quantile: 14 species for Mrich, and 3 for exclusively grazing species; MGrich). We validated the distribution of these areas with fossil evidence (22 sites) from the Last Glacial Maximum and mid-Holocene (see Supplementary Methods and Supplementary Table 9). For this, we also used information about the present dominant vegetation type in the fossil sites, extracted from the reference sources (see Supplementary Table 9), to segregate savanna-forest shifts from data coming from stable savanna patches within forest or long-term savanna regions. We also contrasted the predicted patterns with the present location of savanna patches within the Amazon Forest region from ref. 60.All statistical analyses and data handling were carried out in the R (v.4.0.2) environment, using the previously mentioned packages, in addition to FSA, gridExtra, grid, lattice, lmtest, latticeExtra, olsrr, raster, rgdal, rgeos, sandwich, spatialreg, spdep and vegan, using codes provided in ref. 47.Reporting SummaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More