Metabolic exchange and siderophore sharing underlie emergent biofilm synergism

Metabolic exchange and siderophore sharing underlie emergent biofilm synergism | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Metabolic exchange and siderophore sharing underlie emergent biofilm synergism Ákos T Kovács, Jiyu Xie, Xinli Sun, Viktor Hesselberg-Thomsen, and 8 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-8956555/v1 This work is licensed under a CC BY 4.0 License Status: Under Review Version 1 posted You are reading this latest preprint version Abstract Biofilm communities exhibit emergent properties that exceed the sum of contributions from individual members of the community. Here, we describe a multilayered metabolic interaction that drives enhanced biofilm formation among three bacterial species from the plant rhizosphere. Comparative metatranscriptomic and metabolomic analyses reveal that Bacillus velezensis-secreted 5-aminovaleric acid promotes the growth of the other community members, Burkholderia contaminans and Acinetobacter baumannii. In return, B. contaminans supplies branched-chain amino acids for B. velezensis. Branched-chain amino acids and cell–cell signaling acyl-homoserine lactones from B. contaminans induce biosynthesis of the siderophore bacillibactin in B. velezensis, that is further enhanced by A. baumannii. In exchange, the B. velezensis-secreted siderophore promotes the growth of B. contaminans in iron-limited conditions, which benefits the multispecies biofilm community in vitro and promotes plant growth performance in iron-depleted soil. Our study reveals the molecular mechanisms underlying an emergent rhizosphere biofilm community function and demonstrates its importance in plant–microbe interactions. Biological sciences/Microbiology/Microbial communities/Microbial ecology Biological sciences/Ecology/Microbial ecology multispecies biofilm microbial community siderophore cross-feeding plant growth promotion emerging function Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Introduction Biofilms are complex aggregates of microorganisms attached to surfaces or to each other, composed of microbial cells embedded in extracellular polymeric substances. Biofilms are one of the primary modes of microbial existence in natural environments, including soil 1,2 . This lifestyle provides numerous advantages to microorganisms, including improved nutrient acquisition, enhanced resistance to environmental stresses, and facilitation of horizontal gene transfer. In nature, biofilms typically contain multiple species; they have emergent community properties that cannot be predicted from the individual species alone and exceed the sum of contributions of individual member species of the community 3,4 . Understanding and ultimately predicting emergent properties of biofilm communities remains a key challenge in microbial ecology 5 . Although it is evident that emergent properties arise from interactions within the community, how these emergent features develop in biofilm communities and how they contribute to the advancement of biofilm populations or even influence their hosts remain poorly understood 1,6 . For example, multispecies biofilms in the rhizosphere can provide plants with enhanced drought tolerance and pathogen suppression 7,8 . These experimental observations suggest that rhizosphere biofilm communities can influence plant growth and health, yet a definitive link between biofilm community emergent properties and plant performance is yet to be established. Metabolic cross-feeding is a widespread interaction in microbial communities and is particularly important in resource-limited environments 9,10 . In multispecies biofilms, spatial proximity creates favorable conditions for metabolic cross-feeding, whereas the matrix-embedded structure minimizes diffusion losses and decreases the risk of metabolite exploitation by free-living microorganisms 11 . Amino acid cross-feeding, particularly of branched-chain amino acids (BCAAs), is a common form of metabolic cross-feeding among bacteria 12 . BCAAs not only serve as nutritional resources but also act as regulatory molecules of nonribosomal peptide synthesis 13–15 . Other compounds, including vitamins, heme, organic acids, and sugars, have been implicated in microbial cross-feeding 16–19 . Metabolite cross-feeding underpins the stability of biofilm communities. However, their rapid development and expansion may necessitate the management of other limiting resources. Iron is an essential element for nearly all microbial life activities 20 , participating in important processes such as electron transfer, DNA synthesis, and cellular respiration 21,22 . However, in most natural environments, iron exits primarily as insoluble Fe (III) oxides with extremely low bioavailability 23 . This limitation is particularly problematic in agricultural soils, especially in alkaline conditions, where iron deficiency significantly constrains both microbial activity and plant growth 24 . To address this challenge, microorganisms produce and secrete high-affinity iron chelators called siderophores as a widespread, effective iron acquisition strategy 25,26 . In the rhizosphere, microbial siderophores not only support bacterial iron acquisition but also improve plant iron nutrition 27,28 . For example, bacillibactin, the siderophore produced by Bacillus , promotes plant iron acquisition in alkaline soil 29 . Bacillibactin not only affects iron uptake but also influences the secretion of extracellular matrix and biofilm formation by Bacillus 30 . For wild-type (WT) strains, the iron content in the extracellular matrix is 5-to-10-times higher than that inside the cells. In contrast, mutants with disrupted bacillibactin synthesis genes exhibit a significant decrease in extracellular matrix secretion capability and lose the ability to enrich iron in the extracellular matrix 31,32 . Our latest research indicates that bacillibactin also acts as a signaling molecule that activates plant iron uptake, thereby enhancing iron absorption in iron-limited conditions 29 . Traditionally, siderophores have primarily been regarded as mediators of interspecies competition 33 , where microorganisms compete for limited iron resources through the production of siderophores with different affinities, or by exploiting the siderophores produced by competitors. This competitive paradigm has been well-documented in pathogens, such as Pseudomonas aeruginosa and Staphylococcus aureus , which engage in fierce competition for limited iron resources 34,35 . Similar competition occurs in soil, where siderophore-producing Pseudomonas can inhibit competitors through iron sequestration 36 . Additionally, some microorganisms can directly use siderophores produced by other species, which has been observed in interactions between P. aeruginosa and Mycobacterium species 37 . The interspecies exploitation of siderophores raises intriguing questions: Can siderophore sharing promote cooperation among microorganisms? Additionally, if siderophores function as a public good induced within biofilm communities, can these specialized metabolites also act as signals in the rhizosphere to modulate plant–microbe interactions? This study explored the community-level emergent properties of a trispecies community of soil bacteria consisting of Acinetobacter baumannii XL380 (Ab), Burkholderia contaminans XL73 (Bc), and Bacillus velezensis SQR9 (Bv), that exhibits synergistic biofilm formation. We discovered a signaling and metabolite sharing network, including cross-feeding between Bv and the other two species through 5-aminovaleric acid (Bv to Bc/Ab) and BCAAs (Bc to Bv). Simultaneously, Bc secretes signaling molecules, including acyl-homoserine lactones (AHLs), that upregulate bacillibactin biosynthesis in Bv, which is exploited by Bc, explaining the enhanced biofilm community productivity. The synergistic interaction among the three species promotes plant growth performance in iron-deficient soil. We reveal a reciprocal mutualism that challenges traditional competitive interaction models in bacteria, and demonstrate how metabolic complementation can transform resource limitation into community-wide benefits. Results Synergistic biofilm formation by the trispecies community We have previously collected bacterial isolates from cucumber rhizosphere soil and performed static cultivation, which demonstrated enrichment of Ab, Bc, and Bv in pellicle biofilms 38,39 . Biomass quantification revealed that a trispecies biofilm formed by Ab, Bc, and Bv (AbBcBv) had significantly higher fresh weight than any mono- or dual-species biofilms (Fig 1A). Multispecies biofilms exhibited more robust and wrinkled morphology than single-species pellicle biofilms. To determine whether the enhanced biofilm in the coculture resulted from mutualism or competition, we quantified the cell numbers of each species. In the trispecies biofilm, Bv was the most abundant species (~2.2×10 8 ), followed by Bc (~7.5×10 7 ), and Ab (~1.3×10 7 ). The observation was confirmed by scanning electron microscopy (SEM), which showed a low abundance of Ab closely associated with abundant Bc and Bv (Fig 1B, Fig S1). The trispecies combination (AbBcBv) demonstrated the highest total cell number among all groups, significantly exceeding those in any other culture condition. The cell numbers of Bc (Fig 1D) and Bv (Fig 1E) were significantly higher than in their respective monocultures, whereas Ab decreased compared with its monoculture (Fig 1C). Pairwise comparisons revealed variable interaction patterns. Both AbBv and BcBv showed significantly higher cell numbers than their respective monocultures, whereas AbBc showed no significant difference compared with Bc alone (Fig 1A, right axis, p < 0.05). These findings suggest the enhanced trispecies biofilm results from complex interactions involving both mutualism and competition, where Bc and Bv benefit from the multispecies environment whereas Ab, though decreased in amount, remains an integral member of the community. Transcriptomic analysis of the trispecies biofilm To elucidate the molecular details governing synergistic interactions, global gene expression profiles were systematically investigated. Compared with their monocultures, Bc and Bv differentially expressed >1000 and hundreds of genes, respectively, in BcBv dual- and AbBcBv trispecies biofilms (Fig 2A). Lists of the differentially expressed genes in BcBv and AbBcBv biofilms showed substantial overlap, including 1643 genes in Bc and 248 genes in Bv (Fig 2B,C). In contrast, only a few genes were differentially expressed in AbBc and AbBv combinations compared with monocultures (Fig 2A). Therefore, Bc and Bv are the primary interacting species that drive major transcriptomic changes, whereas Ab exerts limited effects on the transcriptional landscapes of the other community members. Functional analysis of multispecies biofilms compared with monocultures revealed that 24 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were differentially regulated ( p < 0.05) (Fig 2D). Several KEGG pathways of Bc were significantly upregulated both in BcBv and the trispecies biofilm, including xylene degradation, two-component systems, sulfur metabolism, and ABC transporters. Conversely, genes encoding the biosynthesis of antimicrobials were significantly downregulated. In Bv, the gene cluster encoding siderophore biosynthesis was significantly upregulated in both BcBv and trispecies biofilms compared with monoculture. The genes encoding the ribosomes of Bv were significantly upregulated in all coculture groups, suggesting enhanced protein synthesis and faster growth, which is consistent with the higher cell number observed in mixed species biofilms. Conversely, genes encoding the biosynthesis of BCAAs (valine, leucine, and isoleucine) were significantly downregulated in BcBv and trispecies biofilms, suggesting decreased BCAAs production by Bv. Given that Bv emerged as the predominant species in biofilm formation, we specifically focused on two pathways that were significantly regulated in multispecies biofilms compared with Bv monoculture: the upregulated siderophore pathway, and the downregulated BCAA synthesis pathway. Real-time quantitative -PCR (RTq-PCR) experiments confirmed that transcription of all genes for biosynthesis of bacillibactin (the siderophore of Bv) was upregulated in the triculture biofilm (Fig 2F, Fig S2A). In contrast, the transcription of all genes for biosynthesis of pyochelin and ornibactin (siderophores of Bc) was downregulated (Fig 2F, Fig S2B). Regarding the siderophores of Ab, 75% of the biosynthesis genes encoding acinetobactin synthesis were upregulated in the triculture biofilm, whereas the transcription of genes encoding the baoumannoferrin biosynthesis pathway remained unaffected. In the BCAA biosynthesis pathway, the ilvC and ilvH genes of Ab and Bc were upregulated in the trispecies biofilm but downregulated in Bv compared with monocultures (Fig 2E, Fig S2C). These expression patterns suggest that Bv may provide siderophore for the biofilm community, whereas Ab and Bc provide BCAAs in return, establishing a potential metabolic complementation relationship. Bc induces bacillibactin production in Bv and uses it for iron acquisition To explore the importance of iron and siderophores in synergistic biofilm formation, we used three mutants of Bv 30 : Δdhb ,which is unable to produce bacillibactin and its precursor 2,3-dihydroxybenzoate (DHBA); ΔdhbF , which produces DHBA but not bacillibactin; and ΔfeuB , which lacks Fe–bacillibactin transport into the cells 40 (Fig S3). Using trispecies cocultures with either WT or mutant strains of Bv, we tested three different iron levels by adding extra FeCl 3 or the iron chelator 2-dipyridyl (2DP) to standard tryptic soy broth (TSB) culture medium (Fig 3A, Fig S3). In low-iron (2DP-supplemented) medium, Co- Δdhb , Co- ΔdhbF , and Co- ΔfeuB exhibited severely impaired growth compared with Co-WT (Co- denotes coculture), including significantly decreased cell numbers (Fig 3A,B). This suggests that siderophore produced by Bv supported the growth of the biofilm community in iron-limited conditions. In high-iron (FeCl 3 -supplemented) medium, all cocultures showed similar biofilm phenotypes, because abundant iron was directly accessible without a requirement for siderophores (Fig 3A). In TSB (intermediate level of iron), the biofilms of Co- Δdhb , Co- ΔdhbF , and Co- ΔfeuB appeared thinner than Co-WT, and they displayed decreased cell numbers (Fig 3A,B), confirming the necessity of Bv siderophore in enhancing trispecies biofilm formation. Biofilms of Co-WT formed in high-iron conditions had significantly lower cell numbers than those in normal iron conditions, suggesting that siderophores play an additional role in the interspecies interaction, beyond iron scavenging. Next, we monitored dhbA expression in WT Bv when it was cocultured with Ab or Bc. As the distance between colonies on an agar plate was decreased, Bc induced dhbA expression in Bv, whereas a comparable effect was not observed when Bv was paired with Ab (Fig 3C). The combination of Ab and Bc induced higher expression than did Bc alone, indicating that although Ab did not directly induce dhbA expression in Bv, it enhanced the interaction between Bc and Bv. To validate Bc-induced bacillibactin production in Bv, matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI–TOF) imaging was used, which confirmed an increased level of bacillibactin in the Bv colony in the presence of Bc (Fig 3D). Whereas bacillibactin production was significantly enhanced, the spatial distribution of other Bv-produced secondary metabolites showed varying responses in the presence of Bc: bacillomycin D and surfactin showed moderate increases, whereas the polyketides difficidin and macrolactin A were not changed. Quantitative analysis of coculture biofilms revealed an elevated bacillibactin concentration in the trispecies biofilm (142.09 ± 17.59 ug/well) compared with Bv monoculture (22.44 ± 0.77 ug/well) (Fig 3E, Fig S4). These results confirmed that the presence of Bc increased bacillibactin production in Bv . Next, we evaluated whether Bv siderophores can support Bc growth at different iron levels (Fig 3F). In low-iron medium (TSB + 2DP), Bc showed no growth when cultured alone or when supplemented with spent medium from culture of the Bv Δdhb strain. However, spent medium from culture of WT Bv restored the growth of Bc to a similar level to that in TSB without iron depletion, demonstrating that Bc can use bacillibactin produced by Bv to overcome iron limitation. We then hypothesized that Bc secretes signaling molecules that induce siderophore production in Bv. Previous studies have demonstrated that AHLs produced by Burkholderia cepacia can regulate siderophore production 41 . We thus tested induction of dhbA expression in Bv by eight AHL homologs (Fig S5). At 12 and 18 h, C4-HSL and C6-HSL significantly induced dhbA expression compared with the negative control, with C4-HSL showing induction levels comparable to that in the positive control (Bc supernatant) by 18 h (Fig S5A,B). High-performance liquid chromatography (HPLC) analysis confirmed the presence of C4-HSL in Bc supernatant (Fig S6). In conclusion, a sophisticated chemical interplay is present within the biofilm community: Bc stimulates bacillibactin production in Bv, possibly through specific AHLs; Bv provides an iron-scavenging siderophore that supports Bc growth in iron-limited conditions; and Ab reinforces this mutualism, establishing a synergistic relationship critical for the enhanced trispecies biofilm formation. Ab and Bc provide BCAAs to support the growth of Bv To validate the hypothesis that Bv exploits BCAAs produced by Ab and Bc in the mixed species biofilm, we cocultured Bv mutants Δ ilvA , Δ ilvCH , Δ ilvD (unable to synthesize all three BCAAs) and Δ leuBCD (unable to synthesize leucine) 42 with Ab and Bc in TSB (Fig 4A). The biomass of the resulting coculture biofilms showed a clear increase compared with monoculture biofilms. Compared with the monocultures, the cell numbers of the WT, Δ ilvA ,Δ ilvD ,and Δ ilvCH strains exhibited a significant increase, whereas strain Δ leuBCD showed minimal growth enhancement (Fig 4B). These results demonstrated that Ab and Bc can provide BCAAs, particularly valine and isoleucine, to support the growth of Bv in the trispecies biofilm, confirming the existence of BCAA cross-feeding between these species. BCAAs have been shown to influence secondary metabolite biosynthesis pathways 43 , which is connected with nonribosomal peptide synthetases (NRPSs) 44 . To examine whether BCAA availability affects bacillibactin production in Bv, bacillibactin levels were quantified (Fig 4C). In monocultures, the BCAA-deficient mutant Δ ilvA showed decreased bacillibactin production compared with the WT strain. Furthermore, direct supplementation of individual BCAAs significantly induced dhbA expression in Bv compared with negative controls, with the positive control (Bc spent medium) showing the highest level of induction (Fig 4D). When the BCAA-deficient mutant was cocultured with Ab and Bc in trispecies biofilms, bacillibactin production was significantly restored (Fig 4C). This restoration of bacillibactin production in the mutant coculture demonstrates that BCAAs secreted by Ab and Bc can be used by Bv to support siderophore biosynthesis, Bv growth, or both at the same time. Bv supplies metabolites that promote the growth of Bc and Ab To explore other potential cross-feeding interactions in the trispecies community, the influence of spent medium on growth was tested (Fig 5A) 42 . During cultivation in M9 medium, Bc depletes glucose, whereas residual glucose remains in Ac culture supernatant. The culture supernatant of Bc did not support growth of the other species (Fig S7). In contrast, spent medium from culture of Bv supported the growth of both Ab and Bc, indicating that Bv provides metabolites to the other two species. The top 30 differential compounds were analyzed based on the relative abundance (Fig 5B). Highly abundant compounds in Bv spent medium that were depleted after cultivation of Ab or Bc were designated as cross-fed compounds. These included valeric acid, 5-aminovaleric acid, pentadecanoic acid, levulinic acid, cinnamic acid, valine, 2,3,4,5-tetramethylpyrazine, and tretinoin. Additionally, species-specific use patterns were observed: nine compounds ( N -α-L-acetyl-arginine, agmatine, N -acetyl-DL-norvaline, N -acetylvaline, valylproline, prolylleucine, N -acetyl-D-alloisoleucine, nicotinic acid, and 2-isopropylmalic acid) were used exclusively by Bc, whereas two compounds (phenylethanolamine and synephrine) were consumed specifically by Ab. Testing the impact of identified compounds on bacterial growth revealed that Ab grew on synephrine, 5-aminovaleric acid, and pentadecanoic acid, whereas Bc used 5-aminovaleric acid and nicotinic acid (Fig S8). Plant growth promotion in iron-deficient soil To evaluate the ecological significance of iron-mediated interaction of the trispecies biofilm community, plant-growth-promoting (PGP) traits of the community were tested. The trispecies community demonstrated superior PGP traits, including enhanced phosphate solubilization (Fig 6A,B), ammonia production (Fig 6C,D), siderophore production (Fig 6E), and exopolysaccharide production (Fig 6F), compared with mono- and dual-species treatments. Growth of Arabidopsis thaliana was severely stunted in the iron-deficient soil without inoculation (the control, CTL), confirming the detrimental influence of iron limitation on plant growth (Fig 6G,H). Plants inoculated with the trispecies consortium (AbBcBv) exhibited the most robust growth and highest fresh weight, significantly outperforming all other treatments (Fig 6G,H). Critically, when the bacillibactin-deficient mutant (Bv-Δ dhb ) was used instead of WT Bv, plant growth was substantially decreased, demonstrating the critical role of bacillibactin in plant iron nutrition in iron-deficient conditions. Dual-species combinations showed intermediate growth promotion effects, with BcBv performing better than AbBv and AbBc. This plant growth promotion pattern correlated with the biofilm productivity, where BcBv demonstrated stronger synergistic interactions than the other dual-species combinations. Discussion Here, we reveal emergent properties of a trispecies bacterial community, including enhanced biofilm formation, metabolic cross-feeding, elevated siderophore production, and improved plant growth promotion (Fig 7). Within the community, Bv serves as the primary metabolic supplier of essential metabolites, such as 5-aminovaleric acid, pentadecanoic acid, nicotinic acid, and synephrine, that support the growth of Bc and Ab, whereas Bc and Ab reciprocally deliver BCAAs. Biofilm synergism is further promoted by elevated bacillibactin production by Bv that promotes iron acquisition by the entire community; the increased bacillibactin production is stimulated both by the signaling molecule C4-HSL from Bc and by BCAAs from Bc and Ab. In iron-deficient conditions, this metabolic complementation not only secures community survival but also confers plant growth promotion effects. Our findings demonstrate that siderophores facilitate interspecies synergism in biofilms and enhance plant iron nutrition. The enhancement of the trispecies biofilm observed in this study represents emergent properties arising from complex interspecies interactions. The significantly higher biomass and cell numbers in AbBcBv treatment compared with monoculture demonstrate such emergent benefits, where the community performance surpasses what would be expected from simply combining the capabilities of individual species. Emergent properties in diverse multispecies biofilms are well-documented, including enhanced pathogen suppression by skin bacterial symbionts and improved stress tolerance in soil microbial communities 45 . Skin bacterial symbionts with greater species richness exhibit superior suppression of the amphibian fungal pathogen Batrachochytrium dendrobatidis through combined dominant effects and complementarity mechanisms. Additionally, probiotic bacterial consortia consisting of up to eight Pseudomonas strains demonstrate emergent consortium-level effects in the tomato rhizosphere, where four- and eight- strain consortia reached densities up to 10-times higher than single-strain inoculants, exhibiting transgressive overyielding that could not be predicted from individual strain performance 46 . These examples demonstrate how community-level emergent properties can provide ecological advantages that individual species cannot achieve alone. These emergent characteristics reflect a complex interplay of spatial organization, metabolic interdependencies, and chemical communication networks within structured microbial communities. Metabolic cross-feeding is a fundamental mechanism that drives microbial community assembly and stability. This phenomenon has been extensively documented across diverse systems, from vitamin interdependencies 17 to amino acid exchange networks in synthetic communities 47 . Metabolic interactions can augment biofilm formation in dual-species systems (e.g., heme cross-feeding in Staphylococcus aureus and Enterococcus faecalis biofilms), and multigenome metabolic modeling predicted functional interdependencies in plant microbiomes 18,48 . Among the three dual-species treatments tested here, BcBv exhibited the strongest synergistic effects, a pattern consistently validated across transcriptomic and metabolomic analyses. Metabolome analysis revealed a biochemical basis underlying the mutualism between Bc and Bv, where Bv functions as the primary metabolite provider, secreting compounds including 5-aminovaleric acid, pentadecanoic acid, nicotinic acid, and synephrine that support community growth. Reciprocally, Ab and Bc provide BCAAs that promote the metabolic processes of Bv, particularly siderophore biosynthesis. BCAA cross-feeding has also been reported in other PGP rhizosphere biofilm communities 7,42 suggesting that amino acid cross-feeding represents a widespread strategy in multispecies biofilms. g-Aminobutyric acid (GABA) can function as a cross-kingdom signal, mediating communication between bacteria, influencing gene expression and behaviors, and thereby modulating the composition and dynamics of microbial communities 49 . We revealed that a Bv-secreted GABA homolog, 5-aminovaleric acid, supports the growth of Ab and Bc; however, its potential signaling function remains to be explored. In natural environments, iron limitation represents one of the primary challenges faced by microorganisms, particularly in alkaline soils and rhizosphere environments 16 . Our study reveals that siderophore sharing can facilitate interspecies cooperation in biofilm communities, contrasting with the conventional view of siderophores primarily being competitive tools. Such cooperative siderophore-mediated community performance challenges the traditional competition-focused models and suggests a mutualistic framework where iron acquisition becomes a community-wide benefit rather than an individual advantage. The ecological relevance of this siderophore-mediated cooperation was validated through PGP experiments, where the trispecies community significantly promoted the growth of Arabidopsis in iron-deficient soil. The use of bacillibactin-deficient mutants confirmed that cooperative siderophore production was essential for optimal PGP, directly linking laboratory observations of microbial cooperation to a tangible ecological outcome. These findings demonstrate how interspecies cooperation in biofilm communities can enhance both microbial survival and plant performance. The siderophore-mediated interaction extends beyond microbial communities to encompass cross-kingdom interactions, where bacterial siderophores can directly facilitate plant iron acquisition, ultimately benefiting both microbial survival and plant health in iron-limited conditions 28,29 . The regulatory network underlying the trispecies biofilm represents a highly integrated metabolic system where signaling molecules and nutritional exchanges contribute to community function. The simultaneous upregulation of siderophore biosynthesis and downregulation of BCAA synthesis pathways in Bv suggests a complementary expression pattern that demonstrates metabolic coupling between amino acid metabolism and siderophore synthesis. From a biosynthetic perspective, bacillibactin synthesis requires a substantial of amino acid precursors and energy investment, particularly because its backbone structure formation depends on amino acid-derived building blocks 50,51 . The downregulation of BCAA synthesis pathways accompanied by upregulation of protein synthesis-related genes suggests that Bv optimizes siderophore production through reallocation of amino acid resources 20 . Formation of the trispecies biofilm induces numerous unique gene expression changes that could not be predicted from simple dual-species interactions, suggesting the existence of higher-order regulatory mechanisms in multispecies environments, such as cross-regulation of quorum sensing signals and system-level redistribution of metabolic flux 52,53 . The AHL-mediated interspecies signaling mechanism exemplifies the sophisticated regulatory characteristics of this synergistic system, where Bc-derived C4-HSL and C6-HSL induce expression of the dhbA gene of Bv. Moreover, although Ab does not directly induce dhbA expression in Bv, it enhances the impact of Bc on Bv, suggesting that Ab may exert regulatory control by modulating signal molecule stability or the sensitivity of Bv to these interspecies signals. The decreased bacillibactin production in Bv in the absence of BCAAs further confirms the biochemical coupling between amino acid metabolism and siderophore synthesis, where BCAAs influence siderophore synthesis through multiple pathways: serving as direct substrates for peptide synthesis 54 , functioning as energy sources to supply ATP and reducing power 15 , and acting as metabolic regulatory signaling molecules 13 . These integrated mechanisms create a mutualistic cycle that enhances biofilm formation, increases community biomass, and promotes plant growth in iron-deficient conditions. Beyond the rhizosphere environment, internal metabolic patterns and emergent properties of multispecies biofilm communities are relevant in microbial coexistence across diverse habitats, drug resistance, and the spread of antibiotic resistance genes. Potential human pathogens often reside within multispecies biofilm communities. Therefore, a deeper understanding of the metabolic interactions and key signaling pathways within multispecies biofilms could motivate strategies for efficient disruption of biofilm communities and precise elimination of pathogenic bacteria. Future research could explore the influence of quorum sensing signals (such as AHLs) on siderophore production in microbiomes, as well as the universality of cross-kingdom signal regulation in multispecies biofilm communities. Materials and Methods Strains and growth conditions Bacterial strains and plasmids used in this study are listed in Table S1. Bacillus velezensis SQR9 (China General Microbiology Culture Collection Center, CGMCC no. 5808, NCBI accession no. PRJNA227504), Acinetobacter baumannii XL380 (Agricultural Culture Collection of China, ACCC61689, NCBI accession no. PRJNA593376), and Burkholderia contaminans XL73 (ACCC61690, NCBI accession no. PRJNA593683) were isolated from the cucumber rhizosphere 38 . Starting inoculum was prepared from overnight cultures grown in TSB (Hopebio HB4114; 30°C, 180 rpm), harvested by centrifugation (6,000 x g for 2 min), and resuspended in 0.9% NaCl solution to achieve an optical density at 600 nm of 1.0 (OD 600 = 1.0). When required, medium was supplemented with zeocin (20 μg L −1 ) or chloramphenicol (5 μg L −1 ). M9 medium was prepared using minimal salts (5´, Sigma-Aldrich, M6030, Germany) supplemented with 2 mmol L −1 MgSO 4 and 0.1 mmol L −1 CaCl 2 . This basal medium was further amended with various sole carbon sources as detailed in Table S2. Water-soluble compounds were added at 10 mg mL −1 , whereas compounds with limited solubility were added to their maximum solubility. Biofilm formation and quantification For visualization, colony biofilms were formed by spotting 2 μL of starting inoculum onto solid TSB (1.5% agar) and incubating statically at 30°C for 24 h. Pellicle biofilms were cultivated in 24-well plates by inoculating 20 μL of the starting inoculum into 2 mL of liquid TSB and incubating at 30°C for 24 h. For coculture biofilms, equal volumes of Ab, Bc, and Bv (set to same OD 600 ) starting inoculum were combined. Iron availability was modulated by supplementing TSB with either 0.1% 2DP stock solution (Macklin, B807242, China; 78.10 mg mL −1 ) to establish iron-restricted conditions, or 0.1% FeCl 3 ·6H 2 O stock solution (131.6 mg mL −1 ) to establish iron-rich conditions. Pellicle biofilm biomass was determined gravimetrically as described before 42 . Briefly, 100 μL of starting inoculum was added into 10 mL of TSB in six-well plates containing sterile nylon mesh cell strainers (100-μm pore size). After 24 h of static incubation at 30°C, the strainers were taken out, excess liquid was removed, and the total weight was recorded. Biomass was calculated as the difference between the total weight and the weight of the strainer. For strain-specific quantification of biofilm-associated cells, 100 μL of inoculum was added into 10 mL of TSB in six-well plates equipped with nylon mesh strainers containing sterilized 1.5-cm 2 woven filters. After 24 h of static incubation at 30°C, filters carrying pellicle biofilms were transferred to 1.5-mL tubes and stored at −80°C for DNA extraction. For planktonic growth analysis, 40 μL of inoculum was added to 4 mL of TSB in sterile tubes and cultured at 30°C, 180 rpm, for 24 h, followed by centrifugation and storage at −80°C. Strain-specific primers were designed through comparative genomic analysis, as described before 42 (see Table S3 for primer sequences and standard curves). Absolute quantification was performed in 20-μL reactions with ChamQ SYBR qPCR Master Mix (Vazyme, Q711, China) according to the manufacturer’s instructions. The thermal cycling conditions were: 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 45 s, with standard melt curve analysis. All treatments included six biological replicates. Whole-genome transcriptomic analysis and qRT–PCR Pellicle biofilms were harvested from woven filters as described above. RNA was extracted using the E.Z.N.A. Bacterial RNA Kit (Omega Bio-tek, USA). For transcriptomics, sequencing libraries were constructed with the NEBNext Ultra Directional RNA Library Prep Kit and sequenced on an Illumina HiSeq 2000 platform. Raw data have been deposited in the NCBI SRA database (BioProject accession number PRJNA1281721). Following quality filtration, reads were aligned to reference genomes using Bowtie2. Differential expression was analyzed with DESeq2 using false discovery rate (FDR) correction. Genes with log 2 fold-change > 2 and FDR < 0.05 were considered significantly differentially regulated. Functional categorization and pathway enrichment were performed with EggNOG-mapper v2 for KEGG Orthology terms. For qRT–PCR, RNA was reverse-transcribed into cDNA using the PrimeScript RT Reagent Kit. Expression profiles of selected genes [ dhbA , dhbF , btr , feuB , tasA , and epsD of Bv; pchA , pchG , pchF , pchE , and Ipr1 of Bc; as well as ilvC and ilvH of all three species (Ab, Bc, and Bv)] were determined by qRT–PCR using ChamQ SYBR qPCR Master Mix. The recA gene served as an internal control for Bv and Bc, whereas gyrB was used for Ab. Primer sequences are provided in Table S3. The PCR conditions were: 95°C for 10 s, followed by 40 cycles of 95°C for 5 s and 60°C for 34 s. Relative expression levels were calculated using the 2 −ΔΔCT method. Microscopy Bacterial cell morphology was examined using transmission electron microscopy. Overnight cultures were harvested, washed with phosphate-buffered saline, adjusted to OD 600 = 0.5, applied to formvar/carbon-coated copper grids, and negatively stained with 2% (w/v) uranyl acetate. Imaging was on a Hitachi HT7800 TEM (Hitachi High-Tech, Japan) Biofilm architecture was analyzed by SEM. Samples were fixed in 2.5% glutaraldehyde and 2% osmium tetroxide, dehydrated through graded ethanol and tert -butanol, and sputter-coated, before imaging on an SEM Regulus 8100 (Hitachi High-Tech). Fluorescent reporter expression in colony and pellicle biofilms was visualized with an Axio Zoom V16 stereomicroscope (Carl Zeiss, Jena, Germany). Green-fluorescent protein (GFP) signals were captured using excitation at 470/440 nm and emission at 525/550 nm. Exposure time was 50–100 ms with 15% LED intensity, and acquisition parameters were kept constant across all samples to allow comparison. Fluorescently labeled Bv strain construction Fluorescent transcriptional reporters were created as described by Xu et al . 8 Briefly, promoter regions of target genes were fused to a gfp coding sequence via overlap PCR and cloned into vector pNW33n. The resulting constructs (primer sequences given in Table S3) were introduced into WT Bv and derived mutants. Growth curve assays Bacterial planktonic growth was monitored in 96-well plates using an Agilent Synergy H1 microplate reader (BioTek, USA). Cultures were initiated by putting 2 μL of the starting inoculum into 200 μL of TSB or M9 medium, then incubated at 30°C with continuous shaking. OD 600 was measured every 30 min for up to 96 h. For carbon source use tests, M9 medium was supplemented with individual sole carbon sources (Table S2). For BCAA induction experiments, TSB was supplemented with leucine, isoleucine, valine, or a mixture (100 μg mL −1 ). Both OD 600 and GFP fluorescence values (excitation/emission, 485/528 nm) were recorded every 30 min for 48 h. All treatments included five replicates. MALDI–TOF imaging mass spectrometry (MS) Colony biofilms were incubated as described above. Spatial distribution of metabolites within colonies was analyzed using MALDI–TOF MS. Colonies were excised, transferred to glass slides, desiccated at 30°C for 2 h, and coated with 20 mg mL −1 2,5-dihydroxybenzoic acid containing 1.0% trifluoroacetic acid. The prepared slides were analyzed on an UltraFlextreme MALDI TOF/TOF system (Bruker Daltonics, USA) in reflector negative mode, m/z 100–2000. Data were processed with SCiLS Lab 2014b software with total ion count normalization. Target secondary metabolites (bacillibactin, bacillomycin D, difficidin, macrolactin A, and surfactin in Bv) were analyzed based on specific m/z values. HPLC/MS analysis Standard bacillibactin was obtained from Biophore Research Products (Universität Tübingen, Germany). Biofilm extracellular extracts were prepared by sonication of pellicles in water, followed by centrifugation and ethyl acetate extraction. Dried extracts were reconstituted in methanol and analyzed using an Agilent 6410B Triple Quadrupole LC/MS instrument coupled to an Agilent 1200 HPLC system, equipped with a reverse-phase column (XBridge C18 5 μm, 4.6 × 250 mm). The samples were resolved at 220 nm and 30°C, at a flow-rate of 0.3 mL min −1 , with a linear gradient from 70% to 100% (0–20 min) acetonitrile/water solvent (0.1% formic acid), and at 70% acetonitrile from 20–30 min 55 . Mass spectra were collected in negative electrospray ionization mode. For AHL analysis, Bc culture supernatants were extracted with ethyl acetate, dried, and dissolved in methanol. Samples were analyzed on a Waters e2695 equipped with an XBridge C18 column (5 μm, 4.6 × 250 mm) with isocratic elution (methanol/water 65:35, 1 mL min − 1 , 30°C, 40 min). Cross-feeding assay and metabolome analysis Individual strains were cultured in 100 mL of M9 medium with glucose as the sole carbon source at 30°C and 180 rpm, until glucose depletion by Bc and Bv. However, Ab could not consume all the glucose. Spent media were obtained by centrifugation and filtration, and used as growth substrates for reciprocal cultivation by inoculating 1% (v/v) of each strain into spent media from the other species. Growth was monitored by OD 600 for up to 72 h. Extracellular metabolites were analyzed by untargeted metabolomics using a Vanquish UHPLC system coupled to an Orbitrap Q Exactive HF-X mass spectrometer (Thermo Fisher Scientific, USA). Samples were separated on a Hyperil Gold C18 column (100 × 2.1 mm, 1.9 μm) with a 16-min linear gradient of water and methanol containing either 0.1% formic acid (positive mode) or 5 mmol L −1 ammonium acetate, pH 9.0 (negative mode). The flow-rate was 0.2 mL min −1 at 30°C. Mass spectra were acquired in both positive and negative electrospray ionization modes. Raw data were processed with Compound Discoverer 3.0 (Thermo Fisher Scientific) for peak alignment, feature extraction, and intensity normalization. Metabolites were annotated by matching to the mzCloud and ChemSpider databases, and further analyzed on the Majorbio Cloud Platform. Features with relative standard deviation >30% in quality control samples were excluded. Differential metabolites were identified by orthogonal partial least squares discriminant analysis combined with Student’s t -test (VIP ≥ 1, p < 0.05). Heatmaps of the top 30 differential metabolites were generated. Four biological replicates were used for each treatment. Analysis of PGP traits Phosphate solubilization activity was examined on specific agar medium supplemented with either calcium phytate or Ca 3 (PO 4 ) 2 as insoluble phosphate sources, following an established protocol 56 . Exopolysaccharide content within biofilm was determined using the phenol–sulfuric acid colorimetric assay 57 . Ammonia production was evaluated by assay with Nessler’s reagent in peptone medium 58 . Siderophore production was quantified using the Chrome Azurol S colorimetric assay with cell-free supernatants from cultures grown in modified King’s B medium 33 . All assays were performed with three to six biological replicates, and absorbance was measured at appropriate wavelengths according to standard methods. Plant pot experiment design Arabidopsis thaliana (Col0) seeds were surface-sterilized and germinated on Murashige and Skoog medium with 1% agar (Hopebio, HB8469). Seven-day-old uniform seedlings were transplanted into pots containing 250 g of a sterilized growth substrate composed of vermiculite and perlite (5:1, v/v). Four seedlings were planted per pot, with two replicates per treatment. The starting inoculum was mixed into the substrate to achieve a final density of 10 7 colony-forming units g −1 substrate. Iron deficiency was induced by adjusting ¼ Murashige and Skoog liquid medium to pH 8.0 with KOH before watering. Plants were grown with a long-day photoperiod (16-h light/8-h dark) at 22°C and watered every 4 days. After 4 weeks, plant weight was measured. Uninoculated plants were set as the control. The inoculation treatments were as follows: Ab, Bc, Bv- Δdhb , Bv, AbBc, Ab- Δdhb , AbBv, Bc- Δdhb , BcBv, AbBc- Δdhb , and AbBcBv. Statistical analysis Statistical analyses were performed using R version 4.4.1 or GraphPad Prism 10. Figures were generated using the ggplot2 R package, GraphPad Prism 10, and Adobe Illustrator 2025. Details of specific statistical tests, significance thresholds, and sample sizes are provided in the figure captions. Declarations Data availability Genome data for B. velezensis SQR9, A. baumannii XL380 and B. contaminans XL73 are available under NABI BioProject accession number PRJNA227504, PRJNA593376 and PRJNA593683, respectively. Transcriptome raw data have been deposited in the NCBI SRA database (BioProject accession number PRJNA1281721). All other data generated and analyzed during this study are either available in figure legends or can be requested from the corresponding author. Acknowledgments This work was financially supported by Jiangsu Provincial Key Basic Research Projects (BK20253034), and National Natural Science Foundation of China (32361143785, 42477310 and 32400113). JX was supported by a China Scholarship Council fellowship during his stay in Leiden. ÁTK was funded by the European Union (ERC, MicroClock, 101166968). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the European Research Council Executive Agency. Neither the European Union nor the granting authority can be held responsible for them. MLS is supported by the Danish National Research Foundation (DNRF137). Author contributions J.X., X.S., Á.T.K., Z.X., and R.Z conceived the project. J.X., X.S., K.D., and H.Z. performed the experiments. V.H.T and M.L. performed transcriptome analysis. W.X., N.Z., Á.T.K., Q.S., and R.Z. contributed to experimental design and methodology. X.S. performed the metabolomic analysis. J.X., K.D., and H.Z. conducted biofilm formation assays and microscopy analysis. J.X. performed gene knockout experiments, MALDI-TOF imaging and plant growth experiments. J.X., X.S., W.X., and N.Z. performed data analysis and statistics. J.X., X.S., and Á.T.K wrote the manuscript and corrections from all authors. Competing interests The authors declare that they have no competing interests. References Flemming, H.-C. et al. Biofilms: an emergent form of bacterial life. Nat Rev Microbiol 14 , 563–575 (2016). Hall-Stoodley, L., Costerton, J. W. & Stoodley, P. Bacterial biofilms: from the Natural environment to infectious diseases. Nat Rev Microbiol 2 , 95–108 (2004). Røder, H. L., Olsen, N. M. C., Whiteley, M. & Burmølle, M. Unravelling interspecies interactions across heterogeneities in complex biofilm communities. Environ Microbiol 22 , 5–16 (2020). Sadiq, F. A. et al. Community-wide changes reflecting bacterial interspecific interactions in multispecies biofilms. Crit Rev Microbiol 47 , 338–358 (2021). Van Den Berg, N. I. et al. Ecological modelling approaches for predicting emergent properties in microbial communities. Nat Ecol Evol 6 , 855–865 (2022). Estrela, S. & Brown, S. P. Metabolic and demographic feedbacks shape the emergent spatial structure and function of microbial communities. PLoS Comput Biol 9 , e1003398 (2013). Yang, N. Emergent bacterial community properties induce enhanced drought tolerance in Arabidopsis . npj Biofilms and Microbiomes 7, 82 (2021). Xu, Z. et al. Enhanced control of plant wilt dsease by a xylose-inducible degQ gene engineered into Bacillus velezensis strain SQR9XYQ. Phytopathology 109 , 36–43 (2019). Zelezniak, A. et al. Metabolic dependencies drive species co-occurrence in diverse microbial communities. Proc Natl Acad Sci USA 112 , 6449–6454 (2015). Culp, E. J. & Goodman, A. L. Cross-feeding in the gut microbiome: Ecology and mechanisms. Cell Host Microbe 31 , 485–499 (2023). D’Souza, G. et al. Ecology and evolution of metabolic cross-feeding interactions in bacteria. Nat. Prod. Rep. 35 , 455–488 (2018). Pande, S. et al. Privatization of cooperative benefits stabilizes mutualistic cross-feeding interactions in spatially structured environments. ISME J 10 , 1413–1423 (2016). Lynch, C. J. & Adams, S. H. Branched-chain amino acids in metabolic signalling and insulin resistance. Nat Rev Endocrinol 10 , 723–736 (2014). Gulick, A. M. Nonribosomal peptide synthetase biosynthetic clusters of ESKAPE pathogens. Nat Prod Rep 34 , 981–1009 (2017). Brosnan, J. T. & Brosnan, M. E. Branched-chain amino acids: enzyme and substrate regulation. J Nutr 136 , 207S-211S (2006). Schäfer, M. et al. Metabolic interaction models recapitulate leaf microbiota ecology. Science 381 , eadf5121 (2023). Sokolovskaya, O. M., Shelton, A. N. & Taga, M. E. Sharing vitamins: Cobamides unveil microbial interactions. Science 369 , eaba0165 (2020). Ch’ng, J.-H. et al. Heme cross-feeding can augment Staphylococcus aureus and Enterococcus faecalis dual species biofilms. ISME J 16 , 2015–2026 (2022). Wang, W. Narrow-spectrum resource-utilizing bacteria drive the stability of synthetic communities through enhancing metabolic interactions. Nat Commun 16 , 6088 (2025.) Andrews, S. C., Robinson, A. K. & Rodríguez-Quiñones, F. Bacterial iron homeostasis. FEMS Microbiol Rev 27 , 215–237 (2003). Tangen, K. L., Jung, W. H., Sham, A. P., Lian, T. & Kronstad, J. W. The iron- and cAMP-regulated gene SIT1 influences ferrioxamine B utilization, melanization and cell wall structure in Cryptococcus neoformans. Microbiology 153 , 29–41 (2007). Caza, M. & Kronstad, J. Shared and distinct mechanisms of iron acquisition by bacterial and fungal pathogens of humans. Front. Cell. Infect. Microbiol. 3 , (2013). Neilands, J. B. Iron absorption and transport in microorganisms. Annu. Rev. Nutr. 1 , 27–46 (1981). Saleem, A. et al. Alkaline and acidic soil constraints on iron accumulation by Rice cultivars in relation to several physio-biochemical parameters. BMC Plant Biol 23 , 397 (2023). Neilands, J. B. Siderophores: structure and function of microbial iron transport compounds . J Biol Chem 270 , 26723–26726 (1995). Miethke, M. & Marahiel, M. A. Siderophore-based iron acquisition and pathogen control. Microbiol Mol Biol Rev 71 , 413–451 (2007). Wei, Z. et al. Microbial siderophores for One Health. Trends Microbiol 33 , 1277-1285 (2025). Liu, Y. et al. Plant commensal type VII secretion system causes iron leakage from roots to promote colonization. Nat Microbiol 8 , 1434–1449 (2023). Tan, T. et al. Siderophore-mediated iron enrichment in the biofilm matrix enhances plant iron nutrition. Cell Rep 44 , 116481 (2025). Xu, Z. et al. Antibiotic bacillomycin D affects iron acquisition and biofilm formation in bacillus velezensis through a btr-mediated FeuABC-dependent pathway. Cell Rep 29 , 1192-1202.e5 (2019). Rizzi, A. et al. Bacillus subtilis modulates its usage of biofilm-bound iron in response to environmental iron availability. Appl Environ Microbiol 86 , e00944-20 (2020). Qin, Y. et al. Heterogeneity in respiratory electron transfer and adaptive iron utilization in a bacterial biofilm. Nat Commun 10 , 3702 (2019). Gu, S. et al. Competition for iron drives phytopathogen control by natural rhizosphere microbiomes. Nat Microbiol 5 , 1002–1010 (2020). Harrison, F., Paul, J., Massey, R. C. & Buckling, A. Interspecific competition and siderophore-mediated cooperation in Pseudomonas aeruginosa . ISME J 2 , 49–55 (2008). Jenul, C. et al. Pyochelin biotransformation by Staphylococcus aureus shapes bacterial competition with Pseudomonas aeruginosa in polymicrobial infections. Cell Rep 42 , (2023). Butaitė, E., Baumgartner, M., Wyder, S. & Kümmerli, R. Siderophore cheating and cheating resistance shape competition for iron in soil and freshwater Pseudomonas communities. Nat Commun 8 , 414 (2017). D’Onofrio, A. et al. Siderophores from neighboring organisms promote the growth of uncultured bacteria. Chem Biol 17 , 254–264 (2010). Sun, X. et al. Metabolic interactions affect the biomass of synthetic bacterial biofilm communities. mSystems 8 , e01045-23 (2023). Wang, P. et al. An approach to constructing multispecies biofilm communities from rhizosphere soil. JoVE e66926 (2024) doi:10.3791/66926. Miethke, M. et al. Ferri-bacillibactin uptake and hydrolysis in Bacillus subtilis . Mol Microbiol 61 , 1413–1427 (2006). Lewenza S, Sokol PA. Regulation of ornibactin biosynthesis and N-Acyl-l-Homoserine lactone production by CepR in Burkholderia cepacia . J Bacteriol 183 , 7 (2001). Sun, X. et al. Bacillus velezensis stimulates resident rhizosphere Pseudomonas stutzeri for plant health through metabolic interactions. ISME J 16 , 774–787 (2022). Li, Z.-L., Wang, Y.-H., Chu, J., Zhuang, Y.-P. & Zhang, S.-L. Effect of branched-chain amino acids, valine, isoleucine and leucine on the biosythesis of bitespiramycin 4"-O-acylspiramycins. Braz. J. Microbiol. 40 , 734–746 (2009). Blomstrand, E., Eliasson, J., Karlsson, H. K. & Köhnke, R. Branched-chain amino acids activate key enzymes in protein synthesis after physical exercise. J Nutr 136 , 269S-273S (2006). Piovia-Scott, J. et al. Greater species richness of bacterial skin symbionts better suppresses the amphibian fungal pathogen Batrachochytrium dendrobatidis . Microb Ecol 74 , 217–226 (2017). Hu, J. et al. Introduction of probiotic bacterial consortia promotes plant growth via impacts on the resident rhizosphere microbiome. Proc Biol Sci 13 , 1960 (2021). Mee, M. T., Collins, J. J., Church, G. M. & Wang, H. H. Syntrophic exchange in synthetic microbial communities. Proc Natl Acad Sci USA. 111 , (2014). Mataigne, V., Vannier, N., Vandenkoornhuyse, P. & Hacquard, S. Multi-genome metabolic modeling predicts functional inter-dependencies in the Arabidopsis root microbiome. Microbiome 10 , 217 (2022). Quillin, S. J., Tran, P. & Prindle, A. Potential roles for gamma-aminobutyric acid signaling in bacterial communities. Bioelectricity 3 , 120–125 (2021). Challis, G. L. & Naismith, J. H. Structural aspects of non-ribosomal peptide biosynthesis. Curr Opin Struct Biol 14 , 748–756 (2004). May, J. J., Kessler, N., Marahiel, M. A. & Stubbs, M. T. Crystal structure of DhbE, an archetype for aryl acid activating domains of modular nonribosomal peptide synthetases. Proc Natl Acad Sci USA. 99 , 12120–12125 (2002). Momeni, B., Waite, A. J. & Shou, W. Spatial self-organization favors heterotypic cooperation over cheating. eLife 2 , e00960 (2013). Estrela, S., Morris, J. J. & Kerr, B. Private benefits and metabolic conflicts shape the emergence of microbial interdependencies. Environ Microbiol 18 , 1415–1427 (2016). Harper, A. E. & Miller, R. H. Branched-chain amino acid metabolism. Annu Rev Nutr 4 , 409–54 (1984). Dertz, E. A., Stintzi, A. & Raymond, K. N. Siderophore-mediated iron transport in Bacillus subtilis and Corynebacterium glutamicum . J Biol Inorg Chem 11 , 1087–1097 (2006). Ahmad, F., Ahmad, I. & Khan, M. S. Screening of free-living rhizospheric bacteria for their multiple plant growth promoting activities. Microbiol Res 163 , 173–181 (2008). Masuko, T. et al. Carbohydrate analysis by a phenol–sulfuric acid method in microplate format. Anall Biochem 339 , 69–72 (2005). Nautiyal, C. S. An efficient microbiological growth medium for screening phosphate solubilizing microorganisms. FEMS Microbiol Lett 170 , 265–270 (1999). Additional Declarations There is NO Competing Interest. Supplementary Files Supplementarytables.xlsx Supplementary tables Supplementaryfigures.pdf Supplementary figures Cite Share Download PDF Status: Under Review Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-8956555","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":604036947,"identity":"363b83b0-6024-43ab-8570-8ef521e33c5f","order_by":0,"name":"Ákos T Kovács","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABOUlEQVRIie3RPUvEMBgH8JTCTZWuGaT5Ck8peMpV/BqOKQe9pSeO3axLXQ66ZvMr5BZdUwLe0sGxg0O7dHK4SW64wVy509J6uArmD2nS8PzyQhDS0fmjMZK2M18QRefthEBgtnNCtdMBMA9kFCqC95UdYh0n1pn6HAgyvwr6xE5sYbBYEjsrPqoqxg65pkisbyfOGM9zaaG3PsFC7cIL6bJy/gy0wJ5bUpQzmHkX7IYq0gwOtiN1GhocnzzhIMXBklEkLZABLyPImRr3BNmTK24XzTfZgrw7RmBHlqkfcBSNWvKIFUEgKSgi1kPiSlMtVfhTVobj9i5g1Um+gJnLi3dFYECc1X1dL2J8mWWyqTex75CHaV5tthMCq8hb01g6w2dp36tz1Ncg6f72wQ8hmfi1RkdHR+d/5RO133f+cIPsQwAAAABJRU5ErkJggg==","orcid":"https://orcid.org/0000-0002-4465-1636","institution":"Leiden University","correspondingAuthor":true,"prefix":"","firstName":"Ákos","middleName":"T","lastName":"Kovács","suffix":""},{"id":604036948,"identity":"da81bd84-fd09-4945-a083-86e45f0ea73d","order_by":1,"name":"Jiyu Xie","email":"","orcid":"","institution":"Nanjing Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Jiyu","middleName":"","lastName":"Xie","suffix":""},{"id":604036949,"identity":"a74d4c8c-50e9-4b17-bf8a-3352c3f15ac5","order_by":2,"name":"Xinli Sun","email":"","orcid":"","institution":"Nanjing Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Xinli","middleName":"","lastName":"Sun","suffix":""},{"id":604036950,"identity":"c917e071-49ac-4af1-a195-1d7193c2dd04","order_by":3,"name":"Viktor Hesselberg-Thomsen","email":"","orcid":"","institution":"Technical University of Denmark","correspondingAuthor":false,"prefix":"","firstName":"Viktor","middleName":"","lastName":"Hesselberg-Thomsen","suffix":""},{"id":604036951,"identity":"dad0f13c-4e8a-416a-abb8-edcdf4c0c4c9","order_by":4,"name":"Kun Duan","email":"","orcid":"","institution":"Nanjing Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Kun","middleName":"","lastName":"Duan","suffix":""},{"id":604036952,"identity":"ff626edd-0b07-4bd2-8cb3-7f8d6a86f8dd","order_by":5,"name":"Huihui Zhang","email":"","orcid":"","institution":"Nanjing Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Huihui","middleName":"","lastName":"Zhang","suffix":""},{"id":604036953,"identity":"8df5ae0c-2c76-4429-813b-522f48938286","order_by":6,"name":"Weibing Xun","email":"","orcid":"https://orcid.org/0000-0003-0068-8514","institution":"Nanjing Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Weibing","middleName":"","lastName":"Xun","suffix":""},{"id":604036954,"identity":"77d4fe9e-c18f-4e72-a7d5-764dfa8f4c75","order_by":7,"name":"Nan Zhang","email":"","orcid":"https://orcid.org/0000-0001-8444-7456","institution":"Nanjing Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Nan","middleName":"","lastName":"Zhang","suffix":""},{"id":604036955,"identity":"2dbd548e-2d80-4f66-99c7-d9647623ceb0","order_by":8,"name":"Mikael Strube","email":"","orcid":"https://orcid.org/0000-0003-0905-5705","institution":"Technical University of Denmark","correspondingAuthor":false,"prefix":"","firstName":"Mikael","middleName":"","lastName":"Strube","suffix":""},{"id":604036956,"identity":"a01734aa-537b-4334-a90a-dc8aeeedabac","order_by":9,"name":"Qirong Shen","email":"","orcid":"https://orcid.org/0000-0002-5662-9620","institution":"Nanjing Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Qirong","middleName":"","lastName":"Shen","suffix":""},{"id":604036957,"identity":"bb871bb5-0c24-4d76-98c8-61451147731d","order_by":10,"name":"Ruifu Zhang","email":"","orcid":"https://orcid.org/0000-0002-3334-4286","institution":"Nanjing Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Ruifu","middleName":"","lastName":"Zhang","suffix":""},{"id":604036958,"identity":"996a5ba3-072e-4c7a-b7b0-4d7f0ee5008e","order_by":11,"name":"Zhihui Xu","email":"","orcid":"https://orcid.org/0000-0002-3987-8836","institution":"Nanjing Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Zhihui","middleName":"","lastName":"Xu","suffix":""}],"badges":[],"createdAt":"2026-02-24 11:05:14","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-8956555/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-8956555/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":104480687,"identity":"c3939702-bb1c-4a4e-bd3b-34fd235848a0","added_by":"auto","created_at":"2026-03-12 09:19:20","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":306238,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eSynergistic biofilm formation. (A) \u003c/strong\u003eBiofilm phenotype and biomass quantification. The well diameter was 15.6 mm.\u003cstrong\u003e \u003c/strong\u003eThe\u003cstrong\u003e \u003c/strong\u003eleft \u003cem\u003ey\u003c/em\u003e-axis indicates the fresh weight of biofilms. Data presented are the mean ± s.d. (\u003cem\u003en\u003c/em\u003e= 4). The right \u003cem\u003ey\u003c/em\u003e-axis displays the total cell numbers of biofilms as determined by quantitative PCR (qPCR). Data presented are the mean ± s.d. (\u003cem\u003en\u003c/em\u003e = 6). \u003cstrong\u003e(B) \u003c/strong\u003eScanning electron microscopy of synergic biofilms. Green: \u003cem\u003eAcinetobacter baumannii\u003c/em\u003e(Ab); blue: \u003cem\u003eBurkholderia contaminans\u003c/em\u003e (Bc); red: \u003cem\u003eBacillus velezensis\u003c/em\u003e (Bv)\u003cem\u003e. \u003c/em\u003e\u003cstrong\u003e(C\u003c/strong\u003e–\u003cstrong\u003eE) \u003c/strong\u003eCell numbers of Ab\u003cem\u003e \u003c/em\u003e\u003cstrong\u003e(C)\u003c/strong\u003e, Bc\u003cem\u003e\u003cstrong\u003e \u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e(D)\u003c/strong\u003e,\u003cstrong\u003e \u003c/strong\u003eand Bv\u003cem\u003e \u003c/em\u003e\u003cstrong\u003e(E)\u003c/strong\u003e in different treatments. Data presented are the mean ± s.d. (\u003cem\u003en\u003c/em\u003e = 6). Different letters indicate statistically significant (\u003cem\u003ep\u003c/em\u003e \u0026lt; 0.05) differences according to one-way analysis of variance (ANOVA) followed by Tukey’s honestly significant difference (HSD) \u003cem\u003epost-hoc\u003c/em\u003etest.\u003c/p\u003e","description":"","filename":"fig1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-8956555/v1/17cf05f52c5cf3a7dfb0a5fb.jpg"},{"id":104480686,"identity":"53d9780f-2355-45d0-a809-e4a50cfd712c","added_by":"auto","created_at":"2026-03-12 09:19:20","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":664983,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eTranscriptome response in multispecies biofilms. (A)\u003c/strong\u003e Numbers of differentially regulated genes in coculture compared with monoculture (log\u003csub\u003e2\u003c/sub\u003e fold-change \u0026gt; 2, \u003cem\u003ep\u003c/em\u003e \u0026lt; 0.05, \u003cem\u003et\u003c/em\u003e-test). The overlap of differentially upregulated \u003cstrong\u003e(B) \u003c/strong\u003eand downregulated \u003cstrong\u003e(C) \u003c/strong\u003egenes between different comparison groups. \u003cstrong\u003e(D) \u003c/strong\u003eDifferentially regulated Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in multispecies biofilms compared with monoculture. An asterisk (*) indicates that a pathway was significantly regulated (\u003cem\u003ep\u003c/em\u003e \u0026lt; 0.05, Fisher’s exact test, false discovery rate was adjusted by Benjamini–Hochberg correction). The pathway significance is color-coded, as shown (color represents the \u003cem\u003ep\u003c/em\u003e-value). \u003cstrong\u003e(E \u003c/strong\u003e\u0026amp;\u003cstrong\u003e F) \u003c/strong\u003eTranscriptome results for relative gene expression in branched-chain amino acid (BCAA) biosynthesis pathways \u003cstrong\u003e(E) \u003c/strong\u003eand siderophore biosynthesis pathways \u003cstrong\u003e(F) \u003c/strong\u003ein trispecies coculture compared with monoculture. Because the siderophore biosynthesis pathways of Ab and Bc have not been incorporated into the KEGG database, they have been plotted manually.\u003c/p\u003e","description":"","filename":"fig2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-8956555/v1/2104b524219358ff6c5b9929.jpg"},{"id":104480682,"identity":"b98debca-bcb3-4a99-b7cc-8ca764550502","added_by":"auto","created_at":"2026-03-12 09:19:20","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":362389,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eEnhanced bacillibactin production by Bv in interspecies interaction. (A) \u003c/strong\u003ePhenotype of trispecies biofilms formed by Ab, Bc, and wild-type (WT) or mutants of Bv in tryptic soy broth (TSB) with three different iron concentrations. FeCl\u003csub\u003e3\u003c/sub\u003e or 2-2‘-dipyridyl (2DP) was added at a final concentration of 0.5 mM to adjust the iron concentration. \u003cstrong\u003e(B) \u003c/strong\u003eTotal cell numbers in the biofilms, determined by qPCR. Data presented are the mean ± s.d. (\u003cem\u003en\u003c/em\u003e = 6). \u003cstrong\u003e(C) \u003c/strong\u003eInduction of the \u003cem\u003edhbABCEF\u003c/em\u003e operon promoter in Bv upon contact with Bc.\u003cstrong\u003e (D) \u003c/strong\u003eMatrix-assisted laser desorption/ionization time-of-flight mass (MALDI-TOF) spectrometry imaging of bacillibactin during interactions. \u003cstrong\u003e(E) \u003c/strong\u003eThe concentration of bacillibactin in biofilms. \u003cstrong\u003e(F) \u003c/strong\u003eGrowth curves of Bc\u003cem\u003e \u003c/em\u003ein different media. TSB with normal iron concentration was the control. TSB + 2DP was medium with a low iron concentration. The red line indicates medium containing 30% Bv spent\u003cem\u003e \u003c/em\u003eculture supernatant, which contains bacillibactin. The yellow line indicates medium containing 30% spent culture supernatant of Bv mutant ∆\u003cem\u003edhb\u003c/em\u003e, which does not produce bacillibactin. Data presented are the mean ± s.d. (\u003cem\u003en\u003c/em\u003e = 5). Different letters indicate statistically significant differences (\u003cem\u003ep\u003c/em\u003e \u0026lt; 0.05) according to one-way ANOVA followed by Tukey’s HSD \u003cem\u003epost-hoc\u003c/em\u003e test.\u003c/p\u003e","description":"","filename":"fig3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-8956555/v1/04a9440a38099c5679f3f311.jpg"},{"id":104480684,"identity":"ee8f5741-1fae-48fa-8937-105399b8a242","added_by":"auto","created_at":"2026-03-12 09:19:20","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":155374,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eTriculture promoted the growth of Bv BCAA-deficient mutants in biofilms. (A)\u003c/strong\u003e Phenotype of pellicles formed by BCAA-biosynthesis mutants of Bv in monoculture or in coculture with Ab and Bc. \u003cstrong\u003e(B) \u003c/strong\u003eCell numbers of Bv in monoculture and coculture. Data presented are the mean ± s.d. (\u003cem\u003en\u003c/em\u003e = 6). \u003cstrong\u003e(C) \u003c/strong\u003eThe concentration of bacillibactin in biofilms. Data presented are the mean ± s.d. (\u003cem\u003en\u003c/em\u003e = 3). Different letters indicate statistically significant (\u003cem\u003ep\u003c/em\u003e \u0026lt; 0.05) differences according to one-way ANOVA followed by Tukey’s HSD \u003cem\u003epost-hoc\u003c/em\u003e test. \u003cstrong\u003e(D) \u003c/strong\u003eExpression level of \u003cem\u003edhb \u003c/em\u003eof Bv in TSB supplemented with BCAAs as quantified by green-fluorescent protein (GFP) value (GFP/OD\u003csub\u003e600\u003c/sub\u003e). Data presented are the mean ± s.d. (\u003cem\u003en\u003c/em\u003e = 5). Statistical analysis was performed by \u003cem\u003et\u003c/em\u003e-test. (*\u003cem\u003ep\u003c/em\u003e \u0026lt; 0.05, **\u003cem\u003ep\u003c/em\u003e \u0026lt; 0.01).\u003c/p\u003e","description":"","filename":"fig4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-8956555/v1/58292358683378a2794b3987.jpg"},{"id":104780622,"identity":"51d22ce8-06bc-4b3f-a0ab-c375521766ea","added_by":"auto","created_at":"2026-03-17 07:53:25","extension":"jpg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":395653,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eMetabolic facilitation stabilizes the trispecies cooperation. (A)\u003c/strong\u003eSchematic representation of the spent medium assay. \u003cstrong\u003e(B) \u003c/strong\u003eMetabolic\u003cstrong\u003e \u003c/strong\u003eprofiles of spent media. “Bv” indicates the spent medium of Bv grown on M9 glucose medium; “Bv_Ab” indicates the spent medium of Ab grown on spent medium of Bv; and “Bv_Bc” indicates the spent medium of Bc grown on spent medium of Bv. An orange cycle indicates a substance that can be used by Bc,\u003cem\u003e \u003c/em\u003eas verified by growth assay. A green square indicates a substance that can be used by Ab,\u003cem\u003e \u003c/em\u003eas verified by growth assay.\u003c/p\u003e","description":"","filename":"fig5.jpg","url":"https://assets-eu.researchsquare.com/files/rs-8956555/v1/ec41a085d0c4ea956eb4eae9.jpg"},{"id":104480690,"identity":"f5474a38-4564-41a3-96b5-514d679d1af1","added_by":"auto","created_at":"2026-03-12 09:19:20","extension":"jpg","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":715311,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eMultispecies community enhanced plant growth in low-iron conditions. \u003c/strong\u003eThe\u003cstrong\u003e \u003c/strong\u003ephosphate [Ca\u003csub\u003e3\u003c/sub\u003e(PO\u003csub\u003e4\u003c/sub\u003e)\u003csub\u003e2 \u003c/sub\u003eand calcium phytate] solubilization \u003cstrong\u003e(A \u003c/strong\u003e\u0026amp;\u003cstrong\u003e B)\u003c/strong\u003e, ammonia production \u003cstrong\u003e(C \u003c/strong\u003e\u0026amp;\u003cstrong\u003e D)\u003c/strong\u003e,\u003cstrong\u003e \u003c/strong\u003esiderophore production \u003cstrong\u003e(E)\u003c/strong\u003e,\u003cstrong\u003e \u003c/strong\u003eand exopolysaccharide production \u003cstrong\u003e(F)\u003c/strong\u003e in different treatments. The plate diameter was 9 cm. Data presented are the mean ± s.d. (\u003cem\u003en\u003c/em\u003e = 5–9). \u003cstrong\u003e(G) \u003c/strong\u003eThe growth of \u003cem\u003eArabidopsis thaliana \u003c/em\u003ein an iron-deficient environment with different microbial treatments. \u003cstrong\u003e(H) \u003c/strong\u003eThe fresh weight of \u003cem\u003eA. thaliana\u003c/em\u003e. Data presented are the mean ± s.d. (\u003cem\u003en\u003c/em\u003e = 8). Different letters indicate statistically significant (\u003cem\u003ep\u003c/em\u003e \u0026lt; 0.05) differences according to one-way ANOVA followed by Tukey’s HSD \u003cem\u003epost-hoc\u003c/em\u003e test. CTL: control, not inoculated with bacteria; Bv-Δ\u003cem\u003edhb\u003c/em\u003e: a mutant of Bv\u003cem\u003e \u003c/em\u003ethat cannot produce bacillibactin.\u003c/p\u003e","description":"","filename":"fig6.jpg","url":"https://assets-eu.researchsquare.com/files/rs-8956555/v1/c5aa7e48be692ae01bcb3e60.jpg"},{"id":104781293,"identity":"3e981b66-769f-43fd-ab12-30991b03f964","added_by":"auto","created_at":"2026-03-17 07:55:20","extension":"jpg","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":369387,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eSummary figure. (A) \u003c/strong\u003eMetabolic cross-feeding and iron mobilization mechanism. Mutual metabolic exchange occurs between Bc and Bv. Bc metabolites, including BCAAs (valine, leucine, isoleucine), C4-HSL, and other substances, induced and improved bacillibactin biosynthesis in Bv. Bacillibactin chelates Fe\u003csup\u003e3+\u003c/sup\u003e, and FeuB transports iron complexes. FRO2 reduces Fe\u003csup\u003e3+\u003c/sup\u003e to Fe\u003csup\u003e2+\u003c/sup\u003e for plant uptake via IRT1 in \u003cem\u003eArabidopsis\u003c/em\u003e roots. \u003cstrong\u003e(B)\u003c/strong\u003e AbBcBv consortium inoculation enhanced plant uptake of iron, ammonia and phosphate, promoting the growth of \u003cem\u003eA. thaliana\u003c/em\u003e in iron-deficient soil.\u003c/p\u003e","description":"","filename":"fig7.jpg","url":"https://assets-eu.researchsquare.com/files/rs-8956555/v1/3d09e2d34be961d3071d9ea7.jpg"},{"id":104784455,"identity":"7f2330a0-2c0e-48a9-bd32-7efaaa97bc4d","added_by":"auto","created_at":"2026-03-17 08:07:43","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":4195956,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8956555/v1/3e8873a4-8091-49e1-8a6d-a970a60d9d8f.pdf"},{"id":104780945,"identity":"3295a21f-be40-4b15-946e-8da577541e89","added_by":"auto","created_at":"2026-03-17 07:54:20","extension":"xlsx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":18000,"visible":true,"origin":"","legend":"Supplementary tables","description":"","filename":"Supplementarytables.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-8956555/v1/8cedb54f35ea3e071b33ae1b.xlsx"},{"id":104480689,"identity":"8f7aec38-9799-46fd-9004-6122d5b45d71","added_by":"auto","created_at":"2026-03-12 09:19:20","extension":"pdf","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":3736897,"visible":true,"origin":"","legend":"Supplementary figures","description":"","filename":"Supplementaryfigures.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8956555/v1/b0a06841c7932fe954e30b1b.pdf"}],"financialInterests":"There is \u003cb\u003eNO\u003c/b\u003e Competing Interest.","formattedTitle":"Metabolic exchange and siderophore sharing underlie emergent biofilm synergism","fulltext":[{"header":"Introduction ","content":"\u003cp\u003eBiofilms are complex aggregates of microorganisms attached to surfaces or to each other, composed of microbial cells embedded in extracellular polymeric substances. Biofilms are one of the primary modes of microbial existence in natural environments, including soil\u0026nbsp;\u003csup\u003e1,2\u003c/sup\u003e. This lifestyle provides numerous advantages to microorganisms, including improved nutrient acquisition, enhanced resistance to environmental stresses, and facilitation of horizontal gene transfer. In nature, biofilms typically contain multiple species; they have emergent community properties that cannot be predicted from the individual species alone and exceed the sum of contributions of individual member species of the community\u0026nbsp;\u003csup\u003e3,4\u003c/sup\u003e. Understanding and ultimately predicting emergent properties of biofilm communities remains a key challenge in microbial ecology\u0026nbsp;\u003csup\u003e5\u003c/sup\u003e. Although it is evident that emergent properties arise from interactions within the community, how these emergent features develop in biofilm communities and how they contribute to the advancement of biofilm populations or even influence their hosts remain poorly understood\u0026nbsp;\u003csup\u003e1,6\u003c/sup\u003e.\u0026nbsp;For example, multispecies biofilms in the rhizosphere can provide plants with enhanced drought tolerance and pathogen suppression\u0026nbsp;\u003csup\u003e7,8\u003c/sup\u003e. These experimental observations suggest that rhizosphere biofilm communities can influence plant growth and health, yet a definitive link between biofilm community emergent properties and plant performance is yet to be established.\u003c/p\u003e\n\u003cp\u003eMetabolic cross-feeding is a widespread interaction in microbial communities and is particularly important in resource-limited environments\u0026nbsp;\u003csup\u003e9,10\u003c/sup\u003e. In multispecies biofilms, spatial proximity creates favorable conditions for metabolic cross-feeding, whereas the matrix-embedded structure minimizes diffusion losses and decreases the risk of metabolite exploitation by free-living microorganisms\u0026nbsp;\u003csup\u003e11\u003c/sup\u003e. Amino acid cross-feeding, particularly of branched-chain amino acids (BCAAs), is a common form of metabolic cross-feeding among bacteria\u0026nbsp;\u003csup\u003e12\u003c/sup\u003e.\u0026nbsp;BCAAs not only serve as nutritional resources but also act as regulatory molecules of nonribosomal peptide synthesis\u0026nbsp;\u003csup\u003e13\u0026ndash;15\u003c/sup\u003e. Other compounds, including vitamins, heme, organic acids, and sugars, have been implicated in microbial cross-feeding\u0026nbsp;\u003csup\u003e16\u0026ndash;19\u003c/sup\u003e. Metabolite cross-feeding underpins the stability of biofilm communities. However, their rapid development and expansion may necessitate the management of other limiting resources.\u003c/p\u003e\n\u003cp\u003eIron is an essential element for nearly all microbial life activities\u0026nbsp;\u003csup\u003e20\u003c/sup\u003e, participating in important processes such as electron transfer, DNA synthesis, and cellular respiration\u0026nbsp;\u003csup\u003e21,22\u003c/sup\u003e. However, in most natural environments, iron exits primarily as insoluble Fe (III) oxides with extremely low bioavailability\u0026nbsp;\u003csup\u003e23\u003c/sup\u003e. This limitation is particularly problematic in agricultural soils, especially in alkaline conditions, where iron deficiency significantly constrains both microbial activity and plant growth\u0026nbsp;\u003csup\u003e24\u003c/sup\u003e. To address this challenge, microorganisms produce and secrete high-affinity iron chelators called siderophores as a widespread, effective iron acquisition strategy\u0026nbsp;\u003csup\u003e25,26\u003c/sup\u003e. In the rhizosphere, microbial siderophores not only support bacterial iron acquisition but also improve plant iron nutrition\u0026nbsp;\u003csup\u003e27,28\u003c/sup\u003e. For example, bacillibactin, the siderophore produced by \u003cem\u003eBacillus\u003c/em\u003e, promotes plant iron acquisition in alkaline soil\u0026nbsp;\u003csup\u003e29\u003c/sup\u003e. Bacillibactin not only affects iron uptake but also influences the secretion of extracellular matrix and biofilm formation by \u003cem\u003eBacillus\u0026nbsp;\u003c/em\u003e\u003csup\u003e30\u003c/sup\u003e. For wild-type (WT) strains, the iron content in the extracellular matrix is 5-to-10-times higher than that inside the cells. In contrast, mutants with disrupted bacillibactin synthesis genes exhibit a significant decrease in extracellular matrix secretion capability and lose the ability to enrich iron in the extracellular matrix\u0026nbsp;\u003csup\u003e31,32\u003c/sup\u003e. Our latest research indicates that bacillibactin also acts as a signaling molecule that activates plant iron uptake, thereby enhancing iron absorption in iron-limited conditions\u0026nbsp;\u003csup\u003e29\u003c/sup\u003e.\u003c/p\u003e\n\u003cp\u003eTraditionally, siderophores have primarily been regarded as mediators of interspecies competition\u0026nbsp;\u003csup\u003e33\u003c/sup\u003e, where microorganisms compete for limited iron resources through the production of siderophores with different affinities, or by exploiting the siderophores produced by competitors. This competitive paradigm has been well-documented in pathogens, such as \u003cem\u003ePseudomonas aeruginosa\u003c/em\u003e and \u003cem\u003eStaphylococcus aureus\u003c/em\u003e, which engage in fierce competition for limited iron resources\u0026nbsp;\u003csup\u003e34,35\u003c/sup\u003e.\u0026nbsp;Similar competition occurs in soil, where siderophore-producing \u003cem\u003ePseudomonas\u003c/em\u003e can inhibit competitors through iron sequestration\u0026nbsp;\u003csup\u003e36\u003c/sup\u003e. Additionally, some microorganisms can directly use siderophores produced by other species, which has been observed in interactions between \u003cem\u003eP. aeruginosa\u003c/em\u003e and \u003cem\u003eMycobacterium\u003c/em\u003e species\u0026nbsp;\u003csup\u003e37\u003c/sup\u003e. The interspecies exploitation of siderophores raises intriguing questions: Can siderophore sharing promote cooperation among microorganisms? Additionally,\u0026nbsp;if siderophores function as a public good induced within biofilm communities, can these specialized metabolites also act as signals in the rhizosphere to modulate plant\u0026ndash;microbe interactions?\u003c/p\u003e\n\u003cp\u003eThis study explored the community-level emergent properties of a trispecies community of soil bacteria consisting of \u003cem\u003eAcinetobacter baumannii\u003c/em\u003e XL380 (Ab), \u003cem\u003eBurkholderia contaminans\u003c/em\u003e XL73 (Bc), and \u003cem\u003eBacillus velezensis\u003c/em\u003e SQR9 (Bv), that exhibits synergistic biofilm formation. We discovered a signaling and metabolite sharing network, including cross-feeding between Bv and the other two species through 5-aminovaleric acid (Bv to Bc/Ab) and BCAAs (Bc to Bv). Simultaneously, Bc secretes signaling molecules, including acyl-homoserine lactones (AHLs), that upregulate bacillibactin biosynthesis in Bv, which is exploited by Bc, explaining the enhanced biofilm community productivity. The synergistic interaction among the three species promotes plant growth performance in iron-deficient soil. We reveal a reciprocal mutualism that challenges traditional competitive interaction models in bacteria, and demonstrate how metabolic complementation can transform resource limitation into community-wide benefits.\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cstrong\u003eSynergistic biofilm formation by the trispecies community\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe have previously collected bacterial isolates from cucumber rhizosphere soil and performed static cultivation, which demonstrated enrichment of Ab, Bc, and Bv in pellicle biofilms\u0026nbsp;\u003csup\u003e38,39\u003c/sup\u003e. Biomass quantification revealed that a trispecies biofilm formed by Ab, Bc, and Bv (AbBcBv) had significantly higher fresh weight than any mono- or dual-species biofilms (Fig 1A). Multispecies biofilms exhibited more robust and wrinkled morphology than single-species pellicle biofilms.\u003c/p\u003e\n\u003cp\u003eTo determine whether the enhanced biofilm in the coculture resulted from mutualism or competition, we quantified the cell numbers of each species. In the trispecies biofilm, Bv was the most abundant species (~2.2×10\u003csup\u003e8\u003c/sup\u003e), followed by Bc (~7.5×10\u003csup\u003e7\u003c/sup\u003e), and Ab (~1.3×10\u003csup\u003e7\u003c/sup\u003e).\u0026nbsp;The observation was confirmed by scanning electron microscopy (SEM), which showed a low abundance of Ab closely associated with abundant Bc and Bv (Fig 1B, Fig S1). The trispecies combination (AbBcBv) demonstrated the highest total cell number among all groups, significantly exceeding those in any other culture condition. The cell numbers of Bc (Fig 1D) and Bv (Fig 1E) were significantly higher than in their respective monocultures, whereas Ab decreased compared with its monoculture (Fig 1C). Pairwise comparisons revealed variable interaction patterns. Both AbBv and BcBv showed significantly higher cell numbers than their respective monocultures, whereas AbBc showed no significant difference compared with Bc alone (Fig 1A, right axis, \u003cem\u003ep\u0026nbsp;\u003c/em\u003e\u0026lt; 0.05). These findings suggest the enhanced trispecies biofilm results from complex interactions involving both mutualism and competition, where Bc and Bv benefit from the multispecies environment whereas Ab, though decreased in amount, remains an integral member of the community.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTranscriptomic analysis of the trispecies biofilm\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo elucidate the molecular details governing synergistic interactions, global gene expression profiles were systematically investigated. Compared with their monocultures, Bc and Bv differentially expressed \u0026gt;1000 and hundreds of genes, respectively, in BcBv dual- and AbBcBv trispecies biofilms (Fig 2A). Lists of the differentially expressed genes in BcBv and AbBcBv biofilms showed substantial overlap, including 1643 genes in Bc and 248 genes in Bv (Fig 2B,C). In contrast, only a few genes were differentially expressed in AbBc and AbBv combinations compared with monocultures (Fig 2A). Therefore, Bc and Bv are the primary interacting species that drive major transcriptomic changes, whereas Ab exerts limited effects on the transcriptional landscapes of the other community members.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eFunctional analysis\u0026nbsp;of multispecies biofilms compared with monocultures\u0026nbsp;revealed that 24 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways\u0026nbsp;were differentially regulated (\u003cem\u003ep\u003c/em\u003e \u0026lt; 0.05) (Fig 2D). Several KEGG pathways of Bc were significantly upregulated both in BcBv and the trispecies biofilm, including xylene degradation, two-component systems, sulfur metabolism, and ABC transporters. Conversely, genes encoding the biosynthesis of antimicrobials were significantly downregulated.\u0026nbsp;In Bv, the gene cluster encoding siderophore biosynthesis was significantly upregulated in both BcBv and trispecies biofilms compared with monoculture. The genes encoding the ribosomes of Bv were significantly upregulated in all coculture groups, suggesting enhanced protein synthesis and faster growth, which is consistent with the higher cell number observed in mixed species biofilms. Conversely, genes encoding the biosynthesis of BCAAs (valine, leucine, and isoleucine) were significantly downregulated in BcBv and trispecies biofilms, suggesting decreased BCAAs production by Bv.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eGiven that Bv emerged as the predominant species in biofilm formation, we specifically focused on two pathways that were significantly regulated in multispecies biofilms compared with Bv monoculture: the upregulated siderophore pathway, and the downregulated BCAA synthesis pathway. Real-time quantitative -PCR (RTq-PCR) experiments confirmed that transcription of all genes for biosynthesis of bacillibactin (the siderophore of Bv) was upregulated in the triculture biofilm (Fig 2F, Fig S2A). In contrast, the transcription of all genes for biosynthesis of pyochelin and ornibactin (siderophores of Bc) was downregulated (Fig 2F, Fig S2B). Regarding the siderophores of Ab, 75% of the biosynthesis genes encoding acinetobactin synthesis were upregulated in the triculture biofilm, whereas the transcription of genes encoding the baoumannoferrin biosynthesis pathway remained unaffected. In the BCAA biosynthesis pathway, the \u003cem\u003eilvC\u003c/em\u003e and \u003cem\u003eilvH\u003c/em\u003e genes of Ab and Bc were upregulated in the trispecies biofilm but downregulated in Bv compared with monocultures (Fig 2E, Fig S2C). These expression patterns suggest that Bv may provide siderophore for the biofilm community, whereas Ab and Bc provide BCAAs in return, establishing a potential metabolic complementation relationship.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eBc induces bacillibactin production in Bv and uses it for iron acquisition\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo explore the importance of iron and siderophores in synergistic biofilm formation, we used three mutants of Bv\u0026nbsp;\u003csup\u003e30\u003c/sup\u003e: \u003cem\u003eΔdhb\u003c/em\u003e,which is unable to produce bacillibactin and its precursor 2,3-dihydroxybenzoate (DHBA); \u003cem\u003eΔdhbF\u003c/em\u003e, which produces DHBA but not bacillibactin; and \u003cem\u003eΔfeuB\u003c/em\u003e, which lacks Fe–bacillibactin transport into the cells\u0026nbsp;\u003csup\u003e40\u003c/sup\u003e (Fig S3). Using trispecies cocultures with either WT or mutant strains of Bv, we tested three different iron levels by adding extra FeCl\u003csub\u003e3\u003c/sub\u003e or the iron chelator 2-dipyridyl (2DP) to standard tryptic soy broth (TSB) culture medium (Fig 3A, Fig S3). In low-iron (2DP-supplemented) medium, Co-\u003cem\u003eΔdhb\u003c/em\u003e, Co-\u003cem\u003eΔdhbF\u003c/em\u003e, and Co-\u003cem\u003eΔfeuB\u003c/em\u003e exhibited severely impaired growth compared with Co-WT (Co- denotes coculture), including significantly decreased cell numbers (Fig 3A,B). This suggests that siderophore produced by Bv supported the growth of the biofilm community in iron-limited conditions. In high-iron (FeCl\u003csub\u003e3\u003c/sub\u003e-supplemented) medium, all cocultures showed similar biofilm phenotypes, because abundant iron was directly accessible without a requirement for siderophores (Fig 3A). In TSB (intermediate level of iron), the biofilms of Co-\u003cem\u003eΔdhb\u003c/em\u003e, Co-\u003cem\u003eΔdhbF\u003c/em\u003e, and Co-\u003cem\u003eΔfeuB\u003c/em\u003e appeared thinner than Co-WT, and they displayed decreased cell numbers (Fig 3A,B), confirming the necessity of Bv siderophore in enhancing trispecies biofilm formation. Biofilms of Co-WT formed in high-iron conditions had significantly lower cell numbers than those in normal iron conditions, suggesting that siderophores play an additional role in the interspecies interaction, beyond iron scavenging.\u003c/p\u003e\n\u003cp\u003eNext, we monitored \u003cem\u003edhbA\u0026nbsp;\u003c/em\u003eexpression in WT Bv when it was cocultured with Ab or Bc. As the distance between colonies on an agar plate was decreased, Bc induced \u003cem\u003edhbA\u003c/em\u003e expression in Bv, whereas a comparable effect was not observed when Bv was paired with Ab (Fig 3C). The combination of Ab and Bc induced higher expression than did Bc alone, indicating that although Ab did not directly induce \u003cem\u003edhbA\u003c/em\u003e expression in Bv, it enhanced the interaction between Bc and Bv. To validate Bc-induced bacillibactin production in Bv, matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI–TOF) imaging was used, which confirmed an increased level of bacillibactin in the Bv colony in the presence of Bc (Fig 3D). Whereas bacillibactin production was significantly enhanced, the spatial distribution of other Bv-produced secondary metabolites showed varying responses in the presence of Bc: bacillomycin D and surfactin showed moderate increases, whereas the polyketides difficidin and macrolactin A were not changed. Quantitative analysis of coculture biofilms revealed an elevated bacillibactin concentration in the trispecies biofilm (142.09 ± 17.59 ug/well) compared with Bv monoculture (22.44 ± 0.77 ug/well) (Fig 3E, Fig S4). These results confirmed that the presence of Bc increased bacillibactin production in Bv\u003cem\u003e.\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003eNext, we evaluated whether Bv siderophores can support Bc growth at different iron levels (Fig 3F). In low-iron medium (TSB + 2DP), Bc showed no growth when cultured alone or when supplemented with spent medium from culture of the Bv \u003cem\u003eΔdhb\u003c/em\u003e strain. However, spent medium from culture of WT Bv restored the growth of Bc to a similar level to that in TSB without iron depletion, demonstrating that Bc can use bacillibactin produced by Bv to overcome iron limitation.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eWe then hypothesized that Bc secretes signaling molecules that induce siderophore production in Bv. Previous studies have demonstrated that AHLs produced by \u003cem\u003eBurkholderia cepacia\u003c/em\u003e can regulate siderophore production\u0026nbsp;\u003csup\u003e41\u003c/sup\u003e. We thus tested induction of \u003cem\u003edhbA\u0026nbsp;\u003c/em\u003eexpression in Bv by eight AHL homologs (Fig S5). At 12 and 18 h, C4-HSL and C6-HSL significantly induced \u003cem\u003edhbA\u003c/em\u003e expression compared with the negative control, with C4-HSL showing induction levels comparable to that in the positive control (Bc supernatant) by 18 h (Fig S5A,B). High-performance liquid chromatography (HPLC) analysis confirmed the presence of C4-HSL in Bc supernatant (Fig S6).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eIn conclusion, a sophisticated chemical interplay is present within the biofilm community: Bc stimulates bacillibactin production in Bv, possibly through specific AHLs; Bv provides an iron-scavenging siderophore that supports Bc growth in iron-limited conditions; and Ab reinforces this mutualism, establishing a synergistic relationship critical for the enhanced trispecies biofilm formation.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAb and Bc provide BCAAs to support the growth of Bv\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo validate the hypothesis that Bv exploits BCAAs produced by Ab and Bc in the mixed species biofilm, we cocultured Bv mutants Δ\u003cem\u003eilvA\u003c/em\u003e, Δ\u003cem\u003eilvCH\u003c/em\u003e, Δ\u003cem\u003eilvD\u0026nbsp;\u003c/em\u003e(unable to synthesize all three BCAAs) and Δ\u003cem\u003eleuBCD\u003c/em\u003e (unable to synthesize leucine)\u003csup\u003e42\u003c/sup\u003e with Ab and Bc in TSB (Fig 4A). The biomass of the resulting coculture biofilms showed a clear increase compared with monoculture biofilms. Compared with the monocultures, the cell numbers of the WT, Δ\u003cem\u003eilvA\u003c/em\u003e,Δ\u003cem\u003eilvD\u003c/em\u003e,and Δ\u003cem\u003eilvCH\u0026nbsp;\u003c/em\u003estrains exhibited a significant increase, whereas strain Δ\u003cem\u003eleuBCD\u0026nbsp;\u003c/em\u003eshowed minimal growth enhancement (Fig 4B). These results demonstrated that Ab and Bc can provide BCAAs, particularly valine and isoleucine, to support the growth of Bv in the trispecies biofilm, confirming the existence of BCAA cross-feeding between these species.\u003c/p\u003e\n\u003cp\u003eBCAAs have been shown to influence secondary metabolite biosynthesis pathways\u0026nbsp;\u003csup\u003e43\u003c/sup\u003e, which is connected with nonribosomal peptide synthetases (NRPSs)\u0026nbsp;\u003csup\u003e44\u003c/sup\u003e. To examine whether BCAA availability affects bacillibactin production in Bv, bacillibactin levels were quantified (Fig 4C). In monocultures, the BCAA-deficient mutant Δ\u003cem\u003eilvA\u003c/em\u003e showed decreased bacillibactin production compared with the WT strain. Furthermore, direct supplementation of individual BCAAs significantly induced \u003cem\u003edhbA\u0026nbsp;\u003c/em\u003eexpression in Bv compared with negative controls, with the positive control (Bc spent medium) showing the highest level of induction (Fig 4D). When the BCAA-deficient mutant was cocultured with Ab and Bc in trispecies biofilms, bacillibactin production was significantly restored (Fig 4C). This restoration of bacillibactin production in the mutant coculture\u0026nbsp;demonstrates that BCAAs secreted by Ab and Bc can be used by Bv to support siderophore biosynthesis, Bv growth, or both at the same time.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eBv supplies metabolites that promote the growth of Bc and Ab\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo explore other potential cross-feeding interactions in the trispecies community, the influence of spent medium on growth was tested (Fig 5A)\u0026nbsp;\u003csup\u003e42\u003c/sup\u003e. During cultivation in M9 medium, Bc depletes glucose, whereas residual glucose remains in Ac culture supernatant. The culture supernatant of Bc did not support growth of the other species (Fig S7). In contrast, spent medium from culture of Bv supported the growth of both Ab and Bc, indicating that Bv provides metabolites to the other two species.\u003c/p\u003e\n\u003cp\u003eThe top 30 differential compounds were analyzed based on the relative abundance (Fig 5B). Highly abundant compounds in Bv spent medium that were depleted after cultivation of Ab or Bc were designated as cross-fed compounds. These included valeric acid, 5-aminovaleric acid, pentadecanoic acid, levulinic acid, cinnamic acid, valine, 2,3,4,5-tetramethylpyrazine, and tretinoin. Additionally, species-specific use patterns were observed: nine compounds (\u003cem\u003eN\u003c/em\u003e-α-L-acetyl-arginine, agmatine, \u003cem\u003eN\u003c/em\u003e-acetyl-DL-norvaline, \u003cem\u003eN\u003c/em\u003e-acetylvaline, valylproline, prolylleucine, \u003cem\u003eN\u003c/em\u003e-acetyl-D-alloisoleucine, nicotinic acid, and 2-isopropylmalic acid) were used exclusively by Bc, whereas two compounds (phenylethanolamine and synephrine) were consumed specifically by Ab. Testing the impact of identified compounds on bacterial growth revealed that Ab grew on synephrine, 5-aminovaleric acid, and pentadecanoic acid, whereas Bc used 5-aminovaleric acid and nicotinic acid (Fig S8).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003ePlant growth promotion in iron-deficient soil\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo evaluate the ecological significance of iron-mediated interaction of the trispecies biofilm community, plant-growth-promoting (PGP) traits of the community were tested. The trispecies community demonstrated superior PGP traits, including enhanced phosphate solubilization (Fig 6A,B), ammonia production (Fig 6C,D), siderophore production (Fig 6E), and exopolysaccharide production (Fig 6F), compared with mono- and dual-species treatments. Growth of \u003cem\u003eArabidopsis thaliana\u003c/em\u003e was severely stunted in the iron-deficient soil without inoculation (the control, CTL), confirming the detrimental influence of iron limitation on plant growth (Fig 6G,H). Plants inoculated with the trispecies consortium (AbBcBv) exhibited the most robust growth and highest fresh weight, significantly outperforming all other treatments (Fig 6G,H). Critically, when the bacillibactin-deficient mutant (Bv-Δ\u003cem\u003edhb\u003c/em\u003e) was used instead of WT Bv, plant growth was substantially decreased, demonstrating the critical role of bacillibactin in plant iron nutrition in iron-deficient conditions. Dual-species combinations showed intermediate growth promotion effects, with BcBv performing better than AbBv and AbBc. This plant growth promotion pattern correlated with the biofilm productivity, where BcBv demonstrated stronger synergistic interactions than the other dual-species combinations.\u0026nbsp;\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eHere, we reveal emergent properties of a trispecies bacterial community, including enhanced biofilm formation, metabolic cross-feeding, elevated siderophore production, and improved plant growth promotion (Fig 7). Within the community, Bv serves as the primary metabolic supplier of essential metabolites, such as 5-aminovaleric acid, pentadecanoic acid, nicotinic acid, and synephrine, that support the growth of Bc and Ab, whereas Bc and Ab reciprocally deliver BCAAs. Biofilm synergism is further promoted by elevated bacillibactin production by Bv that promotes iron acquisition by the entire community; the increased bacillibactin production is stimulated both by the signaling molecule C4-HSL from Bc and by BCAAs from Bc and Ab. In iron-deficient conditions, this metabolic complementation not only secures community survival but also confers plant growth promotion effects. Our findings demonstrate that siderophores facilitate interspecies synergism in biofilms and enhance plant iron nutrition.\u003c/p\u003e\n\u003cp\u003eThe enhancement of the trispecies biofilm observed in this study represents emergent properties arising from complex interspecies interactions. The significantly higher biomass and cell numbers in AbBcBv treatment compared with monoculture demonstrate such emergent benefits, where the community performance surpasses what would be expected from simply combining the capabilities of individual species. Emergent properties in diverse multispecies biofilms are well-documented, including enhanced pathogen suppression by skin bacterial symbionts and improved stress tolerance in soil microbial communities\u0026nbsp;\u003csup\u003e45\u003c/sup\u003e. Skin bacterial symbionts with greater species richness exhibit superior suppression of the amphibian fungal pathogen \u003cem\u003eBatrachochytrium dendrobatidis\u003c/em\u003e through combined dominant effects and complementarity mechanisms. Additionally, probiotic bacterial consortia consisting of up to eight \u003cem\u003ePseudomonas\u003c/em\u003e strains demonstrate emergent consortium-level effects in the tomato rhizosphere, where four- and eight- strain consortia reached densities up to 10-times higher than single-strain inoculants, exhibiting transgressive overyielding that could not be predicted from individual strain performance\u0026nbsp;\u003csup\u003e46\u003c/sup\u003e. These examples demonstrate how community-level emergent properties can provide ecological advantages that individual species cannot achieve alone. These emergent characteristics reflect a complex interplay of spatial organization, metabolic interdependencies, and chemical communication networks within structured microbial communities.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eMetabolic cross-feeding is a fundamental mechanism that drives microbial community assembly and stability. This phenomenon has been extensively documented across diverse systems, from vitamin interdependencies\u0026nbsp;\u003csup\u003e17\u003c/sup\u003e to amino acid exchange networks in synthetic communities\u0026nbsp;\u003csup\u003e47\u003c/sup\u003e. Metabolic interactions can augment biofilm formation in dual-species systems (e.g., heme cross-feeding in \u003cem\u003eStaphylococcus aureus\u003c/em\u003e and \u003cem\u003eEnterococcus faecalis\u003c/em\u003e biofilms), and multigenome metabolic modeling predicted functional interdependencies in plant microbiomes\u0026nbsp;\u003csup\u003e18,48\u003c/sup\u003e. Among the three dual-species treatments tested here, BcBv exhibited the strongest synergistic effects, a pattern consistently validated across transcriptomic and metabolomic analyses. Metabolome analysis revealed a biochemical basis underlying the mutualism between Bc and Bv, where Bv functions as the primary metabolite provider, secreting compounds including 5-aminovaleric acid, pentadecanoic acid, nicotinic acid, and synephrine that support community growth. Reciprocally, Ab and Bc provide BCAAs that promote the metabolic processes of Bv, particularly siderophore biosynthesis. BCAA cross-feeding has also been reported in other PGP rhizosphere biofilm communities\u0026nbsp;\u003csup\u003e7,42\u003c/sup\u003e suggesting that amino acid cross-feeding represents a widespread strategy in multispecies biofilms.\u0026nbsp;g-Aminobutyric acid (GABA) can function as a cross-kingdom signal, mediating communication between bacteria, influencing gene expression and behaviors, and thereby modulating the composition and dynamics of microbial communities\u0026nbsp;\u003csup\u003e49\u003c/sup\u003e. We revealed that a Bv-secreted GABA homolog, 5-aminovaleric acid, supports the growth of Ab and Bc; however, its potential signaling function remains to be explored.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eIn natural environments, iron limitation represents one of the primary challenges faced by microorganisms, particularly in alkaline soils and rhizosphere environments\u0026nbsp;\u003csup\u003e16\u003c/sup\u003e. Our study reveals that siderophore sharing can facilitate interspecies cooperation in biofilm communities, contrasting with the conventional view of siderophores primarily being competitive tools. Such cooperative siderophore-mediated community performance challenges the traditional competition-focused models and suggests a mutualistic framework where iron acquisition becomes a community-wide benefit rather than an individual advantage. The ecological relevance of this siderophore-mediated cooperation was validated through PGP experiments, where the trispecies community significantly promoted the growth of \u003cem\u003eArabidopsis\u003c/em\u003e in iron-deficient soil. The use of bacillibactin-deficient mutants confirmed that cooperative siderophore production was essential for optimal PGP, directly linking laboratory observations of microbial cooperation to a tangible ecological outcome. These findings demonstrate how interspecies cooperation in biofilm communities can enhance both microbial survival and plant performance. The siderophore-mediated interaction extends beyond microbial communities to encompass cross-kingdom interactions, where bacterial siderophores can directly facilitate plant iron acquisition, ultimately benefiting both microbial survival and plant health in iron-limited conditions\u0026nbsp;\u003csup\u003e28,29\u003c/sup\u003e.\u003c/p\u003e\n\u003cp\u003eThe regulatory network underlying the trispecies biofilm represents a highly integrated metabolic system where signaling molecules and nutritional exchanges contribute to community function. The simultaneous upregulation of siderophore biosynthesis and downregulation of BCAA synthesis pathways in Bv suggests a complementary expression pattern that demonstrates metabolic coupling between amino acid metabolism and siderophore synthesis. From a biosynthetic perspective, bacillibactin synthesis requires a substantial of amino acid precursors and energy investment, particularly because its backbone structure formation depends on amino acid-derived building blocks\u0026nbsp;\u003csup\u003e50,51\u003c/sup\u003e. The downregulation of BCAA synthesis pathways accompanied by upregulation of protein synthesis-related genes suggests that Bv optimizes siderophore production through reallocation of amino acid resources\u0026nbsp;\u003csup\u003e20\u003c/sup\u003e. Formation of the trispecies biofilm induces numerous unique gene expression changes that could not be predicted from simple dual-species interactions, suggesting the existence of higher-order regulatory mechanisms in multispecies environments, such as cross-regulation of quorum sensing signals and system-level redistribution of metabolic flux\u0026nbsp;\u003csup\u003e52,53\u003c/sup\u003e. The AHL-mediated interspecies signaling mechanism exemplifies the sophisticated regulatory characteristics of this synergistic system, where Bc-derived C4-HSL and C6-HSL induce expression of the \u003cem\u003edhbA\u003c/em\u003e gene of Bv. Moreover, although Ab does not directly induce \u003cem\u003edhbA\u0026nbsp;\u003c/em\u003eexpression in Bv, it enhances the impact of Bc on Bv, suggesting that Ab may exert regulatory control by modulating signal molecule stability or the sensitivity of Bv to these interspecies signals. The decreased bacillibactin production in Bv in the absence of BCAAs further confirms the biochemical coupling between amino acid metabolism and siderophore synthesis, where BCAAs influence siderophore synthesis through multiple pathways: serving as direct substrates for peptide synthesis\u0026nbsp;\u003csup\u003e54\u003c/sup\u003e, functioning as energy sources to supply ATP and reducing power\u0026nbsp;\u003csup\u003e15\u003c/sup\u003e, and acting as metabolic regulatory signaling molecules\u0026nbsp;\u003csup\u003e13\u003c/sup\u003e. These integrated mechanisms create a mutualistic cycle that enhances biofilm formation, increases community biomass, and promotes plant growth in iron-deficient conditions.\u003c/p\u003e\n\u003cp\u003eBeyond the rhizosphere environment, internal metabolic patterns and emergent properties of multispecies biofilm communities are relevant in microbial coexistence across diverse habitats, drug resistance, and the spread of antibiotic resistance genes. Potential human pathogens often reside within multispecies biofilm communities. Therefore, a deeper understanding of the metabolic interactions and key signaling pathways within multispecies biofilms could motivate strategies for efficient disruption of biofilm communities and precise elimination of pathogenic bacteria. Future research could explore the influence of quorum sensing signals (such as AHLs) on siderophore production in microbiomes, as well as the universality of cross-kingdom signal regulation in multispecies biofilm communities.\u003c/p\u003e"},{"header":"Materials and Methods ","content":"\u003cp\u003e\u003cstrong\u003eStrains and growth conditions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eBacterial strains and plasmids used in this study are listed in Table S1.\u0026nbsp;\u003cem\u003eBacillus\u0026nbsp;\u003c/em\u003e\u003cem\u003evelezensis\u003c/em\u003e SQR9 (China General Microbiology Culture Collection Center, CGMCC no. 5808, NCBI accession no. PRJNA227504),\u0026nbsp;\u003cem\u003eAcinetobacter\u0026nbsp;\u003c/em\u003e\u003cem\u003ebaumannii\u003c/em\u003e XL380 (Agricultural Culture Collection of China, ACCC61689, NCBI accession no. PRJNA593376), and\u0026nbsp;\u003cem\u003eBurkholderia\u003c/em\u003e\u003cem\u003e\u0026nbsp;contaminans\u003c/em\u003e XL73 (ACCC61690, NCBI accession no. PRJNA593683) were isolated from the cucumber rhizosphere\u0026nbsp;\u003csup\u003e38\u003c/sup\u003e. Starting inoculum was prepared from overnight cultures grown in TSB (Hopebio HB4114; 30°C, 180 rpm), harvested by centrifugation (6,000 x g for 2 min), and resuspended in 0.9% NaCl solution to achieve an optical density at 600 nm of 1.0 (OD\u003csub\u003e600\u003c/sub\u003e = 1.0). When required, medium was supplemented with zeocin (20 μg L\u003csup\u003e−1\u003c/sup\u003e) or chloramphenicol (5 μg L\u003csup\u003e−1\u003c/sup\u003e).\u003c/p\u003e\n\u003cp\u003eM9 medium was prepared using minimal salts (5´, Sigma-Aldrich, M6030, Germany) supplemented with 2 mmol L\u003csup\u003e−1\u003c/sup\u003e MgSO\u003csub\u003e4\u003c/sub\u003e and 0.1 mmol L\u003csup\u003e−1\u003c/sup\u003e CaCl\u003csub\u003e2\u003c/sub\u003e. This basal medium was further amended with various sole carbon sources as detailed in Table S2. Water-soluble compounds were added at 10 mg mL\u003csup\u003e−1\u003c/sup\u003e, whereas compounds with limited solubility were added to their maximum solubility.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eBiofilm formation and quantification\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFor visualization, colony biofilms were formed by spotting 2 μL of starting inoculum onto solid TSB (1.5% agar) and incubating statically at 30°C for 24 h. Pellicle biofilms were cultivated in 24-well plates by inoculating 20 μL of the starting inoculum into 2 mL of liquid TSB and incubating at 30°C for 24 h. For coculture biofilms, equal volumes of Ab, Bc, and Bv (set to same OD\u003csub\u003e600\u003c/sub\u003e) starting inoculum were combined. Iron availability was modulated by supplementing TSB with either 0.1% 2DP stock solution (Macklin, B807242, China; 78.10 mg mL\u003csup\u003e−1\u003c/sup\u003e) to establish iron-restricted conditions, or 0.1% FeCl\u003csub\u003e3\u003c/sub\u003e·6H\u003csub\u003e2\u003c/sub\u003eO stock solution (131.6 mg mL\u003csup\u003e−1\u003c/sup\u003e) to establish iron-rich conditions.\u003c/p\u003e\n\u003cp\u003ePellicle biofilm biomass was determined gravimetrically as described before\u0026nbsp;\u003csup\u003e42\u003c/sup\u003e. Briefly, 100 μL of starting inoculum was added into 10 mL of TSB in six-well plates containing sterile nylon mesh cell strainers (100-μm pore size). After 24 h of static incubation at 30°C, the strainers were taken out, excess liquid was removed, and the total weight was recorded. Biomass was calculated as the difference between the total weight and the weight of the strainer.\u003c/p\u003e\n\u003cp\u003eFor strain-specific quantification of biofilm-associated cells, 100 μL of inoculum was added into 10 mL of TSB in six-well plates equipped with nylon mesh strainers containing sterilized 1.5-cm\u003csup\u003e2\u003c/sup\u003e woven filters. After 24 h of static incubation at 30°C, filters carrying pellicle biofilms were transferred to 1.5-mL tubes and stored at −80°C for DNA extraction. For planktonic growth analysis, 40 μL of inoculum was added to 4 mL of TSB in sterile tubes and cultured at 30°C, 180 rpm, for 24 h, followed by centrifugation and storage at −80°C. Strain-specific primers were designed through comparative genomic analysis, as described before\u0026nbsp;\u003csup\u003e42\u003c/sup\u003e (see Table S3 for primer sequences and standard curves). Absolute quantification was performed in 20-μL reactions with ChamQ SYBR qPCR Master Mix (Vazyme, Q711, China) according to the manufacturer’s instructions. The thermal cycling conditions were: 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 45 s, with standard melt curve analysis. All treatments included six biological replicates.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eWhole-genome transcriptomic analysis and qRT–PCR\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003ePellicle biofilms were harvested from woven filters as described above. RNA was extracted using the E.Z.N.A. Bacterial RNA Kit (Omega Bio-tek, USA). For transcriptomics, sequencing libraries were constructed with the NEBNext Ultra Directional RNA Library Prep Kit and sequenced on an Illumina HiSeq 2000 platform. Raw data have been deposited in the NCBI SRA database (BioProject accession number PRJNA1281721). Following quality filtration, reads were aligned to reference genomes using Bowtie2. Differential expression was analyzed with DESeq2 using false discovery rate (FDR) correction. Genes with log\u003csub\u003e2\u003c/sub\u003e fold-change \u0026gt; 2 and FDR \u0026lt; 0.05 were considered significantly differentially regulated. Functional categorization and pathway enrichment were performed with EggNOG-mapper v2 for KEGG Orthology terms.\u003c/p\u003e\n\u003cp\u003eFor qRT–PCR, RNA was reverse-transcribed into cDNA using the PrimeScript RT Reagent Kit. Expression profiles of selected genes [\u003cem\u003edhbA\u003c/em\u003e, \u003cem\u003edhbF\u003c/em\u003e, \u003cem\u003ebtr\u003c/em\u003e, \u003cem\u003efeuB\u003c/em\u003e, \u003cem\u003etasA\u003c/em\u003e, and \u003cem\u003eepsD\u0026nbsp;\u003c/em\u003eof Bv; \u003cem\u003epchA\u003c/em\u003e, \u003cem\u003epchG\u003c/em\u003e, \u003cem\u003epchF\u003c/em\u003e, \u003cem\u003epchE\u003c/em\u003e, and \u003cem\u003eIpr1\u003c/em\u003e of Bc; as well as \u003cem\u003eilvC\u003c/em\u003e and \u003cem\u003eilvH\u003c/em\u003e of all three species (Ab, Bc, and Bv)] were determined by qRT–PCR using ChamQ SYBR qPCR Master Mix. The \u003cem\u003erecA\u003c/em\u003e gene served as an internal control for Bv and Bc, whereas \u003cem\u003egyrB\u003c/em\u003e was used for Ab. Primer sequences are provided in Table S3. The PCR conditions were: 95°C for 10 s, followed by 40 cycles of 95°C for 5 s and 60°C for 34 s. Relative expression levels were calculated using the 2\u003csup\u003e−ΔΔCT\u003c/sup\u003e method.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMicroscopy\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eBacterial cell morphology was examined using transmission electron microscopy. Overnight cultures were harvested, washed with phosphate-buffered saline, adjusted to OD\u003csub\u003e600\u003c/sub\u003e = 0.5, applied to formvar/carbon-coated copper grids, and negatively stained with 2% (w/v) uranyl acetate. Imaging was on a Hitachi HT7800 TEM (Hitachi High-Tech, Japan)\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eBiofilm architecture was analyzed by SEM. Samples were fixed in 2.5% glutaraldehyde and 2% osmium tetroxide, dehydrated through graded ethanol and \u003cem\u003etert\u003c/em\u003e-butanol, and sputter-coated, before imaging on an SEM Regulus 8100 (Hitachi High-Tech).\u003c/p\u003e\n\u003cp\u003eFluorescent reporter expression in colony and pellicle biofilms was visualized with an Axio Zoom V16 stereomicroscope (Carl Zeiss, Jena, Germany). Green-fluorescent protein (GFP) signals were captured using excitation at 470/440 nm and emission at 525/550 nm. Exposure time was 50–100 ms with 15% LED intensity, and acquisition parameters were kept constant across all samples to allow comparison.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFluorescently labeled Bv strain construction\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFluorescent transcriptional reporters were created as described by Xu \u003cem\u003eet al\u003c/em\u003e.\u0026nbsp;\u003csup\u003e8\u003c/sup\u003e Briefly, promoter regions of target genes were fused to a \u003cem\u003egfp\u003c/em\u003e coding sequence via overlap PCR and cloned into vector pNW33n. The resulting constructs (primer sequences given in Table S3) were introduced into WT Bv and derived mutants.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eGrowth curve assays\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eBacterial planktonic growth was monitored in 96-well plates using an Agilent Synergy H1 microplate reader (BioTek, USA). Cultures were initiated by putting 2 μL of the starting inoculum into 200 μL of TSB or M9 medium, then incubated at 30°C with continuous shaking. OD\u003csub\u003e600\u003c/sub\u003e was measured every 30 min for up to 96 h. For carbon source use tests, M9 medium was supplemented with individual sole carbon sources (Table S2).\u003c/p\u003e\n\u003cp\u003eFor BCAA induction experiments, TSB was supplemented with leucine, isoleucine, valine, or a mixture (100 μg mL\u003csup\u003e−1\u003c/sup\u003e). Both OD\u003csub\u003e600\u003c/sub\u003e and GFP fluorescence values (excitation/emission, 485/528 nm) were recorded every 30 min for 48 h. All treatments included five replicates.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMALDI–TOF imaging mass spectrometry (MS)\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eColony biofilms were incubated as described above. Spatial distribution of metabolites within colonies was analyzed using MALDI–TOF MS. Colonies were excised, transferred to glass slides, desiccated at 30°C for 2 h, and coated with 20 mg mL\u003csup\u003e−1\u003c/sup\u003e 2,5-dihydroxybenzoic acid containing 1.0% trifluoroacetic acid. The prepared slides were analyzed on an UltraFlextreme MALDI TOF/TOF system (Bruker Daltonics, USA) in reflector negative mode, \u003cem\u003em/z\u003c/em\u003e 100–2000. Data were processed with SCiLS Lab 2014b software with total ion count normalization. Target secondary metabolites (bacillibactin, bacillomycin D, difficidin, macrolactin A, and surfactin in Bv) were analyzed based on specific \u003cem\u003em/z\u003c/em\u003e values.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eHPLC/MS analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eStandard bacillibactin was obtained from Biophore Research Products (Universität Tübingen, Germany). Biofilm extracellular extracts were prepared by sonication of pellicles in water, followed by centrifugation and ethyl acetate extraction. Dried extracts were reconstituted in methanol and analyzed using an Agilent 6410B Triple Quadrupole LC/MS instrument coupled to an Agilent 1200 HPLC system, equipped with a reverse-phase column (XBridge C18 5 μm, 4.6 × 250 mm). The samples were resolved at 220 nm and 30°C, at a flow-rate of 0.3 mL min\u003csup\u003e−1\u003c/sup\u003e, with a linear gradient from 70% to 100% (0–20 min) acetonitrile/water solvent (0.1% formic acid), and at 70% acetonitrile from 20–30 min\u0026nbsp;\u003csup\u003e55\u003c/sup\u003e. Mass spectra were collected in negative electrospray ionization mode.\u003c/p\u003e\n\u003cp\u003eFor AHL analysis, Bc culture supernatants were extracted with ethyl acetate, dried, and dissolved in methanol. Samples were analyzed on a Waters e2695 equipped with an XBridge C18 column (5 μm, 4.6 × 250 mm) with isocratic elution (methanol/water 65:35, 1 mL min\u003csup\u003e−\u003c/sup\u003e\u003csup\u003e1\u003c/sup\u003e, 30°C, 40 min).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCross-feeding assay and metabolome analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIndividual strains were cultured in 100 mL of M9 medium with glucose as the sole carbon source at 30°C and 180 rpm, until glucose depletion by Bc and Bv. However, Ab could not consume all the glucose. Spent media were obtained by centrifugation and filtration, and used as growth substrates for reciprocal cultivation by inoculating 1% (v/v) of each strain into spent media from the other species. Growth was monitored by OD\u003csub\u003e600\u003c/sub\u003e for up to 72 h.\u003c/p\u003e\n\u003cp\u003eExtracellular metabolites were analyzed by untargeted metabolomics using a Vanquish UHPLC system coupled to an Orbitrap Q Exactive HF-X mass spectrometer (Thermo Fisher Scientific, USA). Samples were separated on a Hyperil Gold C18 column (100 × 2.1 mm, 1.9 μm) with a 16-min linear gradient of water and methanol containing either 0.1% formic acid (positive mode) or 5 mmol L\u003csup\u003e−1\u003c/sup\u003e ammonium acetate, pH 9.0 (negative mode). The flow-rate was 0.2 mL min\u003csup\u003e−1\u003c/sup\u003e at 30°C. Mass spectra were acquired in both positive and negative electrospray ionization modes. Raw data were processed with Compound Discoverer 3.0 (Thermo Fisher Scientific) for peak alignment, feature extraction, and intensity normalization. Metabolites were annotated by matching to the mzCloud and ChemSpider databases, and further analyzed on the Majorbio Cloud Platform. Features with relative standard deviation \u0026gt;30% in quality control samples were excluded. Differential metabolites were identified by orthogonal partial least squares discriminant analysis combined with Student’s \u003cem\u003et\u003c/em\u003e-test (VIP ≥ 1, p \u0026lt; 0.05). Heatmaps of the top 30 differential metabolites were generated. Four biological replicates were used for each treatment.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAnalysis of PGP traits\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003ePhosphate solubilization activity was examined on specific agar medium supplemented with either calcium phytate or Ca\u003csub\u003e3\u003c/sub\u003e(PO\u003csub\u003e4\u003c/sub\u003e)\u003csub\u003e2\u003c/sub\u003e as insoluble phosphate sources, following an established protocol\u0026nbsp;\u003csup\u003e56\u003c/sup\u003e. Exopolysaccharide content within biofilm was determined using the phenol–sulfuric acid colorimetric assay\u0026nbsp;\u003csup\u003e57\u003c/sup\u003e. Ammonia production was evaluated by assay with Nessler’s reagent in peptone medium\u0026nbsp;\u003csup\u003e58\u003c/sup\u003e. Siderophore production was quantified using the Chrome Azurol S colorimetric assay\u0026nbsp;with cell-free supernatants from cultures grown in modified King’s B medium\u0026nbsp;\u003csup\u003e33\u003c/sup\u003e. All assays were performed with three to six biological replicates, and absorbance was measured at appropriate wavelengths according to standard methods.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003ePlant pot experiment design\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eArabidopsis thaliana\u003c/em\u003e (Col0) seeds were surface-sterilized and germinated on Murashige and Skoog medium with 1% agar (Hopebio, HB8469). Seven-day-old uniform seedlings were transplanted into pots containing 250 g of a sterilized growth substrate composed of vermiculite and perlite (5:1, v/v). Four seedlings were planted per pot, with two replicates per treatment. The starting inoculum was mixed into the substrate to achieve a final density of 10\u003csup\u003e7\u003c/sup\u003e colony-forming units g\u003csup\u003e−1\u003c/sup\u003e substrate. Iron deficiency was induced by adjusting ¼ Murashige and Skoog liquid medium to pH 8.0 with KOH before watering. Plants were grown with a long-day photoperiod (16-h light/8-h dark) at 22°C and watered every 4 days. After 4 weeks, plant weight was measured. Uninoculated plants were set as the control. The inoculation treatments were as follows: Ab, Bc, Bv-\u003cem\u003eΔdhb\u003c/em\u003e, Bv, AbBc, Ab-\u003cem\u003eΔdhb\u003c/em\u003e, AbBv, Bc-\u003cem\u003eΔdhb\u003c/em\u003e, BcBv, AbBc-\u003cem\u003eΔdhb\u003c/em\u003e, and AbBcBv.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eStatistical analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eStatistical analyses were performed using R version 4.4.1 or GraphPad Prism 10. Figures were generated using the ggplot2 R package, GraphPad Prism 10, and Adobe Illustrator 2025. Details of specific statistical tests, significance thresholds, and sample sizes are provided in the figure captions.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eData availability\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eGenome data for \u003cem\u003eB. velezensis\u003c/em\u003e SQR9, \u003cem\u003eA. baumannii\u003c/em\u003e XL380 and \u003cem\u003eB. contaminans\u003c/em\u003e XL73 are available under NABI BioProject accession number PRJNA227504, PRJNA593376 and PRJNA593683, respectively. Transcriptome raw data have been deposited in the NCBI SRA database (BioProject accession number PRJNA1281721). All other data generated and analyzed during this study are either available in figure legends or can be requested from the corresponding author.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgments\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work was financially supported by Jiangsu Provincial Key Basic Research Projects (BK20253034), and National Natural Science Foundation of China (32361143785, 42477310 and 32400113). JX was supported by a China Scholarship Council fellowship during his stay in Leiden. \u0026Aacute;TK was funded by the European Union (ERC, MicroClock, 101166968). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the European Research Council Executive Agency. Neither the European Union nor the granting authority can be held responsible for them. MLS is supported by the Danish National Research Foundation (DNRF137).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor contributions\u003c/strong\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eJ.X., X.S., \u0026Aacute;.T.K., Z.X., and R.Z conceived the project. J.X., X.S., K.D., and H.Z. performed the experiments. V.H.T and M.L. performed transcriptome analysis. W.X., N.Z., \u0026Aacute;.T.K., Q.S., and R.Z. contributed to experimental design and methodology. X.S. performed the metabolomic analysis. J.X., K.D., and H.Z. conducted biofilm formation assays and microscopy analysis. J.X. performed gene knockout experiments, MALDI-TOF imaging and plant growth experiments. J.X., X.S., W.X., and N.Z. performed data analysis and statistics. J.X., X.S., and \u0026Aacute;.T.K wrote the manuscript and corrections from all authors.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no competing interests.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eFlemming, H.-C. \u003cem\u003eet al.\u003c/em\u003e Biofilms: an emergent form of bacterial life. \u003cem\u003eNat Rev Microbiol\u003c/em\u003e \u003cstrong\u003e14\u003c/strong\u003e, 563\u0026ndash;575 (2016).\u003c/li\u003e\n\u003cli\u003eHall-Stoodley, L., Costerton, J. W. \u0026amp; Stoodley, P. Bacterial biofilms: from the Natural environment to infectious diseases. \u003cem\u003eNat Rev Microbiol\u003c/em\u003e \u003cstrong\u003e2\u003c/strong\u003e, 95\u0026ndash;108 (2004).\u003c/li\u003e\n\u003cli\u003eR\u0026oslash;der, H. L., Olsen, N. M. C., Whiteley, M. \u0026amp; Burm\u0026oslash;lle, M. Unravelling interspecies interactions across heterogeneities in complex biofilm communities. \u003cem\u003eEnviron Microbiol\u003c/em\u003e \u003cstrong\u003e22\u003c/strong\u003e, 5\u0026ndash;16 (2020).\u003c/li\u003e\n\u003cli\u003eSadiq, F. A. \u003cem\u003eet al.\u003c/em\u003e Community-wide changes reflecting bacterial interspecific interactions in multispecies biofilms. \u003cem\u003eCrit Rev Microbiol\u003c/em\u003e \u003cstrong\u003e47\u003c/strong\u003e, 338\u0026ndash;358 (2021).\u003c/li\u003e\n\u003cli\u003eVan Den Berg, N. I. \u003cem\u003eet al.\u003c/em\u003e Ecological modelling approaches for predicting emergent properties in microbial communities. \u003cem\u003eNat Ecol Evol\u003c/em\u003e \u003cstrong\u003e6\u003c/strong\u003e, 855\u0026ndash;865 (2022).\u003c/li\u003e\n\u003cli\u003eEstrela, S. \u0026amp; Brown, S. P. Metabolic and demographic feedbacks shape the emergent spatial structure and function of microbial communities. \u003cem\u003ePLoS Comput Biol\u003c/em\u003e \u003cstrong\u003e9\u003c/strong\u003e, e1003398 (2013).\u003c/li\u003e\n\u003cli\u003eYang, N. Emergent bacterial community properties induce enhanced drought tolerance in \u003cem\u003eArabidopsis\u003c/em\u003e. \u003cem\u003enpj Biofilms and Microbiomes\u003c/em\u003e \u003cstrong\u003e7,\u003c/strong\u003e 82 (2021).\u003c/li\u003e\n\u003cli\u003eXu, Z. \u003cem\u003eet al.\u003c/em\u003e Enhanced control of plant wilt dsease by a xylose-inducible \u003cem\u003edegQ\u003c/em\u003e gene engineered into \u003cem\u003eBacillus velezensis\u003c/em\u003e strain SQR9XYQ. \u003cem\u003ePhytopathology\u003c/em\u003e \u003cstrong\u003e109\u003c/strong\u003e, 36\u0026ndash;43 (2019).\u003c/li\u003e\n\u003cli\u003eZelezniak, A. \u003cem\u003eet al.\u003c/em\u003e Metabolic dependencies drive species co-occurrence in diverse microbial communities. \u003cem\u003eProc Natl Acad Sci USA\u003c/em\u003e \u003cstrong\u003e112\u003c/strong\u003e, 6449\u0026ndash;6454 (2015).\u003c/li\u003e\n\u003cli\u003eCulp, E. J. \u0026amp; Goodman, A. L. Cross-feeding in the gut microbiome: Ecology and mechanisms. \u003cem\u003eCell Host Microbe\u003c/em\u003e \u003cstrong\u003e31\u003c/strong\u003e, 485\u0026ndash;499 (2023).\u003c/li\u003e\n\u003cli\u003eD\u0026rsquo;Souza, G. \u003cem\u003eet al.\u003c/em\u003e Ecology and evolution of metabolic cross-feeding interactions in bacteria. \u003cem\u003eNat. Prod. Rep.\u003c/em\u003e \u003cstrong\u003e35\u003c/strong\u003e, 455\u0026ndash;488 (2018).\u003c/li\u003e\n\u003cli\u003ePande, S. \u003cem\u003eet al.\u003c/em\u003e Privatization of cooperative benefits stabilizes mutualistic cross-feeding interactions in spatially structured environments. \u003cem\u003eISME J\u003c/em\u003e \u003cstrong\u003e10\u003c/strong\u003e, 1413\u0026ndash;1423 (2016).\u003c/li\u003e\n\u003cli\u003eLynch, C. J. \u0026amp; Adams, S. H. Branched-chain amino acids in metabolic signalling and insulin resistance. \u003cem\u003eNat Rev Endocrinol\u003c/em\u003e \u003cstrong\u003e10\u003c/strong\u003e, 723\u0026ndash;736 (2014).\u003c/li\u003e\n\u003cli\u003eGulick, A. M. Nonribosomal peptide synthetase biosynthetic clusters of ESKAPE pathogens. \u003cem\u003eNat Prod Rep\u003c/em\u003e \u003cstrong\u003e34\u003c/strong\u003e, 981\u0026ndash;1009 (2017).\u003c/li\u003e\n\u003cli\u003eBrosnan, J. T. \u0026amp; Brosnan, M. E. Branched-chain amino acids: enzyme and substrate regulation. \u003cem\u003eJ Nutr\u003c/em\u003e \u003cstrong\u003e136\u003c/strong\u003e, 207S-211S (2006).\u003c/li\u003e\n\u003cli\u003eSch\u0026auml;fer, M. \u003cem\u003eet al.\u003c/em\u003e Metabolic interaction models recapitulate leaf microbiota ecology. \u003cem\u003eScience\u003c/em\u003e \u003cstrong\u003e381\u003c/strong\u003e, eadf5121 (2023).\u003c/li\u003e\n\u003cli\u003eSokolovskaya, O. M., Shelton, A. N. \u0026amp; Taga, M. E. Sharing vitamins: Cobamides unveil microbial interactions. \u003cem\u003eScience\u003c/em\u003e \u003cstrong\u003e369\u003c/strong\u003e, eaba0165 (2020).\u003c/li\u003e\n\u003cli\u003eCh\u0026rsquo;ng, J.-H. \u003cem\u003eet al.\u003c/em\u003e Heme cross-feeding can augment \u003cem\u003eStaphylococcus aureus\u003c/em\u003e and \u003cem\u003eEnterococcus faecalis\u003c/em\u003e dual species biofilms. \u003cem\u003eISME J \u003c/em\u003e\u003cstrong\u003e16\u003c/strong\u003e, 2015\u0026ndash;2026 (2022).\u003c/li\u003e\n\u003cli\u003eWang, W. Narrow-spectrum resource-utilizing bacteria drive the stability of synthetic communities through enhancing metabolic interactions. \u003cem\u003eNat Commun\u003c/em\u003e \u003cstrong\u003e16\u003c/strong\u003e, 6088 (2025.)\u003c/li\u003e\n\u003cli\u003eAndrews, S. C., Robinson, A. K. \u0026amp; Rodr\u0026iacute;guez-Qui\u0026ntilde;ones, F. Bacterial iron homeostasis. \u003cem\u003eFEMS Microbiol Rev\u003c/em\u003e \u003cstrong\u003e27\u003c/strong\u003e, 215\u0026ndash;237 (2003).\u003c/li\u003e\n\u003cli\u003eTangen, K. L., Jung, W. H., Sham, A. P., Lian, T. \u0026amp; Kronstad, J. W. The iron- and cAMP-regulated gene SIT1 influences ferrioxamine B utilization, melanization and cell wall structure in Cryptococcus neoformans. \u003cem\u003eMicrobiology \u003c/em\u003e\u003cstrong\u003e153\u003c/strong\u003e, 29\u0026ndash;41 (2007).\u003c/li\u003e\n\u003cli\u003eCaza, M. \u0026amp; Kronstad, J. Shared and distinct mechanisms of iron acquisition by bacterial and fungal pathogens of humans. \u003cem\u003eFront. Cell. Infect. Microbiol.\u003c/em\u003e \u003cstrong\u003e3\u003c/strong\u003e, (2013).\u003c/li\u003e\n\u003cli\u003eNeilands, J. B. Iron absorption and transport in microorganisms. \u003cem\u003eAnnu. Rev. Nutr.\u003c/em\u003e \u003cstrong\u003e1\u003c/strong\u003e, 27\u0026ndash;46 (1981).\u003c/li\u003e\n\u003cli\u003eSaleem, A. \u003cem\u003eet al.\u003c/em\u003e Alkaline and acidic soil constraints on iron accumulation by Rice cultivars in relation to several physio-biochemical parameters. \u003cem\u003eBMC Plant Biol\u003c/em\u003e \u003cstrong\u003e23\u003c/strong\u003e, 397 (2023).\u003c/li\u003e\n\u003cli\u003eNeilands, J. B. Siderophores: structure and function of microbial iron transport compounds . \u003cem\u003eJ Biol Chem\u003c/em\u003e \u003cstrong\u003e270\u003c/strong\u003e, 26723\u0026ndash;26726 (1995).\u003c/li\u003e\n\u003cli\u003eMiethke, M. \u0026amp; Marahiel, M. A. Siderophore-based iron acquisition and pathogen control. \u003cem\u003eMicrobiol Mol Biol Rev\u003c/em\u003e \u003cstrong\u003e71\u003c/strong\u003e, 413\u0026ndash;451 (2007).\u003c/li\u003e\n\u003cli\u003eWei, Z. \u003cem\u003eet al.\u003c/em\u003e Microbial siderophores for One Health. \u003cem\u003eTrends Microbiol\u003c/em\u003e \u003cstrong\u003e33\u003c/strong\u003e, 1277-1285 (2025).\u003c/li\u003e\n\u003cli\u003eLiu, Y. \u003cem\u003eet al.\u003c/em\u003e Plant commensal type VII secretion system causes iron leakage from roots to promote colonization. \u003cem\u003eNat Microbiol\u003c/em\u003e \u003cstrong\u003e8\u003c/strong\u003e, 1434\u0026ndash;1449 (2023).\u003c/li\u003e\n\u003cli\u003eTan, T. \u003cem\u003eet al.\u003c/em\u003e Siderophore-mediated iron enrichment in the biofilm matrix enhances plant iron nutrition. \u003cem\u003eCell Rep\u003c/em\u003e \u003cstrong\u003e44\u003c/strong\u003e, 116481 (2025).\u003c/li\u003e\n\u003cli\u003eXu, Z. \u003cem\u003eet al.\u003c/em\u003e Antibiotic bacillomycin D affects iron acquisition and biofilm formation in bacillus velezensis through a btr-mediated FeuABC-dependent pathway. \u003cem\u003eCell Rep\u003c/em\u003e \u003cstrong\u003e29\u003c/strong\u003e, 1192-1202.e5 (2019).\u003c/li\u003e\n\u003cli\u003eRizzi, A. \u003cem\u003eet al.\u003c/em\u003e Bacillus subtilis modulates its usage of biofilm-bound iron in response to environmental iron availability. \u003cem\u003eAppl Environ Microbiol\u003c/em\u003e \u003cstrong\u003e86\u003c/strong\u003e, e00944-20 (2020).\u003c/li\u003e\n\u003cli\u003eQin, Y. \u003cem\u003eet al.\u003c/em\u003e Heterogeneity in respiratory electron transfer and adaptive iron utilization in a bacterial biofilm. \u003cem\u003eNat Commun\u003c/em\u003e \u003cstrong\u003e10\u003c/strong\u003e, 3702 (2019).\u003c/li\u003e\n\u003cli\u003eGu, S. \u003cem\u003eet al.\u003c/em\u003e Competition for iron drives phytopathogen control by natural rhizosphere microbiomes. \u003cem\u003eNat Microbiol\u003c/em\u003e \u003cstrong\u003e5\u003c/strong\u003e, 1002\u0026ndash;1010 (2020).\u003c/li\u003e\n\u003cli\u003eHarrison, F., Paul, J., Massey, R. C. \u0026amp; Buckling, A. Interspecific competition and siderophore-mediated cooperation in \u003cem\u003ePseudomonas aeruginosa\u003c/em\u003e. \u003cem\u003eISME J\u003c/em\u003e \u003cstrong\u003e2\u003c/strong\u003e, 49\u0026ndash;55 (2008).\u003c/li\u003e\n\u003cli\u003eJenul, C. \u003cem\u003eet al.\u003c/em\u003e Pyochelin biotransformation by \u003cem\u003eStaphylococcus aureus \u003c/em\u003eshapes bacterial competition with \u003cem\u003ePseudomonas aeruginosa\u003c/em\u003e in polymicrobial infections. \u003cem\u003eCell Rep\u003c/em\u003e \u003cstrong\u003e42\u003c/strong\u003e, (2023).\u003c/li\u003e\n\u003cli\u003eButaitė, E., Baumgartner, M., Wyder, S. \u0026amp; K\u0026uuml;mmerli, R. Siderophore cheating and cheating resistance shape competition for iron in soil and freshwater \u003cem\u003ePseudomonas \u003c/em\u003ecommunities. \u003cem\u003eNat Commun\u003c/em\u003e \u003cstrong\u003e8\u003c/strong\u003e, 414 (2017).\u003c/li\u003e\n\u003cli\u003eD\u0026rsquo;Onofrio, A. \u003cem\u003eet al.\u003c/em\u003e Siderophores from neighboring organisms promote the growth of uncultured bacteria. \u003cem\u003eChem Biol\u003c/em\u003e \u003cstrong\u003e17\u003c/strong\u003e, 254\u0026ndash;264 (2010).\u003c/li\u003e\n\u003cli\u003eSun, X. \u003cem\u003eet al.\u003c/em\u003e Metabolic interactions affect the biomass of synthetic bacterial biofilm communities. \u003cem\u003emSystems\u003c/em\u003e \u003cstrong\u003e8\u003c/strong\u003e, e01045-23 (2023).\u003c/li\u003e\n\u003cli\u003eWang, P. \u003cem\u003eet al.\u003c/em\u003e An approach to constructing multispecies biofilm communities from rhizosphere soil. \u003cem\u003eJoVE\u003c/em\u003e e66926 (2024) doi:10.3791/66926.\u003c/li\u003e\n\u003cli\u003eMiethke, M. \u003cem\u003eet al.\u003c/em\u003e Ferri-bacillibactin uptake and hydrolysis in \u003cem\u003eBacillus subtilis\u003c/em\u003e. \u003cem\u003eMol Microbiol\u003c/em\u003e \u003cstrong\u003e61\u003c/strong\u003e, 1413\u0026ndash;1427 (2006).\u003c/li\u003e\n\u003cli\u003eLewenza S, Sokol PA. Regulation of ornibactin biosynthesis and N-Acyl-l-Homoserine lactone production by CepR in \u003cem\u003eBurkholderia cepacia\u003c/em\u003e. \u003cem\u003eJ Bacteriol \u003c/em\u003e\u003cstrong\u003e183\u003c/strong\u003e, 7 (2001). \u003c/li\u003e\n\u003cli\u003eSun, X. \u003cem\u003eet al.\u003c/em\u003e \u003cem\u003eBacillus velezensis\u003c/em\u003e stimulates resident rhizosphere \u003cem\u003ePseudomonas stutzeri \u003c/em\u003efor plant health through metabolic interactions. \u003cem\u003eISME J\u003c/em\u003e \u003cstrong\u003e16\u003c/strong\u003e, 774\u0026ndash;787 (2022).\u003c/li\u003e\n\u003cli\u003eLi, Z.-L., Wang, Y.-H., Chu, J., Zhuang, Y.-P. \u0026amp; Zhang, S.-L. Effect of branched-chain amino acids, valine, isoleucine and leucine on the biosythesis of bitespiramycin 4\u0026quot;-O-acylspiramycins. \u003cem\u003eBraz. J. Microbiol.\u003c/em\u003e \u003cstrong\u003e40\u003c/strong\u003e, 734\u0026ndash;746 (2009).\u003c/li\u003e\n\u003cli\u003eBlomstrand, E., Eliasson, J., Karlsson, H. K. \u0026amp; K\u0026ouml;hnke, R. Branched-chain amino acids activate key enzymes in protein synthesis after physical exercise. \u003cem\u003eJ Nutr \u003c/em\u003e\u003cstrong\u003e136\u003c/strong\u003e, 269S-273S (2006).\u003c/li\u003e\n\u003cli\u003ePiovia-Scott, J. \u003cem\u003eet al.\u003c/em\u003e Greater species richness of bacterial skin symbionts better suppresses the amphibian fungal pathogen \u003cem\u003eBatrachochytrium dendrobatidis\u003c/em\u003e. \u003cem\u003eMicrob Ecol\u003c/em\u003e \u003cstrong\u003e74\u003c/strong\u003e, 217\u0026ndash;226 (2017).\u003c/li\u003e\n\u003cli\u003eHu, J. \u003cem\u003eet al.\u003c/em\u003e Introduction of probiotic bacterial consortia promotes plant growth via impacts on the resident rhizosphere microbiome. \u003cem\u003eProc Biol Sci\u003c/em\u003e \u003cstrong\u003e13\u003c/strong\u003e, 1960 (2021).\u003c/li\u003e\n\u003cli\u003eMee, M. T., Collins, J. J., Church, G. M. \u0026amp; Wang, H. H. Syntrophic exchange in synthetic microbial communities. \u003cem\u003eProc Natl Acad Sci USA.\u003c/em\u003e \u003cstrong\u003e111\u003c/strong\u003e, (2014).\u003c/li\u003e\n\u003cli\u003eMataigne, V., Vannier, N., Vandenkoornhuyse, P. \u0026amp; Hacquard, S. Multi-genome metabolic modeling predicts functional inter-dependencies in the \u003cem\u003eArabidopsis\u003c/em\u003e root microbiome. \u003cem\u003eMicrobiome\u003c/em\u003e \u003cstrong\u003e10\u003c/strong\u003e, 217 (2022).\u003c/li\u003e\n\u003cli\u003eQuillin, S. J., Tran, P. \u0026amp; Prindle, A. Potential roles for gamma-aminobutyric acid signaling in bacterial communities. \u003cem\u003eBioelectricity\u003c/em\u003e \u003cstrong\u003e3\u003c/strong\u003e, 120\u0026ndash;125 (2021).\u003c/li\u003e\n\u003cli\u003eChallis, G. L. \u0026amp; Naismith, J. H. Structural aspects of non-ribosomal peptide biosynthesis. \u003cem\u003eCurr Opin Struct Biol\u003c/em\u003e \u003cstrong\u003e14\u003c/strong\u003e, 748\u0026ndash;756 (2004).\u003c/li\u003e\n\u003cli\u003eMay, J. J., Kessler, N., Marahiel, M. A. \u0026amp; Stubbs, M. T. Crystal structure of DhbE, an archetype for aryl acid activating domains of modular nonribosomal peptide synthetases. \u003cem\u003eProc Natl Acad Sci USA.\u003c/em\u003e \u003cstrong\u003e99\u003c/strong\u003e, 12120\u0026ndash;12125 (2002).\u003c/li\u003e\n\u003cli\u003eMomeni, B., Waite, A. J. \u0026amp; Shou, W. Spatial self-organization favors heterotypic cooperation over cheating. \u003cem\u003eeLife\u003c/em\u003e \u003cstrong\u003e2\u003c/strong\u003e, e00960 (2013).\u003c/li\u003e\n\u003cli\u003eEstrela, S., Morris, J. J. \u0026amp; Kerr, B. Private benefits and metabolic conflicts shape the emergence of microbial interdependencies. \u003cem\u003eEnviron Microbiol\u003c/em\u003e \u003cstrong\u003e18\u003c/strong\u003e, 1415\u0026ndash;1427 (2016).\u003c/li\u003e\n\u003cli\u003eHarper, A. E. \u0026amp; Miller, R. H. Branched-chain amino acid metabolism. \u003cem\u003eAnnu Rev Nutr\u003c/em\u003e \u003cstrong\u003e4\u003c/strong\u003e, 409\u0026ndash;54 (1984).\u003c/li\u003e\n\u003cli\u003eDertz, E. A., Stintzi, A. \u0026amp; Raymond, K. N. Siderophore-mediated iron transport in \u003cem\u003eBacillus subtilis\u003c/em\u003e and \u003cem\u003eCorynebacterium glutamicum\u003c/em\u003e. \u003cem\u003eJ Biol Inorg Chem\u003c/em\u003e \u003cstrong\u003e11\u003c/strong\u003e, 1087\u0026ndash;1097 (2006).\u003c/li\u003e\n\u003cli\u003eAhmad, F., Ahmad, I. \u0026amp; Khan, M. S. Screening of free-living rhizospheric bacteria for their multiple plant growth promoting activities. \u003cem\u003eMicrobiol Res\u003c/em\u003e \u003cstrong\u003e163\u003c/strong\u003e, 173\u0026ndash;181 (2008).\u003c/li\u003e\n\u003cli\u003eMasuko, T. \u003cem\u003eet al.\u003c/em\u003e Carbohydrate analysis by a phenol\u0026ndash;sulfuric acid method in microplate format. \u003cem\u003eAnall Biochem\u003c/em\u003e \u003cstrong\u003e339\u003c/strong\u003e, 69\u0026ndash;72 (2005).\u003c/li\u003e\n\u003cli\u003eNautiyal, C. S. An efficient microbiological growth medium for screening phosphate solubilizing microorganisms. \u003cem\u003eFEMS Microbiol Lett\u003c/em\u003e \u003cstrong\u003e170\u003c/strong\u003e, 265\u0026ndash;270 (1999).\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"nature-portfolio","isNatureJournal":true,"hasQc":false,"allowDirectSubmit":false,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"","title":"Nature Portfolio","twitterHandle":"","acdcEnabled":false,"dfaEnabled":false,"editorialSystem":"ejp","reportingPortfolio":"","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"multispecies biofilm, microbial community, siderophore, cross-feeding, plant growth promotion, emerging function","lastPublishedDoi":"10.21203/rs.3.rs-8956555/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-8956555/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"Biofilm communities exhibit emergent properties that exceed the sum of contributions from individual members of the community. Here, we describe a multilayered metabolic interaction that drives enhanced biofilm formation among three bacterial species from the plant rhizosphere. Comparative metatranscriptomic and metabolomic analyses reveal that Bacillus velezensis-secreted 5-aminovaleric acid promotes the growth of the other community members, Burkholderia contaminans and Acinetobacter baumannii. In return, B. contaminans supplies branched-chain amino acids for B. velezensis. Branched-chain amino acids and cell–cell signaling acyl-homoserine lactones from B. contaminans induce biosynthesis of the siderophore bacillibactin in B. velezensis, that is further enhanced by A. baumannii. In exchange, the B. velezensis-secreted siderophore promotes the growth of B. contaminans in iron-limited conditions, which benefits the multispecies biofilm community in vitro and promotes plant growth performance in iron-depleted soil. Our study reveals the molecular mechanisms underlying an emergent rhizosphere biofilm community function and demonstrates its importance in plant–microbe interactions.","manuscriptTitle":"Metabolic exchange and siderophore sharing underlie emergent biofilm synergism","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2026-03-12 09:19:15","doi":"10.21203/rs.3.rs-8956555/v1","editorialEvents":[],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"nature-microbiology","isNatureJournal":true,"hasQc":false,"allowDirectSubmit":false,"externalIdentity":"nmicrobiol","sideBox":"Learn more about [Nature Microbiology](http://www.nature.com/nmicrobiol/)","snPcode":"","submissionUrl":"","title":"Nature Microbiology","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"ejp","reportingPortfolio":"Nature Research","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"ed7ced57-45c2-4826-8f45-dbae29f0d667","owner":[],"postedDate":"March 12th, 2026","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"under-review","subjectAreas":[{"id":64279952,"name":"Biological sciences/Microbiology/Microbial communities/Microbial ecology"},{"id":64279953,"name":"Biological sciences/Ecology/Microbial ecology"}],"tags":[],"updatedAt":"2026-04-28T10:18:45+00:00","versionOfRecord":[],"versionCreatedAt":"2026-03-12 09:19:15","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-8956555","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-8956555","identity":"rs-8956555","version":["v1"]},"buildId":"XKTyCvWXoU3ODBz1xrDgd","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}