Screening and Characterization of Potentially Suppressive Soils against Gaeumannomyces graminis under Extensive Wheat Cropping by Chilean Indigenous Communities

Wheat production around the world is severely compromised by the occurrence of “take-all” disease, which is caused by the soil-borne pathogen Gaeumannomyces graminis var. tritici (Ggt). In this context, suppressive soils are those environments in which plants comparatively suffer less soil-borne pathogen diseases than expected, owing to native soil microorganism activities. In southern Chile, where 85% of the national cereal production takes place, several studies have suggested the existence of suppressive soils under extensive wheat cropping. Thus, this study aimed to screen Ggt-suppressive soil occurrence in 16 locations managed by indigenous “Mapuche” communities, using extensive wheat cropping for more than 10 years. Ggt growth inhibition in vitro screenings allowed the identification of nine putative suppressive soils. Six of these soils, including Andisols and Ultisols, were confirmed to be suppressive, since they reduced take-all disease in wheat plants growing under greenhouse conditions. Suppressiveness was lost upon soil sterilization, and recovered by adding 1% of the natural soil, hence confirming that suppressiveness was closely associated to the soil microbiome community composition. Our results demonstrate that long-term extensive wheat cropping, established by small Mapuche communities, can generate suppressive soils that can be used as effective microorganism sources for take-all disease biocontrol. Accordingly, suppressive soil identification and characterization are key steps for the development of environmentally-friendly and efficient biotechnological applications for soil-borne disease control.


INTRODUCTION
Currently, the main concerns in modern agriculture are those associated with the effect of climate change on biotic and abiotic parameters in arable soils, and their impact on crop yields and food supply at global level (Soussana et al., 2010). In this context, several authors have pointed out an increase of soil-borne disease incidence in winter cereals as a consequence of climate change (French et al., 2009;Manici et al., 2014). Similarly, these studies have also described that climate change may contribute to certain soil-borne pathogen migration toward niches or regions previously uncolonized by these pathogens (French et al., 2009).
The southern region of Chile produces around 85% of cereals, where 40% is wheat (Triticum aestivum L.; ODEPA, 2016). However, wheat production is frequently reduced by the incidence of "take-all" disease, which is caused by the soil-borne pathogen Gaeumannomyces graminis var tritici (Ggt; Andrade et al., 2011), causing the highest wheat crop losses in Chile (Moya-Elizondo et al., 2015). Soil-borne pathogen incidence in cereal cropping is difficult to control due to their natural persistence in soils and the inefficiency of chemical controls (De Coninck et al., 2015); thus, biological control becomes a very promising alternative to prevent soil diseases. The incidence of take-all disease outbreak is favored under particular soil conditions, called "conducive" soils (Chng et al., 2015). In contrast, "suppressive" soils occur as a natural phenomena, preventing soil-borne pathogen establishment or reducing disease incidence (Jara et al., 2011) even in the presence of a susceptible host plant and favorable soil or climate conditions. In this context, several studies have shown that native soil microorganism activity can be pivotal in Ggt disease suppression (Weller et al., 2002;Cook, 2003;Mendes et al., 2011Mendes et al., , 2013. Ggt (and other soil-borne pathogens) suppressive soils have been reported and characterized around the world (Bull et al., 1991;Bithell et al., 2012;Chng et al., 2013).
In the twentieth century, agrarian policies resulted in the establishment of numerous small farmers practicing extensive agriculture in southern Chile, particularly in the "Region de La Araucania" (38 • 54 ′ 00 ′′ S; 72 • 40 ′ 00 ′′ S) (Clapp, 1998). Most farmers belong to the "Mapuche" ethnic group, which represents 50% of the total population in La Araucania. The Mapuche community is characterized by the use of ancestral agronomic techniques to produce their own agricultural products without the application of chemicals such as, commercial fertilizers and pesticides. The extensive Mapuche agriculture is mainly directed to family group subsistence. Under this scenario, we hypothesize that this long-term land-use could represent a natural and effective source of suppressive soils against soil-borne pathogen diseases. Thus, suppressive soils may also be relevant in the context of alterations related to soil-borne pathogen incidence and migration that have been predicted by climate change, and the deleterious effect on intensive agro-chemical product use in arable soils (Meza and Silva, 2009;Neuenschwander, 2010). In fact, Andrade et al. (2011) detected five Ggt-suppressive soils located in La Araucania with a long history of monoculture and natural pasture. Similarly, Arismendi et al. (2012) reported the presence of Pseudomonas fluorescens strains, which are able to produce 2, 4-DAPG, a known biocontrol compound involved in soil Ggt suppression Weller et al., 2012), in 13 soils from La Araucania and Los Lagos regions. Moreover, we have recently isolated and characterized four endophytic bacteria from wheat plants in this area (Acinetobacter sp. E6.2, Bacillus sp. E8.1, Bacillus sp. E5 and Klebsiella sp. E1), which are able to inhibit Ggt mycelia growth in vitro, promote plant growth, and diminish take-all disease incidence under greenhouse experiments (Durán et al., 2014). Although most studies on Ggt suppression have focussed on bacteria, there are some reports showing that some fungal strains can also reduce Ggt incidence (Macia-Vicente et al., 2008). These examples illustrate the great potential of native microorganisms as soil inoculants able to increase plant growth and prevent soil-borne pathogen disease incidence in cereal cropping. Despite of this potential, studies focusing on Ggt-suppressive soils occurrence, and their characterization, are still very limited.
The main objective of our study was to screen Ggtsuppressive soils in 16 locations managed by indigenous Mapuche communities from La Araucania, southern Chile, where extensive wheat cropping was applied for more than 10 years. Soil chemical properties and microbial community composition were characterized. Furthermore, soils were categorized as Ggt suppressive or conducive based on their potential for pathogen inhibition in vitro, and their efficacy in controlling take-all disease under greenhouse conditions.

Soil Sampling
In southern Chile, acidic volcanic soils (Andisols and Ultisols) are the predominant soil types supporting the bulk of agricultural and forestry production. Andisols include modern and recent ash deposits, and Ultisols correspond to ancient deposits (Andrade et al., 2011). Soil samples were collected from extensive wheat cropping areas managed by Mapuche communities in 16 locations from La Araucania ( Figure 1A; Table 1). Nine samples (number 3, 4, 5, 7, 10, 12, 13, 14, 15, and 16) were taken from soils with a long rotation history between wheat monoculture and natural pasture for more than 10 years, and six samples (number 2, 6, 8, 9, 11, and 12) were taken from soils with wheat monoculture including oat rotation. An additional sample number 1 was taken from a Ggt-conducive soil with wheat, including clover rotation, and used as positive control. Parent material from soils 1 to 10 was classified as Andisol, whereas 11-16 were classified as Ultisols.
Samples were collected from rhizosphere and bulk soil at 0-20 cm depth, and then stored in a cold room at 5 • C until usage. Soil sample chemical properties were determined as follow. Briefly, available P (P Olsen ) was extracted by using 0.5 M Na-bicarbonate and analyzed by using the molybdate method (Murphy and Riley, 1962). Organic matter contents were estimated by wet digestion (Walkley and Black, 1934). Soil pH was measured in 1:2.5 soil/deionized water suspensions. Exchangeable potassium (K + ), calcium (Ca 2+ ), magnesium (Mg 2+ ), and sodium (Na + ) were extracted with 1 M ammonium acetate (CH 3 COONH 4 ) at pH 7.0 and analyzed by flame atomic adsorption spectrophotometry Triangle colors indicate type of soil: andisol (green) and ultisol (red).  (Warncke and Brown, 1998). Exchangeable aluminum (Al 3+ ) was extracted with 1 M KCl and analyzed by FAAS (Bertsch and Bloom, 1996). All samples analyses were made in triplicate. To group and determine significant differences between samples based on their chemical properties, data were imported into the PRIMER 7 software (PRIMER-E Ltd, Ivybridge, UK), transformed and normalized using square-root followed by a log (Xþ1) transformations (Lee et al., 2012). Then, a distance matrix was generated based on Euclidean distances and samples were grouped by hierarchical clustering (group average), and then visualized by principal component analysis (PCA).

Bacterial Community Composition and Ggt Detection in Soil Samples
Bacterial community compositions in rhizosphere and bulk soil samples were examined by denaturing gradient gel electrophoresis (DGGE), according to Iwamoto et al. (2000). Total genomic DNA was extracted from 0.5 to 1 g of soil samples using the PowerSoil R DNA Isolation Kit (Mo-Bio Laboratories, Carlsbad, CA, USA), according to manufacturer's instructions. The 16S rRNA gene fragments were amplified by touchdown PCR, using EUBf933-GC/EUBr1387 primer set (Iwamoto et al., 2000). DGGE analysis was performed using a DCode system (Bio-Rad Laboratories, Inc., USA). The PCR product (20 µL) was loaded onto a 6% (v/v) polyacrylamide gel with a 40-70% gradient (urea and formamide). Electrophoresis was run for 12 h at 100 V. Banding profiles were visualized by staining the gel 1:10.000 (v/v) with SYBR Gold (Molecular Probes, Invitrogen Co., USA) for 30 min, followed by image capture using GelDoc-It TS2 Imager (UVP, Upland, CA, USA). DGGE banding profile clustering, using a dendrogram, was also carried out by using Phoretix 1D analysis software (TotalLab Ltd., UK). The correlation between bacterial communities (biological parameters) and chemical soil properties (ecological parameters) was visualized by non-metric multidimensional scaling (MDS) analysis using Primer 7 + Permanova software (Primer-E Ltd., Ivybridge, UK) (Clarke, 1993). The in silico analysis was also used to estimate the bacterial diversity by richness (S), Shannon-Wiener, and dominance by Simpson Index (D), represented by 1-D or 1-λ (Sagar and Sharma, 2012). Ggt occurrence in soil samples: DNA extracts from soil samples were subjected to PCR using the primer set NS5 (5 ′ -AAC TTA AAG GAA TTG ACG GAA G-3 ′ ) and GGT-RP (5 ′ -TGC AAT GGC TTC GTG AA-3 ′ ) designed by Fouly and Wilkinson (2000) specifically for Ggt. The PCR conditions were as follow: an initial denaturation at 93 • C for 3 min, followed for 93 • C for 1 min, 52 • C for 1 min, and 72 • C for 1 min to 35 cycles, and finally with 72 • C for 5 min. All soil samples were tested in triplicate, and pure G. graminis var tritici (Andrade et al., 2011) and Aspergillus niger DNA (code CCT-UFRO 15.62), obtained from La Frontera University Type Strain Culture Collection (http://ccct.ufro.cl/), were used as positive and negative controls, respectively. The presence of Ggt in soil samples showing positive Ggt reaction was also confirmed by the presence of take-all disease symptoms, chlorosis, and blackening roots in wheat seedlings in the pot assays.

Putative Ggt-Suppressive Soil Screening
Putative Ggt-suppressive soil screening was performed by two in vitro inhibition tests using rhizosphere soil samples as follows:

In vitro Inhibition Test on Solid Media
A first screening was carried out in order to evaluate rhizosphere soil capability of inhibiting Ggt growth on agar plates (Supplementary Figure 1A). Briefly, the Ggt inoculum was produced by growing the fungus on potato dextrose agar (PDA) medium at 25 • C for 1 week. Agar disks (4-mm diameter) containing Ggt were aseptically incised and transferred to the center of agar plates containing fresh Luria Bertani (LB) and PDA medium (proportion 1:1). A hole of 10 mm was performed in the agar medium at a distance of 3 cm from Ggt inoculum, and 0.05 g of rhizosphere soil were placed in the agar hole. Ggt mycelia growth was registered at 3, 5, and 7 days of incubation at 25 • C in the darkness, as described by Liu et al. (2011). A fraction of all soil samples was sterilized and samples included in the agar as negative controls. All tests were carried out in triplicate.

In vitro Inhibition Test in Aqueous Soil Extracts
Because some soils contained elevated loads of microorganisms affecting fungal measurements on agar medium (categorized as un-determined in the in vitro test described in Section In vitro Inhibition Test on Solid Media), a second assay was performed in tubes with rhizosphere soil extracts (Supplementary Figure 1B). Briefly, 1 g of rhizosphere soil sample was suspended in 9 ml of sterile phosphate buffer saline (PBS; pH 7.4), and sonicated at 60% amplitude for 30 s. Then, 1 mL of supernatant was added to an eppendorf, and then inoculated with 1% Ggt inoculum. Soil extract tubes were incubated at room temperature for 3, 5, and 7 days, and fungal growth was estimated by quantification of fungal biomass by crystal violet (CV) staining as follows. After incubation, soil extract samples were washed with distilled water and fixed with 500 µl methanol for 15 min at room temperature. Later, they were centrifuged at 13,000 rpm × 1 min, supernatants were discarded, and tubes were air-dried; 400 µl of CV was added to each tube and incubated for 5 min. Tubes were washed three times with distilled water. Finally, 400 µl of acetic acid (33% v/v) were added and kept in the tubes for 5 min. The absorbance of the obtained solution was determined in triplicate in a multi-plate reader at 590 nm (Silva et al., 2009, Supplementary Figure 1B).
Nine putative suppressive soils from the inhibition tests, in both agar plates and soil extract tubes, were used for the greenhouse assay, using soil 1 as Ggt-conducivepositive control.

Inoculum Preparation
The characterization of the pathogen as Ggt was done based on the sequencing of ribosomal internal transcribed spacer 2 region (ITS2). ITS2 was amplified by touchdown polymerase chain reaction (PCR) with primer sets fITS9 (5 ′ -GAACGCAG CRAAIIGYG-3 ′ ) and ITS4 (5-′ TCCTCCGCTTATTGATATGC-3 ′ ), using the following conditions: an initial denaturation at 95 • C for 3 min, followed by 25 cycles-each at 95 • C for 30 s, with an annealing step with a 0.5 • C decrease-each cycle from 65 • C to 52.5 • C, and extension at 72 • C for 30 s. Twenty-five additional cycles were then carried out at 95 • C denaturation for 30 sec, 55 • C annealing, primer extension at 72 • C for 30 s, and a final extension step of 7 min at 72 • C. The PCR products were purified and sequenced by Austral-Omics (Universidad Austral of Valdivia-Chile). The sequence was compared with those in the GenBank database, fungal identity was confirmed (99%), and then deposited under accession no. KY689233.
The Ggt inoculum was prepared as follows: oat kernels were soaked in water for 24 h and sterilized for 3 consecutive days. Then, Ggt pathogenic isolate discs were grown on PDA for 7 days, put on the sterile oat, and maintained at room temperature for 20 days. Colonized oat kernels were blended, sieved to a particle size of 0.5-1.0 mm, and stored at 5 • C until usage (Andrade et al., 2011).

Greenhouse Assay
Nine putative suppressive soils (number 2, 3, 4, 6, 11, 13, 14, 15, and 16), selected according to results obtained in the in vitro tests (Section Putative Ggt-Suppressive Soil Screening), were tested in their ability to suppress take-all disease in planta under greenhouse conditions. Plastic containers enclosing 200 g of soil were used in quintuplicate. Wheat seeds Otto cv were surface sterilized (15% ethanol plus 1% sodium hypochlorite for 2 min) and 5 seeds were used for each treatment. Plants were watered every 3 days, and Taylor and Foyd nutrient solution (Taylor and Foyd, 1985) was applied each 15 days.
The experimental design consisted in sterile (heat-treated) rhizosphere soil to determine disease level incidence, discarding, or diminishing the effect of soil microorganisms; untreated, air-dried soil (natural) to determine the effect of soil native microorganisms; and, sterile soil +1% untreated, natural soil (ster+1% nat) to assess suppression transferability to sterile soils. Suppression transferability has been shown to occur when natural soils are added in as low as 1% (v:v) to non-suppressive soils, as reported earlier (Shipton et al., 1973;Andrade et al., 1994). Ggt inocula were applied at 0.1% in relation to soil weight (2 g), and all treatments were also performed in soil without Ggt inocula as controls. After 40 days, plants were carefully removed from the soil, weighted and the root blackening percentage was determined.

Ggt Presence in Plant Tissues
The Ggt presence was evaluated in roots of PCR infected and non-infected plants. Roots were individually assessed for infection and root blackening percentage was evaluated against a white background. Shoots were carefully separated from the roots, placed into individual paper envelopes, and dried at 70 • C for 72 h, to obtain shoot dry weight. In order to confirm Ggt infection, total DNA from wheat infected tissues was extracted with soil DNA Isolation Kit (Ultraclean, Mo-Bio Laboratories) according to manufacturer's instructions. Specific Ggt DNA fragments were amplified by using NS5 and GGT-RP primer sets, as described in Section Bacterial Community Composition and Ggt Detection in Soil Samples. Pure DNA extracts from Ggt and A. niger collection strains were used as positive and negative control, respectively.

Presence of 2,4-Diacetylphloroglucinol-Producing Bacteria
The presence of 2,4-DAPG-producing bacteria was also evaluated by PCR. Total DNA from rhizosphere soil was extracted with soil DNA Isolation Kit (Ultraclean, Mo-Bio Laboratories) according to manufacturer's instructions. Specific primer sets B2BF (5 ′ ACC CAC CGC AGC ATC GTT TAT GAG C-3 ′ ) and BPR4 (5 ′ CCG CCG GTA TGG AAG ATG AAA AAG TC-3 ′ ), which target the phlD gene (encoding a polyketide synthase that synthesizes monoacetylphloroglucinol, the precursor to 2,4-DAPG that is essential for the phloroglucinol biosynthesis) were used in the PCR reaction . PCR conditions were: an initial denaturation at 95 • C for 3 min, followed by 35 cycles each at 95 • C denaturation for 1 min, 60 • annealing for 1 min, 72 • C extension for 1 min, and final extension step for 10 ′ at 72 • C. Pure DNA extracts from Pseudomonas spp. (SA 32A) and Enterobacter spp. (RJAL6) strains were used as positive and negative controls, respectively (Mora et al., 2017).

Bacterial Community Structure in Suppressive Soils
The bacterial community composition in the rhizosphere suppressive soils from the greenhouse assay was examined by DGGE as described in Section Bacterial Community Composition and Ggt Detection in Soil Samples. The similarity between bacterial communities was visualized by non-metric multidimensional scaling analysis (MDS), using Primer 7 software (Primer-E Ltd., Ivybridge, UK), which showed a Bray-Curtis similarity index higher than 50% and 0.14 stress values (Clarke, 1993).

Statistical Analyses
Data normality was analyzed according to Kolmogorov's test. Data obtained in Section In vitro Inhibition Test on Solid Media (in vitro plate assay) were analyzed by a one-way analysis of variance (ANOVA) and compared by Tukey test, using SPSS software (SPSS, Inc.). Comparisons between inoculated and noninoculated samples from screening 2 were made, and Student t-test was used for related samples with 95% confidence interval. For the greenhouse assay multivariate analysis of variance (MANOVA) and comparisons were carried out for each set with Tukey's test by SPSS software (SPSS, Inc.). Values were given as means ± standard errors. Differences were considered significant when the P value was lower than or equal to 0.01. The microbial diversity analysis was described above.

Bacterial Community Composition is Related to Chemical Soil Composition and Differs in Bulk and Rhizosphere Soils in Terms of Dominance and Diversity
Rhizosphere community composition is highly related to soil classification as revealed by MDS, in which microbiology and soil chemical/environmental properties were analyzed (Figure 2A). Samples from Andisol and Ultisol were notably different from each other, with 55% similarity. According to the Spearman correlation, Ultisol soils were more significantly related with chemical parameters than Andisol soils (r 2 = 0.90 and r 2 = 0.59, respectively, Table 3). Bulk soils were also separated according to the chemical soil composition, but this difference was not significant, and all samples were grouped at 55% (Figure 2B). In TABLE 2 | Average values for some chemical properties of rhizosphere soils used in this study.
Frontiers in Microbiology | www.frontiersin.org  this sense, both Andisol and Ultisol were equally related with soil chemical parameters (r 2 = 0.83, Table 3). Differences in bacterial community structures between bulk and rhizosphere soil were revealed by MDS analysis, based on DGGE banding profiles. In relation to bulk soils, the nMDS analysis revealed the existence of major groups at 40% similarity formed by rhizosphere soils indistinctly grouped, all correlated among them. In general, bulk soils were separated into two main groups formed mainly by Andisol (soil 1, 2, 3, 4, 5, 6, 7) and Ultisol soils (11,12,13,14,15,16) (Figure 3).
Regarding the microbial diversity in bulk and rhizosphere soils (Supplementary Figure 2), in general, the Shannon index (H') showed values <2.0 for bulk (except soils 5 and 6 from Quilaco) whereas it reached values ∼2.5 for rhizosphere soil, indicating a lower diversity in bulk than in rhizosphere soils. A similar trend was observed in the case of species number (S). The lowest diversities (<1.0) and species richness (S) were observed in bulk soils from Membrillar (soils 9 and 10) whereas highest diversity values were obtained in rhizosphere soils from Quilaco (soil 5 and 6). As for bulk soils, rhizosphere soils taken from Membrillar also showed lower diversity values (≤2.0). In contrast, less dominance values in rhizosphere soils compared with bulk soils were also observed by Simpsons (D) index represented by 1-D, in which the lower values indicate major dominance of species. Thus, samples 9 and 10, which showed less diversity according to the Shannon index, also showed major dominance of species. Therefore, bulk soils showed lower biodiversity and major dominance compared to rhizosphere soils (Supplementary Figure 2), and we observed a direct Pearson correlation between both indexes (P < 0.01, data not shown).

Screening for Putative Ggt-Suppressive Soils
According to the in vitro inhibition test on solid media ( Figure 4A and Supplementary Figure 1A), four soils from Las Cardas (soils 2 and 3), Boyeco (soil 11), and Llufquentue (soil 16) were considered as suppressive against Ggt, since Ggt growth was significantly inhibited when compared with the positive control. However, several soils could not be properly evaluated based on this assay and were classified as undetermined. A second evaluation, using the in vitro inhibition test in aqueous soil extracts ( Figure 4B and Supplementary Figure 1B) suggested the presence of nine Ggt-suppressive soils. They were collected from Las Cardas, Momberg,Quilaco,Boyeco,and Llufquentue (soils 2,3,4,6,11,13,14,15,and 16,respectively). In line with the assay on solid agar, this test also showed that soils collected from Membrillar (soils 7, 9, 10, and 12) were Ggt conducive. This assay also revealed the presence of other Ggt-conducive soils (soils 5 and 8), including the positive control soil 1. In summary, based on in vitro experiments, 9 soils were determined as suppressive, and 7 as conducive ( Figure 4B).

Incidence of Take-All Disease
Firstly, we aimed to establish a positive control for take all decline and soil conduciveness for our suppressive soil assays. Plants growing in soil 1 showed clear disease symptoms including leaf chlorosis and blackening root, distinctive take-all disease symptoms. We confirmed the presence and the identity of Ggt in the rhizospheric soil, extracting the total DNA and using Ggt specific primers. The PCR amplification confirmed that soil 1 was positive in terms of Ggt presence (Supplementary Figure  3). Then, to confirm the ability to suppress take-all disease in the putative nine rhizosphere soils selected through the in vitro assays (Figure 4), an assay with the pathosystem wheat-Ggt was performed under greenhouse conditions by growing wheat in Ggt-inoculated and non-inoculated soils. Take-all disease symptoms were evident in all inoculated soils, except in soil 11 (data not shown). The origin of the symptoms was confirmed by Ggt DNA amplification with Ggt specific primers of infected plants; the amplification band was present in all Ggt-inoculated rootswhereas it was absent in the non-inoculated plants (Supplementary Figure 4). In order to determine whether Ggt suppressiveness was due to microbial community or to soil physico-chemical characteristics, sterilization by heat treatment of the soils was performed. The plants grown in the greenhouse on heat-treated sterilized soils showed higher disease symptoms than in the corresponding natural soils (Figure 5): the percentage of blackening root for plants growing in sterile soils ranged between 10 and 40%, while only between 3 and 10% for natural FIGURE 4 | In vitro assays to determine the inhibition of Gaeumannomyces graminis growth in the presence of rhizosphere soil samples. (A) Shows Ggt mycelial growth on solid media supplemented with the different rhizosphere soils after 3, 5, and 7 days of incubation. (B) Shows Ggt fungal biomass in liquid rhizosphere soil extracts after 3, 5, and 7 days of incubation. Asterisks indicate significant differences (Tukey's test, p ≤0.05) in relation to the positive control (solid media assay) and negative control (liquid assay). soils, except for soils 1, 6, and 14 ( Figure 5). In fact, no significant differences were found between sterile and natural soil in the case of soils 1 and 14 (blackening roots levels ranging from 15 to 35%).
A negative correlation between blackening roots and biomass was also observed. Thus, the lower biomass was found in the sterile Ggt-inoculated soil treatments, where plants showed higher plant infection, with the exception, again, of soils 6 and 14, in which plant biomass was similar in sterile and natural soils (Figure 6). In treatments where a 1% of the natural soils was added to the sterile soils, biomass production increased significantly (P ≤ 0.05). In fact, in soils 4, 13, and 15 the highest biomass for this treatment was found, being similar in the rest of the soils. Only in soil 1 (Ggt-conducive positive control) plant biomass was lower in the 1% natural soil supplemented treatment.
Presence of 2,4-DAPG Producing Bacteria 2,4-DAPG-producing bacteria have been reported as major takeall disease suppressors in soils . Accordingly, we checked for their presence in our selected suppressive soils. However, we found that only soil 15 was positive for the phlD gene, essential in the DAPG biosynthetic pathway . The rest of suppressive soil samples did not show the presence of amplicons for this gene (Supplementary Figure 5), FIGURE 5 | Take all disease symptoms (root blackening) in wheat plants from the suppressive soil bioassay in greenhouse. Treatments were: natural soil (natural), sterile soil (sterile), and sterile soil supplemented with 1% of natural soil (ster+1% nat), inoculated (+) or not (−) by Ggt (n = 5). The disease index scale used is represented by the pictures in the bottom panel. Tukey's test was used to compare treatments means, values followed by the same letter do not differ at P ≤ 0.05 (n = 5). Green bars represent suppressive soils, red bars conducive soils, and yellow bars represent undeterminated soil.
suggesting that other metabolic pathways should be responsible for suppressiveness in those soils.

Rhizosphere Bacterial Community Composition in Suppressive and Conducive Soils
Regarding the composition of the rhizosphere bacterial communities (as revealed by DGGE analysis), the dendrogram showed differences between suppressive and conducive soils in the case ofAndisols (Figure 7). These results were confirmed by using MDS analysis, showing clear differences between the conducive control soil (soil 1) and the rest of treatments in Andisol (55% similarity). However, in the case of Ultisols in which the bacterial communities were more significantly related with the soil chemical parameters than in Andisol soils (Table 3), this tendency was not observed (Figure 7B).

DISCUSSION
Several studies have described the importance of soil microorganisms as suppressive agents against phytopathogens worldwide, including soils with chemical properties similar to Chilean agricultural soils, such as, New Zealand soils (Bithell et al., 2013;Chng et al., 2013Chng et al., , 2015Perez et al., 2016). However, despite of the great potential that Chilean suppressive soils offer in terms of microbial diversity and potential for the development of biocontrol strategies, their studies are extremely limited (Andrade et al., 2011). This topic acquires special relevance FIGURE 6 | Shoot biomass (dry weight) of wheat plants from suppressive soil assay in greenhouse. Treatments were: natural soil (natural), sterile soil (sterile), and sterile soil supplemented with 1% of natural soil (ster+1% nat)-and inoculated (+) or not (−) by Ggt (n = 5). Tukey's test was used to compare treatments means, values followed by the same letter do not differ at P ≤ 0.05 (n = 5). when ancestral extensive agriculture has been applied for a long term by native communities; Mapuche people cultivate in small areas to produce their own agricultural products using rustic metal tools and low inputs (Montalba-Navarro, 2004). This ancestral production could be replaced by more aggressive and intensive modern agriculture techniques, with the possible subsequent loss of diversity and suppressive potential of soil microbial communities. Under this scenario, it is essential to identify suppressive soils in these areas. In this work we screened the ability of different agricultural soils from southern Chile to suppress take-all disease in wheat plants. Suppressive soils were collected from little farms mainly belonging to indigenous Mapuche communities from Southern Chile, who practice monoculture and cultivate small subsistence wheat areas.

Soil Chemical Parameters and Their Relation with Microbial Community Composition
Our results on soil chemical properties showed that 25% of the studied soils presented a low content of available P (< 20 mg kg −1 ), moderate acidity (pH < 5.5), and high Al saturation (>10%), which are main characteristics of agricultural soils from southern Chile (Mora et al., 2006). In general, Andisols grouped together when considering the chemical parameters, whereas Ultisols were more diverse. However, when we compared microbial diversity we observed a direct correlation with soil chemistry mainly in rhizospheric soils, and more significantly for Ultisols. Similar results were found by Smalla et al. (2007) despite the amplified fragments comprised different variable regions and lengths, DGGE, T-RFLP, and SSCP analyses led to clustering of fingerprints, which correlated with soil physico-chemical properties. In Chilean Andisol soils, Jorquera et al. (2014) also showed that soil chemistry influenced the composition of rhizobacterial communities, and in Ultisol soils from China, Li et al. (2017) found through nMDS analyses that microbial communities also correlated with soil chemical parameters and fertilization strategies. Although there is some similarity in microbial composition between rhizosphere and bulk microbial composition (de Ridder-Duine et al., 2005), in our study rhizosphere soils were more diverse in terms of richness and dominance than bulk soils. In fact, it is known that the amount of microorganism around the rhizosphere is 10-to 1,000fold higher than that found in bulk soil due to rhizodeposition (Doughari, 2015;Glick, 2015).

Suppressive Soil Identification in Southern Chile
In terms of potential take-all disease suppression, in vitro Ggtgrowth inhibition tests with 16 different soils allowed to identify 9 potentially suppressive soils, and 6 of them were later confirmed to be suppressive in plant bioassays under greenhouse conditions (soil 2, 3, 4, 13, 15, and 16). In general, these soils were cultivated with wheat monoculture and natural pasture for more than 10 years, a management similar to that described previously for other suppressive soils in South Chile (Andrade et al., 2011). Regarding the timing, take-all suppression appeared after 4-6 years of wheat monoculture (McSpadden Gardener and Weller, 2001), and even later, although they can also occur in soils with 3-4 years of monoculture under relatively high pathogen concentrations (Chng et al., 2015). In fact, early studies by Baker and Cook (1974) showed that 3 years of successive wheat cropping could be sufficient for the development of specific suppression. As an exception, soil 2, withrotation based in wheattriticale and oat for more than 10 years, was also found to be Ggt-suppressive. It is well known that oat roots produce saponin avenacin, a glycosylated triterpenoid secondary metabolite with antifungal properties that has been involved in determining oat resistance to soil fungal pathogens (Osbourn et al., 1994;Freeman and Ward, 2004). On the other hand, wheat plants grown in soil 11, characterized by a very high humidity, showed no symptoms of Ggt infection in roots. This is in agreement with earlier reports by Nish (1973), who studied Ggt survival in the field under controlled conditions, and showed a significant reduction in Ggt incidence in wet cool soils.
In greenhouse assays, wheat plants grown on sterile Ggtinoculated soils presented the highest disease incidence and the lowest biomass production, suggesting the relevance of the autochthonous soil microbial communities on plant growth promotion and disease suppression. The effect observed in sterile soils-increased susceptibility and reduced biomasswas improved when they were supplemented with 1% of the same natural soils, except in the case of the conducive soils 1 and 14. This improvement also confirms the role of the microbial communities originally present in the take-all disease suppressive soils and Ggt control. This characteristic is known as specific suppression (Cook, 2003;Andrade et al., 2011;Chng et al., 2015). It is noteworthy that we did not find any relation between suppressive soils and the presence of 2,4-DAPG-producing bacteria since its presence was only detected in one out of the six confirmed suppressive soils. The low occurrence of 2,4-DAPG determined by phlD gene, suggests that other mechanism(s) or antifungal compound(s) are synthetized by native soil microorganisms that could contribute to the effective biocontrol against Ggt. Thus, strains related to soil suppressiveness seem to be differentially shaped by multiple soil factors (Imperiali et al., 2017). Therefore, further studies are required to identify the mechanisms involved in Ggt disease suppression.

Suppressive Soils and Microbial Community Composition
As mentioned above, the selected suppressive soils did not group together in relation to their chemical properties and geographical origin. However, when bacterial communities were analyzed by DGGE, within Andisols, suppressive soils grouped separately with respect to the control conducive soil 1, but not for Ultisols, in which suppressive soils grouped together with soil 14, also classified as conducive. This could be attributed to the high relation between Ultisol soils and soil chemical parameters when compared to that relation in Andisol soils (r 2 = 0.90 and 0.59, respectively). Future research should explore which specific microbial groups act directly upon Ggt suppression and how rhizosphere microbial communities are selected and regulated by the plant rhizosphere, especially in the presence of the phytopathogens. Moreover, identifying the bacterial groups and their antagonist mechanisms, as well as exploring the potential stimulation of plant defense mechanisms are pivotal for the development of biocontrol strategies based on the use of suppressive soils. With this aim, future research should tackle the multifactor soil-microbiome-plant-pathogen systems, considering not only direct antifungal activities in the rhizosphere, but also potential stimulation of plant defense and microbe-selection mechanisms.

CONCLUSIONS
Suppressive soils represent an important microbial source for the biocontrol of soil-borne pathogens, and their identification and characterization is crucial since many of these soils may be lost by the increase of intensive agriculture practices worldwide. Here, we identified six suppressive soils against take-all disease, which have been managed under ancestral and rudimental agronomic practices by Chilean indigenous Mapuche communities. Then, we showed that suppressive activity in the tested soils correlated with the microbial community composition and not with the chemical properties and geographical origin of the studied soils. The key role of the soil microbial communities in Ggt suppression was confirmed in assays with sterile (heat-treated) suppressive soils where Ggt suppresiveness was completely lost, and recovered again through the addition of 1% of the corresponding non sterile natural suppressive soil. Our understanding of microbial communities in suppressive soils as well as the mechanisms acting in disease suppression in the rhizosphere must be considered as a valuable tool for the development of sustainable control of soil-borne pathogen (such as, take-all disease) in agriculture.

AUTHOR CONTRIBUTIONS
PD, SV, MJ, MM designed the research. MP and VC supervised the study. PD organized the soil sampling and chemical soil analysis. SV performed the greenhouse experiments. VC analyzed the data. PD and MP wrote the manuscript and authors critically revised the manuscript and approved the final version.

FUNDING
This study was supported by the Comisión Nacional de Investigación Científica and Tecnológica (CONICYT), FONDECYT Iniciation Project No. 11150540 from Chilean Government. MJ thanks the support by FONDECYT No. 1160302. We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).