Catch the Best: Novel Screening Strategy to Select Stress Protecting Agents for Crop Plants

Climate change increases stress levels for crops and affects the economic and environmental aspects of agricultural management systems. The application of stress tolerance-mediating microorganisms is an auspicious strategy for improving crop protection, and as such, we developed a direct selection strategy to obtain cultivable microorganisms from promising bioresources using the bait plants, maize, oilseed rape, sorghum and sugar beet. Alpine mosses, lichens and primrose were selected as bioresources, as each is adapted to adverse environmental conditions. A 10% crop-specific selection was found for bait plant rhizosphere communities using cultivation-independent fingerprints, and their potential role as stress protecting agents (SPA) was evaluated following the cultivation of captured bacteria. In addition to assays identifying phytopathogen antagonism and plant growth promotion capacities, our evaluation included those that test the ability to allocate nutrients. Moreover, we developed new assays to measure tolerance in diverse stress conditions. A score scheme was applied to select SPAs with desired properties, and three Pseudomonas species with pronounced antagonistic activity that showed elevated tolerance to desiccation and an improved seed germination rate were subsequently chosen. Screening for environmentally-conditioned and host-adapted microorganisms provides a novel tool for target-oriented exploitation of microbial bioresources for the management of ecofriendly crops facing biotic and abiotic stresses.


Introduction
The impact of climate change has had many dramatic and multifaceted ecological consequences on natural and human-managed systems, apart from the prevalent threat of phytopathogens and nutrient deficiency on monocultured crops.Plants are affected by increasing temperatures, salination, drought and solar radiation, and with these direct biological effects come rapidly evolving pathogens with shifting host and geographic distribution, resulting in reduced crop yield and quality [1][2][3].The most effective targets for improved sustainability are those with tolerance traits for abiotic and biotic stresses that match with the associated microbiome, the so-called second plant genome [4].For biological control of phytopathogens, beneficial microorganisms were developed as biological control agents (BCAs) to act as an alternative to conventional pesticides.However, to overcome stress conditions, stress protecting agents (SPAs) [5], which are microorganisms with the ability to mediate stress protection to the host plant and to survive under adverse abiotic conditions, are required.
The first step in our study was to develop an ecologically-based screening strategy for where to search for SPAs.Plants are associated with highly diverse functional microorganisms, which promote plant growth, health and stress tolerance [6].In natural environments, microorganisms' versatility and ability to cope with harsh conditions is facilitated by various elaborate enzymes and sophisticated mechanisms.Thus, plants from natural environments can be a promising source to screen for bacteria and fungi with properties relevant to almost all fields of biotechnology.Today, numerous microbes initially isolated from plant habitats are employed in biotechnological processes, especially in agriculture [7,8].Microorganisms associated with various plant microenvironments constitute a rich source for applications in biotechnological processes and sustainable agriculture.Thus, the rhizosphere and endosphere have been found as particularly interesting reservoirs for potential microorganisms [6] capable of promoting plant growth and inducing systemic tolerance towards salt and drought [9].Recently, various species of lichens and mosses have emerged as potential bioresources for antagonistic and plant growth promoting (PGP) microorganisms [10][11][12].Additionally, antagonistic microorganisms have been found to be associated with mistletoe species and sclerotia from Rhizoctonia solani [13,14].Although the mentioned bioresources are well-known and interesting, it is still difficult to isolate microorganisms from these sources, because they need specific or currently unknown conditions for their growth in vitro.The development of biological plant protection products usually begins with the enrichment of microorganisms on artificial growth media.Regardless of the medium used for cultivation, two major problems arise.Firstly, only 1% to 10% of microbial species are able to grow on conventional media [15].Secondly, faster growing microbes typically outgrow the slower ones.As a result, current cultivation techniques yield a low microbial diversity.Furthermore, most of the applied enrichment strategies neglect the conditions to which the selected microorganisms will be later applied.Superior growth under laboratory conditions does not necessarily guarantee a good performance ad planta.Therefore, a second question arose as to how we could gather suitable SPAs.
The objective of this study was to use bait plants for a more directed selection of suitable microorganisms from stress-exposed bioresources.We used pathogen-and low pH-tolerant mosses [11], cold-adapted primrose [7,8] and drought-adapted lichens [12] as natural sources for microorganisms.The microorganisms from the original host were captured by priming the seeds of the bait plants, maize, oilseed rape, sorghum and sugar beet, and, subsequently, culturing the microbial community from the rhizosphere of each bait plant.The screening procedure was accompanied by microbial community fingerprinting using the single strand conformation polymorphism (SSCP) of 16S rRNA genes [16], which is useful for a fast and safe evaluation of microbial shifts [17].Re-cultivated bacteria were intensively screened for their potential as biocontrol agents and stress protection mediators.The screening strategy included assays to identify: (1) antagonistic activity against diverse Ascomycota (Alternaria alternata, Botrytis cinerea, Sclerotium sclerotiorum, Verticillium dahliae), Basidiomycota (Rhizoctonia solani AG2-2IIIB, R. solani AG4, Sclerotium rolfsii) and the oomycete, Phytophthora infestans; (2) nutrient allocation (iron, nitrogen, phosphate); (3) exoenzyme activity (chitinase, glucanase, protease); (4) antibiotic genes (pyrrolnitrin prnD, diacetylphloroglucinol phlD); (5) growth in the presence of diverse stress conditions (drought, reactive oxygen species, salt, temperature); (6) seed germination promotion; as well as (7) biosurfactant production.There is a general trend towards greater suppression or enhanced consistency of the antagonistic and stress protecting effect of microorganisms.Therefore, a scoring scheme was applied to specifically select desired properties as an example for the most suitable BCAs and SPAs in potential applications, and the resulting strategy provides a novel tool for target-oriented exploitation of bioresources.

Bait Plant Strategy
All microbial communities extracted from the field-collected moss, lichen and primrose from their natural environments were used to prime seeds of maize, oilseed rape, sorghum and sugar beet.After two weeks of cultivation in pots with standardized potting soil, the structure and function of the microbial community was analyzed with cultivation-dependent and -independent methods.The resulting bacterial isolates were assumed to be environmentally-conditioned and host-adapted (Figure 1).

Insight into Rhizosphere Communities Using Molecular Fingerprints
The bioresources were selected due to their potential to survive under particular adverse abiotic conditions.For example, Sphagnum mosses and their diverse microbiome undergo repetitive desiccation and oxidative stress.The bacterial capacity to tolerate oxidative stress, due to the high abundance of genes responsible for the oxidative stress response may determine the effective and stable colonization of the Sphagnum mosses [18].Figure 2 shows microbial fingerprints of the moss-associated community captured by the bait plants.Interestingly, in each capture experiment, a portion of the moss-associated community was able to colonize the crops, and the bacterial community patterns from moss-treated crops showed host-specific communities.Bacteria derived from the moss community were specifically captured by the different crops; the bacterial communities captured through oilseed rape and sugar beet differed significantly (p ≤ 0.0360, p ≤ 0.0144) from the other crops, and the captured communities of maize and sorghum differed at p ≤ 0.2980 and p ≤ 0.1827, respectively.Comparing bacterial community patterns from moss-microbiome-treated and non-treated control plants (sterile sodium chloride treatment), no significant differences were found.The microbial patterns resembled each other at a similarity level of 90% (Figure 3).The highly similar microbial patterns observed for moss-treated crops compared to the controls showed that the bacteria added to the bait plants in target host treatments need to be both compatible and competitive within the already existing and established microbiome of the unsterile soil and seed community.In this way, microorganisms with antagonistic activity are typically known for their competitive nature and other beneficial properties, and it is thus more likely that they successfully establish in new environments [19].As soil is the main reservoir for rhizosphere microorganisms [20], the plant host habitat and root exudates shape the soil bacterial community structure.A subset of the soil bacterial community is present around the plant roots, and active community members benefit from its basic functions and, in turn, help the plant thrive in specific environments [21,22].Specific enrichment of microorganisms is already a standard in laboratory practices through the application of defined media; however, the application of natural bait material is a rarely applied technique in capturing potential target-adapted antagonists.Similar experiments were previously successfully applied, such as using Rhizoctonia solani sclerotia as bait material deposited in soil to enrich antagonistic counterparts [14].The SSCP microbial fingerprint analysis was not suitable to detect bacterial populations derived from the treatments, as only the first one hundred dominant species are detectable.However, the most matching and stress protection mediating bacteria could be a part of the minor fraction of the microbial community.The importance of minor abundant microorganisms' function is gaining recognition in scientific research and discussion, and a better understanding is becoming possible, due to next generation sequencing methods that allow deeper insight into the microbial communities.Unweighted pair group method with arithmetic mean (UPGMA) dendrogram of bacterial rhizosphere communities of moss-primed and untreated crop plants in two replicates (first number) with three repetitions per plant species (second number).The dendrograms were generated from SSCP community profiles using GelCompar II ® .The following settings were used: dendrogram type: UPGMA; similarity coefficient: curve-based Pearson correlation; optimization: 4%; position tolerance: 1%.The dashed line indicates a similarity level at 90%.Letter code: first letter of the bait plant species; M, maize; R, oilseed rape; S, sorghum; Z, sugar beet; the second letter is the treatment; M, moss; C, control.

Characterization of Selected Bacteria and Score Evaluation for SPAs
In the first step, 283 microorganisms were isolated from the rhizospheres of bait plants.All bacteria antagonistic (105) against at least one of the tested pathogens, Alternaria alternata, Botrytis cinerea, Sclerotium sclerotiorum, Verticillium dahliae, Rhizoctonia solani AG2-2IIIB and AG4, Sclerotium rolfsii and Phytophthora infestans, were further characterized using various assays (Table 1).The screenings was comprised of parameters of: (1) nutrient allocation (iron, nitrogen, phosphate); (2) exoenzyme activity (chitinase, glucanase, protease); (3) antibiotic genes (pyrrolnitrin prnD, diacetylphloroglucinol phlD); (4) growth in the presence of diverse stress conditions (reactive oxygen species, drought, salt, temperature); (5) seed germination promotion; as well as (6) biosurfactant production; which were evaluated and weighted (one to five) according to their intended field application.Each parameter was defined with "one" for a positive reaction or maximum data range comparing values of all tested bacteria and "zero" as a negative or minimum data range.Positive values were multiplied with the defined evaluation weight, and all values were summarized to the final score.The higher the final score, the more the bacteria were assumed to be suitable for particular applications as SPAs.The weighted rating is adjustable in accordance to the target application.We defined a rating in order to yield the most potential antagonists with lytic and phosphate-solubilizing activity, which could efficiently enhance the mobilization of nutrients for themselves and for the associated host plant [12].Additionally, detected genes that decode for antibiotic synthesis were excluded to advance other mechanisms involved in biocontrol, e.g., the production of volatiles.Due to elevated stress levels in crops as a result of climate change, the application of stress tolerance-mediating microorganisms is a new possibility to improve crop protection.Therefore, we highlighted bacteria that were able to grow in the presence of stresses, i.e., desiccation, low temperatures and oxidative agents.We valued protease production less, because it is a general mechanism and does not help differentiate between antagonists.Moreover, the score assists in the selection of bacteria for certain target applications, yet remains adjustable for different focuses.To address certain host plants, we set various values in plant growth promotion assays.Additionally, we differentiated on the basis of seed germination after three days and plant emergence after six days to evaluate the possible positive or negative effects of the applied bacteria.Only treatments comparable or better than the controls were highly weighted.The score of the isolates ranged from 10 (one isolate) to 47 (two isolates), whereas the majority shared the score of 32 (10) (Table 2).The three highest assessed bacteria, RM1-3-1 (KA38), SP2-1-2 (KA77) and MP1-3-2 (KA23), shared the score 47, 47 and 46, respectively.They are each a member of the genus, Pseudomonas, have pronounced antagonistic activity, can be re-cultivated after desiccation and improved the germination rate of seeds.Several of the analyzed parameters, e.g., antagonism, exoenzyme activity and PGP, are well known to characterize BCAs [23].Plants and their microbiome are similarly influenced by environmental conditions with specific biotic and abiotic factors [6].The novelties of this study are the directed selection by bait plants and, particularly, the application of diverse stress assays, which support the potential of BCAs to act as SPAs (parameter 5).Interestingly, the majority of the screened isolates could grow at low temperatures (91.4% to 95.2%) and tolerated desiccation for 24 h (81.9%), yet only a few specialists were able to withstand oxidative stresses (6.7% to 15.2%).To select a diverse group of bacteria and also to maintain quality control for further application processes, bacteria were genetically characterized using fingerprints from ARDRA and BOX-PCR.The antagonists were grouped according to the similarity index at 70% into 44 various BOX groups (Figure S1, Table 2).One isolate from the most dominant group was sequenced, and species identification was determined after BLASTn analysis.Group 13 included the highest number of isolates (29), which were assigned to Pseudomonas kilonensis/P.corrugata.
Registration procedures inhibit the progress of the commercial application of bioproducts, due to their intensive costs and time requirements [24].To improve positive product evaluation, potential human pathogenicity has to be identified and excluded from further work in advance.In an extensive study, a broad range of BCAs, pathogens and plant-associated bacteria were used in a rapid and inexpensive slow killing assay using Caenorhabditis elegans to estimate the risk of bacteria to human health [25].In the current study, we analyzed the potential of the SPAs to grow at 37 °C, the human body temperature, and sequenced the partial 16S rRNA genes of antagonists.A high percentage (81.0 to 81.9%) of the bacteria, including the highest scored bacteria, RM1-3-1 (KA38), SP2-1-2 (KA77) and MP1-3-2 (KA23), were able to grow at 37 °C; however, they were identified as Pseudomonas species in risk group 1. Identified species with potential human pathogenicity (risk group two and higher; Table 2) will be excluded from the next step of analysis in both greenhouse and field trials.Sequence similarities based on the 16S rRNA genes were visualized in a phylogenetic tree (Figure 4).The biggest ARDRA groups (10 with 46 isolates, 16 with 21 isolates) with a similarity index of 80% and BOX groups (Table 2) from the diverse Pseudomonas clusters were dominated by P. kilonensis and P. corrugata, and minorly represented by P. brassicacearum subsp.brassicacearum and P. poae species.Pseudomonas species are one of the primary bacteria that succeed in applications against phytopathogens [26].
To link the dominant bacteria from cultivation to fingerprints, the randomly picked isolates were included within the microbial fingerprint analysis ladder (Figure 1).Particular bands of the ladders were also found within the community patterns, e.g., Rhodanobacter fulvus strain ZM2-3-2 (KA95, score 26) isolated from the moss treated sugar beet samples was found within the sugar beet fingerprints with the same migration distance.In the moss-treated oilseed rape fingerprints, the bands of the cultivated bacteria, Arthrobacter nicotinovorans MM1-3-4 (KA14, score 27), Pseudomonas brassicacearum subsp.brassicacearum RM1-1-3 (KA36, score 26), Pseudomonas otitidis RM1-3-1 (KA38, score 47) and Pantoea sp.RM2-2-4 (KA42, score 36), were detected at the same migration level.The bacteria with the highest scores, i.e., Pseudomonas otitidis RM1-3-1 (KA38) or Pantoea sp.RM2-2-4 (KA42), were not part of the dominant community (weak SSCP band signal), but were successfully captured by the bait plants and selected through the appropriate screenings.However, in addition to the directed selection, those bacteria already present in the host-associated microbial community could also be isolated, e.g., Pseudomonas veronii strain SM1-1-2 (KA63, score 24).Directed selection is the first step in SPA isolation; however, more importantly, it is the evaluation of isolates in regards to their role as mediators of stress protection.The first results from evaluation under greenhouse conditions are promising.The tree is comprised of 47 partial 16S rRNA gene sequences.Neighbor-joining phylogenetic tree (maximum composite likelihood method) was constructed with both one representative sequence per major species and the nearest reference sequences (accession numbers in brackets).The reliability of the tree topology was evaluated using bootstrap analysis with 10,000 re-samplings.Color code: rhizosphere samples from maize, yellow/orange/red; oilseed rape, green shades; sorghum, blue shades; sugar beet, purple shades.Sequences of MF2-1-4 (KA3) and ZM1-1-3 (KA90) were excluded from MEGA5 analysis due to their short sequence length.R ._f u lv u s _ K A 0 9 5 R ._ f u lv u s _ ( A B 1 0 0 6 0 8 . 1 )

Experimental Design for Isolation and Characterization of Stress Protecting Bacteria
The experimental design is comprised of the moss, Sphagnum magellanicum (M), the lichen, Solorina crocea (F), and the Alpine Primula vulgaris (P) for the screening of stress tolerance mediating bacteria (Tables S1 and S2).To ensure that we captured the most diverse microbial community possible, entire individuals of the two collected plant and lichen species were suspended in sterile 0.85% sodium chloride for three min in a stomacher laboratory blender (BagMixer; Interscience, St. Nom, France) in four replicates (P.vulgaris, only 3 replicates).Aliquots of the suspensions were diluted and spread on R2A (the minimal medium to isolate universal bacteria bacteria) for determination of abundances and isolation after incubation for five days at room temperature.For directed selection of microorganisms over the rhizosphere of bait plants (maize, M; oilseed rape, R; sorghum, S; sugar beet, Z), the seeds were primed for four h in two different suspensions per material with three replicates at room temperature under agitation.Seeds were planted in 200-mL pots filled with soil and subsequently cultivated (Gramoflor Profi-Substrat, Topfpikier M + Clay + Fe, pH (CaCl 2 ) 5.8, Garten Bau Centrum, Kalsdorf, Austria, containing white and black peat with moist clay, 90 kg m 3 , PG-Mix 1.0 kg m 3 (Greenworld, Wels, Austria), Radigen ® 50 g m 3 (Terraflor, Iserlohn, Germany)) and incubated under greenhouse conditions (12 h light/dark at 15 °C/10 °C, 60% humidity).Two weeks after germination, when plants developed the first leaves, the rhizospheres were harvested, and microbial communities were suspended as described above.Dilutions were spread on different media as mentioned above.Host-adapted and stress protecting agents (SPAs) were randomly selected (30 per plant and treatment) and characterized.For cultivation-independent analysis, four milliliters of each suspension were centrifuged (16,500× g, 20 min, 4 °C) and stored at −70 °C.Total community DNA was extracted using the FastDNA SPIN Kit for Soil (MP Biomedicals, Solon, USA) according to the manufacturer's protocol and used for fingerprint analysis.Bacterial DNA from single isolates was extracted according to Zachow et al. [27].Bacterial isolates were labeled according to their host, treatment, replicate and continuous numbering; e.g., the third selected bacterial isolate selected from maize treated with the first extract of lichens in the second replicate is MF1-2-3; for practical laboratory work, a consecutive numbering was used, e.g., KA2.

Molecular Analysis
Profiling of microbial communities was carried out with single strand conformational polymorphism analysis (SSCP) on a PhorU2 apparatus (Ingeny, Goes, The Netherlands) and as described by Köberl et al. [7].16S rRNA genes of bacteria were amplified using a universal primer according to the protocol described in Zachow et al. [27], and three independent replicates were included for each sample.Bacterial isolates were screened for genes encoding for nitrogenase (nifH); fragments were amplified using a nested approach with primer pairs nifH3/19F and nifH11/nifH22 [28].Genes for potential antibiotic production (2,4-diacetylphloroglucinol phlD, pyrrolnitrin prnD) were amplified according to Raaijmakers et al. [29] and de Souza and Raaijmakers [30].Amplicons were compared with positive control strains containing the corresponding genes.Amplified ribosomal RNA gene restriction analysis (ARDRA) by applying the restriction endonuclease, HhaI (MP Biomedicals, Eschwege, Germany), to the 1492r and 27f fragment was performed according to Zachow et al. [27].ARDRA patterns were used to group isolates at the genus level, and isolates displaying similar patterns were further analyzed using the PCR-dependent fingerprinting method with repetitive BOX-elements, according to Martin et al. [31].Antagonists with unique BOX patterns and a high evaluation score (see part 2.3) were identified by partial 16S rRNA gene sequence analysis, according to Berg et al. [32].The neighbor-joining phylogenetic tree was constructed using Molecular Evolutionary Genetics Analysis version 5 (MEGA5) [33].

Stress Confrontation Assays
Overnight cultures (ONC) grown in nutrient broth II (NBII, Sifin, Germany) were used as inoculum (10 µL) for all bacterial assays.For osmolarity stress, bacteria were cultivated in 1.5 mL tubes filled with 700 µL NBII with various sodium chloride concentrations (1% steps, 0%-20%) and incubated at 30 °C for 48 h under agitation.For desiccation assays, 20 µL of the ONCs were dried under sterile conditions.After 24 h, bacteria were re-suspended with the same volume of NBII and dropped on nutrient broth II agar (NAII).For growth assays under different temperatures, bacteria were streaked onto NAII and R2A plates and incubated at 13 °C and 37 °C.In reactive oxygen species stress assays, bacteria were cultivated in NBII amended with different tellurite (1, 3, 5, 7, 9 and 10 mg mL −1 ) and H 2 O 2 (100, 300, 500, 700, 900, 1000, 1300 and 1500 mM) concentrations.The assays were performed in 96-well-plates filled with 195 µL medium per well and 5 µL ONC.In all assays, bacterial growth was visually evaluated after 24 h and 48 h.Bacteria were cultivated in 40-mL Luria-Bertani medium (LB, per liter: tryptone 10 g, yeast extract 5 g, sodium chloride 10 g, pH 7.0) at 30 °C to measure antioxidant activity in a radical stress assay with tetramethoxy azobismethylene quinone (TMAMQ).After 24 h, cells were harvested (2500× g, 5 min, 4 °C) and washed three times with the same volume of sterile 0.85% sodium chloride.The antioxidant activity of unstressed (control) and stressed isolates (50 °C for 6 min and 20 min) was measured with TMAMQ developed by Prasetyo et al. [34].

Bioassays
The ability to produce indole-3-acetic acid (IAA) was tested for selected bacterial strains in a plate test according to Sawar and Kremer [38].Nitrogen fixation and production of the ethylene precursor, 1-aminocyclopropane-1-carboxylic acid (ACC) deaminase, was tested according to Fürnkranz et al. [37].Germination of primed seeds (seed number: M, 3; O, 20; S, 15; Z, 20; 1 mL ONC per well of a 24-well plate, 4 h, room temperature, with agitation) was tested on wet filter paper in a 9-cm plate at 13 °C.Germination was evaluated after three and six days in comparison to tap water, distilled water, sodium chloride and unprimed control.The experiment evaluated the first positive or negative impact of the bacteria on seed germination.Only treatments comparable or better than the controls were highly weighted.The best isolates from the pre-screening were tested in a large batch under the same conditions with four repetitions (25 seeds each plant, 10 mL ONC).

Statistics
To compare the SSCP fingerprints of the microbial communities, a computer-assisted cluster analysis was carried out using the GelCompar II ® software (Applied Math, Kortrijk, Belgium).Computer-assisted pattern analyses of molecular fingerprints of the communities resulted in correlation matrices that were then subjected to a significance test by applying 10 4 random permutations of sample elements [39].The relationships between samples taken from different treated crops were investigated.

Conclusions
The bait plant-based selection with subsequent stress protection screening assays resulted in target-adapted microorganisms.The resulting SPAs from this selection strategy are taxonomically highly diverse and belong primarily to the Pseudomonas group (e.g., P. asplenii, P. brassicacearum, P. corrugata, P. kilonensis, P. otitidis, P. poae).They (except plant and human pathogenic species) can be further implemented into commercial crop production for enhanced plant health and stress tolerance, as the first field results are promising, and the applied evaluation score allows for a flexible and adjustable constitution of SPAs independent of their supposed biotechnological application.In using different bait plants, environmental conditions or microenvironments other than the rhizosphere, new implementations can be discovered in regards to the changing climate conditions or different soil types, e.g., desert soil, soils with extreme pH values or contaminated soil.Furthermore, the system can be used to "fine tune" highly specific host plant species and even cultivars.

Figure 3 .
Figure3.Unweighted pair group method with arithmetic mean (UPGMA) dendrogram of bacterial rhizosphere communities of moss-primed and untreated crop plants in two replicates (first number) with three repetitions per plant species (second number).The dendrograms were generated from SSCP community profiles using GelCompar II ® .The following settings were used: dendrogram type: UPGMA; similarity coefficient: curve-based Pearson correlation; optimization: 4%; position tolerance: 1%.The dashed line indicates a similarity level at 90%.Letter code: first letter of the bait plant species; M, maize; R, oilseed rape; S, sorghum; Z, sugar beet; the second letter is the treatment; M, moss; C, control.

Figure 4 .
Figure 4. Phylogenetic composition of the antagonistic community selected by bait plants.The tree is comprised of 47 partial 16S rRNA gene sequences.Neighbor-joining phylogenetic tree (maximum composite likelihood method) was constructed with both one representative sequence per major species and the nearest reference sequences (accession numbers in brackets).The reliability of the tree topology was evaluated using bootstrap analysis with 10,000 re-samplings.Color code: rhizosphere samples from maize, yellow/orange/red; oilseed rape, green shades; sorghum, blue shades; sugar beet, purple shades.Sequences of MF2-1-4 (KA3) and ZM1-1-3 (KA90) were excluded from MEGA5 analysis due to their short sequence length.

Table 2 .
Score scheme results with maximum numbers of positive reactions for each characterization group.a Strain number with color code: rhizosphere samples from maize, yellow/orange/red; oilseed rape, green shades; sorghum, blue shades; sugar beet, purple shades.The letter code is the first letter of the crop: M, maize; R, oilseed rape; S, sorghum; Z, sugar beet.The second letter is the treatment: F, lichen; M, moss; P, primrose.AB, antibiotic; n.d., not determined; ROS, reactive oxygen species.ARDRA, amplified ribosomal RNA gene restriction analysis.

Table S2 .
Sampling details: treated plants.Treatments: C, control; F, lichen; M, moss; P, primrose.Figure S1.BOX-PCR fingerprints of single isolates selected after priming seeds of maize, oilseed rape, sorghum and sugar beet with moss, lichen and primrose communities.UPGMA dendrogram obtained by GelCompar II ® ; Pearson correlation with 5% optimization.