Microbial inoculants modulate the rhizosphere microbiome, alleviate plant stress responses and enhance maize growth at field scale

doi:10.21203/rs.3.rs-4993260/v1

Download PDF

Research Article

Microbial inoculants modulate the rhizosphere microbiome, alleviate plant stress responses and enhance maize growth at field scale

https://doi.org/10.21203/rs.3.rs-4993260/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background

Field inoculation of crops with beneficial microbes is a promising sustainable strategy to enhance plant fitness and nutrient acquisition. However, effectiveness can vary due to environmental factors, microbial competition, and methodological challenges, while their precise modes of action remain uncertain. This underscores the need for further research to optimize inoculation strategies for consistent agricultural benefits.

Results

Using a comprehensive, multidisciplinary approach, we investigated the effects of a consortium of beneficial microbes (BMc) (Pseudomonas sp. RU47, Bacillus atrophaeus ABi03, Trichoderma harzianum OMG16) on maize (cv. Benedictio) through an inoculation experiment conducted within a long-term field trial across intensive and extensive farming practices. Additionally, an unexpected early drought stress emerged as a climatic variable, offering further insight into the effectiveness of the microbial consortium. Our findings demonstrate that BMc inoculation primarily enhanced plant growth and fitness, particularly by increasing iron uptake, which is crucial for drought adaptation. Inoculated maize plants showed improved shoot growth and fitness compared to non-inoculated plants, regardless of farming practices. Specifically, BMc modulated plant hormonal balance, enhanced the detoxification of reactive oxygen species (ROS), and increased root exudation of iron-chelating metabolites. Amplicon sequencing revealed shifts in rhizosphere bacterial and fungal communities mediated by the consortium, while metagenomic shotgun sequencing indicated enrichment of genes related to antimicrobial lipopeptides and siderophores.

Conclusions

Our findings highlight the multifaceted benefits of BMc inoculation on plant fitness, significantly influencing metabolism, stress responses, and the rhizosphere microbiome. These improvements are crucial for advancing sustainable agricultural practices by enhancing plant resilience and productivity.

Modern agriculture relies on intensive tillage, fertilization, and pesticide applications to achieve high crop yields. However, these practices carry potential for environmental risks, including promotion of soil erosion, ecosystem eutrophication, biodiversity decline, and impairment of soil integrity^1–3. The implementation of sustainable farming strategies to maintain agro-ecosystem functions and services has been proposed to address these concerns. The integration of reduced tillage and fertilization intensity, while considering the beneficial interaction between plants and rhizosphere microorganisms, has potential to stabilize crop yield, plant fitness and resilience, and soil health and fertility^4–7. Indeed, such practices also entail risks for nutrient limitation or increased pathogen/pest pressure that may, due to still limited sustainable options, require conventional agricultural measures^8–10. Considering the positive effects of beneficial microorganisms (BMs) on plant performance^11–15, their use in farming offers a chance to exploit the interplay between plants and beneficial rhizosphere microbes ^16–18, which also have the potential to alleviate plant responses to biotic and abiotic stresses^19–22. Host plants typically recruit beneficial microorganisms primarily through root exudates that generally serve as chemo-attractants or repellants and deliver nutrients that facilitate microbial rhizosphere colonization. However, the efficacy of BM inoculants at field scale was often inconsistent^23–25. This remains a major challenge for successful agricultural application. Insufficient root colonization (either root surface, intercellular space or intracellular) and unfavorable environmental conditions are major factors that may hinder the interaction between crops and BMs^26,27. For optimal rhizosphere competence, inoculated BMs must compete with the indigenous soil microbial community and adapt to diverse and occasionally harsh environmental conditions.

BMs often show synergisms of their potential functions in the rhizosphere with enhanced promotion of plant performance. Consequently, the application of BM consortia (BMc) was recognized as an opportunity to exploit such synergisms. Especially Bacillus and Pseudomonas strains with plant-beneficial activity often show synergisms, and their combined inoculation can promote plant growth²⁸. In addition, Trichoderma strains can interact positively with Bacillus strains or show synergistic effects upon co-inoculation despite belonging to different kingdoms²⁹. In a previous field experiment we investigated the ability of a BM consortium (BMc) to establish in the winter rye rhizosphere and assessed its impact on plant performance under various farming practices¹³. Composed of Pseudomonas sp. (RU47), Bacillus atrophaeus (ABi03), and Trichoderma harzianum (OMG16), the BMc demonstrated a high capacity for winter rye rhizosphere colonization. These microorganisms persisted through winter dormancy and significantly enhanced plant nutrient status and performance, particularly under organic farming conditions, which provide less nutrients to the plants than intensive conventional fertilization. Reflecting a broad host range, plant growth-promoting properties of the selected strains have additionally been demonstrated for various crops including lettuce^30,31, oilseed rape³², tomato³⁶ and maize³⁷. Combining different microbial strains as BM consortia (BMc) with complementary properties is discussed as an approach to increase the efficiency and flexibility of BM-based production strategies under variable environmental conditions³⁵. For example, Trichoderma-Bacillus combinations are well-documented examples in this context^36,37. Consequently, the use of BMc for sustainable farming, particularly emphasizing reduced agrochemical input and tillage, may be crucial for improving crop health and yield stability.

Seeking to comprehensively examine the impact of the above mentioned BMc on plant performance and to test the hypothesis that BMc inoculation mitigates growth limitation associated with reduced N-fertilization intensity, we conducted a systematic inoculation trial in maize within a long-term field experiment (LTE) established in 1992. This LTE employs intensive (recommended N-fertilization intensity and pesticide application) and extensive (reduced N-fertilization intensity without fungicide application) farming practice within a crop rotation (wheat, maize, wheat, rapeseed, barley) using cultivator tillage. To rigorously assess the impact of BMc inoculation on maize performance, we evaluated plant traits (growth, physiological stress indicators, phytohormonal profiles, and nutrient status), as well as root and rhizosphere metabolites. Additionally, we investigated how BMc inoculation correlated with structural and functional changes within rhizosphere bacterial and fungal communities, utilizing amplicon and metagenome shotgun sequencing. Severe drought conditions during the experimental year, especially in the juvenile maize establishment phase, provided a unique opportunity to investigate the BMc effects under substantial abiotic stress. This pioneering study represents the first field investigation of the impact of BMc inoculation on the soil/microbe/plant system under contrasting mineral N-fertilization intensities, employing a multidisciplinary framework encompassing diverse analytical techniques.

Description of the experimental design

An inoculation experiment was conducted within a long-term field trial (LTE) to investigate the effects of a consortium of beneficial microorganisms (BMc; Pseudomonas sp. RU47, Bacillus atrophaeus ABi03, and Trichoderma harzianum OMG16) on maize performance (cv. Benedictio). The experimental plots employed conservation tillage and received either the recommended use of nitrogen (N) and pesticides for maize (intensive, Int) or reduced N use and no fungicide application (extensive, Ext). Maize plants were drenched twice with the BMc at growth stages EC 12 and 14–16 and harvested at maturity (EC 53–63). During the experiment duration, spring precipitation was exceptionally low compared to the average spring rainfall of the past 30 years, with almost no rainfall in April and less in May and July (Fig. S1). To elucidate the effects of the BMc on the maize-microbiome interaction in the rhizosphere, we analyzed various parameters, including root traits (total root length, root hair length), shoot dry mass (SDM), physiological stress indicators, phytohormones, nutrient status, root and rhizosphere metabolites, and bacterial/archaeal and fungal community compositions, as well as soil metagenome profiles.

BMc strains colonized the maize rhizosphere and increased shoot dry mass and iron uptake

The inoculated strains effectively colonized the rhizosphere, as determined by CFU counts through selective plating (Fig. 1A). Among the BMc strains, ABi03 exhibited the most efficient colonization with a density two orders of magnitude higher than RU47 and OMG (Two-Way ANOVA, F = 1889.68, p < 0.05, n = 8; Fig. 1A). N-fertilization intensity influenced the ABi03 rhizo-competence (Two-Way ANOVA, F = 16.45, p < 0.05, n = 8; Fig. 1A), but had no significant effect on RU47 and OMG16 (Two-Way ANOVA, F = 1.2–0.8, p > 0.05, n = 8; Fig. 1A). Maize SDM was significantly higher by approximately 30% in BMc inoculated plants compared to non-inoculated controls (Two-Way ANOVA, F = 149.7, p < 0.05, n = 8; Fig. 1B). The increase in SDM was accompanied by higher concentrations of sulfur, iron, zinc and manganese in maize leaves of inoculated plants (Two-Way ANOVA, F = 0.3-17.14, p < 0.05, n = 8; Table S1 and Table 1). However, only iron concentration showed an increase regardless of N-fertilization intensity (Two-Way ANOVA, F = 17.14, Tukey’s HSD, p < 0.05, Table 1). Moreover, none of the nutrients (except N in extensive fertilized control plants) showed concentrations below the deficiency limit (Table 1). In contrast, we did not observe any significant effects of the experimental variables on most of the soil measured chemical parameters (Table S2). Furthermore, BMc inoculation did not significantly affect root growth (total root and hair length) or arbuscular mycorrhizal fungi (AMF) colonization, processes crucial for the spatial acquisition of mineral nutrients (Fig. S2). However, AMF root colonization was significantly reduced under intensive N-fertilization, supporting the notion that high N-fertilization and pesticide use have detrimental effects on mycorrhizal symbiosis^38,39.

Table 1

Macro- and micro-nutrient concentrations measured in maize shoots (cv. Benedictio) grown under extensive and intensive N-fertilization intensities, without (Ctrl; control plants) and with BMc (beneficial microorganisms consortium) inoculation. DT indicates deficiency threshold based on Campbell (2000). Data represent mean values ± standard deviation. Different letters indicate significant differences (p < 0.05, Tukey’s HSD test, Benjamini-Hochberg correction, n = 4).
		Extensive						Intensive
	DT	Ctrl			BMc			Ctrl			BMc
Macro-nutrients [g kg^− 1]
C		453	± 0		454	± 3.00		453	± 1.00		454	± 3.00
N	30	29.7	± 0.30		31.0	± 0.80		31.1	± 1.40		31.8	± 1.20
P	2.5	2.72	± 0.16		2.83	± 0.14		2.57	± 0.10		2.59	± 0.13
K	20	22.1	± 0.80		22.6	± 1.50		22.0	± 0.90		21.1	± 1.00
Mg	2.5	1.84	± 0.07	a	1.85	± 0.05	a	1.72	± 0.07	ab	1.63	± 0.09	b
Ca	4	5.00	± 0.36		5.03	± 0.32		4.83	± 0.20		5.05	± 0.72
S	1.2	1.95	± 0.06	b	2.13	± 0.10	ab	2.14	± 0.08	ab	2.23	± 0.07	a
Micro-nutrients [mg kg^− 1]
Fe	15	85.0	± 5.0	bc	96.2	± 4.70	a	80.2	± 5.70	c	91.7	± 6.30	ab
Cu	5	8.78	± 0.32		8.41	± 0.48		10.07	± 1.01		9.29	± 0.77
Zn	15	37.5	3.3	b	40.8	± 4.50	ab	40.9	± 3.90	ab	47.0	± 3.60	a
Mn	15	59.1	± 5.0	b	61.7	± 6.20	ab	62.4	± 4.60	ab	77.3	± 16.4	a
Campbell CR (ed). Reference Sufficiency Ranges for Plant Analysis in the Southern Region of the United States. Southern Cooperative Series Bulletin No. 2000, 394.

In response to BMc inoculation, a substantial number of stress-related maize genes, in terms of relative transcript abundances, tended to be repressed in the leaves regardless of N-fertilization level (Fig. 1C). In contrast, stress-related genes encoding for iron-dependent enzymes (Table S3), involved in detoxification of reactive oxygen species (ROS), such as superoxide dismutase (SOD) and ascorbate peroxidase (APX) activity, and genes associated with nutrient acquisition pathways, exhibited distinct upregulation in response to BMc inoculation, particularly those involved in iron uptake. Specifically, NAS3, a gene associated with the synthesis of the iron/zinc chelator nicotianamine (Table S3) likely contributing to enhanced iron concentrations in leaves, was consistently upregulated (t-test, p < 0.05, n = 4; Fig. 1C). The iron transporter IRTa showed a contrasting pattern, albeit statistically significant only in the context of extensive N-fertilization samples (t-test, p < 0.05, n = 4).

BMc inoculation exerted systemic effects on maize hormonal status and improved stress resilience independent of N-fertilization intensity

BMc inoculation significantly increased the concentration of growth-promoting phytohormones (indole acetic acid, gibberellic acid, and zeatin) in the shoots by approximately twofold, while decreasing the concentrations of stress-related hormones (abscisic acid (ABA) and jasmonic acid (JA)) (Two-Way ANOVA, F = 28.20-279.34, p < 0.05, Fig. 2). Salicylic acid (SA) concentrations in shoots remained unaffected by BMc inoculation (Fig. 2). In the roots, a similar response of plant-growth hormones to BMc inoculation was observed, with significantly higher concentrations of indole acetic acid (IAA) and zeatin (cytokinin) in the inoculated plants compared to the control (Fig. S3). Interestingly, BMc inoculation led to a significant increase in the concentration of stress-related hormones involved in plant defense in the roots (JA and SA, t-test, p < 0.05, n = 4, Fig. S3), which could potentially be associated with the plant response to the presence of a high microbial load or the drastic change in microbial community composition in the maize rhizosphere.

Since BMc inoculation enhanced the gene expression and shoot activities of iron-dependent ROS detoxification (APX and SOD, ANOVA, F = 45.09, p < 0.05, Fig. 3), consequently it halved the accumulation of hydrogen peroxide (H₂O₂) in the leaf tissue (ANOVA, F = 285.5, p < 0.05, Fig. 3), and additionally increased leaf concentrations of total antioxidant and glycine betaine (ANOVA, F = 12.54, p < 0.05, Fig. 3), which have functions in non-enzymatic ROS detoxification and osmotic adjustment. The observed decrease in stress-resilience markers within leaf tissues of BMc-inoculated maize plants aligns with the downregulation of stress-related genes in this tissue, indicating a plausible involvement of BMc inoculation in regulating plant stress responses.

The patterns of rhizosphere metabolites are shaped by BMc inoculation

To evaluate the impact of BMc inoculation on rhizosphere metabolites, we established experimental plots with root observation windows that enabled the sampling and profiling of low molecular weight metabolites across distinct maize root zones (Fig. S4). As expected, metabolite concentrations exhibited significant variation across the sampled root zones, reflecting the distinct physiological and functional characteristics of different root types (Fig. S5, Table S4). Given the higher detection of secondary metabolites in the apical root zone, likely exerting a significant effect on microbiota assembly around maize roots in contrast to other root types, our investigation specifically targeted the analysis of rhizosphere metabolites from the apical root zone (Fig. S5). In line with the improved growth performance exhibited by BMc-inoculated plants, a majority of the identified metabolites showed significant increases following BMc inoculation, irrespective of N-fertilization intensity (p < 0.05, Tukey’s HSD test, Fig. 4). These enhancements encompassed organic acid anions (predominantly malate), sugars (such as glucose and trehalose), specific phenolic acids and flavonoids (like caffeic acid, quercetin/naringenin), as well as benzoxazinoid metabolites with antibiotic, allelopathic, and iron-mobilizing properties (e.g., MBOA) (Fig. 4).

BMc inoculation increased Myxococcota, Bacteroidota, and Proteobacteria in the maize rhizosphere

BMc inoculation significantly impacted bacterial β-diversity within the maize rhizosphere, suggesting its potential role as a determinant factor shaping bacterial community assembly as revealed by 16S rRNA gene amplicon sequencing (PERMANOVA, BMc: R² = 0.12, p = 1x10^− 4, n = 8; Fig. 5A). Actinobacteriota, Proteobacteria, and Firmicutes were the most abundant phyla, while Acidobacteria, Bacteroidota, Gemmatimonadota, Chloroflexi, Myxococcota, and Verrucomicrobiota occurred in relative abundances below 5% (Fig. 6A). BMc inoculation specifically decreased the relative abundance of Actinobacteriota in both N-fertilization intensities, while increasing the relative abundance of Acidobacteria, Bacteroidota, Chloroflexi, and Myxococcota under intensive fertilization and Proteobacteria under extensive fertilization (logistic regression, p < 0.05, n = 4; Fig. 6A). Additionally, BMc inoculation reduced the relative abundance of Gemmatimonadota in the rhizosphere of maize grown under intensive fertilization (logistic regression, p < 0.05, n = 4; Fig. 6A).

BMc inoculation significantly increased the relative abundance of 55 and 48 ASVs in intensive and extensive fertilization, respectively, out of a total of 6595 bacterial ASVs (logistic regression, p < 0.05, Benjamini-Hochberg correction, n = 4; Table S5). Concurrently, the relative abundance of 48 and 33 bacterial ASVs decreased due to BMc inoculation under intensive or extensive N-fertilization, respectively (logistic regression, p < 0.05, Benjamini-Hochberg correction, n = 4, Table S5). Out of the 15 most abundant bacterial ASVs identified in the maize rhizosphere, the relative abundance of five of these ASVs, classified as four genera belonging to the Actinobacteriota phylum (Agromyces: ASV3, Pseudarthrobacter: ASV2 and ASV5, unclassified Micrococcaceae ASV24, and Streptomyces: ASV27), decreased in response to BMc inoculation under extensive N-fertilization intensity (logistic regression, p < 0.05, Benjamini-Hochberg correction, n = 4, Fig. 6B). Moreover, the relative abundance of three out of 15 ASVs classified as Actinobacteriota (Agromyces: ASV3, Streptomyces: ASV25 and ASV27), decreased in response to BMc inoculation under extensive N-fertilization intensity (logistic regression, p < 0.05, Benjamini-Hochberg correction, n = 4, Fig. 6B). In contrast, one ASV (ASV20), classified as Bacillus sp., increased due to BMc treatment in samples under extensive N-fertilization (logistic regression, p < 0.05, Benjamini-Hochberg correction, n = 4; Fig. 6B). Notably, ASV20 exhibited high sequence identity with Bacillus atrophaeus and the ABi03 inoculated strain based on pairwise sequence alignment of the amplified region of the 16S rRNA gene (Fig. S6A). Furthermore, ASV20 was consistently among the most abundant ASVs detected in the rhizosphere across all samples. Given its high abundance in control samples (Fig. 6B) and similarity to other B. atrophaeus strains (Fig. S6A), we retained ASV20 in the dataset used to assess BMc-driven modulation of the microbiome. Unlike ABi03, we did not observe ASVs affiliated with Pseudomonas (especially RU47) among the dominant ASVs in the maize rhizosphere, a result that was expected based on our culture-dependent data (Fig. 1A). However, bacterial ASV1083 was phylogenetically associated with RU47, and its relative abundance tended to increase with BMc inoculation (logistic regression, p > 0.05, Benjamini-Hochberg correction; n = 4, data not shown).

Fungal community assembly in the rhizosphere was remarkably affected by BMc inoculation

The fungal community composition in the maize rhizosphere was significantly affected by BMc inoculation (PERMANOVA, R² = 0.33, p = 1x10^− 4, n = 8, Fig. 5). Fertilization intensity also influenced fungal community composition, albeit to a lesser extent (PERMANOVA, R² = 0.09, p = 0.028, n = 8, Fig. 5B). Accordingly, both experimental variables, particularly BMc inoculation, significantly affected the composition of the fungal community across all taxonomic levels in the rhizosphere (Fig. 6C and 6D). For example, at the phylum level, the relative abundance of Ascomycota, the most dominant phylum, increased from 68.2–80.7% in extensively fertilized samples due to BMc inoculation (logistic regression, p < 0.05, n = 4; Fig. 6B). Contrarily, Mortierellomycota and unclassified fungi (Fungi_phy_Incertae) significantly decreased in relative abundance following BMc inoculation under extensive N-fertilization (logistic regression; p < 0.05, n = 4; Fig. 6C).

Among the 1605 fungal ASVs identified, 121 and 124 exhibited increased relative abundance in response to BMc inoculation under intensive and extensive N-fertilization, respectively, while 49 and 110 ASVs displayed decreased relative abundance under intensive and extensive N-fertilization, respectively (Table S6). Notably, eight dominant ASVs responded to BMc inoculation in the extensive treatments, while only three ASVs responded in the intensive N-fertilization context. Specifically, four fungal ASVs, affiliated with the genera Trichoderma, Marquandomyces, Talaromyces, and Penicillium, were significantly enriched under extensive fertilization (logistic regression, p < 0.05, n = 4; Fig. 6D), while Podila, Mortierella, Stachybotrys and Fungi_gen_incertae_sedis showed an opposite trend (logistic regression, p < 0.05, n = 4; Fig. 6D). Interestingly, BMc inoculation also led to increased abundance of Trichoderma and Talaromyces in intensive fertilized treatments (Fig. 6D). Moreover, the shifts in abundance of the highly-abundant fungal taxa under intensive N-fertilization were restricted to ASV1530 (Pseudogymnoascus), which significantly decreased from the explicitly high abundance of 40.23 ± 19.89% to 8.51 ± 3.74% (logistic regression, p < 0.05, n = 4; Fig. 6D). However, none of the enriched Trichoderma ASVs in response to BMc were closely phylogenetically associated with OMG16, based on the pairwise alignment of the amplified ITS region (Fig. S6).

BMc inoculation enhanced relative abundance of microbial lipopeptide and siderophore biosynthesis genes in the maize rhizosphere

Analysis of the metagenomic dataset revealed numerous sequences annotated to OMG16 ABi03, and RU47, corroborating their establishment in the rhizosphere. BMc inoculation significantly increased the relative abundance of reads assigned to B. atrophaeus ABi03 (ANOVA, F = 96.0, p = 4x10^− 8, n = 4; Fig. 7A) and Pseudomonas sp. RU47 (ANOVA, F = 50.6, p = 12x10^− 4, n = 4; Fig. 7A), while no such increase was detected for Trichoderma harzianum OMG16. To investigate the impact of BMc inoculation on microbial functions in the rhizosphere, metagenomic shotgun sequencing reads were mapped against a customized database targeting potential plant-beneficial microbial functions (Table S7). BMc inoculation significantly influenced the functional composition of bacterial communities in the rhizosphere (PERMANOVA, R² = 0.15, p = 0.0014, Fig. 7B) and led to a significant increase in the abundance of genes associated with chemotaxis, quorum sensing, the degradation of aromatic compounds, as well as synthesis of auxin, siderophores, lipopeptides (e.g., surfactin, iturin, athrofactin) and spermidine (potA) (edgeR, p < 0.05, n = 4. Figure 7C). Genes such as mycobactin synthetase (mbtE), involved in the synthesis of iron-chelating siderophores, were highly enriched due to BMc inoculation, particularly under extensive N-fertilization (Fig. 7C). In contrast, lvH, dppA, rbsB, potD and 4,5-diphosphate decarboxylase (involved in type-3-secretion system, chemotaxis, spermidine production and aromatic compound degradation, respectively) decreased with BMc inoculation (edgeR, p < 0.05, n = 4; Fig. 7C).

In order to link microbial taxa to their respective functions, we taxonomically annotated sequences assigned to differentially abundant functions. Approximately 66% of sequences assigned to specific functions met the criteria for reliable taxonomic classification. Eight genes associated with the production of lipopeptides were linked to Bacillus spp., as discerned by sequence similarity. The majority of these genes exhibited 100% identity and 100% alignment with ABi03 (reads from each gene annotated as ABi03: fenD 63%, ituA 75%, ituB 67%, ituC 61%, srfAA 79%, srfAB 54%, srfAC 74%, and srfAD 93%; Fig. 7C). Furthermore, five BMc-enriched genes associated with siderophore production (dhbf 64%, entE 66%, entF 66%, mbtE 63%, and mbtF 59% reads, Table 2) were mapped to bacterial taxa affiliated with the phylum Actinobacteriota, with the highest proportion of attributed sequences annotated as Streptomyces spp. (Fig. S7). Notably, only a limited percentage of these genes was annotated to RU47 and ABi03 genomes, and these instances were exclusively restricted to the BMc samples (Table 2).

Table 2

Annotated siderophore genes in the shotgun sequencing data and their association with *Bacillus atrophaeus* ABi03 and *Pseudomonas* sp. RU47. In addition, the association of these genes to ABi03 and RU47 genomes was evaluated by DIAMOND (e = 10^− 5, identity = 100%).
Siderophore genes	Total counts	Number of counts annotated to genera	Number of genera annotated per gene	Read counts annotated as Pseudomonas	Read counts annotated as RU47 (100% identity)	Read counts annotated as Bacillus	Read counts annotated as ABi03 (100% identity)
2,3-dihydroxybenzoate-AMP-ligase, entE	647	428	185	4	0	48	27
L-serine ligase, entF	3389	2242	59	167	0	30	26
Glycine ligase, dhbF	4889	3138	176	257	1	191	141
Mycobactin peptide synthetase, mbtE	1146	721	98	78	2	4	2
Mycobactin peptide synthetase, mbtF	579	344	56	5	0	8	6

Improved plant micronutrient uptake, rhizosphere metabolite patterns and microbiome modulation are linked to BMc inoculation

To delve deeper into the interconnected shifts within the soil-plant system resulting from BMc application, we conducted an integrated network analysis. Only highly significant correlations (Pearson ρ ≥ |0.8|, p < 0.05, Benjamini-Hochberg correction, n = 16) and distance-based clustering were used (Fig. 8, Table S8). This analysis unveiled two primary clusters, Module 1 and Module 2, incorporating all variables positively or negatively correlated with maize growth, respectively (Fig. 8). Module 1 comprised ASV20, taxonomically linked to Bacillus ABi03 from the applied BMc (Fig. S6A), alongside various fungal ASVs taxonomically categorized as Trichoderma spp., although none were phylogenetically associated with Trichoderma OMG16 (Fig. S6B). All of these were positively correlated with growth-promoting variables (Fig. 8, Table S8). Interestingly, the concentrations of iron, zinc, and manganese, along with the functional genes encoding bacterial siderophore synthetases (mycobactin synthetase MbtE, Table S8) displayed a positive correlation with maize growth and the expression of NAS3, a plant gene involved in iron/zinc translocation (Fig. 1C). These positive correlations suggest that BMc inoculation improves iron uptake by facilitating alterations in bacterial siderophore release, and by stimulating root exudation of iron-mobilizing phenolics. This resulted in an improved iron-nutritional status, promoting ROS detoxification with protective effects under drought stress.

Intensified agricultural practices have revolutionized agricultural production enabling sustained high crop yields but concurrently compromising soil health and plant resilience^1,2,40. Therefore, implementation of sustainable farming practices is required, including environmentally sound methods for maintaining ecosystem functions. Recognizing the pivotal role played by rhizosphere microbiota in supporting plant performance, we hypothesized that the use of a beneficial microorganism consortium (BMc) mitigates growth reductions associated with reduced (extensive) N-fertilization intensity. Contrarily to our expectations, BMc inoculation outperformed plants grown under intensive N-fertilization, demonstrating the potential of BMc inoculation for promoting sustainable agriculture through reduced reliance on mineral N-fertilizer inputs.

The observed increase in maize shoot biomass and concomitant improvement of plant fitness in BMc inoculated plants were likely attributable to mitigate the impact of drought induced by BMc inoculation. This association may be explained by the prolonged period of unusually low rainfall experienced during the early plant development stages of maize in our field experiment (Fig. S1). Notably, no significant differences in shoot macronutrient concentrations were detected between control and BMc inoculated plants. However, the latter exhibited a significantly higher iron concentration compared to control plants, regardless of N-fertilization intensity. This observation is consistent with existing studies, which suggests a strong link between iron uptake and plant performance under drought stress conditions⁴¹. Indeed, recent research highlighted a notable reduction in the uptake of this micronutrient in Poaceae roots under drought conditions. For example, Xu et al. (2021)⁴¹ employed root RNA-seq data to demonstrate that drought-stressed sorghum exhibited downregulated iron transporter genes and upregulated iron storage functionality. Similar observations were found in barley under water limiting conditions⁴². These findings support the notion that plants respond to drought by downregulating iron uptake genes and upregulating iron storage^43,44. This strategic metabolic shift is thought to address the drought-induced reduction in photosynthesis, as the photosynthetic machinery is a major consumer of iron in plants^45,46. Additionally, it might reflect the decreased iron availability in dry soils and the negative impact of drought on iron acquisition mediated by mugineic acid phytosiderophores⁴¹. Accordingly, no mugineic acids were detected the rhizosphere soil solutions collected in our study (Fig. S4). Therefore, the prolonged drought during early maize growth could have limited iron uptake, consequently hindering plant development. Conversely, BMc-inoculated maize plants exhibited a significant increase in the relative abundance of genes associated with microbial siderophore production (specifically, dhbf, entE, entF, mbtE, and mbtF). This finding supports the presumption of enhanced iron availability in BMc inoculated plants, potentially contributing to the mitigation of negative effects of severe drought conditions. In contrast to the stimulation of shoot growth, root development and root AMF colonization remained unaffected by BMc inoculation (Fig. S2), suggesting that BMc inoculation promoted iron availability rather than spatial iron acquisition.

Interestingly, metagenomic analysis highlighted Actinobacteriota, particularly Streptomyces spp., as harboring genes for siderophore production. This finding aligns with recent research emphasizing the pivotal role of Actinobacteriota in mitigating drought stress in plants^41,47–49, particularly by secreting various iron-chelating organic molecules that sequester and immobilize soluble iron (Fe³⁺) for plant uptake in iron-deficient soils^50,51. More specifically, the genus Streptomyces, renowned for its capacity to produce diverse hydroxamate siderophores via divergent biosynthetic gene clusters⁵², has demonstrated an ability to attain a competitive growth advantage over other soil bacteria in environments characterized by limited iron availability^53,54 and during periods of drought⁵⁵. Network analysis further supported the potential benefits of Actinobacteriota, revealing a significant positive correlation between shoot dry mass and the presence of genes associated with Actinobacteriota-derived siderophores, underscoring their role in enhancing plant fitness. Control maize plants surprisingly also displayed a notable enrichment in Actinobacteriota, concomitant with a decrease in the relative abundance of taxa affiliated with the Bacteroidota and Proteobacteria phyla. This finding is consistent with observations commonly reported in moisture-limited soils, where an increase in the ratio of Gram-positive to Gram-negative bacteria is widely noted under drought conditions^56–58. Specifically, an increase in the abundance of Actinobacteriota, coupled with a concurrent reduction of Proteobacteria and Bacteroidotes in the plant rhizosphere under drought, has been frequently documented^59–63. The fact that BMc inoculated plants exhibited an opposite pattern, with a decrease in Actinobacteriota and an enrichment of Gram-negative taxa, could further substantiate the drought mitigation effect promoted by the inoculated BMc. However, the decrease in Actinobacteriota abundance, along with the increase in the relative abundance of Actinobacteriota-related siderophore genes, could also indicate a direct interaction between the inoculated microbes and the resident rhizosphere microbiota. This interaction may selectively enrich Actinobacteriota taxa possessing siderophore synthesis capabilities while depleting antagonistic members within the same phylum. Unfortunately, the employed short-read metagenomic sequencing approach failed to yield fully assembled genomes, hindering the estimation of genomic content for Actinobacteriota members positively associated with iron uptake. To address this limitation, further sequencing endeavors utilizing a combination of short and long-read and Hi-C sequencing technologies could potentially enhance metagenome assemblies. This approach holds promise in facilitating a more accurate assessment of the genomic potential of the identified Actinobacteriota members, particularly those linked to iron uptake.

The maize rhizosphere microbiota exhibited significant responses to both N-fertilization intensities and BMc inoculation. Notably, these factors impacted fungal community assembly more than bacterial communities. BMc inoculation specifically led to an enrichment of fungal amplicon sequence variants (ASVs) with potential plant growth-promoting capabilities. This included Talaromyces spp., previously documented to enhance plant growth under controlled drought conditions^64–66. Additionally, the rhizosphere of inoculated maize plants harbored a greater abundance of ASVs classified as Trichoderma spp., although not exhibiting close phylogenetic relationship to the inoculated OMG16. Several Trichoderma strains are known for producing phytohormones and siderophores^67–71. Therefore, it is plausible that these Trichoderma taxa contributed to the observed improvement in micronutrient uptake following BMc inoculation. It is worthwhile noting that the BMc strains used are commonly found in commercial biofertilizers and are known for their beneficial effects on plants. These strains possess documented beneficial traits, including modulation of plant hormones, antagonism against pathogens, and promotion of plant nutrition¹³. Among them, ABi03 displayed the highest density and relative abundance in the maize rhizosphere, indicating superior rhizo-competence compared to the other strains. Integrated network analysis identified ABi03 as a key driver of improved nutrient uptake, leading to increased shoot dry mass and concomitant alterations in plant gene transcription, hormonal levels, stress factors, and root exudate composition. However, the possibility that RU47 and OMG16 persisted in sufficient densities during initial inoculation stages and subsequently declined over time cannot be excluded.

The establishment of rhizosphere microorganisms is highly influenced by root exudates^72–74. To comprehensively assess the impact of BMc inoculation on the composition of root exudates and rhizosphere soil solution, root windows were installed. This methodology facilitated a nuanced understanding of the intricate interplay between BMc and the specific root exudates that shape the rhizosphere microbial community. Plants inoculated with BMc exhibited higher levels of organic acids, sugars, and amino acids in their rhizosphere solution, indicating enhanced physiological activity of roots and rhizosphere microbes. The increase in trehalose within the rhizosphere metabolite profile of BMc inoculated plants could be associated with BMc-induced adaptations of maize plants to drought conditions. Indeed, trehalose was extensively documented for its role in improving adaptation to water deficit in various crops, including maize^75–77. Additionally, BMc promoted rhizosphere accumulation of several secondary metabolites known to influence rhizosphere microbiota, including phenolic acids, flavonoids, and benzoxazinoids. The latter, a fascinating class of secondary metabolites mainly known from the Poaceae, exhibit a multifaceted role in plant physiology. Benzoxazinoids are renowned for their allelopathic effects, herbivore deterrence, and antimicrobial activity, significantly contributing to plant resilience^65–68. In addition to their role as defense mechanism activators, recent research has revealed that benzoxazinoids exhibit siderophore-like properties, demonstrating a high affinity to chelate metals, particularly iron ^78–80. This property also applies to phenolic acids and flavonoids, which accumulated in the rhizosphere in response to BMc inoculation^81,82. This siderophore capability has implications for nutrient acquisition, potentially influencing plant-microorganism interactions and nutrient cycling within the rhizosphere. The detected levels of these metabolites (10–15 µg cm^− 1 root length, Fig. 7) translate to rhizosphere concentrations in the millimolar range, which are sufficiently high to mediate the mobilization of mineral nutrients⁸³. Furthermore, benzoxazinoids can act as chemo-attractants for beneficial microorganisms, including Pseudomonas and Bacillus. The tolerance of these bacteria to benzoxazinoids may explain the observed advantage of ABi03 in achieving high densities in the maize rhizosphere. Bacterial tolerance to benzoxazinoids, particularly MBOA, likely influences the benzoxazinoid-dependent community structure of the maize rhizosphere microbiome^84–86. Moreover, recent studies have reported increased accumulation of jasmonic acid (JA) in oilseed rape inoculated with Bacillus spp. and Trichoderma OMG16³², similar to the increased JA levels observed in the root tissue of BMc inoculated maize plants in this study (Fig. S3). Interestingly, JA triggers the biosynthesis of benzoxazinoids⁸⁷. Therefore, the BMc may have induced the production and release of benzoxazinoids by stimulating JA accumulation. The intricate interplay between benzoxazinoids and siderophore dynamics underscores the multifaceted and integral role these compounds play within the broader ecological and physiological context. These observations also suggest a potential direct effect of BMc on maize metabolic patterns, leading to a shift in root exudate production and, consequently, shaping bacterial and fungal community compositions in the rhizosphere.

Our findings suggest that BMc inoculation potentially exerted systemic mitigation effects on the prolonged and severe drought stress experienced by the investigated maize plants during their early growth stages. Supporting this hypothesis, a significant decrease in the expression of various maize genes associated with metabolic stress was observed in the leaf tissue of BMc inoculated plants. Furthermore, a reduction in shoot stress indicators, including hydrogen peroxide accumulation, ABA, and JA levels, was accompanied by elevated gene expression and activities of key iron-dependent enzymes involved in hydrogen peroxide detoxification, such as ascorbate peroxidase and superoxide dismutase. Additionally, there was an increase in glycine betaine, an essential organic osmolyte that enhances plant responses to water stress. Apart from alleviating metabolic stress, BMc inoculation also led to increased concentrations of plant growth-promoting phytohormones (IAA and cytokinins). While the precise mechanisms underlying the observed alterations in plant hormone and metabolic stress profiles remain elusive, potentially involving enhanced micronutrient uptake or direct cross-signaling with the applied BMc, our results demonstrate that BMc inoculation promoted plant health and mitigated the negative impacts of drought conditions.

Our study demonstrates beneficial effects of BMc inoculation in enhancing maize growth under both conventional and reduced N-fertilization intensities under field conditions. We elucidate the intricate interplay between BMc strains, iron metabolism, drought conditions, the rhizosphere microbiome, and root exudate profiles. This highlights the efficacy of a multidisciplinary approach in unraveling these complex relationships. Furthermore, we demonstrate the potential application of BMc alongside reduced N-fertilization intensity to promote soil health and sustainable agriculture without compromising crop productivity. BMc inoculation not only modulates the rhizosphere microbiome but also concurrently alters root exudation patterns, thereby enhancing iron uptake and overall plant health. Our observations reveal that BMc inoculation significantly influences the assembly of maize rhizosphere microbiota while simultaneously increasing the abundance of gene related with microbial siderophores and the root exudation of benzoxazinoids, which are potential iron chelators. These findings hold particular relevance for alleviating iron deficiency, a common consequence of drought conditions in the rhizosphere. This indicates the potential of BMc to mitigate drought effects in long-term applications and foster sustainable agriculture in regions suffering from extended drought periods.

Experimental design at long-term field trial

An experiment was conducted within an ongoing long-term field trial (LTE Westerfeld, https://doi.org/10.20387/bonares-9ne9-a828) in Bernburg, Germany (51° 82’ N, 11° 70’ E) to assess the impact of a beneficial microorganisms consortium (BMc: Pseudomonas sp. RU47 (DSM177411), Bacillus atrophaeus ABi03 (DSMZ 32285), and Trichoderma harzianum OMG16 (DSMZ 32722) on maize performance (cv. Benedictio, KWS Saat SE & Co. KGaA) during the 2020 growing season. The LTE, with the soil type loess chernozem over limestone (8% sand, 70% silt, 22% clay; pH 7.0-7.4) was described in detail by Sommermann et al. (2018)⁸⁸. The treatments were arranged within single field plots (0.6 ha, subdivided into four subplots/replicates) managed with cultivator tillage (12–15 cm depth), which were further subdivided into strips with reduced nitrogen (N) fertilization intensity (40 kg ha^− 1 N) without fungicide use (extensive; Ext) compared to conventional N-fertilization intensity (100 kg ha^− 1 N, intensive; Int) with recommended pesticide (including fungicide) application. Within the respective replicate subplots (size 2.25 m²), controls (Ctrl) and maize BMc treated plants were arranged at a distance of 2 m. Each treatment (Ctrl-Ext, BMc-Ext, Ctrl-Int, and BMc-Int) included four replicates (or subplots) with 33 maize plants each, with an intra-row distance of 14 cm between plants and inter-row distance of 60 cm.

Inoculum preparation and BMc application

An endospore suspension of a rifampicin-resistant B. atrophaeus ABi03 selectant was provided by ABiTEP GmbH (Berlin, Germany). The inoculum of T. harzianum OMG16 (strain collection of Anhalt University of Applied Sciences, Bernburg, Germany) was prepared as described by Hafiz et al. (2022)³². The inoculum of the rifampicin-resistant Pseudomonas sp. RU47 selectant (DSMZ117411) was grown in nutrient broth (Sifin diagnostics GmbH, Germany) supplemented with rifampicin (75 µg ml^− 1; Th. Geyer GmbH & Co. KG, Germany) on a rotary shaker (200 rpm) at 28°C for 24 h to obtain fresh viable cells. RU47, ABi03 and OMG16 were mixed in tap water to a final concentration of 10⁸ CFU ml^− 1 per strain immediately before drenching the plants. Each maize plant was drenched two and five weeks after emergence (EC 12 and EC14) with 50 ml BMc-suspension. The control plants were drenched with 50 ml tap water at the respective inoculation times.

Installation of root windows

To study effects of the experimental variables on root exudation and rhizosphere metabolite profiles, root windows were installed in each control and BMc treated subplots/replicates prior to the second drenching, according to Neumann et al. (2009)⁸⁹. Specifically, soil profiles measuring 50x50 cm were excavated adjacent to maize plants at EC 14–16 using steel plates. After removing the plates, the profiles were covered with plexiglass sheets stabilized with wooden timbers inserted into the soil. Subsequently, the root windows were thermally insulated with two layers of Styrofoam® panels and finally covered with black plastic foil and soil to prevent light and rain exposure (Fig. S4).

Sampling and sample processing

For molecular and physiological analyses, a sample (3x2 cm from the middle portion of the leaf lamina) from the youngest fully developed leaf was collected at the flowering stage (EC 53–63, three maize plants per treatment replicate). The leaf material from each treatment was pooled and immediately immersed in a total of 10 mL of RNAlater (Thermo Fisher Scientific, Darmstadt, Germany). The samples were incubated at 4°C overnight and then stored at -80°C. After grinding the leaf material in liquid nitrogen, total RNA was extracted from 100 mg of the homogenized sample using the RNeasy Plant Mini Kit (Qiagen GmbH, Hilden, Germany) and quantified spectrophotometrically. Single-stranded cDNA synthesis from 2 µg of total RNA was performed with the High Capacity cDNA Reverse Transcription Kit with RNase Inhibitor (Applied Biosystems, Foster City, CA, USA), and qPCR was performed as previously described⁹⁰. The targeted genes and corresponding primer pairs are listed in Table S3. An aliquot of the leaf material preserved in RNAlater was used for the analysis of plant hormones and stress-associated metabolites and enzymes.

Following leaf sampling, these maize plants were excavated, the shoots were separated from the roots, and soil samples were collected. Shoots, roots and soil samples were stored immediately at 4°C. Maize shoots were dried at 60°C until and shoot dry mass (SDM) was measured. Shoot and soil nutrient content was analyzed as described by Behr et al. (2023)¹³.

Sampling of maize rhizosphere for microbial analysis

Rhizosphere soil, defined as soil loosely adhering to roots, was obtained by shaking roots and subsequently sieving through a 2 mm mesh before being stored at − 20°C. The rhizosphere competence of the individual strains of BMc was assessed via serial dilution plating counts following the method described by Behr et al. (2023)¹³. Total community DNA was extracted from the collected rhizosphere soil using the FastPrep-24 bead-beating system and FastDNA Spin Kit for Soil (MP Biomedicals, Santa Ana, CA, USA) following the manufacturer’s recommendations. DNA samples were further purified with the GeneClean Spin Kit (MP Biomedicals, Santa Ana, CA, USA).

Analysis of root traits and mycorrhization

The WinRHIZO root analysis system (Regent Instruments, Quebec, Canada) was used to measure total root length (TRL) of the collected roots by optical scanning. Root hair length was determined after magnification of high-resolution digital photographs taken from roots growing along the observation windows (Zeiss Axiovision software, Oberkochen, Germany). For analysis of mycorrhization, a modified staining method following Vierheilig et al. (1998) was used. Collected roots were washed and cut into segments of 1–2 cm and incubated at 90◦C for 45 min in 10% (w/v) KOH. Bleached roots were acidified with 1–3% HCl and then stained in a 5% ink-vinegar solution [5% (v/v) ink in 5% (v/v) acetic acid] at 90◦C for 10 min. Destaining was performed with acidified tap water. The grid–line intersection method of Giovanetti and Mosse (1980)⁹¹ was employed to determine the rate of mycorrhization.

Sampling rhizosphere solution from root window plots for metabolite profiling

The collection of rhizosphere solutions from the root windows was conducted on the same day as the collection of soil and plant materials described previously, by following the procedure outlined by Neumann et al., (2014)⁹² (Fig. S4). After root window opening, moist filter discs (5 mm diameter) were placed onto the surface of first order lateral roots in subapical root zones (1–2 cm behind the root tip), mature root zones (8–10 cm behind the root tip), crown roots, and soil areas without visible roots. Two filters were applied for each root zone in five replicates, covering root segments of at least 1 cm per replicate. After a sampling period of four hours, the filters were removed and pooled in 2 mL of 80% v/v methanol solution. The sorption filters were precipitated by centrifugation for extraction, and the supernatants were filtered (Chromafil R O-20/15MS) and stored at -80°C.

Amplicon sequencing of bacterial and fungal rhizosphere communities

For prokaryotic community analyses, DNA from the maize rhizosphere was used for amplification of 16S rRNA gene fragments using the 16S primers Uni341F (5’-CCTAYGGGRBGCASCAG-3’) and Uni806R (5’- GGACTACNNGGGTATCTAAT − 3’) targeting Bacteria and Archaea⁹³. Library construction and sequencing was carried out by Novogene (Cambridge, UK) on NovaSeq 6000 PE250. Primers and adapters were removed from sequences using cutadapt⁹⁴ (v3.7). Paired-end reads were processed using the DADA2 pipeline⁹⁵ (v1.26.0) using R⁹⁶ (v.4.2.2).. The obtained amplicon sequence variants (ASVs) were taxonomically classified to the lowest possible taxonomic level by using a Naive Bayesian Classifier⁹⁷, trained on the SILVA SSU Reference Taxonomy database⁹⁸ (v138.1). Sequences originating from chloroplasts or mitochondria, and sequences with less than five reads were removed. Following these steps, no sequence classified as Archaea was present in the final dataset. A total of 1.35x10⁶ out of 4.31x10⁶ reads (31.8 ± 5.2%) were retained after processing, generating 6595 bacterial ASVs in total with an average of 42114 reads per sample.

In parallel, we investigated the fungal community in the rhizosphere by performing sequencing of the Internal Transcribed Spacer (ITS2) region. Briefly, three PCRs per sample were carried out at three annealing temperatures (56°C ± 2°C) with 24 cycles in 20 µl volume using Q5® High-Fidelity 2x Master Mix (New England Biolabs). Each PCR included 10 ng TC-DNA and bovine serum albumin (final concentration 0.5 mg ml^− 1), and PCRs were performed using the standard 8-nucleotide Illumina barcoded ITS86F/ITS4 primer pair for PCR amplification⁸⁸. High-throughput amplicon sequencing was carried out on the Illumina® MiSeq® platform in paired-end mode (2x 300 bp) as previously described⁹⁹. ASV generation and taxonomic assignment based on database-dependent strategy¹⁰⁰ using the GALAXY platform and UNITE v9.0 database^101,102 was performed as described in Behr et al. (2024)¹⁰³. A total of approximately 3.43x10⁶ out of 3.82x10⁶ reads (89.8%±1.5%) were retained, and after singleton removal the total number of fungal ITS2-ASVs was 1965 (average of 107251 reads per sample). Rarefaction curves were generated for estimating the read coverage of each sample (Fig. S8).

Library preparation and shotgun sequencing for metagenome analysis

Shotgun metagenomic libraries for functional analysis of rhizosphere DNA samples were prepared using NEBNext Ultra II FS DNA Library Prep Kit® (New England Biolabs, Frankfurt, Germany) according to the manufacturer’s protocol with modifications as described previously¹⁰³. The libraries were diluted to 1 nM and sequenced on a NextSeq550 sequencer (Illumina, San Diego, CA, USA) using the NextSeq 500/550 High Output Kit v2.5 (300 Cycles) after equimolar pooling. The shotgun metagenomic sequencing generated 186,817,794 forward and 186,817,794 reverse high-quality raw reads with a sequencing depth of 16x10⁶ reads per sample. Trimming and merging resulted in a total of 170,500,215 merged sequences (library size range 7,168,494 − 18,104,610) representing the majority of reads (91.3%). Nonpareil output showed coverage between 5.4% − 13.9%, which indicated a sufficient coverage for read-based analyses (Fig. S9). Taxonomic and functional annotation were conducted following the methodology outlined in our prior study¹⁰³. For detection of the applied BM, sequences taxonomically assigned to Trichoderma harzianum, Bacillus atrophaeus and Pseudomonas sp. RU47 were extracted. In order to confirm whether these sequences actually represent the applied BM, sequences were compared with the respective BM genomes using NCBI BLAST (blastn¹⁰⁴) with 100% sequence identity and 100% sequence alignment. Additionally, sequences of differentially abundant genes were compared with the genomes of the BM to estimate their potential contribution to the respective function, such as iturin, fengycin, surfactin and genes associated with siderophore production.

Determination of plant hormones and stress metabolites

The plant hormones (indole acetic acid, gibberellic acid, cytokinin, abscisic acid, jasmonic acid and salicylic acid) were analyzed from shock-frozen maize shoot and root material. UHPLC-MS analysis of phytohormones was carried out on a Velos LTQ System (Thermo Fisher Scientific, Waltham, MA, United States) utilizing a Synergi Polar column, 4µ, 150×3.0 mm, (Phenomenex, Torrance, CA, United States) as described by Moradtalab et al. (2020)³⁴. Part of the same material, was also used for determination of selected stress metabolites via spectrophotometry (Spectrophotometer, Hitachi, Tokyo, Japan). Total phenolics were determined, after methanol extraction (80% v/v methanol) at 750 nm, using the Folin method¹⁰⁵. Proline analysis was conducted at 520 nm after acetic acid and acid ninhydrin derivatization¹⁰⁶. Glycine betaine determination was performed according to the method of Grieve and Grattan (1983)¹⁰⁷ with modifications described by Valadez-Bustos et al. (2016)¹⁰⁸. The 1,1-diphenyl-2-picrylhydrazyl radical (DPPH) method was used to evaluate the free radical scavenging activity of antioxidants in the plant tissue³⁴. Hydrogen peroxide levels were determined at 390 nm as described previously¹⁰⁶. Ascorbate peroxidase (APX, EC 1.11.1.11) activity was recorded according to the method described by Boominathan and Doran (2002)⁹⁸. Extraction and determination of superoxide dismutase (SOD, EC 1.15.1.1) was optimized for shoot tissue of maize according to the method of Moradtalab et al. (2018)¹⁰⁶.

High performance liquid chromatography profiles of rhizosphere soil solutions

HPLC-profiling of organic acids, sugars and amino acids in the rhizosphere soil solutions in 80% methanol extracts of the root windows sorption filters was conducted as described by Windisch et al. (2021)⁹⁰. For phenolic compounds and benzoxazinoids, identification was performed with positive/negative switching LC-MS on a QExactive Plus Electrospray Mass Spectrometer (Thermo Fisher Scientific) coupled to an Agilent 1290 Ultra Performance Liquid Chromatography System. Details of the HPLC conditions are presented in Table S9. A list of investigated phenolics is reported in Table S4.

Data analysis and statistics

Statistical analysis of SDM and nutrient data was performed using R⁹⁶ (v.4.2.2). Main and interaction effects between fertilization intensity and BMc inoculation were analyzed by two-way ANOVA for SDM, shoot/soil nutrients and BMc rhizosphere competence. Homoscedasticity and normal distribution of residuals was inspected visually using the performance package¹⁰⁹ (v.0.10.2). Pairwise comparison was conducted using the Tukey‘s HSD test of the agricolae package¹¹⁰ (v.1.3.5). Data were visualized using the ggplot2¹¹¹ (v.3.4.1) and ggpattern packages¹¹² (v.1.0.1). For gene expression data we performed t-test (Ctrl vs BMc).

For handling the sequencing data, we mainly used the “tidyverse” set of packages¹¹³ (v.1.3.1⁠) and vegan¹¹⁴ (v.2.6.1⁠). Samples were analysed in relative abundance (% of reads) for sequencing data. Bray-Curtis distance was used for estimating the β-diversity with log₁₀ transformed datasets and pseudocount addition (only amplicon sequencing data). PERMANOVA tests were applied to evaluate how the N-fertilization intensity or BMc inoculation affected bacterial and fungal β-diversity by using the “vegan” package.

For investigating the effect of BMc inoculation on different bacterial and fungal taxa we performed logistic regression and p-values were corrected with Benjamini-Hochberg. The differential abundance testing focused on the rhizosphere and was performed separately for intensive and extensive fertilization. The tests were applied for several taxonomic groups, starting from the Phylum and reaching up to the ASV level. For metagenomic shotgun sequencing data, we performed differential abundance testing via edgeR¹¹⁵ (v2.1.2⁠) with maximum likelihood test and p-values corrected with Benjamini-Hochberg method. For each differential abundance test the minimum cut-off was 75%. Comparisons with p-values below 0.05 were considered as statistically significant (α = 0.05). To estimate whether an ASV corresponds to the 16S rRNA gene or ITS from our BMc strains we performed phylogenetic association of responder ASVs to BMc inoculation with the 16S rRNA gene or ITS2 DNA sequence from our BMc strains. In addition, we also downloaded 16S rRNA gene and ITS2 sequences from the NCBI type strains for comparison. We performed multiple sequencing alignment (packages: “DECIPHER”¹¹⁶, v2.26.0).

An integrated network analysis was performed with the variables (e.g. SDM, plant gene expression, bacterial/fungal ASVs) that passed the differential testing. To create the network, we performed multiple Pearson correlation (p-values were corrected with Benjamini-Hochberg method). We included only correlation with ρ>|0.8| and p < 0.05. Clustering was carried out by using the positive correlations and the Louvain method to divide the network into different modules (package “igraph”, v1.5.1¹¹⁷). The network was visualized with the packages “igraph” and “ggraph”¹¹⁸ (v2.1).

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors have read the final version of the manuscript and have given their explicit consent for publication.

Competing interests

The authors declare no competing interests.

Funding

This work was supported by the German Federal Ministry of Education and Research (BMBF, Germany) within the framework of the project DiControl (http://dicontrol.igzev.de/de/; grant number 031B0514A- E) as part of the BonaRes funding initiative ‘Soil as a sustainable resource for the bioeconomy’ (https://www.bonares.de).

Author Contribution

D.F. (Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft), I.K. (Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft), T.K.N. (Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft), D.B. (Data curation, Formal analysis, Investigation, Methodology, Project administration, Writing – review & editing), L.S. (Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization), J.H.B. (Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization), S.P.C. (Data curation, Formal analysis, Investigation, Methodology), R.Z. (Formal analysis, Investigation, Methodology), N.G. (Data curation, Formal analysis, Investigation, Methodology), M.S. (Supervision, Writing – review & editing), J.G. (Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Writing – review & editing), U.W. (Project administration, Supervision, Writing – review & editing), G.N. (Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Validation, Writing – review & editing), K.S. (Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing), and R.G. (Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Validation, Writing – review & editing).

Acknowledgement

We thank Peggy Gabor, Jasmin Schmidt, Sabine Grabner and Maryna Lakishyk for excellent technical assistance.

Data Availability

Sequencing data (amplicon and shotgun) were deposited at the Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra) under the BioProject accession PRJNA1045550. Fungal ITS2 sequences can be found under PRJEB74508. All scripts including detailed parameters and data to reproduce the figures for this study have been deposited in github: https://github.com/JonKampouris/BMc_inoculants_on_promoting_plant_growth

Emmerson M, et al. Chapter Two - How Agricultural Intensification Affects Biodiversity and Ecosystem Services. In: Dumbrell AJ, Kordas RL, Woodward G, editors. Advances in Ecological Research. Volume 55. Academic; 2016. pp. 43–97.
Liu Y, Pan X, Li J. Current Agricultural Practices Threaten Future Global Food Production. J Agric Environ Ethics. 2015;28:203–16.
Hoekstra AY, Wiedmann TO. Humanity’s unsustainable environmental footprint. Science. 2014;344:1114–7.
Crowder DW, Northfield TD, Strand MR, Snyder WE. Organic agriculture promotes evenness and natural pest control. Nature. 2010;466:109–12.
Tuck SL, et al. Land-use intensity and the effects of organic farming on biodiversity: a hierarchical meta-analysis. J Appl Ecol. 2014;51:746–55.
MacLaren C, et al. Long-term evidence for ecological intensification as a pathway to sustainable agriculture. Nat Sustain. 2022;5:770–9.
Jat ML, et al. Conservation agriculture for sustainable intensification in South Asia. Nat Sustain. 2020;3:336–43.
Seufert V, Ramankutty N, Foley JA. Comparing the yields of organic and conventional agriculture. Nature. 2012;485:229–32.
Hart MM, Antunes PM, Chaudhary VB, Abbott LK. Fungal inoculants in the field: Is the reward greater than the risk? Funct Ecol. 2018;32:126–35.
Meier MS, et al. Environmental impacts of organic and conventional agricultural products – Are the differences captured by life cycle assessment? J Environ Manage. 2015;149:193–208.
Lutz S, et al. Soil microbiome indicators can predict crop growth response to large-scale inoculation with arbuscular mycorrhizal fungi. Nat Microbiol. 2023;8:2277–89.
Zhou Y, et al. Crop rotation and native microbiome inoculation restore soil capacity to suppress a root disease. Nat Commun. 2023;14:8126.
Behr JH et al. Beneficial microbial consortium improves winter rye performance by modulating bacterial communities in the rhizosphere and enhancing plant nutrient acquisition. Front Plant Sci 14, (2023).
Zhou X, et al. Cross-kingdom synthetic microbiota supports tomato suppression of Fusarium wilt disease. Nat Commun. 2022;13:7890.
Wang M, et al. Dynamic root microbiome sustains soybean productivity under unbalanced fertilization. Nat Commun. 2024;15:1668.
Trivedi P, Batista BD, Bazany KE, Singh BK. Plant–microbiome interactions under a changing world: responses, consequences and perspectives. New Phytol. 2022;234:1951–9.
Busby PE, et al. Research priorities for harnessing plant microbiomes in sustainable agriculture. PLOS Biol. 2017;15:e2001793.
Qiu Z, Egidi E, Liu H, Kaur S, Singh BK. New frontiers in agriculture productivity: Optimised microbial inoculants and in situ microbiome engineering. Biotechnol Adv. 2019;37:107371.
Schmitz L, et al. Synthetic bacterial community derived from a desert rhizosphere confers salt stress resilience to tomato in the presence of a soil microbiome. ISME J. 2022;16:1907–20.
Porter SS, et al. Beneficial microbes ameliorate abiotic and biotic sources of stress on plants. Funct Ecol. 2020;34:2075–86.
Vorholt JA, Vogel C, Carlström CI, Müller DB. Establishing Causality: Opportunities of Synthetic Communities for Plant Microbiome Research. Cell Host Microbe. 2017;22:142–55.
Emmenegger B, et al. Identifying microbiota community patterns important for plant protection using synthetic communities and machine learning. Nat Commun. 2023;14:7983.
French E, Kaplan I, Iyer-Pascuzzi A, Nakatsu CH, Enders L. Emerging strategies for precision microbiome management in diverse agroecosystems. Nat Plants. 2021;7:256–67.
Haskett TL, Tkacz A, Poole PS. Engineering rhizobacteria for sustainable agriculture. ISME J. 2021;15:949–64.
Saad MM, Eida AA, Hirt H. Tailoring plant-associated microbial inoculants in agriculture: a roadmap for successful application. J Exp Bot. 2020;71:3878–901.
Khare E, Arora NK. Effects of Soil Environment on Field Efficacy of Microbial Inoculants. in Plant Microbes Symbiosis: Applied Facets (ed. Arora, N. K.) 353–381Springer India, (2015). 10.1007/978-81-322-2068-8_19
Chouhan GK, et al. Phytomicrobiome for promoting sustainable agriculture and food security: Opportunities, challenges, and solutions. Microbiol Res. 2021;248:126763.
Sun X, et al. Bacillus velezensis stimulates resident rhizosphere Pseudomonas stutzeri for plant health through metabolic interactions. ISME J. 2022;16:774–87.
Hafiz FB, et al. Synergistic Effects of a Root-Endophytic Trichoderma Fungus and Bacillus on Early Root Colonization and Defense Activation against Verticillium longisporum in Rapeseed. Mol Plant Microbe Interact. 2022;35:380–92.
Schreiter S, Sandmann M, Smalla K, Grosch R. Soil Type Dependent Rhizosphere Competence and Biocontrol of Two Bacterial Inoculant Strains and Their Effects on the Rhizosphere Microbial Community of Field-Grown Lettuce. PLoS ONE. 2014;9:e103726.
Windisch S, et al. Rhizoctonia solani and bacterial inoculants stimulate root exudation of antifungal compounds in lettuce in a soil-type specific manner. Agronomy. 2017;7:1–17.
Hafiz FB, et al. Tissue-Specific Hormone Signalling and Defence Gene Induction in an In Vitro Assembly of the Rapeseed Verticillium Pathosystem. Int J Mol Sci. 2023;24:10489.
Mpanga IK et al. The form of n supply determines plant growth promotion by p-solubilizing microorganisms in maize. Microorganisms 7, (2019).
Moradtalab N et al. Synergisms of Microbial Consortia, N Forms, and Micronutrients Alleviate Oxidative Damage and Stimulate Hormonal Cold Stress Adaptations in Maize. Front Plant Sci 11, (2020).
Woo SL, Pepe O. Microbial consortia: Promising probiotics as plant biostimulants for sustainable agriculture. Front Plant Sci. 2018;9:428205.
Zaim S, Bekkar AA, Belabid L. Efficacy of Bacillus subtilis and Trichoderma harzianum combination on chickpea Fusarium wilt caused by F. oxysporum f. sp. ciceris. Arch Phytopathol Plant Prot. 2018;51:217–26.
Karuppiah V, Sun J, Li T, Vallikkannu M, Chen J. Co-cultivation of Trichoderma asperellum GDFS1009 and Bacillus amyloliquefaciens 1841 Causes Differential Gene Expression and Improvement in the Wheat Growth and Biocontrol Activity. Front Microbiol. 2019;10:441938.
Edlinger A, et al. Agricultural management and pesticide use reduce the functioning of beneficial plant symbionts. Nat Ecol Evol. 2022;6:1145–54.
Hage-Ahmed K, Rosner K, Steinkellner S. Arbuscular mycorrhizal fungi and their response to pesticides. Pest Manag Sci. 2019;75:583–90.
Evenson RE, Gollin D. Assessing the Impact of the Green Revolution, 1960 to 2000. Science. 2003;300:758–62.
Xu L, et al. Genome-resolved metagenomics reveals role of iron metabolism in drought-induced rhizosphere microbiome dynamics. Nat Commun. 2021;12:3209.
Cantalapiedra CP et al. Large Differences in Gene Expression Responses to Drought and Heat Stress between Elite Barley Cultivar Scarlett and a Spanish Landrace. Front Plant Sci 8, (2017).
Tripathi DK et al. Acquisition and Homeostasis of Iron in Higher Plants and Their Probable Role in Abiotic Stress Tolerance. Front Environ Sci 5, (2018).
D’Oria A, et al. Drought specifically downregulates mineral nutrition: Plant ionomic content and associated gene expression. Plant Direct. 2022;6:e402.
Schmidt SB, Eisenhut M, Schneider A. Chloroplast Transition Metal Regulation for Efficient Photosynthesis. Trends Plant Sci. 2020;25:817–28.
Müller B. Iron transport mechanisms and their evolution focusing on chloroplasts. J Plant Physiol. 2023;288:154059.
Metze D, et al. Microbial growth under drought is confined to distinct taxa and modified by potential future climate conditions. Nat Commun. 2023;14:5895.
Selim S, Hassan YM, Saleh AM, Habeeb TH, AbdElgawad H. Actinobacterium isolated from a semi-arid environment improves the drought tolerance in maize (Zea mays L). Plant Physiol Biochem. 2019;142:15–21.
Martins BR, et al. Cultivar-specific dynamics: unravelling rhizosphere microbiome responses to water deficit stress in potato cultivars. BMC Microbiol. 2023;23:377.
Shivlata L, Tulasi S. Thermophilic and alkaliphilic Actinobacteria: biology and potential applications. Front Microbiol 6, (2015).
Thoenen L et al. Bacterial tolerance to host-exuded specialized metabolites structures the maize root microbiome. Proc. Natl. Acad. Sci. U. S. A. 120, e2310134120.
Cruz-Morales P, et al. Actinobacteria phylogenomics, selective isolation from an iron oligotrophic environment and siderophore functional characterization, unveil new desferrioxamine traits. FEMS Microbiol Ecol. 2017;93:fix086.
Xu L et al. Drought delays development of the sorghum root microbiome and enriches for monoderm bacteria. Proc. Natl. Acad. Sci. 115, E4284–E4293 (2018).
Fitzpatrick CR et al. Assembly and ecological function of the root microbiome across angiosperm plant species. Proc. Natl. Acad. Sci. 115, E1157–E1165 (2018).
Jones SE et al. Streptomyces Volatile Compounds Influence Exploration and Microbial Community Dynamics by Altering Iron Availability. mBio 10, 10.1128/mbio.00171 – 19 (2019).
Acosta-Martínez V, et al. Predominant bacterial and fungal assemblages in agricultural soils during a record drought/heat wave and linkages to enzyme activities of biogeochemical cycling. Appl Soil Ecol. 2014;84:69–82.
Fuchslueger L, Bahn M, Fritz K, Hasibeder R, Richter A. Experimental drought reduces the transfer of recently fixed plant carbon to soil microbes and alters the bacterial community composition in a mountain meadow. New Phytol. 2014;201:916–27.
Chodak M, Gołębiewski M, Morawska-Płoskonka J, Kuduk K, Niklińska M. Soil chemical properties affect the reaction of forest soil bacteria to drought and rewetting stress. Ann Microbiol. 2015;65:1627–37.
Naylor D, DeGraaf S, Purdom E, Coleman-Derr D. Drought and host selection influence bacterial community dynamics in the grass root microbiome. ISME J. 2017;11:2691–704.
Taketani RG, et al. Dry Season Constrains Bacterial Phylogenetic Diversity in a Semi-Arid Rhizosphere System. Microb Ecol. 2017;73:153–61.
Bouskill NJ, et al. Pre-exposure to drought increases the resistance of tropical forest soil bacterial communities to extended drought. ISME J. 2013;7:384–94.
Barnard RL, Osborne CA, Firestone MK. Responses of soil bacterial and fungal communities to extreme desiccation and rewetting. ISME J. 2013;7:2229–41.
Naylor D, Coleman-Derr D. Drought Stress and Root-Associated Bacterial Communities. Front Plant Sci 8, (2018).
Halo BA, Al-Yahyai RA, Al-Sadi AM. An endophytic Talaromyces omanensis enhances reproductive, physiological and anatomical characteristics of drought-stressed tomato. J Plant Physiol. 2020;249:153163.
Kaur R, Saxena S. Evaluation of drought-tolerant endophytic fungus Talaromyces purpureogenus as a bioinoculant for wheat seedlings under normal and drought-stressed circumstances. Folia Microbiol (Praha). 2023;68:781–99.
Ou T, et al. Study on the Potential for Stimulating Mulberry Growth and Drought Tolerance of Plant Growth-Promoting Fungi. Int J Mol Sci. 2023;24:4090.
Macías-Rodríguez L, Contreras-Cornejo HA, Adame-Garnica SG, del-Val E, Larsen J. The interactions of Trichoderma at multiple trophic levels: inter-kingdom communication. Microbiol Res. 2020;240:126552.
Contreras-Cornejo HA, Macías-Rodríguez L, del-Val E, Larsen J. Interactions of Trichoderma with Plants, Insects, and Plant Pathogen Microorganisms: Chemical and Molecular Bases. in Co-Evolution of Secondary Metabolites (eds. Merillon, J.-M. & Ramawat, K. G.) 1–28Springer International Publishing, (2018). 10.1007/978-3-319-76887-8_23-1
Woo SL, Hermosa R, Lorito M, Monte E. Trichoderma: a multipurpose, plant-beneficial microorganism for eco-sustainable agriculture. Nat Rev Microbiol. 2023;21:312–26.
Abdenaceur R, et al. Effective biofertilizer Trichoderma spp. isolates with enzymatic activity and metabolites enhancing plant growth. Int Microbiol. 2022;25:817–29.
Vinale F et al. Trichoderma Secondary Metabolites Active on Plants and Fungal Pathogens. Open Mycol J 8, (2014).
Chaparro JM, Badri DV, Vivanco JM. Rhizosphere microbiome assemblage is affected by plant development. ISME J. 2014;8:790–803.
Jing H, Wang H, Wang G, Liu G, Cheng Y. The mechanism effects of root exudate on microbial community of rhizosphere soil of tree, shrub, and grass in forest ecosystem under N deposition. ISME Commun. 2023;3:1–13.
Haichar F et al. el Z. Plant host habitat and root exudates shape soil bacterial community structure. ISME J. 2, 1221–1230 (2008).
Nuccio ML, et al. Expression of trehalose-6-phosphate phosphatase in maize ears improves yield in well-watered and drought conditions. Nat Biotechnol. 2015;33:862–9.
Acosta-Pérez P, et al. Characterization of Trehalose-6-phosphate Synthase and Trehalose-6-phosphate Phosphatase Genes and Analysis of its Differential Expression in Maize (Zea mays) Seedlings under Drought Stress. Plants. 2020;9:315.
Ali Q, Ashraf M. Induction of Drought Tolerance in Maize (Zea mays L.) due to Exogenous Application of Trehalose: Growth, Photosynthesis, Water Relations and Oxidative Defence Mechanism. J Agron Crop Sci. 2011;197:258–71.
Zhou S, Richter A, Jander G. Beyond Defense: Multiple Functions of Benzoxazinoids in Maize Metabolism. Plant Cell Physiol. 2018;59:1528–37.
Hu L, et al. Plant iron acquisition strategy exploited by an insect herbivore. Science. 2018;361:694–7.
Bigler L, Baumeler A, Werner C, Hesse M. Detection of Noncovalent Complexes of Hydroxamic-Acid Derivatives by means of electrospray mass spectrometry. Helv Chim Acta. 1996;79:1701–9.
Wang L, et al. Multifaceted roles of flavonoids mediating plant-microbe interactions. Microbiome. 2022;10:233.
Neumann G, Ludewig U. Rhizosphere chemistry influencing plant nutrition. in Marschner’s Mineral Nutrition of Plants (Fourth Edition) (eds. Rengel, Z., Cakmak, I. & White, P. J.) 545–585Academic Press, (2023). 10.1016/B978-0-12-819773-8.00013-7
Bar-Yosef B. Root excretions and their environmental effects. Influence on availability of phosphorus. Plant Roots. The Hidden Half. Marcel Dekker, Inc.; 1991.
Cadot S, et al. Specific and conserved patterns of microbiota-structuring by maize benzoxazinoids in the field. Microbiome. 2021;9:103.
Hu L, et al. Root exudate metabolites drive plant-soil feedbacks on growth and defense by shaping the rhizosphere microbiota. Nat Commun. 2018;9:2738.
Gfeller V et al. Plant secondary metabolite-dependent plant-soil feedbacks can improve crop yield in the field. eLife 12, e84988 (2023).
Sue M, Fujii M, Fujimaki T. Increased benzoxazinoid (Bx) levels in wheat seedlings via jasmonic acid treatment and etiolation and their effects on Bx genes including Bx6. Biochem. Biophys. Rep. 27, 101059 (2021).
Sommermann L, et al. Fungal community profiles in agricultural soils of a long-term field trial under different tillage, fertilization and crop rotation conditions analyzed by high-throughput ITS-amplicon sequencing. PLoS ONE. 2018;13:e0195345.
Neumann G, George TS, Plassard C. Strategies and methods for studying the rhizosphere—the plant science toolbox. Plant Soil. 2009;321:431–56.
Windisch S et al. Impact of Long-Term Organic and Mineral Fertilization on Rhizosphere Metabolites, Root–Microbial Interactions and Plant Health of Lettuce. Front Microbiol 11, (2021).
Giovannetti M, Mosse B. An Evaluation of Techniques for Measuring Vesicular Arbuscular Mycorrhizal Infection in Roots. New Phytol. 1980;84:489–500.
Neumann G et al. Root exudation and root development of lettuce (Lactuca sativa L. cv. Tizian) as affected by different soils. Front Microbiol 5, (2014).
Sundberg C, et al. 454 pyrosequencing analyses of bacterial and archaeal richness in 21 full-scale biogas digesters. FEMS Microbiol Ecol. 2013;85:612–26.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10–2.
Callahan BJ, et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.
R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing; 2020.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naïve Bayesian Classifier for Rapid Assignment of rRNA Sequences into the New Bacterial Taxonomy. Appl Environ Microbiol. 2007;73:5261–7.
Quast C, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.
Babin D, et al. Distinct rhizomicrobiota assemblages and plant performance in lettuce grown in soils with different agricultural management histories. FEMS Microbiol Ecol. 2021;97:fiab027.
Antweiler K, et al. Statistical test for tolerability of effects of an antifungal biocontrol strain on fungal communities in three arable soils. Microb Biotechnol. 2017;10:434–49.
Kõljalg U, et al. UNITE: a database providing web-based methods for the molecular identification of ectomycorrhizal fungi. New Phytol. 2005;166:1063–8.
Abarenkov K, et al. The UNITE database for molecular identification and taxonomic communication of fungi and other eukaryotes: sequences, taxa and classifications reconsidered. Nucleic Acids Res. 2024;52:D791–7.
Behr JH, et al. Long-term conservation tillage with reduced nitrogen fertilization intensity can improve winter wheat health via positive plant–microorganism feedback in the rhizosphere. FEMS Microbiol Ecol. 2024;100:fiae003.
McGinnis S, Madden TL. BLAST: at the core of a powerful and diverse set of sequence analysis tools. Nucleic Acids Res. 2004;32:W20–5.
Swain T, Hillis WE. The phenolic constituents of Prunus domestica. I.—The quantitative analysis of phenolic constituents. J Sci Food Agric. 1959;10:63–8.
Moradtalab N et al. Silicon Improves Chilling Tolerance During Early Growth of Maize by Effects on Micronutrient Homeostasis and Hormonal Balances. Front Plant Sci 9, (2018).
Grieve CM, Grattan SR. Rapid assay for determination of water soluble quaternary ammonium compounds. Plant Soil. 1983;70:303–7.
Valadez-Bustos MG, et al. A reliable method for spectrophotometric determination of glycine betaine in cell suspension and other systems. Anal Biochem. 2016;498:47–52.
Lüdecke D, et al. see: An R Package for Visualizing Statistical Models. J Open Source Softw. 2021;6:3393.
Felipe de Mendiburu & Muhammad Yaseen. agricolae: Statistical Procedures for Agricultural Research. (2020).
Wickham H. ggplot2: Elegant Graphics for Data Analysis. New York: Springer-; 2016.
FC M, Davis TL. ggpattern: ‘ggplot2’ Pattern Geoms. (2022).
Wickham H, et al. Welcome to the Tidyverse. J Open Source Softw. 2019;4:1686.
Oksanen J et al. vegan: Community Ecology Package. (2008).
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
Wright ES, Using. DECIPHER v2.0 to Analyze Big Biological Sequence Data in R. R J. 2016;8:352–9.
Csardi G, Nepusz T. The igraph software package for complex network research. InterJournal Complex Syst, 1695 (2006).
Pedersen TL. ggraph: An Implementation of Grammar of Graphics for Graphs and Networks. (2024).

No competing interests reported.

MS3SITablesfinal.docx

Download PDF

Reviewers agreed at journal
14 Nov, 2024
Reviews received at journal
21 Oct, 2024
Reviewers agreed at journal
29 Sep, 2024
Reviewers invited by journal
25 Sep, 2024
Editor assigned by journal
05 Sep, 2024
Submission checks completed at journal
29 Aug, 2024
First submitted to journal
28 Aug, 2024

You are reading this latest preprint version

Microbial inoculants modulate the rhizosphere microbiome, alleviate plant stress responses and enhance maize growth at field scale

Status:

Version 1

Abstract

Background

Results

Conclusions

Figures

Background

Results

Description of the experimental design

BMc strains colonized the maize rhizosphere and increased shoot dry mass and iron uptake

The patterns of rhizosphere metabolites are shaped by BMc inoculation

BMc inoculation increased Myxococcota, Bacteroidota, and Proteobacteria in the maize rhizosphere

Fungal community assembly in the rhizosphere was remarkably affected by BMc inoculation

Improved plant micronutrient uptake, rhizosphere metabolite patterns and microbiome modulation are linked to BMc inoculation

Discussion

Conclusions

Methods

Experimental design at long-term field trial

Inoculum preparation and BMc application

Installation of root windows

Sampling and sample processing

Sampling of maize rhizosphere for microbial analysis

Analysis of root traits and mycorrhization

Sampling rhizosphere solution from root window plots for metabolite profiling

Amplicon sequencing of bacterial and fungal rhizosphere communities

Library preparation and shotgun sequencing for metagenome analysis

Determination of plant hormones and stress metabolites

High performance liquid chromatography profiles of rhizosphere soil solutions

Data analysis and statistics

Declarations

Funding

Author Contribution

Acknowledgement

Data Availability

References

Additional Declarations

Supplementary Files

Status:

Version 1