Philippa Kaur

Basic information of intervieweesResults of the descriptive analysis summarized in Table 2 show that more than half of the respondents were males (69%) and were on average 41.3 years old while more than 32 years of farming experience. The study area is comprised of multiple ethnic groups (Han, Tibetan, Yugur, Mongolian, Hui, etc.). In most cases, the main livelihood activity of the Ethnic Minorities (Tibetan, Yugur, Mongolian, Hui, etc.) is livestock, while Han people main livelihood activity is farming. The majority of respondents (64%) were minority nationality. The vast majority of the agro-pastoralists (86%) have a primary school education or above, even though only 1% of them have Undergraduate education or Above. The results also reveal that 92% of respondents have access to weather information. The average cultivated land Per household is 10.23 Mu and Grassland is 156.21 Mu, respectively. The average per household income is RMB78000, and agricultural income is RMB52000.Table 2 Descriptive statistics of agro-pastoralist characteristics.Full size tableDue to their long-term farming experience, the agro-pastoralists were expected to have a high-level of understanding of local climate knowledge. Also contributing to this could be the information they receive about climate change and for some, the associated training through agro-pastoralists’ associations. Therefore, they also have a propensity to adapt to adverse conditions resulting from climate change impacts. In addition, the high-level of farming experience, the cultivated-land size, grassland size, Credit loan, Insurance, Village cadres all have a positive impact on the level of agro-pastoralists’ adaptation to new climate scenarios.However, the education level and cadres experience may be the major limiting factors for adopting specific long-term adaptation strategies. Ethnicity and gender are also expected to be key factors influencing awareness and adaptation to climate change. There are differences in relative perception intensity between Ethnic Minority and Han because of their cultural ecology (the main livelihood activity of minorities nationality is livestock, while Han main livelihood activity is farming.). In terms of gender, women in rural areas are less mobile and have less access to information and rights. They are also heavily involved in domestic work. However, men may have easier access to information (socializing, going out to work, etc.) Therefore, male headed households are expected to be more likely to adapt to the impact of climate change.Climate change trend in the study areaFigure 2 shows the trend of annual precipitation, annual rainfall and annual snow at different meteorological stations in the study area. As shown in the Fig. 2, precipitation, rainfall and snow show an increasing trend, but the increase range of snow (0.0325–0.375/a) is significantly lower than that of precipitation (1.22–3.1/a) and rainfall (1.04–2.81/a). Similarly, through the inspection, it is found that the multi-collinearity among precipitation, rainfall and snow at each meteorological station is obvious (most R2  > 0.5, and p  More

Hanski, I. Habitat fragmentation and species richness. J. Biogeogr. 42, 989–993. https://doi.org/10.1111/jbi.12478 (2015).Article 

Google Scholar 
Hoffmann, M. et al. The impact of conservation on the status of the world’s vertebrates. Science 330, 1503–1509. https://doi.org/10.1126/science.1194442 (2010).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Pimm, S. L. et al. The biodiversity of species and their rates of extinction, distribution, and protection. Science 344, 1246752 (2014).CAS 
Article 

Google Scholar 
Crooks, K. R. Relative sensitivities of mammalian carnivores to habitat fragmentation. Conserv. Biol. 16, 488–502. https://doi.org/10.1046/j.1523-1739.2002.00386.x (2002).Article 

Google Scholar 
Elliot, N. B., Cushman, S. A., Macdonald, D. W. & Loveridge, A. J. The devil is in the dispersers: Predictions of landscape connectivity change with demography. J. Appl. Ecol. 51, 1169–1178. https://doi.org/10.1111/1365-2664.12282 (2014).Article 

Google Scholar 
Carroll, C. Interacting effects of climate change, landscape conversion, and harvest on carnivore populations at the range margin: Marten and lynx in the northern Appalachians. Conserv. Biol. 21, 1092–1104. https://doi.org/10.1111/j.1523-1739.2007.00719.x (2007).Article 
PubMed 

Google Scholar 
Ripple, W. J. et al. Status and ecological effects of the world’s largest carnivores. Science 343, 1241484. https://doi.org/10.1126/science.1241484 (2014).CAS 
Article 
PubMed 

Google Scholar 
Farris, Z. J. et al. Hunting, exotic carnivores, and habitat loss: Anthropogenic effects on a native carnivore community, Madagascar. PLOS ONE 10, e0136456. https://doi.org/10.1371/journal.pone.0136456 (2015).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Farris, Z. J. et al. Threats to a rainforest carnivore community: A multi-year assessment of occupancy and co-occurrence in Madagascar. Biol. Cons. 210, 116–124. https://doi.org/10.1016/j.biocon.2017.04.010 (2017).Article 

Google Scholar 
Swihart, R. K., Gehring, T. M., Kolozsvary, M. B. & Nupp, T. E. Responses of “resistant” vertebrates to habitat loss and fragmentation: The importance of niche breadth and range boundaries. Divers. Distrib. 9, 1–18. https://doi.org/10.1046/j.1472-4642.2003.00158.x (2003).Article 

Google Scholar 
Caryl, F. M., Quine, C. P. & Park, K. J. Martens in the matrix: The importance of nonforested habitats for forest carnivores in fragmented landscapes. J. Mammal. 93, 464–474. https://doi.org/10.1644/11-MAMM-A-149.1 (2012).Article 

Google Scholar 
Pereboom, V. et al. Movement patterns, habitat selection, and corridor use of a typical woodland-dweller species, the European pine marten (Martes martes), in fragmented landscape. Can. J. Zool. 86, 983–991. https://doi.org/10.1139/Z08-076 (2008).Article 

Google Scholar 
Fleschutz, M. M. et al. Response of a small felid of conservation concern to habitat fragmentation. Biodivers. Conserv. 25, 1447–1463. https://doi.org/10.1007/s10531-016-1118-6 (2016).Article 

Google Scholar 
Gálvez, N. et al. Forest cover outside protected areas plays an important role in the conservation of the Vulnerable guiña Leopardus guigna. Oryx 47, 251–258. https://doi.org/10.1017/S0030605312000099 (2013).Article 

Google Scholar 
Belcher, C. A. Demographics of tiger quoll (Dasyurus maculatus maculatus) populations in south-eastern Australia. Aust. J. Zool. 51, 611–626. https://doi.org/10.1071/ZO02051 (2003).Article 

Google Scholar 
Maxwell, S., Burbidge, A. & Morris, K. Spotted-tailed Quoll (SE mainland and Tas); recovery outline. (1996).Jones, M. E., Rose, R. K. & Burnett, S. Dasyurus maculatus. Mammalian Species 676, 1–9 (2001).Article 

Google Scholar 
Long, K. & Nelson, J. National recovery plan for the spotted-tailed Quoll Dasyurus maculatus. Victorian Department of Sustainability and Environment (2010).Claridge, A. W. et al. Home range of the spotted-tailed quoll (Dasyurus maculatus), a marsupial carnivore, in a rainshadow woodland. Wildl. Res. 32, 7–14. https://doi.org/10.1071/WR04031 (2005).Article 

Google Scholar 
Glen, A. S. & Dickman, C. R. Home range, denning behaviour and microhabitat use of the carnivorous marsupial Dasyurus maculatus in eastern Australia. J. Zool. 268, 347–354. https://doi.org/10.1111/j.1469-7998.2006.00064.x (2006).Article 

Google Scholar 
Körtner, G. et al. Population structure, turnover and movement of spotted-tailed quolls on the New England Tablelands. Wildl. Res. 31, 475–484. https://doi.org/10.1071/WR03041 (2004).Article 

Google Scholar 
Belcher, C. The Largest Surviving Marsupial Carnivore on Mainland Australia: The Tiger or Spotted-Tailed Quoll Dasyurus maculatus, A Nationally Threatened, Forest-Dependent Species 612–623 (Royal Zoological Society of New South Wales, Sydney, 2004).
Google Scholar 
Henderson, T., Fancourt, B. A., Rajaratnam, R., Vernes, K. & Ballard, G. Spatial and temporal interactions between endangered spotted-tailed quolls and introduced red foxes in a fragmented landscape. J. Zool. https://doi.org/10.1111/jzo.12919 (2021).Article 

Google Scholar 
Troy, S. N. Spatial Ecology of the Tasmanian Spotted-Tailed Quoll. Ph.D. Thesis, University of Tasmania, (2014).Jones, M. E. et al. Research supporting restoration aiming to make a fragmented landscape ‘functional’ for native wildlife. Ecol. Manag. Restor. 22, 65–74. https://doi.org/10.1111/emr.12504 (2021).Article 

Google Scholar 
Andersen, G. E., Johnson, C. N., Barmuta, L. A. & Jones, M. E. Use of anthropogenic linear features by two medium-sized carnivores in reserved and agricultural landscapes. Scientific Reports 7, 1–11. https://doi.org/10.1038/s41598-017-11454-z (2017).CAS 
Article 

Google Scholar 
Nichols, J. D. in Applied Ecology and Human Dimensions in Biological Conservation (eds L. M. Verdade, M.C. Lyra-Jorge, & C.I. Pina) 117–131 (Springer, 2014).Royle, J. A., Chandler, R. B., Sollmann, R. & Gardner, B. in Spatial Capture-Recapture (eds J. Andrew Royle, Richard B. Chandler, Rahel Sollmann, & Beth Gardner) 3–19 (Academic Press, 2014).Sollmann, R., Gardner, B. & Belant, J. L. How does spatial study design influence density estimates from spatial capture-recapture models?. PLoS ONE 7, e34575 (2012).ADS 
CAS 
Article 

Google Scholar 
Kalle, R., Ramesh, T., Qureshi, Q. & Sankar, K. Density of tiger and leopard in a tropical deciduous forest of Mudumalai Tiger Reserve, southern India, as estimated using photographic capture–recapture sampling. Acta Theriol. 56, 335–342. https://doi.org/10.1007/s13364-011-0038-9 (2011).Article 

Google Scholar 
Vissia, S., Wadhwa, R. & van Langevelde, F. Co-occurrence of high densities of brown hyena and spotted hyena in central Tuli, Botswana. J. Zool. 314, 143–150. https://doi.org/10.1111/jzo.12873 (2021).Article 

Google Scholar 
Henderson, T., Fancourt, B. A. & Ballard, G. The importance of species-specific survey designs: Prey camera trap surveys significantly underestimate the detectability of endangered spotted-tailed quolls. Aust. Mammalogy https://doi.org/10.1071/AM21039 (2022).Gorta, S. B. Z., Alting, B., Claridge, A. & Henderson, T. Apparent piebald variants in quolls (Dasyurus): Examples of three recent cases in the spotted-tailed quoll Dasyurus maculatus. Aust. Mammalogy 43, 373–377. https://doi.org/10.1071/AM20058 (2021).Article 

Google Scholar 
Kowalksi, M. (https://exifpro.informer.com/2.1/, 2011).Efford, M. in R package version 4.5.3 (2022).R Core Team. (R Foundation for Statistical Computing, Vienna, Austria, 2022).Rovero, F. & Zimmermann, F. Camera Trapping for Wildlife Research (Pelagic Publishing Ltd, London, 2016).
Google Scholar 
Efford, M. Density estimation in live-trapping studies. Oikos 106, 598–610. https://doi.org/10.1111/j.0030-1299.2004.13043.x (2004).Article 

Google Scholar 
Niedballa, J., Sollmann, R., Courtiol, A. & Wilting, A. camtrapR: An R package for efficient camera trap data management. Methods Ecol. Evol. 7, 1457–1462. https://doi.org/10.1111/2041-210X.12600 (2016).Article 

Google Scholar 
Burnham, K. P. & Anderson, D. R. A practical information-theoretic approach. Model Sel. Multimodel Inference 2, 70–71 (2002).MATH 

Google Scholar 
Hamer, R. P. et al. Differing effects of productivity on home-range size and population density of a native and an invasive mammalian carnivore. Wildlife Res. 49, 158–168. https://doi.org/10.1071/WR20134 (2021).Article 

Google Scholar 
Glen, A. S. & Dickman, C. R. Complex interactions among mammalian carnivores in Australia, and their implications for wildlife management. Biol. Rev. 80, 387–401. https://doi.org/10.1017/s1464793105006718 (2005).Article 
PubMed 

Google Scholar 
Glen, A. S., Pennay, M., Dickman, C. R., Wintle, B. A. & Firestone, K. B. Diets of sympatric native and introduced carnivores in the Barrington Tops, eastern Australia. Austral Ecol. 36, 290–296. https://doi.org/10.1111/j.1442-9993.2010.02149.x (2011).Article 

Google Scholar 
Glen, A. S. & Dickman, C. R. Population viability analysis shows spotted-tailed quolls may be vulnerable to competition. Aust Mammalogy 35, 180–183. https://doi.org/10.1071/AM12045 (2013).Article 

Google Scholar 
Graham, C. A., Maron, M. & McAlpine, C. A. Influence of landscape structure on invasive predators: Feral cats and red foxes in the brigalow landscapes, Queensland Australia. Wildl. Res. 39, 661–676. https://doi.org/10.1071/WR12008 (2012).Article 

Google Scholar 
Glen, A. S. Population attributes of the spotted-tailed quoll (Dasyurus maculatus) in north-eastern New South Wales. Aust. J. Zool. 56, 137–142. https://doi.org/10.1071/ZO08025 (2008).Article 

Google Scholar 
Chua, M. A., Sivasothi, N. & Meier, R. Population density, spatiotemporal use and diet of the leopard cat (Prionailurus bengalensis) in a human-modified succession forest landscape of Singapore. Mammal Res. 61, 99–108 (2016).Article 

Google Scholar 
Lorica, M. & Heaney, L. Survival of a native mammalian carnivore, the leopard cat Prionailurus bengalensis Kerr, 1792 (Carnivora: Felidae), in an agricultural landscape on an oceanic Philippine island. J. Threatened Taxa, 4451–4460 (2013).Rajaratnam, R., Sunquist, M., Rajaratnam, L. & Ambu, L. Diet and habitat selection of the leopard cat (Prionailurus bengalensis borneoensis) in an agricultural landscape in Sabah, Malaysian Borneo. J. Trop. Ecol. 23, 209–217 (2007).Article 

Google Scholar 
Belcher, C. A. & Darrant, J. P. Den use by the spotted-tailed quoll Dasyurus maculatus in south-eastern Australia. Aust Mammalogy 28, 59–64. https://doi.org/10.1071/AM06007 (2006).Article 

Google Scholar 
Glen, A. & Dickman, C. Why are there so many spotted-tailed Quolls Dasyurus maculatus in parts of north-eastern New South Wales?. Aust Zool 35, 711–718. https://doi.org/10.7882/az.2011.023 (2011).Article 

Google Scholar 
Hanski, I. Metapopulation ecology (Oxford University Press, Oxford, 1999).
Google Scholar 
Pulliam, H. R. Sources, sinks, and population regulation. Am. Nat. 132, 652–661 (1988).Article 

Google Scholar 
Belcher, C. A. Susceptibility of the tiger quoll, Dasyurus maculatus, and the eastern quoll, D. viverrinus, to 1080-poisoned baits in control programmes for vertebrate pests in eastern Australia. Wildl. Res. 25, 33–40. https://doi.org/10.1071/WR95077 (1998).Article 

Google Scholar 
Schmidt, G. M., Graves, T. A., Pederson, J. C. & Carroll, S. L. Precision and bias of spatial capture–recapture estimates: A multi-site, multi-year Utah black bear case study. Ecological Applications 32, e2618. https://doi.org/10.1002/eap.2618 (2022).Article 
PubMed 
PubMed Central 

Google Scholar 
White, G. C. Capture-Recapture and Removal Methods for Sampling Closed Populations (Los Alamos National Laboratory, New Mexico, 1982).
Google Scholar 
Thornton, D. H. & Pekins, C. E. Spatially explicit capture–recapture analysis of bobcat (Lynx rufus) density: Implications for mesocarnivore monitoring. Wildl. Res. 42, 394–404. https://doi.org/10.1071/WR15092 (2015).Article 

Google Scholar 
Sollmann, R. et al. Improving density estimates for elusive carnivores: Accounting for sex-specific detection and movements using spatial capture–recapture models for jaguars in central Brazil. Biol. Cons. 144, 1017–1024 (2011).Article 

Google Scholar 
Green, A. M., Chynoweth, M. W. & Şekercioğlu, Ç. H. Spatially explicit capture-recapture through camera trapping: A review of benchmark analyses for wildlife density estimation. Front. Ecol. Evol. 8, 473. https://doi.org/10.3389/fevo.2020.563477 (2020).Article 

Google Scholar 
du Preez, B. D., Loveridge, A. J. & Macdonald, D. W. To bait or not to bait: a comparison of camera-trapping methods for estimating leopard Panthera pardus density. Biol. Cons. 176, 153–161 (2014).Article 

Google Scholar 
Zimmermann, F., Breitenmoser-Würsten, C., Molinari-Jobin, A. & Breitenmoser, U. Optimizing the size of the area surveyed for monitoring a Eurasian lynx (Lynx lynx) population in the Swiss Alps by means of photographic capture–recapture. Integr. Zool. 8, 232–243 (2013).Article 

Google Scholar 
Dupont, P., Milleret, C., Gimenez, O. & Bischof, R. Population closure and the bias-precision trade-off in spatial capture–recapture. Methods Ecol. Evol. 10, 661–672. https://doi.org/10.1111/2041-210X.13158 (2019).Article 

Google Scholar 
Mergey, M., Helder, R. & Roeder, J. J. Effect of forest fragmentation on space-use patterns in the European pine marten (Martes martes). J. Mammal. 92, 328–335. https://doi.org/10.1644/09-MAMM-A-366.1 (2011).Article 

Google Scholar 
Silmi, M. et al. Activity and ranging behavior of leopard cats (Prionailurus bengalensis) in an oil palm landscape. Frontiers in Environmental Science 9, 651939. https://doi.org/10.3389/fenvs.2021.651939 (2021).Article 

Google Scholar  More

Ecosystems and Human Well-being: Synthesis (Millennium Ecosystem Assessment, 2005).Huang, J. et al. Dryland climate change: recent progress and challenges. Rev. Geophys. 55, 719–778 (2017).Article 

Google Scholar 
Fu, B. et al. The Global-DEP conceptual framework — research on dryland ecosystems to promote sustainability. Curr. Opin. Environ. Sustain. 48, 17–28 (2021).Article 

Google Scholar 
He, C. et al. Detecting global urban expansion over the last three decades using a fully convolutional network. Environ. Res. Lett. 14, 034008 (2019).Article 

Google Scholar 
Güneralp, B., Reba, M., Hales, B. U., Wentz, E. A. & Seto, K. C. Trends in urban land expansion, density, and land transitions from 1970 to 2010: a global synthesis. Environ. Res. Lett. 15, 044015 (2020).Article 

Google Scholar 
McDonald, R. I. et al. Research gaps in knowledge of the impact of urban growth on biodiversity. Nat. Sustain. 3, 16–24 (2019).Article 

Google Scholar 
Liu, X. et al. High-spatiotemporal-resolution mapping of global urban change from 1985 to 2015. Nat. Sustain. 3, 564–570 (2020).Article 

Google Scholar 
Güneralp, B. & Seto, K. C. Futures of global urban expansion: uncertainties and implications for biodiversity conservation. Environ. Res. Lett. 8, 014025 (2013).Article 

Google Scholar 
McDonald, R. I., Kareiva, P. & Forman, R. T. T. The implications of current and future urbanization for global protected areas and biodiversity conservation. Biol. Conserv. 141, 1695–1703 (2008).Article 

Google Scholar 
McDonald, R. I., Marcotullio, P. J. & Güneralp, B. Urbanization, Biodiversity and Ecosystem Services: Challenges and Opportunities (Springer, 2013).van Vliet, J. Direct and indirect loss of natural area from urban expansion. Nat. Sustain. 2, 755–763 (2019).Article 

Google Scholar 
Sharp, R. et al. InVEST 3.2.0 User’s Guide (The Natural Capital Project, Stanford Univ., Univ. Minnesota, The Nature Conservancy and World Wildlife Fund, 2015).Terrado, M. et al. Model development for the assessment of terrestrial and aquatic habitat quality in conservation planning. Sci. Total Environ. 540, 63–70 (2016).CAS 
Article 

Google Scholar 
Bai, Y. et al. Developing China’s Ecological Redline Policy using ecosystem services assessments for land use planning. Nat. Commun. 9, 3034 (2018).Article 
CAS 

Google Scholar 
McDonald, R. I. et al. Urban effects, distance, and protected areas in an urbanizing world. Landsc. Urban Plan. 93, 63–75 (2009).Article 

Google Scholar 
Mirzabaev, A. et al. in Climate Change and Land (eds Shukla, P. R. et al.) 249–343 (IPCC, 2019).Friis, C. & Nielsen, J. Telecoupling. Exploring Land-use Change in a Globalised World (Palgrave Macmillan, 2019).Maestre, F. et al. Structure and functioning of dryland ecosystems in a changing world. Annu. Rev. Ecol. Evol. Syst. 47, 215–237 (2016).Article 

Google Scholar 
Leh, M. D. K., Matlock, M. D., Cummings, E. C. & Nalley, L. L. Quantifying and mapping multiple ecosystem services change in West Africa. Agric. Ecosyst. Environ. 165, 6–18 (2013).Article 

Google Scholar 
Xie, W., Huang, Q., He, C. & Zhao, X. Projecting the impacts of urban expansion on simultaneous losses of ecosystem services: a case study in Beijing, China. Ecol. Indic. 84, 183–193 (2018).Article 

Google Scholar 
Whitford, W. & Wade, E. L. Ecology of Desert Systems (Academic Press, 2002).Brito, J. C. et al. Conservation biogeography of the Sahara‐Sahel: additional protected areas are needed to secure unique biodiversity. Divers. Distrib. 22, 371–384 (2016).Article 

Google Scholar 
Jenkins, C. N., Pimm, S. L. & Joppa, L. N. Global patterns of terrestrial vertebrate diversity and conservation. Proc. Natl Acad. Sci. USA 110, E2602–E2610 (2013).CAS 
Article 

Google Scholar 
Salafsky, N. et al. A standard lexicon for biodiversity conservation: unified classifications of threats and actions. Conserv. Biol. 22, 897–911 (2008).Article 

Google Scholar 
Oliver, T. H. et al. Declining resilience of ecosystem functions under biodiversity loss. Nat. Commun. 6, 10122 (2015).Article 
CAS 

Google Scholar 
Powers, R. P. & Jetz, W. Global habitat loss and extinction risk of terrestrial vertebrates under future land-use-change scenarios. Nat. Clim. Change 9, 323–329 (2019).Article 

Google Scholar 
Díaz, S. M. et al. The Global Assessment Report on Biodiversity and Ecosystem Services: Summary for Policy Makers (IPBES, 2019).Seto, K. C., Guneralp, B. & Hutyra, L. R. Global forecasts of urban expansion to 2030 and direct impacts on biodiversity and carbon pools. Proc. Natl Acad. Sci. USA 109, 16083–16088 (2012).CAS 
Article 

Google Scholar 
Pautasso, M. Scale dependence of the correlation between human population presence and vertebrate and plant species richness. Ecol. Lett. 10, 16–24 (2007).Article 

Google Scholar 
Luck, G. W. A review of the relationships between human population density and biodiversity. Biol. Rev. Camb. Phil. Soc. 82, 607–645 (2007).Article 

Google Scholar 
McDonald, R. I., Güneralp, B., Huang, C.-W., Seto, K. C. & You, M. Conservation priorities to protect vertebrate endemics from global urban expansion. Biol. Conserv. 224, 290–299 (2018).Article 

Google Scholar 
The IUCN Red List of Threatened Species Version 2017-3 (IUCN, 2017); https://www.iucnredlist.org/resources/spatial-data-downloadTucker, M. A. et al. Moving in the Anthropocene: global reductions in terrestrial mammalian movements. Science 359, 466–469 (2018).CAS 
Article 

Google Scholar 
Howard, C., Flather, C. H. & Stephens, P. A. A global assessment of the drivers of threatened terrestrial species richness. Nat. Commun. 11, 993 (2020).CAS 
Article 

Google Scholar 
Guidelines for Geoconservation in Protected and Conserved Areas (IUCN, 2020).Gao, J. How China will protect one-quarter of its land. Nature 569, 457 (2019).CAS 
Article 

Google Scholar 
Chen, C. et al. China and India lead in greening of the world through land-use management. Nat. Sustain. 2, 122–129 (2019).Article 

Google Scholar 
Gao, B., Huang, Q., He, C., Sun, Z. & Zhang, D. How does sprawl differ across cities in China? A multi-scale investigation using nighttime light and census data. Landsc. Urban Plan. 148, 89–98 (2016).Article 

Google Scholar 
Mace, G. M. et al. Aiming higher to bend the curve of biodiversity loss. Nat. Sustain. 1, 448–451 (2018).Article 

Google Scholar 
Lambin, E. A. & Meyfroidt, P. Global land use change, economic globalization, and the looming land scarcity. Proc. Natl Acad. Sci. USA 108, 3465–3472 (2011).CAS 
Article 

Google Scholar 
Arlidge, W. et al. A global mitigation hierarchy for nature conservation. Bioscience 68, 336–347 (2018).Article 

Google Scholar 
Moallemi, E. A., Kwakkel, J., de Haan, F. J. & Bryan, B. A. Exploratory modeling for analyzing coupled human-natural systems under uncertainty. Glob. Environ. Change 65, 102186 (2020).Article 

Google Scholar 
Luck, M. A., Jenerette, G. D., Wu, J. & Grimm, N. B. The urban funnel model and the spatially heterogeneous ecological footprint. Ecosystems 4, 782–796 (2001).Article 

Google Scholar 
Ramaswami, A. et al. A social‐ecological‐infrastructural systems framework for interdisciplinary study of sustainable city systems. J. Ind. Ecol. 16, 801–813 (2012).Article 

Google Scholar 
Boerema, A. et al. Soybean trade: balancing environmental and socio-economic impacts of an intercontinental market. PLoS ONE 11, e0155222 (2016).Article 
CAS 

Google Scholar 
Garrett, R. D., Lambin, E. F. & Naylor, R. L. Land institutions and supply chain configurations as determinants of soybean planted area and yields in Brazil. Land Use Policy 31, 385–396 (2013).Article 

Google Scholar 
Friess, D. A., Rogers, K., Lovelock, C. E., Krauss, K. W. & Shi, S. The state of the world’s mangrove forests: past, present, and future. Annu. Rev. Environ. Resour. 44, 89–115 (2019).Article 

Google Scholar 
Ferreira, A. C. & Lacerda, L. D. Degradation and conservation of Brazilian mangroves, status and perspectives. Ocean Coast. Manage. 125, 38–46 (2016).Article 

Google Scholar 
Richards, D. R. & Friess, D. A. Rates and drivers of mangrove deforestation in Southeast Asia, 2000–2012. Proc. Natl Acad. Sci. USA 113, 201510272 (2016).
Google Scholar 
García-Vega, D. & Newbold, T. Assessing the effects of land use on biodiversity in the world’s drylands and Mediterranean environments. Biodivers. Conserv. 29, 393–408 (2020).Article 

Google Scholar 
Martínez-Valderrama, J., Guirado, E. & Maestre, F. Desertifying deserts. Nat. Sustain. 3, 572–575 (2020).Article 

Google Scholar 
Maestre, F. et al. Biogeography of global drylands. New Phytol. 231, 540–558 (2021).Article 

Google Scholar 
United Nations Environment World Conservation Monitoring Centre. World dryland areas according to UNCCD and CBD definitions. https://resources.unep-wcmc.org/products/789fcac8959943ab9ed7a225e5316f08 (2022).Olson, D. M. et al. Terrestrial ecoregions of the world: a new map of life on Earth. Bioscience 51, 933–938 (2001).Article 

Google Scholar 
Goldewijk, K. K., Beusen, A., Doelman, J. & Stehfest, E. Anthropogenic land use estimates for the Holocene – HYDE 3.2. Earth Syst. Sci. Data 9, 927–953 (2017).Article 

Google Scholar 
Revision of World Urbanization Prospects (United Nations, 2018); https://esa.un.org/unpd/wupLand Cover CCI—Product User Guide Version 2.0. (European Space Agency, 2017); http://maps.elie.ucl.ac.be/CCI/viewer/index.phpGrekousis, G., Mountrakis, G. & Kavouras, M. An overview of 21 global and 43 regional land-cover mapping products. Int. J. Remote Sens. 36, 5309–5335 (2015).Article 

Google Scholar 
Xu, X., Jain, A. K. & Calvin, K. V. Quantifying the biophysical and socioeconomic drivers of changes in forest and agricultural land in South and Southeast Asia. Glob. Change Biol. 25, 2137–2151 (2019).Article 

Google Scholar 
Gong, P. et al. Annual maps of global artificial impervious area (GAIA) between 1985 and 2018. Remote Sens. Environ. 236, 111510 (2020).Article 

Google Scholar 
Huang, Q. et al. The occupation of cropland by global urban expansion from 1992 to 2016 and its implications. Environ. Res. Lett. 15, 084037 (2020).Article 

Google Scholar 
He, C., Liu, Z., Tian, J. & Ma, Q. Urban expansion dynamics and natural habitat loss in China: a multiscale landscape perspective. Glob. Change Biol. 20, 2886–2902 (2014).Article 

Google Scholar 
Di Febbraro, M. et al. Expert-based and correlative models to map habitat quality: which gives better support to conservation planning? Glob. Ecol. Conserv. 16, e00513 (2018).Article 

Google Scholar 
Anselin, L. Local Indicators of Spatial Association—LISA. Geogr. Anal. 27, 93–115 (2010).Article 

Google Scholar  More

Experimental designThis experiment was set up containing two levels of soil biodiversity (high and low soil biodiversity) and seven treatments considering the number of global change factors (GCFs) (0, 1, 2, 4, 6, 8, 10) (Table 1, Supplementary Fig. 1 and Supplementary Data 1). We used the dilution-to-extinction approach to create the high and low soil biodiversity treatments (Supplementary Methods). Soil dilution can lead to a gradual loss of rare soil microbes, which can simulate a realistic loss of soil biodiversity, because rare soil microbes are more sensitive to anthropogenic pressures, e.g., warming, nitrogen addition and drought15. We note that the low soil biodiversity treatment is a subset of the high biodiversity, as many rare species have been eliminated through the dilution; this approach will likely lead to relatively more tolerant microbes in the resulting communities.An increasing number of GCFs was created inspired by the experimental design of the studies on biodiversity-ecosystem function relationships, based on random sampling from a species pool5,6,14. The combination of multiple GCFs was replicated 15 times at each level by randomly selecting GCFs from a pool with 10 GCFs for each replicate (Table 1 and Supplementary Data 1). For each replicate of combined GCFs, there were identical GCF combinations between the high and low soil biodiversity treatments to avoid a confounding effect of GCF combination and soil biodiversity treatments. The pool of 10 GCFs included: warming, nitrogen deposition, drought, heavy metal pollution, plastic mulching film residues, salinity, agricultural fungicide, bactericide application, surfactant contaminant and soil compaction (Supplementary Methods). These GCFs frequently occur in intensively managed agroecosystems and are treated as anthropogenic pressures10,13,14,15.MicrocosmsThis experiment was conducted using 50 ml conical Mini Bioreactors (Product Number 431720, Corning Inc., NY) as experimental units (Supplementary Fig. 2). Each Mini Bioreactor has four vents in the cap, where a hydrophobic membrane avoids microbial contamination but allows gas exchange. We filled each Mini Bioreactor with 40.0 g (dry weight, d.w.) of soil in total, which received the appropriate treatments.Soil sterilization and inoculum preparationWe collected the field soil from the top 10 cm of an intensive farming system in Berlin (52.466°N, 13.303°E). Field soil was passed through a 2 mm mesh to remove large roots and stones. We sterilized 20 kg of soil for 90 min at 121 °C, and stored 2 kg of fresh soil at 4 °C. The dilution-to-extinction approach38,39,40,41,42 was used to create high and low soil biodiversity (Supplementary Fig. 1). A parent inoculum suspension was prepared by mixing 100 g of fresh soil with 200 ml of sterilized VE water. The sediment settled for 1 min. The upper 200 ml of soil suspension was treated as parent inoculum suspension. 50 ml of parent inoculum suspension was added to 500 g of sterilized soil in a plastic bag, and homogenized by turning the bag up and down 30 times to obtain the inoculum of high soil biodiversity. Another 5 ml of parent inoculum suspension was mixed with 45 ml of sterilized parent inoculum suspension to create the 10-1 dilution. This procedure was repeated five times to reach the 10-6 dilution. 50 ml of the 10-6 dilution was mixed and homogenized with 500 g of sterilized soil in a plastic bag to obtain the inoculum of low soil biodiversity. This whole dilution procedure was repeated five times to obtain 10 bags of soil inoculum (five bags for each soil biodiversity inoculum).Sterile water was added to each plastic bag to reach the water content of the fresh soil in the field. All bags were closed with a sterilized cotton plug and a rubber band to avoid microbial contamination but allow gas exchange42. All bags were incubated in a dark room at 20 °C until similar microbial abundance was observed between the high and low soil biodiversity inoculum. Soil inoculum was homogenized by shaking and turning the bags once a week. After incubation, 2.0 g of soil in each bag was collected and stored at −80 °C for DNA extraction. Quantitative real-time PCR (qPCR) was used to determine fungal and bacterial abundance. In the present study, it took two months to recover soil microbial biomass (Supplementary Fig. 3).The implementation of GCFs and harvestAgroecosystems, some of the most intensively managed ecosystems, are affected by the co-occurrence of multiple GCFs13,14,15. This study focused on GCFs that frequently occur in agroecosystems, including warming, nitrogen deposition, drought, heavy metal pollution, plastic film residues, salinity, agricultural fungicide and bactericide application, surfactant contaminant and soil compaction. We present the rationale for the 10 tested GCFs in the Supplementary Methods.Loading soils were used to achieve an effective mixing of chemical agents into 40.0 g soil in each Mini Bioreactor. We created separate ‘loading soil’ for each GCF with chemical addition by mixing an appropriate dose of a chemical agent with sterilized soil through careful homogenization. This was done to avoid exaggerated effects of more concentrated chemicals when mixed with soil. For each chemical, 1.0 g (d.w.) of loading soil contained an appropriate dose for 40.0 g soil in a Mini Bioreactor. For instance, 1 634 mg of NH4NO3 was mixed with 100 g (d.w.) of sterilized soil, to ensure that there was about 16.34 mg of NH4NO3 in 1.0 g of sterilized loading soil. We weighed 40.0 g (d.w.) of soils, including 1.0 g (d.w.) of each loading soil, 5.0 g (d.w.) of soil inoculum (high or low soil biodiversity), an appropriate amount of film (0 or 0.16 g plastic film) and sterilized soil, according to GCF combination for each experimental unit. We put 40.0 g of mixed soils into a clean and sterilized cup (200 ml) with a cap, and then homogenized it by turning the cup up and down for 5 min using a shaking machine (Heidolph Reax 2, Heidolph Instruments GmbH & CO. KG, Schwabach, Germany). After homogenization, 40.0 g of mixed soils was transferred to a Mini Bioreactor, and a mesh bag containing about 100 mg of dry Medicago lupulina leaves (65 °C for 72 h) was buried 1 cm below the soil surface. We used a stick to press soils in each Mini Bioreactor to simulate an ambient condition (1.3 g cm−3) or mechanical compaction (1.7 g cm−3) in farmland.For the warming treatment with an increment of 5.0 °C over the ambient temperature (20 °C), we wrapped heating cables (Exo Terra PT-2012; Hagen Deutschland GmbH & Co. KG, Holm, Germany) around the outside of the bioreactors. A set temperature was maintained by temperature controllers (Voltcraft ETC-902; Conrad Electronic SE, Hirschau, Germany), which can switch off and on heating cables depending on the real-time temperature of the outside surface of Mini Bioreactors. Mini Bioreactors were placed in beakers filled with sand to reduce thermal radiation from neighboring units. At the start of the experiment, we added suitable amounts of sterilized water to each experimental unit to reach 60% of water holding capacity for the non-drought treatment and 30% water holding capacity for the drought treatment.All Mini Bioreactors were incubated at 20.0 °C in the dark for six weeks before the final harvest. Because there was 2.0 g of weight loss on average in each Mini Bioreactor in the first three weeks, we added 2 ml of sterilized water to each Mini Bioreactor on the first day of the fourth week. During the final harvest, soil in each Mini Bioreactor was gently homogenized using a spoon, and then 2.0 g of fresh soil was collected and stored at −80 °C for DNA extraction; 5.0 g was stored at 4 °C for the determination of soil enzyme activity; the leftover was oven-dried at 40 °C for other measurements. DNA of each soil sample was extracted from 250 mg soil, using DNeasy PowerSoil Pro Kit (QIAGEN GmbH, Hilden, Germany), following manufacturer’s instructions. Soil DNA extraction was stored at −80 °C for further analysis.The measurement of response variablesWe measured the following response variables: microbial activity (soil respiration), microbial abundance (bacterial and fungal abundance), nutrient cycling (litter decomposition rate and soil enzyme activity), physical properties (water-stable soil aggregates and soil water repellency), bacterial and fungal biodiversity (richness and microbial network features) (See details in the Supplementary Methods). We measured soil respiration as CO2 concentration (ppm h−1 g−1 soil) in the third and sixth week as an indicator for soil microbial activity. Bacterial and fungal abundance was estimated by qPCR. The proportional loss of litter (Medicago lupulina leaves) dry weight during soil incubation was used as an indication of decomposition rate. We measured the activity of β-glucosidase (cellulose degradation), β-D-celluliosidase (cellulose degradation), N-acetyl-β-glucosaminidase (chitin degradation) and phosphatase (organic phosphorus mineralization) using high throughput microplate assay43,44. A modified protocol by Kemper and Rosenau was used to measure water-stable soil aggregates45. Soil water repellency was measured using the water drop penetration time method46. High throughput sequencing (Illumina MiSeq) was used to measure the taxonomic composition of soil fungal and bacterial communities with the primers fITS7 and ITS4 for fungi and 515F-Y and 806 R for bacteria47,48 (Supplementary Methods).Statistical analysesFor diversity and community composition analysis, we excluded the samples with less than 1% of the observations of the largest sample in the ASV table. For network analysis, we then removed ASVs with low prevalence, which presented less than 20% of samples across all experimental units to reduce the high percentage of zero counts. A co-occurrence network was constructed based on both fungal and bacterial ASV tables. The PLNnetwork function in the R package PLNmodels was employed to infer the network, using a sparse multivariate Poisson log-normal (PLN) model49. According to the Extended Bayesian Information Criterion (EBIC), the best model was extracted with the function getBestModel. The network was compartmentalized into different modules using the cluster_fast_greedy function in the igraph package and visualized with partial correlations with |ρ| > 0.05. We focused on the response of the relative abundance of modules, also known as clusters, which represent the closely associated microbes, e.g., groups of coexisting or co-evolving microbes27. The relative abundance of modules was calculated by summing relative abundances for individual ASVs in modules. We used the package FUNGuildR50 to taxonomically parse fungal trait information, using the FUNGuild database51.For each single GCF treatment, we took the average response from the 10 replicates before analysis. To confirm how the effect of soil biodiversity treatment can change along with the increasing number of GCFs, we quantified the effect size of soil biodiversity treatment for each response variable using Hedges’ g (mean and 95% CIs) at each level of the number of GCFs, using the R package effsize52. Hedges’ g is calculated as the mean difference between the high and low soil biodiversity treatments in units of the pooled standard deviation as a paired-samples because there were identical GCFs and GCF combinations for both high and low biodiversity conditions, with the exception of the zero and 10 GCF treatments.To evaluate how each of the response variables changes along with the number of GCFs, we applied a generalized additive model (GAM)53. GAM is a penalized generalized linear model that fits a nonparametric, nonlinear smooth curve54. The degree of smoothness of model terms is estimated as part of fitting, using the generalized cross validation. We reasonably assume that the curve shapes are different between high and low soil biodiversity treatments. Therefore, we included biodiversity conditions (low/high) as the model intercept and as the “factor smooth” smoothing class, where a smooth function is created for each factor level independently55. For GAM modeling, we used the mgcv package55. The dimension of the basis used to represent the smooth term was set as k = 5 so that the model does not overfit to the data. For this, we compared some other values (from 3 to 8) and confirmed that the results are essentially the same within the tested range. The other parameters were set as default.The relationships between soil microbial indices and other soil indices were tested using Spearman correlation in the package microbiome, and the adjustment method “fdr” was employed to control the false discovery rate for multiple testing correction56. For the further multivariate integration of soil functions/properties and composition of modules, the DIABLO (Data Integration Analysis for Biomarker discovery using a Latent component method for Omics studies) was employed to detect correlation (Pearson’s correlation |r| > 0.5) among variables using the package mixOmics57.The Z-scores for each of the eight soil functions (as shown in Fig. 1, with the exception of soil water repellency) were evaluated, and then we computed an improved weighted multifunctionality metric to represent soil multifunctionality (Supplementary Methods)58. Structural equation models (SEMs) were used to reveal the direct and indirect effects of an increasing number of GCFs on soil multifunctionality within each soil biodiversity treatment using the package lavaan59. We assumed that an increasing number of GCFs influences soil multifunctionality by regulating the bacterial and fungal abundance and the relative abundance of modules. All response variables were standardized to the same comparison scale using the z-score transformation before constructing SEMs. Models with optimal fitting indices were reported (Supplementary Fig. 11).The permutational multivariate analysis (ADONIS) and non-metric multidimensional scaling (NMDS) ordination based on the Bray-Curtis distance were conducted to test the effect of soil biodiversity and GCF treatments on the community composition of bacteria and fungi using the R package vegan60. For the data handling, processing, and visualization, we used the packages tidyverse61, reshaping62, cowplot63, RColorBrewer64, qgraph65, igraph66, factoextra67, phyloseq68 and itsadug69. These data manipulation and analyses were conducted using R version 4.1.370. The R script and data are available in a publicly accessible database71.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

The pre-evaporitic, evaporitic, and post-evaporitic phases are recognized for the late Aptian. These phases are recorded within the K40–K50 sequences (Fig. 2A), and show an average maximum thickness of approximately 650 m in the studied basins. The pre-evaporitic phase is represented by carbonate and siliciclastic deposits formed in fluvial and lacustrine deltaic environments within a large proto-oceanic gulf28 (Fig. 2A). The peak of the evaporitic deposition is recorded in the K50 sequence, with widespread occurrences in the Brazilian equatorial margin. The origin of these deposits is the heat intensification associated with the widening of the Atlantic Ocean. These conditions caused strong evaporation leading to a wide distribution of evaporites (mainly halite and anhydrite gypsum) in the South Atlantic basins. The eastern continental margin of Brazil contains a restricted marine section characterized by evaporites, which are particularly prominent in thickness and occurrence in the Espírito Santo Basin (Itaúnas Member of the Mariricu Formation) and the Sergipe Basin (the Ibura Member of the Muribeca Formation)28. Evaporites form the most prominent evidence of dry climates in the South Atlantic basins11, with evaporation exceeding precipitation. The post-evaporitic phase is characterized by fully marine conditions evidenced by rich assemblages of marine fossils. During this phase, carbonates were deposited, followed by muddy and sandy sediments in shallow-marine and slope environments.Figure 2Paleoclimatic phases scheme and principal component analysis for paleoclimatic phases. (A) Paleoclimatic phases scheme for the late Aptian and the main depositional environments. (B) Principal component plot of bioclimatic groups. (C) Principal component for the pre-evaporitic phase (N = 92), evaporitic phase (N = 78), and post-evaporitic phase (N = 385); see Supplementary Fig. 9 for individual basins.Full size imagePaleovegetationWe identified a rich plant community with 139 spore and pollen genera/morphotypes representing all plant groups: bryophytes (five genera), ferns (58 genera), lycophytes (18 genera), pteridosperms (one genus), gymnosperms (27 genera), and angiosperms (30 genera) (Supplementary Table 2). The inferred systematic affinities at the family level reached 100% in bryophytes, 56.9% in ferns, 100% in lycophytes, 100% in pteridosperms, 92.6% in gymnosperms, and 40.0% in angiosperms, totaling 67.6% of the recorded genera (Supplementary Table 2). Marine elements (e.g., dinoflagellate cysts and microforaminiferal linings) were identified, in particular from the Sergipe and Araripe basins (Fig. 1). Pollen grains from gymnosperms were most abundant, represented mainly by the conifer families Cheirolepidiaceae, Araucariaceae, and Podocarpaceae, although representing different climatic settings. Classopollis (Cheirolepidiaceae) is the most abundant genus in all sections studied, followed by Araucariacites (Araucariaceae). Gymnosperms showed low diversity. Spore-producing plants are the most diverse in the assemblages of all basins (82 genera) and represented by several families of bryophytes, ferns, and lycophytes (e.g., Sphagnaceae, Anemiaceae, Cyatheaceae, Marsileaceae, Selaginellaceae, and Lycopodiaceae). These plant groups depend on water to reproduce and are therefore associated with humid settings.Cicatricosisporites (Anemiaceae) is the third most abundant palynomorph in all the basins, but especially in the northeastern basins (e.g., Sergipe Basin). Angiosperms are among the least abundant; however, they are diverse and include the most abundant and controversial genus Afropollis, herein attributed to angiosperms. In the most recent publication that addressed this question, ref.29 suggest that Afropollis should be treated as an angiosperm genus, although without more precise systematic assignment. The 30 genera/morphotypes of angiosperms are assigned to 8 families, viz., Arecaceae, Chloranthaceae, Euphorbiaceae, Flacourtiaceae, Illiciaceae, Liliaceae, Solanaceae and Trimeniaceae. The second most abundant genus is Stellatopollis also without precise systematic assignment.Spatio-temporal distribution of bioclimatic groupsOn the basis of their botanical affinities, most taxa were classified into five bioclimatic groups [see “Methods” section and Supplementary information], viz., hydrophytes, hygrophytes, tropical lowland flora, upland flora, and xerophytes (Supplementary Table 2) (Fig. 3).Figure 3Relevant palynomorphs of bioclimatic groups: (1) Aequitriradites sp.; (2) Crybelosporites sp.; (3) Perotriletes sp.; (4) Cicatricosisporites sp.; (5) Echinatisporis sp.; (6) Verrucosisporites sp.; (7) Bennettitaepollenites sp.; (8) Stellatopollis sp.; (9) Afropollis sp.; (10) Dejaxpollenites microfoveolatus; (11) Classopollis classoides; (12) Equisetosporites ovatus; (13) Gnetaceaepollenites jansonii; (14) Regalipollenites sp.; (15) Araucariacites sp.; (16) Callialasporites dampieri; (17) Complicatissacus cearensis; (18) Cyathidites sp.. Scale bar 20 µm.Full size imageOverall, the vegetation is dominated by the xerophytic bioclimatic group on account of the very high abundance of Classopollis (Cheirolepidiaceae) (general mean of 60.5%). However, the stratigraphic distribution of the bioclimatic groups in the sections studied (Supplementary Figs. 1–6) indicates wet phases confirmed by the curves of the other bioclimatic groups (hygrophytes, hydrophytes, tropical lowland flora, and upland flora). We used Pearson correlation analysis (Supplementary Fig. 7) to assess the correlation between the bioclimatic groups. The analysis revealed positive correlations between the bioclimatic groups of hygrophytes, hydrophytes, tropical lowland flora, and upland flora, and a negative correlation between these groups and the xerophyte group (Supplementary Fig. 7). The positive correlation between upland flora and hygrophytes confirms previous studies for the Sergipe Basin6,7, suggesting a relation between these groups and the hot and humid climate. The weak negative correlation between tropical lowland flora and upland flora is presumably related to elevation.The upland flora forms the second most abundant bioclimatic group, with an average of 18.9%. The large number of specimens of Araucariacites (Araucariaceae) in this group is notable. The hydrophytes are the least abundant group, with an average of only 1.4%. In this group, the highest values are attributed to the genus Crybelosporites (Marsileaceae).Principal component analyses (PCA) were used to reduce the multidimensional dataset, based on the percent abundance of the bioclimatic groups to a smaller number of dimensions for interpretive analysis. For all sections, two components or axes explain 97.6% of the observed variability (Fig. 2B). Hygrophytes, hydrophytes, tropical lowland flora, and upland flora show positive correlation (positive loading, 0.320, 0.029, 0.006, and 0.468, respectively), whereas xerophytes show a negative relationship (negative loading, − 0.823) on the first axis, which alone explains 83.0% of the variability. In summary, the first axis of the PCA reveals a separation of two major climatic conditions (wet and dry) along the axis (Fig. 2B). The wet conditions include the associations of hygrophytes, hydrophytes, tropical lowland flora, and upland flora, with dry conditions associated with taxa from the xerophyte group. The second axis explains 14.6%, in which hygrophytes, hydrophytes, and tropical lowland flora show a positive correlation relationship (positive loading, 0.719, 0.037, 0.036, respectively), whereas upland flora and xerophytes show a negative relationship (negative loading, − 0.684 and − 0.108, respectively). With respect to the second axis, a polarization between the hygrophytes (positive loading, 0.719) and the upland flora (negative loading, − 0.684) can be interpreted as a lowland–upland trend. The same pattern was recorded for all paleoclimatic phases (Fig. 2C) and sections (Supplementary Fig. 8), that is, the first axis is related to humidity vs. aridity, and the second axis to elevation (lowland vs. upland). This suggests that these two factors, particularly the first one, controlled the vegetation distribution in the late Aptian of the region. As all bioclimatic groups occurred in the three evaporitic phases, these trends in abundance reflect expansion and contraction of the recorded vegetation.Parallel increasing trends of bioclimatic groups mark the pre-evaporitic phase: hygrophytes and upland flora in the Bragança-Viseu, São Luís, Parnaíba, Ceará, Potiguar, and Araripe basins (Supplementary Figs. 1–3 and 5), suggesting that there was a certain amount of moisture in these areas. The xerophytes show the lowest average of this phase (44.1%) (Table 1), whereas hygrophytes show the highest average (27.0%). These humid conditions are confirmed by the highest mean of the Fs/X ratio (Fs/X = 0.4), representing the predominance of spore-producing plants [see Methods section and Supplementary information]. Despite the low abundance of hydrophytes in the sections, a prominent feature is the highest average (2.5%) of this group (Table 1), which is assigned to aquatic environments, confirming relatively wet conditions in this phase. There are no pre-evaporitic samples available from the Sergipe and Espírito Santo basins.Table 1 Average abundance of bioclimatic groups, diversity, Fs/X and marine elements for the paleoclimatic phases.Full size tableThe evaporitic phase is characterized by the highest abundance of the xerophyte bioclimatic group (76.4%) (Table 1), represented mainly by Classopollis (Supplementary Figs. 1–6). A high abundance of xerophytes occurred widely distributed in all basins studied. In this phase, tropical lowland flora is notable, showing an average higher than the overall average (3.3%), particularly in the Bragança-Viseu, São Luís, Parnaíba, and Ceará basins (Supplementary Figs. 1 and 2). This result is related to the moderate to high abundance of the genus Afropollis in these basins. The evaporitic phase is also characterized by the lowest average Fs/X ratio (Fs/X = 0.1) (Table 1), confirming the dominance of xerophytes.The post-evaporitic phase is characterized by the upland flora bioclimatic group (mean = 24.4%) (Table 1). The moderate to high abundance of upland flora in this phase is represented, in particular, by pollen grains of Araucariacites, which represent the high-relief family Araucariaceae. This bioclimatic group is associated with more humid conditions, as confirmed by an Fs/X ratio higher than the overall average (Fs/X = 0.2). The upland flora is significant in all basins, except the Espírito Santo Basin, where xerophytes predominate in both studied phases in this basin.Latitudinal biome distributionsBiome change is a fundamental biological response to climate change. In the study area, the predominance of a specific biome is mainly related to humidity, since all five recorded bioclimatic groups are related to a warm climate (Supplementary Table 2) representing two biomes: tropical xerophytic shrubland and tropical rainforest. In the rainforest biome two phytophysiognomies are recognized: lowland and montane rainforest. The tropical xerophytic shrubland biome predominates in the three paleoclimatic phases, with a wide latitudinal range from the Bragança-Viseu, São Luís, and Parnaíba basins (1° S) to the Espírito Santo Basin (20° S). This wide distribution is compatible with a predominantly arid climate in South America in the late Aptian, as indicated by paleoclimatic maps8,9,15 (Fig. 4A). Most arid and semi-arid ecosystems are mainly controlled by precipitation. Other climate parameters are less important, a condition that simplifies cause-effect interpretations. The PCA (Fig. 2B) demonstrated that the wet–dry trend, which reflects high–low precipitation, was the main determinant in the distribution of the biomes. However, considering all phases, an increasing trend in humidity was observed from the southeast (Espírito Santo Basin) to the northeast (e.g., Potiguar Basin) (Fig. 4B), coinciding with the hot and wet belt attributed to the ITCZ (Fig. 4A)15. The latitudinal distribution of diversity also follows this trend. Diversity increased significantly towards in the basins near the equator. Diversity indices (Shannon – H’) peaked in the Sergipe Basin (H’ = 3.5, CL-47 section) at 11° S. Conversely, the lowest average diversity is recorded in the Espírito Santo Basin (H’ = 1.1) at 20° S. Additionally, there is a clear correlation between high diversity (H’) and humidity (Fs/X ratio) (r = 0.691), regardless of paleoclimatic phase, as evidenced by the synchronicity of the H’ and Fs/X curves (Fig. 5). After data normalization between humidity (Fs/X) and marine elements (dinoflagellate cysts and microforaminifer linings), we performed linear correlation analyses, which showed a weak but positive correlation (r = 0.137). This is due to the fact that pre- evaporitic deposits contain only 19 occurrences of dinoflagellate cysts in 90 samples. Despite this, the curves of Fs/X, marine elements and diversity are synchronous (Fig. 5), suggesting a relation between humidity, diversity, and marine incursions.Figure 4Latitudinal changes in late Aptian biomes from southeast to center-north. (A) Paleoclimatic belts of the late Aptian in South America (climatic belts modified from refer.14). Reconstruction map at 116 Ma modified from ODSN Plate Tectonic Reconstruction Service. The Reconstruction map at 116 Ma was generated by ODSN Plate Tectonic Reconstruction Service (https://www.odsn.de/odsn/services/paleomap/paleomap.html). (B) Late Aptian latitudinal distribution of the tropical xerophytic biome in Brazil. (C) Stratigraphic distribution of biomes for individual basins. (D) Relative Importance of biomes for paleoclimatic phases.Full size imageFigure 5Biome trends in relation to paleoclimatic phases. Change in biomes, diversity, Fs/X ratio and marine elements shown by changepoint analysis plotted against paleoclimatic phases.Full size imageThe pre-evaporitic phase is marked by a certain balance between the biomes (Fig. 4C,D). In the lowlands, the tropical xerophytic shrubland biome predominated in the Bragança- Viseu, São Luís, Parnaíba, and Ceará basins, but in the Potiguar Basin it is co-dominant with the lowland rainforest. The montane rainforest was relatively extensive in this phase, although with several areal changes, and reached its widest extent in the Araripe (7° S) and Potiguar (5° S) basins in response to the deterioration of the tropical xerophytic shrubland biome. These conditions demonstrate that humidity was relatively high at this stage. The pre-evaporitic deposits were characterized by the highest diversity average (H’ = 1.8).The method of indicator species analysis (IndVal) was used to identify the key species of each paleoclimatic phase (Supplementary Table 15). The species identified for the pre-evaporitic phase, Deltoidospora spp. (Cyatheaceae-Dicksoniaceae) related to the montane rainforest, are indicator species for the Bragança-Viseu, São Luís, Parnaíba, and Ceará basins. The Gnetaceaepollenites spp. (Gnetaceae) of the Potiguar Basin and Equisetosporites spp. (Ephedraceae) of the Araripe Basin are related to the tropical xerophytic shrubland biome (Supplementary Table 15). Even for the pre-evaporitic phase, a progressive increase in the tropical xerophytic shrubland biome was observed and interpreted as the start of a climatic deterioration stage (Fig. 4C), which culminated in the evaporitic phase. Shifts in vegetation types may occur when precipitation reaches a threshold value, which means that a regionally synchronous gradual climate change can cause abrupt vegetation shifts. The change from humid to warm and arid conditions (evaporitic phase) is directly related to a decrease in precipitation. This aridization process coincides with the appearance of marine elements (e.g., dinoflagellate cysts). The threshold effect (intense evaporation) is reflected in an abrupt decrease in the abundance of lowland and montane rainforest and a sharp increase to a very high abundance of the tropical xerophytic shrubland biome (Supplementary Figs. 4C and 5). The threshold effect was not detected in the Espírito Santo Basin, where the arid conditions remained stable with minimal shift (expansion and contraction) of the biome. The main representatives of this biome are conifers of the family Cheirolepidiaceae (Classopollis), which were most abundant in lagoons and coastal environments and are often associated with evaporates30,31,32,33,34,35. Even under xeric or water-stressed conditions there was a slight increase in biomes related to a humid climate (lowland and montane rainforest phytophysiognomies) towards the equatorial region, suggesting influence of the ITCZ (Fig. 4A,B).The evaporitic phase was characterized by the lowest diversity average (H’ = 1.2). With modest rainfall, arid regions are generally characterized by fewer species than moister biomes36. However, diversity indices peaked in the Bragança-Viseu, São Luís, and Parnaíba basins (H’ = 2.6, RL-01 section) and along the equatorial margin (2° S) (Supplementary Fig. 1).IndVal emphasizes the xeric conditions in the evaporitic phase by association with the species from the tropical xerophytic shrubland biome: Classopollis spp. (Ceará and Potiguar basins), Classopollis classoides (Sergipe Basin), Classopollis intrareticulatus (Araripe Basin), and Gnetaceaepollenites spp. (Espírito Santo Basin). For the Bragança-Viseu, São Luís, Parnaíba, and Ceará basins, where xeric restrictions are milder, the indicator taxon is Afropollis spp. from the lowland rainforest. This genus shows the weakest negative correlation with xerophytes.After the end of evaporite deposition, all sections indicate climatic stability, which kept the climate hot and arid even in the post-evaporitic phase, although the response was not linear.The shift in the biomes, especially the tropical xerophytic shrubland in the Bragança-Viseu, São Luís, Parnaíba, Ceará, and Araripe basins, occurred in the transition between the evaporitic and post-evaporitic phases, whereas in the Potiguar and Sergipe basins it occurred within the post-evaporitic phase. As indicated in the dendrograms of each section (Supplementary Figs. 1–6), the shift occurred abruptly in all basins, except the Espírito Santo Basin. The tropical rainforest biome (lowland and montane rainforests) replaced the tropical xerophytic shrubland in almost all basins (Fig. 4C). Even the Espírito Santo Basin, far from the influence of the ITCZ, shows a slight increase in lowland rainforest. The changes in the biomes are attributable to threshold effects caused by gradual climate change related to the ITCZ intensification shift and progressive increase in marine influence, indicated by an increase in marine microplankton from an average of 3.9% in the evaporitic phase to 44.1%. The increase in marine influence is reflected in the first major flooding surface observed in the Cretaceous succession27. Thus, a climate amelioration stage was established in the post- evaporitic phase (Fig. 5). In combination with published paleotopographic information25, the bioclimatic groups associated to the humid conditions (hygrophytes, hydrophytes, tropical lowland flora, and upland flora) were combined and visualized to create Fig. 6.Figure 6Reconstruction of the transitional gradient between marine to terrestrial environment (uplands) under ITCZ influence. The illustration is based on paleoflora and environmental information from palynological data from studied sections. Original size illustration: 18 × 24 cm, by Julio Lacerda.Full size imageAccording to refs.7,37, arid conditions are characterized by sea-level lowstands, whereas warm and humid conditions are correlated with sea levels rise, which explains the increase in the tropical rainforest biome (lowland and montane rainforests). The more intense humidity is supported by the results of IndVal for the post-evaporitic phase, with all species related to humid climate: Deltoidospora spp. (Bragança-Viseu, São Luís and Parnaíba basins), Araucariacites limbatus (Ceará Basin), Cicatricosisporites spp. (Potiguar Basin), Cicatricosisporites spp. and Araucariacites australis (Sergipe Basin), Inaperturopollenites spp. (Araripe Basin) and Inaperturopollenites simplex (Espírito Santo Basin).Our results show that the ITCZ combined with the opening of the South Atlantic Ocean during the late Aptian altered vegetation dynamics. As today, the ITCZ influence is stronger in the northeastern and north-central regions of South America. It is notable that the late Aptian climate evolution in the South Atlantic, culminating in higher humidity, was accompanied by an intrinsic relation between plant diversity, humidity, and marine influence. More

Chen, L. et al. Geographic variation in traits of fruit stones and seeds of Melia azedarach. J. Beijing For. Univ. 36, 15–20 (2014).CAS 

Google Scholar 
Angamuthu, D., Purushothaman, I., Kothandan, S. & Swaminathan, R. Antiviral study on Punica granatum L., Momordica charantia L., Andrographis paniculata Nees, and Melia azedarach L., to human herpes virus-3. Eur. J. Integr. Med. 28, 98–108. https://doi.org/10.1016/j.eujim.2019.04.008 (2019).Article 

Google Scholar 
Wang, N. et al. Selective ERK1/2 agonists isolated from Melia azedarach with potent anti-leukemic activity. BMC Cancer 19, 1–9. https://doi.org/10.1186/s12885-019-5914-8 (2019).CAS 
Article 

Google Scholar 
Khoshraftar, Z., Safekordi, A., Shamel, A. & Zaefizadeh, M. Evaluation of insecticidal activity of nanoformulation of Melia azedarach (leaf) extract as a safe environmental insecticide. Int. J. Environ. Sci. Technol. 17, 1159–1170. https://doi.org/10.1007/s13762-019-02448-7 (2020).CAS 
Article 

Google Scholar 
Sivaraj, I., Nithaniyal, S., Bhooma, V., Senthilkumar, U. & Parani, M. Species delimitation of Melia dubia Cav. from Melia azedarach L. complex based on DNA barcoding. Botany 96, 329–336. https://doi.org/10.1139/cjb-2017-0148 (2018).CAS 
Article 

Google Scholar 
Liao, B. et al. Population structure and genetic relationships of Melia Taxa in China assayed with sequence-related amplified polymorphism (SRAP) markers. Forests 7, 81. https://doi.org/10.3390/f7040081 (2016).Article 

Google Scholar 
Wu, L., Kaewmano, A., Fu, P., Wang, W. & Fan, Z. Intra-annual radial growth of Melia azedarach in a tropical moist seasonal forest and its response to environmental factors in Xishuangbanna Southwest China. Acta Ecol. Sin. 40, 6831–6840. https://doi.org/10.5846/stxb202003120508 (2020).Article 

Google Scholar 
Hoegh-Guldberg, O. et al. The human imperative of stabilizing global climate change at 1.5 C. Science 365, eaaw6974. https://doi.org/10.1126/science.aaw6974 (2019).CAS 
Article 
PubMed 

Google Scholar 
López-Tirado, J., Vessella, F., Schirone, B. & Hidalgo, P. J. Trends in evergreen oak suitability from assembled species distribution models: Assessing climate change in south-western Europe. New For. 49, 471–487. https://doi.org/10.1007/s11056-018-9629-5 (2018).Article 

Google Scholar 
Xu, Y. et al. Modelling the effects of climate change on the distribution of endangered Cypripedium japonicum in China. Forests 12, 429. https://doi.org/10.3390/f12040429 (2021).Article 

Google Scholar 
Booth, T. H. Species distribution modelling tools and databases to assist managing forests under climate change. For. Ecol. Manag. 430, 196–203. https://doi.org/10.1016/j.foreco.2018.08.019 (2018).Article 

Google Scholar 
Dyderski, M. K., Paź, S., Frelich, L. E. & Jagodziński, A. M. How much does climate change threaten European forest tree species distributions?. Glob. Change Biol. 24, 1150–1163. https://doi.org/10.1111/gcb.13925 (2018).ADS 
Article 

Google Scholar 
Zhong, Y. et al. A generalized linear mixed model approach to assess emerald ash Borer diffusion. ISPRS Int. J. Geo Inf. 9, 414. https://doi.org/10.3390/ijgi9070414 (2020).Article 

Google Scholar 
Chang, Z., Meng, J., Shi, Y. & Mo, F. Lnc RNA recognition by fusing multiple features and its function prediction. CAAI Trans. Intell. Syst. 13, 928–934. https://doi.org/10.11992/tis.201806008 (2018).Article 

Google Scholar 
Shiferaw, H., Bewket, W. & Eckert, S. Performances of machine learning algorithms for mapping fractional cover of an invasive plant species in a dryland ecosystem. Ecol. Evol. 9, 2562–2574. https://doi.org/10.1002/ece3.4919 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Tang, X., Yuan, Y., Li, X. & Zhang, J. Maximum entropy modeling to predict the impact of climate change on pine wilt disease in China. Front. Plant Sci. 12, 764. https://doi.org/10.3389/fpls.2021.652500 (2021).Article 

Google Scholar 
Chhogyel, N., Kumar, L., Bajgai, Y. & Jayasinghe, L. S. Prediction of Bhutan’s ecological distribution of rice (Oryza sativa L.) under the impact of climate change through maximum entropy modelling. J. Agric. Sci. 158, 25–37. https://doi.org/10.1017/S0021859620000350 (2020).Article 

Google Scholar 
Ahmad, Z. et al. Melia Azedarach impregnated Co and Ni zero-valent metal nanoparticles for organic pollutants degradation: Validation of experiments through statistical analysis. J. Mater. Sci. Mater. Electron. 31, 16938–16950. https://doi.org/10.1007/s10854-020-04250-5 (2020).CAS 
Article 

Google Scholar 
Hijmans, R. J., Huaccho, L. & Zhang, D. In I International Conference on Sweetpotato. Food and Health for the Future 583, 41–49.Luo, M., Wang, H. & Lyu, Z. Evaluating the performance of species distribution models Biomod2 and MaxEnt using the giant panda distribution data. J. Appl. Ecol. 28, 4001–4006. https://doi.org/10.13287/j.1001-9332.201712.011 (2017).Article 

Google Scholar 
Wang, T., Wang, G., Innes, J. L., Seely, B. & Chen, B. ClimateAP: An application for dynamic local downscaling of historical and future climate data in Asia Pacific. Front. Agric. Sci. Eng. 4, 448–458. https://doi.org/10.15302/J-FASE-2017172 (2017).CAS 
Article 

Google Scholar 
Yang, X.-Q., Kushwaha, S., Saran, S., Xu, J. & Roy, P. Maxent modeling for predicting the potential distribution of medicinal plant, Justicia adhatoda L Lesser Himalayan foothills. Ecol. Eng. 51, 83–87. https://doi.org/10.1016/j.ecoleng.2012.12.004 (2013).CAS 
Article 

Google Scholar 
Pepe, M. S., Cai, T. & Longton, G. Combining predictors for classification using the area under the receiver operating characteristic curve. Biometrics 62, 221–229. https://doi.org/10.1111/j.1541-0420.2005.00420.x (2006).MathSciNet 
Article 
PubMed 
MATH 

Google Scholar 
McHugh, M. L. Interrater reliability: The kappa statistic. Biochem. Med. 22, 276–282. https://hrcak.srce.hr/89395 (2012).Article 

Google Scholar 
Lu, C. Y., Gu, W., Dai, A. H. & Wei, H. Y. Assessing habitat suitability based on geographic information system (GIS) and fuzzy: A case study of Schisandra sphenanthera Rehd. et Wils. in Qinling Mountains, China. Ecol. Model. 242, 105–115. https://doi.org/10.1016/j.ecolmodel.2012.06.002 (2012).Article 

Google Scholar 
Zhang, L. et al. The basic principle of random forest and its applications in ecology: A case study of Pinus yunnanensis. Acta Ecol. Sin. 34, 650–659. https://doi.org/10.5846/stxb201306031292 (2014).Article 

Google Scholar 
Williams, J. N. et al. Using species distribution models to predict new occurrences for rare plants. Divers. Distrib. 15, 565–576. https://doi.org/10.1111/j.1472-4642.2009.00567.x (2009).Article 

Google Scholar 
Akpoti, K., Kabo-Bah, A. T., Dossou-Yovo, E. R., Groen, T. A. & Zwart, S. J. Mapping suitability for rice production in inland valley landscapes in Benin and Togo using environmental niche modeling. Sci. Total Environ. 709, 136165. https://doi.org/10.1016/j.scitotenv.2019.136165 (2020).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Dutra Silva, L., de Brito, A. E., Vieira Reis, F., Bento Elias, R. & Silva, L. Limitations of species distribution models based on available climate change data: a case study in the Azorean forest. Forests 10, 575. https://doi.org/10.3390/f10070575 (2019).Article 

Google Scholar 
Lin, H. Y. et al. Climate-based approach for modeling the distribution of montane forest vegetation in Taiwan. Appl. Veg. Sci. 23, 239–253. https://doi.org/10.1111/avsc.12485 (2020).Article 

Google Scholar 
Zhang, L. et al. Consensus forecasting of species distributions: The effects of niche model performance and niche properties. PLoS ONE 10, e0120056. https://doi.org/10.1371/journal.pone.0120056 (2015).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Zhang, H. The optimality of naive Bayes. Am. Assoc. Artif. Intell. 1, 3 (2004).
Google Scholar 
Wang, Q., Nguyen, T.-T., Huang, J. Z. & Nguyen, T. T. An efficient random forests algorithm for high dimensional data classification. Adv. Data Anal. Classif. 12, 953–972. https://doi.org/10.1007/s11634-018-0318-1 (2018).MathSciNet 
Article 
MATH 

Google Scholar 
Zheng-tao, Y., Bin, D., Bo, H., Lu, H. & Jian-yi, G. Word sense disambiguation based on bayes model and information gain. Proc. Int. J. Adv. Sci. Technol. 2, 153–157. https://doi.org/10.1109/FGCN.2008.188 (2009).Article 

Google Scholar 
Yu, B. et al. SubMito-XGBoost: Predicting protein submitochondrial localization by fusing multiple feature information and eXtreme gradient boosting. Bioinformatics 36, 1074–1081. https://doi.org/10.1093/bioinformatics/btz734 (2020).CAS 
Article 
PubMed 

Google Scholar 
Hailu, B. T., Siljander, M., Maeda, E. E. & Pellikka, P. Assessing spatial distribution of Coffea arabica L. in Ethiopia’s highlands using species distribution models and geospatial analysis methods. Ecol. Inf. 42, 79–89. https://doi.org/10.1016/j.ecoinf.2017.10.001 (2017).Article 

Google Scholar 
Ramirez-Reyes, C. et al. Embracing ensemble species distribution models to inform at-risk species status assessments. J. Fish Wildl. Manag. 12, 98–111. https://doi.org/10.3996/JFWM-20-072 (2021).Article 

Google Scholar 
Wisz, M. S. et al. Effects of sample size on the performance of species distribution models. Divers. Distrib. 14, 763–773. https://doi.org/10.1111/j.1472-4642.2008.00482.x (2008).Article 

Google Scholar 
Feng, L., Sun, J., Shi, Y., Wang, G. & Wang, T. Predicting suitable habitats of camptotheca acuminata considering both climatic and soil variables. Forests 11, 891. https://doi.org/10.3390/f11080891 (2020).Article 

Google Scholar 
Wang, T., Campbell, E. M., O’Neill, G. A. & Aitken, S. N. Projecting future distributions of ecosystem climate niches: Uncertainties and management applications. For. Ecol. Manag. 279, 128–140. https://doi.org/10.1016/j.foreco.2012.05.034 (2012).Article 

Google Scholar 
Wang, T., Hamann, A., Spittlehouse, D. L. & Murdock, T. Q. ClimateWNA—high-resolution spatial climate data for western North America. J. Appl. Meteorol. Climatol. 51, 16–29. https://doi.org/10.1175/JAMC-D-11-043.1 (2012).ADS 
Article 

Google Scholar 
Feng, L. et al. Predicting suitable habitats of ginkgo biloba L. fruit forests in China. Clim. Risk Manag. 34, 100364. https://doi.org/10.1016/j.crm.2021.100364 (2021).Article 

Google Scholar 
Wang, T., Hamann, A., Spittlehouse, D. & Carroll, C. Locally downscaled and spatially customizable climate data for historical and future periods for North America. PLoS ONE 11, e0156720. https://doi.org/10.1371/journal.pone.0156720 (2016).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Guo, Y. et al. Spatial prediction and delineation of Ginkgo biloba production areas under current and future climatic conditions. Ind. Crops Prod. 166, 113444. https://doi.org/10.1016/j.indcrop.2021.113444 (2021).Article 

Google Scholar 
Jiao, C., Lan, G., Sun, Y., Wang, G. & Sun, Y. Dopamine alleviates chilling stress in watermelon seedlings via modulation of proline content, antioxidant enzyme activity, and polyamine metabolism. J. Plant Growth Regul. 40, 2. https://doi.org/10.1007/s00344-020-10096-2 (2021).CAS 
Article 

Google Scholar 
Thakur, S., Thakur, I. & Sankanur, M. Assessment of genetic diversity in drek (Melia azedarach) using molecular markers. J. Tree Sci. 36, 78–85. https://doi.org/10.5958/2455-7129.2017.00011.5 (2017).Article 

Google Scholar 
Sivasubramaniam, K. et al. Seed priming: Triumphs and tribulations. The Madras Agricultural Journal 98, 197–209. https://www.researchgate.net/publication/267298497 (2011).
Google Scholar 
Xu, L. et al. Effect of salt stress on growth and physiology in Melia azedarach seedlings of six provenances. Int. J. Agric. Biol. 20, 471–480. https://doi.org/10.17957/IJAB/15.0618 (2018).CAS 
Article 

Google Scholar 
Lenoir, J., Gégout, J.-C., Marquet, P., De Ruffray, P. & Brisse, H. A significant upward shift in plant species optimum elevation during the 20th century. Science. 320, 1768–1771. https://doi.org/10.1126/science.1156831 (2008).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Ou-Yang, C.-F. et al. Impact of equatorial and continental airflow on primary greenhouse gases in the northern South China Sea. Environ. Res. Lett. 10, 065005. https://doi.org/10.1088/1748-9326/10/6/065005 (2015).ADS 
CAS 
Article 

Google Scholar 
Liu, B., Zhu, C., Su, J., Ma, S. & Xu, K. Record-breaking northward shift of the western North Pacific subtropical high in July 2018. J. Meteorol. Soc. Japan. 97, 913–925. https://doi.org/10.2151/jmsj.2019-047 (2019).ADS 
Article 

Google Scholar 
Huang, J. et al. Dryland climate change: Recent progress and challenges. Rev. Geophys. 55, 719–778. https://doi.org/10.1002/2016RG000550 (2017).ADS 
Article 

Google Scholar 
Thuiller, W., Lavorel, S., Araújo, M. B., Sykes, M. T. & Prentice, I. C. Climate change threats to plant diversity in Europe. Proc. Natl. Acad. Sci. 102, 8245–8250. https://doi.org/10.1073/pnas.0409902102 (2005).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Waldvogel, A. M. et al. Evolutionary genomics can improve prediction of species’ responses to climate change. Evol. Lett. 4, 4–18. https://doi.org/10.1002/evl3.154 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Vilà-Cabrera, A., Coll, L., Martínez-Vilalta, J. & Retana, J. Forest management for adaptation to climate change in the Mediterranean basin: A synthesis of evidence. For. Ecol. Manag. 407, 16–22. https://doi.org/10.1016/j.foreco.2017.10.021 (2018).Article 

Google Scholar 
He, X., Li, J., Wang, F., Zhang, J. & Chen, X. Variation and selection of Melia azedarach provenances and families. J. Northeast For. Univ. 47, 1–7. https://doi.org/10.13332/j.1000-1522.20170321 (2019).CAS 
Article 

Google Scholar 
Smith, A. B., Alsdurf, J., Knapp, M., Baer, S. G. & Johnson, L. C. Phenotypic distribution models corroborate species distribution models: A shift in the role and prevalence of a dominant prairie grass in response to climate change. Glob. Change Biol. 23, 4365–4375. https://doi.org/10.1111/gcb.13666 (2017).Article 

Google Scholar 
Bellon, M. R., Dulloo, E., Sardos, J., Thormann, I. & Burdon, J. J. In situ conservation—harnessing natural and human-derived evolutionary forces to ensure future crop adaptation. Evol. Appl. 10, 965–977. https://doi.org/10.1111/eva.12521 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Bidak, L. M., Heneidy, S. Z., Halmy, M. W. A. & El-Kenany, E. T. Sustainability potential for Ginkgo biloba L. plantations under climate change uncertainty: An ex-situ conservation perspective. Acta Ecol. Sin. 42, 101–114. https://doi.org/10.1016/j.chnaes.2021.09.012 (2021).Article 

Google Scholar 
Qin, F., Liu, S. & Yu, S. Effects of allelopathy and competition for water and nutrients on survival and growth of tree species in Eucalyptus urophylla plantations. For. Ecol. Manag. 424, 387–395. https://doi.org/10.1016/j.foreco.2018.05.017 (2018).Article 

Google Scholar 
Zabel, F. et al. Global impacts of future cropland expansion and intensification on agricultural markets and biodiversity. Nat. Commun. 10, 1–10. https://doi.org/10.1038/s41467-019-10775-z (2019).ADS 
CAS 
Article 

Google Scholar  More

Pimentel D, McNair S, Janecka J, Wightman J, Simmonds C, O’Connell C, et al. Economic and environmental threats of alien plant, animal, and microbe invasions. Agric Ecosyst Environ. 2001;84:1–20.Article 

Google Scholar 
Diagne C, Leroy B, Vaissière AC, Gozlan RE, Roiz D, Jarić I, et al. High and rising economic costs of biological invasions worldwide. Nature. 2021;592:571–6.Article 
CAS 

Google Scholar 
Catford JA, Jansson R, Nilsson C. Reducing redundancy in invasion ecology by integrating hypotheses into a single theoretical framework. Divers Distrib. 2009;15:22–40.Article 

Google Scholar 
Pearson DE, Ortega YK, Eren Ö, Hierro JL. Community assembly theory as a framework for biological invasions. Trends Ecol Evol. 2018;33:313–25.PubMed 
Article 

Google Scholar 
Inderjit, van der Putten WH. Impacts of soil microbial communities on exotic plant invasions. Trends Ecol Evol. 2010;25:512–9.PubMed 
Article 
CAS 

Google Scholar 
Keane RM, Crawley MJ. Exotic plant invasions and the enemy release hypothesis. Trends Ecol Evol. 2002;17:164–70.Article 

Google Scholar 
Stinson KA, Campbell SA, Powell JR, Wolfe BE, Callaway RM, Thelen GC, et al. Invasive plant suppresses the growth of native tree seedlings by disrupting belowground mutualisms. PLoS Biol. 2006;4:e140.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Hierro JL, Callaway RM. Allelopathy and exotic plant invasion. Plant Soil. 2003;256:29–39.Article 
CAS 

Google Scholar 
Reinhart KO, Callaway RM. Soil biota and invasive plants. N Phytol. 2006;170:445–57.Article 

Google Scholar 
Waller LP, Allen WJ, Barratt BIP, Condron LM, França FM, Hunt JE, et al. Biotic interactions drive ecosystem responses to exotic plant invaders. Science. 2020;368:967–72.PubMed 
Article 
CAS 

Google Scholar 
McLeod ML, Cleveland CC, Lekberg Y, Maron JL, Philippot L, Bru D, et al. Exotic invasive plants increase productivity, abundance of ammonia-oxidizing bacteria and nitrogen availability in intermountain grasslands. J Ecol. 2016;104:994–1002.Article 
CAS 

Google Scholar 
Saul WC, Jeschke JM. Eco-evolutionary experience in novel species interactions. Ecol Lett. 2015;18:236–45.PubMed 
Article 

Google Scholar 
Desprez-Loustau M, Robin C, Buee M, Courtecuisse R, Garbaye J, Suffert F, et al. The fungal dimension of biological invasions. Trends Ecol Evol. 2007;22:472–80.PubMed 
Article 

Google Scholar 
Hierro JL, Maron JL, Callaway RM. A biogeographical approach to plant invasions: the importance of studying exotics in their introduced and native range. J Ecol. 2005;93:5–15.Article 

Google Scholar 
Callaway RM, Thelen GC, Rodriguez A, Holben WE. Soil biota and exotic plant invasion. Nature. 2004;427:731–3.PubMed 
Article 
CAS 

Google Scholar 
Maron JL, Klironomos J, Waller L, Callaway RM. Invasive plants escape from suppressive soil biota at regional scales. J Ecol. 2014;102:19–27.Article 

Google Scholar 
Brundrett MC. Mycorrhizal associations and other means of nutrition of vascular plants: understanding the global diversity of host plants by resolving conflicting information and developing reliable means of diagnosis. Plant Soil. 2009;320:37–77.Article 
CAS 

Google Scholar 
Smith SE, Read DJ. Mycorrhizal symbiosis. London: Academic Press; 2008.O’Neill EG, O’Neill RV, Norby RJ. Hierarchy theory as a guide to mycorrhizal research on large-scale problems. Environ Pollut. 1991;73:271–84.PubMed 
Article 

Google Scholar 
Johnson NC, Wilson GWTT, Wilson JA, Miller RM, Bowker MA. Mycorrhizal phenotypes and the Law of the Minimum. N Phytol. 2015;205:1473–84.Article 
CAS 

Google Scholar 
Lekberg Y, Arnillas CA, Borer ET, Bullington LS, Fierer N, Kennedy PG, et al. Nitrogen and phosphorus fertilization consistently favor pathogenic over mutualistic fungi in grassland soils. Nat Commun. 2021;12:3484.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Richardson DM, Allsopp N, D’Antonio CM, Milton S, Rejmanek M. Plant invasions – the role of mutualisms. Biol Rev. 2000;75:65–93.PubMed 
Article 
CAS 

Google Scholar 
Marler MJ, Zabinski CA, Callaway RM. Mycorrhizae indirectly enhance competitive effects of an invasive forb on a native bunchgrass. Ecology. 1999;80:1180–6.Article 

Google Scholar 
Soti PG, Jayachandran K, Purcell M, Volin JC, Kitajima K. Mycorrhizal symbiosis and Lygodium microphyllum invasion in South Florida—a biogeographic comparison. Symbiosis. 2014;62:81–90.Article 

Google Scholar 
Fumanal B, Plenchette C, Chauvel B, Bretagnolle F. Which role can arbuscular mycorrhizal fungi play in the facilitation of Ambrosia artemisiifolia L. invasion in France? Mycorrhiza. 2006;17:25–35.PubMed 
Article 
CAS 

Google Scholar 
Hart MM, Reader RJ. Taxonomic basis for variation in the colonization strategy of arbuscular mycorrhizal fungi. N Phytol. 2002;153:335–44.Article 

Google Scholar 
Maherali H, Klironomos JN. Influence of phylogeny on fungal community assembly and ecosystem functioning. Science. 2007;316:1746–8.PubMed 
Article 
CAS 

Google Scholar 
Kivlin SN, Hawkes CV, Treseder KK. Global diversity and distribution of arbuscular mycorrhizal fungi. Soil Biol Biochem. 2011;43:2294–303.Article 
CAS 

Google Scholar 
Davison J, Moora M, Öpik M, Adholeya A, Ainsaar L, Bâ A, et al. Global assessment of arbuscular mycorrhizal fungus diversity reveals very low endemism. Science. 2015;349:970–3.PubMed 
Article 
CAS 

Google Scholar 
Mitchell CE, Agrawal AA, Bever JD, Gilbert GS, Hufbauer RA, Klironomos JN, et al. Biotic interactions and plant invasions. Ecol Lett. 2006;9:726–40.PubMed 
Article 

Google Scholar 
Ehrenfeld JG. Effects of exotic plant invasions on soil nutrient cycling processes. Ecosystems. 2003;6:503–23.Article 
CAS 

Google Scholar 
Rout ME, Chrzanowski TH. The invasive Sorghum halepense harbors endophytic N2-fixing bacteria and alters soil biogeochemistry. Plant Soil. 2009;315:163–72.Article 
CAS 

Google Scholar 
Sardans J, Bartrons M, Margalef O, Gargallo-Garriga A, Janssens IA, Ciais P, et al. Plant invasion is associated with higher plant-soil nutrient concentrations in nutrient-poor environments. Glob Change Biol. 2017;23:1282–91.Article 

Google Scholar 
Bossdorf O, Auge H, Lafuma L, Rogers WE, Siemann E, Prati D. Phenotypic and genetic differentiation between native and introduced plant populations. Oecologia. 2005;144:1–11.PubMed 
Article 

Google Scholar 
Lankau RA. Resistance and recovery of soil microbial communities in the face of Alliaria petiolata invasions. N Phytol. 2011;189:536–48.Article 

Google Scholar 
Blossey B, Nötzold R. Evolution of increased competitive ability in invasive nonindigenous plants: a hypothesis. J Ecol. 1995;83:887–9.Article 

Google Scholar 
van Kleunen M, Bossdorf O, Dawson W. The ecology and evolution of alien plants. Annu Rev Ecol Evol Syst. 2018;49:25–47.Article 

Google Scholar 
Rosche C, Hensen I, Schaar A, Zehra U, Jasieniuk M, Callaway RM, et al. Climate outweighs native vs. nonnative range‐effects for genetics and common garden performance of a cosmopolitan weed. Ecol Monogr. 2019;89:e01386.Article 

Google Scholar 
Weaver SE. The biology of Canadian weeds. 115. Conyza canadensis. Can J Plant Sci. 2001;81:867–75.Article 

Google Scholar 
Gange AC, Ayres RL. On the relation between arbuscular mycorrhizal colonization and plant ’ benefit. Oikos. 1999;87:615–21.Article 

Google Scholar 
Řezáčová V, Konvalinková T, Řezáč M. Decreased mycorrhizal colonization of Conyza canadensis (L.) Cronquist in invaded range does not affect fungal abundance in native plants. Biologia. 2020;75:693–9.Article 

Google Scholar 
Zhang Q, Sun Q, Koide RT, Peng Z, Zhou J, Gu X, et al. Arbuscular mycorrhizal fungal mediation of plant-plant onteractions in a marshland plant community. Sci World J. 2014;2014:1–10.
Google Scholar 
Zhang HY, Goncalves P, Copeland E, Qi SS, Dai ZC, Li GL, et al. Invasion by the weed Conyza canadensis alters soil nutrient supply and shifts microbiota structure. Soil Biol Biochem. 2020;143:107739.Article 
CAS 

Google Scholar 
Shah MA, Callaway RM, Shah T, Houseman GR, Pal RW, Xiao S, et al. Conyza canadensis suppresses plant diversity in its nonnative ranges but not at home: a transcontinental comparison. N Phytol. 2014;202:1286–96.Article 

Google Scholar 
Colautti RI, Lau JA. Contemporary evolution during invasion: evidence for differentiation, natural selection, and local adaptation. Mol Ecol. 2015;24:1999–2017.PubMed 
Article 

Google Scholar 
Rosche C, Hensen I, Lachmuth S. Local pre-adaptation to disturbance and inbreeding-environment interactions affect colonisation abilities of diploid and tetraploid Centaurea stoebe. Plant Biol. 2018;20:75–84.PubMed 
Article 
CAS 

Google Scholar 
Hart SC, Start JM, Davidson EA, Firestone MK. Nitrogen mineralization, immobilization, and nitrification. In: Weaver RW, Angle J., Bottomley P., editors. Methods of soil analysis, part 2 microbiological and biochemical properties. Madison, WI: Soil Science Society of America; 1994. p. 985–1018.Brundrett M, Bougher N, Dell B, Grove T, Malajczuk N. Working with mycorrhizas in forestry and agriculture. ACIAR Monogr. 1996;32:1–374.
Google Scholar 
McGonigle TP, Miller MH, Evans DG, Fairchild GL, Swan JA. A new method which gives an objective measure of colonization of roots by vesicular—arbuscular mycorrhizal fungi. N Phytol. 1990;115:495–501.Article 
CAS 

Google Scholar 
Fick SE, Hijmans RJ. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int J Climatol. 2017;4315:4302–15.Article 

Google Scholar 
Hijmans RJ. raster: Geographic data analysis and modeling. R package version 3.3-13. 2020. https://cran.r-project.org/package=raster.R Core Team. R: A language and environment for statistical computing [https://www.r-project.org/]. Vienna, Austria: R Foundation for Statistical Computing; 2019.Oksanen J, Guillaume BF, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: community ecology package. R package. 2019. https://cran.r-project.org/package=vegan.Dumbrell AJ, Ashton PD, Aziz N, Feng G, Nelson M, Dytham C, et al. Distinct seasonal assemblages of arbuscular mycorrhizal fungi revealed by massively parallel pyrosequencing. N Phytol. 2011;190:794–804.Article 
CAS 

Google Scholar 
Lee J, Lee S, Young JPW. Improved PCR primers for the detection and identification of arbuscular mycorrhizal fungi. FEMS Microbiol Ecol. 2008;65:339–49.PubMed 
Article 
CAS 

Google Scholar 
Bullington LS, Lekberg Y, Larkin BG. Insufficient sampling constrains our characterization of plant microbiomes. Sci Rep. 2021;11:3645.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37:852–7.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.PubMed 
PubMed Central 
CAS 

Google Scholar 
Öpik M, Vanatoa A, Vanatoa E, Moora M, Davison J, Kalwij JM, et al. The online database MaarjAM reveals global and ecosystemic distribution patterns in arbuscular mycorrhizal fungi (Glomeromycota). N Phytol. 2010;188:223–41.Article 
CAS 

Google Scholar 
Chen J. GUniFrac: generalized UniFrac distances. R package version 1.1. 2018. https://cran.r-project.org/package=GUniFrac.Webb CO. Exploring the phylogenetic structure of ecological communities: an example for rain forest trees. Am Nat. 2000;156:145–55.PubMed 
Article 

Google Scholar 
Kembel SW, Cowan PD, Helmus MR, Cornwell WK, Morlon H, Ackerly DD, et al. Picante: R tools for integrating phylogenies and ecology. Bioinformatics. 2010;26:1463–4.PubMed 
Article 
CAS 

Google Scholar 
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Price MN, Dehal PS, Arkin AP. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS ONE. 2010;5:e9490.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Yu G, Smith DK, Zhu H, Guan Y, Lam TT. GGTREE: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017;8:28–36.Article 

Google Scholar 
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2014;67.Borcard D, Gillet F, Legendre P. Numerical ecology with R. New York: Springer; 2011.Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995;57:289–300.
Google Scholar 
Anderson MJ, Walsh DCI. PERMANOVA, ANOSIM, and the Mantel test in the face of heterogeneous dispersions: What null hypothesis are you testing? Ecol Monogr. 2013;83:557–74.Article 

Google Scholar 
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Felker-Quinn E, Schweitzer JA, Bailey JK. Meta-analysis reveals evolution in invasive plant species but little support for Evolution of Increased Competitive Ability (EICA). Ecol Evol. 2013;3:739–51.PubMed 
PubMed Central 
Article 

Google Scholar 
Pal RW, Maron JL, Nagy DU, Waller LP, Tosto A, Liao H, et al. What happens in Europe stays in Europe: apparent evolution by an invader does not help at home. Ecology 2020;101:e03072.PubMed 
Article 

Google Scholar 
Matesanz S, Sultan SE. High-performance genotypes in an introduced plant: insights to future invasiveness. Ecology. 2013;94:2464–74.PubMed 
Article 

Google Scholar 
Hart M, Reader R. Host plant benefit from association with arbuscular mycorrhizal fungi: variation due to differences in size of mycelium. Biol Fertil Soils. 2002;36:357–66.Article 

Google Scholar 
Yang H, Zhang Q, Koide RT, Hoeksema JD, Tang J, Bian X, et al. Taxonomic resolution is a determinant of biodiversity effects in arbuscular mycorrhizal fungal communities. J Ecol. 2017;105:219–28.Article 
CAS 

Google Scholar 
Moora M, Berger S, Davison J, Öpik M, Bommarco R, Bruelheide H, et al. Alien plants associate with widespread generalist arbuscular mycorrhizal fungal taxa: evidence from a continental-scale study using massively parallel 454 sequencing. J Biogeogr. 2011;38:1305–17.Article 

Google Scholar 
Policelli N, Bruns TD, Vilgalys R, Nuñez MA. Suilloid fungi as global drivers of pine invasions. N Phytol. 2019;222:714–25.Article 

Google Scholar 
Jia Y, Heijden MGA, Wagg C, Feng G, Walder F. Symbiotic soil fungi enhance resistance and resilience of an experimental grassland to drought and nitrogen deposition. J Ecol. 2021;109:3171–81.Article 
CAS 

Google Scholar 
Van Der Heijden MGAA, Klironomos JN, Ursic M, Moutoglis P, Streitwolf-Engel R, Boller T, et al. Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity. Nature. 1998;396:69–72.Article 
CAS 

Google Scholar 
Zhang Q, Yang R, Tang J, Yang H, Hu S, Chen X. Positive feedback between mycorrhizal fungi and plants influences plant invasion success and resistance to invasion. PLoS ONE. 2010;5:e12380.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Shah MA, Reshi ZA, Khasa DP. Arbuscular mycorrhizas: drivers or passengers of alien plant invasion. Bot Rev. 2009;75:397–417.Article 

Google Scholar 
Valverde-Barrantes OJ, Horning AL, Smemo KA, Blackwood CB. Phylogenetically structured traits in root systems influence arbuscular mycorrhizal colonization in woody angiosperms. Plant Soil. 2016;404:1–12.Article 
CAS 

Google Scholar 
Wilson GWT, Hartnett DC. Interspecific variation in plant responses to mycorrhizal colonization in tallgrass prairie. Am J Bot. 1998;85:1732–8.PubMed 
Article 
CAS 

Google Scholar 
Seifert EK, Bever JD, Maron JL. Evidence for the evolution of reduced mycorrhizal dependence during plant invasion. Ecology 2009;90:1055–62.PubMed 
Article 

Google Scholar 
Deveautour C, Donn S, Power SA, Bennett AE, Powell JR. Experimentally altered rainfall regimes and host root traits affect grassland arbuscular mycorrhizal fungal communities. Mol Ecol. 2018;27:2152–63.PubMed 
Article 

Google Scholar 
Osborne OG, De-Kayne R, Bidartondo MI, Hutton I, Baker WJ, Turnbull CGN, et al. Arbuscular mycorrhizal fungi promote coexistence and niche divergence of sympatric palm species on a remote oceanic island. N Phytol. 2018;217:1254–66.Article 
CAS 

Google Scholar 
Tian B, Pei Y, Huang W, Ding J, Siemann E. Increasing flavonoid concentrations in root exudates enhance associations between arbuscular mycorrhizal fungi and an invasive plant. ISME J. 2021;15:1919–30.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Pimprikar P, Gutjahr C. Transcriptional regulation of arbuscular mycorrhiza development. Plant Cell Physiol. 2018;59:673–90.PubMed 
Article 
CAS 

Google Scholar 
Wendlandt CE, Helliwell E, Roberts M, Nguyen KT, Friesen ML, Wettberg E, et al. Decreased coevolutionary potential and increased symbiont fecundity during the biological invasion of a legume‐rhizobium mutualism. Evolution. 2021;75:731–47.PubMed 
Article 

Google Scholar 
Callaway RM, Bedmar EJ, Reinhart KO, Silvan CG, Klironomos J. Effects of soil biota from different ranges on Robinia invasion: acquiring mutualists and escaping pathogens. Ecology. 2011;92:1027–35.PubMed 
Article 

Google Scholar 
Shelby N, Duncan RP, Putten WH, McGinn KJ, Weser C, Hulme PE. Plant mutualisms with rhizosphere microbiota in introduced versus native ranges. J Ecol. 2016;104:1259–70.Article 
CAS 

Google Scholar 
Yang Q, Carrillo J, Jin H, Shang L, Hovick SM, Nijjer S, et al. Plant–soil biota interactions of an invasive species in its native and introduced ranges: Implications for invasion success. Soil Biol Biochem. 2013;65:78–85.Article 
CAS 

Google Scholar 
Bronstein JL. The exploitation of mutualisms. Ecol Lett. 2001;4:277–87.Article 

Google Scholar 
Kiers ET, Duhamel M, Beesetty Y, Mensah JA, Franken O, Verbruggen E, et al. Reciprocal rewards stabilize cooperation in the mycorrhizal symbiosis. Science. 2011;333:880–2.PubMed 
Article 
CAS 

Google Scholar 
Koziol L, Bever JD. Mycorrhizal feedbacks generate positive frequency dependence accelerating grassland succession. J Ecol. 2019;107:622–32.Article 

Google Scholar 
Yang H, Yuan Y, Zhang Q, Tang J, Liu Y, Chen X. Changes in soil organic carbon, total nitrogen, and abundance of arbuscular mycorrhizal fungi along a large-scale aridity gradient. Catena. 2011;87:70–7.Article 
CAS 

Google Scholar 
Zhang J, Wang F, Che R, Wang P, Liu H, Ji B, et al. Precipitation shapes communities of arbuscular mycorrhizal fungi in Tibetan alpine steppe. Sci Rep. 2016;6:23488.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Read DJ. Mycorrhizas in ecosystems. Experientia. 1991;47:376–91.Article 

Google Scholar 
Delavaux CS, Smith-Ramesh LM, Kuebbing SE. Beyond nutrients: a meta-analysis of the diverse effects of arbuscular mycorrhizal fungi on plants and soils. Ecology. 2017;98:2111–9.PubMed 
Article 

Google Scholar  More

Charlesworth, D., Barton, N. H. & Charlesworth, B. The sources of adaptive variation. Proc. R. Soc. B Biol. Sci. 284(1855), 20162864 (1855).Article 
CAS 

Google Scholar 
Whitlock, M. C. Fixation of new alleles and the extinction of small populations: Drift load, beneficial alleles, and sexual selection. Evolution 54(6), 1855–1861 (2000).CAS 
PubMed 
Article 

Google Scholar 
Hamilton, W. D. & Zuk, M. Heritable true fitness and bright birds: a role for parasites?. Science 218(4570), 384 (1982).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Hadany, L. & Beker, T. Sexual selection and the evolution of obligatory sex. BMC Evol. Biol. 7(1), 245 (2007).PubMed 
PubMed Central 
Article 

Google Scholar 
Clutton-Brock, T. Reproductive competition and sexual selection. Philos. Trans. R. Soc. B Biol. Sci. 372(1729), 20160310 (2017).Article 

Google Scholar 
Taddei, F. et al. Role of mutator alleles in adaptive evolution. Nature 387(6634), 700–702 (1997).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Agrawal, A. F. & Wang, A. D. Increased transmission of mutations by low-condition females: Evidence for condition-dependent DNA repair. PLoS Biol. 6(2), e30 (2008).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Petrie, M. & Roberts, G. Sexual selection and the evolution of evolvability. Heredity 98(4), 198–205 (2007).CAS 
PubMed 
Article 

Google Scholar 
Dugand, R. J., Kennington, W. J. & Tomkins, J. L. Evolutionary divergence in competitive mating success through female mating bias for good genes. Sci. Adv. 4(5), eaaq0369 (2018).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Siller, S. Sexual selection and the maintenance of sex. Nature 411(6838), 689–692 (2001).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Agrawal, A. F. Sexual selection and the maintenance of sexual reproduction. Nature 411(6838), 692–695 (2001).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Lehtonen, J., Jennions, M. D. & Kokko, H. The many costs of sex. Trends Ecol. Evol. 27(3), 172–178 (2012).PubMed 
Article 

Google Scholar 
Maynard Smith, J. What use is sex?. J. Theor. Biol. 30(2), 319–335 (1971).ADS 
MathSciNet 
Article 

Google Scholar 
Trivers, R. L. Parental investment and sexual selection. In Sexual Selection and the Descent of Man 1871–1971 (ed. Campbell, B.) 136–179 (Aldone, 1972).
Google Scholar 
Petrie, M. & Lipsitch, M. Avian polygyny is most likely in populations with high variability in heritable male fitness. Proc. R. Soc. Lond. Ser. B Biol. Sci. 256(1347), 275–280 (1994).ADS 
Article 

Google Scholar 
Lumley, A. J. et al. Sexual selection protects against extinction. Nature 522(7557), 470–473 (2015).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Andersson, M. Sexual Selection (Princeton University Press, 1994).Book 

Google Scholar 
Petrie, M. Improved growth and survival of offspring of peacocks with more elaborate trains. Nature 371(6498), 598–599 (1994).ADS 
CAS 
Article 

Google Scholar 
Møller, A. P. & Alatalo, R. V. Good-genes effects in sexual selection. Proc. R. Soc. Lond. Ser. B Biol. Sci. 266(1414), 85–91 (1999).Article 

Google Scholar 
David, P. et al. Condition-dependent signalling of genetic variation in stalk-eyed flies. Nature 406(6792), 186–188 (2000).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Hale, M. L. et al. Is the peacock’s train an honest signal of genetic quality at the major histocompatibility complex?. J. Evol. Biol. 22(6), 1284–1294 (2009).CAS 
PubMed 
Article 

Google Scholar 
Prokop, Z. M. et al. Meta-analysis suggests choosy females get sexy sons more than “good genes”. Evolution 66(9), 2665–2673 (2012).PubMed 
Article 

Google Scholar 
Kokko, H. et al. The sexual selection continuum. Proc. R. Soc. Lond. Ser. B Biol. Sci. 269(1498), 1331–1340 (2002).Article 

Google Scholar 
Drake, J. W. et al. Rates of Spontaneous Mutation. Genetics 148(4), 1667 (1998).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Keightley, P. D. Rates and fitness consequences of new mutations in humans. Genetics 190(2), 295–304 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
Haag-Liautard, C. et al. Direct estimation of per nucleotide and genomic deleterious mutation rates in Drosophila. Nature 445(7123), 82–85 (2007).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Metzgar, D. & Wills, C. Evidence for the adaptive evolution of mutation rates. Cell 101(6), 581–584 (2000).CAS 
PubMed 
Article 

Google Scholar 
Janetos, A. C. Strategies of female mate choice: A theoretical analysis. Behav. Ecol. Sociobiol. 7(2), 107–112 (1980).Article 

Google Scholar 
Johnstone, R. A. & Earn, D. J. D. Imperfect female choice and male mating skew on leks of different sizes. Behav. Ecol. Sociobiol. 45(3), 277–281 (1999).Article 

Google Scholar 
Petrie, M., Halliday, T. & Sanders, C. Peahens prefer peacocks with elaborate trains. Anim. Behav. 41(2), 323–331 (1991).Article 

Google Scholar 
Cally, J. G., Stuart-Fox, D. & Holman, L. Meta-analytic evidence that sexual selection improves population fitness. Nat. Commun. 10(1), 2017 (2019).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Kotiaho, J. S. et al. On the resolution of the lek paradox. Trends Ecol. Evol. 23(1), 1–3 (2008).PubMed 
Article 

Google Scholar 
Parker, G. A., Baker, R. R. & Smith, V. G. F. The origin and evolution of gamete dimorphism and the male-female phenomenon. J. Theor. Biol. 36(3), 529–553 (1972).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Parker, G. A. The sexual cascade and the rise of pre-ejaculatory (Darwinian) sexual selection, sex roles, and sexual conflict. Cold Spring Harb. Perspect. Biol. 6(10), a017509–a017509 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
Rowe, L. & Houle, D. The lek paradox and the capture of genetic variance by condition dependent traits. Proc. R. Soc. Lond. Ser. B Biol. Sci. 263(1375), 1415–1421 (1996).ADS 
Article 

Google Scholar 
Petrie, M. Evolution by sexual selection. Front. Ecol. Evol. 9, 950 (2021).Article 

Google Scholar 
Petrie, M. & Kempenaers, B. Extra-pair paternity in birds: Explaining variation between species and populations. Trends Ecol. Evol. 13(2), 52–58 (1998).CAS 
PubMed 
Article 

Google Scholar 
Møller, A. P. & Cuervo, J. J. Minisatellite mutation rates increase with extra-pair paternity among birds. BMC Evol. Biol. 9(1), 100 (2009).PubMed 
PubMed Central 
Article 

Google Scholar 
Anmarkrud, J. A. et al. Factors affecting germline mutations in a hypervariable microsatellite: A comparative analysis of six species of swallows (Aves: Hirundinidae). Mutat. Res. Fundam. Mol. Mech. Mutagen. 708(1), 37–43 (2011).CAS 
Article 

Google Scholar 
Ellegren, H. Characteristics, causes and evolutionary consequences of male-biased mutation. Proc. R. Soc. B Biol. Sci. 274(1606), 1–10 (2007).CAS 
Article 

Google Scholar 
Baur, J. & Berger, D. Experimental evidence for effects of sexual selection on condition-dependent mutation rates. Nat. Ecol. Evol. 4, 737–744 (2020).PubMed 
Article 

Google Scholar 
Vrijenhoek, R. C. & Parker, E. D. Geographical parthenogenesis: General purpose genotypes and frozen niche variation. In Lost Sex (eds Schön, I. et al.) (Springer, Dordrecht, 2009).
Google Scholar 
Reudink, M. W. et al. Evolution of song and color in island birds. Wilson J. Ornithol. 133(1), 1–10 (2021).Article 

Google Scholar 
Iglesias-Carrasco, M. et al. Sexual selection, body mass and molecular evolution interact to predict diversification in birds. Proc. R. Soc. B Biol. Sci. 2019(286), 20190172 (1899).
Google Scholar 
Earl, D. J. & Deem, M. W. Evolvability is a selectable trait. Proc. Natl. Acad. Sci. U. S. A. 101(32), 11531 (2004).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar  More