Beneficial microbial consortium improves winter rye performance by modulating bacterial communities in the rhizosphere and enhancing plant nutrient acquisition

The beneficial effect of microbial consortium application on plants is strongly affected by soil conditions, which are influenced by farming practices. The establishment of microbial inoculants in the rhizosphere is a prerequisite for successful plant-microorganism interactions. This study investigated whether a consortium of beneficial microorganisms establishes in the rhizosphere of a winter crop during the vegetation period, including the winter growing season. In addition, we aimed for a better understanding of its effect on plant performance under different farming practices. Winter rye plants grown in a long-time field trial under conventional or organic farming practices were inoculated after plant emergence in autumn with a microbial consortium containing Pseudomonas sp. (RU47), Bacillus atrophaeus (ABi03) and Trichoderma harzianum (OMG16). The density of the microbial inoculants in the rhizosphere and root-associated soil was quantified in autumn and the following spring. Furthermore, the influence of the consortium on plant performance and on the rhizosphere bacterial community assembly was investigated using a multidisciplinary approach. Selective plating showed a high colonization density of individual microorganisms of the consortium in the rhizosphere and root-associated soil of winter rye throughout its early growth cycle. 16S rRNA gene amplicon sequencing showed that the farming practice affected mainly the rhizosphere bacterial communities in autumn and spring. However, the microbial consortium inoculated altered also the bacterial community composition at each sampling time point, especially at the beginning of the new growing season in spring. Inoculation of winter rye with the microbial consortium significantly improved the plant nutrient status and performance especially under organic farming. In summary, the microbial consortium showed sufficient efficacy throughout vegetation dormancy when inoculated in autumn and contributed to better plant performance, indicating the potential of microbe-based solutions in organic farming where nutrient availability is limited.

The beneficial effect of microbial consortium application on plants is strongly affected by soil conditions, which are influenced by farming practices.The establishment of microbial inoculants in the rhizosphere is a prerequisite for successful plant-microorganism interactions.This study investigated whether a consortium of beneficial microorganisms establishes in the rhizosphere of a winter crop during the vegetation period, including the winter growing season.In addition, we aimed for a better understanding of its effect on plant performance under different farming practices.Winter rye plants grown in a long-time field trial under conventional or organic farming practices were inoculated after plant emergence in autumn with a microbial consortium containing Pseudomonas sp.(RU47), Bacillus atrophaeus (ABi03) and Trichoderma harzianum (OMG16).The density of the microbial inoculants in the rhizosphere and root-associated soil was quantified in autumn and the following spring.Furthermore, the influence of the consortium on plant performance and on the rhizosphere bacterial community assembly was investigated using a multidisciplinary approach.Selective plating showed a high colonization density of individual microorganisms of the consortium in the rhizosphere and root-associated soil of winter rye throughout its early growth cycle.16S rRNA gene amplicon sequencing showed that the farming practice affected mainly the rhizosphere bacterial communities in autumn and spring.However, the microbial consortium inoculated altered also the bacterial community composition at each sampling time point, especially at the beginning of the new growing season in spring.Inoculation of winter rye with the microbial consortium significantly improved

Introduction
Soil microorganisms play a central role in almost all soil processes, as they do not only contribute to soil formation (Watteau and Villemin, 2018) but also have a significant impact on important ecosystem functions such as nutrient cycling and plant performance (Latz et al., 2016;Lori et al., 2017).Therefore, managing the soil microbial community is essential to maintain soil health and productivity.Microbial diversity has been positively associated with the resilience of microbial communities and soil fertility (van Elsas et al., 2012;Usero et al., 2021).Crops benefit from diverse soil microbial communities as they recruit and enrich beneficial microorganisms (BM) via root exudates in the rhizosphere (Berg and Smalla, 2009;Liu et al., 2021).BM can stimulate plant growth through various mechanisms, such as modulating phytohormone levels (Zubair et al., 2019) and increasing nutrient availability (Jacoby et al., 2017).Furthermore, beneficial rhizosphere microorganisms can trigger induced systemic resistance and protect the plant against below-and aboveground pathogens (Wei et al., 2020;Vlot et al., 2021).Indeed, plant pathogens can be controlled by BM either directly by producing lytic enzymes and bioactive compounds (Fira et al., 2018;Ali et al., 2020) or indirectly by competing for resources in the same ecological niche like the rhizosphere (Rai et al., 2016).
The inoculation of microorganisms with known beneficial functions for plants is a practical approach that is able to modulate the composition of indigenous microbial communities, particularly in the rhizosphere (Deng et al., 2019;Berg et al., 2021).Although many BM show promising functions on plant growth and protection under controlled conditions in vitro and in vivo (Mendes et al., 2013;Mazzola and Freilich, 2017), their efficacy in the field is often limited due to insufficient colonization of the host rhizosphere and/or unfavorable conditions (Parnell et al., 2016;Batista and Singh, 2021).To successfully colonize the rhizosphere, inoculated BM have to bypass the barrier of the indigenous soil microbial community, adapt to variable environmental conditions, and interact with the host plant using mainly root exudates as a chemoattractant and substrate (Finkel et al., 2017;Berg et al., 2021).Different plant developmental stages and changing environmental conditions shift the root exudation pattern and thereby also the composition of the microbial community in the rhizosphere (Chaparro et al., 2014;Windisch et al., 2021), which may affect the rhizosphere competence of BM.Thus, even if inoculated BM successfully establish in the rhizosphere at early plant development, persistence over an extended period is not guaranteed.Inoculation with multi-strain consortia containing two or more BM has shown increased efficacy compared to single BM species (Sun et al., 2022).Consortia of BM, comprising members with diverse functions, can occupy different ecological niches, making them more resilient to variable environments (Bradaćǒváet al., 2019;Tosi et al., 2020).Furthermore, the interaction among consortium members increases their ecological fitness, and complementary beneficial functions generate synergies that promote their interaction with the plant host (Moradtalab et al., 2020;Pascale et al., 2020).
Different farming practices, such as conventional and organic farming, drastically affect the structure of microbial communities (Francioli et al., 2016;Windisch et al., 2021).Chowdhury et al. (2019) showed that conventional farming using synthetic pesticides and mineral fertilizers shaped the rhizosphere microbial community composition differently from organic farming, which was associated with varying states of plant health.There is currently a limited understanding of how the farming practice impacts the efficacy of a beneficial microorganism consortium (BMc).For instance, crops can benefit more from increased nutrient bioavailability through BMc in nutrient-deficient soils (Eltlbany et al., 2019), while in soils with high nutrient level the rhizosphere colonization by the inoculated BMc might be reduced (Lopes et al., 2021).
The main goal of this study was to assess the rhizosphere competence of the inoculated BMc members in winter rye (Secale cereale cv.Conduct) during the growing season in dependence on different farming practices.Winter cereal crops are exposed to extreme environmental changes during the growing season, raising the question of whether the members of a BMc inoculated in autumn can establish in the rhizosphere of winter rye and maintain a high colonization density over the winter period.Furthermore, the impact of the applied BMc on plant performance may be defined by the physicochemical properties of the soil and the farming practice dependent soil microbial communities.In this context, we hypothesized that (i) early inoculation of winter rye enables sufficient colonization of each BMc member at early plant developmental stage, supporting its persistence in the rhizosphere throughout the vegetation period; (ii) the application of BMc shapes the composition of the rhizosphere bacterial community depending on the farming practice and thus differentially affects the plant performance.
To test our hypotheses, we first characterized in vitro the plant growth-promoting traits and the ability of each BMc member to inhibit the growth of soil-borne pathogens.Then, we used a longtime field experiment to evaluate the rhizosphere competence and efficacy of the tested BMc on plant performance under field conditions with different farming practices (organic vs. conventional farming practice).Winter rye plants were inoculated shortly after emergence and BMc colonization and plant performance were evaluated in autumn and spring.The dynamics of the bacterial community throughout vegetation dormancy were characterized by 16S rRNA amplicon sequencing.

In vitro characterization of consortium members
To provide a more comprehensive characterization of the plantbeneficial traits associated with each consortium member (Pseudomonas sp.RU47, B. atrophaeus ABi03, and T. harzianum OMG16) a series of in vitro tests were conducted to assess their antagonistic capabilities and their capacity for plant growth promotion.The modulation of plant growth hormones by indole-3-acetic acid (IAA) production and 1-aminocyclopropane-1carboxylate (ACC) deaminase activity (Koo et al., 2010), as well as the increased nutrient availability by siderophore production (Schwyn and Neilands, 1987) was tested for each strain.In addition, the solubilization of potassium feldspar (KAlSi 3 O 8 ) (Breitkreuz et al., 2021), calcium phosphate (Ca 3 (PO 4 ) 2 ) (Nautiyal, 1999), zinc oxide (ZnO) (Mumtaz et al., 2017), manganese dioxide (MnO 2 ) (Sanket et al., 2017) and ammonia production (Cappuccino and Sherman, 2014) by the individual strains was analyzed.For nutrient solubilization assays, the medium was modified for OMG16 using Waksman-agar (Huber et al., 1987).
Cellulase (Berg et al., 2005), chitinase (Berg et al., 2001), b-1,3glucanase (Weinert et al., 2010), protease (Berg et al., 2005), and hydrogen cyanide (HCN) (Donate-Correa et al., 2005) formation were tested as antagonistic functions.In vitro tests for enzyme secretion with antagonistic functions were not feasible for OMG16 due to cultivation incompatibilities.A dual culture assay modified after Lin and Ho (2021) was used to evaluate the ability of each strain to inhibit pathogen growth in vitro.For the two bacterial strains RU47 and ABi03, 20 µl of bacterial suspension (10 5 colony forming units (CFU) mL -1 ) was dispensed in a centered line of a nutrientagar I plate (Sifin diagnostics, Berlin, Germany).Freshly grown mycelia plugs (ø 6 mm) of Fusarium culmorum (Isolate F247), F. graminearum (Isolate FG66), and Rhizoctonia solani Isolate AG3 (Isolate Ben3) were placed at a distance of 2.5 cm on both sides of the line.For the fungal strain OMG16, precultured mycelia plugs of OMG16 and the pathogen were placed on opposite sides of a potato dextrose agar (PDA, Carl Roth, Karlsruhe, Germany) plate.Dual culture assays were cultivated in the dark at 25°C for 10 days.Antagonistic properties, i.e. formation of an inhibition zone, of the individual BMc strains were checked every second day (Supplementary Figure 1).

Field site description
The experiment was conducted from 2020-2021 in an experimental field with plots under conventional and organic farming practices as part of the long-time field trial in Thyrow (Thy_ABS; 52°15' N, 13°14' E, 44 m a.s.l.), which was established in 2005 at the Agricultural Research Institute of the Humboldt University of Berlin.Pallid brown earth (Luvisols) of sand over deep loam with low humus contents is the predominant soil type (83% sand, 14% silt, and only 3% clay) at this site.The average annual temperature (1981-2020) is 9.2°C with a mean annual precipitation of 509.8 mm.
Winter rapeseed (Brassica napus ssp.napus), winter wheat (Triticum aestivum), and winter rye (Secale cereale) were rotated in conventional farming and lucerne (Medicago sativa), lucerne, potato (Solanum tuberosum), triticale (Triticale), forage pea (Pisum sativum ssp.arvense), spring barley (Hordeum vulgare), and winter rye (Secale cereale) in organic farming.Mineral fertilization was used in conventional farming practice according to the agricultural standard of the region (N = 120 kg ha −1 , P = 21 kg ha −1 , K = 120 kg ha −1 , Mg = 12 kg ha −1 , S = 14 kg ha −1 ).Nitrogen was applied as calcium ammonium nitrate (KAS) and potassium, magnesium and sulfur as equal ratios of Patentkali ® (K + S Minerals and Agriculture, Kassel, Germany) and Korn-Kali ® (K + S Minerals and Agriculture).Triple superphosphate (46% P 2 O 5 ) was used for phosphate fertilization.Before growing winter rye and triticale, cattle manure was incorporated into the soil in organic farming (15,000 kg ha -1 ), supported by the cultivation of lucerne and forage peas.Plant protection was conducted following the agricultural standards of conventional cultivation [herbicides: 2.0 L ha -1 Picona (BASF, Ludwigshafen, Germany), 0.25 L ha -1 Cadou SC (Bayer, Leverkusen, Germany); fungicides: 2.0 L ha -1 Ceriax (BASF)].In organic farming, no pesticides were used, and emerging weeds were removed with a mechanical cultivator.Before sowing of winter rye (sowing density: 200 seeds per m 2 in conventional farming and 300 seeds per m 2 in organic farming), the soil was tilled, including stubble clearing with a disc harrow, plowing to a depth of about 20-25 cm and seedbed preparation with a reciprocating harrow.

Preparation of the microbial consortium and its field application
The T. harzianum OMG16 strain was cultivated on PDA medium (Carl Roth) supplemented with 100 mg L -1 penicillin (Carl Roth), 50 mg L -1 streptomycin sulfate (Sigma-Aldrich, St. Louis, USA), and 10 mg L -1 tetracycline hydrochloride (AppliChem, Darmstadt, Germany).To prepare the field inoculum, millet grains (Panicum miliaceum) were soaked overnight in tap water.The grains were washed intensively with cold tap water afterwards.Clean grains (500 g per bag) were transferred to bags designed for fungal culture (sun bag; Sigma-Aldrich) and autoclaved on three subsequent days.A disc with freshly grown OMG16 mycelium was added to the grains and incubated for 30 days in the dark at room temperature.After OMG16 entirely colonized the millet grains, the inoculum was freeze-dried and ground.Shortly before winter rye sowing, the powder (100 mg m -2 ) was incorporated into the soil to a depth of 20 cm.Freeze-dried millet grain powder without OMG16 was used as a control.
The B. atrophaeus ABi03 and Pseudomonas sp.RU47 strains were used as spontaneous rifampicin-resistant mutants provided as a prepared inoculum by ABiTEP (Schreiter et al., 2018).Two weeks after winter rye emergence, each plot was drenched with 4 L of the bacterial ABi03 and RU47 suspension, mixed with tap water (each 7.5 x 10 7 CFU mL -1 ).Control plots were drenched with tap water.Each treatment included four replicates (1.5 x 2 m) arranged in a randomized block design for both farming practices.

Sampling and verification of rhizosphere competence
Only plants from the inner plot core (0.5 x 1m) of each replicate were considered for the sampling, which was performed in autumn at seven weeks post inoculation (WPI) at EC 21, and in spring of the following year (22 WPI) at EC 25-29.A total of 20 winter rye plants were sampled per plot and treated as one sample per treatment replicate.After plant excavation, the shoots were separated from the roots.To obtain root-associated soil, the soil fraction loosely adhering to the root system was shaken off.Roots, shoots, and root-associated soil were stored immediately at 4°C.After shoot fresh mass determination, shoots were dried by lyophilization, and shoot dry mass (SDM) was measured.
Roots with the remaining tightly attached soil were washed with sterile tap water and cut into 1 cm pieces for rhizosphere sampling.After pooling the root fragments, five grams of root material were transferred into a sterile stomacher bag containing 15 mL of sterile 0.3% NaCl.The rhizosphere samples were obtained by a Stomacher 400 Circulator (Seward Ltd., Worthing, UK) for 60 s at 300 rpm.The supernatant was transferred into a falcon tube, and the Stomacher blending steps were repeated twice.The combined supernatants was used for determining bacterial BMc counts by immediately plating serial dilutions onto nutrient agar I plates (Sifin diagnostics) containing 75 µg mL -1 rifampicin (Serva Electrophoresis, Heidelberg, Germany) and 100 µg mL -1 cycloheximide (Serva Electrophoresis) and incubated in the dark for 48 h at 28°C.Given the distinct colony morphology of ABi03 and RU47, the CFU were counted and CFU per g root dry mass, used for stomacher processing, were calculated.The remaining Stomacher supernatant was further processed as described by Schreiter et al. (2014) for bacterial rhizosphere community analysis based on total community DNA extraction.
Root-associated soil was used to determine OMG16 counts.Therefore, 5 g of root-associated soil were extracted in a total volume of 50 mL NaCl (0.3%) by shaking for 30 min at room temperature on an orbital shaker (HSM-10; Hettich, Tuttlingen, Germany) at 70 rpm.Subsequently, the supernatant was immediately plated on Trichoderma selection medium as suggested by Williams et al. (2003) containing 250 mg L -1 chloramphenicol (Carl Roth), 90 mg L -1 streptomycin (Sigma-Aldrich), 200 mg L -1 quintozene (Sigma-Aldrich), 926 mg L -1 propamocarb (ProPlant; Arysta LifeScience, Paris, France) and 150 mg L -1 rose bengal (AppliChem).Following an incubation of OMG16 at 28°C in the dark for 10 days, the CFU were quantified per gram soil dry mass.Soil dry mass was obtained by drying five grams of fresh soil at 110°C until a constant weight was reached.
Winter rye ears of the whole plot were sampled at full ripeness.After counting the ears, the grain was threshed and the weight was determined.Yield per hectare was extrapolated for better comparison with literature data.

Nutrient analysis
The nutrient analysis of shoot material was conducted according to the certified protocols of the Association of German Agricultural Analytic and Research Institutes, VDLUFA, Germany (Verband Deutscher Landwirtschaftlicher Untersuchungs-und Forschungsanstalten, 2011).200-500 mg of dry plant material were solubilized in 5 mL HNO 3 (65%) and 3 mL H 2 O 2 (30%) by microwave digestion at 210°C for 25 min.The extract was filtered and diluted to a final volume of 100 mL.K, P, Mg, S, Ca, Mn, Cu, and Zn concentrations were determined via inductively coupled plasma optical emission spectrometry (ICP-OES, Thermo Fisher Scientific, Dreieich, Germany); total C and N were determined via elemental analysis (Elementary Vario El cube; Langenselbold, Germany).The content of macro-and micronutrients in the rootassociated soil were analyzed according to the certified protocols of VDLUFA by AGROLAB (Leinefelde-Worbis, Germany).

DNA extraction and high-throughput amplicon sequencing of the 16S rRNA gene
The total community DNA was extracted from total rhizosphere pellets using the FastPrep-24 bead-beating system and FastDNA Spin Kit for Soil (MP Biomedicals, Santa Ana, USA) following the manufacturer's recommendations.Samples were purified with the GeneClean Spin Kit (MP Biomedicals).Library construction and sequencing of the 16S rRNA gene was carried out by Novogene (Cambridge, UK) on NovaSeq PE250 using the primers Uni341F (5'-CCTAYGGGRBGCASCAG-3') and Uni806R (5'-GGACTACNNGGGTATCTAAT-3') (Sundberg et al., 2013).Primers and adapters were removed with the software Cutadapt (Martin, 2011).Paired-end reads were processed using the DADA2 pipeline (Callahan et al., 2016).For taxonomic assignment of the obtained amplicon sequence variants (ASVs), the representative sequences were taxonomically classified down to the lowest possible taxonomic level in a Galaxy workflow (Cock et al., 2013) with an e-value cut-off of 0.001 and a percent identity cut-off of 80% against the SILVA 138.1 SSU Ref NR99 database (v.138.1, Quast et al., 2012;v138.1).Sequences identified as plastids (chloroplast or mitochondria) and sequences with less than five reads were removed.No ASV classified as archaeal taxa was present in our dataset.Rarefaction curves were generated to estimate the read coverage, and sequencing depth sufficiently covered the diversity of each sample.Amplicon sequence data has been deposited at the Sequence Read Archive (https:// www.ncbi.nlm.nih.gov/sra)under the BioProject accession number PRJNA975889.

Data analysis
Statistical analysis of winter rye growth, yield, nutrient concentration, and BMc abundance was performed using R (v.4.2.2, R Core Team, 2022).The main and interaction effects between longtime farming practices and BMc use were analyzed by two-way ANOVA for SDM and shoot nutrients.Differences in the BMc abundance in the rhizosphere and root-associated soil were tested by two-way ANOVA analyzing the main and interaction effects of the sampling time and different long-time farming practices.All analyses were inspected visually for homoscedasticity and normal distribution of residuals using the Durbin-Watson test from the "car" package (v.3.1.1,Fox and Weisberg, 2019) and normal Q-Q plots.If ANOVA assumptions failed, data were transformed using the R package "rcompanion " (v.2.4.26, Mangiafico, 2023) according to Tukey's ladder of power approach.The statistical evaluation was conducted using the Sidak-test algorithm and the pairwise comparison combined with a compact letter display of the packages "emmeans " (v.1.8.3, Lenth, 2022), "multcomp" (v.1.4, Hothorn et al., 2008) and "multcompView" (v.0.1.8,Graves et al., 2019).Data were visualized using the "ggplot2" (v.3.4.1, Wickham, 2016) and "ggpattern" packages (v.1.0.1, FC and Davis, 2022).Comparisons with p-values below 0.05 were considered statistically significant.
Bray-Curtis distance was used for estimating the b-diversity by normalizing the rarefied datasets to relative abundance (%) and applying log 10 transformation with pseudo-count addition.PERMANOVA tests were applied to evaluate how farming practice and BMc inoculation affected the bacterial community composition.A pairwise PERMANOVA test with Benjamini-Hochberg correction was performed after combining the two treatments into four groups using the "pairwiseadonis" package (v.0.4,Martinez, 2022).Significant differences in a-diversity (ASV Richness, Simpson, Shannon Index, and Evenness) were estimated with two-way ANOVA tests.If the data failed to fulfill the normality criteria based on the Shapiro test, the non-parametric aligned-rank ANOVA was performed ("ARTool" package, v.0.11.1, Kay et al., 2021).
For investigating the effect of farming practice type and BMc inoculation on the abundance of ASVs, ANCOM-BC2 (v2.1.2,Lin and Peddada, 2020) was performed as a differential abundance test with pvalue correction via Benjamini-Hochberg method.Following the differential abundance test, we performed multiple logistic regressions with Benjamini-Hochberg correction as post hoc tests.Logistic regression was selected because of its similar predictive power to machine learning algorithms (e.g., random forest) and its inherent interpretability (Topcuoglu et al., 2020).We used a logistic regression equation to predict the probability of the treatment to occur (i.e., Organic-BMc, Organic-Control, Conventional-BMc, or Conventional-Control) based on the relative abundance of the ASVs.Consequently, the logistic regression tests enabled us to evaluate whether the differentially abundant ASVs responded to a specific treatment (ASV responders).For each particular treatment, the fitted model's positive or negative model coefficient indicated higher or lower relative abundance, respectively.ASVs were assigned as positive or negative responders to each treatment based on the model coefficient.
To evaluate whether the biomass yield of rye was associated with higher proportions of the inoculated bacterial strains in comparison to the rest of the bacterial community members in the different treatments, ASV responder sequences classified as Pseudomonas spp. or Bacillus spp.were aligned to herein constructed databases.For Pseudomonas spp.we created a database that contained the 16S rRNA gene sequences from the six 16S rRNA gene copies of the inoculant RU47 (Kuzmanovićet al., 2018).For the strain ABi03, we performed cloning and sequencing of the complete 16S rRNA region using primers U8-27 and R1494-1514U8-27 and R1494- (Heuer et al., 2009)), where we identified six different 16S rRNA genes.In addition, we retrieved sequences of Pseudomonas spp. or Bacillus spp.from GenBank 16S rRNA gene restricted to the type strain collection (Federhen, 2015).The V3-V4 region of the 16S rRNA gene, used for Illumina amplicon sequencing, was extracted from the database sequences using the Cutadapt software (Martin, 2011).Pairwise alignment was performed for Bacillus spp.and Pseudomonas spp., using the ASVs and the constructed databases with the package "msa" (v.1.30.1, Bodenhofer et al., 2015).A distance matrix of the alignment was calculated via the package "seqinr" (v.4.2.23, Charif and Lobry, 2007).A phylogenetic distance was calculated with the neighbor-joining algorithm (package "ape", v.5.7,Paradis and Schliep, 2019).The phylogenetic distance was visualized as principal coordinates analysis.

In vitro traits of the microbial consortium members
In vitro, the individual members of the used microbial consortium were tested for antagonistic and growth-promoting functions (Table 1).ABi03 showed cellulase, chitinase, b-1,3glucanase, and protease activity but no HCN production.In dual culture assays, ABi03 reduced the growth of all tested pathogens by inhibiting the mycelium growth (Supplementary Figure 1).Despite the inability of RU47 to inhibit the growth of F. culmorum or F. graminearum, it is noteworthy that the RU47 colonies were not overgrown by these pathogens.Additionally, we could confirm a minor inhibition of R. solani growth by RU47 as well as protease and HCN production.OMG16 also reduced the growth of F. culmorum and F. graminearum by spreading faster than the pathogen on PDA.No lytic degradation of the mycelium was observed in the border region between OMG16 and the tested pathogens (Supplementary Figure 1).
RU47 dissolved K-feldspar and Ca 3 (PO 4 ) 2 in vitro in nutrient availability assays and showed siderophore as well as IAA production (Table 1).For ABi03, ammonia and siderophores were detected, but no increase in solubility of the tested nutrients or plant hormone modulation was observed.Furthermore, IAA, low siderophore production, and ZnO solubility could be shown for OMG16.

Density of the consortium members in the rhizosphere and soil
To evaluate the persistence of the individual BM in the rhizosphere and root-associated soil of winter rye at field scale, we estimated their density after inoculation with two subsequent samplings in autumn and spring of the following year.Both bacterial strains ABi03 and RU47 sufficiently colonized the rhizosphere of winter rye with a density of more than six Log 10 CFU per gram root dry mass seven WPI at both farming practices (Figure 1).The fungal strain OMG16 also showed a sufficient density in the root-associated soil of more than five Log 10 CFU per gram soil dry mass.No differences in the density of the individual BM in the respective habitat depending on farming practice were observed.At the second sampling in spring (22 WPI), no changes in the density of ABi03 in the rhizosphere of winter rye were found independent of farming practice (Figure 1, Supplementary Table 1).RU47 density decreased to approximately five Log 10 CFU per gram of root dry mass at both farming practices compared to the first sampling in autumn.Plants grown under organic farming practice showed a significantly reduced density of OMG16 in the root-associated soil compared to the first sampling in autumn.The two different farming practices had no significant impact on the ability of the individual BM to establish in the rye rhizosphere or root-associated soil at both sampling time points.

Plant performance
The mean SDM of winter rye plants at the first sampling (EC 21) in autumn ranged from 0.07 to 0.10 g plant -1 , while the mean SDM in the second sampling (EC 25-29) in spring was between 0.26 and 0.48 g plant -1 (Figure 2A, Supplementary Table 2A).Control plants without BMc inoculation showed no significant difference in SDM in response to farming practice (conventional vs. organic) at both sampling time points.Organically cultivated winter rye plants treated with the BMc showed a significant increase in SDM compared to those without BMc application in autumn (7 WPI, +20%) and particularly in spring (22 WPI, +45%).In contrast to organic farming, the BMc treatment did not affect the SDM of plants grown under conventional farming at both sampling time points.Two-way ANOVA further revealed a significant interaction between farming practice and BMc inoculation on SDM at both sampling time points (Supplementary Table 2B) indicating farming practice-dependent effects of BMc on SDM.At harvest, there were no significant differences in yield between organic and conventional farming, although organically grown plants had slightly lower grain weight per hectare (Figure 2B).The BMc treatment had no significant effect on the yield in either conventional or organic farming.However, a slight tendency towards higher yield per unit area was observed in organic farming after the BMc treatment.
To determine whether the nutrient status of winter rye was associated with farming practice and BMc inoculation, the concentrations of several micro-and macronutrients were measured in plants grown across all studied plots and sampling time points (Table 2).At the first sampling time in autumn, most nutrients were at the same level for both farming practices, regardless of BMc inoculation.Two-way ANOVA revealed an overall significant effect (p ≤ 0.021) of BMc inoculation on the nitrogen concentration (Supplementary Table 3), which was slightly increased with BMc inoculation for both farming practices.However, a pairwise comparison could not verify a significant difference between the four treatments.In contrast to the autumn sampling, the concentrations of almost all nutrients except carbon (C), copper (Cu), calcium (Ca), and manganese (Mn) were reduced in organic farming practice (Table 2) at the second sampling in spring of the following year.Based on a two-way ANOVA, the BMc treatment showed a significant effect on the majority of nutrients in both farming practices, with the exception of zinc (Zn) (Supplementary Table 3).However, a pairwise comparison across all four treatments revealed that the BMc treatment significantly increased only the phosphorus (P) concentration in organically grown plants and the sulphur (S) concentration in conventionally grown plants, in comparison to the control (Table 2).The soil properties were not affected by BMc treatment, but were mainly influenced by the farming practice (Supplementary Table 4).

Bacterial community in the rhizosphere
To elucidate whether the choice of farming practice and BMc inoculation influenced the composition of the winter rye rhizosphere bacterial communities, we performed 16S rRNA gene amplicon sequencing.The farming practice strongly influenced the composition of rhizosphere bacterial communities in the autumn (PERMANOVA, R 2 = 0.29, p < 0.001, Figure 3A) and in the spring sampling (PERMANOVA, R 2 = 0.23, p < 0.001, Figure 3B).Similarly, BMc inoculation significantly influenced the bacterial community composition of the rhizosphere in both samplings (PERMANOVA, R 2 = 0.10-0.13,p < 0.001, Figure 3).However, we did not detect any interaction effect between the choice of farming practice and the BMc inoculation (PERMANOVA, R 2 = 0.04-0.05).Combining the two variables and performing pairwise PERMANOVA tests showed that all treatment combinations significantly differed in the autumn sampling (Table 3).Similarly, most groups significantly differed in the spring sampling, except for the Organic-BMc and Organic-Control combination (Table 3).While both farming practices and BMc inoculation influenced bacterial community composition, they did not consistently affect a-diversity.Specifically, farming practice significantly affected ASV richness, Shannon index, and Evenness in the spring but not in the autumn sampling (Supplementary Figure 2, Supplementary Table 5).Irrespective of the sampling time, BMc inoculation did not influence bacterial a-diversity and no interaction effect between the two experimental factors was detected (Supplementary Figure 2).
We aimed to evaluate the ASVs that were affected by either the choice of the farming practice and/or the BMc inoculation.In total, 371 ASVs in the autumn and 273 ASVs in the spring sampling significantly differed between the two farming practices (Supplementary Table 6, ANCOM-BC2).Out of these differential abundant ASVs, 38 ASVs in the autumn and 16 ASVs in the spring sampling were representative of the bacterial communities investigated as they showed a relative abundance above 0.5% (Figure 4).Dominant discriminant ASVs with the highest differential abundance due to choices of farming practice Abundance of beneficial microorganisms in winter rye (cv.Conduct) cultivated under different long-time farming practices (conventional vs. organic) in autumn and spring of the same growing season.Bars represent means ± standard deviation of four replicates.A two-way ANOVA was conducted to test the effect of the farming practice and the sampling time on the rhizosphere competence of each member of the consortium (Supplementary Table 1).Asterisk indicates a significant difference of the CFU between the sampling time points for each consortium member separately ( * p ≤ 0.05).

B A FIGURE 2
Effect of beneficial microorganisms consortium (BMc) on autumn and spring shoot dry mass (SDM) (A) and yield (B) of winter rye (cv.Conduct) cultivated under conventional and organic long-time farming practices.Bars represent means ± standard deviation of four replicates (Supplementary Table 2).Different letters indicate significant differences calculated by the Sidak-test (p ≤ 0.05).The type of letter indicates a separate statistical analysis (lowercase: autumn, uppercase: spring).6).Specifically, the relative abundance of ASV44 (Rhizobium spp.) increased, while the relative abundance of ASV40 (Pedobacter spp.) decreased due to organic farming practice (Figure 4).In contrast, the ASV14 (Pedobacter spp.) showed opposite trends between the two samplings under organic farming, with an increase in abundance in autumn and a decline in spring.The BMc inoculation altered the relative abundance of 21 ASVs in the autumn and 14 ASVs in the spring sampling (Supplementary Table 6), with five ASVs in the autumn and three ASVs in the spring accounting for more than 0.5% of relative abundance (Figure 4).Two ASVs classified as Luteibacter spp.(ASV48 and ASV61) consistently occurred in higher relative abundance in the BMc inoculated samples.Other ASVs that increased due to BMc inoculation belonged to Mucilaginibacter spp.(ASV37) and Pedobacter spp.(ASV78, ASV73).Moreover, one dominant ASVs, which belonged to Pedobacter spp.(ASV31) was negatively affected by BMc inoculation in the spring sampling.
We were then interested to elucidate ASVs that significantly differed between each treatment combination.In total, 304 and 236 ASVs were found significantly different between the treatment combinations in the autumn and spring samplings, respectively (Supplementary Table 7), with 43 ASVs significantly differing independent from the sampling time (Figure 4, Supplementary Table 8).The logistic regression models revealed 277 (autumn sampling) and 210 (spring sampling) ASVs as representative predictors for each specific treatment combination (Supplementary Table 9).Of these ASVs, the relative abundance of 30 and 46 ASVs increased in the BMc-inoculated samples under organic fertilization in autumn and spring, respectively (Supplementary Table 9), where significantly higher SDM was detected (Figure 2A).In contrast, in the organic farming, BMc inoculation resulted in a decreased relative abundance of 57 ASVs at the autumn sampling and 19 ASVs at the spring sampling (Supplementary Table 9).Furthermore, three ASVs consistently responded to the BMc-inoculation in organic treatment in each sampling campaign (Supplementary Table 10).Out of these differentially abundant ASV between the four treatment combinations, 28 and 12 ASVs exceeded the relative abundance of 0.5% in at least one sample in the autumn (Figure 5A) and the spring samplings, respectively (Figure 5B). ) was the single top abundant, significantly different ASV that negatively responded to BMc inoculation in the organic treatment at both sampling time points (Figure 5).Some of these ASVs were also associated with other treatments, especially in the control of conventional farming.For example, the ASV88 (Chitinophaga spp.) was enriched in the control samples under conventional farming and depleted with BMc inoculation in the organic treatment at both sampling time points.However, in the spring sampling, ASV61 (Luteibacter spp.) responded positively to BMc inoculation in the organic treatment and negatively to control in conventional farming (Figure 5).
Finally, we performed a phylogenetic clustering analysis to estimate whether ASV responders to the BMc inoculation in the organic treatment were associated phylogenetically with the inoculated BM strains (Figure 6).ASV responders classified as Bacillus spp.clustered more closely to other Bacillus spp.than the 16S rRNA genes of B. atrophaeus ABi03 (Figure 6A).In contrast, two out of the three ASVs classified as Pseudomonas spp.(ASV1123 and ASV1576), were clustered together with the six 16S rRNA gene copies of Pseudomonas sp.RU47 (Figure 6B).These two ASVs were positively associated with the BMc inoculation in the organic treatment in the autumn sampling but not in the spring sampling (Supplementary Figure 3).

Discussion
Successful establishment of the microbial consortium throughout winter dormancy In this study, we hypothesized that each consortium member would colonize the rhizosphere of winter rye plants and maintain a high density over the winter dormancy independent of farming practices (organic vs. conventional).The three consortium members were selected based on their previously observed positive effects on plant performance.Pseudomonas sp.RU47 strain enhanced growth of tomato and maize plants in P-deficient soils (Eltlbany et al., 2019) and suppressed bottom rot (Rhizoctonia solani) in lettuce very efficiently under greenhouse and field conditions, although no strong antagonistic properties were found in in vitro assays (Adesina et al., 2009;Schreiter et al., 2018).We assumed that the disease suppression effect was mainly based on the production of HCN (Kuzmanovićet al., 2018).In contrast, the Bacillus atrophaeus ABi03 strain showed various antagonistic functions in vitro.Plant growth promotion was also observed in consortia containing the Trichoderma harzianum OMG16 strain in tomato (Mpanga et al., 2018) and maize (Mpanga et al., 2019;Moradtalab et al., 2020).Furthermore, OMG16 reduced root infection by Verticillium longisporum in rapeseed as well (Hafiz et al., 2022).In this study, the strains with different modes of action were applied as a consortium in the field.Several studies have shown that using a consortium of two or more microorganisms can have a stronger and more sustainable effect on plant health and performance than a single beneficial microorganism (Sharma et al., 2020;Hafiz et al., 2022).Additive or synergistic effects can, in part, be attributed to separate habitats such as the rhizosphere, the rootassociated soil, and the root cortex but also to different beneficial functions such as improved nutrient availability, phytohormone production, antagonistic functions against plant pathogens and (a)biotic stress mitigation (Jha and Saraf, 2012;De Vrieze et al., 2018;Gu et al., 2020;Santoyo et al., 2021).
A long-term coexistence between the members of a BMc in the rhizosphere or root-associated soil is crucial to ensure beneficial effects on plants in short and long periods.Seven weeks after inoculation (autumn sampling), all members of the consortium colonized the winter rye roots in a similar density compared to previous studies (Schreiter et al., 2018;Jamil et al., 2021), indicating that the combined application of the beneficial microorganisms as a consortium did not impair their rhizosphere competence.The bacterial strains ABi03 and RU47 were mainly detected in the rhizosphere and did not inhibit each other as observed in dual culture (data not shown), whereas OMG16 inhabits especially the root-associated soil but was also previously reported as an endophyte in rapeseed (Hafiz et al., 2022).Because of this spatial Besides direct competition, the establishment of applied microorganisms is highly dependent on the composition and structure of the indigenous microbial community, as demonstrated by previous studies that linked a low rhizosphere competence of applied microorganisms to a high microbial diversity, which acts as a barrier against external invaders (Schierstaedt et al., 2020;Mawarda et al., 2022).The successful establishment of inoculated microorganisms within the rhizosphere Log 2 fold-change (FC) of ASVs that significantly differed and occurred in relative abundance higher than 0.5%.The comparisons were performed between farming practices (Organic vs. Conventional) or BMc inoculation (BMc vs. Ctrl) in the autumn or the spring sampling campaign during the growth of winter rye plants.Differential abundance testing was performed on rarefied datasets via ANCOM-BC2 with Benjamini-Hochberg correction (adjusted p-value<0.05).
has been linked to rapid nutrient utilization, biofilm formation, but also antagonistic capabilities (Adam et al., 2016;Kumawat et al., 2019).In vitro demonstrated antagonistic functions of all consortium members may facilitate their selective occupation of niches within the rhizosphere, thus promoting their integration into the microbial community (Berg et al., 2021).
Although our results (Figure 3) confirm that farming practice drives the composition of microbial communities in the soil and Heatmaps with the log 10 transformed mean relative abundance of differential abundant ASVs in the autumn (A) and spring (B) sampling.The presented ASVs occurred in relative abundance higher than 0.5% in at least one sample.Significantly different ASVs were identified via ANCOM-BC2 (adjusted p-value via Benjamini-Hochberg <0.05, n=4).The colored legend (left side) shows whether these ASVs responded positively or negatively to a particular treatment combination (Organic-Ctrl, Organic-BMc, Conventional-Ctrl, Conventional-BMc).The responder ASVs were identified via logistic regression (adjusted p-value via Benjamini-Hochberg < 0.05).RA, Relative abundance.
rhizosphere (Schrama et al., 2018;Babin et al., 2019;Bziuk et al., 2021;Fernandez-Gnecco et al., 2022), the BMc rhizosphere competence was not affected, indicating a robust competence of the consortium members at different environmental conditions.Since only a few responder ASVs consistently differed between the farming practices at both sampling time points, we presume that the assembly of rhizosphere microbiota was a stochastic process influenced by environmental conditions and plant developmental stages.The continuous time-dependent stochastic assembly of bacterial communities in the rhizosphere could have contributed to the high rhizosphere competence of the consortium members independent of farming practices.In spring, the colonization density of RU47 in the rhizosphere decreased overall by a factor of 10, whereas the densities of ABi03 and OMG16 (root-associated soil) remained almost the same compared to the autumn sampling.The higher resilience of ABi03 and OMG16 can be attributed to their potential for sporulation, which allows higher survival under adverse winter conditions compared to RU47 (Fernańdez-Sandoval et al., 2012;Sella et al., 2014).Nevertheless, all three microorganisms colonized the rhizosphere of winter rye sufficiently, even 22 weeks after inoculation.This indicates a strong positive feedback between the plant and the members of the consortium since the persistence of a high colonization density depends on the one side on the root exudation profile of the plant and on the other side on the chemotaxis of each consortium member towards these metabolites and their use as a substrate (Malgioglio et al., 2022).

Consortium improved plant performance and nutrient acquisition under organic farming
Shoot growth of winter rye was promoted by BMc inoculation in autumn (+20%), but especially during regrowth after winter (+45%) exclusively under organic farming.These findings are in line with the robust rhizosphere competence shown by the investigated inoculants, which remained unaffected even during unfavorable winter conditions.Moreover, under long-time organic farming practice, BMc inoculation was associated with the enrichment of bacterial taxa with documented plant growth promoting properties.As a consequence, organic farming practice may support soil biological processes that allow more complex taxonomic and functional communities (Bender et al., 2016).Numerous studies have consistently reported positive effects on beneficial plant-microbe interactions when using manure-based fertilizers alongside inoculants such as Bacillus, Pseudomonas, and Trichoderma (Thonar et al., 2017;Mpanga et al., 2018;Bradaćǒvaé t al., 2019).A beneficial plant-microbe interaction under organic farming might improve the carbon supply of fast-growing copiotrophic inoculants as well as indigenous plant growthpromoting microorganisms, thereby supporting the establishment of a beneficial microbial community.This effect may be of particular importance especially on light sandy soils characterized by low organic carbon content, which were observed in this study even after long-time organic farming practice (Supplementary Table 4).Accordingly, a meta-analysis reported that responsiveness to microbial inoculants decreases with increasing soil organic carbon content (Schütz et al., 2018).Furthermore, the high availability of N and P, which is characteristic of manure-based fertilizers, can serve as a starter fertilization for the host plant, facilitating root growth and the establishment of microbial inoculants in the rhizosphere (Bittman et al., 2006;Chekanai et al., 2018).The beneficial effects of inoculants and their impact on the rhizosphere microbiome can in part be attributed to improved plant nutrient supply.During autumn, no significant differences in shoot nutrient concentrations were initially observed among the treatments reflecting a comparable nutritional status (Table 2).However, ANOVA revealed a general impact of the BMc inoculation on the N concentration of the shoots (Supplementary Table 3).Nutrient deficiencies were only apparent for Ca (Bergmann, 1988) and may be attributed to low Ca levels of the light sandy soil at the field site (Supplementary Table 4).In contrast, during regrowth after winter, critical nutrient concentrations were recorded, particularly under organic fertilization, including N, P, K, Mg, Ca, Cu and Zn.BMc inoculation significantly improved the P nutritional status, which is a critical nutrient in many organic farming systems (Cooper et al., 2018).However, similar trends were also recorded for all remaining nutrients except Zn and Mn (Table 2).An improved plant nutritional status due to BMc inoculation is also supported by ANOVA, which shows a general effect of BMc inoculation on the concentration of these nutrients independent of farming practice (Supplementary Table 3).As a consequence, nutrient concentrations of BMc inoculated plants reached or even exceeded the sufficiency thresholds particularly for N and P.These findings suggest a general positive effect of BMc on the nutrient acquisition, possibly mediated through the stimulation of root growth, which has been well-documented for the inoculant strains (Mpanga et al., 2018;Eltlbany et al., 2019;Moradtalab et al., 2020).However, improved plant nutrient acquisition through direct nutrient mobilizing by the BMc is a less likely scenario, although in vitro tests did detect solubilization of P and K. Analysis of the root-associated nutrient pools in the soil did not show any changes following the BMc inoculation, suggesting no specific effects on solubilization of available soil nutrients in spring (Supplementary Table 4).According to current literature, P-solubilizing microorganisms contribute to the host plant nutrition primarily through their long-term impact on nutrient cycling, rather than direct nutrient solubilization (Raymond et al., 2021).Interestingly, similar trends of an improved nutrient status after BMc inoculation were also shown for plants grown under conventional farming.However, compared to organic farming, conventional farming practice generally resulted in a higher nutritional status of the plants, reaching sufficiency levels for N, P, and Mg, even in non-inoculated plants.This may explain the absence of additional growth responses by BMc inoculation under conventional farming.
Although there were no significant differences among the treatments, the final yield data (Figure 2B) corresponded with the effects of the shoot biomass during spring vegetative growth (Figure 2A).The grain yield under long-time conventional farming practices was slightly higher than the 2015-2020 average rye yields (4.1 t ha -1 ) in the district of Brandenburg (Germany).Organic farming practice tended to decrease yield by approx.22%, similar to the reported average yield losses determined by metastudies (de Ponti et al., 2012;Seufert et al., 2012).Although not significant, our results indicate that this yield decline by organic farming may partially be compensated by BMc application.

Inoculation-dependent modulation of rhizosphere bacterial community persists over winter dormancy
We hypothesized that the application of the BMc would affect the bacterial community composition in the rhizosphere depending on farming practice.Consequently, we anticipated that these variations would exert differential influences on the performance of winter rye.An alteration of the bacterial community composition in the rhizosphere of winter rye was found in the BMc inoculated treatments, regardless of farming practices.This effect persisted until spring, indicating a prolonged impact of the BMc (Figure 3).Similar results were observed in a study by Deng et al. (2019), where BMc inoculation affected the bacterial community composition of strawberry roots over various time points, irrespective of farming practices.This indicates that repeated application of the inoculants might not be necessary, once their effects have been established.However, further time-series experiments are needed to clarify whether the effects of BMc on the rhizosphere microbiota and plant performance persist throughout the growing season at the field scale.
In addition to BMc inoculation, we observed an effect of the farming practices on the soil microbiota, similar to our previous studies (Bziuk et al., 2021;Windisch et al., 2021;Fernandez-Gnecco et al., 2022;Sommermann et al., 2022).However, this is the first time we report that this effect persists across different plant developmental stages.While both farming practices and BMc inoculation significantly affected bacterial community composition, only a few ASVs exhibited consistent differential abundance patterns over time (Figure 4).This indicates a) that the recruitment of microorganisms in the rhizosphere differed over the developmental stages of the plant and b) that the BMc inoculation and the farming practices influenced this recruitment.A similar influence of farming practices on the recruitment of microorganisms in the rhizosphere has been demonstrated previously (Windisch et al., 2021;Sommermann et al., 2022).In summary, the BMc inoculation influence on rhizosphere microbiota persisted over the winter dormancy, indicating specific strategies such as spore formation by the consortium members to survive during winter time.In addition, changes in the root exudation during spring could have further activated dormant BMs (Windisch et al., 2017).
According to ecological theories on plant-microbe interactions, the mutualistic selection of plant beneficial microorganisms often emerges under nutrient limitations (Sańchez-Cañizares et al., 2017).Since chemical fertilizers provide a high availability of nutrients to plant and rhizosphere microbiota, they could interfere with the plant's selection processes for beneficial microorganisms (Sańchez-Cañizares et al., 2017).Interestingly, we observed increased growth performance of plants inoculated with the BMc only under organic farming, where nutrient availability was low (Supplementary Table 4).We believe that this additive effect could be associated with a higher selection pressure on both the plants and soil microorganisms due to nutrient limitations, in contrast to conventional farming with high nutrient availability.However, CFU data indicated that the farming practices negligibly altered the density of the applied BMc in the soil (Figure 1).Herein, the differential abundance testing revealed an enrichment of ASVs closely related to Pseudomonas spp. or Bacillus spp., which may correspond to our BMc inoculants.By applying phylogenetic clustering, we identified two Pseudomonas ASVs closely related to Pseudomonas sp.RU47 (Figure 6).The higher proportions (relative abundance) of Pseudomonas sp.RU47 under organic farming and BMc inoculation, might indicate less competition from other bacterial taxa inhabiting the rhizosphere in organic farming.The interaction effect between farming practices and BMc inoculation on the rhizosphere bacterial community profiles was weak.Thus, the two farming practices only slightly influenced the ability of BMc to modify the rhizosphere bacterial community.Nevertheless, we attempted to identify ASVs associated with BMc inoculation under organic farming to determine which fraction of bacterial taxa are associated with the inoculated BMc.Several ASV responders were detected, but most had a relative abundance below 0.5%.
Among the dominant responders (>0.5% relative abundance), several ASVs closely classified as known potential plant-beneficial bacteria increased in relative abundance in the rhizosphere of winter rye (Figure 5).For instance, members of the genus Luteibacter (which belong to g-Proteobacteria) are known for their plant-growth-promoting traits (Guglielmetti et al., 2013;Hoffman et al., 2013) and were enriched in the rhizosphere of the organic treatment.Interestingly, Luteibacter spp.can act as symbionts to fungal plant endophytes and trigger the production of IAA, which promotes plant growth.Similarly, one ASV classified as Pedobacter spp.(ASV95), a genus that includes several plantbeneficial species, increased in relative abundance due to BMc inoculation under organic farming (Morais et al., 2019).While the combination of BMc inoculation and organic farming might not substantially affect the overall profile of the bacterial community, it increased the relative abundance of bacterial taxa that potentially act as plant-beneficial microorganisms promoting plant growth.However, further field studies are needed to confirm whether the microbial responders due to BMc inoculation also contributed to increased plant growth.
Among the dominant responders with decreasing relative abundance due to BMc inoculation, primarily in organic farming, we identified an ASV classified as Chitinophaga sp. or Duganella sp.These ASVs have been previously associated with the onset of plant diseases (Carrioń et al., 2019;Li et al., 2021).Organic farming practice and the BMc combination probably reduced bacterial taxa associated with plant diseases and stressful conditions.In contrast, some of these ASVs had higher relative abundance in conventional farming practice without BMc inoculation, indicating that the high input of nutrients through chemical fertilization might have led to the accumulation of bacterial taxa associated with plant diseases.Alternatively, plants often recruit several closely taxonomicallyrelated bacteria in response to a pathogen (e.g., Duganella spp.) (Haack et al., 2016), so direct association of these taxa with plant pathogenicity is not possible, especially since no plant disease phenotype was detected.

Conclusion
This study illustrated the effects of a beneficial microorganism consortium (BMc) on the performance of a winter rye depending on different farming strategies, i.e. conventional and organic farming.The main goal of our study was to assess the ability of the consortium members to survive in a sufficient density in the rhizosphere/root-associated soil of winter rye as prerequisite for successful plantmicroorganism interaction.Our findings demonstrated that the consortium members persisted in the rye rhizosphere over the vegetation period, maintaining a high density even after winter dormancy independent of the farming practice.As expected, the BMc had a positive effect on the rye performance, increasing shoot dry biomass and plant nutritional status especially under organic fertilization.These findings further confirmed the in vitro tested plant-promoting features that characterized the microbial members inoculated in the rhizosphere of winter rye under field conditions.Moreover, the inoculated BMc had a significant effect on the bacterial community dynamics as large shifts in rhizosphere bacterial assembly were observed after inoculation in both organic and conventional treatments, and such community shifts were also detected across different plant developmental stages.Interestingly, under organic farming we observed that in the rye rhizosphere, treated with BMc, several bacterial taxa previously associated with plant diseases were depleted, while several putative plant beneficial bacterial taxa were enriched in their relative abundance, thus further highlighting the positive effects of the BMc on the plant health by modulating the rye rhizosphere microbiome.It is noteworthy that although our field experiment indicated a positive effect of BMc on rye performance, especially under organic fertilization, the beneficial impact on rye yield was only marginal.Nevertheless, our results indicate that BMc inoculation might have the potential to compensate yield losses caused by nutrient limitation in organic farming practice.Furthermore, BMc treatment may contribute to a better yield stability by improving plant nutrient status and promoting beneficial microbiota, thus representing a sustainable approach especially in combination with organic farming.Our research emphasizes the importance of conducting additional field experiments to understand the efficacy of microbial inoculants in various cropping systems, promoting the development of efficient microbe-based solutions for sustainable agriculture.

FIGURE 4
FIGURE 4 FIGURE 6(A) Ordination plots based on phylogenetic distance following pairwise local alignment of 16S rRNA gene sequence from the ASV responders to Organic-BMc that were classified as Bacillus spp., 16S rRNA genes from Bacillus atrophaeus ABi03 and Bacillus spp.from the type strain collection of Genbank (Red: ASV responders, Blue: Bacillus atrophaeus ABi03, Green: Bacillus spp.type strains).The plot shows no indication that these ASV responders classified as Bacillus spp.could belong to Bacillus ABi03.(B) Ordination plots based on phylogenetic distance following pairwise local alignment of 16S rRNA gene sequence from the ASV responders to Organic-BMc that were classified as Pseudomonas spp., the strain Pseudomonas sp.RU47 (six 16S rRNA genes) and Pseudomonas spp.from the type strain collection of Genbank (Red: ASV responders, Blue: Pseudomonas sp.RU47, Green: Pseudomonas spp.type strains).The plot indicates that two ASV responders (ASV1123 and ASV1576) probably belong to the inoculum strain Pseudomonas sp.RU47, since they cluster with the 16S rRNA genes of the inoculum strain Pseudomonas sp.RU47.

TABLE 1
In vitro characterization of beneficial functions of the inoculated consortium members.

TABLE 2
Nutrient status of winter rye (cv.Conduct) grown under different long-time farming practices and use of beneficial microorganisms (control vs. BMc) in autumn and spring of the same growing season.

TABLE 3
Pairwise PERMANOVA tests of the combined treatments (farming practices: conventional or organic, use of beneficial microorganisms: BMc or Control) with Benjamini-Hochberg correction (n=4).
separation, potential antagonistic functions demonstrated in vitro for RU47 and ABi03 (Table1) can likely not act against OMG16.