Philippa Kaur

1.Mineau, P. & Palmer, C. Neonicotinoid Insecticides and Birds: The Impact of the Nation’s Most Widely Used Insecticides on Birds. (American Bird Conservancy, USA, 2013).2.Simon-Delso, N. et al. Systemic insecticides (Neonicotinoids and fipronil): Trends, uses, mode of action and metabolites. Environ. Sci. Pollut. Res. 22, 5–34. https://doi.org/10.1007/s11356-014-3470-y (2015).CAS 
Article 

Google Scholar 
3.Jeschke, P., Nauen, R., Schindler, M. & Elbert, A. Overview of the status and global strategy for neonicotinoids. J. Agric. Food Chem. 59, 2897–2908. https://doi.org/10.1021/jf101303g (2011).CAS 
Article 
PubMed 

Google Scholar 
4.Tomizawa, M. & Casida, J. E. Selective toxicity of neonicotinoids attributable to specificity of insect and mammalian nicotining receptors. Annu. Rev. Entomol. 48, 339–364. https://doi.org/10.1146/annurev.ento.48.091801.112731 (2003).CAS 
Article 
PubMed 

Google Scholar 
5.Woodcock, B. A. et al. Impacts of neonicotinoid use on long-term population changes in wild bees in England. Nat. Commun. 7, 1–8. https://doi.org/10.1038/ncomms12459 (2016).CAS 
Article 

Google Scholar 
6.Pisa, L. et al. An update of the Worldwide Integrated Assessment (WIA) on systemic insecticides. Part 2: Impacts on organisms and ecosystems. Environ. Sci. Pollut. Res. 28, 1–49. https://doi.org/10.1007/s11356-017-0341-3 (2017).CAS 
Article 

Google Scholar 
7.Li, Y., Miao, R. & Khanna, M. Neonicotinoids and decline in bird biodiversity in the United States. Nat. Sustain. 3, 1027–1035. https://doi.org/10.1038/s41893-020-0582-x (2020).Article 

Google Scholar 
8.Eng, M. L., Stutchbury, B. J. & Morrissey, C. A. Imidacloprid and chlorpyrifos insecticides impair migratory ability in a seed-eating songbird. Sci. Rep. 7, 1. https://doi.org/10.1038/s41598-017-15446-x (2017).CAS 
Article 

Google Scholar 
9.Eng, M. L., Stutchbury, B. J. & Morrissey, C. A. A neonicotinoid insecticide reduces fueling and delays migration in songbirds. Science 80(365), 1177–1180. https://doi.org/10.1126/science.aaw9419 (2019).ADS 
CAS 
Article 

Google Scholar 
10.Lopez-Antia, A., Ortiz-Santaliestra, M. E., Mougeot, F. & Mateo, R. Imidacloprid-treated seed ingestion has lethal effect on adult partridges and reduces both breeding investment and offspring immunity. Environ. Res. 136, 97–107. https://doi.org/10.1016/j.envres.2014.10.023 (2015).CAS 
Article 
PubMed 

Google Scholar 
11.Pandey, S. P. & Mohanty, B. The neonicotinoid pesticide imidacloprid and the dithiocarbamate fungicide mancozeb disrupt the pituitary-thyroid axis of a wildlife bird. Chemosphere 122, 227–234. https://doi.org/10.1016/j.chemosphere.2014.11.061 (2015).ADS 
CAS 
Article 
PubMed 

Google Scholar 
12.Tokumoto, J. et al. Effects of exposure to clothianidin on the reproductive system of male quails. J. Vet. Med. Sci. 75, 755–760. https://doi.org/10.1292/jvms.12-0544 (2013).CAS 
Article 
PubMed 

Google Scholar 
13.Addy-Orduna, L. M., Brodeur, J. C. & Mateo, R. Oral acute toxicity of imidacloprid, thiamethoxam and clothianidin in eared doves: A contribution for the risk assessment of neonicotinoids in birds. Sci. Total Environ. 650, 1216–1223. https://doi.org/10.1016/j.scitotenv.2018.09.112 (2019).ADS 
CAS 
Article 
PubMed 

Google Scholar 
14.Berheim, E. H. et al. Effects of Neonicotinoid Insecticides on Physiology and Reproductive Characteristics of Captive Female and Fawn White-tailed Deer. Sci. Rep. 9, 1–10. https://doi.org/10.1038/s41598-019-40994-9 (2019).CAS 
Article 

Google Scholar 
15.Wang, Y. et al. Unraveling the toxic effects of neonicotinoid insecticides on the thyroid endocrine system of lizards. Environ. Pollut. 258, 113731. https://doi.org/10.1016/j.envpol.2019.113731 (2020).CAS 
Article 
PubMed 

Google Scholar 
16.Khalil, S. R., Awad, A., Mohammed, H. H. & Nassan, M. A. Imidacloprid insecticide exposure induces stress and disrupts glucose homeostasis in male rats. Environ. Toxicol. Pharmacol. 55, 165–174. https://doi.org/10.1016/j.etap.2017.08.017 (2017).CAS 
Article 
PubMed 

Google Scholar 
17.Abou-Donia, M. B. et al. Imidacloprid induces neurobehavioral deficits and increases expression of glial fibrillary acidic protein in the motor cortex and hippocampus in offspring rats following in utero exposure. J. Toxicol. Environ. Heal. – Part A Curr. Issues 71, 119–130. https://doi.org/10.1080/15287390701613140 (2008).CAS 
Article 

Google Scholar 
18.Gawade, L., Dadarkar, S. S., Husain, R. & Gatne, M. A detailed study of developmental immunotoxicity of imidacloprid in Wistar rats. Food Chem. Toxicol. 51, 61–70. https://doi.org/10.1016/j.fct.2012.09.009 (2013).CAS 
Article 
PubMed 

Google Scholar 
19.Mohanty, B., Pandey, S. P. & Tsutsui, K. Thyroid disrupting pesticides impair the hypothalamic-pituitary-testicular axis of a wildlife bird. Amandava amandava. Reprod. Toxicol. 71, 32–41. https://doi.org/10.1016/j.reprotox.2017.04.006 (2017).CAS 
Article 
PubMed 

Google Scholar 
20.Sun, Q. et al. Imidacloprid Promotes High Fat Diet-Induced Adiposity in Female C57BL/6J Mice and Enhances Adipogenesis in 3T3-L1 Adipocytes via the AMPK(alpha)-Mediated Pathway. J. Agric. Food Chem. 65, 6572–6581. https://doi.org/10.1021/acs.jafc.7b02584 (2017).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
21.Sun, Q. et al. Imidacloprid promotes high fat diet-induced adiposity and insulin resistance in male C57BL/6J mice. J. Agric. Food Chem. 64, 9293–9306. https://doi.org/10.1021/acs.jafc.6b04322 (2016).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
22.Park, Y. et al. Imidacloprid, a neonicotinoid insecticide, potentiates adipogenesis in 3T3-L1 adipocytes. J. Agric. Food Chem. 61, 255–259. https://doi.org/10.1021/jf3039814 (2013).CAS 
Article 
PubMed 

Google Scholar 
23.Ricklefs, R. E., Stark, J. M. & Konarzewski, M. Internal constraints on growth in birds. in Avian Growth and Development. Evolution within the Altricial-Precocial Spectrum (eds Starck, J. M. & Ricklefs, R.E.) 266–287 (Oxford Ornithology Series, Oxford, 1998).
Google Scholar 
24.Bobek, S., Jastrzebski, M. & Pietras, M. Age-related changes in oxygen consumption and plasma thyroid hormone concentration in the young chicken. Gen. Comput. Endocrinol. 31, 169–174. https://doi.org/10.1016/0016-6480(77)90014-4 (1977).CAS 
Article 

Google Scholar 
25.Metcalfe, N. B. & Monaghan, P. Compensation for a bad start: Grow now, pay later?. Trends Ecol. Evol. 16, 254–260. https://doi.org/10.1016/S0169-5347(01)02124-3 (2001).Article 
PubMed 

Google Scholar 
26.Criscuolo, F., Monaghan, P., Nasir, L. & Metcalfe, N. B. Early nutrition and phenotypic development: “catch-up” growth leads to elevated metabolic rate in adulthood. Proc. Biol. Sci. 275(1642), 1565–1570. https://doi.org/10.1098/rspb.2008.0148 (2008).Article 
PubMed 
PubMed Central 

Google Scholar 
27.Monaghan, P. Early growth conditions, phenotypic development and environmental change. Philos. Trans. R. Soc. B Biol. Sci. 363, 1635–1645. https://doi.org/10.1098/rstb.2007.0011 (2008).Article 

Google Scholar 
28.Lee, W. S., Monaghan, P. & Metcalfe, N. B. The pattern of early growth trajectories affects adult breeding performance. Ecology 93, 902–912. https://doi.org/10.1890/11-0890.1 (2012).Article 
PubMed 

Google Scholar 
29.Zera, A. J. & Harshman, L. G. The Physiology of Life History Trade-Offs in Animals. Annu. Rev. Ecol. Syst. 32, 95–126. https://doi.org/10.1146/annurev.ecolsys.32.081501.114006 (2001).Article 

Google Scholar 
30.Botías, C., David, A., Hill, E. M. & Goulson, D. Quantifying exposure of wild bumblebees to mixtures of agrochemicals in agricultural and urban landscapes. Environ. Pollut. 222, 73–82. https://doi.org/10.1016/j.envpol.2017.01.001 (2017).CAS 
Article 
PubMed 

Google Scholar 
31.Hladik, M. L. & Kolpin, D. W. First national-scale reconnaissance of neonicotinoid insecticides in streams across the USA. Environ. Chem. 13, 12. https://doi.org/10.1071/EN15061 (2016).CAS 
Article 

Google Scholar 
32.Morrissey, C. A. et al. Neonicotinoid contamination of global surface waters and associated risk to aquatic invertebrates: A review. Environ. Int. 74, 291–303. https://doi.org/10.1016/j.envint.2014.10.024 (2015).CAS 
Article 
PubMed 

Google Scholar 
33.McNabb, F. M. A. The hypothalamic-pituitary-thyroid (HPT) axis in birds and its role in bird development and reproduction. Crit. Rev. Toxicol. 37(1–2), 163–193. https://doi.org/10.1080/10408440601123552 (2007).CAS 
Article 
PubMed 

Google Scholar 
34.Gobeli, A., Crossley, D., Johnson, J. & Reyna, K. The effects of neonicotinoid exposure on embryonic development and organ mass in northern bobwhite quail (Colinus virginianus). Comp. Biochem. Physiol. Part – C Toxicol. Pharmacol. 195, 9–15. https://doi.org/10.1016/j.cbpc.2017.02.001 (2017).CAS 
Article 

Google Scholar 
35.Mineau, P. & Callaghan, C. Neonicotinoid insecticides and bats: an assessment of the direct and indirect risks. (Canadian Wildlife Federation, 2018).36.Wilson, J. D., Morris, A. J., Arroyo, B. E., Clark, S. C. & Bradbury, R. B. A review of the abundance and diversity of invertebrate and plant foods of granivorous birds in northern Europe in relation to agricultural change. Agric. Ecosyst. Environ. 75, 13–30. https://doi.org/10.1016/S0167-8809(99)00064-X (1999).Article 

Google Scholar 
37.Peig, J. & Green, A. J. New perspectives for estimating body condition from mass/length data: the scaled mass index as an alternative method. Oikos 118, 1883–1891. https://doi.org/10.1111/j.1600-0706.2009.17643.x (2009).Article 

Google Scholar 
38.Spencer, K., Buchanan, K., Goldsmith, A. & Catchpole, C. Song as an honest signal of developmental stress in the zebra finch (Taeniopygia guttata). Horm. Behav. 44, 132–139. https://doi.org/10.1016/S0018-506X(03)00124-7 (2003).CAS 
Article 
PubMed 

Google Scholar 
39.Ayyanath, M.-M., Cutler, G. C., Scott-Dupree, C. D. & Sibley, P. K. Transgenerational Shifts in Reproduction Hormesis in Green Peach Aphid Exposed to Low Concentrations of Imidacloprid. PLoS One 8, e74532. https://doi.org/10.1371/journal.pone.0074532 (2013).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
40.Calabrese, E. J. & Baldwin, L. A. Toxicology rethinks its central belief. Nature 421, 691–692. https://doi.org/10.1038/421691a (2003).ADS 
CAS 
Article 
PubMed 

Google Scholar 
41.Lopez-Antia, A., Ortiz-Santaliestra, M. E., Mougeot, F. & Mateo, R. Experimental exposure of red-legged partridges (Alectoris rufa) to seeds coated with imidacloprid, thiram and difenoconazole. Ecotoxicology 22, 125–138. https://doi.org/10.1007/s10646-012-1009-x (2013).CAS 
Article 
PubMed 

Google Scholar 
42.Rix, R. R., Ayyanath, M. M. & Christopher Cutler, G. Sublethal concentrations of imidacloprid increase reproduction, alter expression of detoxification genes, and prime Myzus persicae for subsequent stress. J. Pest Sci. (2004) 89, 581–589. https://doi.org/10.1007/s10340-015-0716-5 (2016).Article 

Google Scholar 
43.von Engelhardt, N. & Groothuis, T. G. G. Maternal hormones in avian eggs. in Hormones and Reproduction of Vertebrates: Birds, 1st edn. (eds Norris, D. & Lopez, K.) 91–127. https://doi.org/10.1016/C2009-0-01697-3 (Academic Press, 2011).Chapter 

Google Scholar 
44.Hulbert, A. J. Thyroid hormones and their effects: A new perspective. Biol. Rev. Camb. Philos. Soc. 75, 519–631. https://doi.org/10.1017/s146479310000556x (2000).CAS 
Article 
PubMed 

Google Scholar 
45.Darras, V. M. et al. Partial Food Restriction Increases Hepatic Inner Ring Deiodinating Activity in the Chicken and the Rat. Gen. Comp. Endocrinol. 100, 334–338. https://doi.org/10.1006/gcen.1995.1164 (1995).CAS 
Article 
PubMed 

Google Scholar 
46.Klandorf, H. & Harvey, S. Food intake regulation of circulating thyroid hormones in domestic fowl. Gen. Comp. Endocrinol. 60, 162–170. https://doi.org/10.1016/0016-6480(85)90310-7 (1985).CAS 
Article 
PubMed 

Google Scholar 
47.Reyns, G. E., Janssens, K. A., Buyse, J., Kühn, E. R. & Darras, V. M. Changes in thyroid hormone levels in chicken liver during fasting and refeeding. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 132(1), 239–245. https://doi.org/10.1016/s1096-4959(01)00528-0.CAS 
Article 
PubMed 

Google Scholar 
48.Harvey, S. & Klandorf, H. Reduced adrenocortical function and increased thyroid function in fasted and refed chickens. J. Endocrinol. 98, 129–135. https://doi.org/10.1677/joe.0.0980129 (1983).CAS 
Article 
PubMed 

Google Scholar 
49.Rimbach, R., Pillay, N. & Schradin, C. Both thyroid hormone levels and resting metabolic rate decrease in African striped mice when food availability decreases. J. Exp. Biol. 220, 837–843. https://doi.org/10.1242/jeb.151449 (2017).Article 
PubMed 

Google Scholar 
50.Scott, I. & Evans, P. R. The metabolic output of avian (Sturnus vulgaris, Calidris alpina) adipose tissue liver and skeletal muscle: Implications for BMR/body mass relationships. Comp. Biochem. Physiol. Comp. Physiol. 103(2), 329–332. https://doi.org/10.1016/0300-9629(92)90589-I (1992).CAS 
Article 
PubMed 

Google Scholar 
51.Mesnage, R., Biserni, M., Genkova, D., Wesolowski, L. & Antoniou, M. N. Evaluation of neonicotinoid insecticides for oestrogenic, thyroidogenic and adipogenic activity reveals imidacloprid causes lipid accumulation. J. Appl. Toxicol. 38, 1483–1491. https://doi.org/10.1002/jat.3651 (2018).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
52.Lindström, J. Early development and fitness in birds and mammals. Trends Ecol. Evol. 14(9), 343–348. https://doi.org/10.1016/S0169-5347(99)01639-0 (1999).Article 
PubMed 

Google Scholar 
53.Vézina, F., Love, O. P., Lessard, M. & Williams, T. D. Shifts in metabolic demands in growing altricial nestlings illustrate context-specific relationships between basal metabolic rate and body composition. Physiol. Biochem. Zool. 82, 248–257. https://doi.org/10.1086/597548 (2009).Article 
PubMed 

Google Scholar 
54.Swanson, D. L., Mckechnie, A. E. & Vézina, F. How low can you go ? An adaptive energetic framew ork for interpreting basal metabolic rate variation in endotherms. J. Comp. Physiol. B 187, 1039–1056. https://doi.org/10.1007/s00360-017-1096-3 (2017).Article 
PubMed 

Google Scholar 
55.Hao, C., Eng, M. L., Sun, F. & Morrissey, C. A. Part-per-trillion LC-MS/MS determination of neonicotinoids in small volumes of songbird plasma. Sci. Total Environ. 644, 1080–1087. https://doi.org/10.1016/j.scitotenv.2018.06.317 (2018).ADS 
CAS 
Article 
PubMed 

Google Scholar 
56.Taliansky-Chamudis, A., Gómez-Ramírez, P., León-Ortega, M. & García-Fernández, A. J. Validation of a QuECheRS method for analysis of neonicotinoids in small volumes of blood and assessment of exposure in Eurasian eagle owl (Bubo bubo) nestlings. Sci. Total Environ. 595, 93–100. https://doi.org/10.1016/j.scitotenv.2017.03.246 (2017).ADS 
CAS 
Article 
PubMed 

Google Scholar 
57.Lemon, W. C. The energetics of lifetime reproductive success in the zebra finch Taeniopygia guttata. Physiol. Zool. 66, 946–963. https://doi.org/10.1086/physzool.66.6.30163748 (1993).Article 

Google Scholar 
58.Chastel, O., Lacroix, A. & Kersten, M. Pre-breeding energy requirements: thyroid hormone, metabolism and the timing of reproduction in house sparrows (Passer domesticus). J. Avian Biol. 34, 298–306. https://doi.org/10.1034/j.1600-048X.2003.02528.x (2003).Article 

Google Scholar 
59.Hicks, O. et al. The role of parasitism in the energy management of a free-ranging bird. J. Exp. Biol. 221(24), jeb190066. https://doi.org/10.1242/jeb.190066 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
60.Guglielmo, C. G., McGuire, L. P., Gerson, A. R. & Seewagen, C. L. Simple, rapid, and non-invasive measurement of fat, lean, and total water masses of live birds using quantitative magnetic resonance. J. Ornithol. 152, 75–85. https://doi.org/10.1007/s10336-011-0724-z (2011).Article 

Google Scholar 
61.Le Pogam, A. et al. Wintering snow buntings elevate cold hardiness to extreme levels but show no changes in maintenance costs. Physiol. Biochem. Zool. 93, 417–433. https://doi.org/10.1086/711370 (2020).Article 
PubMed 

Google Scholar 
62.Lighton, J. R. B. Measuring Metabolic Rates, 2nd edn. https://doi.org/10.1093/oso/9780198830399.001.0001 (Oxford University Press, Oxford, 2018).Book 

Google Scholar 
63.Gessaman, J. A. & Nagy, K. A. Energy metabolism: Errors in gas-exchange conversion factors. Physiol. Zool. 61, 507–513. https://doi.org/10.1086/physzool.61.6.30156159 (1988).Article 

Google Scholar 
64.R Core Team. R: A Language and Environment for Statistical Computing. https://www.R-project.org/ (R Foundation for Statistical  Computing, Vienna, Austria, 2017).65.Zuur, A. F., Ieno, E. N. & Elphick, C. S. A protocol for data exploration to avoid common statistical problems. Methods Ecol. Evol. 1, 3–14. https://doi.org/10.1111/j.2041-210X.2009.00001.x (2010).Article 

Google Scholar  More

1.Gray, M., Gordon, J. & Brown, E. Geodiversity and the ecosystem approach: The contribution of geoscience in delivering integrated environmental management. Proc. Geol. Assoc. 124, 659–673 (2013).Article 

Google Scholar 
2.Gray, M. Valuing geodiversity in an ‘ecosystem services’ context. Scott. Geogr. J. 128, 177–194 (2012).Article 

Google Scholar 
3.Warren, A. & French, J. R. Habitat Conservation: Managing the Physical Environment (Wiley, Hoboken, 2001).
Google Scholar 
4.Gordon, J. E., Barron, H. F., Hansom, J. D. & Thomas, M. F. Engaging with geodiversity—Why it matters. Proc. Geol. Assoc. 123, 1–6 (2012).Article 

Google Scholar 
5.Hjort, J., Gordon, J. E., Gray, M. & Hunter, M. L. Why geodiversity matters in valuing nature’s stage. Conserv. Biol. 29, 630–639 (2015).PubMed 
Article 
PubMed Central 

Google Scholar 
6.Gray, M. Geodiversity: Valuing and Conserving Abiotic Nature 448 (Wiley, 2004).
Google Scholar 
7.Serrano, E. & Ruiz-Flano, P. Geodiversity. A theoretical and applied concept. Geogr. Helv. Jg 62, 140–147 (2007).Article 

Google Scholar 
8.Comer, P. J. et al. Incorporating geodiversity into conservation decisions. Conserv. Biol. 29, 692–701 (2015).PubMed 
Article 
PubMed Central 

Google Scholar 
9.Pătru-Stupariu, I. et al. Integrating geo-biodiversity features in the analysis of landscape patterns. Ecol. Indic. 80, 363–375 (2017).Article 

Google Scholar 
10.Chakraborty, A. & Gray, M. A call for mainstreaming geodiversity in nature conservation research and praxis. J. Nat. Conserv. 56, 125862 (2020).Article 

Google Scholar 
11.Poesen, J., Torri, D. & Bunte, K. Effects of rock fragments on soil erosion by water at different spatial scales: A review. CATENA 23, 141–166 (1994).Article 

Google Scholar 
12.Zhang, Y., Zhang, M., Niu, J., Li, H. & Xiao, R. Rock fragments and soil hydrological processes: Significance and progress. CATENA 147, 153–166 (2016).Article 

Google Scholar 
13.Xia, L. et al. Effects of rock fragment cover on hydrological processes under rainfall simulation in a semi-arid region of China. Hydrol. Process. 32, 792–804 (2018).ADS 
Article 

Google Scholar 
14.Lavee, H. & Poesen, J. W. A. Overland flow generation and continuity on stone-covered soil surfaces. Hydrol. Process. 5, 345–360 (1991).ADS 
Article 

Google Scholar 
15.Agassi, M. & Levy, G. Stone cover and rain intensity—Effects on infiltration, erosion and water splash. Soil Res. 29, 565–575 (1991).Article 

Google Scholar 
16.Mandal, U. K. et al. Soil infiltration, runoff and sediment yield from a shallow soil with varied stone cover and intensity of rain. Eur. J. Soil Sci. 56, 435–443 (2005).Article 

Google Scholar 
17.Cerdà, A. Effects of rock fragment cover on soil infiltration, interrill runoff and erosion. Eur. J. Soil Sci. 52, 59–68 (2001).Article 

Google Scholar 
18.Jury, W. A. & Bellantuoni, B. Heat and water movement under surface rocks in a field soil: I. Thermal effects. Soil Sci. Soc. Am. J. 40, 505–509 (1976).ADS 
Article 

Google Scholar 
19.Yuan, C., Lei, T., Mao, L., Liu, H. & Wu, Y. Catena soil surface evaporation processes under mulches of different sized gravel. CATENA 78, 117–121 (2009).CAS 
Article 

Google Scholar 
20.Poesen, J. & Lavee, H. Rock fragments in top soils: Significance and processes. CATENA 23, 1–28 (1994).Article 

Google Scholar 
21.Yizhaq, H., Stavi, I., Shachak, M. & Bel, G. Geodiversity increases ecosystem durability to prolonged droughts. Ecol. Complex. 31, 96–103 (2017).Article 

Google Scholar 
22.Stavi, I., Rachmilevitch, S. & Yizhaq, H. Geodiversity effects on soil quality and geo-ecosystem functioning in drylands. CATENA 176, 372–380 (2019).CAS 
Article 

Google Scholar 
23.Preisler, Y. et al. Mortality versus survival in drought-affected Aleppo pine forest depends on the extent of rock cover and soil stoniness. Funct. Ecol. 33, 901–912 (2019).Article 

Google Scholar 
24.Sauer, T. J. & Logsdon, S. D. Hydraulic and Physical Properties of Stony Soils in a Small Watershed. Soil Sci. Soc. Am. J. 66, 1947–1956 (2002).25.Arnau-Rosalén, E., Calvo-Cases, A., Boix-Fayos, C., Lavee, H. & Sarah, P. Analysis of soil surface component patterns affecting runoff generation. An example of methods applied to Mediterranean hillslopes in Alicante (Spain). Geomorphology 101, 595–606 (2008).ADS 
Article 

Google Scholar 
26.Ceacero, C. J., Díaz-Hernández, J. L., de Campo, A. D. & Navarro-Cerrillo, R. M. Soil rock fragment is stronger driver of spatio-temporal soil water dynamics and efficiency of water use than cultural management in holm oak plantations. Soil Tillage Res. 197, 104495 (2020).Article 

Google Scholar 
27.Burnett, M. R., August, P. V., Brown, J. H. & Killingbeck, K. T. The influence of geomorphological heterogeneity on biodiversity I. A patch-scale perspective. Conserv. Biol. 12, 363–370 (2008).Article 

Google Scholar 
28.Engelbrecht, B. et al. Drought sensitivity shapes species distribution patterns in tropical forests. Nature 447, 80–82 (2007).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
29.Stavi, I., Rachmilevitch, S. & Yizhaq, H. Small-scale geodiversity regulates functioning, connectivity, and productivity of shrubby, semi-arid rangelands. L. Degrad. Dev. 29, 205–209 (2018).Article 

Google Scholar 
30.Dubinin, V., Stavi, I., Svoray, T., Dorman, M. & Yizhaq, H. Hillslope geodiversity improves the resistance of shrubs to prolonged droughts in semiarid ecosystems. J. Arid Environ. 188, 104462 (2021).ADS 
Article 

Google Scholar 
31.Ochoa-Hueso, R. et al. Soil fungal abundance and plant functional traits drive fertile island formation in global drylands. J. Ecol. 106, 242–253 (2018).CAS 
Article 

Google Scholar 
32.Stavi, I., Rachmilevitch, S., Hjazin, A. & Yizhaq, H. Geodiversity decreases shrub mortality and increases ecosystem tolerance to droughts and climate change. Earth Surf. Process. Landforms 43, 2808–2817 (2018).ADS 
Article 

Google Scholar 
33.Suggitt, A. J. et al. Extinction risk from climate change is reduced by microclimatic buffering. Nat. Clim. Change 8, 713–717 (2018).ADS 
Article 

Google Scholar 
34.Bailey, J. J., Boyd, D. S. & Field, R. Models of upland species’ distributions are improved by accounting for geodiversity. Landscape Ecol. https://doi.org/10.1007/s10980-018-0723-z (2018).Article 

Google Scholar 
35.Lenoir, J. et al. Local temperatures inferred from plant communities suggest strong spatial buffering of climate warming across Northern Europe. Glob. Change Biol. 19, 1470–1481 (2013).ADS 
Article 

Google Scholar 
36.Lawler, J. J. et al. The theory behind, and the challenges of, conserving nature’s stage in a time of rapid change. Conserv. Biol. 29, 618–629 (2015).PubMed 
Article 
PubMed Central 

Google Scholar 
37.Nichols, W. F., Killingbeck, K. T. & August, P. V. The influence biodiversity of geomorphological heterogeneity: II. A landscape perspective. Soc. Conserv. Biol. 12, 371–397 (1998).Article 

Google Scholar 
38.Alahuhta, J., Toivanen, M. & Hjort, J. Geodiversity–biodiversity relationship needs more empirical evidence. Nat. Ecol. Evol. 4, 2–3 (2020).PubMed 
Article 

Google Scholar 
39.Evenari, M., Shanan, L., Tadmor, N. & Shkolnik, A. The Negev: The Challenge of a Desert (Harvard University Press, 1982).Book 

Google Scholar 
40.Akttani, H., Trimborn, P. & Ziegler, H. Photosynthetic pathways in Chenopodiaceae from Africa, Asia and Europe with their ecological, phytogeographical and taxonomical importance. Plant Syst. Evol. 206, 187–221 (1997).Article 

Google Scholar 
41.Harley, J. The Biology of Mycorrhiza (Leonard Hill, 1969).
Google Scholar 
42.Mejsti, V. K. & Cudlin, P. Mycorrhiza in some plant desert species in Algeria. Plant Soil 71, 363–366 (1983).Article 

Google Scholar 
43.Segoli, M., Ungar, E. D. & Shachak, M. Shrubs enhance resilience of a semi-arid ecosystem by engineering and regrowth. Ecohydrology 1, 330–339 (2008).Article 

Google Scholar 
44.Gilad, E., Von Hardenberg, J., Provenzale, A., Shachak, M. & Meron, E. Ecosystem engineers: From pattern formation to habitat creation. Phys. Rev. Lett. 93, 098105 (2004).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
45.Wright, J. P., Jones, C. G., Boeken, B. & Shachak, M. Predictability of ecosystem engineering effects on species richness across environmental variability and spatial scales. J. Ecol. 94, 815–824 (2006).Article 

Google Scholar 
46.Katra, I., Blumberg, D. G., Lavee, H. & Sarah, P. Spatial distribution dynamics of topsoil moisture in shrub microenvironment after rain events in arid and semi-arid areas by means of high-resolution maps. Geomorphology 86, 455–464 (2007).ADS 
Article 

Google Scholar 
47.Hoffman, O., de Falco, N., Yizhaq, H. & Boeken, B. Annual plant diversity decreases across scales following widespread ecosystem engineer shrub mortality. J. Veg. Sci. 27, 578–586 (2016).Article 

Google Scholar 
48.Shachak, M. et al. Woody species as landscape modulators and their effect on biodiversity patterns. Bioscience 58, 209–221 (2008).Article 

Google Scholar 
49.Madrigal-González, J., García-Rodríguez, J. A. & Alarcos-Izquierdo, G. Testing general predictions of the stress gradient hypothesis under high inter- and intra-specific nurse shrub variability along a climatic gradient. J. Veg. Sci. 23, 52–61 (2012).Article 

Google Scholar 
50.Boeken, B. & Shachak, M. The dynamics of abundance and incidence of annual plant species during colonization in a desert. Ecography (Cop.) 21, 63–73 (1998).Article 

Google Scholar 
51.Golodets, C. & Boeken, B. Moderate sheep grazing in semiarid shrubland alters small-scale soil surface structure and patch properties. CATENA 65, 285–291 (2006).Article 

Google Scholar 
52.Boeken, B. & Shachak, M. Desert plant communities in human-made patches-implications for management. Ecol. Appl. 4, 702–716 (1994).Article 

Google Scholar 
53.Hoffman, O., Yizhaq, H. & Boeken, B. Small-scale effects of annual and woody vegetation on sediment displacement under field conditions. CATENA 109, 157–163 (2013).Article 

Google Scholar 
54.Zaady, E., Arbel, S., Barkai, D. & Sarig, S. Long-term impact of agricultural practices on biological soil crusts and their hydrological processes in a semiarid landscape. J. Arid Environ. 90, 5–11 (2013).ADS 
Article 

Google Scholar 
55.Zaady, E., Stavi, I. & Yizhaq, H. Hillslope geodiversity effects on properties and composition of biological soil crusts in drylands. Eur. J. Soil Sci. https://doi.org/10.1111/ejss.13097 (2021).Article 

Google Scholar 
56.Feinbrun-Dothan, N. & Danin, A. Analytical Flora of Eretz-Israel (Cana Publishing Ltd, 1991).
Google Scholar 
57.Solovchenko, A., Merzlyak, M. N., Khozin-Goldberg, I., Cohen, Z. & Boussiba, S. Coordinated carotenoid and lipid syntheses induced in parietochloris incisa (chlorophyta, trebouxiophyceae) mutant deficient in δ5 desaturase by nitrogen starvation and high light. J. Phycol. 46, 763–772 (2010).CAS 
Article 

Google Scholar 
58.Shannon, C. A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423 (1948).MathSciNet 
MATH 
Article 

Google Scholar 
59.Simpson, E. Measurement of diversity. Nature 163, 688 (1949).ADS 
MATH 
Article 

Google Scholar 
60.R Core Team. R: A language and environment for statistical Computing. R Foundation for Statistical  Computing, Vienna, Austria (2020).61.Richerson, P. J. & Lum, K. Patterns of plant species diversity in California: Relation to weather and topography. Am. Nat. 116, 504–536 (1980).Article 

Google Scholar 
62.Kerr, J. T. & Packer, L. Habitat heterogeneity as a determinant of mammal species richness in high-energy regions. Nature 385, 252–254 (1997).ADS 
CAS 
Article 

Google Scholar 
63.Alahuhta, J. et al. The role of geodiversity in providing ecosystem services at broad scales. Ecol. Indic. 91, 47–56 (2018).Article 

Google Scholar 
64.Zarnetske, P. L. et al. Towards connecting biodiversity and geodiversity across scales with satellite remote sensing. Wiley Online Libr. 28, 548–556 (2019).
Google Scholar 
65.Schrodt, F. et al. To advance sustainable stewardship, we must document not only biodiversity but geodiversity. Proc. Natl. Acad. Sci. U. S. A. 116, 16155–161658 (2019).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
66.Read, Q. D. et al. Beyond counts and averages: Relating geodiversity to dimensions of biodiversity. Wiley Online Libr. 29, 696–710 (2020).
Google Scholar 
67.Antonelli, A. et al. Geological and climatic influences on mountain biodiversity. Nat. Geosci. 11, 718–725 (2018).ADS 
CAS 
Article 

Google Scholar 
68.Knudson, C., Kay, K. & Fisher, S. Appraising geodiversity and cultural diversity approaches to building resilience through conservation. Nat. Clim. Change 8, 678–685 (2018).ADS 
Article 

Google Scholar 
69.Beier, P., Hunter, M. L. & Anderson, M. Special section: Conserving nature’s stage. Conserv. Biol. 29, 613–617 (2015).PubMed 
Article 
PubMed Central 

Google Scholar 
70.Dubinin, V., Svoray, T., Stavi, I. & Yizhaq, H. Using LANDSAT 8 and VENµS data to study the effect of geodiversity on soil moisture dynamics in a semiarid shrubland. Remote Sens. 12, 3377 (2020).ADS 
Article 

Google Scholar 
71.Renne, R. R. et al. Soil and stand structure explain shrub mortality patterns following global change–type drought and extreme precipitation. Ecology 100, e02889 (2019).PubMed 
PubMed Central 

Google Scholar 
72.Gutterman, Y., Golan, T. & Garsani, M. Porcupine diggings as a unique ecological system in a desert environment. Oecologia 85, 122–127 (1990).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
73.Armas, C., Pugnaire, F. I. & Sala, O. E. Patch structure dynamics and mechanisms of cyclical succession in a Patagonian steppe (Argentina). J. Arid Environ. 72, 1552–1561 (2008).ADS 
Article 

Google Scholar 
74.Pickett, S. & White, P. The Ecology of Natural Disturbance and Patch Dynamics (Academic Press, 1985). https://doi.org/10.1016/C2009-0-02952-3.Book 

Google Scholar 
75.Segoli, M., Ungar, E. D., Giladi, I., Arnon, A. & Shachak, M. Untangling the positive and negative effects of shrubs on herbaceous vegetation in drylands. Landsc. Ecol. 27, 899–910 (2012).Article 

Google Scholar 
76.Rodríguez, F., Mayor, Á. G., Rietkerk, M. & Bautista, S. A null model for assessing the cover-independent role of bare soil connectivity as indicator of dryland functioning and dynamics. Ecol. Indic. 94, 512–519 (2018).Article 

Google Scholar 
77.Zelnik, Y. R., Kinast, S., Yizhaq, H., Bel, G. & Meron, E. Regime shifts in models of dryland vegetation. Philos. Trans. R. Soc. A https://doi.org/10.1098/rsta.2012.0358 (2013).Article 
MATH 

Google Scholar 
78.Walker, M. D. et al. Plant community responses to experimental warming across the tundra biome. PNAS 103, 1342–1346 (2006).79.Kardol, P. et al. Climate change effects on plant biomass alter dominance patterns and community evenness in an experimental old-field ecosystem. Glob. Change Biol. 16, 2676–2687 (2010).ADS 
Article 

Google Scholar 
80.Hillebrand, H., Bennett, D. M. & Cadotte, M. W. Consequences of dominance: A review of evenness effects on local and regional ecosystem processes. Ecology 89, 1510–1520 (2008).PubMed 
Article 
PubMed Central 

Google Scholar 
81.Stavi, I., Yizhaq, H., Szitenberg, A. & Zaady, E. Patch-scale to hillslope-scale geodiversity alleviates susceptibility of dryland ecosystems to climate change: Insights from the Israeli Negev. Curr. Opin. Environ. Sustain. 50, 129–137 (2021).Article 

Google Scholar 
82.Loarie, S. R. et al. Climate change and the future of California’s endemic flora. PLoS ONE 3, 2502 (2008).ADS 
Article 
CAS 

Google Scholar 
83.Loarie, S. R. et al. The velocity of climate change. Nature 462, 1052–1055 (2009).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
84.Ashcroft, M. B., Chisholm, L. A. & French, K. O. Climate change at the landscape scale: Predicting fine-grained spatial heterogeneity in warming and potential refugia for vegetation. Glob. Change Biol. 15, 656–667 (2009).ADS 
Article 

Google Scholar 
85.Correa-Metrio, A., Meave, J. A., Lozano-García, S. & Bush, M. B. Environmental determinism and neutrality in vegetation at millennial time scales. J. Veg. Sci. 25, 627–635 (2014).Article 

Google Scholar 
86.Baumgartner, J., Esperon-Rodriguez, M. & Beaumont, L. Identifying in situ climate refugia for plant species. Ecography (Cop.) 41, 1850–1863 (2018).Article 

Google Scholar 
87.Alahuhta, J. et al. The Role of Geodiversity in Providing Ecosystem Services at Broad Scales (Elsevier, 2018).Book 

Google Scholar 
88.Parks, K. E. & Mulligan, M. On the relationship between a resource based measure of geodiversity and broad scale biodiversity patterns. Biodivers. Conserv. 19, 2751–2766 (2010).Article 

Google Scholar 
89.Keppel, G. et al. The capacity of refugia for conservation planning under climate change. Front. Ecol. Environ. 13, 106–112 (2015).Article 

Google Scholar 
90.Mokany, K. et al. Past, present and future refugia for Tasmania’s palaeoendemic flora. J. Biogeogr. 44, 1537–1546 (2017).Article 

Google Scholar  More

1.Clavero, M. & Garcia-Berthou, E. Invasive species are a leading cause of animal extinctions. Trends Ecol. Evol. 20, 110–110. https://doi.org/10.1016/j.tree.2005.01.003 (2005).Article 
PubMed 

Google Scholar 
2.Grinnell, J. The niche-relationships of the California thrasher. Auk 34, 427–433 (1917).Article 

Google Scholar 
3.Elton, C. S. Animal Ecology and Evolution (Clarendon Press, 1930).
Google Scholar 
4.Arim, M., Abades, S. R., Neill, P. E., Lima, M. & Marquet, P. A. Spread dynamics of invasive species. Proc. Natl. Acad. Sci. U.S.A. 103, 374–378. https://doi.org/10.1073/pnas.0504272102 (2006).ADS 
CAS 
Article 
PubMed 

Google Scholar 
5.Kopf, R. K. et al. Confronting the risks of large-scale invasive species control. Nat. Ecol. Evol. 1, 0172. https://doi.org/10.1038/s41559-017-0172 (2017).Article 

Google Scholar 
6.Lonsdale, W. M. Global patterns of plant invasions and the concept of invasibility. Ecology 80, 1522–1536. https://doi.org/10.2307/176544 (1999).Article 

Google Scholar 
7.Pyšek, P. et al. MAcroecological framework for invasive aliens (MAFIA): Disentangling large-scale context dependence in biological invasions. Neobiota https://doi.org/10.3897/neobiota.62.52787 (2020).Article 

Google Scholar 
8.Donaldson, J. E. et al. Invasion trajectory of alien trees: the role of introduction pathway and planting history. Glob. Change Biol. 20, 1527–1537. https://doi.org/10.1111/gcb.12486 (2014).ADS 
Article 

Google Scholar 
9.Sax, D. F. & Brown, J. H. The paradox of invasion. Glob. Ecol. Biogeogr. 9, 363–371. https://doi.org/10.1046/j.1365-2699.2000.00217.x (2000).Article 

Google Scholar 
10.Valido, A., Rodriguez-Rodriguez, M. C. & Jordano, P. Honeybees disrupt the structure and functionality of plant-pollinator networks. Sci. Rep. 9, a4711. https://doi.org/10.1038/s41598-019-41271-5 (2019).ADS 
CAS 
Article 

Google Scholar 
11.Morales, C. L. & Aizen, M. A. Invasive mutualisms and the structure of plant–pollinator interactions in the temperate forests of north-west Patagonia, Argentina. J. Ecol. 94, 171–180. https://doi.org/10.1111/j.1365-2745.2005.01069.x (2006).Article 

Google Scholar 
12.Richardson, D. M., Allsopp, N., D’Antonio, C. M., Milton, S. J. & Rejmanek, M. Plant invasions—The role of mutualisms. Biol. Rev. 75, 65–93. https://doi.org/10.1017/S0006323199005435 (2000).CAS 
Article 
PubMed 

Google Scholar 
13.Simberloff, D. & Von Holle, B. Positive interactions of nonindigenous species: Invasional meltdown?. Biol. Invasions 1, 21–32. https://doi.org/10.1023/A:1010086329619 (1999).Article 

Google Scholar 
14.Vazquez, D. P. & Aizen, M. A. Asymmetric specialization: A pervasive feature of plant-pollinator interactions. Ecology 85, 1251–1257. https://doi.org/10.1890/03-3112 (2004).Article 

Google Scholar 
15.Shigesada, N. & Kawasaki, K. Biological Invasions: Theory and Practice (Oxford University Press, 1997).
Google Scholar 
16.Liebhold, A. M., Keitt, T. H., Goel, N. & Bertelsmeier, C. Scale invariance in the spatial-dynamics of biological invasions. Neobiota https://doi.org/10.3897/neobiota.62.53213 (2020).Article 

Google Scholar 
17.Mainali, K. P. et al. Projecting future expansion of invasive species: Comparing and improving methodologies for species distribution modeling. Glob. Change Biol. 21, 4464–4480. https://doi.org/10.1111/gcb.13038 (2015).ADS 
Article 

Google Scholar 
18.Barbet-Massin, M., Rome, Q., Villemant, C. & Courchamp, F. Can species distribution models really predict the expansion of invasive species?. PLoS ONE 13, e0193085. https://doi.org/10.1371/journal.pone.0193085 (2018).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
19.Dafni, A., Kevan, P., Gross, C. L. & Goka, K. Bombus terrestris, pollinator, invasive and pest: An assessment of problems associated with its widespread introductions for commercial purposes. Appl. Entomol. Zool. 45, 101–113. https://doi.org/10.1303/aez.2010.101 (2010).Article 

Google Scholar 
20.Velthuis, H. H. W. & van Doorn, A. A century of advances in bumblebee domestication and the economic and environmental aspects of its commercialization for pollination. Apidologie 37, 421–451. https://doi.org/10.1051/apido:2006019 (2006).Article 

Google Scholar 
21.Medel, R., González-Browne, C., Salazar, D. A., Ferrer, P. & Ehrenfeld, M. The most effective pollinator principle applies to new invasive pollinators. Biol. Lett. https://doi.org/10.1098/rsbl.2018.0132 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
22.Smith-Ramirez, C. et al. The reasons why Chile should stop importing commercial bumblebee Bombus terrestris (Linnaeus) and to start controlling it. Gayana 82, 118–127. https://doi.org/10.4067/S0717-65382018000200118 (2018).Article 

Google Scholar 
23.Aizen, M. A. et al. Coordinated species importation policies are needed to reduce serious invasions globally: The case of alien bumblebees in South America. J. Appl. Ecol. 56, 100–106. https://doi.org/10.1111/1365-2664.13121 (2019).Article 

Google Scholar 
24.Acosta, A. L., Giannini, T. C., Imperatriz-Fonseca, V. L. & Saraiva, A. M. Worldwide alien invasion: A methodological approach to forecast the potential spread of a highly invasive pollinator. PLoS ONE 11, e0148295. https://doi.org/10.1371/journal.pone.0148295 (2016).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
25.Geslin, B. & Morales, C. L. New records reveal rapid geographic expansion of Bombus terrestris Linnaeus, 1758 (Hymenoptera: Apidae), an invasive species in Argentina. CheckList 11, a1620. https://doi.org/10.15560/11.3.1620 (2015).Article 

Google Scholar 
26.Montalva, J., Sepulveda, V., Vivallo, F. & Silva, D. P. New records of an invasive bumble bee in northern Chile: Expansion of its range or new introduction events?. J. Insect Conserv. 21, 657–666. https://doi.org/10.1007/s10841-017-0008-x (2017).Article 

Google Scholar 
27.González-Varo, J. P. et al. Combined effects of global change pressures on animal-mediated pollination. Trends Ecol. Evol. 28, 524–530. https://doi.org/10.1016/j.tree.2013.05.008 (2013).Article 
PubMed 

Google Scholar 
28.Knapp, J. L., Becher, M. A., Rankin, C. C., Twiston-Davies, G. & Osborne, J. L. Bombus terrestris in a mass-flowering pollinator-dependent crop: A mutualistic relationship?. Ecol. Evol. 9, 609–618. https://doi.org/10.1002/ece3.4784 (2019).Article 
PubMed 

Google Scholar 
29.Nilsen, E. B., Pedersen, S. & Linnell, J. D. C. Can minimum convex polygon home ranges be used to draw biologically meaningful conclusions?. Ecol. Res. 23, 635–639. https://doi.org/10.1007/s11284-007-0421-9 (2008).Article 

Google Scholar 
30.Kadoya, T. & Washitani, I. Predicting the rate of range expansion of an invasive alien bumblebee (Bombus terrestris) using a stochastic spatio-temporal model. Biol. Conserv. 143, 1228–1235. https://doi.org/10.1016/j.biocon.2010.02.030 (2010).Article 

Google Scholar 
31.Kadoya, T., Ishii, H. S., Kikuchi, R., Suda, S. & Washitani, I. Using monitoring data gathered by volunteers to predict the potential distribution of the invasive alien bumblebee Bombus terrestris. Biol. Conserv. 142, 1011–1017. https://doi.org/10.1016/j.biocon.2009.01.012 (2009).Article 

Google Scholar 
32.Murúa, M., Espindola, A., González, A. & Medel, R. Pollinators and crossability as reproductive isolation barriers in two sympatric oil-rewarding Calceolaria (Calceolariaceae) species. Evol. Ecol. 31, 421–434. https://doi.org/10.1007/s10682-017-9894-3 (2017).Article 

Google Scholar 
33.Valdivia, C. E., Carroza, J. P. & Orellana, J. I. Geographic distribution and trait-mediated causes of nectar robbing by the European bumblebee Bombus terrestris on the Patagonian shrub Fuchsia magellanica. Flora 225, 30–36. https://doi.org/10.1016/j.flora.2016.09.010 (2016).Article 

Google Scholar 
34.Herbertsson, L. et al. Long-term data shows increasing dominance of Bombus terrestris with climate warming. Basic Appl. Ecol. 53, 116–123. https://doi.org/10.1016/j.baae.2021.03.008 (2021).Article 

Google Scholar 
35.Aizen, M. A. et al. When mutualism goes bad: Density-dependent impacts of introduced bees on plant reproduction. New Phytol. 204, 322–328. https://doi.org/10.1111/nph.12924 (2014).Article 

Google Scholar 
36.Esterio, G. et al. Assessing the impact of the invasive buff-tailed bumblebee (Bombus terrestris) on the pollination of the native Chilean herb Mimulus luteus. Arthropod-Plant Interact. 7, 467–474. https://doi.org/10.1007/s11829-013-9264-1 (2013).Article 

Google Scholar 
37.Morales, C. L., Arbetman, M. P., Cameron, S. A. & Aizen, M. A. Rapid ecological replacement of a native bumble bee by invasive species. Front. Ecol. Environ. 11, 529–534. https://doi.org/10.1890/120321 (2013).Article 

Google Scholar 
38.Plischuk, S., Antunez, K., Haramboure, M., Minardi, G. M. & Lange, C. E. Long-term prevalence of the protists Crithidia bombi and Apicystis bombi and detection of the microsporidium Nosema bombi in invasive bumble bees. Environ. Microbiol. Rep. 9, 169–173. https://doi.org/10.1111/1758-2229.12520 (2017).CAS 
Article 
PubMed 

Google Scholar 
39.Plischuk, S. & Lange, C. E. Invasive Bombus terrestris (Hymenoptera: Apidae) parasitized by a flagellate (Euglenozoa: Kinetoplastea) and a neogregarine (Apicomplexa: Neogregarinorida). J. Invertebr. Pathol. 102, 261–263. https://doi.org/10.1016/j.jip.2009.08.005 (2009).Article 

Google Scholar 
40.Plischuk, S., Meeus, I., Smagghe, G. & Lange, C. E. Apicystis bombi (Apicomplexa: Neogregarinorida) parasitizing Apis mellifera and Bombus terrestris (Hymenoptera: Apidae) in Argentina. Environ. Microbiol. Rep. 3, 565–568. https://doi.org/10.1111/j.1758-2229.2011.00261.x (2011).Article 
PubMed 

Google Scholar 
41.Ruz, L. & Herrera, R. Preliminary observations on foraging activities of Bombus dahlbomii and Bombus terrestris (Hymenoptera: Apidae) on native and non-native vegetation in Chile. Acta Hortic. 561, 165–169. https://doi.org/10.17660/ActaHortic.2001.561.24 (2000).Article 

Google Scholar 
42.Sáez, A., Morales, C. L., Garibaldi, L. A. & Aizen, M. A. Invasive bumble bees reduce nectar availability for honey bees by robbing raspberry flower buds. Basic Appl. Ecol. 19, 26–35. https://doi.org/10.1016/j.baae.2017.01.001 (2017).Article 

Google Scholar 
43.Sáez, A., Morales, J. M., Morales, C. L., Harder, L. D. & Aizen, M. A. The costs and benefits of pollinator dependence: Empirically based simulations predict raspberry fruit quality. Ecol. Appl. 28, 1215–1222. https://doi.org/10.1002/eap.1720 (2018).Article 
PubMed 

Google Scholar 
44.Schmid-Hempel, R. et al. The invasion of southern South America by imported bumblebees and associated parasites. J. Anim. Ecol. 83, 823–837. https://doi.org/10.1111/1365-2656.12185 (2014).Article 
PubMed 

Google Scholar 
45.Torretta, J. P., Medan, D. & Arahamovich, A. H. First record of the invasive bumblebee Bombus terrestris (L.) (Hymenoptera, Apidae) in Argentina. Trans. Am. Entomol. Soc. 132, 285–289 (2006).Article 

Google Scholar 
46.Arismendi, N., Bruna, A., Zapata, N. & Vargas, M. Molecular detection of the tracheal mite Locustacarus buchneri in native and non-native bumble bees in Chile. Insect Soc 63, 629–633. https://doi.org/10.1007/s00040-016-0502-2 (2016).Article 

Google Scholar 
47.Polidori, C. & Nieves-Aldrey, J. Comparative flight morphology in queens of invasive and native Patagonian bumblebees (Hymenoptera: Bombus). Crit. Biol. 338, 126–133. https://doi.org/10.1016/j.crvi.2014.11.001 (2015).Article 

Google Scholar 
48.Vieli, L., Davis, F. W., Kendall, B. E. & Altieri, M. Landscape effects on wild Bombus terrestris (Hymenoptera: Apidae) queens visiting highbush blueberry fields in south-central Chile. Apidologie 47, 711–716. https://doi.org/10.1007/s13592-015-0422-6 (2016).Article 

Google Scholar 
49.Sáez, A., Morales, C. L., Ramos, L. Y. & Aizen, M. A. Extremely frequent bee visits increase pollen deposition but reduce drupelet set in raspberry. J. Appl. Ecol. 51, 1603–1612. https://doi.org/10.1111/1365-2664.12325 (2014).Article 

Google Scholar 
50.Montalva, J., Dudley, L., Arroyo, M. K., Retamales, H. & Abrahamovich, A. H. Geographic distribution and associated flora of native and introduced bumble bees (Bombus spp.) in Chile. J. Apicult Res. 50, 11–21. https://doi.org/10.3896/Ibra.1.50.1.02 (2011).Article 

Google Scholar 
51.GBIF.org. GBIF Occurrence Download (April 15, 2020). https://doi.org/10.15468/dl.f15467jezh.52.R: A Language and Environment for Statistical Computing, Version 3.6.3 (Foundation for Statistical Computing, 2020). More

Data collectionWe systematically searched all peer-reviewed publications that were published prior to May 2021, which investigated the effects of plant diversity on terrestrial C:N:P ratios (i.e., plants, soils, soil microbial biomass, and extracellular enzymes) using the Web of Science (Core Collection; http://www.webofknowledge.com), Google Scholar (http://scholar.google.com), and the China National Knowledge Infrastructure (CNKI; https://www.cnki.net) using the search term: “C:N or C:P or N:P or C:N:P AND plant OR soil OR microbial biomass OR extracellular enzyme OR exoenzyme AND plant diversity OR richness OR mixture OR pure OR polyculture OR monoculture OR overyielding”, and also searched for references within these papers. Our survey also included studies summarized in previously published diversity-ecosystem functioning meta-analyses15,17,20,33. The literature search was performed following the guidelines of PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) (Moher, Liberati44; Supplementary Fig. 5).We employed the following criteria to select the studies: (i) they were purposely designed to test the effects of plant diversity on C:N:P ratios, (ii) they had at least one species mixture treatment and corresponding monocultures, (iii) they had the same initial climatic and soil properties in the monoculture and mixture treatment plots. In thirteen publications, several experiments, each with independent controls, were conducted at different locations and were considered to be distinct studies. In total, 169 studies met these criteria (Supplementary Fig. 5 and Supplementary Table 3). When different publications included the same data, we recorded the data only once. When a study included plant species mixtures of different numbers of species, we considered them as distinct observations.For each site, we extracted the means, the number of replications, and standard deviations of the C:N, N:P, and C:P ratios of plants (including leaves, shoots, fine roots, total roots), soils, soil enzymes as well as soil microbial biomass C:N ratios, if reported. Similar to Zhou and Staver45, we collected nine types of soil enzymes and integrated individual soil enzymes into combined enzymes to represent proxies targeting specific resource acquisitions: C-acquisition (average of Invertase, α-Glucosidase, β-1,4-Glucosidase, Cellobiohydrolase, β-1,4-Xylosidase), N-acquisition (average of β-1,4-N-acetylglucosaminidase, Leucine-aminopeptidase, Urease), and P-acquisition (phosphatase). The ratios of each type of enzyme were subsequently calculated, referred to as soil enzyme C:N, C:P, and N:P. When an original study reported the results graphically, we used Plot Digitizer version 2.0 (Department of Physics at the University of South Alabama, Mobile, AL, USA) to extract data from the figures. This resulted in 52 studies for plant C:N ratios, 35 studies for plant N:P ratios, 17 studies for plant C:P ratios, 83 studies for soil C:N ratios, 42 studies for soil N:P ratios, 19 studies for soil C:P ratios, 33 studies for soil microbial biomass C:N ratios, 41 studies for soil enzyme C:N ratios, 40 studies for soil enzyme N:P ratios and 34 studies for soil enzyme C:P ratios (Supplementary Table 3).We also extracted species compositions in mixtures, latitude, longitude, stand age, ecosystem type (i.e, forest, grassland, cropland, pot), mean annual temperature (MAT, °C), management practice (fertilization or not), soil type (FAO classification) and sampled soil depth from original or cited papers, or cited data sources. The mean annual aridity index and solar radiation data were retrieved from the CGIAR-CSI Global Aridity Index data set46 and WorldClim Version 247 using location information. The annual aridity index was calculated as the ratio of the mean annual precipitation to mean annual potential evapotranspiration48. Stand age (SA) was recorded as the number of years since stand establishment following stand-replacing disturbances in forests, and the number of years between the initiation and measurements of the experiments in grasslands, croplands, and pots. Observations were averaged if multiple measurements were conducted during different seasons within a year. The species proportions in plant mixtures were based on the stem density in forests and pots, coverage in croplands, and sown seeds in grasslands. Soil depth was recorded as the midpoint of each soil depth interval49. We employed the weighted averages of soil C:N, C:P, and N:P ratios of monocultures in each study as proxies for the status of background nutrients. For studies that did not report soil C:N, C:P, and N:P ratios of monocultures, we used the initial soil C:N, C:P, and N:P ratios (before experiment establishment, if reported) as proxies for the status of background nutrients. When a study reported the soils, soil microbial biomass or soil enzyme C:N:P data from multiple soil depths, we used the soil C:N, C:P, and N:P ratios of the corresponding depths as background nutrient proxies. For plant C:N:P data, we used the uppermost soil layer C:N, C:P, and N:P ratios as background nutrient proxies, since it contains the majority of the available nutrients essential for plant growth50. We compared the estimates for the data sets with and without pot studies and found that both data sets yielded qualitatively similar results (Supplementary Tables 2 and 4). Thus, we reported results based on the whole data set.We employed two key functional traits to describe the functional composition: ‘leaf nitrogen content per leaf dry mass’ (Nmass, mg g−1), and “specific leaf area” (SLA, mm2 mg−1; i.e., leaf area per leaf dry mass), as they are expected to be related to plant growth rate, resource uptake and use efficiency27, and are available for large numbers of species. We obtained the mean trait values of Nmass and SLA data by using all available measurements for each plant species from the TRY Plant Trait Database51 except for two studies that included the data in their original publication52, or related publications in the same sites53. Functional diversity (FDis) was calculated as functional dispersion, which is the mean distance of each species to the centroid of all species in the functional trait space, based on the two traits together54. The calculation of FDis was conducted using the FD package54.Data analysisThe natural log-transformed response ratio (lnRR) was employed to quantify the effects of plant mixture following Hedges, Gurevitch55:$${{{{{{mathrm{ln}}}}}}}{RR}={{{{{{mathrm{ln}}}}}}}({bar{X}}_{{{{{{mathrm{t}}}}}}}/{bar{X}}_{{{{{{mathrm{c}}}}}}})={{{{{{{mathrm{ln}}}}}}}bar{X}}_{{{{{{mathrm{t}}}}}}}-{{{{{{{mathrm{ln}}}}}}}bar{X}}_{{{{{{mathrm{c}}}}}}}$$
(1)
where ({overline{X}}_{{{{{{rm{t}}}}}}}) and ({overline{X}}_{{{{{{rm{c}}}}}}}) are the observed values of a selected variable in the mixture and the expected value of the mixture in each study, respectively. If a study has multiple richness levels in mixtures (for example, 1, 4, 8, and 16), lnRR was calculated for the species richness levels 4, 8, and 16, respectively. We calculated ({overline{X}}_{{{{{{rm{c}}}}}}}) based on weighted values of the component species in monocultures following Loreau and Hector39:$$overline{{X}_{{{{{{mathrm{c}}}}}}}}=sum ({p}_{i}times {m}_{i})$$
(2)
where mi is the observed value of the selected variable of the monoculture of species i and pi is the proportion of species i density in the corresponding mixture. When a study reported multiple types of mixtures (species richness levels) and experimental years, ({overline{X}}_{{{{{{rm{t}}}}}}}) and ({overline{X}}_{{{{{{rm{c}}}}}}}) were calculated separately for each mixture type and experimental year.In our data set, sampling variances were not reported in 37 of the 169 studies, and no single control group mean estimate is present with standard deviation or the standard error reported. Like the previous studies6,56, we employed the number of replications for weighting:$${W}_{{{{{{mathrm{r}}}}}}}=({N}_{{{{{{mathrm{c}}}}}}}times {N}_{{{{{{mathrm{t}}}}}}})/({N}_{{{{{{mathrm{c}}}}}}}+{N}_{{{{{{mathrm{t}}}}}}})$$
(3)
where Wr is the weight associated with each lnRR observation, and Nc and Nt are the number of replications in monocultures and corresponding mixtures, respectively.The C:N, N:P, and C:P ratios of plants, soils, and soil enzymes, as well as soil microbial biomass C:N ratios were considered as response variables and analyzed separately. To validate the linearity assumption for the continuous predictors, we initially graphically plotted the lnRR vs. individual predictors, including FDis, SA, and background nutrient status (N, i.e., C:N, C:P, and N:P ratios of soil) and identified logarithmic functions as an alternative to linear functions. We also statistically compared the linear and logarithmic functions with the predictor of interest as the fixed effect, and “study” and measured plant parts (i.e., leaves, shoots, fine roots, total roots) or soil depth as the random effects, using Akaike information criterion (AIC). The random factors were used to account for the autocorrelation among observations within each “Study”, and potential influences of variation in measured plant parts and soil depth. We found that the linear FDis, SA, and N resulted in lower, or similar AIC values (∆AIC  More

AOA kinetic propertiesIn this study we investigated the kinetic properties of 12 AOA strains, including representatives from all four described AOA phylogenetic lineages: Nitrosopumilales (Group I.1a), ‘Ca. Nitrosotaleales’ (Group I.1a-associated), Nitrososphaerales (Group I.1b), and ‘Ca. Nitrosocaldales’ (thermophilic AOA clade) [58, 59] (Fig. 1). These AOA isolates and enrichments were obtained from a variety of habitats (marine, soil, sediment, hot spring) and have optimal growth pH and temperatures ranging from 5.3–7.8 to 25–72 °C, respectively (Table S2). The substrate-dependent oxygen consumption rates for all AOA tested followed Michaelis–Menten kinetics. Below, the kinetic properties of these AOA are put into a broader context with comparisons to previously characterized AOM. It is important to note that the whole cell kinetic properties, such as substrate competitiveness, detailed here were generated from instantaneous activity measurements in the absence of growth. It is unknown how the substrate competitiveness of nitrifiers may or may not differ from their competitiveness when cellular processes such as growth, division, stress, and repair are involved.Fig. 1: Phylogenetic reconstruction of ammonia oxidizing archaea (AOA) rooted on closely related non-AOA members of the “Thaumarchaeota”.Black taxon labels correspond to AOA from cultures or enrichments. Gray taxon labels correspond to representative metagenome assembled genomes from release 05-RS95 of the genome taxonomy database [41]. AOA that were kinetically characterized in the current study are highlighted in gray and AOA that were previously characterized are indicated with an asterisk (*). The phylogeny was calculated with IQ-TREE under model LG + F + R6 using an alignment of 34 universal genes (43 markers) produced by CheckM [42]. Support values (UFboot) greater than 95% for bipartitions are shown with a black circle and support values between 80% and 95% are shown with a gray circle. Order designations reflect lineages proposed by Alves et al. [59]. The scale bar indicates amino acids changes per site.Full size imageNitrosopumilales (Group I.1a)From this lineage, three mesophilic marine (N. piranensis D3C, N. adriaticus NF5, and N. maritimus SCM1) [3, 60], two agricultural soil (N. koreense MY1 and ‘Ca. N. chungbukensis’ MY2) [61, 62] and one thermal spring isolate (‘Ca. N. uzonensis’ N4) [40] were kinetically characterized (Fig. S1). These AOA all displayed a high substrate affinity for NH3, ranging from ~2.2 to 24.8 nM. Thus, all characterized Nitrosopumilales, and not just marine isolates, are adapted to oligotrophic conditions. All possess substrate affinities several orders of magnitude higher (lower Km(app)) than any characterized AOB, with the exception of the recently characterized acidophilic gammaproteobacterial AOB ‘Ca. Nitrosacidococcus tergens’ [55] (Fig. 2a). This finding appears to support the widely reported hypothesis that regardless of the environment, AOA in general are adapted to lower substrate concentrations than AOB [22, 29, 30]. However, as described later, this trend does not apply to all AOA.Fig. 2: Substrate-dependent oxidation kinetics of ammonia-oxidizing microorganisms.The (a) apparent substrate affinity (Km(app)) for NH3, (b) specific substrate affinity (a°) for NH3, (c) Km(app) for total ammonium, (d) a° for total ammonium, and (e) maximum oxidation rate (Vmax), of AOA (red), comammox (blue), and AOB (black) are provided. Symbols filled with light gray represent previously published values from reference studies (references provided in Materials and Methods). The four different gradations of red differentiate the four AOA phylogenetic lineages: (I) Nitrosopumilales, (II) ‘Ca. Nitrosotaleales’, (III) Nitrososphaerales, and (IV) ‘Ca. Nitrosocaldales’. Measurements were performed with either pure (circles) or enrichment (diamonds) cultures. Multiple symbols per strain represent independent measurements performed in this study and/or in the literature. The individual Michaelis–Menten plots for each AOM determined in this study are presented in Figs. S1, S3–5, and S8. Note the different scales.Full size imageAs the substrate oxidation kinetics of the marine AOA strain, N. maritimus SCM1, originally characterized by Martens-Habbena et al. [29] have recently been disputed [63], they were revisited in this study (Fig. S2). With the same strain of N. maritimus used in Hink et al. [63] (directly obtained by the authors), we were able to reproduce (Figs. S1 and S2) the original kinetic properties of N. maritimus SCM1 reported in Martens-Habbena et al. [29] ruling out strain domestication during lab propagation as cause for the observed discrepancy. Therefore, the reported differences in the literature possibly reflect the measurements of two distinct cellular properties, Km(app) [29] and Ks [63], representing the half saturation of activity and growth, respectively. In addition, differences in pre-measurement cultivation and growth conditions could also contribute to these unexpected differences [63, 64]. More details are provided in the Supplementary Results and Discussion.‘Ca. Nitrosotaleales’ (Group I.1a-associated)The only isolated AOA strains in this lineage ‘Ca. Nitrosotalea devanaterra’ Nd1 and ‘Ca. Nitrosotalea sinensis’ Nd2, are highly adapted for survival in acidic environments and grow optimally at pH 5.3 [25, 65]. Both display a relatively low affinity for total ammonium (Km(app) = 3.41–11.23 μM), but their affinity for NH3 is among the highest calculated of any AOA characterized (Km(app) = ~0.6–2.8 nM) (Fig. 2a,c, and Fig. S3). This seemingly drastic difference in substrate affinity for total ammonium versus NH3 is due to the combination of the high acid dissociation constant of ammonium (pKa = 9.25) and the kinetic properties of these strains being carried out at pH 5.3. The very limited availability of NH3 under acidic conditions has led to the hypothesis that these acidophilic AOA should be highly adapted to very low NH3 concentrations and possess a high substrate affinity (low Km(app)) for NH3 [66, 67]. Our data corroborate this hypothesis.Nitrososphaerales (Group I.1b)The AOA strains ‘Ca. N. nevadensis’ GerE (culture information provided in Supplementary Results and Discussion), ‘Ca. N. oleophilus’ MY3 [68] and ‘Ca. N. franklandus’ C13 [69] were kinetically characterized, and contextualized with the previously published kinetic characterization of Nitrososphaera viennensis EN76 and ‘Ca. Nitrososphaera gargensis’ [5]. Together, the Nitrososphaerales AOA possess a wide range of affinities for NH3 (Km(app) = ~0.14–31.5 µM) (Fig. 2a and Fig. S4). Although this range of NH3 affinities spans more than two orders of magnitude, none of the Nitrososphaerales AOA possess an affinity for NH3 as high as any Nitrosopumilales or ‘Ca. Nitrosotaleales’ AOA (Fig. 2a).The moderately thermophilic enrichment culture ‘Ca. N. nevadensis’ GerE displayed a higher substrate affinity (lower Km(app)) for NH3 (0.17 ± 0.03 µM) than the other characterized AOA strains within the genus Nitrososphaera (Fig. 2a). In contrast, ‘Ca. N. oleophilus’ MY3 and ‘Ca. N. franklandus’ C13, which belong to the genus Nitrosocosmicus, had the lowest affinity (highest Km(app)) for NH3 (12.37 ± 6.78 μM and 16.32 ± 14.11 μM, respectively) of any AOA characterized to date. In fact, their substrate affinity is comparable to several characterized AOB (Fig. 2a). In this context it is interesting to note that several Nitrosocosmicus species have been shown to tolerate very high ammonium concentrations [68,69,70], a trait usually associated with AOB [24, 54]. The low substrate affinity observed in Nitrosocosmicus AOA correlates with the absence of a putative Amt-type high affinity ammonium transporter in the genome of any sequenced Nitrosocosmicus species to date [68, 69, 71].
‘Ca. Nitrosocaldales’ (Thermophilic AOA lineage)The thermophilic AOA enrichment cultures, ‘Ca. Nitrosocaldus yellowstonensis’ HL72 [72] and ‘Ca. N. tenchongensis’ DRC1 (culture information provided in Supplementary Results and Discussion), possess affinities for NH3 (Km(app) = ~1.36 ± 0.53 μM and ~0.83 ± 0.01 μM; respectively comparable to AOA within the genus Nitrososphaera (Fig. 2a). Notably, the substrate oxidation rate of these two AOA quickly dropped with increasing substrate concentrations after Vmax was reached (Fig. S5). This trend was not observed with any other AOA tested here and may reflect an increased susceptibly to NH3 stress at high temperatures, as the free NH3 concentration increases with increasing temperatures [33]. It should be noted that both of these AOA cultures are enrichment cultures, as no member of the ‘Ca. Nitrosocaldales’ has been isolated to date.Together, these results highlight that the substrate affinity for NH3 among AOA species is much more variable than previously hypothesized, spanning several orders of magnitude and in some cases overlapping with the substrate affinity values of characterized non-oligotrophic AOB. In addition, the substrate affinity of AOA is related, to a certain degree, to their phylogenetic placement within each of the four AOA phylogenetic lineages mentioned above (Fig. 2). Although the substrate affinity ranges of these AOA lineages overlap, the link between AOA phylogeny and kinetic properties provides deeper insights into the physiological and evolutionary differences among AOA species. As a limited number of AOA have been isolated and characterized to date, the continued isolation and characterization of AOA from underrepresented phylogenetic lineages and new habitats is needed. While substrate affinity is certainly one of multiple factors that contribute to niche differentiation between AOM in general, it may also present a previously under acknowledged factor in AOA niche differentiation.Maximum substrate oxidation rates (V
max)The normalized maximum substrate oxidation rate of all the AOA characterized to date only span about one order of magnitude from 4.27 to 54.68 μmol N mg protein−1 h−1. These normalized AOA Vmax values are in the same range as the recorded Vmax for the comammox N. inopinata (~12 μmol N mg protein−1 h−1) and the marine AOB strain Nitrosococcus oceani ATCC 19707 (~38 μmol N mg protein−1 h−1) but are lower than the normalized Vmax of the AOB Nitrosomonas europaea ATCC 19718 (average of 84.2 μmol N mg protein−1 h−1; Fig. 2e). The high Vmax value for N. europaea is the only real outlier among the AOM characterized to date and it remains to be determined whether other AOB related to N. europaea also possess such a high Vmax or if members of the Nitrosomonadales possess a broad range of Vmax values. Similarly, as additional comammox strains become available as pure cultures their kinetic characterization will be vital in understanding the variability of these ecologically important parameters within this guild.Specific substrate affinity (a°)Although the Km(app) and Vmax of AOM can be compared by themselves and provide useful information on cellular properties, the ability of an AOM to scavenge (and compete for) substrate from a dilute solution is most appropriately represented by the a°, which takes into account both the cellular Km(app) and Vmax [28]. In previous studies, the a° of AOM has been calculated using the Km(app) value for total ammonium (NH3 + NH4+) and not the Km(app) value for NH3 [5, 29]. Calculating the a° based on the Km(app) value for total ammonium allows for the a° of AOM to be compared with the a° of microorganisms that do not use NH3 as a sole energy generating substrate, such as ammonium assimilating heterotrophic bacteria or diatoms [29]. While this is useful when evaluating competition for total ammonium in mixed communities or environmental settings, an a° calculated using the Km(app) value for NH3 may be more useful when directly comparing the interspecies competitiveness of AOM for the following reasons: (i) our data support the hypothesis that the substrate for all AOM is NH3 and not NH4+ (see below) and (ii) the Km(app) value for total ammonium is more dependent on the environmental factors it was measured at (e.g., pH, temperature, salinity) than the Km(app) for NH3.All characterized AOA (with the exception of representatives of the genus Nitrosocosmicus) and the comammox bacterium N. inopinata possess much higher a° for total ammonium or NH3 (~10–3000×) than the AOB, N. oceani or N. europaea (Fig. 2b–d), indicating that they are highly competitive in environments limited in either total ammonium or only NH3. However, due to the low number of published normalized Vmax values for AOB, a° could only be calculated for these two AOB representatives. Thus, extrapolations to the a° of all AOB species, based solely on these observations should be approached with caution.The low variation in experimentally measured Vmax values (Fig. 2e) across all measured AOM in combination with the high variation in Km(app) values leads to a strong relationship between cellular a° and the reciprocal of Km(app) (Fig. 3) according to Eq. 2 (see Materials and Methods). AOM adapted to oligotrophic (low substrate) conditions should possess both a high substrate affinity (low Km(app)) and a high ao [28]. Therefore, the AOM best suited for environments limited in total ammonium are the AOA belonging to the Nitrosopumilales and the comammox isolate N. inopinata, (top right corner of Fig. 3a). Overall, when looking at solely NH3 or total ammonium, the separation of species in these plots remains similar, with the exception that the acidophilic AOA belonging to the ‘Ca. Nitrosotaleales’ are predicted to be best suited for life in environments limited in NH3 (Fig. 3b). The adaptation correlates well with the fact the AOA ‘Ca. Nitrosotalea devanaterra’ Nd1 and ‘Ca. Nitrosotalea sinensis’ Nd2 were isolated from acidic soils with a pH of 4.5 and 4.7, respectively [25, 65], where the NH3:NH4+ equilibrium is heavily shifted toward NH4+.Fig. 3: The reciprocal relationship between the substrate affinity (Km(app)) and specific substrate affinity (a°) of ammonia-oxidizing microorganisms (AOM).Reciprocal plots for both (a) total ammonium and (b) NH3 are depicted. The Km(app) and a° values correspond to the values presented for pure AOM isolates in Fig. 2. Data for AOA (red), comammox (blue), and AOB (black) are shown. The correlation (R2) indicates the linear relationship between the logarithmically transformed data points.Full size imageIn either case, when looking at NH3 or total ammonium, the AOA belonging to the genus Nitrosocosmicus (‘Ca. N. oleophilus’ MY3 and ‘Ca. N. franklandus’ C13) and AOB populate the lower left section of these plots, indicating that they are not strong substrate competitors in NH3 or total ammonium limited environments (Fig. 3). Here, the Vmax of all the AOM reported spans ~10×, whereas the difference in Km(app) spans about five orders of magnitude. If the cellular kinetic property of Vmax really is so similar across all AOB, AOA, and comammox species (Fig. 2e) compared to the large differences in Km(app) values, then substrate competitiveness can be predicted from an AOMs Km(app) for either NH3 or total ammonium (Fig. 2a–c). This may prove especially helpful when characterizing enrichment cultures, where normalizing ammonia-oxidizing activity to cellular protein in order to obtain a comparable Vmax value is not possible. However, there is also a need for more kinetically characterized AOB and comammox species to confirm this hypothesis. In addition, when comparing AOM, differences in the Vmax cellular property will play a larger role, the closer the Km(app) values of the AOM strains are. This is important to consider when comparing AOM from similar habitats and likely adapted to similar substrate concentrations.The effect of environmental and cellular factors on AOA kinetic propertiesThe concentration of NH3 present in a particular growth medium or environment can vary by orders of magnitude, based solely on the pH, temperature, or salinity of the system [73]. This is notable because at a given total ammonium concentration, the concentration of NH3 is ~10 times higher at 70 °C versus 30 °C and ~1000 times lower at pH 5.3 versus pH 8.4 (representative of maximum ranges tested). While it should be recognized that in our dataset no AOM were included that have a pH optimum between 5.3 and 7.0, the effect of pH and temperature on the ammonia oxidation kinetics of AOM must be considered in order to understand their ecophysiological niches. However, there was no correlation between the kinetic properties of AOM (Km(app), Vmax, and a°) measured in this study and their optimal growth temperature or pH. This lack of correlation between AOM species kinetic properties and growth conditions does not imply that the cellular kinetic properties of an individual AOM species will remain the same over a range of pH and temperature conditions. Therefore, we investigated the effect of pH and temperature variation on the substrate-dependent kinetic properties of the AOA strain ‘Ca. N. oleophilus’ MY3, and the effect of pH on the comammox strain N. inopinata. Here, the AOA ‘Ca. N. oleophilus’ MY3 was selected based on the fact that it is a non-marine, mesophilic, pure culture, that does not require external hydrogen peroxide scavengers for growth. These traits are shared with the previously characterized AOB, N. europaea [35], and the comammox organism, N. inopinata (this study) and thus facilitate comparison.The effect of temperatureThe effects of short-term temperature changes on the substrate-dependent kinetic properties of ‘Ca. N. oleophilus’ MY3 were determined. Temperature shifts of 5 °C above and below the optimal growth temperature (30 °C) had no effect on the Km(app) for total ammonium. However, the Km(app) for NH3, Vmax, and a° of ‘Ca. N. oleophilus’ MY3 all increased with increasing temperatures (Fig. S6). Therefore, as temperature increased, ‘Ca. N. oleophilus’ MY3 displayed a lower substrate affinity (higher Km(app) for NH3) but would be able to turnover substrate with a higher Vmax and better compete for substrate with a higher a°. Increasing AOA Km(app) values for NH3 with increasing temperatures have also been observed across studies with N. viennensis EN76 (Fig. S2), and this is discussed in more detail in the Supplementary Results and Discussion. In addition, similar observations have previously been made for AOB strains belonging to the genus Nitrosomonas [33, 34]. The increase in Vmax and a° can be explained in terms of the Van’t Hoff rule (reaction velocity increases with temperature) [74], or in terms of a temperature sensitivity coefficient (Q10; change in reaction velocity over 10 °C) [75]. Here, the maximal reaction velocity of ‘Ca. N. oleophilus’ MY3, displays a relative Q10 of 2.17 between 25 and 35 °C, which is in line with more general microbial respiration measurements [75, 76].The increase in Km(app) for NH3 (lower NH3 affinity) with increasing temperature is less straightforward to interpret. As this is a whole cell measurement, the observed differences may result from either broad cellular changes or from changes in individual enzymes involved in the ammonia oxidation pathway specifically. At the cellular level, changes in the proteinaceous surface layer (S-layer) or lipid cell membrane could affect substrate movement/transport and enzyme complex stability. It has been suggested that the negatively charged AOA S-layer proteins act as a substrate reservoir, trapping NH4+ and consequently increasing the NH3 concentration in the AOA pseudo-periplasmic space [77]. It is interesting to note that sequenced representatives from the genus ‘Ca. Nitrosocosmicus’ lack the main S-layer protein (slp1) found in all Nitrosopumilales, Nitrososphaerales, and ‘Ca. Nitrosotaleales’ sequenced isolates [71], although it remains to be demonstrated whether ‘Ca. Nitrosocosmicus’ members actually lack a S-layer or form S-layers composed of other proteins. In addition, it has been demonstrated that elevated temperatures significantly alter the lipid composition in the AOA cell membrane [78, 79]. However, it is unclear how differences in the cell membrane or S-layer composition between AOA species may affect the observed kinetic properties. In this context it is important to note that on the single enzyme level, previous studies have shown the same trend of decreasing substrate affinity and increasing maximal reaction velocity with increasing temperatures, due to altered protein structures and an increased enzyme-substrate dissociation constant [80, 81].Notably, differing optimum growth and activity conditions were previously determined for the marine AOB strain Nitrosomonas cryotolerans [34]. These observations raise interesting, albeit unanswered, questions about why the growth and activity temperature optima are or can be uncoupled in AOM, and what this means for AOM niche differentiation and their competitiveness in-situ. Moving forward, investigations into the growth and cellular kinetic properties of AOM across a range of environmental factor gradients will be essential in understanding competition between AOM in engineered and environmental systems.The effect of pHThe effects of short-term pH changes on the substrate-dependent kinetics of ‘Ca. N. oleophilus’ MY3 and N. inopinata were determined. The Vmax of both ‘Ca. N. oleophilus MY3’ and N. inopinata were stable at 37.3 ± 6.6 μmol N mg protein−1 h−1 and 11.2 ± 2.5 μmol N mg protein−1 h−1, respectively, in medium with a pH between ~6.5 and ~8.5 (Table S3). The Km(app) for total ammonium of ‘Ca. N. oleophilus MY3’ and N. inopinata decreased by more than an order of magnitude (~11×) across this pH range, while the Km(app) for NH3 remained more stable, increasing only 3–4 times (Fig. 4). This stability of the Km(app) for NH3 compared with the larger change in the Km(app) for total ammonium across this pH range suggests that the actual substrate used by AOA and comammox is indeed the undissociated form (NH3) rather than the ammonium ion (NH4+), as previously demonstrated for AOB [34, 35, 54, 82]. As these kinetic measurements were performed with whole cells, the change in Km(app) for NH3 across this pH range may be due to cellular effects of the differing pH values unrelated to the direct ammonia oxidation pathway. The changes in Km(app) for NH3 and Km(app) for total ammonium demonstrated here for ‘Ca. N. oleophilus’ MY3 and N. inopinata are similar to what has been observed for AOB. That AOA and AOB utilize the NH3 as a substrate, aligns with the fact that both are competitively inhibited by the non-polar acetylene compound [83, 84].Fig. 4: The effect of medium pH on the substrate affinity of ‘Ca. N. oleophilus MY3’ and N. inopinata.The substrate affinities for both (a,b) NH3 and (c,d) total ammonium (NH3 + NH4+) are provided. Individual substrate affinity values determined at each pH are shown as single points (circles). The boxes represent the first and third quartiles (25–75%) of the substrate affinity range under each condition. The median (line within the boxes) and mean substrate affinity (black diamonds) values are also indicated. The whiskers represent the most extreme values within 1.58× of quartile range. The variation of the substrate affinity values across the entire tested pH range are indicated in each panel. In all four instances there was a significant difference between the affinity at the lowest pH and the highest pH, as determined by a Student’s t test (p  More

Biostimulation of ureolytic communitiesEnrichmentThe native communities of the soils were successfully grown in the enrichment media (NB5U). The cultivated communities after two subcultures were serially diluted (10–2 to 10–6) and were spread with a sterile loop over Nutrient agar plates supplemented with 2% urea. Later, 36 morphologically distinct single colonies were obtained on the urea agar base plates based on visual observation.Isolation, identification, and characterization of the ureolytic isolatesOut of 36 isolated bacteria, six isolates (BS1, BS2, BS3, BS4, LS1, and LS2) were selected after checking for the urease activity test on the urea agar base (UAB) plate. These selected isolates turned the color of the UAB plate from orange to pink within 12 h. The 16 s rRNA sequence revealed the isolates as relatives of Sporosarcina pasteurii (SP). The details of the identified isolates are provided by NCMR (details in Supplementary Table 1). The biochemical characterization (details in Supplementary Table 2) of the isolates revealed that all the isolates are Gram-positive. All the isolates were rod-shaped, endospore-forming, urease, and oxidase positive. All the isolates were not able to utilize the Lysine and ONPG, contrary to SP.Further investigation of the isolated sequence was done via the NCBI database. The sequences were submitted to the GenBank database of the NCBI (National Center for Biotechnology Information) under the accession number MW024144 to MW024149. The BLAST analysis suggested that these strains are close relatives and indicate the possibility of being novel strains of the Sporosarcina family. We found that the isolate BS1 and BS2 had 96.62% (coverage 100%) and 96.22% (coverage 99%) identity with Sporosarcina siberinisis (NCBI accession number NR 134188). BS3 had 98.8% identity (coverage 97%) with Sporosarcina pasteurii (NCBI accession number NR 104923). BS4 and LS1 had 97.4% (coverage 99%) and 97.37% identity (coverage 100%) with Sporosarcina soli (NCBI accession number NR 043527). Contrarily, LS2 was found to be related to the Pseudogracilibacillus family. LS2 was observed to be closely related to Pseudogracilibacillus auburnensis P-207 with 97.06% identity (96% coverage). Based on these findings, the Phylogenetic tree was constructed with bootstrap (1000 replicates) considering the reference sequences obtained from the BLAST analysis, as shown in Fig. 3. The threshold criteria to differentiate two species is defined as 98.65% similarity score with the reference culture from databank40, while another study has suggested that in case of similarity index is  > 99%, the unknown isolate should be assigned to a species, and if the unknown isolates have similarity score between 95 to 99% to a reference sequence, the isolate should be assigned to the genus41. However, further investigation is suggested to conclude if the reported strains are novel or merely mutants of the reference strains of the databank. Similar observations were made at Graddy et al.22, where the majority of the isolated strains (47 out of 57) from bio-stimulation soil tanks were found to be strains of the Sporosarcina genus. It is worth noting that the soil enrichment media for stimulation was rich in urea, similar to Gomez et al.42 and Graddy et al.22, which is conditional stress for selective stimulation of ureolytic microorganisms. Moreover, the isolated strains were screened based on morphology and qualitative urease activity.Figure 3Neighbor-joining phylogenetic tree based on the 16S rRNA sequence of the isolates and reference sequence from the GenBank database (NCBI).Full size imageEvaluation of biocementation potential of the isolated strainsGrowth and pHThe various parameters of the biocementation potential of the isolated strains in comparison with Sporosarcina pasteurii (SP) have been plotted in Fig. 4. The growth characteristics of the isolates in NBU media and pH during growth have been represented in Fig. 4a and b. The initial pH of the growth media is kept at 7.5. It was observed that the pH of the growth media rises to 9.5 within 24 h of growth, indicating that these strains favor an alkaline environment to grow similar to SP43. All the isolates start growing when the pH of the media rises to 8.5 or above. Isolate LS2 was observed to have slower growth when compared with other isolates. This can be explained as LS2 belongs to different genera (Pseudogracilibacillus).Figure 4(a) Growth characteristics, (b) pH, (c) specific urease activity, (d) calcium utilization rate, and (e) carbonate precipitation rate of the isolates and consortia.Full size imageSpecific urease activityThe specific urease activities of the isolates were found to be comparable with SP (shown in Fig. 4c). Based on the provided NBU media and growth condition, the specific urease activity of SP is found to be 173.44 mM urea hydrolyzed h−1 (OD600)−1, which is around 2.9 mM urea hydrolyzed min−1 (OD600)−1. The specific urease activity of the isolate BS3 was observed to be maximum as 186.6 mM urea hydrolyzed h−1 (OD600)−1 during a growth period of 24 h and pH  > 9. Consortia also demonstrated significant urease activity as 160 mM urea hydrolyzed h−1 (OD600)−1 at a growth period of 48 h. The maximum ureolytic activity in BS1 was observed after 72 h of growth with a value of 106.67 mM urea hydrolyzed h−1 (OD600)−1. Maximum specific urease activity of the isolate BS2, BS4, and LS1 was observed to be 160.2, 120, and 173.4 mM urea hydrolyzed h−1 (OD600)−1 respectively after a growth duration of 48 h. LS2 demonstrated the maximum specific urease activity of 146.4 mM urea hydrolyzed h−1 (OD600)−1. The observed order of specific urease activity at 24 h of growth period is BS3  > SP  > Consortia  > LS1  > BS2  > BS4  > LS2  > BS1. As the urease activity of the strains depends on the growth media, urea content, and environmental conditions such as pH and Temperature44, we considered the conditions at the riverbank at the time of isolation, and the pH and temperature of the growth media were set at 7.5 and 37 degrees Celsius. The specific urease measured by the electrical conductivity method is reported to be between 3 to 9.7 mM urea hydrolyzed min−1 (OD600)−1 in yeast-extract urea media at pH 7 and temperature 30 degrees Celsius43. It is reported around 5 mM urea hydrolyzed min−1 (OD600)−1 in the nutrient broth urea (2%) media at a temperature of 25 degrees Celsius44. The comparative analysis of the urease activity (measured by electrical conductivity method) was done considering SP as positive control in this study. The maximum specific urease activities of all isolates were found to be in a range of 106.67 to 186.67 mM urea hydrolyzed h−1 (OD600)−1 (1.78 to 3.11 mM urea hydrolyzed min−1 OD600–1), which indicates that all of the isolated strains are capable of biocementation43,45.Calcium utilization and carbonate precipitation potentialIt was experimentally observed that the depletion of the supplemented soluble calcium in the precipitation media (PM) was corresponding to the ureolytic activities of the isolated strains. Within 48 h of introducing 1% bacteria (OD600 = 1) in the precipitation media, the soluble calcium chloride (50 mM) was utilized to precipitate carbonate crystals, as illustrated in Fig. 4d. Within 12 h of the inoculation period, BS3 was able to utilize 75% of the supplied calcium, while SP was able to utilize only 62.5% of the soluble calcium. The order of the calcium utilization potential in the isolates was observed as BS3 ≥ LS2  > L.S.1  > Consortia  > SP  > BS4  > BS2  > BS1 during the inoculation period. Contrarily, LS2, despite being a comparatively slow urease-producing bacteria, was able to utilize calcium ions at par with other isolates. Negligible changes were observed in the soluble calcium concentration of the control group eliminating the possibility of abiotic precipitation.The carbonate precipitation rate for each isolate (1% at OD600 = 1) for the 50 mM cementation media is plotted in Fig. 4e. The isolate BS3 with maximum ureolytic activity (specific urease activity 186.6 mM urea hydrolyzed min−1 OD600–1) precipitated the highest carbonate crystals after 96 h of the incubation period. BS3 precipitated 438 mg/100 ml of carbonate crystals, which is around 87.66% precipitation from the total supplied CaCl2, while precipitation with SP was quantified as 389 mg/100 ml (78%). The precipitation in consortia was observed to be 407 mg/100 ml (81%), which is slightly higher than SP. Precipitation in other isolates was found to be significantly lower than isolate BS3. Isolate BS1and BS2 precipitated 334 mg/100 ml (67%) and 343 mg/100 ml (69%) of carbonate crystals respectively, whereas isolate LS1 and LS2 precipitated around 357 mg/100 (71%) ml of carbonate crystals each. Isolate BS4 precipitated minimum carbonate crystals 292 mg/100 ml (58%). No precipitation was observed in the negative control set. Low concentrations of bacterial cells (1%) were considered in this experiment to slow down the urea hydrolysis in order to differentiate the calcium utilization potential of the isolated strains. This approach was modified from Dhami et al.46, and our results show agreement with their finding where 1% of SP cells depletes the 25 mM of CaCl2 within 24 h. It was observed that all the isolates took approximately 48 h to deplete the 50 mM CaCl2. The depletion of soluble calcium concentration was rapid in the initial 24 h in all the isolates. After 48 h, the residual soluble calcium was observed to be in the range of 2.5–5 mM in all the isolates (except BS1 and BS2), which might be due to loss of super-saturation caused by the unavailability of nutrient for bacterial cells to continue urea hydrolysis in the precipitation media13,43,47. The mineralogy of precipitated carbonate polymorphs (calcite, aragonite, vaterite) and the residual calcium are also influenced by pH, temperature, saturation index, dissolved organic carbon concentration, and the Ca2+ /CO32−ratio along with the presence of metabolites in the precipitation media13,47,48. As maximum precipitation was recovered with the isolate BS3, the isolate BS3 was selected for further investigation on soil improvement.Microstructure analysis of the precipitatesThe FESEM images of the carbonate crystal precipitated from BS3 were investigated further. The shape of the precipitated crystals was observed to be rhombohedral and trigonal (Fig. 5a). The average size of the crystals was observed in a range of 25 to 50 microns. The entrapped bacteria and rod-shaped bacterial imprints were identified (Fig. 5b), indicating that the bacteria acted as a nucleation site14. The smaller crystals were observed to coagulate in layers to develop larger calcite crystals. The entrapped bacteria were noticed on the grown and coagulated calcite crystal in Fig. 5c. After taking the FESEM image (Fig. 5a) of the precipitate, EDX analysis was conducted, and the elemental composition suggested an abundance of calcium, carbon, and oxygen, which indicates the presence of calcium carbonate crystals (Supplementary Fig. 1). XRD analysis was conducted to confirm the mineralogy of the precipitates, and the majority of the observed peaks of the XRD plot belonged to calcite, which is consistent with the observation of rhombohedral crystal shapes in the FESEM image. The XRD analysis also suggested an insignificant presence of aragonite in the precipitates.Figure 5FESEM images of the calcite precipitated from BS3 (a) Coagulated crystals (b) Bacterial imprints, (c) Entrapped bacteria on the precipitates.Full size imageApplication of native communities on riverbank soil and its influence on soil strengthNeedle penetration resistance of treated soilThe average NPI (N/mm) for different cases has been shown in Fig. 6a. No notable resistance was observed in the loose untreated sand (control) against the needle penetration. With one bio cementation cycle treatment, the consortia-treated soil sample (Consortia-BC1) demonstrated a higher value of NPI (5.15 N/mm) than SP-BC1 (4.19 N/mm) and BS3-BC1 (4.64 N/mm). The increase in the biocementation cycle treatment significantly improved the needle penetration resistance. Sample BS3-BC2 showed 116% improvement with the NPI value of 10.03 N/mm when compared to one cycle treated sample BS3-BC1. A similar trend was observed in the sample BS3-BC3 (NPI = 16.12 N/mm), which showed around 347% improvement in NPI when compared to sample BS3-BC1. From the needle penetration test, it was evident that the penetration resistance of treated soil improves significantly with the increased level of biocementation cycles, indirectly indicating an improvement in the soil erodibility resilience. Since non-uniformity is one of the undesired traits of MICP, a contour was plotted corresponding to the 25 points NPI, as shown in Fig. 6b. The contrasting color difference in the contours of the samples BS3-BC1, BS3-BC2, and BS-BC3 clearly demonstrates the stark difference in the strength of treated samples. The non-uniformity in the strength of treated soil crust of sample BS3-BC2 and BS3-BC3 can also be realized with the contrasting color gradient of the NPI contours.Figure 6Comparison of the Needle penetration resistance (N/mm) of the treated soil specimen (a) average values and (b) the contours.Full size imageSince the rate of penetration has an insignificant influence on the test results, the needle penetration test is recommended by the International Society of Rock Mechanics (ISRM) for quick, non-destructive testing of the strength of the stabilized soils and soft rocks49. As a large number of tests can be conducted due to the small diameter of the needle without destroying the sample, the needle penetration test is a better alternative to evaluate the local grain bonding in the biocemented soil than bulk strength properties like unconfined compressive strength and calcite content. Another rationale for choosing needle penetration test over conventional soil strength evaluation tests was that a pocket type penetrometer could be developed with the configuration in the present study for non-destructive monitoring of the soil strength improvement with biocementation application in the field. The response of the needle penetration resistance in terms of nominal strain (ratio of penetration to rod diameter) also indicated that the measured responses are independent of needle diameter for a small range, i.e., 1 to 3 mm49,50. A portable penetrometer of Maruto. Co. ltd. (needle maximum diameter 0.84 mm at 12 mm from the tip) have been correlated with high confidence value to conventional physicochemical parameters such as unconfined compressive strength (UCS), elasticity modulus, and elastic wave velocity in several studies50. In our setup, we have utilized a similar configuration chenille 22 needle with (maximum diameter 0.86 mm at 9 mm from the tip) and a penetration rate of 15 mm/minute for measuring the strength properties of cemented soil. Adopting the UCS and NPI correlation suggested by Ulusay et al.51, the UCS of samples BS3-BC1, BS3-BC2, and BS3-BC3 are around 1.67 MPa, 3.4 Mpa, and 5.3 Mpa.It is worth noting that in the needle penetration resistance tests, the boundary of the Petri dish can influence the test results. Therefore, trials were conducted and based on the findings, all penetrations were conducted at points at least five times the diameter of the needle away from the boundary to negate the influence of boundary conditions. The maximum penetration was conducted only up to 50% of the depth of prepared biocemented soil samples in the Petri dishes to avoid inference from the bottom of the Petri dish.Erodibility test in the hydraulic flumeTo investigate the influence of hydraulic current on different levels of biocementation, all the treated samples were exposed to hydraulic current gradually varying from gentle flow (0.06 m/s) to five times the critical velocity (0.75 m/s) in a 45-min duration test and soil mass loss percentage by the initial dry mass of the treated sample is presented as a measure of soil erodibility in Fig. 7. As expected, with an increase of biocementation cycles, i.e., calcite content, the soil erodibility reduced substantially. The initial dry weight of the samples control, BC1, BC2, and BC3, were measured as 398, 403, 406, and 410 g, respectively. Approximately 7.3% of calcite content resulted in a drastic reduction in erodibility (12% mass loss), while 56% soil mass loss was recorded for control (untreated sand). One biocementation cycle treatment (sample BS3-BC1) produced an average of 2.5% of calcite, reducing the soil loss to 31%. Sample BS3-BC2 with 4.93% calcite content resulted in 22% soil mass loss during the hydraulic flume test. It is worth noting that higher precipitation in the soil pores may hinder the flow of water in the soil matrix and increase the pore water pressure resulting in catastrophic failures. However, the MICP technique is reported to be a great tool to improve soil strength, maintaining an adequate hydraulic conductivity to prevent the build-up of the excess pore water pressure11. Theoretically, the percentage pore volume filled with precipitates for samples BS3-BC1, BS3-BC2, and BS3-BC3 considering the observed calcite contents and pore volume (100 ml) is around 3.7, 7.14, and 11.08%. The influence of pore water pressure on erodibility has not been established in the present study, and it certainly is one of the exciting parameters to consider for future studies.Figure 7Weight of the eroded soil (%) after the hydraulic flume test.Full size imageFrom the visual observation of the soil specimen after the flume test, it was evident that the soil particles start bonding with an increased level of bio cementation. A tough crust was formed on the top of BS3-BC2 and BS3-BC3 the samples, which got eroded with the fluvial current. Insignificant aggregation was observed in the sample BS3-BC1. However, with two and three cycles of biocementation treatment (BS3-BC2 and BS3-BC3), the biocemented soil particles (BCS) were evidently noticed (photos are shown in Supplementary Fig. 4).Clarà Saracho et al.27 addressed the erosion due to tangential flow (similar to river current) by treating the soil specimen with ten pore volume of low concentration of cementation media (0.02 M to 0.1 M) by injection strategy and tested the specimen in the flow velocity ranging from 0.035 to 0.185 m/sec for 120 min in a modified erosion function apparatus (EFA). The study concluded that the treatment with 0.08 M cementation media (calcite content varying from 1.2 to 4%) resulted in negligible erosion in the stated test conditions, and with the increase in MICP treatment, a shift in the erosion mode from particulate mode to block failure was observed indicating that with the increase in calcite content and needle penetration resistance, there might be a threshold for biocemented soil, where the soil erosion might be catastrophic due to block failure. However, in this study, a consistent decrease in soil erodibility is observed with the increase in needle penetration resistance. We found that 7.3% of calcite content was required to control the soil erodibility substantially in the test flow range (0.06 m/s to 0.75 m/s). A similar trend was observed by Kou et al.52 and Chung et al.53, where consistent reduction in wave-induced erosion and rainfall-induced erosion was observed with an increase in needle penetration resistance for biocemented fine sand treated with the exogenous bacteria.Another aspect to note in the context of the applied treatment is the produced ammonia which can be toxic for riparian flora and fauna15. From stochiometric calculations, for each biocementation cycle, the produced ammonia is evaluated as 8.5 kg per metric ton of soil treated, and for the best performing treatment approach, i.e., BS3-BC3, the ammonia generated is evaluated as 2.63% by weight of the retained soil. The acceptable limit of ammonia in the surface water was recommended as 17 mg per liter for acute exposure and 1.9 mg/l for chronic exposure for protecting the aquatic life in the freshwater as per the environmental protection agency54. With the MICP technique, the threat of produced ammonia crossing the maximum acceptable quantity is highly plausible; however, for the field application, the ammonia generated can be reduced by reducing the quantity of reagents and increasing the period of applications. It is to be noted that the produced ammonia will also be subjected to dilution in the river stream. The average discharge of the Brahmaputra river is around 19.8 megaliters per second in the Assam valley33. Therefore, the area of the riverbank to be biocemented must be decided judiciously with the context of the produced ammonia quantity and its possible dilution to non-toxic levels.Microstructure and mineralogical analysis of the biocemented samplesTo investigate the influence of different biocementation levels on the erodibility of the treated sand grains, FESEM imaging was conducted for light biocemented samples (BS3-BC1) and heavy biocemented samples (BS3-BC3). While precipitated crystals were observed to be growing on the grooves of sand grains in the light biocemented sample (BS3-BC1), bridging of sand grains with rhombohedral crystals was observed in the heavy biocemented sample (BS3-BC3), as shown in Fig. 8. The effective calcium carbonate bridging between sand grains increases the frictional and cohesive property of sand grains55,56, leading to a substantial reduction in the erodibility of the soil. Bacterial imprints were observed in both cases, suggesting the precipitation to be biodgenic14. Further EDX analysis on a bridged sand grain (Supplementary Fig. 5) suggested an abundance of silicon and oxygen on the sand grains with a trace amount of chlorine and calcium. This indicates the presence of residual calcium chloride on silica grains. The EDX analysis on the grain bridge indicated the presence of calcium, carbon, and oxygen, suggesting CaCO3 precipitation. XRD analysis on treated and untreated sand confirmed the precipitation of calcite. Most of the peaks correspond to quartz (silica). In the biocemented sand sample, a visible peak of calcite was observed at around 29 degrees of 2Ɵ (Details in Supplementary Fig. 6). Therefore, the incorporation of microbial calcite as a binding agent for loose grain silica soil was found to reduce the soil erodibility.Figure 8FESEM images of the treated sand grains (a). Calcite crystal growing on the sand grains on BS3- BC1 sample (b). Calcite bridging in BS3-BC3 samples.Full size image More

1.Pagano, M. C. & Miransari, M. The importance of soybean production worldwide. In Abiotic and Biotic Stresses in Soybean Production Vol. 1 (ed. Miransari, M.) 1–26 (Academic Press, 2016).
Google Scholar 
2.FAOSTAT. Food and Agriculture Organisation Statistical Database http://www.apps.fao.org/faostat. Accessed 23 May 2021.3.MacDonald, R. S. et al. Environmental influences on isoflavones and saponins in soybeans and their role in colon cancer. J. Nutr. 135, 1239–1242 (2005).CAS 
PubMed 
Article 

Google Scholar 
4.Tidke, S. A. et al. Assessment of anticancer, anti-inflammatory and antioxidant properties of isoflavones present in soybean. Res. J. Phytochem. 12, 35–42 (2018).ADS 
CAS 
Article 

Google Scholar 
5.Hill, J. H. & Whitham, S. A. Control of virus diseases in soybeans. Adv. Virus Res. 90, 355–390 (2014).PubMed 
Article 

Google Scholar 
6.Tian, B. et al. Host adaptation of soybean dwarf virus following serial passages on pea (Pisum sativum) and soybean (Glycine max). Viruses 9, 155 (2017).PubMed Central 
Article 
CAS 
PubMed 

Google Scholar 
7.Wang, R. Y., Kritzman, A., Hershman, D. E. & Ghabrial, S. A. Aphis glycines as a vector of persistently and nonpersistently transmitted viruses and potential risks for soybean and other crops. Plant Dis. 90, 920–926 (2006).CAS 
PubMed 
Article 

Google Scholar 
8.Hesler, L. S., Dashiell, E., Jonathan, A. E. & Lundgren, G. Characterization of resistance to Aphis glycines in soybean accessions. Euphytica 154, 91–99 (2007).Article 

Google Scholar 
9.Baldin, E. L. L. et al. Feeding behavior of Aphis glycines (Hemiptera: Aphididae) on soybeans exhibiting antibiosis, antixenosis, and tolerance resistance. Fla. Entomol. 101, 223–228 (2018).Article 

Google Scholar 
10.Chang, H.-X. & Hartman, G. L. Characterization of insect resistance loci in the USDA soybean germplasm collection using genome-wide association studies. Front. Plant Sci. 8, 670. https://doi.org/10.3389/fpls.2017.00670 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
11.Bansal, R., Mian, M. A. R. & Michel, A. Characterizing resistance to soybean aphid (Hemiptera: Aphididae): Antibiosis and antixenosis assessment. J. Econ. Entomol. https://doi.org/10.1093/jee/toab038 (2021).Article 
PubMed 

Google Scholar 
12.Klein, A. T. et al. Investigation of the chemical interface in the soybean−aphid and rice−bacteria interactions using MALDI-Mass Spectrometry Imaging. Anal. Chem. 87, 5294–5301 (2015).CAS 
PubMed 
Article 

Google Scholar 
13.Hohenstein, J. D. et al. Transcriptional and chemical changes in soybean leaves in response to long-term aphid colonization. Front. Plant Sci. 10, 310 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
14.Blackman, R. L. & Eastop, V. F. Taxonomic issues. In Aphids as Crop Pests (eds van Emden, H. F. & Harrington, R.) 1–36 (CABI, 2017).
Google Scholar 
15.Wale, M., Jembere, B. & Seyoum, E. Occurrence of the pea aphid, Acyrthosiphon pisum (Harris) (Homoptera: Aphididae) on wild leguminous plants in West Gojam, Ethiopia, Sinet. Ethiopian J. Sci. 26, 83–87 (2003).
Google Scholar 
16.Chan, C. K., Forbes, A. R. & Raworth, D. A. Aphid-transmitted viruses and their vectors of the world. Agric. Can. Tech. Bull. 3E, 1–216 (1991).
Google Scholar 
17.Rashed, A. et al. Vector-borne viruses of pulse crops, with a particular emphasis on North American cropping system. Ann. Entomol. Soc. Am. 111, 205–227 (2018).CAS 
Article 

Google Scholar 
18.Stavrinides, J., McCloskey, J. K. & Ochman, H. Pea Aphid as both host and vector for the phytopathogenic bacterium Pseudomonas syringae. Appl. Environ. Microbiol. 75, 2230–2235 (2009).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
19.Peccoud, J., Ollivier, A., Plantegenest, M. & Simon, J.-C. A continuum of genetic divergence from sympatric host races to species in the pea aphid complex. PNAS 106, 7495–7500 (2009).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
20.Caillaud, M. C. & Via, S. Specialized feeding behavior influences both ecological specialization and assortative mating in sympatric host races of pea aphids. Am. Nat. 156, 606–621 (2000).PubMed 
Article 

Google Scholar 
21.Ferrari, J., Godfray, H. C., Faulconbridge, A. S., Prior, K. & Via, S. Population differentiation and genetic variation in host choice among pea aphids from eight host plant genera. Evolution 60, 1574–1584 (2006).PubMed 
Article 

Google Scholar 
22.Mitku, G. & Damte, T. Development, reproduction, and host preference of Acyrthosiphon pisum (Harris) (Homoptera: Aphididae) on selected lentil genotypes and resistance index of these selected lentil genotypes to pea aphid. Int. J. Entomol. Res. 4, 16–22 (2019).
Google Scholar 
23.Powell, G., Tosh, C. R. & Hardie, J. Host plant selection by aphids: behavioral, evolutionary, and applied perspectives. Annu. Rev. Entomol. 51, 309–330 (2006).CAS 
PubMed 
Article 

Google Scholar 
24.Kordan, B. et al. European yellow lupine Lupinus luteus and narrow-leaf lupine Lupinus angustifolius as hosts for the pea aphid Acyrthosiphon pisum. Entomol. Exp. Appl. 128, 139–146 (2008).Article 

Google Scholar 
25.Kordan, B. et al. Susceptibility of forage legumes to infestation by the pea aphid Acyrthosiphon pisum (Harris) (Hemiptera: Aphididae). Crop Pasture Sci. 69, 775–784 (2018).Article 

Google Scholar 
26.Kordan, B. et al. Antixenosis potential in pulses against the pea aphid (Hemiptera: Aphididae). J. Econ. Entomol. 112, 465–474 (2019).PubMed 
Article 

Google Scholar 
27.Pettersson, J., Tjallingii, W. F. & Hardie, J. Host-plant selection and feeding. In Aphids as Crop Pests (eds van Emden, H. F. & Harrington, R.) 173–195 (CABI, 2017).Chapter 

Google Scholar 
28.Martin, B., Collar, J. L., Tjallingi, W. F. & Fereres, A. Intracellular ingestion and salivation by aphids may cause the acquisition and inoculation of non-persistently transmitted plant viruses. J. Gen. Virol. 78, 2701–2705 (1997).CAS 
PubMed 
Article 

Google Scholar 
29.Garzo, E., Moreno, A., Plaza, M. & Fereres, A. Feeding behavior and virus-transmission ability of insect vectors exposed to systemic insecticides. Plants 9, 895. https://doi.org/10.3390/plants9070895 (2020).CAS 
Article 
PubMed Central 
PubMed 

Google Scholar 
30.Onstad, D. W. & Knolhoff, L. Arthropod resistance to crops. In Insect Resistance Management (ed. Onstad, D. W.) 293–326 (Academic Press, 2014).Chapter 

Google Scholar 
31.Smith, C. M. & Clement, S. L. Molecular bases of plant resistance to arthropods. Annu. Rev. Entomol. 57, 309–328 (2012).CAS 
PubMed 
Article 

Google Scholar 
32.Stout, M. J. Reevaluating the conceptual framework for applied research on host-plant resistance. Insect Sci. 20, 263–272 (2013).PubMed 
Article 

Google Scholar 
33.Smith, C. M. & Chuang, W. P. Plant resistance to aphid feeding: behavioral, physiological, genetic and molecular cues regulate aphid host selection and feeding. Pest Manag. Sci. 70, 528–540 (2014).PubMed 
Article 
CAS 

Google Scholar 
34.Dogimont, C., Bendahmane, A., Chovelon, V. & Boissot, N. Host plant resistance to aphids in cultivated crops: Genetic and molecular bases, and interactions with aphid populations. C. R. Biol 333, 566–573 (2010).CAS 
PubMed 
Article 

Google Scholar 
35.Chandran, P. et al. Feeding behavior comparison of soybean aphid (Hemiptera: Aphididae) biotypes on different soybean genotypes. J. Econ. Entomol. 106, 2234–2240 (2013).PubMed 
Article 

Google Scholar 
36.Simmonds, M. S. J. Flavonoid-insect interactions: Recent advances in our knowledge. Phytochemistry 64, 21–30 (2003).CAS 
Article 

Google Scholar 
37.Mai, V. C. et al. Differential induction of Pisum sativum defense signaling molecules in response to pea aphid infestation. Plant Sci. 221–222, 1–12 (2014).PubMed 
Article 
CAS 

Google Scholar 
38.Morkunas, I. et al. Pea aphid infestation induces changes in flavonoids, antioxidative defence, soluble sugars and sugar transporter expression in leaves of pea seedlings. Protoplasma 253, 1063–1079 (2016).CAS 
PubMed 
Article 

Google Scholar 
39.Woźniak, A. et al. The dynamics of the defense strategy of pea induced by exogenous nitric oxide in response to aphid infestation. Int. J. Mol. Sci. 18, 329. https://doi.org/10.3390/ijms18020329 (2017).CAS 
Article 
PubMed Central 

Google Scholar 
40.Buer, C. S., Muday, G. K. & Djordjevic, M. A. Implications of long-distance flavonoid movement in Arabidopsis thaliana. Plant Signal. Behav. 3, 415–417 (2008).PubMed 
PubMed Central 
Article 

Google Scholar 
41.Petrussa, E. et al. Plant flavonoids: Biosynthesis, transport and involvement in stress responses. Int. J. Mol. Sci. 14, 14950–14973 (2013).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
42.Zhao, J. Flavonoid transport mechanisms: How to go, and with whom. Trends Plant Sci. 20, 576–585 (2015).CAS 
PubMed 
Article 

Google Scholar 
43.Alseekh, S., de Souza, L. P., Benina, M. & Fernie, A. L. The style and substance of plant flavonoid decoration; Towards defining both structure and function. Phytochemistry 174, 112347. https://doi.org/10.1016/j.phytochem.2020.112347 (2020).CAS 
Article 
PubMed 

Google Scholar 
44.Klingauf, F. A. Host plant finding and acceptance. In Aphids, Their Biology, Natural Enemies and Control Vol. 2 (eds Minks, A. K. & Harrewijn, P.) 209–223 (Elsevier, 1987).
Google Scholar 
45.Tjallingii, W. F. & Mayoral, A. M. Criteria for host plant acceptance by aphids. In Proceeding 8th International Symposium Insect–Plant Relationships (eds Menken, S. B. J. et al.) 280–282 (Kluwer Academic Publishers, 1992).Chapter 

Google Scholar 
46.Wensler, R. J. & Filshie, B. K. Gustatory sense organs in the food canal of aphids. J. Morph. 129, 473–492 (1969).Article 

Google Scholar 
47.Gabryś, B. & Tjallingii, W. F. The role of sinigrin in host plant recognition by aphids during initial plant penetration. Entomol. Exp. Appl. 104, 89–93 (2002).Article 

Google Scholar 
48.Philippi, J. et al. Correlation of the alkaloid content and composition of narrow-leafed lupins (Lupinus angustifolius L.) to aphid susceptibility. J. Pest Sci. 89, 359–373 (2016).Article 

Google Scholar 
49.Van Hoof, H. A. An investigation of the biological transmission of a non-persistent virus. Doctoral thesis (Van Putten and Oortmijer, 1958).50.Dancewicz, K., Szumny, A., Wawrzeńczyk, C. & Gabryś, B. Repellent and antifeedant activities of citral-derived lactones against the peach potato aphid. Int. J. Mol. Sci. 21, 8029. https://doi.org/10.3390/ijms21218029 (2020).CAS 
Article 
PubMed Central 
PubMed 

Google Scholar 
51.Kordan, B. et al. Variation in susceptibility of rapeseed cultivars to the peach potato aphid. J. Pest. Sci. 94, 435–449 (2021).Article 

Google Scholar 
52.Pritchard, J. & Vickers, L. H. Aphids and stress. In Aphids as Crop Pests (eds Van Emden, H. F. & Harrington, R.) 132–147 (CABI, 2017).Chapter 

Google Scholar 
53.Pompon, J. & Pelletier, Y. Changes in aphid probing behaviour as a function of insect age and plant resistance level. Bull. Entomol. Res. 102, 550–557 (2012).CAS 
PubMed 
Article 

Google Scholar 
54.van Emden, H. F. Host-plant resistance. In Aphids as Crop Pests (eds van Emden, H. F. & Harrington, R.) 515–532 (CABI, 2017).Chapter 

Google Scholar 
55.Gould, K. S. & Lister, C. Flavonoid functions in plants. In Flavonoids, Chemistry, Biochemistry and Applications (eds Andersen, Ø. M. & Markham, K. R.) 397–442 (CRC Press, 2006).
Google Scholar 
56.Goławska, S., Kapusta, I., Łukasik, I. & Wójcicka, A. Effect of phenolics on the pea aphid, Acyrthosiphon pisum (Harris) population on Pisum sativum L. (Fabaceae). Pestycydy. 3–4, 71–77 (2008).
Google Scholar 
57.Goławska, S. & Łukasik, I. Antifeedant activity of luteolin and genistein against the pea aphid, Acyrthosiphon pisum. J. Pest Sci. 85, 443–450 (2012).Article 

Google Scholar 
58.Goławska, S. et al. Alfalfa (Medicago sativa L.) apigenin glycosides and their effect on the pea aphid (Acyrthosiphon pisum). Polish J. Environ. Stud. 19, 913–919 (2010).
Google Scholar 
59.Johnson, A. D. & Singh, A. Larvicidal activity and biochemical effects of apigenin against filarial vector Culex quinquefasciatus. Int. J. Life. Sci. Sci. Res. 3, 1315–1321 (2017).
Google Scholar 
60.Boué, S. M. & Raina, A. K. Effects of plant flavonoids on fecundity, survival, and feeding of the formosan subterranean termite. J. Chem. Ecol. 29, 2575–2584 (2003).PubMed 
Article 

Google Scholar 
61.Xu, D. et al. Antifeedant activities of secondary metabolites from Ajuga nipponensis against adult of striped flea beetles, Phyllotreta striolata. J. Pest Sci. 82, 195–202 (2009).Article 

Google Scholar 
62.Goławska, S., Sprawka, I. & Łukasik, I. Effect of saponins and apigenin mixtures on feeding behavior of the pea aphid, Acyrthosiphon pisum Harris. Biochem. Syst. Ecol. 55, 137–144 (2014).Article 
CAS 

Google Scholar 
63.Zavala, J. A., Scopel, A. L. & Ballaré, C. L. Effects of ambient UV-B radiation on soybean crops: Impact on leaf herbivory by Anticarsia gemmatalis. Plant Ecol. 156, 121–130 (2001).Article 

Google Scholar 
64.Bentivenha, J. P. F. et al. Role of the rutin and genistein flavonoids in soybean resistance to Piezodorus guildinii (Hemiptera: Pentatomidae). Arthropod Plant Interact. 12, 311–320 (2018).Article 

Google Scholar 
65.Hoffmann-Campo, C. B., Harborne, J. B. & McCaffery, A. R. Pre-ingestive and post-ingestive effects of soya bean extracts and rutin on Trichoplusia ni growth. Entomol. Exp. Appl. 98, 181–194 (2001).Article 

Google Scholar 
66.Yuan, E. et al. Increases in genistein in Medicago sativa confer resistance against the Pisum host race of Acyrthosiphon pisum. Insects. 10, 97. https://doi.org/10.3390/insects10040097 (2019).Article 
PubMed Central 
PubMed 

Google Scholar 
67.Meng, F. et al. QTL underlying the resistance to soybean aphid (Aphis glycines Matsumura) through isoflavone-mediated antibiosis in soybean cultivar ‘Zhongdou 27’. Theor Appl. Genet. 123, 1459–1465 (2011).CAS 
PubMed 
Article 

Google Scholar 
68.Murakami, S. et al. Insect-induced daidzein, formononetin and their conjugates in soybean leaves. Metabolites 4, 532–546 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
69.Lattanzio, V. et al. Role of endogenous flavonoids in resistance mechanism of Vigna to aphids. J. Agric. Food Chem. 48, 5316–5320 (2000).CAS 
PubMed 
Article 

Google Scholar 
70.Delgado-Núñez, E. J. et al. Isorhamnetin: A nematocidal flavonoid from Prosopis laevigata leaves against Haemonchus contortus eggs and larvae. Biomolecules 10, 773 (2020).PubMed Central 
Article 
CAS 
PubMed 

Google Scholar 
71.Gómez, J. D., Vital, C. E., Oliveira, M. G. A. & Ramos, H. J. O. Broad range flavonoid profiling by LC/MS of soybean genotypes contrasting for resistance to Anticarsia gemmatalis (Lepidoptera: Noctuidae). PLoS ONE https://doi.org/10.1371/journal.pone.0205010 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
72.Khan, M. A. M., Ulrichs, C. & Mewis, I. Effect of water stress and aphid herbivory on flavonoids in broccoli (Brassica oleracea var. italica Plenck). J. Appl. Bot. Food Qual. 84, 178–182 (2011).CAS 

Google Scholar 
73.Bale, J. S., Ponder, K. L. & Pritchard, J. Coping with stress. In Aphids as Crop Pests (eds van Emden, H. F. & Harrington, R.) 287–309 (CABI, 2007).Chapter 

Google Scholar 
74.Atteyat, M., Abu-Romman, S., Abu-Darwish, M. & Ghabeish, I. Impact of flavonoids against woolly apple aphid, Eriosoma lanigerum (Hausmann) and its sole parasitoid, Aphelinus mali (Hald). J. Agric. Sci. 4, 227–236 (2012).
Google Scholar 
75.Tjallingii, W. F. & Hogen Esch, T. Fine structure of aphid stylet routes in plant tissues in correlation with EPG signals. Physiol. Entomol. 18, 317–328 (1993).Article 

Google Scholar 
76.Cherqui, A. & Tjallingii, W. F. Salivary proteins of aphids, a pilot study on identification, separation and immunolocalisation. J. Insect Physiol. 46, 1177–1186 (2000).CAS 
PubMed 
Article 

Google Scholar 
77.Silva-Sanzana, C., Estevez, J. M. & Blanco-Herrera, F. Influence of cell wall polymers and their modifying enzymes during plant–aphid interactions. J. Exp. Bot. 71, 3854–3864 (2020).CAS 
PubMed 
Article 

Google Scholar 
78.Alvarez, A. E. et al. Location of resistance factors in the leaves of potato and wild tuber-bearing Solanum species to aphid Myzus persicae. Entomol. Exp. Appl. 121, 145–157 (2006).Article 

Google Scholar 
79.Alvarez, A. E. et al. Infection of potato plants with potato leafroll virus changes attraction and feeding behaviour of Myzus persicae. Entomol. Exp. Appl. 125, 135–144 (2007).Article 

Google Scholar 
80.Machado-Assefh, C. R. & Alvarez, A. E. Probing behavior of aposymbiotic green peach aphid (Myzus persicae) on susceptible Solanum tuberosum and resistant Solanum stoloniferum plants. Insect Sci. 25, 127–136 (2018).CAS 
PubMed 
Article 

Google Scholar 
81.Common Catalogue of Varieties of Agricultural Plant Species [CCA]. 37th complete edition. Official Journal of the European Union C 13/1 (2019). Accessed 23 May 2021.82.Porejestrowe doświadczalnictwo odmianowe. Charakterystyka odmian. http://www.coboru.gov.pl/Polska/Rejestr/odm_w_rej.aspx?kodgatunku=SOS. Accessed 23 May 2021.83.Meier, U. Growth stages of mono- and dicotyledonous plants: BBCH. Monograph (Julius Kühn-Institut, 2018).84.Beer, K., Joschinski, J., Sastre, A. A., Kraus, J. & Helfrich-Forster, C. A damping circadian clock drives weak oscillations in metabolism and locomotor activity of aphids (Acyrthosiphon pisum). Sci. Rep. 7, 1–5. https://doi.org/10.1038/s41598-017-15014-3 (2017).CAS 
Article 

Google Scholar 
85.Joschinski, J., Beer, K., Helfrich-Forster, C. & Krauss, J. Pea aphids (Hemiptera: Aphididae) have diurnal rhythms when raised independently of a host plant. J. Insect. Sci. 16, 1–5 (2016).Article 

Google Scholar 
86.Graham, T. L. Flavonoid and isoflavonoid distribution in developing soybean seedling tissues and in seed and root exudates. Plant Physiol. 95, 594–603 (1991).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
87.Lee, J. H. et al. Characterization of isoflavones accumulation in developing leaves of soybean (Glycine max) cultivars. J. Korean Soc. Appl. Biol. Chem. 52(2), 139–143 (2009).CAS 
Article 

Google Scholar 
88.Magarelli, G. et al. Rutin and total isoflavone determination in soybean at different growth stages by using voltammetric methods. Microchem. J. 117, 149–155 (2014).CAS 
Article 

Google Scholar 
89.Perlatti, B. et al. Application of a quantitative HPLC-ESI-MS/MS method for flavonoids in different vegetables matrices. J. Braz. Chem. Soc. 27(3), 475–483 (2016).CAS 

Google Scholar 
90.Biesaga, M. & Pyrzyńska, K. Stability of bioactive polyphenols from honey during different extraction methods. Food Chem. 136, 46–54 (2013).CAS 
PubMed 
Article 

Google Scholar 
91.Sergiel, I., Pohl, P. & Biesaga, M. Characterisation of honeys according to their content of phenolic compounds using high performance liquid chromatography/tandem mass spectrometry. Food Chem. 145, 404–408 (2014).CAS 
PubMed 
Article 

Google Scholar  More

1.Arbour, J. H. & Santana, S. E. A major shift in diversification rate helps explain macroevolutionary patterns in primate species diversity. Evolution 71, 1600–1613 (2017).PubMed 
Article 

Google Scholar 
2.Groves, C. Primates (Taxonomy) in The International Encyclopedia of Primatology (ed Augustin Fuentes) (John Wiley & Sons, Inc., 2016).3.Cotton, A., Clark, F., Boubli, J. & Schwitzer, C. IUCN red list of threatened primate species in An Introduction to Primate Conservation 31–18 (Oxford University Press, 2016).
Google Scholar 
4.Stumpf, R. M. et al. Microbiomes, metagenomics, and primate conservation: New strategies, tools, and applications. Biol. Conserv. 199, 56–66 (2016).Article 

Google Scholar 
5.West, A. G. et al. The microbiome in threatened species conservation. Biol. Conserv. 229, 85–98 (2019).Article 

Google Scholar 
6.Cunningham, A. A., Daszak, P. & Wood, J. L. N. One Health, emerging infectious diseases and wildlife: two decades of progress?. Philos. Trans. R. Soc. B: Biol. Sci. 372, 20160167 (2017).Article 

Google Scholar 
7.Ramey, A. M. & Ahlstrom, C. A. Antibiotic resistant bacteria in wildlife: Perspectives on trends, acquisition and dissemination, data gaps, and future directions. J. Wildl. Dis. 56, 1–15 (2020).PubMed 
Article 

Google Scholar 
8.Clayton, J. B. et al. Captivity humanizes the primate microbiome. Proc. Natl. Acad. Sci. 113, 10376–10381 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
9.Hale, V. L. et al. Gut microbiota in wild and captive Guizhou snub-nosed monkeys. Rhinopithecus brelichi. Am. J. Primatol. 81, e22989 (2019).CAS 
PubMed 

Google Scholar 
10.Kriss, M., Hazleton, K. Z., Nusbacher, N. M., Martin, C. G. & Lozupone, C. A. Low diversity gut microbiota dysbiosis: drivers, functional implications and recovery. Curr. Opin. Microbiol. 44, 34–40 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
11.Mahnert, A. et al. Man-made microbial resistances in built environments. Nat. Commun. 10, 1–12 (2019).CAS 
Article 

Google Scholar 
12.Amato, K. R. et al. Using the gut microbiota as a novel tool for examining colobine primate GI health. Glob. Ecol. Conserv. 7, 225–237 (2016).Article 

Google Scholar 
13.Zhu, H. et al. Diarrhea-associated intestinal microbiota in captive Sichuan golden snub-nosed monkeys (Rhinopithecus roxellana). Microbes Environ. ME17163 (2018).14.Campbell, T. P. et al. The microbiome and resistome of chimpanzees, gorillas, and humans across host lifestyle and geography. ISME J. 14, 1584–1599 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
15.Buzzard, P. J. Ecological partitioning of Cercopithecus campbelli, C. petaurista, and C. diana in the Taï Forest. Int. J. Primatol. 27, 529–558 (2006).Article 

Google Scholar 
16.Chapman, C. A. et al. The guenons: diversity and adaptation in African monkeys. 325–350 (Springer, 2004).17.Krishnadas, M., Chandrasekhara, K. & Kumar, A. The response of the frugivorous lion-tailed macaque (Macaca silenus) to a period of fruit scarcity. Am. J. Primatol. 73, 1250–1260 (2011).PubMed 
Article 
PubMed Central 

Google Scholar 
18.Swedell, L., Hailemeskel, G. & Schreier, A. Composition and seasonality of diet in wild hamadryas baboons: preliminary findings from Filoha. Folia Primatol. 79, 476–490 (2008).Article 

Google Scholar 
19.Basabose, A. K. Diet composition of chimpanzees inhabiting the montane forest of Kahuzi, Democratic Republic of Congo. Am. J. Primatol. 58, 1–21 (2002).PubMed 
Article 
PubMed Central 

Google Scholar 
20.McLennan, M. R. & Ganzhorn, J. U. Nutritional characteristics of wild and cultivated foods for chimpanzees (Pan troglodytes) in agricultural landscapes. Int. J. Primatol. 38, 122–150 (2017).Article 

Google Scholar 
21.Newton-Fisher, N. E. The diet of chimpanzees in the Budongo Forest Reserve Uganda. Afr. J. Ecol. 37, 344–354 (1999).Article 

Google Scholar 
22.Bach, T. H., Chen, J., Hoang, M. D., Beng, K. C. & Nguyen, V. T. Feeding behavior and activity budget of the southern yellow-cheeked crested gibbons (Nomascus gabriellae) in a lowland tropical forest. Am. J. Primatol. 79, e22667 (2017).Article 

Google Scholar 
23.Fan, P.-F., Fei, H.-L., Scott, M. B., Zhang, W. & Ma, C.-Y. Habitat and food choice of the critically endangered cao vit gibbon (Nomascus nasutus) in China: implications for conservation. Biol. Conserv. 144, 2247–2254 (2011).Article 

Google Scholar 
24.Fan, P. F., Fei, H. L. & Ma, C. Y. Behavioral responses of cao vit gibbon (Nomascus nasutus) to variations in food abundance and temperature in Bangliang, Jingxi China. Am. J. Primatol. 74, 632–641 (2012).PubMed 
Article 
PubMed Central 

Google Scholar 
25.McConkey, K. R., Ario, A., Aldy, F. & Chivers, D. J. Influence of forest seasonality on gibbon food choice in the rain forests of Barito Ulu Central Kalimantan. Int. J. Primatol. 24, 19–32 (2003).Article 

Google Scholar 
26.Amora, T. D., BeltrÃO-Mendes, R. & Ferrari, S. F. Use of alternative plant resources by common marmosets (Callithrix jacchus) in the semi-arid Caatinga scrub forests of northeastern Brazil. Am. J. Primatol. 75, 333–341 (2013).PubMed 
Article 
PubMed Central 

Google Scholar 
27.Dietz, J. M., Peres, C. A. & Pinder, L. Foraging ecology and use of space in wild golden lion tamarins (Leontopithecus rosalia). Am. J. Primatol. 41, 289–305 (1997).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
28.Garber, P. A. Feeding ecology and behaviour of the genus Saguinus. Marmosets and tamarins: systematics behaviour and ecology (1993).29.Heymann, E. W., Knogge, C. & Tirado Herrera, E. R. Vertebrate predation by sympatric tamarins, Saguinus mystax and Saguinus fuscicollis. Am. J. Primatol. 51, 153–158 (2000).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
30.Porter, L. M. Dietary differences among sympatric Callitrichinae in northern Bolivia: Callimico goeldii, Saguinus fuscicollis and S. labiatus. Int. J. Primatol. 22, 961–992 (2001).Article 

Google Scholar 
31.Anapol, F. & Lee, S. Morphological adaptation to diet in platyrrhine primates. Am. J. Phys. Anthropol. 94, 239–261 (1994).CAS 
PubMed 
Article 

Google Scholar 
32.Nash, L. T. Dietary, behavioral, and morphological aspects of gummivory in primates. Am. J. Phys. Anthropol. 29, 113–137 (1986).Article 

Google Scholar 
33.Abreu, F., De la Fuente, M. F. C., Schiel, N. & Souto, A. Feeding ecology and behavioral adjustments: flexibility of a small neotropical primate (Callithrix jacchus) to survive in a semiarid environment. Mammal Res. 61, 221–229 (2016).Article 

Google Scholar 
34.Cunha, A. A., Vieira, M. V. & Grelle, C. E. V. Preliminary observations on habitat, support use and diet in two non-native primates in an urban Atlantic forest fragment: the capuchin monkey (Cebus sp.) and the common marmoset (Callithrix jacchus) in the Tijuca forest Rio de Janeiro. Urban Ecosyst. 9, 351–359 (2006).Article 

Google Scholar 
35.Passamani, M. & Rylands, A. B. Feeding behavior of Geoffroy’s marmoset (Callithrix geoffroyi) in an Atlantic forest fragment of south-eastern Brazil. Primates 41, 27–38 (2000).PubMed 
Article 

Google Scholar 
36.Veracini, C. Habitat use and ranging behavior of the silvery marmoset (Mico argentatus) at Caxiuanã National Forest (eastern Brazilian Amazonia) in The smallest anthropoids 221–240 (Springer, 2009).37.Yépez, P., De La Torre, S. & Snowdon, C. T. Interpopulation differences in exudate feeding of pygmy marmosets in Ecuadorian Amazonia. Am. J. Primatol. 66, 145–158 (2005).PubMed 
Article 

Google Scholar 
38.Hale, V. L. et al. Diet versus phylogeny: a comparison of gut microbiota in captive colobine monkey species. Microb. Ecol. 75, 515–527 (2018).PubMed 
Article 

Google Scholar 
39.Amato, K. R. et al. The gut microbiota appears to compensate for seasonal diet variation in the wild black howler monkey (Alouatta pigra). Microb. Ecol. 69, 434–443 (2015).CAS 
PubMed 
Article 

Google Scholar 
40.Frankel, J. S., Mallott, E. K., Hopper, L. M., Ross, S. R. & Amato, K. R. The effect of captivity on the primate gut microbiome varies with host dietary niche. Am. J. Primatol. 81, e23061 (2019).PubMed 
Article 

Google Scholar 
41.McKenzie, V. J. et al. The effects of captivity on the mammalian gut microbiome. Integr. Comp. Biol. 57, 690–704 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
42.Lugli, G. A. et al. Evolutionary development and co‐phylogeny of primate‐associated bifidobacteria. Environ. Microbiol. (2020).43.Milani, C. et al. Unveiling bifidobacterial biogeography across the mammalian branch of the tree of life. ISME J. 11, 2834–2847 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
44.Lugli, G. A. et al. Comparative genomic and phylogenomic analyses of the Bifidobacteriaceae family. BMC Genom. 18, 568 (2017).Article 
CAS 

Google Scholar 
45.Pokusaeva, K., Fitzgerald, G. F. & van Sinderen, D. Carbohydrate metabolism in Bifidobacteria. Genes Nutr. 6, 285–306 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
46.Stewart, C. J. et al. Temporal development of the gut microbiome in early childhood from the TEDDY study. Nature 562, 583–588 (2018).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
47.Orkin, J. D. et al. Seasonality of the gut microbiota of free-ranging white-faced capuchins in a tropical dry forest. ISME J. 13, 183–196 (2019).CAS 
PubMed 
Article 

Google Scholar 
48.Neuzil-Bunesova, V. et al. Five novel bifidobacterial species isolated from faeces of primates in two Czech zoos: Bifidobacterium erythrocebi sp. nov., Bifidobacterium moraviense sp. nov., Bifidobacterium oedipodis sp. nov., Bifidobacterium olomucense sp. nov. and Bifidobacterium panos sp. nov. Int. J. Syst. Evol. Microbiol. (2020).49.Duranti, S. et al. Characterization of the phylogenetic diversity of two novel species belonging to the genus Bifidobacterium: Bifidobacterium cebidarum sp. Nov. and Bifidobacterium leontopitheci sp. nov.. Int. J. Syst. Evol. Microbiol. 70, 2288–2297 (2020).CAS 
PubMed 
Article 

Google Scholar 
50.Modesto, M. et al. Bifidobacterium primatium sp. nov., Bifidobacterium scaligerum sp. nov., Bifidobacterium felsineum sp. nov. and Bifidobacterium simiarum sp. nov.: Four novel taxa isolated from the faeces of the cotton top tamarin (Saguinus oedipus) and the emperor tamarin (Saguinus imperator). Syst. Appl. Microbiol. (2018).51.Neuzil-Bunesova, V. et al. Bifidobacterium canis sp nov a novel member of the Bifidobacterium pseudolongum phylogenetic group isolated from faeces of a dog (Canis lupus f. familiaris). Int. J. Syst. Evol. Microbiol. 70, 5040–5047 (2020).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
52.Vlková, E. et al. A new medium containing mupirocin, acetic acid, and norfloxacin for the selective cultivation of bifidobacteria. Anaerobe 34, 27–33 (2015).PubMed 
Article 
CAS 
PubMed Central 

Google Scholar 
53.Carding, S., Verbeke, K., Vipond, D. T., Corfe, B. M. & Owen, L. J. Dysbiosis of the gut microbiota in disease. Microb. Ecol. Health Dis. 26, 26191 (2015).
Google Scholar 
54.WagnerMackenzie, B. et al. Bacterial community collapse: a meta-analysis of the sinonasal microbiota in chronic rhinosinusitis. Environ. Microbiol. 19, 381–392 (2017).CAS 
Article 

Google Scholar 
55.Arboleya, S., Watkins, C., Stanton, C. & Ross, R. P. Gut bifidobacteria populations in human health and aging. Front. Microbiol. 7 (2016).56.Binda, C. et al. Actinobacteria: a relevant minority for the maintenance of gut homeostasis. Dig. Liver Dis. 50, 421–428 (2018).MathSciNet 
PubMed 
Article 
PubMed Central 

Google Scholar 
57.Tojo, R. et al. Intestinal microbiota in health and disease: role of bifidobacteria in gut homeostasis. World J. Gastroenterol. 20, 15163 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
58.Rodriguez, C. I. & Martiny, J. B. H. Evolutionary relationships among bifidobacteria and their hosts and environments. BMC Genom. 21, 1–12 (2020).Article 

Google Scholar 
59.Sharma, V., Mobeen, F. & Prakash, T. Exploration of survival traits, probiotic determinants, host interactions, and functional evolution of bifidobacterial genomes using comparative genomics. Genes 9, 477 (2018).PubMed Central 
Article 
CAS 

Google Scholar 
60.Sun, Z. et al. Comparative genomic analysis of 45 type strains of the genus Bifidobacterium. a snapshot of its genetic diversity and evolution. PLoS One 10, 0117912 (2015).
Google Scholar 
61.Frey, J. C. et al. Fecal bacterial diversity in a wild gorilla. Appl. Environ. Microbiol. 72, 3788–3792 (2006).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
62.Makovska, M., Modrackova, N., Bolechova, P., Drnkova, B. & Neuzil-Bunesova, V. Antibiotic susceptibility screening of primate-associated Clostridium ventriculi. Anaerobe, 102347 (2021).63.Ushida, K. et al. Draft genome sequences of Sarcina ventriculi strains isolated from wild Japanese macaques in Yakushima Island. Genome announcements 4 (2016).64.Owens, L. A. et al. A Sarcina bacterium linked to lethal disease in sanctuary chimpanzees in Sierra Leone. Nat. Commun. 12, 1–16 (2021).ADS 
Article 
CAS 

Google Scholar 
65.Vlková, E., Rada, V., Šmehilová, M. & Killer, J. Auto-aggregation and co-aggregation ability in bifidobacteria and clostridia. Folia Microbiol. 53, 263–269 (2008).Article 
CAS 

Google Scholar 
66.Wang, L. et al. Adhesive Bifidobacterium induced changes in cecal microbiome alleviated constipation in mice. Front. Microbiol. 10, 1721 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
67.Wei, Y. et al. Protective effects of bifidobacterial strains against toxigenic Clostridium difficile. Front. Microbiol. 9, 888 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
68.Guittar, J., Shade, A. & Litchman, E. Trait-based community assembly and succession of the infant gut microbiome. Nature Commun. 10, 1–11 (2019).Article 
CAS 

Google Scholar 
69.Moore, R. E. & Townsend, S. D. Temporal development of the infant gut microbiome. Open Biol. 9, 190128 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
70.Korpela, K. et al. Probiotic supplementation restores normal microbiota composition and function in antibiotic-treated and in caesarean-born infants. Microbiome 6, 1–11 (2018).Article 

Google Scholar 
71.Timperio, A. M., Gorrasi, S., Zolla, L. & Fenice, M. Evaluation of MALDI-TOF mass spectrometry and MALDI BioTyper in comparison to 16S rDNA sequencing for the identification of bacteria isolated from Arctic sea water. PloS One 12, 0181860 (2017).Article 
CAS 

Google Scholar 
72.Bäckhed, F. et al. Dynamics and stabilization of the human gut microbiome during the first year of life. Cell Host Microbe 17, 690–703 (2015).PubMed 
Article 
CAS 

Google Scholar 
73.Brown, C. J. et al. Comparative genomics of Bifidobacterium species isolated from marmosets and humans. Am. J. Primatol. 81, e983 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
74.Killer, J. et al. Gene encoding the CTP synthetase as an appropriate molecular tool for identification and phylogenetic study of the family Bifidobacteriaceae. MicrobiologyOpen 7, e00579 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
75.Milani, C. et al. Evaluation of bifidobacterial community composition in the human gut by means of a targeted amplicon sequencing (ITS) protocol. FEMS Microbiol. Ecol. 90, 493–503 (2014).CAS 
PubMed 

Google Scholar 
76.Srinivasan, R. et al. Use of 16S rRNA gene for identification of a broad range of clinically relevant bacterial pathogens. PloS One 10, e0117617 (2015).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
77.Maiden, M. C. J. et al. MLST revisited: the gene-by-gene approach to bacterial genomics. Nature Rev. Microbiol. 11, 728–736 (2013).CAS 
Article 

Google Scholar 
78.Lugli, G. A. et al. Phylogenetic classification of six novel species belonging to the genus Bifidobacterium comprising Bifidobacterium anseris sp. nov., Bifidobacterium criceti sp. nov., Bifidobacterium imperatoris sp. nov., Bifidobacterium italicum sp. nov., Bifidobacterium margollesii sp. nov. and Bifidobacterium parmae sp. nov. Syst. Appl. Microbiol. 41, 173–183 (2018).PubMed 
Article 

Google Scholar 
79.Malukiewicz, J. et al. The effects of host taxon, hybridization, and environment on the gut microbiome of Callithrix marmosets. BioRxiv, 708255 (2019).80.Amato, K. R. et al. Phylogenetic and ecological factors impact the gut microbiota of two Neotropical primate species. Oecologia 180, 717–733 (2016).ADS 
PubMed 
Article 
PubMed Central 

Google Scholar 
81.Hernández‐Rodríguez, D., Vásquez‐Aguilar, A. A., Serio‐Silva, J. C., Rebollar, E. A. & Azaola‐Espinosa, A. Molecular detection of Bifidobacterium spp. in faeces of black howler monkeys (Alouatta pigra). J. Med. Primatol. 48, 99–105 (2019).82.Zhu, L. et al. Sex bias in gut microbiome transmission in newly paired marmosets (Callithrix jacchus). Msystems 5, e00910-00919 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
83.Kap, Y. S. et al. Targeted diet modification reduces multiple sclerosis–like disease in adult marmoset monkeys from an outbred colony. J. Immunol. 201, 3229–3243 (2018).CAS 
PubMed 
Article 

Google Scholar 
84.Ren, T., Grieneisen, L. E., Alberts, S. C., Archie, E. A. & Wu, M. Development, diet and dynamism: longitudinal and cross-sectional predictors of gut microbial communities in wild baboons. Environ. Microbiol. 18, 1312–1325 (2016).PubMed 
Article 

Google Scholar 
85.Xu, B. et al. Metagenomic analysis of the Rhinopithecus bieti fecal microbiome reveals a broad diversity of bacterial and glycoside hydrolase profiles related to lignocellulose degradation. BMC Genom. 16, 1–11 (2015).Article 
CAS 

Google Scholar 
86.Baumann, P. Biology of bacteriocyte-associated endosymbionts of plant sap-sucking insects. Annu. Rev. Microbiol. 59, 155–189 (2005).CAS 
PubMed 
Article 

Google Scholar 
87.Killer, J. et al. Bifidobacterium actinocoloniiforme sp. nov. and Bifidobacterium bohemicum sp. nov., from the bumblebee digestive tract. Int. J. Syst. Evol. Microbiol. 61, 1315–1321 (2011).88.Amato, K. R. et al. Evolutionary trends in host physiology outweigh dietary niche in structuring primate gut microbiomes. ISME J. 13, 576–587 (2019).CAS 
PubMed 
Article 

Google Scholar 
89.Garber, P. A., Mallott, E. K., Porter, L. M. & Gomez, A. The gut microbiome and metabolome of saddleback tamarins (Leontocebus weddelli): Insights into the foraging ecology of a small‐bodied primate. Am. J. Primatol. 81, e23003 (2019).90.Gralka, M., Szabo, R., Stocker, R. & Cordero, O. X. Trophic interactions and the drivers of microbial community assembly. Curr. Biol. 30, R1176–R1188 (2020).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
91.Clayton, J. B. et al. Associations between nutrition, gut microbiome, and health in a novel nonhuman primate model. Sci. Rep. 8, 1–16 (2018).CAS 
Article 

Google Scholar 
92.Koo, B. S. et al. Idiopathic chronic diarrhea associated with dysbiosis in a captive cynomolgus macaque (Macaca fascicularis). J. Med. Primatol. 49, 56–59 (2020).PubMed 
Article 
PubMed Central 

Google Scholar 
93.Krynak, K. L., Burke, D. J., Martin, R. A. & Dennis, P. M. Gut microbiome composition is associated with cardiac disease in zoo-housed western lowland gorillas (Gorilla gorilla gorilla). FEMS Microbiol. Lett. 364 (2017).94.Modrackova, N. et al. Prebiotic potential of natural gums and starch for bifidobacteria of variable origins. Bioact. Carbohydr. Diet. Fibre 20, 100199 (2019).95.McKenzie, V. J., Kueneman, J. G. & Harris, R. N. Probiotics as a tool for disease mitigation in wildlife: insights from food production and medicine. Ann. N. Y. Acad. Sci. 1429, 18–30 (2018).ADS 
PubMed 
Article 
PubMed Central 

Google Scholar 
96.Hicks, A. L. et al. Gut microbiomes of wild great apes fluctuate seasonally in response to diet. Nat. Commun. 9, 1–18 (2018).CAS 
Article 

Google Scholar 
97.Hungate, R. E. & Macy, J. The roll-tube method for cultivation of strict anaerobes. Bulletins from the ecological research committee, 123–126 (1973).98.Rada, V. & Petr, J. A new selective medium for the isolation of glucose non-fermenting bifidobacteria from hen caeca. J. Microbiol. Methods 43, 127–132 (2000).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
99.Orban, J. I. & Patterson, J. A. Modification of the phosphoketolase assay for rapid identification of bifidobacteria. J. Microbiol. Methods 40, 221–224 (2000).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
100.Kim, B. J., Kim, H.-Y., Yun, Y.-J., Kim, B.-J. & Kook, Y.-H. Differentiation of Bifidobacterium species using partial RNA polymerase β-subunit (rpoB) gene sequences. Int. J. Syst. Evol. Microbiol. 60, 2697–2704 (2010).PubMed 
Article 
PubMed Central 

Google Scholar 
101.Hall, T. A. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. 41 edn 95–98 ([London]: Information Retrieval Ltd., c1979-c2000.).102.Thompson, J. D., Higgins, D. G. & Gibson, T. J. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22, 4673–4680 (1994).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
103.Callahan, B. J. et al. DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
104.Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Ress 41, D590–D596 (2012).Article 
CAS 

Google Scholar 
105.Shannon, C. E. & Weaver, W. The mathematical theory of information. Urbana: University of Illinois Press 97 (1949).106.Pielou, E. C. The measurement of diversity in different types of biological collections. J. Theor. Biol. 13, 131–144 (1966).ADS 
Article 

Google Scholar 
107.Mandal, S. et al. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb. Ecol. Health Dis. 26, 27663 (2015).PubMed 
PubMed Central 

Google Scholar 
108.fundamental algorithms for scientific computing in Python. Virtanen, P. et al. SciPy 1.0. Nat. Methods 17, 261–272 (2020).Article 
CAS 

Google Scholar 
109.Seabold, S. & Perktold, J. Statsmodels: Econometric and statistical modeling with python in Proceedings of the 9th Python in Science Conference 57 (Austin, TX, 2010).110.MacKinnon, J. G. & White, H. Some heteroskedasticity-consistent covariance matrix estimators with improved finite sample properties. J. Econom. 29, 305–325 (1985).Article 

Google Scholar  More