Cofunctioning of bacterial exometabolites drives root microbiota establishment
Widespread Production of Specialized Exometabolites among Root-Associated Bacteria. To investigate the prevalence of interbacterial competition mediated by secreted metabolites (hereafter referred to as exometabolites), we used a m odified B urkholder plate-based a ssay (mBA; ref. and SI Appendix , Fig. S1 A ) and tested 39,204 binary interbacterial interactions, i.e., 198 producer versus 198 target isolates. Independent validation of 7,470 randomly selected interactions (i.e., 19%) revealed 95% reproducibility of interaction phenotypes ( SI Appendix , Table S1 ). We detected an inhibition halo in 1,011 interactions (i.e., 2.6% of tested pairwise interactions involving 66% of the isolates, SI Appendix , Fig. S1 B and Table S1 ), suggesting antibiosis due to specialized exometabolites. Antagonistic interactions were detected between all bacterial classes tested, indicating that the production of exometabolites is common in the At -RSphere culture collection of commensals ( ). Actinobacteria isolates were most sensitive to all other classes, especially to Gammaproteobacteria, which showed the highest aggregated frequency and average intensity of inhibitory activities ( ). Since we observed taxonomic signals at the class level, we calculated the average halo of inhibition size for all target and producer isolates across all bacterial classes tested (see sensitivity scores in and inhibition scores in , see also SI Appendix , Fig. S1 B ). This revealed that only a few bacteria – mainly belonging to Pseudomonadaceae (R9, R68, R71, R329, R401, R562, and R569) – exhibited broad inhibition of phylogenetically diverse isolates ( ). Extensive strain-specific variation in inhibition scores across closely related bacteria was observed, suggesting large standing genetic variation for the production of bacterial exometabolites among root-derived isolates ( ). Finally, we observed a 2.7× higher inhibitory activity for root-derived bacteria compared to those originating from soil ( P = 0.011; Kruskal–Wallis followed by Dunn’s post-hoc test and Benjamini–Hochberg correction, ), suggesting that the production of exometabolites might be advantageous for bacterial root colonization. Ralstonia spp . are core members of the root microbiota of healthy Arabidopsis plants in nature ( ). We hypothesized that the bacterial root microbiota of A. thaliana contributes to preventing disease in natural environments and tested the inhibitory activities of 167 of the aforementioned bacteria against pathogenic Ralstonia solanacearum GMI1000 (hereafter referred to as Rs ; ref. and ). A subset of root- and soil-derived bacteria inhibited the growth of Rs in mBA experiments (10.9% and 10.3%, respectively). These inhibitory activities were mainly manifested by bacteria from three genera: Pseudomonas , Streptomyces, and Bacillus ( ). Genomic Capacity for Specialized Metabolite Production Explains Pronounced Inhibitory Activity. We hypothesized that the ability to produce inhibitory halos is correlated with the genome-encoded potential for the biosynthesis of specialized metabolites. We predicted BGCs for all strains tested in mBA experiments ( SI Appendix , Table S2 ) and examined whether halo producer strains encode more BGCs than nonproducers. The total number of BGCs was significantly increased in the antagonistic isolates ( P = 0.0166). BGCs, which are thought to be involved in the biosynthesis of nonribosomal peptides (NRPs), aryl polyene, and redox cofactors, were significantly enriched in halo-producing isolates compared with nonproducers ( ), suggesting that these compounds may be important for interbacterial competition. We next investigated the diversity of bacterial metabolites that could explain the observed inhibitory phenotypes. Metabolites were extracted from individual strains (n = 198) grown on the same agar medium used in mBA experiments with two organic solvents with different polarities, ethyl acetate and methanol, to capture greater chemical diversity. Liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) analysis of the resulting 396 samples yielded ~200,000 mass spectra, which were analyzed using the Global Natural Products Social Molecular Networking (GNPS) workflow ( ). Network analysis of the resulting mass spectra revealed 247 families of molecules with similar fragmentation patterns, each representing at least two structurally related analogs ( SI Appendix , Fig. S2 ). While 2,220 nodes were shared between multiple classes, the remaining 1,094 nodes were found to be produced only by individual bacterial classes. Here, Gammaproteobacteria showed the greatest number of class-specific metabolites, indicating that their pronounced inhibitory activity ( ) can be explained by the production of an enormous diversity of specialized metabolites. To determine whether bacteria secrete an increased diversity of specialized metabolites when interacting with competitors ( ), we additionally analyzed the metabolomes of 16 inhibitory zones from 10 producer isolates (R63, R71, R68, R342, R401, R562, R569, R690, R920, and R1310) with broad inhibitory activity against three target isolates (R472D3, R480, and R553) using the same UPLC-MS/MS GNPS workflow. Spectra comparison with those of individually cultured strains revealed 298 additional ions detected exclusively in inhibition zones ( ). Dereplication of the network revealed that the biosynthesis of polyketides, such as two additional congeners of the nactin antibiotics and peptidic compounds, such as cyclo(Trp-Pro), and congeners of a unique lipopeptide family, were specifically triggered by interactions with other bacterial strains, pointing to a possible role of these molecules in interbacterial competition. DAPG Contributes to the Inhibitory Activity of P. brassicacearum R401. P. brassicacearum R401 showed the greatest number of inhibitory interactions (>17× greater than average) and the largest average halo size (>10× greater than average; ). Except for Bacilli and Flavobacteria, 50 isolates belonging to all other tested bacterial classes were sensitive to the inhibitory activity of R401, particularly Actinobacteria ( ), indicating the production of exometabolites with broad spectrum inhibitory activity. antiSMASH-based analysis of a resequenced, circular R401 genome predicted 16 BGCs ( SI Appendix , Fig. S3 and Table S2 ), one of which exactly matches the phl operon of Pseudomonas protegens Pf-5 , which has been shown to encode the enzymatic machinery for the production of the Pseudomonas -specific polyketide (DAPG; ) ( ). Since no other genome from our culture collection harbors the phl operon ( SI Appendix , Fig. S3 ), this points to a role for DAPG in mediating the inhibitory activity of R401. We generated a marker-free deletion mutant of the key biosynthetic gene phlD in R401 ( ). In mBA experiments, the R401 Δ phld mutant was significantly—yet only partly—impaired in its antagonistic activity ( ), retaining 71% of its inhibitory activity toward Rs . As DAPG production was completely lacking in the Δ phld mutant ( ), this suggests that DAPG alone is insufficient to explain the full inhibitory activity of R401 toward Rs. DAPG and Pyoverdine Act Additively to Inhibit Taxonomically Distinct Root Microbiota Members. We performed a forward genetic screen to reveal additional determinants mediating the residual inhibitory activity of R401 Δ phld . First, we generated a R401 mini- Tn5 transposon library with >6,000 insertion mutants. We adopted a fluorescence-based bacterial co-culture system in liquid medium to test all the R401 insertion mutants individually for their ability to suppress growth of a GFP-expressing Rs GMI1000 strain ( Rs GMI1600; ref. ). Of the 230 candidates identified in the primary screen, 38 R401 Tn5 mutants were robustly impaired in Rs suppression after two independent rounds of validation ( and SI Appendix , Table S3 ). One of these candidate R401 mutants was also significantly impaired (on average 26.6%) in its inhibitory activity against Rs in mBA experiments on solid agar, suggesting that this mutant is impaired in the production of an inhibitory exometabolite ( SI Appendix , Table S3 ). This R401 mutant carries the Tn5 transposon insertion within the gene of a putative acyltransferase ( pvdY ) involved in the biosynthesis of the siderophore pyoverdine in Pseudomonas aeruginosa ( ). We validated the contribution of R401 pvdY to pyoverdine biosynthesis and Rs inhibition by generating an independent pvdY deletion mutant, Δ pvdy , and a R401 deletion mutant of the gene encoding NRP synthetase pvdL, located downstream of pvdY (Δ pvdl , ). In mBA experiments with Rs GMI1000, the R401 Δ pvdy strain phenocopied the R401 tn5::pvdy mutant and the mutant phenotype was complemented by the expression of pvdY under its native promoter in the R401 Δ pvdy background (Δ pvdy::pvdY ). R401 Δ pvdl exhibited a slightly weaker impairment of halo production compared to R401 Δ pvdy . It is likely that R401 PvdY is involved in hydroxyornithine acetylation—a component of R401 pyoverdine—while PvdL is involved in the initial amino acid condensation, forming the basis for the peptidic backbone of pyoverdines. Acetyl hydroxyornithine could be involved in the biosynthesis of a third, unknown exometabolite by R401, explaining the difference between both pyoverdine mutants observed in mBA experiments. We also generated a R401 double deletion strain, Δ pvdy Δ pvdl , which is impaired in its inhibitory activity against pathogenic Rs to a similar extent as the R401 Δ pvdy single mutant ( ). Using mass spectrometry, we confirmed that all generated R401 pyoverdine mutants had lost their ability to produce pyoverdine (R401 pyoverdine has the sequence: Glu-Q-Lys-AcOHOrn-Ala-Gly-Ser-Ser-OHAsp-Thr, with Q being the fluorophore moiety; ). All tested Pseudomonas isolates in the At -RSphere contain pyoverdine biosynthetic genes; however, for R9, no characteristic pyoverdine fluorescence was detected ( SI Appendix , Fig. S3 ), potentially explaining its low inhibitory activity in mBAs ( ). Therefore, we hypothesized that pyoverdines contribute more broadly to explaining the unusually high inhibitory activity of the Pseudomonadaceae detected in the mBA experiments and isolated pyoverdine mutants from another Pseudomonas root commensal in our collection. We generated a mini- Tn5 mutant library of approximately 6,000 insertion mutants of pyoverdine-producing Pseudomonas fluorescens R569 and assayed about 2,000 mutants for loss of pyoverdine-mediated fluorescence by fluorimetry ( SI Appendix , Fig. S4 A ). Characterization of the Tn5 integration sites revealed one mutant in each of the R401 pvdY and pvdL homologs ( SI Appendix , Table S3 and Fig. S4 B and C ). Unlike the partial loss of inhibitory activity of the R401 Δ pvdy and Δ pvdl single and Δ pvdy Δ pvdl double mutants against Rs , the R569 Δ pvdy and Δ pvdl single mutants both showed a complete loss of inhibitory activity to Rs in mBA experiments ( SI Appendix , Fig. S4 D ). Although the genomes of root commensals R401 and R569 are assigned to different Pseudomonas sublineages ( ) and the arrangement of genes in the pvd operon is different ( and SI Appendix , Fig. S4 C ), it is likely that the corresponding pyoverdines are directly responsible for the observed inhibitory activity. Our findings indicate that pyoverdine might either be the sole or one of the several exometabolites produced by root commensals that limit the growth of Rs . Since pyoverdines function as iron chelators, we tested whether their inhibitory activity can be explained by interbacterial competition for iron by supplementing mBA agar medium with excess ferric iron (100 µM FeCl 3 ). This resulted in significantly impaired Rs inhibition when confronted with WT R401 or undetectable inhibition when confronted with WT R569, thus phenocopying results obtained with the corresponding strain-specific pyoverdine mutants ( and SI Appendix , Fig. S4 D ). Given that the mutants also exhibit drastically reduced iron mobilization capacity compared with the corresponding WT isolates ( and SI Appendix , Fig. S4 E ), we conclude that their inhibitory activities against root commensals and pathogenic Rs are likely mediated through their iron chelator function. However, we cannot exclude a possible additional pyoverdine activity under limited iron conditions. To test whether these two classes of compounds account for the full inhibitory activity of R401, we generated two R401 double pyoverdine and DAPG mutants, Δ pvdy Δ phld and Δ pvdl Δ phld , and conducted mBA assays using Rs and a representative set of the aforementioned commensals that we previously observed to be sensitive to R401 ( SI Appendix , Table S1 ). The double pyoverdine and DAPG mutants showed severely reduced inhibitory activity against Rs, and all other tested isolates compared to the single mutants, suggesting a cumulative impact of the two metabolites. Although the general sensitivity patterns were again isolate specific, Actinobacteria were typically inhibited by DAPG alone. DAPG and pyoverdine collectively explained >70% of R401 inhibitory activity, but a residual halo of inhibition was still observed, pointing to at least a third—yet to be defined—exometabolite ( ). Beyond the known inhibitory activities of DAPG and pyoverdine, our results provide evidence that metabolites with distinct modes of action jointly act to inhibit a taxonomically broad range of bacteria. DAPG and Pyoverdine Modulate Root Microbiota Assembly and Restrict Bacterial Diversity. We conducted root microbiota reconstitution experiments with germ-free A. thaliana Col-0 in Flowpots to study the contribution of DAPG and pyoverdine to root microbiota structure and to limiting the growth of Rs ( ). We tested heat-killed (HK) or live WT R401, R401 Δ pvdy, and Δ pvdl single as well as Δ pvdy Δ phld and Δ pvdl Δ phld double mutants on a phylogenetically diverse synthetic community (SynCom) that comprises 18 bacterial isolates ( SI Appendix , Fig. S5 A and Table S1 ), in both “soil” (peat matrix) and root compartments. After 3 weeks of plant—microbe cocultivation, DNA was extracted from soil and root samples and was subjected to 16S rRNA amplicon sequencing at isolate-specific resolution. Shoot fresh weight did not differ between conditions, and no wilting was observed ( SI Appendix , Fig. S5 B and C ), indicating that Rs was unable to cause disease on A. thaliana , possibly due to the presence of the SynCom. This is consistent with very low Rs relative abundances found in root samples ( SI Appendix , Fig. S5 D ) and a previous report that direct immersion of A. thaliana roots with a very high Rs inoculum of 10 8 cells/mL was needed to induce wilting symptoms in the crucifer ( ). The addition of live, WT R401 to the SynCom significantly reduced bacterial alpha diversity (Shannon index) compared to the (HK) R401 condition in the root compartment ( P < 0.001). This major impact on the community is gradually lost when R401 exometabolite mutants were co-inoculated with the SynCom ( ). Importantly, no such effects were observed in the soil compartment ( ). Bacterial beta diversity (Bray—Curtis dissimilarity) was also drastically affected by WT R401 inoculation, explaining most of the variation and resulting in a clear separation along axis 1 (HK versus WT, P < 0.001; R2 = 0.63). All mutant samples fall in between these extremes and follow a clear trajectory, WT > single mutants > double mutants > HK, suggesting that DAPG and pyoverdine have a cumulative influence on bacterial community structure in the root compartment ( ). This is likely due to direct exometabolite activities, as the production of either metabolite did not affect plant phenotypes ( SI Appendix , Fig. S5 C ). Furthermore, in silico depletion of R401 16S rRNA sequence reads leads to similar changes in beta diversity with even higher significance levels in some cases (Δ pvdl versus WT; P = 0.002; R2 = 0.14; SI Appendix , Fig. S5 E ), indicating that, irrespective of R401 abundance, the bacterial community is altered by DAPG and pyoverdine production. While in the soil compartment, inoculation of WT R401 leads to a similar shift as in the root compartment (HK versus WT, P < 0.001; R2 = 0.66), genetic depletion of either DAPG or pyoverdine biosynthetic abilities alone or together does not alter community structure ( ), which indicates niche-specific activity of these exometabolites in the root compartment. The impact of R401 on community structure was still significant upon in silico depletion of R401 16S rRNA sequences (HK versus WT , P < 0.001; R2 = 0.35; SI Appendix , Fig. S5 F ), indicating that R401 uses other, unknown mechanisms to influence community structure in soil. To determine whether the mBA results obtained in vitro have physiological relevance in planta , we inspected relative isolate abundances across conditions, reasoning that isolates insensitive to at least one of the R401 exometabolites in vitro would not benefit from genetic disruption of either BGC in a community context in the root compartment. This analysis revealed that only DAPG- and/or pyoverdine-sensitive isolates benefited from disruption of the corresponding BGCs present in R401 ( ), suggesting that binary interaction data in vitro can explain the impact on individual commensal isolates in a community context in the root compartment. This observation prompted us to test whether mBA data could predict community-scale effects as well as effects on single isolates. We computed the average reduction in inhibitory activity for each R401 mutant ( ) and tested whether the lack of competitiveness could inform changes in alpha- and beta-diversity. Linear regression analyses revealed that binary interaction data largely explained the observed effect size in the root compartment for both alpha- and beta-diversity indices but had no predictive power in the soil compartment ( , respectively). In conclusion, two exometabolites produced by a single isolate have large effects on key ecological indices of a taxonomically diverse SynCom and can be linked to pairwise in vitro interaction experiments. DAPG and Pyoverdine Act as Root Competence Determinants in a Community Context. To examine whether R401-induced modulation of bacterial assembly and diversity promoted R401 competitiveness, we determined the relative abundance of R401 in SynCom samples collected from root and soil compartments. In both root and soil samples, the 16S rRNA reads of live R401 by far exceeded the barely detectable HK R401 reads, indicating that live R401 inoculum proliferates in both compartments under all conditions ( ). However, R401 accumulated in association with roots at >2x higher abundance than in the soil compartment. R401 accumulation in the root compartment was gradually reduced in the SynCom context: Single or double mutations of DAPG and pyoverdine biosynthetic genes were sufficient to reduce the R401 abundance by up to approximately 39% compared to WT R401 ( ). Importantly, the capacity of these deletion mutants to colonize the soil compartment remained unaffected ( ). We also investigated the root colonization capacity of WT and all the five R401 deletion strains in mono-associations on axenic A. thaliana plants in an agar-based system ( ) and found no significant differences in live R401 cell counts ( ). Neither the growth of R401 DAPG or pyoverdine mutants was impaired in axenic culture media ( SI Appendix , Fig. S6 B and C ). Similarly, colonization experiments with WT R401 or the Δ pvdy Δ phld double mutant in the Flowpot system revealed no difference in root colonization capacity in the absence of bacterial competitors ( SI Appendix , Fig. S6 A ). Thus, DAPG and pyoverdine cofunction as R401 root competence determinants specifically in competition with other members of this SynCom in the root compartment. To assess the prevalence of Pseudomonas strains capable of producing DAPG and pyoverdine in association with plants when grown in natural soils, we analyzed the genomes of several commensal culture collections of 1,567 isolates from roots or leaves of A. thaliana or roots of the legume Lotus japonicus ( , , , ). These plants had been grown in the same or different soils on different continents, and our genome analysis revealed increased abundance of DAPG and pyoverdine BGCs in root-derived Pseudomonas isolates ( SI Appendix , Fig. S6 D and E ). These results not only support our conclusions obtained with a defined core commensal community and gnotobiotic A. thaliana but also suggest a broader yet root-specific role of these two exometabolites in natural environments beyond the model crucifer. 

To investigate the prevalence of interbacterial competition mediated by secreted metabolites (hereafter referred to as exometabolites), we used a m odified B urkholder plate-based a ssay (mBA; ref. and SI Appendix , Fig. S1 A ) and tested 39,204 binary interbacterial interactions, i.e., 198 producer versus 198 target isolates. Independent validation of 7,470 randomly selected interactions (i.e., 19%) revealed 95% reproducibility of interaction phenotypes ( SI Appendix , Table S1 ). We detected an inhibition halo in 1,011 interactions (i.e., 2.6% of tested pairwise interactions involving 66% of the isolates, SI Appendix , Fig. S1 B and Table S1 ), suggesting antibiosis due to specialized exometabolites. Antagonistic interactions were detected between all bacterial classes tested, indicating that the production of exometabolites is common in the At -RSphere culture collection of commensals ( ). Actinobacteria isolates were most sensitive to all other classes, especially to Gammaproteobacteria, which showed the highest aggregated frequency and average intensity of inhibitory activities ( ). Since we observed taxonomic signals at the class level, we calculated the average halo of inhibition size for all target and producer isolates across all bacterial classes tested (see sensitivity scores in and inhibition scores in , see also SI Appendix , Fig. S1 B ). This revealed that only a few bacteria – mainly belonging to Pseudomonadaceae (R9, R68, R71, R329, R401, R562, and R569) – exhibited broad inhibition of phylogenetically diverse isolates ( ). Extensive strain-specific variation in inhibition scores across closely related bacteria was observed, suggesting large standing genetic variation for the production of bacterial exometabolites among root-derived isolates ( ). Finally, we observed a 2.7× higher inhibitory activity for root-derived bacteria compared to those originating from soil ( P = 0.011; Kruskal–Wallis followed by Dunn’s post-hoc test and Benjamini–Hochberg correction, ), suggesting that the production of exometabolites might be advantageous for bacterial root colonization. Ralstonia spp . are core members of the root microbiota of healthy Arabidopsis plants in nature ( ). We hypothesized that the bacterial root microbiota of A. thaliana contributes to preventing disease in natural environments and tested the inhibitory activities of 167 of the aforementioned bacteria against pathogenic Ralstonia solanacearum GMI1000 (hereafter referred to as Rs ; ref. and ). A subset of root- and soil-derived bacteria inhibited the growth of Rs in mBA experiments (10.9% and 10.3%, respectively). These inhibitory activities were mainly manifested by bacteria from three genera: Pseudomonas , Streptomyces, and Bacillus ( ). 

We hypothesized that the ability to produce inhibitory halos is correlated with the genome-encoded potential for the biosynthesis of specialized metabolites. We predicted BGCs for all strains tested in mBA experiments ( SI Appendix , Table S2 ) and examined whether halo producer strains encode more BGCs than nonproducers. The total number of BGCs was significantly increased in the antagonistic isolates ( P = 0.0166). BGCs, which are thought to be involved in the biosynthesis of nonribosomal peptides (NRPs), aryl polyene, and redox cofactors, were significantly enriched in halo-producing isolates compared with nonproducers ( ), suggesting that these compounds may be important for interbacterial competition. We next investigated the diversity of bacterial metabolites that could explain the observed inhibitory phenotypes. Metabolites were extracted from individual strains (n = 198) grown on the same agar medium used in mBA experiments with two organic solvents with different polarities, ethyl acetate and methanol, to capture greater chemical diversity. Liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) analysis of the resulting 396 samples yielded ~200,000 mass spectra, which were analyzed using the Global Natural Products Social Molecular Networking (GNPS) workflow ( ). Network analysis of the resulting mass spectra revealed 247 families of molecules with similar fragmentation patterns, each representing at least two structurally related analogs ( SI Appendix , Fig. S2 ). While 2,220 nodes were shared between multiple classes, the remaining 1,094 nodes were found to be produced only by individual bacterial classes. Here, Gammaproteobacteria showed the greatest number of class-specific metabolites, indicating that their pronounced inhibitory activity ( ) can be explained by the production of an enormous diversity of specialized metabolites. To determine whether bacteria secrete an increased diversity of specialized metabolites when interacting with competitors ( ), we additionally analyzed the metabolomes of 16 inhibitory zones from 10 producer isolates (R63, R71, R68, R342, R401, R562, R569, R690, R920, and R1310) with broad inhibitory activity against three target isolates (R472D3, R480, and R553) using the same UPLC-MS/MS GNPS workflow. Spectra comparison with those of individually cultured strains revealed 298 additional ions detected exclusively in inhibition zones ( ). Dereplication of the network revealed that the biosynthesis of polyketides, such as two additional congeners of the nactin antibiotics and peptidic compounds, such as cyclo(Trp-Pro), and congeners of a unique lipopeptide family, were specifically triggered by interactions with other bacterial strains, pointing to a possible role of these molecules in interbacterial competition. 

P. brassicacearum R401. P. brassicacearum R401 showed the greatest number of inhibitory interactions (>17× greater than average) and the largest average halo size (>10× greater than average; ). Except for Bacilli and Flavobacteria, 50 isolates belonging to all other tested bacterial classes were sensitive to the inhibitory activity of R401, particularly Actinobacteria ( ), indicating the production of exometabolites with broad spectrum inhibitory activity. antiSMASH-based analysis of a resequenced, circular R401 genome predicted 16 BGCs ( SI Appendix , Fig. S3 and Table S2 ), one of which exactly matches the phl operon of Pseudomonas protegens Pf-5 , which has been shown to encode the enzymatic machinery for the production of the Pseudomonas -specific polyketide (DAPG; ) ( ). Since no other genome from our culture collection harbors the phl operon ( SI Appendix , Fig. S3 ), this points to a role for DAPG in mediating the inhibitory activity of R401. We generated a marker-free deletion mutant of the key biosynthetic gene phlD in R401 ( ). In mBA experiments, the R401 Δ phld mutant was significantly—yet only partly—impaired in its antagonistic activity ( ), retaining 71% of its inhibitory activity toward Rs . As DAPG production was completely lacking in the Δ phld mutant ( ), this suggests that DAPG alone is insufficient to explain the full inhibitory activity of R401 toward Rs. 

We performed a forward genetic screen to reveal additional determinants mediating the residual inhibitory activity of R401 Δ phld . First, we generated a R401 mini- Tn5 transposon library with >6,000 insertion mutants. We adopted a fluorescence-based bacterial co-culture system in liquid medium to test all the R401 insertion mutants individually for their ability to suppress growth of a GFP-expressing Rs GMI1000 strain ( Rs GMI1600; ref. ). Of the 230 candidates identified in the primary screen, 38 R401 Tn5 mutants were robustly impaired in Rs suppression after two independent rounds of validation ( and SI Appendix , Table S3 ). One of these candidate R401 mutants was also significantly impaired (on average 26.6%) in its inhibitory activity against Rs in mBA experiments on solid agar, suggesting that this mutant is impaired in the production of an inhibitory exometabolite ( SI Appendix , Table S3 ). This R401 mutant carries the Tn5 transposon insertion within the gene of a putative acyltransferase ( pvdY ) involved in the biosynthesis of the siderophore pyoverdine in Pseudomonas aeruginosa ( ). We validated the contribution of R401 pvdY to pyoverdine biosynthesis and Rs inhibition by generating an independent pvdY deletion mutant, Δ pvdy , and a R401 deletion mutant of the gene encoding NRP synthetase pvdL, located downstream of pvdY (Δ pvdl , ). In mBA experiments with Rs GMI1000, the R401 Δ pvdy strain phenocopied the R401 tn5::pvdy mutant and the mutant phenotype was complemented by the expression of pvdY under its native promoter in the R401 Δ pvdy background (Δ pvdy::pvdY ). R401 Δ pvdl exhibited a slightly weaker impairment of halo production compared to R401 Δ pvdy . It is likely that R401 PvdY is involved in hydroxyornithine acetylation—a component of R401 pyoverdine—while PvdL is involved in the initial amino acid condensation, forming the basis for the peptidic backbone of pyoverdines. Acetyl hydroxyornithine could be involved in the biosynthesis of a third, unknown exometabolite by R401, explaining the difference between both pyoverdine mutants observed in mBA experiments. We also generated a R401 double deletion strain, Δ pvdy Δ pvdl , which is impaired in its inhibitory activity against pathogenic Rs to a similar extent as the R401 Δ pvdy single mutant ( ). Using mass spectrometry, we confirmed that all generated R401 pyoverdine mutants had lost their ability to produce pyoverdine (R401 pyoverdine has the sequence: Glu-Q-Lys-AcOHOrn-Ala-Gly-Ser-Ser-OHAsp-Thr, with Q being the fluorophore moiety; ). All tested Pseudomonas isolates in the At -RSphere contain pyoverdine biosynthetic genes; however, for R9, no characteristic pyoverdine fluorescence was detected ( SI Appendix , Fig. S3 ), potentially explaining its low inhibitory activity in mBAs ( ). Therefore, we hypothesized that pyoverdines contribute more broadly to explaining the unusually high inhibitory activity of the Pseudomonadaceae detected in the mBA experiments and isolated pyoverdine mutants from another Pseudomonas root commensal in our collection. We generated a mini- Tn5 mutant library of approximately 6,000 insertion mutants of pyoverdine-producing Pseudomonas fluorescens R569 and assayed about 2,000 mutants for loss of pyoverdine-mediated fluorescence by fluorimetry ( SI Appendix , Fig. S4 A ). Characterization of the Tn5 integration sites revealed one mutant in each of the R401 pvdY and pvdL homologs ( SI Appendix , Table S3 and Fig. S4 B and C ). Unlike the partial loss of inhibitory activity of the R401 Δ pvdy and Δ pvdl single and Δ pvdy Δ pvdl double mutants against Rs , the R569 Δ pvdy and Δ pvdl single mutants both showed a complete loss of inhibitory activity to Rs in mBA experiments ( SI Appendix , Fig. S4 D ). Although the genomes of root commensals R401 and R569 are assigned to different Pseudomonas sublineages ( ) and the arrangement of genes in the pvd operon is different ( and SI Appendix , Fig. S4 C ), it is likely that the corresponding pyoverdines are directly responsible for the observed inhibitory activity. Our findings indicate that pyoverdine might either be the sole or one of the several exometabolites produced by root commensals that limit the growth of Rs . Since pyoverdines function as iron chelators, we tested whether their inhibitory activity can be explained by interbacterial competition for iron by supplementing mBA agar medium with excess ferric iron (100 µM FeCl 3 ). This resulted in significantly impaired Rs inhibition when confronted with WT R401 or undetectable inhibition when confronted with WT R569, thus phenocopying results obtained with the corresponding strain-specific pyoverdine mutants ( and SI Appendix , Fig. S4 D ). Given that the mutants also exhibit drastically reduced iron mobilization capacity compared with the corresponding WT isolates ( and SI Appendix , Fig. S4 E ), we conclude that their inhibitory activities against root commensals and pathogenic Rs are likely mediated through their iron chelator function. However, we cannot exclude a possible additional pyoverdine activity under limited iron conditions. To test whether these two classes of compounds account for the full inhibitory activity of R401, we generated two R401 double pyoverdine and DAPG mutants, Δ pvdy Δ phld and Δ pvdl Δ phld , and conducted mBA assays using Rs and a representative set of the aforementioned commensals that we previously observed to be sensitive to R401 ( SI Appendix , Table S1 ). The double pyoverdine and DAPG mutants showed severely reduced inhibitory activity against Rs, and all other tested isolates compared to the single mutants, suggesting a cumulative impact of the two metabolites. Although the general sensitivity patterns were again isolate specific, Actinobacteria were typically inhibited by DAPG alone. DAPG and pyoverdine collectively explained >70% of R401 inhibitory activity, but a residual halo of inhibition was still observed, pointing to at least a third—yet to be defined—exometabolite ( ). Beyond the known inhibitory activities of DAPG and pyoverdine, our results provide evidence that metabolites with distinct modes of action jointly act to inhibit a taxonomically broad range of bacteria. 

We conducted root microbiota reconstitution experiments with germ-free A. thaliana Col-0 in Flowpots to study the contribution of DAPG and pyoverdine to root microbiota structure and to limiting the growth of Rs ( ). We tested heat-killed (HK) or live WT R401, R401 Δ pvdy, and Δ pvdl single as well as Δ pvdy Δ phld and Δ pvdl Δ phld double mutants on a phylogenetically diverse synthetic community (SynCom) that comprises 18 bacterial isolates ( SI Appendix , Fig. S5 A and Table S1 ), in both “soil” (peat matrix) and root compartments. After 3 weeks of plant—microbe cocultivation, DNA was extracted from soil and root samples and was subjected to 16S rRNA amplicon sequencing at isolate-specific resolution. Shoot fresh weight did not differ between conditions, and no wilting was observed ( SI Appendix , Fig. S5 B and C ), indicating that Rs was unable to cause disease on A. thaliana , possibly due to the presence of the SynCom. This is consistent with very low Rs relative abundances found in root samples ( SI Appendix , Fig. S5 D ) and a previous report that direct immersion of A. thaliana roots with a very high Rs inoculum of 10 8 cells/mL was needed to induce wilting symptoms in the crucifer ( ). The addition of live, WT R401 to the SynCom significantly reduced bacterial alpha diversity (Shannon index) compared to the (HK) R401 condition in the root compartment ( P < 0.001). This major impact on the community is gradually lost when R401 exometabolite mutants were co-inoculated with the SynCom ( ). Importantly, no such effects were observed in the soil compartment ( ). Bacterial beta diversity (Bray—Curtis dissimilarity) was also drastically affected by WT R401 inoculation, explaining most of the variation and resulting in a clear separation along axis 1 (HK versus WT, P < 0.001; R2 = 0.63). All mutant samples fall in between these extremes and follow a clear trajectory, WT > single mutants > double mutants > HK, suggesting that DAPG and pyoverdine have a cumulative influence on bacterial community structure in the root compartment ( ). This is likely due to direct exometabolite activities, as the production of either metabolite did not affect plant phenotypes ( SI Appendix , Fig. S5 C ). Furthermore, in silico depletion of R401 16S rRNA sequence reads leads to similar changes in beta diversity with even higher significance levels in some cases (Δ pvdl versus WT; P = 0.002; R2 = 0.14; SI Appendix , Fig. S5 E ), indicating that, irrespective of R401 abundance, the bacterial community is altered by DAPG and pyoverdine production. While in the soil compartment, inoculation of WT R401 leads to a similar shift as in the root compartment (HK versus WT, P < 0.001; R2 = 0.66), genetic depletion of either DAPG or pyoverdine biosynthetic abilities alone or together does not alter community structure ( ), which indicates niche-specific activity of these exometabolites in the root compartment. The impact of R401 on community structure was still significant upon in silico depletion of R401 16S rRNA sequences (HK versus WT , P < 0.001; R2 = 0.35; SI Appendix , Fig. S5 F ), indicating that R401 uses other, unknown mechanisms to influence community structure in soil. To determine whether the mBA results obtained in vitro have physiological relevance in planta , we inspected relative isolate abundances across conditions, reasoning that isolates insensitive to at least one of the R401 exometabolites in vitro would not benefit from genetic disruption of either BGC in a community context in the root compartment. This analysis revealed that only DAPG- and/or pyoverdine-sensitive isolates benefited from disruption of the corresponding BGCs present in R401 ( ), suggesting that binary interaction data in vitro can explain the impact on individual commensal isolates in a community context in the root compartment. This observation prompted us to test whether mBA data could predict community-scale effects as well as effects on single isolates. We computed the average reduction in inhibitory activity for each R401 mutant ( ) and tested whether the lack of competitiveness could inform changes in alpha- and beta-diversity. Linear regression analyses revealed that binary interaction data largely explained the observed effect size in the root compartment for both alpha- and beta-diversity indices but had no predictive power in the soil compartment ( , respectively). In conclusion, two exometabolites produced by a single isolate have large effects on key ecological indices of a taxonomically diverse SynCom and can be linked to pairwise in vitro interaction experiments. 

To examine whether R401-induced modulation of bacterial assembly and diversity promoted R401 competitiveness, we determined the relative abundance of R401 in SynCom samples collected from root and soil compartments. In both root and soil samples, the 16S rRNA reads of live R401 by far exceeded the barely detectable HK R401 reads, indicating that live R401 inoculum proliferates in both compartments under all conditions ( ). However, R401 accumulated in association with roots at >2x higher abundance than in the soil compartment. R401 accumulation in the root compartment was gradually reduced in the SynCom context: Single or double mutations of DAPG and pyoverdine biosynthetic genes were sufficient to reduce the R401 abundance by up to approximately 39% compared to WT R401 ( ). Importantly, the capacity of these deletion mutants to colonize the soil compartment remained unaffected ( ). We also investigated the root colonization capacity of WT and all the five R401 deletion strains in mono-associations on axenic A. thaliana plants in an agar-based system ( ) and found no significant differences in live R401 cell counts ( ). Neither the growth of R401 DAPG or pyoverdine mutants was impaired in axenic culture media ( SI Appendix , Fig. S6 B and C ). Similarly, colonization experiments with WT R401 or the Δ pvdy Δ phld double mutant in the Flowpot system revealed no difference in root colonization capacity in the absence of bacterial competitors ( SI Appendix , Fig. S6 A ). Thus, DAPG and pyoverdine cofunction as R401 root competence determinants specifically in competition with other members of this SynCom in the root compartment. To assess the prevalence of Pseudomonas strains capable of producing DAPG and pyoverdine in association with plants when grown in natural soils, we analyzed the genomes of several commensal culture collections of 1,567 isolates from roots or leaves of A. thaliana or roots of the legume Lotus japonicus ( , , , ). These plants had been grown in the same or different soils on different continents, and our genome analysis revealed increased abundance of DAPG and pyoverdine BGCs in root-derived Pseudomonas isolates ( SI Appendix , Fig. S6 D and E ). These results not only support our conclusions obtained with a defined core commensal community and gnotobiotic A. thaliana but also suggest a broader yet root-specific role of these two exometabolites in natural environments beyond the model crucifer. 

Our identification of a subset of specialized metabolites, primarily polyketides, and NRPs, that are specifically produced upon competition sensing, suggests that interbacterial competition can activate cryptic BGCs in the A. thaliana root microbiota. A related study using an A. thaliana phyllosphere bacterial culture collection also demonstrated widespread production of specialized metabolites and identified 725 binary inhibitory interactions in vitro, with only 3.6% of 224 strains mediating most of these interactions ( ). The A. thaliana leaf and root microbiota overlap extensively at higher taxonomic ranks, and Pseudomonas species are core members of both communities ( ). However, the comparison between root- and leaf-derived Pseudomonas genomes collected from natural environments reveals complete or partial niche specificity, respectively, for Pseudomonas production of DAPG and pyoverdine in the root compartment ( SI Appendix , Fig. S6 D ). Furthermore, the total number of BGCs is significantly higher (approximately 1.6×) in root-derived Pseudomonas isolates ( SI Appendix , Fig. S6 E ). Likewise, the number of antagonistic interbacterial interactions within the root microbiota significantly exceeds those between abundant soil-derived bacteria ( ). These findings could be explained by the fact that natural, unplanted soils and leaves represent oligotrophic habitats ( , ), whereas the rhizoplane is a microenvironment with a continuous supply of root exudation-derived nonstructural carbohydrates to support bacterial proliferation ( ). It is likely that the large number of antagonistic interactions we found between root microbiota members in vitro is also influenced by the agar medium, which is rich in nonstructural organic carbon to mimic the nutrient-rich rhizoplane, thereby inducing the costly biosynthesis of exometabolites. We provided genetic evidence for a root niche-specific cofunction of the exometabolites DAPG and pyoverdine as root competence determinants of commensal R401 that scales at the community level by influencing alpha- and beta-diversity indices. R401 also proliferates in the soil compartment, but its growth there and its influence on the structure of the bacterial soil community are independent of these exometabolites ( and ). These observations, together with the fact that plant roots are a major sink for the uptake of rhizospheric mineral iron ( – ), imply that the public good iron becomes rate limiting in the root compartment. This could explain why the high-affinity iron chelator pyoverdine and antimicrobial DAPG cofunction and maximize R401 growth at the expense of its commensal competitors in the root niche, despite their different modes of action. Consistent with this model, the expression of core biosynthetic genes involved in the production of both DAPG and pyoverdine is induced under iron-limiting conditions in rhizospheric P. fluorescens and human pathogenic P. aeruginosa PAO1 ( , ). The indistinguishable root colonization capacity of WT R401, pyoverdine, or DAPG single, or double mutants in mono-associations with A. thaliana is also consistent with our model and further supports the specific role of these exometabolites in interspecies competition during root microbiota establishment. Likewise, P. fluorescens C7R12 pyoverdine mutants were unaffected in rhizosphere competence under axenic conditions ( ). All the seven A. thaliana root-derived Pseudomonas spp. strains in the At -RSphere culture collection have the genetic capacity to produce pyoverdines ( SI Appendix , Fig. S3 ). Together with the essential function of pyoverdines produced by multiple Pseudomonas root commensals in limiting the growth of pathogenic Rs ( , SI Appendix , Fig. S4 D , and ref. ), this suggests a widespread role for these siderophores in determining Pseudomonas competitiveness in the root microbiota, even in soils that contain replete bioavailable inorganic iron for plant growth ( , , , ). Thus, it is possible that in natural soils, regulated secretion of pyoverdines is an adaptive trait of the genus Pseudomonas to simultaneously compete against the root iron sink and enable pervasiveness of the Pseudomonadaceae lineage via interbacterial competition ( , , ). Our genetic data show that a sequential reduction in the diversity of specialized exometabolites in one species is sufficient to increase alpha-diversity indices and significantly alter the beta diversity of a synthetic root microbiota, thereby establishing a causal link between within-species genetic diversity and interspecies diversity changes. This suggests that in the root microenvironment, DAPG and pyoverdine produced by WT R401 cofunction to locally inhibit the growth of multiple SynCom members, which in turn increases the abundance of the producer R401. Our results obtained with the SynCom in coculture with gnotobiotic A. thaliana grown in peat matrix bear striking similarity to the changes in alpha diversity reported for this plant when grown in natural soils, where bacterial alpha diversity in bulk soil, rhizosphere, and root compartments gradually decreases toward the root ( , , ). This suggests that exometabolite-mediated antagonistic interactions of commensals underpin at least part of the reduction of alpha diversity consistently observed in natural environments at the soil—root interface. Given the prevalence of operons predicted to encode pyoverdines among root-derived Pseudomonas spp. isolates in the At -RSphere culture collection (ref. and SI Appendix , Fig. S3 ), we consider the moderately complex 18-member SynCom employed here as necessary to overcome genetic redundancy, at least for pyoverdine production and probably also for the DAPG exometabolite in the root microbiota, when A. thaliana plants are grown in natural soils ( SI Appendix , Table S2 and refs. , ). It remains to be tested whether genetic depletion of individual bacterial exometabolites can have similarly striking effects on highly complex natural microbial communities, as the assessment of microbial antagonism in vitro is often not transferable to field conditions ( – ). However, DAPG- and phenazine-overproducing mutant strains of Pseudomonas putida WCS358r were shown to differentially shape root-associated fungal communities of field-grown wheat over an experimental period of up to 139 d compared to the WCS358r WT strain ( , ). This suggests that bacterial exometabolites—such as DAPG—can have long-lasting effects on the structure of complex microbial communities. The approach described here to remove genetic redundancy by using annotated genomes of cultured microbiota members may be more generally applicable to explore community functions of other microbial genetic determinants with SynComs in gnotobiotic plant growth systems. Taken together, our study suggests that high-throughput binary interaction experiments, combined with genome mining for BGCs of root microbiota culture collections, can be applied to identify strains with broad-spectrum antagonistic activities that are likely robust root colonizers. This might have relevance for future interventions in the root microbiota with rational biologicals that confer beneficial traits on the host, including indirect pathogen protection and mineral nutrition. 

Detailed descriptions of all utilized methods and data analysis workflows can be found in SI Appendix . Screen for Antagonistic Interbacterial Interactions. For all mBA experiments, bacterial strains were cultured axenically. Target strains were embedded in molten 25% tryptic soy agar (TSA) and producer strains were then dropped out on top. After up to 96 h of cultivation, pictures were taken and halos were quantified. Metabolomic Analyses. Untargeted metabolomic analysis was either conducted with axenically grown strains or on the mBA inhibition zones of ten strains (R63, R68, R71, R342, R401, R562, R569, R690, R920, and R1310), each tested against three target strains (R472D3, R480, and R553). Detection of R401 DAPG and Pyoverdine. Using metabolite analyses, R401 DAPG and pyoverdine were detected in WT extracts and lack of metabolite production was confirmed in the respective mutants. BGC Prediction Using antiSMASH. antiSMASH 6.0 ( ) was used to predict BGCs for all the tested strains. Mutant Generation. Using homologous recombination, marker-free targeted R401 knockout mutants were generated lacking DAPG and/or pyoverdine biosynthetic genes. Establishment of mini- Tn5 Transposon Mutant Collections in R401 and R569. R401 and R569 were transformed with pUTmTn5Km2 ( ), carrying a mini- Tn5 transposon, and individual colonies were then picked into 96-well plates. mini- Tn5 Transposon Mutant Screen for Loss of R401-Mediated Growth Inhibition of Rs GMI1600. R401 mini- Tn5 transposon mutants were evaluated in parallel for their WT-like growth and loss of inhibitory activity against GFP-expressing Rs GMI1600. mini- Tn5 Transposon Mutant Screen for Lack of Pyoverdine Fluorescence of R569. Axenically grown R569 mini- Tn5 transposon mutants were evaluated for their lack of fluorescence characteristic of pyoverdine at λ excitation = 395 nm and λ emission = 470 nm. Identification of mini- Tn5 Transposon Integration Sites in the Genomes of R401 and R569. The chromosomal mini- Tn5 transposon integration sites in the R401 or R569 genomes were determined similarly as described before ( ). Complementation of R401 Δ pvdy . Complementation of R401 Δ pvdy was conducted by expressing the coding region of R401 pvdY with its native 5′ regulatory sequences from the low-copy pSEVA22l plasmid ( ) in the Δ pvdy mutant background. In vitro Iron Mobilization Assay. The capability of R401 and R569 mutants to solubilize inaccessible ferric iron was tested using a previously described photometric assay ( ). Validation of Bacterial Growth Rates. The growth of each R401 mutant and wild type was assessed by continuously measuring the OD 600 of actively growing bacterial cultures. Microbiota Reconstitution in the Gnotobiotic Flowpot System. Flowpots were assembled according to ref. with minor modifications. A final bacterial OD 600 of 0.0025 was inoculated in ½ MS, and five surface-sterilized A. thaliana Col-0 seeds per Flowpot were sown. After 21 dpi, the fresh weight of shoots was determined, and the roots and peat matrix were harvested for bacterial community profiling. Monoassociation Experiment of R401 on A. thaliana Seedlings. Colony-forming units of WT R401 or its mutants colonizing A. thaliana seedlings grown in ½ MS agar were determined as described ( ). DNA Isolation. DNA was isolated from A. thaliana roots and Flowpot peat using a modified high-throughput version of the FastDNA SPIN kit for Soil (MP Biomedicals). Library Preparation for Bacterial 16S rRNA Gene Profiling. The v5v7 variable regions of the bacterial 16S rRNA gene were amplified from DNA template derived from roots and Flowpot peat. The resulting amplicons were tagged with sample-specific barcodes using a dual-indexing approach. Illumina paired-end sequencing was then performed in-house with the MiSeq benchtop sequencer. 

For all mBA experiments, bacterial strains were cultured axenically. Target strains were embedded in molten 25% tryptic soy agar (TSA) and producer strains were then dropped out on top. After up to 96 h of cultivation, pictures were taken and halos were quantified. 

Untargeted metabolomic analysis was either conducted with axenically grown strains or on the mBA inhibition zones of ten strains (R63, R68, R71, R342, R401, R562, R569, R690, R920, and R1310), each tested against three target strains (R472D3, R480, and R553). 

Using metabolite analyses, R401 DAPG and pyoverdine were detected in WT extracts and lack of metabolite production was confirmed in the respective mutants. 

antiSMASH 6.0 ( ) was used to predict BGCs for all the tested strains. 

Using homologous recombination, marker-free targeted R401 knockout mutants were generated lacking DAPG and/or pyoverdine biosynthetic genes. 

Tn5 Transposon Mutant Collections in R401 and R569. R401 and R569 were transformed with pUTmTn5Km2 ( ), carrying a mini- Tn5 transposon, and individual colonies were then picked into 96-well plates. 

Tn5 Transposon Mutant Screen for Loss of R401-Mediated Growth Inhibition of Rs GMI1600. R401 mini- Tn5 transposon mutants were evaluated in parallel for their WT-like growth and loss of inhibitory activity against GFP-expressing Rs GMI1600. 

Tn5 Transposon Mutant Screen for Lack of Pyoverdine Fluorescence of R569. Axenically grown R569 mini- Tn5 transposon mutants were evaluated for their lack of fluorescence characteristic of pyoverdine at λ excitation = 395 nm and λ emission = 470 nm. 

Tn5 Transposon Integration Sites in the Genomes of R401 and R569. The chromosomal mini- Tn5 transposon integration sites in the R401 or R569 genomes were determined similarly as described before ( ). 

pvdy . Complementation of R401 Δ pvdy was conducted by expressing the coding region of R401 pvdY with its native 5′ regulatory sequences from the low-copy pSEVA22l plasmid ( ) in the Δ pvdy mutant background. 

The capability of R401 and R569 mutants to solubilize inaccessible ferric iron was tested using a previously described photometric assay ( ). 

The growth of each R401 mutant and wild type was assessed by continuously measuring the OD 600 of actively growing bacterial cultures. 

Flowpots were assembled according to ref. with minor modifications. A final bacterial OD 600 of 0.0025 was inoculated in ½ MS, and five surface-sterilized A. thaliana Col-0 seeds per Flowpot were sown. After 21 dpi, the fresh weight of shoots was determined, and the roots and peat matrix were harvested for bacterial community profiling. 

thaliana Seedlings. Colony-forming units of WT R401 or its mutants colonizing A. thaliana seedlings grown in ½ MS agar were determined as described ( ). 

DNA was isolated from A. thaliana roots and Flowpot peat using a modified high-throughput version of the FastDNA SPIN kit for Soil (MP Biomedicals). 

16S rRNA Gene Profiling. The v5v7 variable regions of the bacterial 16S rRNA gene were amplified from DNA template derived from roots and Flowpot peat. The resulting amplicons were tagged with sample-specific barcodes using a dual-indexing approach. Illumina paired-end sequencing was then performed in-house with the MiSeq benchtop sequencer. 

Appendix 01 (PDF) Click here for additional data file. 

