Philippa Kaur

Geisen S, Wall DH, van der Putten WH. Challenges and opportunities for soil biodiversity in the anthropocene. Curr Biol. 2019;29:R1036–44.Article 
CAS 
PubMed 

Google Scholar 
Fierer N. Embracing the unknown: disentangling the complexities of the soil microbiome. Nat Rev Microbiol. 2017;15:579–90.Article 
CAS 
PubMed 

Google Scholar 
Gossner MM, Lewinsohn TM, Kahl T, Grassein F, Boch S, Prati D, et al. Land-use intensification causes multitrophic homogenization of grassland communities. Nature. 2016;540:266–9.Article 
CAS 
PubMed 

Google Scholar 
Leff JW, Jones SE, Prober SM, Barberán A, Borer ET, Firn JL, et al. Consistent responses of soil microbial communities to elevated nutrient inputs in grasslands across the globe. Proc Natl Acad Sci USA. 2015;112:10967–72.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Alberti M, Correa C, Marzluff JM, Hendry AP, Palkovacs EP, Gotanda KM, et al. Global urban signatures of phenotypic change in animal and plant populations. Proc Natl Acad Sci USA. 2017;114:8951–6.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
El-Sabaawi R. Trophic structure in a rapidly urbanizing planet. Funct Ecol. 2018;32:1718–28.Article 

Google Scholar 
Yu S, Wu Z, Xu G, Li C, Wu Z, Li Z, et al. Inconsistent patterns of soil fauna biodiversity and soil physicochemical characteristic along an urbanization gradient. Front Ecol Evol. 2022;9:824004.Article 

Google Scholar 
Zambrano L, Aronson MFJ, Fernandez T. The consequences of landscape fragmentation on socio-ecological patterns in a rapidly developing urban area: a case study of the National Autonomous University of Mexico. Front. Environ Sci. 2019;7:152.
Google Scholar 
Wilson MC, Chen XY, Corlett RT, Didham RK, Ding P, Holt RD, et al. Habitat fragmentation and biodiversity conservation: key findings and future challenges. Landsc Ecol. 2016;31:219–27.Article 

Google Scholar 
Guilland C, Maron PA, Damas O, Ranjard L. Biodiversity of urban soils for sustainable cities. Environ Chem Lett. 2018;16:1267–82.Article 
CAS 

Google Scholar 
Dou Y, Kuang W. A comparative analysis of urban impervious surface and green space and their dynamics among 318 different size cities in China in the past 25 years. Sci. Total Environ. 2020;706:135828.Article 
CAS 
PubMed 

Google Scholar 
Francini G, Hui N, Jumpponen A, Kotze D, Romantschuk M, Allen J, et al. Soil biota in boreal urban greenspace: responses to plant type and age. Soil Biol Biochem. 2018;118:145–55.Article 
CAS 

Google Scholar 
Corline NJ, Peek RA, Montgomery J, Katz JVE, Jeffres CA. Understanding community assembly rules in managed floodplain food webs. Ecosphere. 2021;12:e03330.Article 

Google Scholar 
Tripathi BM, Stegen JC, Kim M, Dong K, Adams JM, Lee YK. Soil pH mediates the balance between stochastic and deterministic assembly of bacteria. ISME J. 2018;12:1072–83.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Liu W, Graham EB, Dong Y, Zhong L, Zhang J, Qiu C, et al. Balanced stochastic versus deterministic assembly processes benefit diverse yet uneven ecosystem functions in representative agroecosystems. Environ Microbiol. 2021;23:391–404.Article 
CAS 
PubMed 

Google Scholar 
Thakur MP, Phillips HR, Brose U, De Vries FT, Lavelle P, Loreau M, et al. Towards an integrative understanding of soil biodiversity. Biol Rev. 2020;95:350–64.Article 
PubMed 

Google Scholar 
Bahram M, Kohout P, Anslan S, Harend H, Abarenkov K, Tedersoo L. Stochastic distribution of small soil eukaryotes resulting from high dispersal and drift in a local environment. ISME J. 2016;10:885–96.Article 
PubMed 

Google Scholar 
Luan L, Jiang Y, Cheng M, Dini-Andreote F, Sui Y, Xu Q, et al. Organism body size structures the soil microbial and nematode community assembly at a continental and global scale. Nat Commun. 2020;11:1–11.Article 

Google Scholar 
Isabwe A, Yang JR, Wang Y, Wilkinson DM, Graham EB, Chen H, et al. Riverine bacterioplankton and phytoplankton assembly along an environmental gradient induced by urbanization. Limnol Oceanogr. 2022;67:1943–58.Article 
CAS 

Google Scholar 
Nemergut DR, Schmidt SK, Fukami T, O’Neill SP, Bilinski TM, Stanish LF, et al. Patterns and processes of microbial community assembly. Microbiol Mol Biol Rev. 2013;77:342–56.Article 
PubMed 
PubMed Central 

Google Scholar 
Zinger L, Taberlet P, Schimann H, Bonin A, Boyer F, De Barba M, et al. Body size determines soil community assembly in a tropical forest. Mol Ecol. 2019;28:528–43.Article 
CAS 
PubMed 

Google Scholar 
Jiao S, Yang Y, Xu Y, Zhang J, Lu Y. Balance between community assembly processes mediates species coexistence in agricultural soil microbiomes across eastern China. ISME J. 2020;14:202–16.Article 
PubMed 

Google Scholar 
Jiao S, Chen W, Wei G. Biogeography and ecological diversity patterns of rare and abundant bacteria in oil‐contaminated soils. Mol Ecol. 2017;26:5305–17.Article 
CAS 
PubMed 

Google Scholar 
Wu W, Lu H-P, Sastri A, Yeh Y-C, Gong G-C, Chou W-C, et al. Contrasting the relative importance of species sorting and dispersal limitation in shaping marine bacterial versus protist communities. ISME J. 2018;12:485–94.Article 
PubMed 

Google Scholar 
Farjalla VF, Srivastava DS, Marino NA, Azevedo FD, Dib V, Lopes PM, et al. Ecological determinism increases with organism size. Ecology. 2012;93:1752–9.Article 
PubMed 

Google Scholar 
Carscadden KA, Emery NC, Arnillas CA, Cadotte MW, Afkhami ME, Gravel D, et al. Niche breadth: causes and consequences for ecology, evolution, and conservation. Q Rev Biol. 2020;95:179–214.Article 

Google Scholar 
Beissinger SR. Ecological mechanisms of extinction. Proc Natl Acad Sci USA. 2000;97:11688–9.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Poiani KA, Richter BD, Anderson MG, Richter HE. Biodiversity conservation at multiple scales: functional sites, landscapes, and networks. Bioscience. 2000;50:133–46.Article 

Google Scholar 
Yang J, Zhang X, Jin X, Seymour M, Richter C, Logares R, et al. Recent advances in environmental DNA-based biodiversity assessment and conservation. Divers Distrib. 2021;27:1876–9.Article 

Google Scholar 
Breed MF, Harrison PA, Blyth C, Byrne M, Gaget V, Gellie NJC, et al. The potential of genomics for restoring ecosystems and biodiversity. Nat Rev Genet. 2019;20:615–28.Article 
CAS 
PubMed 

Google Scholar 
Department of Economic and Social Affairs (DESA). World Urbanization Prospects. The 2018 Revision. United Nations. 2019. https://population.un.org/wup/publications/Files/WUP2018-Report.pdf. Accessed 13 Mar 2022.Qiao Z, Wang B, Yao H, Li Z, Scheu S, Zhu Y-G, et al. Urbanization and greenspace type as determinants of species and functional composition of collembola communities. Geoderma. 2022;428:116175.Article 

Google Scholar 
Shrestha S, Cui S, Xu L, Wang L, Manandhar B, Ding S. Impact of land use change due to urbanisation on surface runoff using GIS-based SCS–CN Method: a case study of Xiamen City, China. Land. 2021;10:839.Article 

Google Scholar 
R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing. 2022. Vienna, Austria. https://www.R-project.org/.Wickham. H ggplot2: elegant graphics for data analysis. Springer-Verlag New York, 2016.Kassambara A. ggpubr: ‘ggplot2’ based publication ready plots. 2020. https://CRAN.R-project.org/package=ggpubr.Morlon H, Chuyong G, Condit R, Hubbell S, Kenfack D, Thomas D, et al. A general framework for the distance–decay of similarity in ecological communities. Ecol Lett. 2008;11:904–17.Article 
PubMed 
PubMed Central 

Google Scholar 
Goslee S, D. Urban, Goslee, MS. ecodist: dissimilarity-based functions for rcological analysis. 2020. https://cran.r-project.org/web/packages/ecodist/index.html.Ofiţeru ID, Lunn M, Curtis TP, Wells GF, Criddle CS, Francis CA, et al. Combined niche and neutral effects in a microbial wastewater treatment community. Proc Natl Acad Sci USA. 2010;107:15345–50.Article 
PubMed 
PubMed Central 

Google Scholar 
Burns AR, Stephens WZ, Stagaman K, Wong S, Rawls JF, Guillemin K, et al. Contribution of neutral processes to the assembly of gut microbial communities in the zebrafish over host development. ISME J. 2016;10:655–64.Article 
CAS 
PubMed 

Google Scholar 
Chen W, Ren K, Isabwe A, Chen H, Liu M, Yang J. Stochastic processes shape microeukaryotic community assembly in a subtropical river across wet and dry seasons. Microbiome. 2019;7:138.Article 
PubMed 
PubMed Central 

Google Scholar 
Chase JM, Kraft NJ, Smith KG, Vellend M, Inouye BD. Using null models to disentangle variation in community dissimilarity from variation in α‐diversity. Ecosphere. 2011;2:1–11.Article 

Google Scholar 
Pandit SN, Kolasa J, Cottenie K. Contrasts between habitat generalists and specialists: an empirical extension to the basic metacommunity framework. Ecology. 2009;90:2253–62.Article 
PubMed 

Google Scholar 
Salazar G. EcolUtils: utilities for community ecology analysis. 2019. https://github.com/GuillemSalazar/EcolUtils.Kraft NJB, Adler PB, Godoy O, James EC, Fuller S, Levine JM. Community assembly, coexistence and the environmental filtering metaphor. Funct Ecol. 2015;29:592–9.Article 

Google Scholar 
Cadotte MW, Tucker CM. Should environmental filtering be abandoned? Trends Ecol Evol. 2017;32:429–37.Article 
PubMed 

Google Scholar 
Leibold MA, McPeek MA. Coexistence of the niche and neutral perspectives in community ecology. Ecology. 2006;87:1399–410.Article 
PubMed 

Google Scholar 
Evans S, Martiny JB, Allison SD. Effects of dispersal and selection on stochastic assembly in microbial communities. ISME J. 2017;11:176–85.Article 
PubMed 

Google Scholar 
Jiang Y, Liu M, Zhang J, Chen Y, Chen X, Chen L, et al. Nematode grazing promotes bacterial community dynamics in soil at the aggregate level. ISME J. 2017;11:2705–17.Article 
PubMed 
PubMed Central 

Google Scholar 
Douhan GW, Vincenot L, Gryta H, Selosse M-A. Population genetics of ectomycorrhizal fungi: from current knowledge to emerging directions. Fungal Biol. 2011;115:569–97.Article 
PubMed 

Google Scholar 
Granot I, Belmaker J. Niche breadth and species richness: correlation strength, scale and mechanisms. Glob Ecol Biogeogr. 2020;29:159–70.Article 

Google Scholar 
Sexton JP, Montiel J, Shay JE, Stephens MR, Slatyer RA. Evolution of ecological niche breadth. Annu Rev Ecol Evol Syst Annu Rev Ecol Evol S. 2017;48:183–206.Article 

Google Scholar 
Fraaije RGA, ter Braak CJF, Verduyn B, Verhoeven JTA, Soons MB. Dispersal versus environmental filtering in a dynamic system: drivers of vegetation patterns and diversity along stream riparian gradients. J Ecol. 2015;103:1634–46.Article 

Google Scholar 
Soininen J, McDonald R, Hillebrand H. The distance decay of similarity in ecological communities. Ecography. 2007;30:3–12.Article 

Google Scholar 
Zhang K, Delgado-Baquerizo M, Zhu Y-G, Chu H. Space is more important than season when shaping soil microbial communities at a large spatial scale. mSystems. 2020;5:e00783–19.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ma B, Dai Z, Wang H, Dsouza M, Liu X, He Y, et al. Distinct biogeographic patterns for archaea, bacteria, and fungi along the vegetation gradient at the continental scale in Eastern China. mSystems. 2017;2:e00174–16.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Wang J, Zhang T, Li L, Li J, Feng Y, Lu Q. The patterns and drivers of bacterial and fungal β-diversity in a typical dryland ecosystem of northwest China. Front Microbiol. 2017;8:2126.Article 
PubMed 
PubMed Central 

Google Scholar 
Kang L, Chen L, Zhang D, Peng Y, Song Y, Kou D, et al. Stochastic processes regulate belowground community assembly in alpine grasslands on the Tibetan Plateau. Environ Microbiol. 2021;24:179–94.Article 
PubMed 

Google Scholar 
Chen Q-L, Hu H-W, Yan Z-Z, Li C-Y, Nguyen B-AT, Sun A-Q, et al. Deterministic selection dominates microbial community assembly in termite mounds. Soil Biol Biochem. 2021;152:108073.Article 
CAS 

Google Scholar 
Huang S, Tucker MA, Hertel AG, Eyres A, Albrecht J. Scale-dependent effects of niche specialisation: the disconnect between individual and species ranges. Ecol Lett. 2021;24:1408–19.Article 
PubMed 

Google Scholar 
Rapacciuolo G, Blois JL. Understanding ecological change across large spatial, temporal and taxonomic scales: integrating data and methods in light of theory. Ecography. 2019;42:1247–66.
Google Scholar 
van der Gast CJ. Microbial biogeography: the end of the ubiquitous dispersal hypothesis? Environ Microbiol. 2015;17:544–6.Article 
PubMed 

Google Scholar 
Levy-Booth DJ, Giesbrecht IJW, Kellogg CTE, Heger TJ, D’Amore DV, Keeling PJ, et al. Seasonal and ecohydrological regulation of active microbial populations involved in DOC, CO2, and CH4 fluxes in temperate rainforest soil. ISME J. 2019;13:950–63.Article 
CAS 
PubMed 

Google Scholar 
De Gannes V, Bekele I, Dipchansingh D, Wuddivira MN, De Cairies S, Boman M, et al. Microbial community structure and function of soil following ecosystem conversion from native forests to teak plantation forests. Front Microbiol. 2016;7:1976.Article 
PubMed 
PubMed Central 

Google Scholar 
Männistö M, Vuosku J, Stark S, Saravesi K, Suokas M, Markkola A, et al. Bacterial and fungal communities in boreal forest soil are insensitive to changes in snow cover conditions. FEMS Microbiol. 2018;94:fiy123.
Google Scholar 
Sakarika M, Spanoghe J, Sui Y, Wambacq E, Grunert O, Haesaert G, et al. Purple non‐sulphur bacteria and plant production: benefits for fertilization, stress resistance and the environment. Microb Biotechnol. 2020;13:1336–65.Article 
CAS 
PubMed 

Google Scholar 
Kernaghan G, Patriquin G. Diversity and host preference of fungi co-inhabiting Cenococcum mycorrhizae. Fungal Ecol. 2015;17:84–95.Article 

Google Scholar 
Lumibao CY, Kimbrough ER, Day RH, Conner WH, Krauss KW, Van Bael SA. Divergent biotic and abiotic filtering of root endosphere and rhizosphere soil fungal communities along ecological gradients. FEMS Microbiol. 2020;96:fiaa124.Article 
CAS 

Google Scholar 
Rueckert S, Betts EL, Tsaousis AD. The symbiotic spectrum: where do the gregarines fit? Trends Parasitol. 2019;35:687–94.Article 
PubMed 

Google Scholar 
Butaeva F, Paskerova G, Entzeroth R. Ditrypanocystis sp.(Apicomplexa, Gregarinia, Selenidiidae): the mode of survival in the gut of Enchytraeus albidus (Annelida, Oligochaeta, Enchytraeidae) is close to that of the coccidian genus Cryptosporidium. Tsitologiia. 2006;48:695–704.CAS 
PubMed 

Google Scholar 
Pavao-Zuckerman MA, Coleman DC. Urbanization alters the functional composition, but not taxonomic diversity, of the soil nematode community. Appl Soil Ecol. 2007;35:329–39.Article 

Google Scholar 
Gaspar C, Borges PA, Gaston KJ. Diversity and distribution of arthropods in native forests of the Azores archipelago. Arquipelago: Life Mar Sci. 2008;25:1–30.
Google Scholar 
Suter RB, Doyle G, Shane CM. Oviposition site selection by Frontinella pyramitela (Araneae, Linyphiidae). J Arachnol. 1987;15:349–54.Tian T, Ren Q, Fan J, Haseeb M, Zhang R. Too dry or too wet soils have a negative impact on larval pupation of fall armyworm. J Appl Entomol. 2022;146:196–202.Article 

Google Scholar 
Marczylo EL, Macchiarulo S, Gant TW. Metabarcoding of soil fungi from different urban greenspaces around Bournemouth in the UK. EcoHealth. 2021;18:315–30.Article 
PubMed 
PubMed Central 

Google Scholar 
Corline NJ, Peek RA, Montgomery J, Katz JVE, Jeffres CA. Understanding community assembly rules in managed floodplain food webs. Ecosphere. 2021;12:e03330.Article 

Google Scholar 
Schlägel UE, Grimm V, Blaum N, Colangeli P, Dammhahn M, Eccard JA, et al. Movement-mediated community assembly and coexistence. Biol Rev Camb Philos Soc. 2020;95:1073–96.Article 
PubMed 

Google Scholar 
Stubner S. Enumeration of 16S rDNA of desulfotomaculum lineage 1 in rice field soil by real-time PCR with SybrGreen™ detection. J Microbiol Methods. 2002;50:155–64.Article 
CAS 
PubMed 

Google Scholar 
DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Toju H, Tanabe AS, Yamamoto S, Sato H. High-coverage ITS primers for the DNA-based identification of ascomycetes and basidiomycetes in environmental samples. PloS One. 2012;7:e40863.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Abarenkov K, Henrik Nilsson R, Larsson KH, Alexander IJ, Eberhardt U, Erland S, et al. The UNITE database for molecular identification of fungi–recent updates and future perspectives. New Phytol. 2010;186:281–5.Article 
PubMed 

Google Scholar 
Stoeck T, Bass D, Nebel M, Christen R, Jones MD, Breiner H-W, et al. Multiple marker parallel tag environmental DNA sequencing reveals a highly complex eukaryotic community in marine anoxic water. Mol Ecol. 2010;19:21–31.Article 
CAS 
PubMed 

Google Scholar 
Guillou L, Bachar D, Audic S, Bass D, Berney C, Bittner L, et al. The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote small sub-unit rRNA sequences with curated taxonomy. Nucleic Acids Res. 2012;41:D597–604.Article 
PubMed 
PubMed Central 

Google Scholar 
Porazinska DL, Giblin‐Davis RM, Faller L, Farmerie W, Kanzaki N, Morris K, et al. Evaluating high‐throughput sequencing as a method for metagenomic analysis of nematode diversity. Mol Ecol Res. 2009;9:1439–50.Article 
CAS 

Google Scholar 
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2012;41:D590–96.Article 
PubMed 
PubMed Central 

Google Scholar 
Leray M, Yang JY, Meyer CP, Mills SC, Agudelo N, Ranwez V, et al. A new versatile primer set targeting a short fragment of the mitochondrial COI region for metabarcoding metazoan diversity: application for characterizing coral reef fish gut contents. Front Zool. 2013;10:1–14.Article 

Google Scholar 
Porter TM, Hajibabaei M. Over 2.5 million COI sequences in GenBank and growing. PloS One. 2018;13:e0200177.Article 
PubMed 
PubMed Central 

Google Scholar  More

This study quantitatively assessed waterbody loss due to urban expansion of large Chinese cities. We first extracted multi-temporal urban boundaries to determine the expansion of cities of over one million in population from 1990 to 2018. The monthly surface-water dataset was then used to identify surface waterbodies in the study period. Depending on the ratio of surface waterbody area to urban area, cities were further divided into three categories (i.e. water-abundant, water-medium, water-deficient). Finally, we quantified the rate of waterbody loss and evaluated the spatial and temporal variation of waterbody loss as a function of urban expansion and according to city type.GUB datasetThe Global Urban Boundary (GUB) dataset (http://data.ess.tsinghua.edu.cn) was used to determine urban expansion. GUB provides data on built-up areas over 30 years, with a spatial resolution of 30 m. In the GUB dataset, nonurban areas (such as green space and water space) surrounded by artificial impervious areas are filled within the urban boundary and removed by the algorithm, which is consistent with global mapping methods. The continuous urban boundary was demarcated by morphological image processing methods, which have an overall accuracy of over 90%. In this dataset, extensive water and forests are excluded, and the impervious surface within the urban boundaries accounts for about 60% of the total surface area47. Compared with urban boundaries obtained from night-time light, GUB better separates urban areas from surrounding nonurban areas.Monthly waterbody datasetWe selected the JRC Monthly Water History V1.3 dataset(https://global-surface-water.appspot.com/), which is available from the Google Earth Engine, as the basis for representing surface waterbodies48. This data collection, which was produced by using images from the Landsat series, contains 442 images of global monthly waterbody area from March 1984 to December 2020. In this dataset, the validation confirmed that fewer than 1% of waterbodies were incorrectly detected, and fewer than 5% of waterbodies were missed altogether. We chose this dataset due to the long-term spatial distribution of waterbodies and due to mountain shadows and urban-constructions masking, which reflects the real changes in waterbodies.Theoretical backgroundIt is well known that cities have high concentrations of population and resources and expand spatially during development. There are many different perspectives on the size of cities, and studies have mostly used urban density and population to characterize them. However, because it is challenging to standardize data sources and quality, there is no unified quantitative standard49. Urban construction has concentrated human activity and brought about changes in land types. Cities are also identified as physical spaces, which can be defined as the built environment50,51. The built environment, which includes structures like buildings, roads, and other artificial constructions, is sometimes referred to as a non-natural environment52.Rural is the antithesis of urban. As large cities have spread outward in developing nations like Asia, a transitional fringe has been created by the gradual blurring of the line separating urban and rural areas53. According to McGee, good locations, easy access, and sizable agricultural land all contribute to the development potential of large cities. Thus, between urban and rural areas, there are transitional areas of active spatial morphological change known as desakota33,54. The peri-urban areas, like desakota, are gradually developed and incorporated into original built-up urban areas in urbanization. The original landscape, which included agricultural land, vegetation, and waterbodies, gradually changed into an urban land use type, i.e. impervious surface, and thus the city continues to expand outwards. Waterbody, an essential ecological element, has been heavily developed or filled in during urbanization, which may present dangerous ecological risks. In this paper, we identified the urban boundaries based on physical space to explore the encroachment activities on waterbodies during the urbanization of large cities. We determined whether existing waterbodies were transformed into urban waterbodies or encroached upon and whether waterbodies were increased in the expansion of urban boundaries, thus proposing strategies for protecting waterbodies in the future.Extracting the extent of large Chinese cities from GUB datasetTo characterize urban expansion, GUB data are selected as the original data for urban boundary selection. The Chinese administrative scale of municipalities is not exclusively urban, but also includes rural areas. In our study, cities were defined as municipal districts excluding the vast countryside within the administrative boundaries of prefecture-level cities. We identified urban areas based on the physical boundaries from the perspective of remote sensing, which can precisely track urban expansion51.In this work, we selected 159 cities with a population of over one million in 2018 based on the average annual population of urban districts from the 2019 China City Statistical Yearbook (Fig. S1). Taiwan, Hong Kong, and Macau are omitted. According to statistics, China had 160 cities with populations exceeding one million in 2018. However, due to the lack of data for the built-up area in 1990, Guang’an was not included in the study. We thus obtained 159 cities from the GUB dataset. Due to numerous fragmented patches within the administrative boundary, the population identified the main urban areas, and max patch areas were comprehensively based on the urban boundaries. Through manual detection and adjustment of the map, we determined that the location of the extracted urban area was consistent with that of the municipal government, and the boundary was extracted for each period. We took the growth area as the expansion area, with the original area being the city at the onset of each period (Fig. S3).We used the average annual urban growth (AUG) rate to characterize the rate of urban expansion, as is widely done to evaluate urban expansion55,56. It is calculated as$${text{AUG}} = left[ {frac{{Land_{t1} }}{{Land_{t0} }}^{{frac{1}{t1 – t0}}} – 1} right] times 100% ,$$
where (Land_{t0}) and (Land_{t1}) represent the urban land area at time t0 and t1, where t0 and t1 are the start and end of the given study period.Identification of urban waterbodiesUrban waterbodies contain all the components of urban flow networks above the ground and include natural waterbodies such as lakes, rivers, streams, and wetlands and artificial waterbodies such as parks and ponds48. We identified all waterbodies existing within the urban boundary as urban waterbody. Considering urban expansion, urban waterbodies vary as urban boundary shift at different stages. Our study explored how the original waterbodies changed under urban expansion, including whether they were kept as urban waterbodies or encroached upon. Considering the dryness or wetness of each year, we used the data for 3 years (36 months) around each period (1990, 1995, 2000, 2005, 2010, 2015, and 2018) to describe the waterbody. Not all waterbodies could be detected for each month of the year; for example, freezing may prevent waterbodies from being detected. To cover seasonal and permanent waterbodies, we used the waterbody frequency index (WFI), which is calculated as the fraction of waterbody months within the 3 years to identify stable waterbodies pixel by pixel57. The spatial distribution of each waterbody was then mapped comprehensively for each period. By comparing the extracted waterbody with the long-time-series high-resolution remote-sensing images from Google Earth, we found that the extracted waterbodies fit the actual waterbody distribution quite well (Fig. S2):$$WFIleft( i right) = frac{WMleft( i right)}{{DMleft( i right)}}$$
where WFI(i) is the water occurrence for pixel i in the images before and after the given year, and i is the pixel number for the study area. WM(i) is the number of months during which the waterbody is detected in i pixel over the 3 years. DM(i) is the number of months during which the data are available in pixel i. If the waterbody frequency index of a pixel is greater than 25%, this pixel is considered as a waterbody; otherwise, it is not.City classification based on surface waterbodyCities with over one million in population may not be short of waterbodies, but significant differences remain in surface waterbody abundance. Due to large differences in city size, it is inappropriate to use waterbody area as a criterion. Considering the influence of urban expansion, we ranked 159 cities according to the indicator of waterbody fraction (WF), namely the fraction of the original surface water within the urban boundary in 2018. Waterbodies not impacted by urbanization were taken as the original surface waterbody, which used the average surface waterbody from 1985 to 1991 as baseline. We used the natural break method to divide cities into abundant, moderate, and deficient levels (referred to as Type I, Type II, and Type III, respectively) and evaluate the abundance of waterbodies in cities. Based on the waterbody fraction (WF) value, which is calculated as follows:$${text{WF}} = frac{{Water_{origin} }}{{Land_{2018} }}$$
where WF is used to judge the urban waterbody abundance in cities. (Water_{1990}) is the origin surface waterbody area (used the year in 1985–1991) in the urban boundary of 2018, (Land_{2018}) the urban land area in the urban boundary of 2018.Temporal characteristic of waterbody loss and gainTo understand the spatial–temporal features of surface waterbodies, we used five normalized indicators to compare waterbody variations between cities during urban expansion from the overall perspective and from the city perspective.The variation in original natural waterbodies reflects the intensity of the natural resource development in urban expansion. We summarized the reduction and preservation of original waterbodies in urban expansion areas with a population of over one million to represent the encroachment of urban expansion on waterbodies:$$WL = frac{{sum NWL_{t0_t1} left( i right)}}{{sum W_{t0} left( i right)}} times 100%$$$$WP = frac{{sum (W_{t0} left( i right) – NWL_{t0_t1} left( i right))}}{{sum W_{t0} left( i right)}} times 100%$$
where i labels the city within the 159 cities, WL and WP are the fractions of waterbody loss and preservation in urban expansion areas of all cities, (NWL_{t0_t1}) is the net waterbody loss during period t0–t1 (, and W_{t0}) is the natural waterbody in the urban expansion area at time t0.To estimate the net waterbody loss caused by urban expansion at various stages, we used the standardized indicator, annual average net waterbody loss rate (ANWL), to compare waterbody loss speeds over time. This indicator is independent of the difference in waterbody abundance and can be compared over time. Waterbody loss is one part of the impact of urbanization; the other is waterbody gain. We used the same method to evaluate the annual average net waterbody gain rate (ANWG). The formulas are$$A{text{NWL}} = frac{{NWL_{t0_t1} }}{{W_{t0} left( {t1 – t0} right)}} times 100%$$$$ANWG = frac{{NWG_{t0 – t1} }}{{W_{t0} left( {t1 – t0} right)}} times 100%$$
where NWL and NWG are the net waterbody loss and gain, respectively, and the other abbreviations are the same as above.Considering the direct impact of urban expansion, we used a normalized indicator, the average net waterbody loss velocity of urban expansion ((AWLV)), which refers to the amount of waterbody encroachment per unit urban expansion area. It quantifies the time-heterogeneity of waterbody loss due to urban expansion and is calculated as follows:$$AWLV = frac{{NWL_{t0_t1} }}{{Land_{t1} – Land_{t0} }}$$We calculated these indicators for the six expansion periods (1990–1995, 1995–2000, 2000–2005, 2005–2010, 2010–2015, and 2015–2018) (Fig. 3). In the study, if the waterbody pixel count is zero at the onset of the period, the indicator for the period is abnormal and thus excluded. More

Retrieval of E. coli plasmid sequencesAll E. coli sequences were downloaded from the NCBI FTP server in May 2020. To establish an initial collection of plasmids, only complete genomes with an associated plasmid were retained. All genomes were verified for belonging to the species E. coli using kmerfinder (https://cge.cbs.dtu.dk/services/KmerFinder/). Sequence type (ST) was determined via multi-locus sequence typing (MLST) based on the 7-gene Achtman scheme using pubMLST (https:/github.com/tseemann/mlst). Only genomes with exact matches were assigned for each ST and used for subsequent analysis. To ensure our sequences were sufficiently representative of E. coli pathogens expected in nature, a systematic literature search (see description below and Fig. S1) was conducted to establish an expected distribution of STs (Table S1). This information was used to update our initial collection to match the top 4 most prevalent STs (131, 11, 73, and 95). Specifically, to identify supplementary plasmid sequences, genome accession IDs were chosen from EnteroBase based on the following criteria: the strain was matched to the correct ST and had a high-quality genome sequence (based on N50  > 20,000 and the number of contigs  0.1, 2-tailed student t test). For the second method, all kanR plasmids were used, and instead changed the hosts such that DH5αPro cells were in competition with DH5αPro containing a spontaneous rifampicin-resistant mutant (rifR). Any rifR strain was quantified on rifampicin-containing plates, and the second strain was quantified by rifampicin CFU minus CFU obtained on blank plates. We established that rifR exhibited no fitness defects by (1) growth rates between the wild-type (WT) strain (W) and rifR (M) (Fig. S5D), and (2) directly competing the two control strains (Fig. S5E). In both cases, results were statistically indistinguishable (p  > 0.1, two-tailed student t test). KanR/cmR and WT/rifR experiments were each conducted in LB or M9CAG, respectively. In all cases, experiments were repeated with at least three independent biological replicates.Time-kill measurements in the presence of carbenicillinAll strains were grown as previously described. Time-kill experiments entailed hourly measurements of CFU in presence of carbenicillin at either 3.75 μg/mL (3x IC50) or 5 μg/mL (4x IC50) over a span of 2 or 3 h, including time 0. Specifically, overnight cultures were first diluted 1:100 into LB media containing 1 mM IPTG and 50 μg/mL kanamycin and sub-cultured for two hours in a 37 °C incubator with shaking at 250 rpm. Following this, cell density was adjusted as necessary to achieve a starting OD600 of ~0.15 in all cases. Adjusted subcultures were then aliquoted into a 96-well plate and the appropriate carbenicillin treatments were added directly to the well. Plates were sealed with a paper film and placed in a 37 °C incubator with shaking at 250 rpm. Initial collection for time=0 was acquired before carbenicillin treatment. Thereafter, 10 μL of culture was removed from the well every hour, 10-fold serial dilutions were performed and 10 μL was plated on blank LB agar with three technical replicates at each time point. Colonies were counted after plates were grown for 16 h in a 37 °C incubator to determine CFU. This procedure utilized 14 strains of DH5αPro transformed with kanR plasmids of interest – ctrl, katG, lpxM, yfbR, aroH, pld, fdtC, agp, eptC, arcA, argF, mmuM, ahr, and fabG. CFUs were averaged for all technical replicates, and experiments were conducted with at least three independent biological replicates.Oxygen consumption rateOxygen consumption rates (OCR) were obtained with the Resipher device from Lucid Scientific. The selected strains were grown overnight as previously described. Overnight cultures were resuspended in M9CAG media with 1 mM IPTG and 50 μg/mL kanamycin, and placed in 25 °C for one hour to initiate gene expression. Following this, cells were diluted 10x into M9CAG media containing kanamycin and IPTG, and 100 μL was aliquoted per well into a 96-well microtiter plate according to the manufacturer’s instructions. Plates were placed at 30 °C to minimize growth, and oxygen concentration (μM) was measured immediately thereafter. 24 wells were measured consisting of 6 technical replicates for each strain. Given the clear well-well variability (Fig. S8B, C), data shown are for one biological replicate. However, qualitative trends were consistently reproduced in multiple independent experiments.StatisticsIn all cases where t tests and ANOVA’s were used, data was first verified to be normally distributed using Kolmogorov test for normality. Otherwise, Mann-Whitney U-tests were conducted. For panels with multiple tests, Bonferroni correction was used to adjust the p values. To determine whether any metabolic category was significantly dependent on incompatibility groups, we implemented logistic regressions in MATLAB with the function fitglm. Random forest classification was used to establish the relative importance of prevalent metabolic genes and gene categories predicting the presence of antibiotic resistance genes. Chi-square tests were conducted to determine significant co-occurrence of individual antibiotic resistant and metabolism genes. Dissociative relationships were distinguished by the odds ratios from the chi-square tests. To investigate whether the strong associations and disassociations were driven by evolutionary constraints, or simply artifacts of a common ancestor, we re-ran our statistical analysis using Coinfinder [29] to take in our gene presence-absence data, along with the genome phylogeny, and compute the Bonferroni-corrected statistical likelihood of coincidence (either associations or dissociations), thereby accounting for evolutionary relatedness. More

Sequencing and assembly of the H. schachtii genomeWe measured (Supplemental Fig. 1), sequenced (BioProject PRJNA722882), and assembled the genome of H. schachtii (population Bonn) using a combination of flow cytometry, Pacific Biosciences sequencing, and Illumina sequencing. H. schachtii has the largest genome (160–170 Mb) of any cyst nematode measured/sequenced to date (Supplementary Table 1). It was sequenced to 192-fold coverage using Pacific Biosciences sequencing (fragment n50 of 16 kb), and 144-fold coverage using Illumina sequencing (150 bp Paired-end reads). The final, polished, contamination-free (Supplemental Fig. 2), assembly (v1.2) included ~179 Mbp contained within 395 scaffolds: 90% of the sequence is contained on scaffolds longer than 281,463 bp (n = 154). The assembly is a largely complete haploid representation of the diploid genome, as evidenced by core eukaryotic genes being largely present, complete and single copy (CEGMA 93.15% complete with an average of 1.12 copies each, and BUSCO (Eukaryota odb9) 79% complete with 8.2% duplicated—Supplementary Table 2). Over three million variants were phased into haplotypes (2029 blocks, N50 239.5 kb, covering 94.7% of the reference) which can be used to predict true protein variants (Supplementary data 1), and 601 larger structural variants were identified (Supplementary data 2).The trans-kingdom, lifestage-specific, transcriptomes of H. schachtii and A. thaliana provide a holistic view of parasitismWe devised a sampling procedure to cover all major life stages/transitions of the parasitic life cycle to generate a simultaneous, chronological, and comprehensive picture of nematode gene expression, and infection-site-specific plant gene expression patterns. We sampled cysts and pre-infective second-stage juveniles (J2s), as well as infected segments of A. thaliana root and uninfected adjacent control segments of root at 10 hours post infection (hpi – migratory J2s, pre-establishment of the feeding site), 48 hpi (post establishment of the feeding site), 12 days post infection females (dpi – virgin), 12 dpi males (differentiated, pre-emergence, most if not all stopped feeding), and 24 dpi females (post mating), each in biological triplicate (Fig. 1A). We generated approximately nine billion pairs of 150 bp strand-specific RNAseq reads (Supplementary data 3) covering each stage in biological triplicate (for the parasite and the host): in the early stages of infection we generated over 400 million reads per replicate, to provide sufficient coverage of each kingdom.Fig. 1: Trans-kingdom, lifestage-specific, transcriptome of H. schachtii and A. thaliana.A Schematic representation of the life cycle of H. schachtii infecting A. thaliana, highlighting the 7 stages sampled in this study. For each stage, the average number of trimmed RNAseq read pairs per replicate is shown, with the proportion of reads mapping to either parasite or host in parentheses. B Principle components 1 and 2 for H. schachtii and A. thaliana expression data are plotted. Arrows indicate progression through the life cycle/real-time. Hours post infection (hpi), days post infection (dpi).Full size imageStrand-specific RNAseq reads originating from host and parasite were deconvoluted by mapping to their respective genome assemblies (H. schachtii v.1.2 and TAIR10). For the parasite, ~500 million Illumina RNAseq read pairs uniquely mapping to the H. schachtii genome were used to generate a set of 26,739 gene annotations (32,624 transcripts – detailed further in the next section), ~77% of which have good evidence of transcription in at least one lifestage (≥10 reads in at least one rep). Similarly for the host, ~2.8 billion Illumina RNAseq read pairs uniquely mapping to the A. thaliana genome show that ~77% of the 32,548 gene models have good evidence of transcription in at least one stage (≥10 reads in at least one rep, even though we only sampled roots). A principal component analysis of the host and parasite gene expression data offers several insights into the parasitic process. Principle component 1 (60% of the variance) and 2 (19% of the variance) of the parasite recapitulate the life cycle in PCA space (Fig. 1B). The 12 dpi female transcriptome is more similar to the 24 dpi female transcriptome than to the 12 dpi male transcriptome. Principle components 1 (75% of the variance) and 2 (10% of the variance) of the host show that the greatest difference between infected and uninfected plant tissue is at the early time points (10 hpi), and that the transcriptomes of infected and uninfected plant material converge over time, possibly due to systemic effects of infection. A 12 dpi male syncytium transcriptome is roughly intermediate between a control root transcriptome and a 12 dpi female syncytium transcriptome. Given that at this stage most if not all of the males will have ceased feeding, this could be due to inadequate formation of the feeding site, or regression of the tissue. In any case, by comparing both principal component analyses, we can see that what is a relatively small difference in the transcriptomes of the feeding sites of males and females is amplified to a relatively large difference in the transcriptomes of the males and females themselves (Fig. 1B).The consequences, and possible causes, of large-scale segmental duplication in the Heterodera lineageTo understand the evolutionary origin(s) of the relatively large number of genes in H. schachtii in particular, and Heterodera spp. in general, we analysed the abundance and categories of gene duplication in the predicted exome. Compared to a related cyst nematode, Globodera pallida (derived using comparable methodology and of comparable contiguity) the exomes of H. schachtii and H. glycines are characterised by a relatively smaller proportion of single-copy genes (as classified by MCSanX toolkit17, and a relatively greater proportion of segmental duplications (at least five co-linear genes with no >25 genes between them), with relatively similar proportions of dispersed duplications (two similar genes with >20 other genes between them), proximal duplications (two similar genes with  +0.5 or  More

Dorning, M. & Cipollini, D. Leaf and root extracts of the invasive shrub, Lonicera maackii, inhibit seed germination of three herbs with no autotoxic effects. Plant Ecol. 184, 287–296 (2006).Article 

Google Scholar 
Greer, M. J., Wilson, G. W., Hickman, K. R. & Wilson, S. M. Experimental evidence that invasive grasses use allelopathic biochemicals as a potential mechanism for invasion: Chemical warfare in nature. Plant Soil 385, 165–179 (2014).Article 
CAS 

Google Scholar 
Möhler, H., Diekötter, T., Herrmann, J. D. & Donath, T. W. Allelopathic vs. autotoxic potential of a grassland weed-evidence from a seed germination experiment. Plant Ecol. Divers. 11, 539–549 (2018).Article 

Google Scholar 
Callaway, R. M. & Aschehoug, E. T. Invasive plants versus their new and old neighbors: A mechanism for exotic invasion. Science 290, 521–523 (2000).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Niu, H. B., Liu, W. X., Wan, F. H. & Liu, B. An invasive aster (Ageratina adenophora) invades and dominates forest understories in China: Altered soil microbial communities facilitate the invader and inhibit natives. Plant Soil 294, 73–85 (2007).Article 
CAS 

Google Scholar 
Wardle, D. A., Karban, R. & Callaway, R. M. The ecosystem and evolutionary contexts of allelopathy. Trends Ecol. Evol. 26, 655–662 (2011).Article 
PubMed 

Google Scholar 
Meiners, S. J., Kong, C. H., Ladwig, L. M., Pisula, N. L. & Lang, K. A. Developing an ecological context for allelopathy. Plant Ecol. 213, 1221–1227 (2012).Article 

Google Scholar 
Liebhold, A. M., Brockerhoff, E. G., Kalisz, S., Nunez, M. A. & Wardle, D. A. Biological invasions in forest ecosystems. Biol. Invasions 19, 3437–3458 (2017).Article 

Google Scholar 
Liao, H. X. et al. Soil microbes regulate forest succession in a subtropical ecosystem in China: Evidence from a mesocosm experiment. Plant Soil 430, 277–289 (2018).Article 
CAS 

Google Scholar 
Wardle, D. A., Nilsson, M. C., Gallet, C. & Zackrisson, O. An ecosystem-level perspective of allelopathy. Biol. Rev. 73, 305–319 (2010).Article 

Google Scholar 
Hierro, J. L. & Callaway, R. M. Allelopathy and exotic plant invasion. Plant Soil 256, 29–39 (2003).Article 
CAS 

Google Scholar 
Uddin, M. N., Robinson, R. W., Buultjens, A., Harun, M. A. & Shampa, S. H. Role of allelopathy of Phragmites australis in its invasion processes. J. Exp. Mar. Biol. Ecol. 486, 237–244 (2017).Article 

Google Scholar 
Thiébaut, G., Tarayre, M. & Rodríguez-Pérez, H. Allelopathic effects of native versus invasive plants on one major invader. Front. Plant Sci. 2, 854 (2019).Article 

Google Scholar 
Smith, M., Cecchi, L., Skjøth, C. A., Karrer, G. & Šikoparijae, B. Common ragweed: A threat to environmental health in Europe. Environ. Int. 61, 115–126 (2013).Article 
CAS 
PubMed 

Google Scholar 
Montagnani, C., Gentili, R., Smith, M., Guarino, M. F. & Citterio, S. The worldwide spread, success, and impact of ragweed (Ambrosia spp.). Crit. Rev. Plant Sci. 36, 1–40 (2017).Article 

Google Scholar 
Zeng, K., Zhu, Y. Q. & Liu, J. X. Research progress on ragweed (Ambrosia). Acta Prataculturae Sin. 19, 212–219 (2010).
Google Scholar 
Jacobs, R. L. et al. Responses to ragweed pollen in a pollen challenge chamber versus seasonal exposure identify allergic rhinoconjunctivitis endotypes. J. Allergy Clin. Immun. 130, 122-127.e8 (2012).Article 
PubMed 

Google Scholar 
Lake, R. I. et al. Climate change and future pollen allergy in Europe. Environ. Health Perspect. 125, 385–391 (2017).Article 
PubMed 

Google Scholar 
Wang, J. J., Zhao, B. Y., Li, M. T. & Li, R. Ecological invasion plant-bitter weed (Ambrosia artemisiifolia) and integrated control strategy. Pratacultural Sci. 023, 71–75 (2006).CAS 

Google Scholar 
Deng, Z. Z., Bai, J. D., Zhao, C. Y. & Li, J. S. Advance in invasion mechanisms of Ambrosia artemisiifolia. Pratacultural Sci. 32, 54–63 (2015).
Google Scholar 
Dong, H. G. et al. Diffusion and intrusion features of Ambrosia artemisiifolia and Ambrosia trifida in Yili River Valley. J. Arid Land Resour. Environ. 31, 175–180 (2017).
Google Scholar 
Vink, J. P. et al. Glyphosate-resistant giant ragweed (Ambrosia trifida) control in dicamba-tolerant soybean. Weed Technol. 26, 422–428 (2012).Article 
CAS 

Google Scholar 
Simard, M. J. & Benoit, D. L. Effect of repetitive mowing on common ragweed (Ambrosia artemisiifolia L.) pollen and seed production. Ann. Agric. Environ. Med. 18, 55–62 (2011).PubMed 

Google Scholar 
Goplen, J. J. et al. Seedbank depletion and emergence patterns of giant ragweed (Ambrosia trifida) in Minnesota cropping systems. Weed Sci. 65, 52–60 (2017).Article 

Google Scholar 
Jurik, T. W. Population distributions of plant size and light environment of giant ragweed (Ambrosia trifida L.) at three densities. Oecologia 87, 539–550 (1991).Article 
ADS 
PubMed 

Google Scholar 
Patracchini, C., Vidotto, F. & Ferrero, A. Common ragweed (Ambrosia artemisiifolia) growth as affected by plant density and clipping. Weed Technol. 25, 268–276 (2011).Article 

Google Scholar 
Kazinczi, G. Ragweed seed bank in the soils of arable fields of Transdanubia, Hungary. Hung. Weed Res. Technol. 19(1), 21–36 (2018).
Google Scholar 
Essl, F. et al. Biological flora of the British Isles: Ambrosia artemisiifolia. J. Ecol. 103, 1069–1098 (2015).Article 

Google Scholar 
Goplen, J. J. Giant Ragweed (Ambrosia trifida) Seed Bank Dynamics and Management. (Master’s dissertation, University of Minnesota.) Retrieved from https://hdl.handle.net11299174767 (2015).Yoda, K. Self-thinning in overcrowded pure stands under cultivated and natural conditions. J. Biol. 14, 107–129 (1963).
Google Scholar 
Friedman, J. & Waller, G. R. Allelopathy and autotoxicity. Trends Biochem. Sci. 10, 47–50 (1985).Article 
CAS 

Google Scholar 
Weller, D. E. The interspecific size-density relationship among crowded plant stands and its implications for the −3/2 power rule of self-thinning. Am. Nat. 133, 20–41 (1989).Article 

Google Scholar 
Deng, J. et al. Autotoxicity of phthalate esters in tobacco root exudates: Effects on seed germination and seedling growth. Pedosphere 27, 1073–1082 (2017).Article 
CAS 

Google Scholar 
Sudatti, D. B., Duarte, H. M., Soares, A. R., Salgado, L. T. & Pereira, R. C. New ecological role of seaweed secondary metabolites as autotoxic and allelopathic. Front. Plant Sci. 11, 347 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Singh, H. P., Batish, D. & Kohil, R. Autotoxicity: Concepts, organisms, and ecological significance. Plant Sci. 18, 757–772 (1999).CAS 

Google Scholar 
Chon, S. U. et al. Effects of alfalfa leaf extracts and phenolic allelochemicals on early seedling growth and root morphology of alfalfa and barnyard grass. Crop Prot. 21, 1077–1082 (2002).Article 
CAS 

Google Scholar 
Chen, B. M., D’Antonio, C. M., Molinari, N. & Peng, S. L. Mechanisms of influence of invasive grass litter on germination and growth of coexisting species in California. Biol. Invasions 20, 1881–1897 (2018).Article 

Google Scholar 
Chen, L. C., Wang, S. L., Wang, P. & Kong, C. H. Autoinhibition and soil allelochemical (cyclic dipeptide) levels in replanted Chinese fir (Cunninghamia lanceolata) plantations. Plant Soil 374, 793–801 (2014).Article 
CAS 

Google Scholar 
Perry, L. G. et al. Retracted: Dual role for an allelochemical: catechin from Centaurea maculosa root exudates regulates conspecific seedling establishment. J. Ecol. 93, 1126–1135 (2005).Article 
CAS 

Google Scholar 
Yu, J. Q., Ye, S. F., Zhang, M. F. & Hu, W. H. Effects of root exudates and aqueous root extracts of cucumber (Cucumis sativus) and allelochemicals, on photosynthesis and antioxidant enzymes in cucumber. Biochem. Syst. Ecol. 31, 129–139 (2003).Article 
CAS 

Google Scholar 
Kong, C. H., Wang, P. & Xu, X. H. Allelopathic interference of Ambrosia trifida with wheat (Triticum aestivum). Agric. Ecosyst. Environ. 119, 416–420 (2007).Article 
CAS 

Google Scholar 
Béres, I., Kazinczi, G. & Narwal, S. S. Allellopathic plants. 4. Common ragweed (Ambrosia elatior L. syn. A. artemisiifolia). Allelopathy J. 9, 27–34 (2002).
Google Scholar 
Bauer, J. T., Shannon, S. M., Stoops, R. E. & Reynolds, H. L. Context dependency of the allelopathic effects of Lonicera maackii on seed germination. Plant Ecol. 213, 1907–1916 (2012).Article 

Google Scholar 
Renne, I. J., Sinn, B. T., Shook, G. W., Sedlacko, D. M. & Hierro, J. L. Eavesdropping in plants: Delayed germination via biochemical recognition. J. Ecol. 102, 86–94 (2014).Article 

Google Scholar 
Loydi, A., Donath, T. W., Eckstein, R. L. & Otte, A. Non-native species litter reduces germination and growth of resident forbs and grasses: Allelopathic, osmotic or mechanical effects?. Biol. Invasions 17, 581–595 (2014).Article 

Google Scholar 
Bais, H. P., Weir, T. L., Perry, L. G., Gilroy, S. & Vivanco, J. M. The role of root exudates in rhizosphere interactions with plants and other organisms. Annu. Rev. Plant Biol. 57, 233–266 (2006).Article 
CAS 
PubMed 

Google Scholar 
Bonea, D., Bonciu, E., Niculescu, M. & Olaru, A. L. The allelopathic, cytotoxic and genotoxic effect of Ambrosia artemisiifolia on the germination and root meristems of Zea mays. Caryologia 71, 24–28 (2017).Article 

Google Scholar 
Dadkhah, A. Allelopathic effect of sugar beet (Beta vulgaris) and eucalyptus (Eucalyptus camaldulensis) on seed germination and growth of Portulaca oleracea. Russ. Agric. Sci. 39, 117–123 (2013).Article 

Google Scholar 
Zheng, L. & Feng, Y. L. Allelopathic effects of Eupatorium adenophorum Spreng on. seed germination and seedling growth in ten herbaceous species. Acta Ecol. Sin. 25, 2782–2787 (2005).CAS 

Google Scholar 
Brückner, D. J. The allelopathic effect of ragweed (Ambrosia artemisiifolia L.) on the germination of cultivated plants. Novenytermeles 47, 635–644 (1998).
Google Scholar 
Qin, R. M. et al. The evolution of increased competitive ability, innate competitive advantages, and novel biochemical weapons act in concert for a tropical invader. New Phytol. 197, 979–988 (2012).Article 
PubMed 

Google Scholar 
Zheng, Y. L. et al. Integrating novel chemical weapons and evolutionarily increased competitive ability in success of a tropical invader. New Phytol. 205, 1350–1359 (2015).Article 
PubMed 

Google Scholar 
Kaushal, R., Verma, K. S. & Singh, K. N. Effect of Grewia optiva and Populus deltoides leachatesv on field crops. Allelopathy J. 11, 229–234 (2003).
Google Scholar 
Kumari, A. & Kohli, R. Autotoxicity of ragweed parthenium (Parthenium hysterophorus). Weed Sci. 35, 629–632 (1987).Article 

Google Scholar 
Einhellig, F. A. Allelopathy: Current status and future goals. In Allelopathy: Organisms, processes and applications (ed. Inderjit Dakshini, K. M. M.) 1–24 (Am Chem. Soc, Washington, 1995).
Google Scholar 
Hadack, F. Secondary metabolites as plant traits: Current assessment and future perspectives. Crit. Rev. Plant Sci. 21, 273–322 (2002).Article 

Google Scholar 
Rice, E. L. Biological Control of Weeds and Plant Diseases (Oklahomka Press, 1995).
Google Scholar 
Choi, B. et al. Common ragweed-derived phenolic compounds and their effects on germination and seedling growth of weed species. Weed Turfgrass Sci. 30, 396–404 (2010).
Google Scholar 
Friedman, J. & Waller, G. R. Seeds as allelopathic agents. Chem. Ecol. 9, 1107–1117 (1983).Article 
CAS 

Google Scholar 
Canals, R. M., Emeterio, L. S. & Peralta, J. Autotoxicity in Lolium rigidum: Analyzing the role of chemically mediated interactions in annual plant populations. J. Theor. Biol. 235, 402–407 (2005).Article 
ADS 
CAS 
PubMed 

Google Scholar 
San Emeterio, L., Damgaard, C. & Canals, R. M. Modelling the combined effect of chemical interference and resource competition on the individual growth of two herbaceous populations. Plant Soil 292, 95–103 (2007).Article 
CAS 

Google Scholar 
Dickerson, C. T. Studies on the germination, growth, development and control of Common Ragweed (Ambrosia artemisiifolia L.). PhD thesis, Cornell University, Ann Arbor (1968).Nuutinen, V. & Butt, K. R. Homing ability widens the sphere of influence of the earthworm Lumbricus terrestris L. Soil Biol. Biochem. 37, 805–807 (2005).Article 
CAS 

Google Scholar 
Favaretto, A., Scheffer-basso, S. M. & Perez, N. B. Autotoxicity in tough lovegrass (Eragrostis plana). Planta Daninha 35(35), e017164046 (2017).
Google Scholar 
Sinkkonen, A. Modelling the effect of autotoxicity on density-dependent phytotoxicity. J. Theor. Biol. 244, 218–227 (2007).Article 
ADS 
MathSciNet 
CAS 
PubMed 
MATH 

Google Scholar 
Zhang, S. S., Shi, F. Q., Yang, W. Z., Xiang, Z. Y. & Duan, Z. L. Autotoxicity as a cause for natural regeneration failure in Nyssa yunnanensis and its implications for conservation. Isr. J. Plant Sci. 62, 187–197 (2015).Article 

Google Scholar 
Liu, Y. et al. Relationship between seed germination and invasion of Ambrosia artemisiifolia and A. trifida at different positions. Acta Ecol. Sin. 39, 9079–9088 (2019).

Google Scholar  More

Wurst, S., Sonnemann, I. & Zaller, J. G. Soil Macro-Invertebrates: Their Impact on Plants and Associated Aboveground Communities in Temperate Regions. in Aboveground-Belowground Community Ecology. Ecological Studies 234 (eds. Ohgushi, T., Wurst, S. & Johnson, S. N.) 175–200 (Springer International Publishing, 2018).Wagg, C., Bender, S. F., Widmer, F. & van der Heijden, M. G. A. Soil biodiversity and soil community composition determine ecosystem multifunctionality. PNAS 111, 5266–5270 (2014).ADS 
CAS 
Article 

Google Scholar 
Potapov, A. M. et al. Feeding habits and multifunctional classification of soil-associated consumers from protists to vertebrates. Biol. Rev. 97, 1057–1117 (2022).Article 

Google Scholar 
Potapov, A., Tiunov, A. V. & Scheu, S. Uncovering trophic positions and food resources of soil animals using bulk natural stable isotope composition. Biol. Rev. 94, 37–59 (2019).Article 

Google Scholar 
Voroney, R. P. & Heck, R. J. The Soil Habitat. Soil Microbiology, Ecology and Biochemistry (Elsevier Inc., 2015).Wardle, D. A. The influence of biotic interactions on soil biodiversity. Ecol. Lett. 9, 870–886 (2006).Article 

Google Scholar 
De Deyn, G. B. & Kooistra, L. The role of soils in habitat creation, maintenance and restoration. Philos. Trans. R. Soc. Lond., B, Biol. Sci. 376 (2021).Phillips, H. R. P. et al. Global distribution of earthworm diversity. Science 366, 480–485 (2019).ADS 
CAS 
Article 

Google Scholar 
Geitner, C. et al. Soil and Land Use in the Alps—Challenges and Examples of Soil-Survey and Soil-Data Use to Support Sustainable Development. In Soil Mapping and Process Modelling for Sustainable Land Use Management 221–292 (Elsevier, 2017).FAO. Understanding Mountain Soils: A contribution from mountain areas to the International Year of Soils 2015. (2015).Praeg, N. et al. The role of land management and elevation in shaping soil microbial communities: Insights from the Central European Alps. Soil Biol. Biochem. 150, 107951 (2020).CAS 
Article 

Google Scholar 
Somme, L. Adaptations of terrestrial arthropods to the alpine environment. Biol. Rev. 64, 367–407 (1989).Article 

Google Scholar 
Meyer, E. & Thaler, K. Animal Diversity at High Altitudes in the Austrian Central Alps. In Arctic and Alpine Biodiversity (eds. Chapin, F. S. & Körner, C.) (Springer Berlin Heidelberg, 1995).Guerra, C. A. et al. Blind spots in global soil biodiversity and ecosystem function research. Nat. Commun. 11, 1–13 (2020).Article 

Google Scholar 
Kempson, D., Lloyd, M. & Ghelardi, R. A new extractor for woodland litter. Pedobiologia 3, 1–30 (1963).
Google Scholar 
Tasser, E., Ruffini, F. V. & Tappeiner, U. An integrative approach for analysing landscape dynamics in diverse cultivated and natural mountain areas. Landsc. Ecol. 24, 611–628 (2009).Article 

Google Scholar 
European Environment Agency, (EAA). European Union, Copernicus Land Monitoring Service. https://land.copernicus.eu/pan-european/corine-land-cover/clc2018 (2018).IUSS Working Group WRB. World Reference Base for Soil Resources 2014, update 2015 International soil classification system for naming soils and creating legends for soil maps. World Soil Resources Reports No. 106 (2015).Christian, E. & Zicsi, A. A synoptic key to the earthworms of Austria (Oligochaeta: Lumbricidae). Die Bodenkultur 50, 121–131 (1999).Czusdi, C. & Zicsi, A. Earthworms of Hungary (Annelida, Oligochaeta, Lumbricidae). (Hungarian Natural History Museum, 2003).Schaefer, M. Brohmer – Fauna von Deutschland: Ein Bestimmungsbuch unserer heimischen Tierwelt. (Quelle & Meyer Verlag, 2009).Klausnitzer, B. Exkursionsfauna von Deutschland. (Springer Berlin Heidelberg, 2011).Ellis, D. Taxonomic sufficiency in pollution assessment. Mar. Pollut. Bull. 16, 459 (1985).Article 

Google Scholar 
Caruso, T. & Migliorini, M. Micro-arthropod communities under human disturbance: is taxonomic aggregation a valuable tool for detecting multivariate change? Evidence from Mediterranean soil oribatid coenoses. Acta Oecol. 30, 46–53 (2006).ADS 
Article 

Google Scholar 
Steinwandter, M., Schlick-Steiner, B. C., Seeber, G. U. H., Steiner, F. M. & Seeber, J. Effects of Alpine land-use changes: Soil macrofauna community revisited. Ecol. Evol. 7, 5389–5399 (2017).Article 

Google Scholar 
de Oliveira, S. S. et al. Higher taxa are sufficient to represent biodiversity patterns. Ecol. Indic. 111, 105994 (2020).Article 

Google Scholar 
Lavelle, P. et al. Soil macroinvertebrate communities: A world‐wide assessment. Glob. Ecol. Biogeogr. 31, 1261–1276 (2022).Article 

Google Scholar 
Pik, A. J., Oliver, I. & Beattie, A. J. Taxonomic sufficiency in ecological studies of terrestrial invertebrates. Austral. Ecol. 24, 555–562 (1999).Article 

Google Scholar 
Parisi, V., Menta, C., Gardi, C., Jacomini, C. & Mozzanica, E. Microarthropod communities as a tool to assess soil quality and biodiversity: A new approach in Italy. Agric. Ecosyst. Environ. 105, 323–333 (2005).Article 

Google Scholar 
Ruiz, N. et al. IBQS: A synthetic index of soil quality based on soil macro-invertebrate communities. Soil Biol. Biochem. 43, 2032–2045 (2011).
Google Scholar 
Seeber, J. et al. A 30-years collection of soil macro-invertebrate abundance data from the European Alps. PANGAEA https://doi.org/10.1594/PANGAEA.944405 (2022).de Jong, Y. et al. Fauna Europaea – All European animal species on the web. Biodivers. Data J. 2 (2014).Shannon, C. E. A Mathematical Theory of Communication. Bell Syst. Tech. J. 27, 379–423 (1948).MathSciNet 
Article 

Google Scholar 
Steinberger, K.-H. & Meyer, E. The Spider Fauna of the Nature Reserve Rheindelta (Vorarlberg, Western Austria) (Arachnida: Araneae). Ber. nat.-med. Verein Innsbruck 82, 195–215 (1995).
Google Scholar 
Kopeszki, H. & Meyer, E. Species Composition and Abundance of Collembola in Forest Soils in the Provinces of Bozen and Trient (Italy). Ber. nat.-med. Verein Innsbruck 83, 221–237 (1996).
Google Scholar 
Geitner, C., Mätzler, A., Bou-Vinals, A., Meyer, E. & Tusch, M. Soil characteristics and colonization by earthworms (Lumbricidae) on pastures and hay meadows in the Brixenbach Valley (Kitzbühel Alps, Tyrol). Die Bodenkultur 65, 39–51 (2014).
Google Scholar 
Meyer, E. & Steinberger, K.-H. Über die Bodenfauna in Wäldern Vorarlbergs (Österreich) Bestand und Auswirkungen von Gesteinsmehlapplikationen. Verhandlungen der Gesellschaft für Ökologie 23, 149–164 (1994).
Google Scholar 
Peham, T. & Meyer, E. Kommentierte Artenlisten ausgewählter Bodentiergruppen aus der Erhebung des SoilDiv-Projektes in Südtirol. Gredleriana 14, 227–262 (2014).
Google Scholar 
Steinwandter, M., Rief, A., Scheu, S., Traugott, M. & Seeber, J. Structural and functional characteristics of high alpine soil macro-invertebrate communities. Eur. J. Soil Biol. 86, 72–80 (2018).Article 

Google Scholar 
Steinwandter, M. et al. Does green manuring positively affect the soil macro-invertebrates in vineyards? A case study from Kaltern/Caldaro (South Tyrol, Italy). Gredleriana 18, 17–26 (2018).
Google Scholar 
Steinwandter, M. et al. Raw data from: Does green manuring positively affect the soil macro-invertebrates in vineyards? A case study from Kaltern/ Caldaro (South Tyrol, Italy). PANGAEA https://doi.org/10.1594/PANGAEA.900632 (2019).Damisch, K., Steinwandter, M., Tappeiner, U. & Seeber, J. Soil Macroinvertebrate Distribution along a Subalpine Land Use Transect. Mt. Res. Dev. 40, R1–R10 (2020).Article 

Google Scholar 
Damisch, K., Steinwandter, M., Tappeiner, U. & Seeber, J. Abundance data from soil macro-invertebrates along a subalpine land-use transect. PANGAEA https://doi.org/10.1594/PANGAEA.918958 (2020).Schneider, E., Steinwandter, M. & Seeber, J. A comparison of Alpine soil macro-invertebrate communities from European larch and Swiss pine forests in the LTSER area “Val Mazia/Matschertal”, South Tyrol. Gredleriana 19, 217–228 (2019).
Google Scholar 
Schneider, E., Steinwandter, M. & Seeber, J. Raw data from: A comparison of Alpine soil macro-invertebrate communities from European larch and Swiss pine forests in the LTSER area “Val Mazia/Matschertal”, South Tyrol. PANGAEA https://doi.org/10.1594/PANGAEA.910666 (2020).Seeber, J. et al. Soil invertebrate abundance, diversity, and community composition across steep high elevation snowmelt gradients in the European Alps. Arct. Antarct. Alp. Res. 53, 288–299 (2021).Article 

Google Scholar 
Seeber, J. et al. Abundance data from soil macro- and mesofauna of alpine snowbeds in the European Alps (summer 2017). PANGAEA https://doi.org/10.1594/PANGAEA.935737 (2021).Guariento, E. et al. Management intensification of hay meadows and fruit orchards alters soil macro-invertebrate communities differently. Agronomy 10, 767 (2020).CAS 
Article 

Google Scholar  More

Similar non-native and invasive flora, such as the fever tree (pictured) are found in regions previously occupied by the same European empire.Credit: Alamy

In 1860, a British expedition raided the highland forests of South America, looking for a hot commodity: Cinchona seeds. The bark of these ‘fever’ trees produces the anti-malarial compound quinine, and the British Empire sought a stable source of the drug for its soldiers and civil service in India. After cultivation in the United Kingdom, young Cinchona trees were planted across southern India and what is now Sri Lanka.The British quinine scheme failed — instead, a species introduced to Java, now part of Indonesia, by the Dutch Empire later dominated the global market — but Cinchona trees are still common in parts of India.Such botanical legacies of imperial rule are common, finds a study published on 17 October in Nature Ecology & Evolution1. Regions that were once occupied by the same European colonial power — such as India and Sri Lanka — tend to have similar species of non-native and invasive plants. The longer the regions were occupied, the more their populations of invasive species resemble each other, the research found.Alien floraThe link between European colonialism and invasive species is intuitive, and has been noted by other researchers, says Bernd Lenzner, a macro-ecologist at the University of Vienna who led the study. To test the association, his team turned to the Global Naturalized Alien Flora database, which maps the distribution of nearly 14,000 invasive plant species.
The imperial roots of climate science
Across more than 1,100 regions, including 404 islands, the researchers found that regions once occupied by the British Empire had more similarities in their invasive flora than did ‘artificial’ empires that the team assembled from random regions. This was also the case for regions once part of the Dutch Empire (former Spanish and Portuguese colonies had alien-plant compositions similar to those of the artificial empires).Climate and geography play an important part in explaining the overlap in the diversity of invasive species, modelling by Lenzner’s team found, but so does the length of time regions were occupied by an imperial power. Regions that were central to trade, such as southern India for the British Empire and Indonesia for the Dutch Empire, formed clusters with considerable overlap in invasive-plant composition.The analysis did not look at when individual plant species were introduced or why. But anecdotally, many of the plants that were commonly taken to former empires were once of economic value and their populations were probably established on purpose, says Lenzner.Global trade impactsThe study’s conclusions might be “super obvious”, but they have important implications for conservation, says Nussaïbah Raja, a palaeontologist at Friedrich-Alexander University of Erlangen–Nürnberg in Erlangen, Germany. “We should be taking this history into consideration when we think about management of species.” Appreciating the history of introduced plants — as well as their place in today’s ecosystems — could help conservationists to handle future changes in biodiversity, such as those driven by climate change, Raja adds.Global trade is beginning to overwrite the colonial legacy of introduced plants. For example, the analysis showed similarities between invasive plant populations in Fujian, China, and some parts of Australia. Although both places were once connected by the British Empire, more recent global trade might also be partly responsible for the overlap.“We are still seeing these imprints of the colonial-empire legacies from centuries ago,” Lenzner says. “So what we’re doing and the species we’re redistributing today will be visible far into the future.” More

Kelemen PB, Manning CE. Reevaluating carbon fluxes in subduction zones, what goes down, mostly comes up. Proc Natl Acad Sci USA. 2015;112:E3997–4006.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Vitale Brovarone A, Sverjensky DA, Piccoli F, Ressico F, Giovannelli D, Daniel I. Subduction hides high-pressure sources of energy that may feed the deep subsurface biosphere. Nat Commun. 2020;11:1–1.Article 

Google Scholar 
Harris RN, Wang K. Thermal models of the middle America trench at the Nicoya Peninsula, Costa Rica. Geophys Res Lett. 2002;29:6–1.Article 

Google Scholar 
Plümper O, King HE, Geisler T, Liu Y, Pabst S, Savov IP, et al. Subduction zone forearc serpentinites as incubators for deep microbial life. Proc Natl Acad Sci USA. 2017;114:4324–9.PubMed 
PubMed Central 
Article 

Google Scholar 
Lee H, Fischer TP, de Moor JM, Sharp ZD, Takahata N, Sano Y. Nitrogen recycling at the Costa Rican subduction zone: the role of incoming plate structure. Sci Rep. 2017;7:1–10.
Google Scholar 
Stern RJ. Subduction zones. Rev Geophys. 2002;40:3–38.Article 

Google Scholar 
Fullerton KM, Schrenk MO, Yücel M, Manini E, Basili M, Rogers TJ, et al. Effect of tectonic processes on biosphere–geosphere feedbacks across a convergent margin. Nat Geosci. 2021;14:301–6.CAS 
Article 

Google Scholar 
Barry PH, de Moor JM, Giovannelli D, Schrenk M, Hummer DR, Lopez T, et al. Forearc carbon sink reduces long-term volatile recycling into the mantle. Nature. 2019;568:487–92.CAS 
PubMed 
Article 

Google Scholar 
Moore EK, Jelen BI, Giovannelli D, Raanan H, Falkowski PG. Metal availability and the expanding network of microbial metabolisms in the Archaean eon. Nat Geosci. 2017;10:629–36.CAS 
Article 

Google Scholar 
Barnes JD, Cullen J, Barker S, Agostini S, Penniston-Dorland S, Lassiter JC, et al. The role of the upper plate in controlling fluid-mobile element (Cl, Li, B) cycling through subduction zones: Hikurangi forearc, New Zealand. Geosphere. 2019;15:642–58.Article 

Google Scholar 
Clift P, Vannucchi P. Controls on tectonic accretion versus erosion in subduction zones: Implications for the origin and recycling of the continental crust. Rev Geophys. 2004;42:1–31.Article 

Google Scholar 
Rüpke LH, Morgan JP, Hort M, Connolly JA. Serpentine and the subduction zone water cycle. Earth Planet Sci Lett. 2004;223:17–34.Article 

Google Scholar 
Carr MJ, Feigenson MD, Bennett EA. Incompatible element and isotopic evidence for tectonic control of source mixing and melt extraction along the Central American arc. Contrib Miner Pet. 1990;105:369–80.CAS 
Article 

Google Scholar 
Gazel E, Carr MJ, Hoernle K, Feigenson MD, Szymanski D, Hauff F, et al. Galapagos‐OIB signature in southern Central America: mantle refertilization by arc–hot spot interaction. Geochem Geophys Geosyst. 2009;10:1–32.Article 

Google Scholar 
Trembath-Reichert E, Butterfield DA, Huber JA. Active subseafloor microbial communities from Mariana back-arc venting fluids share metabolic strategies across different thermal niches and taxa. ISME J. 2019;13:2264–79. https://doi.org/10.1038/s41396-019-0431-y.Power JF, Carere CR, Lee CK, Wakerley GL, Evans DW, Button M, et al. Microbial biogeography of 925 geothermal springs in New Zealand. Nat Commun. 2018;9:1–2.CAS 
Article 

Google Scholar 
Acocella V, Spinks K, Cole J, Nicol A. Oblique back arc rifting of Taupo Volcanic zone. NZ Tecton. 2003;22:1–18.
Google Scholar 
Curtis AC, Wheat CG, Fryer P, Moyer CL. Mariana forearc serpentinite mud volcanoes harbor novel communities of extremophilic archaea. Geomicrobiol J. 2013;30:430–41.Article 

Google Scholar 
Inskeep WP, Jay ZJ, Herrgard MJ, Kozubal MA, Rusch DB, Tringe SG, et al. Phylogenetic and functional analysis of metagenome sequence from high-temperature archaeal habitats demonstrate linkages between metabolic potential and geochemistry. Front Microbiol. 2013;4:1–21.Article 

Google Scholar 
Colman DR, Lindsay MR, Amenabar MJ, Boyd ES. The intersection of geology, geochemistry, and microbiology in continental hydrothermal systems. Astrobiology. 2019;19:1505–22.CAS 
PubMed 
Article 

Google Scholar 
Inskeep WP, Jay ZJ, Tringe SG, Herrgård MJ, Rusch DB, YNP Metagenome Project Steering Committee and Working Group Members. The YNP metagenome project: environmental parameters responsible for microbial distribution in the Yellowstone geothermal ecosystem. Front Microbiol. 2013;4:1–15.Article 

Google Scholar 
Hou W, Wang S, Dong H, Jiang H, Briggs BR, Peacock JP, et al. A comprehensive census of microbial diversity in hot springs of Tengchong, Yunnan Province China using 16S rRNA gene pyrosequencing. PloS One. 2013;8:1–15.
Google Scholar 
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017;27:824–34.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bowers RM, Kyrpides NC, Stepanauskas R, Harmon-Smith M, Doud D, Reddy TB, et al. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotechnol. 2017;35:725–31.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Probst AJ, Castelle CJ, Singh A, Brown CT, Anantharaman K, Sharon I, et al. Genomic resolution of a cold subsurface aquifer community provides metabolic insights for novel microbes adapted to high CO2 concentrations. Environ Microbiol. 2017;19:459–74.CAS 
PubMed 
Article 

Google Scholar 
Probst AJ, Ladd B, Jarett JK, Geller-McGrath DE, Sieber CM, Emerson JB, et al. Differential depth distribution of microbial function and putative symbionts through sediment-hosted aquifers in the deep terrestrial subsurface. Nat Microbiol. 2018;3:328–36.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
He C, Keren R, Whittaker M, Farag IF, Doudna J, Cate JH, et al. Genome-resoled metagenomics reveals site-specific diversity of episymbiotic CPR bacteria and DPANN archaea in groundwater ecosystems. Nat. Microbiol. 2021;6:354–65.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Grettenberger CL, Hamilton TL. Metagenome-assembled genomes of novel taxa from an acid mine drainage environment. Appl Environ Microbiol. 2021;87:e0077221. https://doi.org/10.1101/2020.07.02.185728.Uritskiy GV, DiRuggiero J, Taylor J. MetaWRAP–a flexible pipeline for genome-resolved metagenomic data analysis. Microbiome. 2018;6:1–3.Article 

Google Scholar 
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30:2068–9.CAS 
PubMed 
Article 

Google Scholar 
Garber AI, Nealson KH, Okamoto A, McAllister SM, Chan CS, Barco RA, et al. FeGenie: a comprehensive tool for the identification of iron genes and iron gene neighborhoods in genome and metagenome assemblies. Front Microbiol. 2020;11:37. https://doi.org/10.3389/fmicb.2020.00037.Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J Mol Biol. 2016;428:726–31.CAS 
PubMed 
Article 

Google Scholar 
Graham ED, Heidelberg JF, Tully BJ. Potential for primary productivity in a globally distributed bacterial phototroph. ISME J. 2018;12:1861–6.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Chaumeil PA, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics. 2020;36:1925–27.CAS 

Google Scholar 
Berg IA, Kockelkorn D, Ramos-Vera WH, Say RF, Zarzycki J, Hügler M, et al. Autotrophic carbon fixation in archaea. Nat Rev Microbiol. 2010;8:447–60.CAS 
PubMed 
Article 

Google Scholar 
Berg IA. Ecological aspects of the distribution of different autotrophic CO2 fixation pathways. Appl Environ Microbiol. 2011;77:1925–36.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Youssef NH, Farag IF, Hahn CR, Jarett J, Becraft E, Eloe-Fadrosh E, et al. Genomic characterization of candidate division LCP-89 reveals an atypical cell wall structure, microcompartment production, and dual respiratory and fermentative capacities. Appl Environ Microbiol. 2019;85:1–19.Article 

Google Scholar 
Nigro LM, King GM. Disparate distributions of chemolithotrophs containing form IA or IC large subunit genes for ribulose-1, 5-bisphosphate carboxylase/oxygenase in intertidal marine and littoral lake sediments. FEMS Microbiol Ecol. 2007;60:113–25.CAS 
PubMed 
Article 

Google Scholar 
Aminuddin M, Nicholas DJ. Electron transfer during sulphide and sulphite oxidation in Thiobacillus denitrificans. Microbiology. 1974;82:115–23.
Google Scholar 
Giovannelli D, Sievert SM, Hügler M, Markert S, Becher D, Schweder T, et al. Insight into the evolution of microbial metabolism from the deep-branching bacterium, Thermovibrio ammonificans. eLife. 2017;6:1–31.Article 

Google Scholar 
Nakagawa S, Shataih Z, Banta A, Beveridge TJ, Sako Y, Reysenbach AL. Sulfurihydrogenibium yellowstonense sp. nov., an extremely thermophilic, facultatively heterotrophic, sulfur-oxidizing bacterium from Yellowstone National Park, and emended descriptions of the genus Sulfurihydrogenibium, Sulfurihydrogenibium subterraneum. Int J Syst Evol Microbiol. 2005;55:2263–8.CAS 
PubMed 
Article 

Google Scholar 
Leclerque A, Kleespies RGA. Rickettsiella bacterium from the hard tick, Ixodes woodi: molecular taxonomy combining multilocus sequence typing (MLST) with significance testing. PLoS One. 2012;7:e38062. https://doi.org/10.1371/journal.pone.0038062.Quatrini R, Johnson DB. Acidithiobacillus ferrooxidans. Trends Microbiol. 2019;27:282–3.CAS 
PubMed 
Article 

Google Scholar 
Spang A, Poehlein A, Offre P, Zumbrägel S, Haider S, Rychlik N, et al. The genome of the ammonia‐oxidizing Candidatus Nitrososphaera gargensis: insights into metabolic versatility and environmental adaptations. Environ Microbiol. 2012;14:3122–45.CAS 
PubMed 
Article 

Google Scholar 
Chen CY, Chen PC, Weng FC, Shaw GT, Wang D. Habitat and indigenous gut microbes contribute to the plasticity of gut microbiome in oriental river prawn during rapid environmental change. PLoS One. 2017;12:e0181427. https://doi.org/10.1371/journal.pone.0181427.Garcia R, Müller R. The family Myxococcaceae. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. The prokaryotes: Deltaproteobacteria and Epsilonproteobacteria. Berlin: Springer; 2014. p. 191–212.Garcia R, Müller R. Simulacricoccus ruber gen. nov., sp. nov., a microaerotolerant, non-fruiting, myxospore-forming soil myxobacterium and emended description of the family Myxococcaceae. Int J Syst Evol Microbiol. 2018;68:3101–10.CAS 
PubMed 
Article 

Google Scholar 
Iino T. The family Ignavibacteriaceae. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. The prokaryotes: other major lineages of bacteria and the archaea. New York, NY: Springer Science + Business Media; 2014. p. 701–3.Petrie L, North NN, Dollhopf SL, Balkwill DL, Kostka JE. Enumeration and characterization of iron (III)-reducing microbial communities from acidic subsurface sediments contaminated with uranium (VI). Appl Environ Microbiol. 2003;69:7467–79.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Fincker M, Huber JA, Orphan VJ, Rappé MS, Teske A, Spormann AM. Metabolic strategies of marine subseafloor Chloroflexi inferred from genome reconstructions. Environ Microbiol. 2020;22:3188–204.CAS 
PubMed 
Article 

Google Scholar 
Chen YJ, Leung PM, Wood JL, Bay SK, Hugenholtz P, Kessler AJ, et al. Metabolic flexibility allows bacterial habitat generalists to become dominant in a frequently disturbed ecosystem. ISME J. 2021;15:2986–3004.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Flieder M, Buongiorno J, Herbold CW, Hausmann B, Rattei T, Lloyd KG, et al. Novel taxa of Acidobacteriota implicated in seafloor sulfur cycling. ISME J. 2021;15:3159–80.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kim M, Wilpiszeski RL, Wells M, Wymore AM, Gionfriddo CM, Brooks SC, et al. Metagenome-assembled genome sequences of novel prokaryotic species from the mercury-contaminated East Fork Poplar Creek, Oak Ridge, Tennessee, USA. Microbiol Resour Announc. 2021;10:e00153–21. https://doi.org/10.1128/MRA.00153-21.Santos‐Júnior CD, Logares R, Henrique‐Silva F. Microbial population genomes from the Amazon River reveal possible modulation of the organic matter degradation process in tropical freshwaters. Mol Ecol. 2022;31:206–19.PubMed 
Article 

Google Scholar 
Yamada T, Sekiguchi Y. Cultivation of uncultured Chloroflexi subphyla: significance and ecophysiology of formerly uncultured Chloroflexi ‘subphylum I’ with natural and biotechnological relevance. Microbes Environ. 2009;24:205–16.PubMed 
Article 

Google Scholar 
Sheik CS, Reese BK, Twing KI, Sylvan JB, Grim SL, Schrenk MO, et al. Identification and removal of contaminant sequences from ribosomal gene databases: lessons from the census of deep life. Front Microbiol. 2018;9:840. https://doi.org/10.3389/fmicb.2018.00840.Doughari HJ, Ndakidemi PA, Human IS, Benade S. The ecology, biology and pathogenesis of Acinetobacter spp.: an overview. Microbes Environ. 2011;26:101–12.PubMed 
Article 

Google Scholar 
Han XY, Han FS, Segal J. Chromobacterium haemolyticum sp. nov., a strongly haemolytic species. Int J Syst Evol Microbiol. 2008;58:1398–403.CAS 
PubMed 
Article 

Google Scholar 
Lau MC, Kieft TL, Kuloyo O, Linage-Alvarez B, Van Heerden E, Lindsay MR, et al. An oligotrophic deep-subsurface community dependent on syntrophy is dominated by sulfur-driven autotrophic denitrifiers. Proc Natl Acad Sci USA. 2016;113:E7927–36.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Momper L, Jungbluth SP, Lee MD, Amend JP. Energy and carbon metabolisms in a deep terrestrial subsurface fluid microbial community. ISME J. 2017;11:2319–33.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Worzewski T, Jegen M, Kopp H, Brasse H, Taylor Castillo W. Magnetotelluric image of the fluid cycle in the Costa Rican subduction zone. Nat Geosci. 2011;4:108–11.CAS 
Article 

Google Scholar 
Hensen C, Wallmann K, Schmidt M, Ranero CR, Suess E. Fluid expulsion related to mud extrusion off Costa Rica—a window to the subducting slab. Geology. 2004;32:201–4.CAS 
Article 

Google Scholar 
Simpson DR. Aluminum phosphate variants of feldspar. Am Miner. 1977;62:351–5.CAS 

Google Scholar 
London DA, Cerny P, Loomis J, Pan JJ. Phosphorus in alkali feldspars of rare-element granitic pegmatites. Can Miner. 1990;28:771–86.CAS 

Google Scholar 
Petrillo C, Castaldi S, Lanzilli M, Selci M, Cordone A, Giovannelli D, et al. Genomic and physiological characterization of Bacilli isolated from salt-pans with plant growth promoting features. Front Microbiol. 2021;12:715678. https://doi.org/10.3389/fmicb.2021.715678.Ghiorse WC, Wilson JT. Microbial ecology of the terrestrial subsurface. Adv Appl Microbiol. 1988;33:107–72.CAS 
PubMed 
Article 

Google Scholar 
Barker WW, Welch SA, Chu S, Banfield JF. Experimental observations of the effects of bacteria on aluminosilicate weathering. Am Miner. 1998;83:1551–63.CAS 
Article 

Google Scholar 
Bennett PC, Rogers JR, Choi WJ, Hiebert FK. Silicates, silicate weathering, and microbial ecology. Geomicrobiol J. 2001;18:3–19.CAS 
Article 

Google Scholar 
Hügler M, Sievert SM. Beyond the Calvin cycle: autotrophic carbon fixation in the ocean. Ann Rev Mar Sci. 2011;3:261–89.PubMed 
Article 

Google Scholar 
Markert S, Arndt C, Felbeck H, Becher D, Sievert SM, Hügler M, et al. Physiological proteomics of the uncultured endosymbiont of Riftia pachyptila. Science. 2007;315:247–50.CAS 
PubMed 
Article 

Google Scholar 
Bar-Even A, Noor E, Milo R. A survey of carbon fixation pathways through a quantitative lens. J Exp Bot. 2012;63:2325–42.CAS 
PubMed 
Article 

Google Scholar 
Stevens TO, McKinley JP. Lithoautotrophic microbial ecosystems in deep basalt aquifers. Science. 1995;270:450–5.CAS 
Article 

Google Scholar 
Barker WW, Welch SA, Banfield JF. Biogeochemical weathering of silicate minerals. Rev Miner Geochem. 1997;35:391–428.CAS 

Google Scholar 
Frank YA, Kadnikov VV, Lukina AP, Banks D, Beletsky AV, Mardanov AV, et al. Characterization and genome analysis of the first facultatively alkaliphilic Thermodesulfovibrio isolated from the deep terrestrial subsurface. Front Microbiol. 2016;7:2000. https://doi.org/10.3389/fmicb.2016.02000.Woycheese KM, Meyer-Dombard DA, Cardace D, Argayosa AM, Arcilla CA. Out of the dark: transitional subsurface-to-surface microbial diversity in a terrestrial serpentinizing seep (Manleluag, Pangasinan, the Philippines). Front Microbiol. 2015;6:1–12.Article 

Google Scholar 
Brazelton WJ, Morrill PL, Szponar N, Schrenk MO. Bacterial communities associated with subsurface geochemical processes in continental serpentinite springs. Appl Environ Microbiol. 2013;79:3906–16.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Moser DP, Gihring TM, Brockman FJ, Fredrickson JK, Balkwill DL, Dollhopf ME, et al. Desulfotomaculum and Methanobacterium spp. dominate a 4-to 5-kilometer-deep fault. Appl Environ Microbiol. 2005;71:8773–83.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Schwarzenbach EM, Gill BC, Gazel E, Madrigal P. Sulfur and carbon geochemistry of the Santa Elena peridotites: comparing oceanic and continental processes during peridotite alteration. Lithos. 2016;252:92–108.Article 

Google Scholar 
Sánchez‐Murillo R, Gazel E, Schwarzenbach EM, Crespo‐Medina M, Schrenk MO, Boll J, et al. Geochemical evidence for active tropical serpentinization in the Santa Elena Ophiolite, Costa Rica: an analog of a humid early Earth? Geochem Geophys Geosyst. 2014;15:1783–800.Article 

Google Scholar 
Crespo-Medina M, Twing KI, Sánchez-Murillo R, Brazelton WJ, McCollom TM, Schrenk MO. Methane dynamics in a tropical serpentinizing environment: the Santa Elena Ophiolite, Costa Rica. Front Microbiol. 2017;8:916. https://doi.org/10.3389/fmicb.2017.00916.DeShon HR, Schwartz SY. Evidence for serpentinization of the forearc mantle wedge along the Nicoya Peninsula, Costa Rica. Geophys Res Lett. 2004;31. https://doi.org/10.1029/2004GL021179.Delmelle P, Stix J. Volcanic gases. In: Sigurdsson H, Houghton B, McNutt S, Rymer H, Stix J, editors. Encyclopedia of volcanoes. New York, NY: Elsevier; 2000. p 803–15.Kharaka YK, Mariner RH. Geothermal systems. In: Sigurdsson H, Houghton B, McNutt S, Rymer H, Stix J, editors. Encyclopedia of volcanoes. New York, NY: Elsevier; 2000. p. 817–34.Badger MR, Bek EJ. Multiple Rubisco forms in proteobacteria: their functional significance in relation to CO2 acquisition by the CBB cycle. J Exp Bot. 2008;59:1525–41.CAS 
PubMed 
Article 

Google Scholar 
West-Roberts JA, Carnevali PB, Scholmerich MC, Al-Shayeb B, Thomas A, Sharrar AM, et al. The Chloroflexi supergroup is metabolically diverse and representatives have novel genes for non-photosynthesis based CO2 fixation. bioRxiv [Preprint]. 2021. Available from: https://doi.org/10.1101/2020.05.14.094862.Lloyd KG, Steen AD, Ladau J, Yin J, Crosby L. Phylogenetically novel uncultured microbial cells dominate earth microbiomes. mSystems. 2018;3:1–12.Article 

Google Scholar 
Colman DR, Lindsay MR, Boyd ES. Mixing of meteoric and geothermal fluids supports hyperdiverse chemosynthetic hydrothermal communities. Nat Commun. 2019;10:1–3.Article 

Google Scholar  More