Philippa Kaur

Smith, S. E. & Read, D. J. Mycorrhizal Symbiosis (Academic Press, 2008).Tedersoo, L., Bahram, M. & Zobel, M. How mycorrhizal associations drive plant population and community biology. Science 367, eaba1223 (2020).CAS 
Article 

Google Scholar 
Connell, J. H. & Lowman, M. D. Low-diversity tropical rain forests: some possible mechanisms for their existence. Am. Nat. 134, 88–119 (1989).Article 

Google Scholar 
Brundrett, M. in Advances in Ecological Research, Vol. 21 (eds Begon M. et al.) 171–313 (Academic Press, 1991).Allen, E. B. et al. Patterns and regulation of mycorrhizal plant and fungal diversity. Plant Soil 170, 47–62 (1995).CAS 
Article 

Google Scholar 
Laliberté, E., Lambers, H., Burgess, T. I. & Wright, S. J. Phosphorus limitation, soil-borne pathogens and the coexistence of plant species in hyperdiverse forests and shrublands. New Phytol. 206, 507–521 (2015).Article 

Google Scholar 
Bennett, J. A. et al. Plant-soil feedbacks and mycorrhizal type influence temperate forest population dynamics. Science 355, 181–184 (2017).CAS 
Article 

Google Scholar 
Teste, F. P. et al. Plant-soil feedback and the maintenance of diversity in Mediterranean-climate shrublands. Science 355, 173–176 (2017).CAS 
Article 

Google Scholar 
van der Heijden, M. G. A., Martin, F. M., Selosse, M.-A. & Sanders, I. R. Mycorrhizal ecology and evolution: the past, the present, and the future. New Phytol. 205, 1406–1423 (2015).CAS 
Article 

Google Scholar 
Burrill, E. A. et al. The Forest Inventory and Analysis Database: Database Description and User Guide Version 8.0 for Phase 2 (U.S. Department of Agriculture, Forest Service, 2018).Jo, I., Fei, S., Oswalt, C. M., Domke, G. M. & Phillips, R. P. Shifts in dominant tree mycorrhizal associations in response to anthropogenic impacts. Sci. Adv. 5, eaav6358 (2019).CAS 
Article 

Google Scholar 
Jo, I. & Fei, S. Responses of Dominant Tree-Mycorrhizal Associations to Anthropogenic Impacts in the USA (Purdue Univ. Research Repository, 2019).Read, D. J. Mycorrhizas in ecosystems. Experientia 47, 376–391 (1991).Article 

Google Scholar 
Laliberté, E., Zemunik, G. & Turner, B. L. Environmental filtering explains variation in plant diversity along resource gradients. Science 345, 1602–1605 (2014).Article 

Google Scholar 
Ricklefs, R. E. Community diversity: relative roles of local and regional processes. Science 235, 167–171 (1987).CAS 
Article 

Google Scholar 
Cleland, D. T. et al. Ecological Subregions: Sections and Subsections for the Conterminous United States. General Technical Report WO-76D (U.S. Department of Agriculture, 2007).Tedersoo, L. & Bahram, M. Mycorrhizal types differ in ecophysiology and alter plant nutrition and soil processes. Biol. Rev. Camb. Philos. Soc. 94, 1857–1880 (2019).Article 

Google Scholar 
Taylor, D. L. et al. A first comprehensive census of fungi in soil reveals both hyperdiversity and fine-scale niche partitioning. Ecol. Monogr. 84, 3–20 (2014).Article 

Google Scholar 
Mariotte, P. et al. Plant–soil feedback: bridging natural and agricultural sciences. Trends Ecol. Evol. 33, 129–142 (2018).Article 

Google Scholar 
Zemunik, G., Turner, B. L., Lambers, H. & Laliberté, E. Diversity of plant nutrient-acquisition strategies increases during long-term ecosystem development. Nat. Plants 1, 15050 (2015).CAS 
Article 

Google Scholar 
van der Heijden, M. G. A. et al. Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity. Nature 396, 69–72 (1998).CAS 
Article 

Google Scholar 
Maherali, H. & Klironomos, J. N. Influence of phylogeny on fungal community assembly and ecosystem functioning. Science 316, 1746–1748 (2007).CAS 
Article 

Google Scholar 
Lin, G., McCormack, M. L. & Guo, D. Arbuscular mycorrhizal fungal effects on plant competition and community structure. J. Ecol. 103, 1224–1232 (2015).CAS 
Article 

Google Scholar 
Hartnett, D. C. & Wilson, G. W. T. The role of mycorrhizas in plant community structure and dynamics: lessons from grasslands. Plant Soil 244, 319–331 (2002).CAS 
Article 

Google Scholar 
Mangan, S. A. et al. Negative plant–soil feedback predicts tree-species relative abundance in a tropical forest. Nature 466, 752–755 (2010).CAS 
Article 

Google Scholar 
Peh, K. S.-H., Lewis, S. L. & Lloyd, J. Mechanisms of monodominance in diverse tropical tree-dominated systems. J. Ecol. 99, 891–898 (2011).Article 

Google Scholar 
ter Steege, H. et al. Hyperdominance in the Amazonian tree flora. Science 342, 1243092 (2013).Article 

Google Scholar 
Newman, E. I. & Reddell, P. Relationship between mycorrhizal infection and diversity in vegetation: evidence from the Great Smoky Mountains. Funct. Ecol. 2, 259–262 (1988).Article 

Google Scholar 
Bahram, M. et al. Plant nutrient-acquisition strategies drive topsoil microbiome structure and function. New Phytol. 227, 1189–1199 (2020).CAS 
Article 

Google Scholar 
Beisner, B. E., Haydon, D. T. & Cuddington, K. Alternative stable states in ecology. Front. Ecol. Environ. 1, 376–382 (2003).Article 

Google Scholar 
Zhu, K., Woodall, C. W. & Clark, J. S. Failure to migrate: lack of tree range expansion in response to climate change. Glob. Change Biol. 18, 1042–1052 (2012).Article 

Google Scholar 
Phillips, R. P., Brzostek, E. & Midgley, M. G. The mycorrhizal-associated nutrient economy: a new framework for predicting carbon–nutrient couplings in temperate forests. New Phytol. 199, 41–51 (2013).CAS 
Article 

Google Scholar 
Jenkins, C. N., Van Houtan, K. S., Pimm, S. L. & Sexton, J. O. US protected lands mismatch biodiversity priorities. Proc. Natl Acad. Sci. USA 112, 5081–5086 (2015).CAS 
Article 

Google Scholar 
Chao, A., Chiu, C.-H. & Jost, L. Unifying species diversity, phylogenetic diversity, functional diversity, and related similarity and differentiation measures through Hill numbers. Annu. Rev. Ecol. Evol. Syst. 45, 297–324 (2014).Article 

Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2018).Bürkner, P.-C. brms: an R package for Bayesian multilevel models using Stan. J. Stat. Softw. 80, 1–28 (2017).Article 

Google Scholar 
Dowle, M. & Srinivasan, A. data.table: Extension of ‘data.frame’ R package version 1.14.0 https://CRAN.R-project.org/package=data.table (2017).Wickham, H., Francois, R., Henry, L. & Müller, K. dplyr: A grammar of data manipulation. R package version 1.0.2 https://CRAN.R-project.org/package=dplyr (2017).Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer-Verlag, 2016).Book 

Google Scholar 
Kassambara, A. ggpubr: ‘ggplot2’ Based Publication Ready Plots. R package version 3.3.3 https://CRAN.R-project.org/package=ggpubr (2018).Dunnington, D. ggspatial: Spatial Data Framework for ggplot2. R package version 1.1.1 https://CRAN.R-project.org/package=ggspatial (2018).Hijmans, R. J. raster: Geographic Data Analysis and Modeling. R package version 3.4-5 https://CRAN.R-project.org/package=raster (2019).Wickham, H. Reshaping data with the reshape Package. J. Stat. Softw. 21, 1–20 (2007).Article 

Google Scholar 
Pebesma, E. Simple features for R: standardized support for spatial vector data. R. J. 10, 439–446 (2018).Article 

Google Scholar 
Wickham, H. & Henry, L. tidyr: Tidy messy data. R package version 1.1.0 https://CRAN.R-project.org/package=tidyr (2019).Oksanen, J. et al. vegan: Community ecology package. R package version 2.5-6 https://cran.r-project.org/web/packages/vegan/index.html (2017).Forest Inventory and Analysis; https://apps.fs.usda.gov/fia/datamart/CSV/ENTIRE.zipCarteron, A. alexiscarter/mycorrhiza_tree_diversity: custom code https://doi.org/10.5281/ZENODO.5713273 (2021). More

Anderegg, W. R. L. AGU Advances 2, e2021AV000490 (2021).Article 

Google Scholar 
Goldstein, A. et al. Nat. Clim. Chang. 10, 287–295 (2020).CAS 
Article 

Google Scholar 
Tanneberger, F. et al. Adv. Sustain. Syst. 5, 2000146 (2021).CAS 
Article 

Google Scholar 
Paustian, K. et al. Nature 532, 49–57 (2016).CAS 
Article 

Google Scholar 
Magill, B. Biden climate plan to save forests pivots on swamps, wetland (1). Bloomberg Law, https://go.nature.com/3gzaYKM (2 November 2021).Meier, J. Designating Dark Sky Areas: Actors and Interests (Routledge, 2014).Climate Action Reserve. Key Accounting Principles for Improved Forest Management Projects within the Forest Protocol (Climate Action Reserve, 2019).Verified Carbon Standard. VM0007 REDD+ Methodology Framework (REDD+ MF), v1.6. verra.org, https://go.nature.com/3GCEKbZ (8 September 2020).PCAF. The Global GHG Accounting and Reporting Standard for the Financial Industry (Partnership for Carbon Accounting Financials, 2021).Xu, J., Morris, P. J., Liu, J. & Holden, J. Catena 160, 134–140 (2018).Article 

Google Scholar 
Yu, Z., Loisel, J., Brosseau, D. P., Beilman, D. W. & Hunt, S. J. Geophys. Res. Lett. 37, L13402 (2010).
Google Scholar 
Loisel, J. et al. Nat. Clim. Chang. 11, 70–77 (2021).Article 

Google Scholar 
Leifeld, J., Wüst-Galley, C. & Page, S. Nat. Clim. Chang. 9, 945–947 (2019).CAS 
Article 

Google Scholar 
Leifeld, J. & Menichetti, L. Nat. Commun. 9, 1071 (2018).CAS 
Article 

Google Scholar 
Juffe-Bignoli, D. et al. Protected Planet Report 2014 (UNEP-WCMC: 2014).Hoyos-Santillan, J., Miranda, A., Lara, A., Rojas, M. & Sepulveda-Jauregui, A. Science 366, 1207–1208 (2019).Article 

Google Scholar 
von Unger, M. & Emmer, I. Carbon Market Incentives to Conserve, Restore and Enhance Soil Carbon. (The Nature Conservancy, 2018).Bonn, A. et al. (eds). Peatland Restoration and Ecosystem Services: Science, Policy and Practice (Cambridge Univ. Press, 2016).Buotte, P. C., Law, B. E., Ripple, W. J. & Berner, L. T. Ecol. Appl. 30, e02039 (2020).Article 

Google Scholar 
Soto-Navarro, C. et al. Phil. Trans. R. Soc. Lond. B 375, 20190128 (2020).CAS 
Article 

Google Scholar  More

Chemical and pathogenic hazards in aquaculture supply chains threaten the provision of safe aquatic food. The Seafood Risk Tool is an integrated, semi-quantitative system that develops bespoke supply chain and risk management strategies.Although wild fish catches have plateaued globally, the aquaculture sector continues to expand in response to increasing demand for fish and other aquatic foods. Economic progress and the growing consumer awareness of aquatic foods in sustainable, healthy and nutritious diets are contributing to sector expansion1. However, rapid expansion must be achieved in a socially responsible and environmentally sustainable manner. Aquaculture produces over 400 different species across marine and freshwater fish, crustaceans, molluscs, plants and algae2 — all presenting complex and unique risk profiles to the environment, the industry, investors and consumers. Measures to mitigate these risks, including aquaculture certification and legislation, are inconsistent across nations and regions, and cohesive risk management has been difficult to manage and implement. More

Barcoding of Pseudomonas isolates and experimental designTo test possible host–commensal–pathogen dynamics in a local population, we spray inoculated six A. thaliana accessions with synthetic bacterial communities composed of pathogenic and commensal Pseudomonas candidates. Because we wanted to study interactions that are likely to occur in nature, we used A. thaliana genotypes that originated from the same plant populations near Tübingen, Germany28, from which the Pseudomonas strains had been isolated (Fig. 1a). Classification of Pseudomonas lineages as pathogenic or commensal was based on observed effects in axenic infections11. Only one lineage, previously named OTU5, which dominated local plant populations, was associated with pathogenicity, both based on negative impact on rosette weight and visible disease symptoms11. We henceforth call this lineage ATUE5 (isolates sampled from ‘Around TUEbingen, group 5’) and all other Pseudomonas lineages from the Karasov collection non-ATUE5. We interchangeably use the terms ‘pathogens’ or ‘ATUE5’, and ‘commensals’ or ‘non-ATUE5’.Fig. 1: Study system.a, Location of original a. thaliana and Pseudomonas sampling sites around Tübingen. b, Taxonomic representation of the 14 Pseudomonas isolates used and the prevalence of closely related strains (divergence 0.99 and P value |± 0.2| shown. Node colours indicate the bacterial isolate classification, ATUE5 or non-ATUE5. e, In planta abundance change of the seven ATUE5 isolates in non-ATUE5 inclusive treatments in comparison with PathoCom. Abundance mean difference was estimated with the model [log10(isolate load) ~ treatment × experiment + treatment + experiment + error] for each individual strain. Thus, the treatment coefficient was estimated per isolate. Dots indicate the median estimates, and vertical lines represent 95% Bayesian credible intervals of the fitted parameter. ‘Combi’ indicates combination of the isolates C3,C4,C5 and C7, and n = 23.Full size imageEach of the 14 isolates was examined for growth inhibition against all other isolates, covering all possible combinations of binary interactions. In total, three strains out of the 14 had inhibitory activity; all were non-ATUE5 (Fig. 4c). Specifically, C4 and C5 showed the same pattern: both inhibited all pathogenic isolates but P1, and both inhibited the same two commensals, C6 and weakly C3. C3 inhibited three ATUE5 isolates: P5, P6 and P7. In summary, the in vitro assay provides evidence that among the tested Pseudomonas isolates, direct inhibition was a trait unique to commensals, and susceptible bacteria were primarily pathogens. This supports the notion that ATUE5 and non-ATUE5 isolates employ divergent competition strategies, or that if they use the same mechanism, they differ in the effectiveness of such a mechanism.The in vitro results recapitulated the general trend of pathogen inhibition found among treatments in planta. Nevertheless, we observed major discrepancies between the two assays. First, P1 was not inhibited by any isolate in the host-free assay (Fig. 4c), though it was the most inhibited member in planta among the communities (Fig. 4b). Second, no commensal isolate was inhibited in planta among communities (Fig. 4b), while two commensals, C3 and C6, were inhibited in vitro (Fig. 4c). Both observations are compatible with an effect of the host on microbe–microbe interactions. To explore such effects, we analysed all pairwise microbe–microbe abundance correlations within MixedCom-infected hosts. When we used absolute abundances, all pairwise correlations were positive, also in CommenCom and PathoCom (Extended Data Fig. 8a), consistent with there being a positive correlation between absolute abundance of individual isolates and total abundance of the entire community (Supplementary Fig. 7), that is, no isolate was less abundant in highly colonized plants than in sparsely colonized plants. This indicates that there does not seem to be active killing of competitors in planta in the CommenCom, which is probably not surprising. With relative abundances, however, a clear pattern emerged with a cluster of commensals that were positively correlated, possibly reflecting mutual growth promotion, and several commensal strains being negatively correlated with both P6 and C7, possibly reflecting unidirectional growth inhibition (Fig. 4d). We did not observe the same correlations within CommenCom among commensals and within PathoCom among pathogens as we did for either subgroup in MixedCom, reflecting higher-order interactions. Thus, interactions among pathogens were constrained by the presence of commensals and vice versa (Extended Data Fig. 8b).The in planta patterns measured in complex communities did not fully recapitulate what we had observed in vitro with pairwise interactions. We therefore investigated individual commensal isolates for their ability to suppress pathogens in planta and also tested the entourage effect. We focused on the three commensals C3, C4 and C5, which had directly inhibited pathogens in vitro, and as a control C7, which had not shown any inhibition activity in vitro. We infected plants with mixtures of PathoCom and each of the four individual commensals and also with PathoCom mixed with all four commensals. Because pathogen inhibition seemed to be independent of the host genotype, we arbitrarily chose HE-1. Regardless of the commensal isolate, only P1 was suppressed with high probability in all commensal-including treatments (Fig. 4e), with P2, P3 and P4 being substantially inhibited only by the mixture of all four commensals. Together with the lack of meaningful differences between individual commensals, this indicates that pathogen inhibition is either a function of commensal dose or a result of interaction among commensals.An important finding was that four commensal strains had much more similar inhibitory activity in planta than in vitro and that the combined action was greater than the individual effects. Together, this suggested that the host contributes to the observed interactions between commensal and pathogenic Pseudomonas isolates. To begin to investigate this possibility, we next studied potential host immune responses with RNA sequencing.Defensive response elicited by non-ATUE5For the RNA-sequencing experiment, we treated plants of the genotype Lu3-30 with the three synthetic communities and also used a bacteria-free control treatment. We sampled the treated plants 3 DPI and 4 DPI, thus increasing the ability to pinpoint differentially expressed genes (DEGs) between treatments that are not highly time specific. Exploratory analysis indicated that the two time points behaved similarly, and they were combined for further in-depth analysis.We first looked at DEGs in a comparison between infected plants and control (Supplementary Table 5); with PathoCom, there were only 14 DEGs; with CommenCom, there were 1,112 DEGs; and with MixedCom, there were 1,949 DEGs, suggesting that the CommenCom isolates, which are also present in the MixedCom, elicited a stronger host response than the PathoCom members. Furthermore, the high number of DEGs in MixedCom, higher than both PathoCom and CommenCom together, suggested a synergistic response derived from inclusion of both PathoCom and CommenCom members. Alternatively, this could also be a consequence of the higher initial inoculum in the 14-member MixedCom than either the 7-member PathoCom or 7-member CommenCom, or a combination of the two effects (Fig. 5a,b and Extended Data Fig. 9). The genes induced by the MixedCom fell into two classes: Group 5 (Fig. 5a,b) was also induced, albeit more weakly, by the CommenCom but not by the PathoCom. This group was overrepresented for non-redundant gene ontology (GO) categories linked to defence (Fig. 5c) and most likely explains the protective effects of commensals in the MixedCom. Specifically, among the top ten enriched GO categories in the shared MixedCom and CommenCom set, eight relate to immune response or response to another organism (‘defence response’, ‘multi-organism process’, ‘immune response’, ‘response to stimulus’, ‘response to biotic stimulus’, ‘response to other organism’, ‘immune system process’, ‘response to stress’; Fig. 5c).Fig. 5: Only commensal members elicit a strong host-defensive response.a, Relative expression (RE) pattern of 2,727 DEGs found in at least one of the comparisons of CommenCom, PathoCom and MixedCom with Control. DEGs were hierarchically clustered. b, Euler diagram of DEGs in PathoCom-, CommenCom- and MixedCom-treated plants compared with Control (log2[fold change]  >|± 1|; false discovery rate (FDR) 0.05); n = 4.Full size imageGroup 4 was only induced in MixedCom, either indicating synergism between commensals and pathogens or reflecting a consequence of the higher initial inoculum. This group included a small number of redundant GO categories indicative of defence, such as ‘salicylic acid mediated signalling pathway’, ‘multi-organism process’, ‘response to other organism’ and ‘response to biotic stimulus’ (Supplementary Table 6). Moreover, the MixedCom response cannot simply be explained by synergistic effects or commensals suppressing pathogen effects because there was a prominent class, Group 2, which included genes that were induced in the CommenCom but to a much lesser extent in the PathoCom or MixedCom. From their annotation, it was unclear how they can be linked to infection (Fig. 5c). About 500 genes (Group 1) that were downregulated by all bacterial communities are unlikely to contain candidates for commensal protection (Fig. 5a).Cumulatively, these results imply that the CommenCom members elicited a defensive response in the host regardless of PathoCom members, while the mixture of both led to additional responses. To better understand if selective suppression of ATUE5 in MixedCom infections may have resulted from the recognition of both non-ATUE5 and ATUE5 (reflected by a unique MixedCom set of DEGs) or solely non-ATUE5 (a set of DEGs shared by MixedCom and CommenCom), we examined the expression of key genes related to the salicylic acid pathway and downstream immune responses. Activation of the salicylic acid pathway was previously related to increased fitness of A. thaliana in the presence of wild bacterial pathogens, a phenomenon which was attributed to an increased systemic acquired resistance32.We observed a general trend of higher expression in MixedCom- and CommenCom-infected hosts for several such genes (Fig. 5d). Examples are PR1 and PR5, marker genes for systemic acquired resistance and resistance execution. Therefore, according to the marker genes we tested, non-ATUE5 elicited a defensive response in the host, regardless of ATUE5 presence.We conclude that the expression profiles of non-ATUE5-infected Lu3-30 plants point to an increased defensive status, supporting our hypothesis regarding host-mediated ATUE5 suppression. We note that ATUE5 suppression was not associated with full plant protection and thus control-like weight levels in all plant genotypes. One accession, Ey15-2, was only partially protected in the MixedCom (Fig. 2), despite levels of pathogen inhibition being not very different from other host genotypes (Extended Data Fig. 7).Lack of protection explained by a single pathogenic isolateThe fact that Ey15-2 was only partially protected by MixedCom (Fig. 2) underlines the importance of the host genotype in plant–microbe–microbe interactions, apparently reflecting the dynamics between microbes and plants in wild populations. We therefore wanted to reveal the cause for this differential interaction.Our first aim was to rank compositional variables in MixedCom according to their impact on plant weight, regardless of host genotype. Next, we asked whether any of the top-ranked variables could explain the lack of protection in Ey15-2. With Random Forest analysis, we estimated the weight-predictive power of all individual isolates in MixedCom and three cumulative variables: total bacterial abundance, total ATUE5 abundance and total non-ATUE5 abundance. We found that the best weight-predictive variable was the abundance of pathogenic isolate P6, followed by total bacterial load and total ATUE5 load, which were probably confounded by the abundance of P6 (Fig. 6a). In agreement, P6 was the dominant ATUE5 in MixedCom (Fig. 6b and Extended Data Fig. 10a). We thus hypothesized that the residual pathogenicity in MixedCom-infected Ey15-2 was caused by P6. Although P6 grew best in Ey15-2, the difference to most other genotypes was unlikely to be important (Extended Data Fig. 10b). However, P6 was particularly dominant in Ey15-2 (Fig. 6b).Fig. 6: The effect of isolate P6 on weight in MixedCom-infected hosts and particularly on accession Ey15-2.a, Relative importance (mean decrease accuracy; ‘MSE’) of 20 examined variables in weight prediction of MixedCom-infected hosts as determined by Random Forest analysis. The best predictor was the abundance of isolate P6. ‘Total bacterial’, ‘Total ATUE5’ and ‘Total non-ATUE5’ indicate the cumulative abundances of the 14 isolates, seven ATUE5 isolates and seven non-ATUE5 isolates, respectively. b, Abundance of P6 compared with the other 13 barcoded isolates in MixedCom-infected hosts across the six A. thaliana genotypes used in this study. Dots indicate the median estimates, and vertical lines represent 95% Bayesian credible intervals of the fitted parameter, following the model [log10(isolate load) ~ isolate × experiment + isolate + experiment + error]. Each genotype was analysed individually, thus the model was utilized for each genotype separately. The shaded area denotes the 95% Bayesian credible intervals for the isolate P6. c, Fresh rosette weight of Ey15-2 plants treated with Control, MixedCom and MixedCom without P6 (MixedCom ΔP6). Fresh rosette weight was measured 12 DPI. The top panel presents the raw data, with the breaks in the vertical black lines denoting the mean value of each group, and the vertical lines indicating standard deviation. The lower panel presents the mean difference to control, plotted as bootstrap sampling55,56, indicating the distribution of effect size that is compatible with the data. The 95% confidence intervals are indicated by the black vertical bars, and n = 19.Full size imageGiven that pathogen load in Ey15-2 was driven to a substantial extent by P6, we assumed that this isolate had a stronger impact on the weight of Ey15-2 than on other accessions. We experimentally validated that removal of P6 restored protection when Ey15-2 was infected with the MixedCom (Fig. 6c). To confirm that restored protection was due to the interaction of commensals with the five other pathogenic isolates (P1–P5), rather than simply removal of P6, we also treated Ey15-2 with PathoCom only, but not P6. The removal of P6 did not diminish the negative weight impact of PathoCom (P1–P5, Supplementary Fig. 8), implying that it was indeed the interaction between commensals with five out of six pathogenic isolates that mitigated the harmful effect of pathogens in Ey15-2 plants. More

Positive correlation between biosynthetic gene cluster (BGC) and phylogenetic distance in the genus Bacillus
BGCs are responsible for the synthesis of secondary metabolites involved in microbial interference competition. To investigate the relationship between BGC and phylogenetic distance within the genus Bacillus, we collected 4268 available Bacillus genomes covering 139 species from the NCBI database (Supplementary Data 1). Phylogenetic analysis based on the sequences of 120 ubiquitous single-copy proteins27 showed that the 139 species could be generally clustered into four clades (Fig. 1 and Supplementary Data 2; the phylogenetic tree including all the detailed species information is shown in Supplementary Fig. 1), including a subtilis clade that includes species from diverse niches and can be further divided into the subtilis and pumilus subclades, a cereus clade that contains typical pathogenic species (B. cereus, B. anthracis, B. thuringiensis, etc.), a megaterium clade, and a circulans clade.Fig. 1: Phylogram of the tested Bacillus genomes.The maximum likelihood (ML) phylogram of 4268 Bacillus genomes was based on the sequences of 120 ubiquitous single-copy proteins27. The phylogenetic tree shows that Bacillus species can be generally clustered into the subtilis (light green circle; further includes subtilis (dark green) and pumilus (blue) subclades as shown in the branches), cereus (red), megaterium (yellow), and circulans (gray) clades. For detailed information of the species, please refer to the phylogenetic tree in Supplementary Fig. 1.Full size imagePrediction using the bioinformatic tool antiSMASH15 detected 49,671 putative BGCs in the 4268 genomes, corresponding to an average of 11.6 BGCs per genome (Supplementary Data 3). The subtilis clade had the most BGCs, 13.1 BGCs per genome (Fig. 2a); the subtilis subclade especially accommodates a high abundance of BGCs as 13.6 per genome (Supplementary Fig. 2a), which corresponds to their adaptation in diverse competitive habitats such as plant rhizosphere. The cereus and megaterium clades possessed moderate number of BGCs as 11.7 and 7.4 per genome, respectively; while the circulans clade only had 4.3 BGCs/genome (Fig. 2a and Supplementary Table 1), suggesting a distinct physiological feature and niche adaptation strategy. The two most abundant BGC classes were nonribosomal peptide-synthetase (NRPS) and RiPPs, which had an abundance of 3.7 and 3.1 per genome on average, respectively (Supplementary Fig. 2b and Supplementary Table 1). Interestingly, subtilis clade accommodated significantly higher abundance of BGCs in another polyketide synthase (PKSother; 2.0 per genome vs. 0.0–1.1 per genome) and PKS-NRPS Hybrids (0.7 vs. 0.0–0.2) classes, as compared with the three other clades (Supplementary Table 1); while cereus clade had more BGCs in RiPPs than other clades on average (Supplementary Table 1). Overall, the profile of BGC products and classification was generally consistent with the phylogenetic tree (Supplementary Fig. 3).Fig. 2: Biosynthetic gene cluster (BGC) distribution is correlated with phylogeny in the genus Bacillus.a The numbers of BGCs in the 4268 Bacillus genomes from different clades as defined by antiSMASH15. In the violin plot, the centre line represents the median, violin edges show the 25th and 75th percentiles, and whiskers extend to 1.5× the interquartile range. b Hierarchal clustering among the 545 representative Bacillus genomes based on the abundance of the different biosynthesis gene cluster families (GCFs). Each column represents a GCF, which was classified through BiG-SCAPE by calculating the Jaccard index (JI), adjacency index (AI), and domain sequence similarity (DSS) of each BGC28; the color bar on the top of the heatmap represents the BGC class of each GCF, where PKS includes classes of PKSother and PKSI, PKS-NRPS means PKS-NRPS Hybrids, Others includes classes of saccharides, terpene, and others. Each row represents a Bacillus genome, and the abundance of each GCF in different genomes is shown in the heatmap. The left tree was constructed based on the distribution pattern of GCFs, which showes a similar pattern to the phylogram in Fig. 1. c The correlation between the BGC and phylogenetic distance of the 545 representative Bacillus genomes (P  More

Meredith, M. et al. Polar regions. IPCC Intergov. Panel Clim. Chang. Geneva, Switz. 3, 203–320 (2019).Raftery, A. E., Zimmer, A., Frierson, D. M. W., Startz, R. & Liu, P. Less than 2 °C warming by 2100 unlikely. Nat. Clim. Chang. 7, 637–641 (2017).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
Ernakovich, J. G. et al. Predicted responses of arctic and alpine ecosystems to altered seasonality under climate change. Glob. Chang. Biol. 20, 3256–3269 (2014).PubMed 
ADS 

Google Scholar 
Parmesan, C. & Yohe, G. A globally coherent fingerprint of climate change impacts across natural systems. Nature 421, 37–42 (2003).CAS 
PubMed 
ADS 

Google Scholar 
Myers-Smith, I. H. & Hik, D. S. Climate warming as a driver of tundra shrubline advance. J. Ecol. 106, 547–560 (2018).
Google Scholar 
Martin, A. C., Jeffers, E. S., Petrokofsky, G., Myers-Smith, I. & Macias-Fauria, M. Shrub growth and expansion in the Arctic tundra: An assessment of controlling factors using an evidence-based approach. Environ. Res. Lett. 12, (2017).Myers-Smith, I. H. et al. Climate sensitivity of shrub growth across the tundra biome. Nat. Clim. Chang. 5, 887–891 (2015).ADS 

Google Scholar 
Myers-Smith, I. H. et al. Complexity revealed in the greening of the Arctic. Nat. Clim. Chang. 10, 106–117 (2020).ADS 

Google Scholar 
Epstein, H. E., Myers-Smith, I. & Walker, D. A. Recent dynamics of arctic and sub-arctic vegetation. Environ. Res. Lett. 8, 015040 (2013).ADS 

Google Scholar 
Ackerman, D., Griffin, D., Hobbie, S. E. & Finlay, J. C. Arctic shrub growth trajectories differ across soil moisture levels. Glob. Chang. Biol. 23, 4294–4302 (2017).PubMed 

Google Scholar 
Carrer, M., Pellizzari, E., Prendin, A. L., Pividori, M. & Brunetti, M. Winter precipitation – not summer temperature – is still the main driver for Alpine shrub growth. Sci. Total Environ. 682, 171–179 (2019).CAS 
PubMed 
ADS 

Google Scholar 
Xu, Y., Ramanathan, V. & Washington, W. M. Observed high-altitude warming and snow cover retreat over Tibet and the Himalayas enhanced by black carbon aerosols. Atmos. Chem. Phys. 16, 1303–1315 (2016).CAS 
ADS 

Google Scholar 
Francon, L. et al. Assessing the effects of earlier snow melt-out on alpine shrub growth: The sooner the better? Ecol. Indic. 115, (2020).López-Blanco, E. et al. Exchange of CO2 in Arctic tundra: impacts of meteorological variations and biological disturbance. Biogeosciences 14, 4467–4483 (2017).ADS 

Google Scholar 
Lund, M. et al. Larval outbreaks in West Greenland: Instant and subsequent effects on tundra ecosystem productivity and CO2 exchange. Ambio 46, 26–38 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
Prendin, A. L. et al. Immediate and carry-over effects of insect outbreaks on vegetation growth in West Greenland assessed from cells to satellite. J. Biogeogr. 47, 87–100 (2020).
Google Scholar 
Hobbie, S. E., Nadelhoffer, K. J. & Högberg, P. A synthesis: The role of nutrients as constraints on carbon balances in boreal and arctic regions. Plant Soil 242, 163–170 (2002).CAS 

Google Scholar 
Bret-Harte, M. S., Shaver, G. R. & Chapin, F. S. Primary and secondary stem growth in arctic shrubs: Implications for community response to environmental change. J. Ecol. 90, 251–267 (2002).
Google Scholar 
Sullivan, P. F., Ellison, S. B. Z., McNown, R. W., Brownlee, A. H. & Sveinbjörnsson, B. Evidence of soil nutrient availability as the proximate constraint on growth of treeline trees in northwest Alaska. Ecology 96, 716–727 (2015).PubMed 

Google Scholar 
Craine, J. M. et al. Global patterns of foliar nitrogen isotopes and their relationships with climate, mycorrhizal fungi, foliar nutrient concentrations, and nitrogen availability. New Phytol. 183, 980–992 (2009).CAS 
PubMed 

Google Scholar 
Shaver, G. R. & Chapin, F. S. Long-term responses to factorial, NPK fertilizer treatment by Alaskan wet and moist tundra sedge species. Ecography (Cop.) 18, 259–275 (1995).
Google Scholar 
Choudhary, S., Blaud, A., Osborn, A. M., Press, M. C. & Phoenix, G. K. Nitrogen accumulation and partitioning in a High Arctic tundra ecosystem from extreme atmospheric N deposition events. Sci. Total Environ. 554–555, 303–310 (2016).PubMed 
ADS 

Google Scholar 
Bergström, A. K. & Jansson, M. Atmospheric nitrogen deposition has caused nitrogen enrichment and eutrophication of lakes in the northern hemisphere. Glob. Chang. Biol. 12, 635–643 (2006).ADS 

Google Scholar 
Wild, B. et al. Plant-derived compounds stimulate the decomposition of organic matter in arctic permafrost soils. Sci. Rep. 6, 25607 (2016).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
Pedersen, E. P., Elberling, B. & Michelsen, A. Foraging deeply: Depth-specific plant nitrogen uptake in response to climate-induced N-release and permafrost thaw in the High Arctic. Glob. Chang. Biol. 26, 6523–6536 (2020).PubMed 
ADS 

Google Scholar 
Mack, M. C., Schuur, E. A. G. & Bret-harte, M. S. Ecosystem carbon storage in arctic tundra reduced by long-term nutrient fertilization. 431, 658–661 (2004).
Google Scholar 
Zamin, T. J. & Grogan, P. Birch shrub growth in the low Arctic: the relative importance of experimental warming, enhanced nutrient availability, snow depth and caribou exclusion. Environ. Res. Lett. 7, 034027 (2012).ADS 

Google Scholar 
DeMarco, J., MacK, M. C., Bret-Harte, M. S., Burton, M. & Shaver, G. R. Long-term experimental warming and nutrient additions increase productivity in tall deciduous shrub tundra. Ecosphere 5, 1–22 (2014).
Google Scholar 
Zamin, T. J., Bret-Harte, M. S. & Grogan, P. Evergreen shrubs dominate responses to experimental summer warming and fertilization in Canadian mesic low arctic tundra. J. Ecol. 102, 749–766 (2014).
Google Scholar 
Fenger-Nielsen, R. et al. Footprints from the past: The influence of past human activities on vegetation and soil across five archaeological sites in Greenland. Sci. Total Environ. 654, 895–905 (2019).CAS 
PubMed 
ADS 

Google Scholar 
Forbes, B. C., Ebersole, J. J. & Strandberg, B. Anthropogenic disturbance and patch dynamics in Circumpolar Arctic ecosystems. Conserv. Biol. 15, 954–969 (2001).
Google Scholar 
Andersen, E. A. S. et al. Nitrogen isotopes reveal high N retention in plants and soil of old Norse and Inuit deposits along a wet-dry arctic fjord transect in Greenland. Plant Soil 455, 241–255 (2020).CAS 

Google Scholar 
Normand, S. et al. Legacies of historical human activities in Arctic woody plant dynamics. Annu. Rev. Environ. Resour. 42, 541–567 (2017).
Google Scholar 
Walker, D. A. et al. The Circumpolar Arctic vegetation map. J. Veg. Sci. 16, 267–282 (2005).
Google Scholar 
Cappelen, J., Vinther, B. M., Kern-Hansen, C., Laursen, E. V. & Jørgensen, P. V. Greenland-DMI Historical Climate Data Collection 1784–2020 (Danish Meteorological Institute, 2021).
Google Scholar 
Hollesen, J. et al. Predicting the loss of organic archaeological deposits at a regional scale in Greenland. Sci. Rep. 9, 9097 (2019).PubMed 
PubMed Central 
ADS 

Google Scholar 
Fenger-Nielsen, R. et al. Arctic archaeological sites threatened by climate change: A regional multi-threat assessment of sites in south-west Greenland. Archaeometry 62, 1280–1297 (2020).CAS 

Google Scholar 
Fettweis, X. et al. Reconstructions of the 1900–2015 Greenland ice sheet surface mass balance using the regional climate MAR model. Cryosphere 11, 1015–1033 (2017).ADS 

Google Scholar 
Berner, L. T. et al. Summer warming explains widespread but not uniform greening in the Arctic tundra biome. Nat. Commun. 11, 1–12 (2020).ADS 

Google Scholar 
Assmann, J. J. et al. Local snow melt and temperature—but not regional sea ice—explain variation in spring phenology in coastal Arctic tundra. Glob. Chang. Biol. 25, 2258–2274 (2019).PubMed 
ADS 

Google Scholar 
Bhatt, U. S. et al. Climate drivers of Arctic tundra variability and change using an indicators framework. Environ. Res. Lett. 16, (2021).Hollesen, J., Matthiesen, H. & Elberling, B. The impact of Climate Change on an archaeological site in the Arctic. Archaeometry 59, 1175–1189 (2017).CAS 

Google Scholar 
Tolvanen, A. & Henry, G. H. R. Responses of carbon and nitrogen concentrations in high arctic plants to experimental warming. Can. J. Bot. 79, 711–718 (2001).CAS 

Google Scholar 
Oppen, J. et al. Annual air temperature variability and biotic interactions explain tundra shrub species abundance. J. Veg. Sci. 32, (2021).Hobbie, S. E. Temperature and plant species control over litter decomposition in Alaskan tundra. Ecol. Monogr. 66, 503–522 (1996).
Google Scholar 
Nadelhoffer, K. J., Giblin, A. E., Shaver, G. R. & Laundre, J. A. Effects of temperature and substrate quality on element mineralization in six Arctic soils. Ecology 72, 242–253 (1991).
Google Scholar 
Arens, S. J. T., Sullivan, P. F. & Welker, J. M. Nonlinear responses to nitrogen and strong interactions with nitrogen and phosphorus additions drastically alter the structure and function of a high Arctic ecosystem. J. Geophys. Res. Biogeosciences 113, 1–10 (2008).
Google Scholar 
Baddeley, J. A., Woodin, S. J. & Alexander, I. J. Effects of increased nitrogen and phosphorus availability on the photosynthesis and nutrient relations of three Arctic dwarf shrubs from Svalbard. Funct. Ecol. 8, 676 (1994).
Google Scholar 
Anadon-Rosell, A. et al. Xylem anatomical and growth responses of the dwarf shrub Vaccinium myrtillus to experimental CO2 enrichment and soil warming at treeline. Sci. Total Environ. 642, 1172–1183 (2018).CAS 
PubMed 
ADS 

Google Scholar 
Dawes, M. A. et al. Soil warming and CO2 enrichment induce biomass shifts in alpine tree line vegetation. Glob. Chang. Biol. 21, 2005–2021 (2015).PubMed 
ADS 

Google Scholar 
Walker, M. D. et al. Plant community responses to experimental warming across the tundra biome. Proc. Natl. Acad. Sci. 103, 1342–1346 (2006).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
Matthiesen, H., Fenger-Nielsen, R. F., Harmsen, H., Madsen, C. K. & Hollesen, J. The impact of vegetation on archaeological sites in the low arctic in light of climate change. Arctic 73, 141–152 (2020).
Google Scholar 
Dahl, M. B. et al. Warming, shading and a moth outbreak reduce tundra carbon sink strength dramatically by changing plant cover and soil microbial activity. Sci. Rep. 7, 1–13 (2017).CAS 

Google Scholar 
Westergaard-Nielsen, A., Karami, M., Hansen, B. U., Westermann, S. & Elberling, B. Contrasting temperature trends across the ice-free part of Greenland. Sci. Rep. 8, 1586 (2018).PubMed 
PubMed Central 
ADS 

Google Scholar 
Schweingruber, F. H., Börner, A. & Schulze, E.-D. Atlas of Stem Anatomy in Herbs, Shrubs and Trees. (Springer, Berlin, 2013). https://doi.org/10.1007/978-3-642-20435-7Pellizzari, E., Camarero, J. J., Gazol, A., Sangüesa-Barreda, G. & Carrer, M. Wood anatomy and carbon-isotope discrimination support long-term hydraulic deterioration as a major cause of drought-induced dieback. Glob. Chang. Biol. 22, 2125–2137 (2016).PubMed 
ADS 

Google Scholar 
Myers-Smith, I. H. et al. Methods for measuring arctic and alpine shrub growth: A review. Earth-Science Rev. 140, 1–13 (2015).ADS 

Google Scholar 
Stokes, M. A. & Smiley, T. L. Introduction to Tree-Ring Dating. (University of Chicago Press, 1968).Cook, E. R., Briffa, K., Shiyatov, S. & Mazepa, V. Methods of Dendrochronology: Applications in the Environmental Sciences. (Kluwer Academic Publisher, 1990).Gärtner, H. & Schweingruber, F. H. Microscopic preparation techniques for plant stem analysis. Kessel 95, 132–150 (2013).
Google Scholar 
von Arx, G., Crivellaro, A., Prendin, A. L., Čufar, K. & Carrer, M. Quantitative wood anatomy—practical guidelines. Front. Plant Sci. 7, 781 (2016).
Google Scholar 
Holmes, R. L. Computer-assisted quality control in tree- ring dating and measurement. Tree-ring Bulletin 43, 69–78 (1983).
Google Scholar 
Belokopytova, L. V, Babushkina, E. A., Zhirnova, D. F., Panyushkina, I. P. & Vaganov, E. A. Pine and larch tracheids capture seasonal variations of climatic signal at moisture-limited sites. Trees 33, 227–242 (2019).Büntgen, U. et al. Temperature-induced recruitment pulses of Arctic dwarf shrub communities. J. Ecol. 103, 489–501 (2015).
Google Scholar 
Fritts., H. C. Dendrochronology and Dendroclimatology. in Tree Rings and Climate 1–54 (1976). https://doi.org/10.1016/B978-0-12-268450-0.50006-9Briffa, K. & Jones, P. Basic chronology statistics and assessment. in Methods of Dendrochronology: Applications in the Environmental Sciences 137–152 (Kluwer Academic Publishers, 1990).Zuur, A. F., Ieno, E. N., Walker, N., Saveliev, A. A. & Smith, G. M. Mixed effects models and extensions in ecology with R. (Springer New York, 2009). https://doi.org/10.1007/978-0-387-87458-6Gazol, A. & Camarero, J. J. Mediterranean dwarf shrubs and coexisting trees present different radial-growth synchronies and responses to climate. Plant Ecol. 213, 1687–1698 (2012).
Google Scholar 
Crawley, M. J. Mixed-Effects Models. in R Book Second edition 681–714 (2007).Zar, J. H. Biostatistical analysis Fifth edition. USA Prentice Hall 4165 4159–4165, (1999).Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting Linear Mixed-Effects Models Using lme4. J. Stat. Softw. 67, arXiv:1406.5823 (2015).Pinheiro, J. C. & Bates, D. M. Linear Mixed-Effects Models: Basic Concepts and Examples. in Mixed-Effects Models in S and S-PLUS 3–56 (Springer-Verlag, 2000). https://doi.org/10.1007/0-387-22747-4_1Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. B. lmerTest Package: Tests in Linear Mixed Effects Models. J. Stat. Softw. 82, (2017).R Core Team (2021). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/. More

Daszak, P., Cunningham, A. A. & Hyatt, A. D. Emerging infectious diseases of wildlife: threats to biodiversity and human health. Science 287, 443–449 (2000).CAS 
PubMed 
ADS 

Google Scholar 
Pedersen, A. B., Jones, K. E., Nunn, C. L. & Altizer, S. Infectious diseases and extinction risk in wild mammals. Conserv. Biol. 21, 1269–1279 (2007).PubMed 
PubMed Central 

Google Scholar 
Smith, K. F., Sax, D. F. & Lafferty, K. D. Evidence for the role of infectious disease in species extinction and endangerment. Conserv. Biol. 20, 1349–1357 (2006).PubMed 

Google Scholar 
Fisher, M. C. et al. Emerging fungal threats to animal, plant and ecosystem health. Nature 484, 186–194 (2012).CAS 
PubMed 
ADS 

Google Scholar 
Paré, J. A. & Sigler, L. An overview of reptile fungal pathogens in the genera Nannizziopsis, Paranannizziopsis, and Ophidiomyces. J. Herpetol. Med. Surg. 26, 46–53 (2016).
Google Scholar 
Warnecke, L. et al. Inoculation of bats with European Geomyces destructans supports the novel pathogen hypothesis for the origin of white-nose syndrome. Proc. Natl. Acad. Sci. 109, 6999–7003 (2012).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
Skerratt, L. F. et al. Spread of chytridiomycosis has caused the rapid global decline and extinction of frogs. EcoHealth 4, 125 (2007).
Google Scholar 
Franklinos, L. H. V. et al. Emerging fungal pathogen Ophidiomyces ophiodiicola in wild European snakes. Sci. Rep. 7, 3844 (2017).PubMed 
PubMed Central 
ADS 

Google Scholar 
Lorch, J. M. et al. Snake fungal disease: an emerging threat to wild snakes. Philos. Trans. R. Soc. B Biol. Sci. 371, 20150457 (2016).
Google Scholar 
Bustos, M. L., Nicolás Sánchez, M., Peichoto, M. E. & Teibler, G. P. First report of fungal disease in a South American snake. Rev. Investig. Vet. Perú 29, 1036–1042 (2018).Sun, P.-L. et al. Infection with Nannizziopsis guarroi and Ophidiomyces ophiodiicola in reptiles in Taiwan. Transbound. Emerg. Dis. https://doi.org/10.1111/tbed.14049 (2021).Article 
PubMed 

Google Scholar 
Haynes, E. et al. First report of ophidiomycosis in a free-ranging California kingsnake (Lampropeltis californiae) in California, USA. J. Wildl. Dis. 57, 246–249 (2021).CAS 
PubMed 

Google Scholar 
Takami, Y. et al. First report of ophidiomycosis in Asia caused by Ophidiomyces ophiodiicola in captive snakes in Japan. J. Vet. Med. Sci. 83, 1234–1239 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Lorch, J. M. et al. Experimental infection of snakes with Ophidiomyces ophiodiicola causes pathological changes that typify snake fungal disease. MBio 6, e01534 (2015).CAS 
PubMed 
PubMed Central 

Google Scholar 
Koenig, J. E. et al. Succession of microbial consortia in the developing infant gut microbiome. Proc. Natl. Acad. Sci. 108, 4578–4585 (2011).CAS 
PubMed 
ADS 

Google Scholar 
Huttenhower, C. et al. Structure, function and diversity of the healthy human microbiome. Nature 486, 207–214 (2012).CAS 
ADS 

Google Scholar 
Grice, E. A. & Segre, J. A. The skin microbiome. Nat. Rev. Microbiol. 9, 244–253 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
Gall, C. A. et al. The bacterial microbiome of Dermacentor andersoni ticks influences pathogen susceptibility. ISME J. 10, 1846–1855 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Gould, A. L. et al. Microbiome interactions shape host fitness. Proc. Natl. Acad. Sci. 115, E11951–E11960 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
Hanning, I. & Diaz-Sanchez, S. The functionality of the gastrointestinal microbiome in non-human animals. Microbiome 3, 51 (2015).PubMed 
PubMed Central 

Google Scholar 
Jani, A. J. & Briggs, C. J. The pathogen Batrachochytrium dendrobatidis disturbs the frog skin microbiome during a natural epidemic and experimental infection. Proc. Natl. Acad. Sci. 111, E5049–E5058 (2014).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
Walker, D. M. et al. Variability in snake skin microbial assemblages across spatial scales and disease states. ISME J. 13, 2209–2222 (2019).PubMed 
PubMed Central 

Google Scholar 
Allender, M. C., Baker, S., Britton, M. & Kent, A. D. Snake fungal disease alters skin bacterial and fungal diversity in an endangered rattlesnake. Sci. Rep. 8, 12147 (2018).PubMed 
PubMed Central 
ADS 

Google Scholar 
Rykiel, E. J. Jr. Towards a definition of ecological disturbance. Aust. J. Ecol. 10, 361–365 (1985).
Google Scholar 
Kong, H. H. et al. Temporal shifts in the skin microbiome associated with disease flares and treatment in children with atopic dermatitis. Genome Res. 22, 850–859 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Ferrenberg, S. et al. Changes in assembly processes in soil bacterial communities following a wildfire disturbance. ISME J. 7, 1102–1111 (2013).PubMed 
PubMed Central 

Google Scholar 
Mackey, R. L. & Currie, D. J. The diversity–disturbance relationship: is it generally strong and peaked?. Ecology 82, 3479–3492 (2001).
Google Scholar 
Connell, J. H. Diversity in tropical rain forests and coral reefs: high diversity of trees and corals is maintained only in a nonequilibrium state. Science 199, 1302–1310 (1978).CAS 
PubMed 
ADS 

Google Scholar 
Guthrie, A. L., Knowles, S., Ballmann, A. E. & Lorch, J. M. Detection of snake fungal disease due to Ophidiomyces ophiodiicola in Virginia, USA. J. Wildl. Dis. 52, 143–149 (2016).PubMed 

Google Scholar 
Chandler, H. C. et al. Ophidiomycosis prevalence in Georgia’s eastern indigo snake (Drymarchon couperi) populations. PLoS ONE 14, e0218351 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Tetzlaff, S. J. et al. Snake fungal disease affects behavior of free-ranging massasauga rattlesnakes (Sistrurus catenatus). Herpetol. Conserv. Biol. 12, 624–634 (2017).
Google Scholar 
Lind, C. M., McCoy, C. M. & Farrell, T. M. Tracking outcomes of snake fungal disease in free-ranging pygmy rattlesnakes (Sistrurus miliarius). J. Wildl. Dis. 54, 352–356 (2018).PubMed 

Google Scholar 
Lind, C. M., Lorch, J. M., Moore, I. T., Vernasco, B. J. & Farrell, T. M. Seasonal sex steroids indicate reproductive costs associated with snake fungal disease. J. Zool. 307, 104–110 (2019).
Google Scholar 
McKenzie, J. M. et al. Field diagnostics and seasonality of Ophidiomyces ophiodiicola in wild snake populations. EcoHealth 16, 141–150 (2019).PubMed 

Google Scholar 
McCoy, C. M., Lind, C. M. & Farrell, T. M. Environmental and physiological correlates of the severity of clinical signs of snake fungal disease in a population of pigmy rattlesnakes, Sistrurus miliarius. Conserv. Physiol. 5, cow077 (2017).Hill, A. J. et al. Common cutaneous bacteria isolated from snakes inhibit growth of Ophidiomyces ophiodiicola. EcoHealth 15, 109–120 (2018).PubMed 

Google Scholar 
Baker, S. et al. Case definition and diagnostic testing for Snake Fungal Disease. Herpetol. Rev. 50, 279–285 (2019).
Google Scholar 
Chase, J. M., Kraft, N. J. B., Smith, K. G., Vellend, M. & Inouye, B. D. Using null models to disentangle variation in community dissimilarity from variation in α-diversity. Ecosphere 2, art24 (2011).Agugliaro, J., Lind, C. M., Lorch, J. M. & Farrell, T. M. An emerging fungal pathogen is associated with increased resting metabolic rate and total evaporative water loss rate in a winter-active snake. Funct. Ecol. 34, 486–496 (2020).
Google Scholar 
Frick, W. F. et al. Pathogen dynamics during invasion and establishment of white-nose syndrome explain mechanisms of host persistence. Ecology 98, 624–631 (2017).PubMed 

Google Scholar 
Gervasi, S. S., Hunt, E. G., Lowry, M. & Blaustein, A. R. Temporal patterns in immunity, infection load and disease susceptibility: understanding the drivers of host responses in the amphibian-chytrid fungus system. Funct. Ecol. 28, 569–578 (2014).
Google Scholar 
Allender, M. C. et al. Development of snake fungal disease after experimental challenge with Ophidiomyces ophiodiicola in cottonmouths (Agkistrodon piscivorous). PLoS ONE 10, e0140193 (2015).PubMed 
PubMed Central 

Google Scholar 
Briggs, C. J., Knapp, R. A. & Vredenburg, V. T. Enzootic and epizootic dynamics of the chytrid fungal pathogen of amphibians. Proc. Natl. Acad. Sci. 107, 9695–9700 (2010).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
Neuman-Lee, L. A. et al. Assessing multiple endpoints of atrazine ingestion on gravid Northern Watersnakes (Nerodia sipedon) and their offspring. Environ. Toxicol. 29, 1072–1082 (2014).CAS 
PubMed 
ADS 

Google Scholar 
Turnbaugh, P. J. et al. A core gut microbiome in obese and lean twins. Nature 457, 480–484 (2009).CAS 
PubMed 
ADS 

Google Scholar 
Kueneman, J. G. et al. Community richness of amphibian skin bacteria correlates with bioclimate at the global scale. Nat. Ecol. Evol. 3, 381–389 (2019).PubMed 

Google Scholar 
Manichanh, C. et al. Reduced diversity of faecal microbiota in Crohn’s disease revealed by a metagenomic approach. Gut 55, 205–211 (2006).CAS 
PubMed 
PubMed Central 

Google Scholar 
Murphy, G. E. P. & Romanuk, T. N. A meta-analysis of declines in local species richness from human disturbances. Ecol. Evol. 4, 91–103 (2014).PubMed 

Google Scholar 
Jani, A. J. et al. The amphibian microbiome exhibits poor resilience following pathogen-induced disturbance. ISME J. 15, 1628–1640 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Zaneveld, J. R., McMinds, R. & Vega Thurber, R. Stress and stability: applying the Anna Karenina principle to animal microbiomes. Nat. Microbiol. 2, 1–8 (2017).
Google Scholar 
Anderson, M. J. et al. Navigating the multiple meanings of β diversity: a roadmap for the practicing ecologist. Ecol. Lett. 14, 19–28 (2011).PubMed 
ADS 

Google Scholar 
Lankau, E. W., Hong, P.-Y. & Mackie, R. I. Ecological drift and local exposures drive enteric bacterial community differences within species of Galápagos iguanas. Mol. Ecol. 21, 1779–1788 (2012).PubMed 

Google Scholar 
Mebert, K. Good species despite massive hybridization: genetic research on the contact zone between the watersnakes Nerodia sipedon and N. fasciata in the Carolinas, USA. Mol. Ecol. 17, 1918–1929 (2008).CAS 
PubMed 

Google Scholar 
Bohuski, E., Lorch, J. M., Griffin, K. M. & Blehert, D. S. TaqMan real-time polymerase chain reaction for detection of Ophidiomyces ophiodiicola, the fungus associated with snake fungal disease. BMC Vet. Res. 11, 95 (2015).PubMed 
PubMed Central 

Google Scholar 
Wiens, J. A. Spatial scaling in ecology. Funct. Ecol. 3, 385–397 (1989).
Google Scholar 
Caporaso, J. G. et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. 108, 4516–4522 (2011).CAS 
PubMed 
ADS 

Google Scholar 
Fadrosh, D. W. et al. An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform. Microbiome 2, 6 (2014).PubMed 
PubMed Central 

Google Scholar 
Schloss, P. D. et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541 (2009).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K. & Schloss, P. D. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the Miseq Illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120 (2013).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucl. Acids Res. 41, D590–D596 (2013).CAS 
PubMed 

Google Scholar 
Yilmaz, P. et al. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucl. Acids Res. 42, D643–D648 (2014).CAS 
PubMed 

Google Scholar 
Weiss, S. et al. Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome 5, 27 (2017).PubMed 
PubMed Central 

Google Scholar 
R Core Team. R: A language and environment for statistical computing. (R Foundation for Statistical Computing, 2021).Bozdogan, H. Model selection and akaike’s information criterion (AIC): the general theory and its analytical extensions. Psychometrika 52, 345–370 (1987).MathSciNet 
MATH 

Google Scholar 
Fox, J. & Weisberg, S. An R Companion to Applied Regression (SAGE Publications, Thousand Oaks, 2011).
Google Scholar 
Zuur, A., Ieno, E. N., Walker, N., Saveliev, A. A. & Smith, G. M. Mixed Effects Models and Extensions in Ecology with R (Springer, 2009).MATH 

Google Scholar 
Heip, C. A new index measuring evenness. J. Mar. Biol. Assoc. UK 54, 555–557 (1974).
Google Scholar  More

Leibniz Institute DSMZ German Collection of Microorganisms and Cell Cultures, Braunschweig, GermanyAmber Hartman Scholz, Rodrigo Sara, Scarlett Sett, Andrew Lee Hufton & Jörg OvermannLeibniz Institute of Plant Genetics and Crop Plant Research (IPK), Seeland, GermanyJens FreitagNatural History Museum, London, UKChristopher H. C. LyalOne Planet Solutions, Montpellier, FranceRodrigo SaraUniversidad de los Andes, Bogotá, ColombiaMartha Lucia CepedaPlentzia Marine Station (PiE-UPV/EHU), European Marine Biological Resource Centre – Spain (EMBRC-Spain), Plentzia, SpainIbon CancioEthiopian Biotechnology Institute, Addis Ababa, EthiopiaYemisrach Abebaw & Kassahun TesfayeNational Academy of Agricultural Science and Global Plant Council, New Delhi, IndiaKailash BansalNational Council of Scientific Research and Technologies (NCSRT), Algiers, AlgeriaHalima BenbouzaMinistry of Agriculture, Livestock, Fisheries and Cooperatives, Nairobi, KenyaHamadi Iddi BogaInstitut Pasteur, Paris, FranceSylvain Brisse, Anne-Caroline Deletoille & Raquel Hurtado-OrtizSchool of Biosciences, Cardiff University, Cardiff, UKMichael W. BrufordWellcome Sanger Institute, Hinxton, UKHayley Clissold & David NicholsonEuropean Molecular Biology Laboratory European Bioinformatics Institute (EMBL-EBI), Hinxton, UKGuy CochraneGlobal Genome Initiative, Smithsonian National Museum of Natural History, Washington, DC, USAJonathan A. CoddingtonAlexander von Humboldt Biological Resources Research Institute, Bogota, ColombiaFelipe García-CardonaSouth African National Biodiversity Institute, Cape Town, South AfricaMichelle Hamer, Jessica da Silva & Krystal A. TolleyUniversity of Nairobi, Nairobi, KenyaDouglas W. MianoInstituto Tecnologico Vale (ITV), Belem, BrazilGuilherme OliveiraMinistry of Environment and Sustainable Development, Bogota, ColombiaCarlos Ospina BravoUniversity of Lethbridge, Lethbridge, CanadaFabian RohdenNatural History Museum of Denmark, Copenhagen, DenmarkOle SebergUniversity of Freiburg, Freiburg, GermanyGernot SegelbacherNational Centre for Cell Science, Pune, IndiaYogesh ShoucheMariano Galvez University, Guatemala City, GuatemalaAlejandra Sierra National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, USAIlene Karsch-MizrachiCentre for Ecological Genomics and Wildlife Conservation, University of Johannesburg, Johannesburg, South AfricaJessica da Silva & Krystal A. TolleyUniversity of the Philippines Los Banos, Laguna, PhilippinesDesiree M. HauteaFundação Oswaldo Cruz (FIOCRUZ), Rio de Janeiro, BrazilManuela da SilvaNational Institute of Genetics, Mishima, JapanMutsuaki SuzukiInstitute of Biotechnology, Addis Ababa University, Addis Ababa, EthiopiaKassahun TesfayeCentre for Tropical Livestock Genetics and Health (CTLGH) – International Livestock Research Institute (ILRI), Nairobi, KenyaChristian Keambou TiamboMurdoch University, Murdoch, AustraliaRajeev VarshneyCorporación CorpoGen, Bogotá, ColombiaMaría Mercedes ZambranoTechnical University of Braunschweig, Braunschweig, GermanyJörg OvermannConceptualization: A.H.S., J.F., C.H.C.L., R.S., M.L.C., I.C., S.S., Y.A., K.B., H.B., H.I.B., S.Y., M.W.B., H.C., G.C., J.A.C., A.D., F.G.C., M.H., R.H.O., D.W.M., G.O., C.O.B., F.B., O.S., G.S., Y.S., A.S., J.d.S., M.d.S., M.S., K.T., K.A.T., M.M.Z., and J.O. Visualization: J.O., I.C., S.S., R.S., C.H.C.L., G.C., and A.H.S. Funding acquisition: A.H.S., J.F., and J.O. Writing—original draft: A.H.S., R.S., M.L.C., C.H.C.L., I.C., and S.S. Writing—review & editing: A.H.S., J.F., C.H.C.L., R.S., M.L.C., I.B., S.S., A.L.H., D.N., M.d.S., S.B., M.M.Z., O.S., K.T., K.A.T., R.H.O., J.d.S., C.K.T., R.V., J.O., D.H., and I.K.M. More