Bacillales Members from the Olive Rhizosphere Are Effective Biological Control Agents against the Defoliating Pathotype of Verticillium dahliae

The use of biological control agents (BCAs) is of interest within an integrated management strategy of Verticillium wilt of olive (VWO) caused by the soil-borne fungus Verticillium dahliae Kleb. Previous studies have shown that the root/rhizosphere of healthy olive plants is an important reservoir of microorganisms displaying biocontrol activity against VWO (i.e., Pseudomonas strains PICF7 and PIC141). Moreover, these BCAs are already adapted to the ecological niche where they are deployed. Three novel bacteria (strains PIC28, PIC73 and PIC167) from nursery-produced olive plants were in-depth characterized using a previously implemented approach consisting of in situ isolation, in vitro antagonism tests, in planta bioassays, phenotypic and metabolic characterization, genome analyses and in silico identification of traits involved in plant-bacteria interactions, and multi-locus sequence analyses. All strains displayed in vitro growth inhibition of different olive pathogens and biocontrol effectiveness against Verticillium dahliae, with strain PIC73 being the most effective BCA. Strains PIC73 and PIC167 were identified as Paenibacillus polymyxa (Prazmowski) Ash et al. and Paenibacillus terrae Yoon et al., respectively. Strain PIC28 belongs to the Bacillus genus. Some of these Bacillales members showed in vitro compatibility with previously characterized BCAs (Pseudomonas spp. strains) also originating from the olive rhizosphere, paving the way for the future development of tailored bacterial consortia effective against VWO.


Introduction
Verticillium wilt of olive (VWO) is one of the most serious diseases affecting olive (Olea europaea L. subsp.europaea var.europaea) cultivation worldwide.Verticillium dahliae Kleb., the causal agent of VWO, is responsible for substantial production losses thereby constituting the main phytopathological issue in many olive-growing regions.Multiple factors explain the current spread and importance of the disease: the long-term endurance of microsclerotia (resistance structures of V. dahliae) in soil, the broad host range of the pathogen, the vascular nature of the disease, inappropriate agronomical practices, etc. [1].
All these factors make it difficult for an efficient control of VWO, and available measures have so far proven unsuccessful when implemented individually.Therefore, an integrated disease management strategy is recommended, highlighting pre-plating control tools since no resistant varieties and/or effective fungicides are available.Within this context, the development and use of biological control agents (BCAs) is an appealing, sustainable and environmentally friendly management approach to control VWO.
Biocontrol experiments are often carried out using sterile soils/substrates; therefore, major determinants influencing results are nutrient availability and/or soil physico-chemical conditions [2].Moreover, these bioassays are frequently performed using plants artificially inoculated via root dipping in (or soil drenching with) the pathogen's conidia suspensions.However, this experimental scenario is quite different to the situation taking place under field conditions.For instance, in sterile soils there is no competence with other microorganisms and the BCA can easily be established.When several microorganisms are used at the same time under these very specific (and artificial) conditions, the final outcome is strongly determined by the inherent capability of each microorganism to multiply and prevail over the others.In contrast, under natural soil conditions, the introduced BCAs face a very complex situation in a harsh environment considered as one of the most microbiologically diverse on earth [3].This is particularly true when the introduced microorganism is not adapted to the target ecological niche.
The use of BCAs under greenhouse or nursery conditions is easier, cheaper, and less time-consuming than large-scale field treatments which, in addition, often suffer from inconsistency or lack of effectiveness.Hence, implementing biological control measures to protect woody plants (e.g., olive) at the propagation stage (i.e., greenhouses or nurseries) against soil-borne pathogens like V. dahliae appears as a preventive and recommendable approach.This 'pretreatment strategy' can be performed before olive plants are transferred to the field, usually cultivated in large growing areas where biocontrol approaches are often unfeasible both from economical and agronomical perspectives [4].
Our current knowledge on fungi, oomycetes and bacteria as effective BCAs against V. dahliae infecting different crops is steadily increasing (reviewed by [5]).Berg et al. [6] already demonstrated that naturally occurring rhizobacteria are effective BCA against V. dahliae.Zhengjun et al. [7] reported on rhizosphere and endophytic bacteria able to control Verticillium wilt of cotton.Berg and Lottmann [8] informed that Verticillium wilt of oilseed rape could be controlled by a Stenotrophomonas maltophilia (Hugh) Palleroni and Bradbury.strain.Tjamos et al. [9] selected two Bacillales strains, Bacillus sp.5-27 and Paenibacillus alvei K-165, from root tips of tomato plants grown in solarized soils able to inhibit V. dahliae's growth and to successfully control Verticillium wilt.The use of strain K-165 also reduced germination of V. dahliae microsclerotia in root tips and in the root elongation zone of eggplants [10].Malandraki et al. [11] selected two strains of the Pseudomonas fluorescens Migula.complex and two isolates of the Fusarium oxysporum Schltdl.complex, all capable to reduce disease severity in eggplants.Likewise, Veloso and Díaz [12] demonstrated that inoculation of pepper plants with the nonpathogenic isolate F. oxysporum 47 reduced Verticillium wilt blight symptoms.Regarding biocontrol of VWO, advances are still limited, although some examples have shown promising results.For instance, olive root-associated fluorescent pseudomonads have been demonstrated to control VWO caused by the highly virulent, defoliating (D) pathotype of V. dahliae in young, nursery-produced plants [13][14][15][16].Paenibacillus alvei K-165 has also been shown to suppress VWO under both greenhouse and field conditions [17].Finally, Acetobacter aceti VIN02 reduced the percentage of mortality in V. dahliae-inoculated plants to the level observed in non-inoculated plants [18].However, even though several microorganisms have been shown to be effective BCAs against Verticillium wilt, hardly any of them are available as biopesticides against Verticillium in Europe [19].
Disease suppressive soils, organic amendments, and healthy plants enduring in pathogen-infested fields are interesting reservoirs of potential BCAs [5,11,16,20].Microbial communities associated with plant roots are complex and diverse, and they strongly influence the health and development of the host [21].Moreover, the rhizosphere microbiome plays a significant role in reprogramming defense responses of plants [22].Therefore, to look for beneficial indigenous microbial inhabitants adapted to the conditions found in the olive root/rhizosphere seems to be a good strategy to identify novel BCAs candidates against VWO.In a previous study, we designed a holistic approach to identify, characterize and evaluate new and effective BCAs originating from the roots of healthy nursery-propagated olive plants [16].This strategy was based on the isolation from the target niche, in vitro antagonism assays against relevant olive pathogens (V.dahliae, Phytophthora cinnamomi Rands., Pseudomonas savastanoi (Janse) Gardan et al.), phenotypic and metabolic characterization in order to identify traits associated with biological control or plant growth promotion (PGP), in planta experiments under non-gnotobiotic conditions to assess their effectiveness against the D pathotype of V. dahliae, genome sequencing and in silico analysis, taxonomic identification by multi-locus sequence analyses (MLSA), and evaluation of root colonization abilities [16].Eventually, we were able to identify and in-depth characterized three novel Gram-negative resident olive rhizobacteria belonging to the Pseudomonas genus (namely, strains PIC25, PIC105 and PICF141) showing good biocontrol ability against VWO.Current research trends in biocontrol aim to investigate the role and effectiveness of microbial consortia, synthetic communities and even tailored microbiomes in promoting plant growth and/or defense against pathogens [23,24].The use of bacterial consortia to control soil-borne pathogens may have the advantage over single microbe-based formulations to enhance biocontrol effectiveness.Indeed, the introduction of diverse beneficial microbes may facilitate the colonization of different microhabitats on/in roots thereby enhancing competition for colonization sites used by the pathogen.Furthermore, a range of biocontrol mechanisms can be deployed by each BCA in the consortium contributing in this manner to improved disease suppressiveness [25].
This study aims to identify and characterize Gram-positive bacteria from the olive root/rhizosphere with potential as BCAs against olive pathogens, with emphasis on V. dahliae.We implemented the same comprehensive strategy that previously yielded successful results when screening and evaluating indigenous beneficial Pseudomonas spp.An additional aim was to in vitro assess whether combinations of these previously verified BCAs and the newly described antagonists here reported, all adapted to the olive rhizosphere, could be formulated as consortia to control VWO.

Olive Roots Sampling and Bacterial Isolation
This study focused on three Gram-positive strains (namely PIC28, PIC73 and PIC167) of an olive rhizobacteria collection obtained from roots of one-year-old 'Picual' plants, a cultivar qualified as highly susceptible to V. dahliae [26], collected from ten commercial nurseries located in Córdoba province (South Spain).Experimental details on olive roots sampling, bacteria isolation and their preliminary identification, and generation of the olive rhizobacteria collection have been previously described [27].Briefly, individual root tissue samples (three plants per nursery) were ground using a mortar and pestle.Aliquots of 10-fold serial dilutions of the macerates obtained from each sample were plated on different culturing media: potato dextrose agar (PDA; Oxoid Ltd., Basingstoke, Hampshire, UK), LuriaBertani agar (LBA; Oxoid Ltd., Basingstoke, Hampshire, UK) and nutrient agar (NA; Oxoid Ltd., Basingstoke, Hampshire, UK).Plates were incubated at 28 • C in the dark for at least 72 h.Single colonies grown in these media showing distinctive morphological characteristics were selected and isolated.Each colony were cryopreserved as described [27] until use.

DNA Isolation, Molecular Identification and Phylogenetic Analysis of Selected Olive Rhizobacteria
Total DNAs from bacteria were extracted using the 'JETFLEX Genomic DNA Purification Kit' (Genomed Ltd., Harrow, Middlesex, UK) according to the manufacturer's instructions.Bacteria were molecularly identified (at least at the genus level) by partial sequencing of the gyrB gene.PCR experiments and purification of DNA fragments were performed as previously described [27].
PCR products were sequenced at Sistemas Genómicos S.L. (Paterna, Valencia, Spain) and assembled in 'contigs' using the 'CLC Main Workbench' (version 7.6.4,CLC Bio, Qiagen, Aarhus, Denmark).The genus/species of the strains under study were identified by comparing the resulting contigs with NCBI Blast using the BLASTN algorithm [28].
An MLSA (5717 nt positions) was conducted using partial sequences of several housekeeping genes: gyrB (1582 nt), atpA (1518 nt), nusA (758 nt), pyrH (732 nt) and dnaJ (1127 nt) in order to refine the taxonomic position of strains PIC28, PIC73 and PIC167 previously obtained by just comparing the gyrB gene sequence.Gene sequences, except that of the gyrB gene, were obtained from the three genomes here sequenced (see below) and then compared with homolog sequences of ad hoc selected Bacillus spp.(13) and Paenibacillus spp.(15) type strains retrieved from different public databases, (i.e., NCBI, EMBL, KEGG, etc.), and that were phylogenetically related to the three olive rhizobacteria under study.A circular dendrogram was then constructed with Molecular Evolutionary Genetics Analysis version 7 (MEGA7) [29], using the Neighbor-Joining algorithm.For this phylogenetic tree Paenibacillus larvae dsm 25,430 was used as out-group species.Bootstrap analyses of 1000 replications were performed.Distance matrixes were generated using the Maximum Composite Likelihood model [30] with MEGA7 [29], and percentages identity of the concatenated sequences were inferred with the Clustal 2.1 version.

In Vitro Antagonism Tests
Strains PIC28, PIC73 and PIC167 were in vitro tested to assess their growth inhibition ability against several olive pathogens (Table 1), following an experimental setup previously described [16,27].In addition, and in order to quantify the growth inhibition degree of V. dahliae, antagonism tests against a representative of the D pathotype (strain Lebrija 1, [16]) were carried out.Individual drops (10 µL) containing V. dahliae biomass (i.e., mycelium plus conidia, 1 × 10 6 conidia/mL) were placed in the middle of PDA plates, and four equidistant 10-µL drops of each bacterial strain (1 × 10 8 ufc/mL) were inoculated.Control plates containing only an individual drop of the pathogen biomass were included in each assay.The relative inhibition index was calculated according to the equation (Rc-Ra)/Rc, where Rc is the average radius of V. dahliae colony in the absence of the antagonist bacterium and Ra is the average radius of V. dahliae colony in the presence of the antagonist bacterium (four equidistant points).This experiment was repeated once including three biological replicates per each pathogen-BCA interaction.Statistical significance was determined by ANOVA analysis of the relative inhibition index.The Tukey HSD Test at p = 0.05 was applied for comparison of mean values by using the analytical software 'Statistix' (Version 10.0 for Windows, 1985-2013, Tallahassee, FL, USA).

Biocontrol Ability of Selected Olive Rhizobacteria Against Verticillium Dahliae D Pathotype
Two independent biocontrol experiments were carried out under non-gnotobiotic conditions to test the ability of the three selected strains, PIC28, PIC73 and PICF167, to control VWO.In these bioassays strains P. fluorescens PICF7 [31] and Pseudomonas sp.PICF141 [16] were included for comparative purposes because of their well-known biocontrol activity.Fifty-four (Experiment I; ExpI) ninety and (Experiment II; ExpII) 'Picual' plants (5-month-old) were grown in pots (11 × 11 × 11 cm; one plant per pot), each containing 300 g of potting substrate used in the nursery, and distributed in three random blocks (nine plants in ExpI and fifteen plants in ExpII per treatment) within a greenhouse under natural lighting and day/night temperature of 27/21 • C. For each bacterial treatment, olive plants were inoculated as described by Gómez-Lama Cabanás and co-workers [16].Control plants were only amended with 150 mL of SO 4 Mg•7H 2 O 10 mM.One week after bacterization, plants were inoculated with isolate V937I of V. dahliae, representative of the D pathotype, as previously described [16].Disease progress was assessed twice a week by scoring the severity of symptoms on a 0-4 scale according to the percentage of affected leaves and twigs (0, no symptoms; 1, 1-33%; 2, 34-66%; 3, 67-100%; 4, dead plant) during 90 days after pathogen inoculation [15].Parameters such as Disease Incidence (DI), Mortality (M), and a Disease Intensity Index (DII) were also calculated for each treatment.This latter parameter was defined as DII = (ΣSi × Ni)/(4 × Nt), where Si is severity of symptoms, Ni is the number of plants with Si symptoms severity, and Nt the total number of plants.A final DI was also established as the percentage of affected plants at the end of the experiments.Finally, the area under the disease progress curve (AUDPC) over time (days; [32]) and the final severity were calculated.Statistical significance was determined by ANOVA.The Fisher's protected LSD Test at p = 0.05 was applied for comparison of mean values by using the analytical software 'Statistix' (Version 10.0 for Windows, 1985-2013, Tallahassee, FL, USA).

Phenotypic and Metabolic Characterization of Olive Rhizobacteria
The ability of strains PIC28, PIC73 and PICF167 to produce HCN, siderophores, the phytohormone indole acetic acid (IAA) and the volatile 2,3-butanediol, and to display enzymatic activities such as protease, catalase, phosphatase, chitinase, phytase, cellulase, xylanase, β-glucosidase, glucanase, pectinase and amylase was evaluated to identify phenotypes traditionally associated with biological control or PGP.For such an aim, the procedures described by Gómez-Lama Cabanás and co-workers ( [16]; and references therein) were implemented.The GEN III MicroPlate™ system (Biolog, Hayward, CA, USA) was used to obtain the phenotypic profiles of the three olive rhizobacteria here under study.These plates contain 71 carbon sources and 23 antimicrobials.This system analyzes the ability of a tested microorganism to metabolize a range of compounds, as well as enables to determine important physiological properties such as pH, salt and lactic acid tolerance, reducing power, and chemical sensitivities.Bacteria suspensions were prepared as recommended by the manufacturer.All these experiments were performed at least twice.

Genome Sequencing and Bioinformatics Analysis
The genomes of strains PIC28, PIC73, and PIC167 were sequenced according to a previously described strategy [16] at Sistemas Genómicos S.L. (Paterna; Valencia, Spain).The sequencing was made with a final size of 300 bp for the paired-end reads.Total initial reads were 3,361,104 paired-ends giving a fold coverage of 193.67 for PIC73, 3,061,700 paired-ends (168.36X) for PIC28, and 2,278,298 paired-ends (130.81X) for PIC167.The quality of the raw data was measured and checked using FASTQC tools [33].Illumina Truseq adaptors were removed with the Fastq mcf [34] (v1.04.803) tool, and then a quality filter was made with Cutadapt [35] (v1.9.1) using a quality window value of 30.With the reads clean, the R1 and R2 paired-end reads were merged using Flash [36] (v1.2.11).With a masking step for the low-quality bases and another for complexity reduction, the reads were ready to be assembled.The assemblers used were mainly Megahit [37] (v1.0.3-29-g707d683) and Velvet [38] (v1.2.10).Several k-mer sizes were used, from 71 to 99.The best assemblies were selected using the best N50 criteria, and over this selected assembly the annotation process was made using Glimmer3 [37][38][39][40] for the ORF detection.Blast V.2.2.30+ [41] with an E-value cutoff of 1e −3 against the latest version (UniProtKB/Swiss-Prot Release 2015_08) of the Uniprot Swissprot protein curated database for bacteria [42] was used for the annotation.The small local alignments were removed from the Blast result file applying a filter.All the sequences without a hit after removing the small local alignments, were annotated using Blast V.2.2.30+ [41] as previously described [16].Identified genes were functionally annotated using functional annotation of Uniprot database according to the three different functional categories; biological process, molecular function and cellular component from Gene Ontology [43].To improve the understanding of the sequenced genomes and their arrangements, each one of them was aligned to the closest relative species available at NCBI [28] using Projector2 [44].The new pseudoassembly sequences were processed again using Glimmer3 [40] and the new ORFs were annotated with Blast V.2.2.30+ [41] against Uniprot (The Universal Protein Resource Database, 2008) [42] for the functional annotation.The associated Kegg Ontology pathways [45], Gene Ontology [43] and PFam [46] terms were obtained to count and sort them in order to create a functional view of the differences.Genome sequences were deposited at Genbank under the accession IDs SAMN08323495, SAMN08323496 and SAMN08323497 for PIC28, PIC73 and PIC167, respectively.

In Silico Identification of Genetic Factors Involved in Plant-Bacteria Interactions
'PIFAR' web-based tool [47] was used to identify genetic factors involved in the interaction of bacteria with host plants.Thus, the presence of genes related to adhesion, antibiotic production, biofilm, detoxification, secretion of exopolysaccharides (EPSs), bacterial lipopolysaccharides (LPSs), microbe-associated molecular patterns (MAMPs), multidrug resistance (MDRs), plant cell wall-degrading enzymes (PCWDEs), phytohormones, phytotoxins, pigmentation, proteases, siderophores, type III effectors and volatiles was investigated in the three olive rhizobacteria genomes under study.Furthermore, by using the 'T346Hunter' web-based tool [48] we also aimed to uncover the presence of three types of bacterial secretion systems (TSS), namely flagellar and non-flagellar type III secretion systems (T3SS), type IV (T4SS) and type VI (T6SS).Finally, since an indispensable prerequisite for the use of BCAs in agriculture is the absence of undesirable traits (i.e., genetic factors potentially deleterious for plants, animals and humans), the annotated genomes of the three selected rhizobacteria were manually examined in order to discard the presence of genes potentially involved in pathogenesis/virulence processes.

In Vitro Antagonism among Olive Rhizobacteria Displaying Biocontrol Against VWO
In order to evaluate whether well-characterized Pseudomonas spp.strains indigenous to olive roots and with demonstrated biocontrol activity against VWO (i.e., strains PIC25, PIC105, PICF141 [16] and PICF7 [31]) could be combined and used in future formulations (consortia of BCAs) with the Bacillales members described in this study, in vitro antagonism tests were carried out.Four drops (5 µL, 10 7 cfu mL −1 ) of each strain (challenger) were placed on the surface of LBA (Luria Bertani Agar) and PDA media plates previously inoculated with 100 µL of a bacterial suspension of each strain (10 6 cfu mL −1 ) (target) in all possible combinations.Two series of plates (two per assayed media) were incubated at 28 • C in the dark, and after 48 h, bacterial growth and presence/absence of inhibition halos around colonies were observed.For each combination, the assay was performed twice.

Molecular Identification and Phylogenetic Analysis of Newly Isolated Olive Rhizobacteria
The gyrB sequences (amplicon size ≈ 1200 bp) of the three olive rhizobacteria were compared to sequences deposited on GenBank using a BLASTn search.The analysis enabled us to identify the three strains as members of the order Bacillales.Strain PIC28 showed high identity (99%, e-value 0.0) with two different species from the Bacillus genus, Bacillus cereus Frankland and Frankland.(strain C1L) and Bacillus thuringiensis serovar canadensis (strain IEBC-T05A001).Therefore, the gyrB gene sequence did not allow the accurate identification (species level) of this strain.In contrast, strain PIC73 could be unambiguously identified as Paenibacillus polymyxa (Prazmowski) Ash et al. (99% identity, e-value 0.0), and strain PIC167 as Paenibacillus terrae Yoon et al. (96% identity, e-value 0.0).
In order to more accurately identify the three novel Bacillales strains, an MLSA (5717 nucleotide positions) of concatenated partial sequences of the housekeeping genes gyrB, atpA, nusA, pyrH and dnaJ, that were obtained from their sequenced genomes (see below) was performed.A phylogenetic tree was then generated including our three olive rhizobacteria and 13 Bacillus spp.type strains belonging to five different species and 15 Paenibacillus spp.type strains belonging to 12 different species (Figure 1).Results enabled to assign strains PIC73 and PIC167 to P. polymyxa and P. terrae species, respectively, in agreement with the preliminary molecular identification only based on gyrB gene sequences comparison.However, strain PIC28 clustered ambiguously among different Bacillus species (Figure 1), B. cereus (89% identity) and B. thuringiensis Berliner.(89.59% identity) being the closest relatives.Therefore, this strain was kept as incertae sedis within the Bacillus genus until further evidence.

Bacillales from Olive Roots Showed In Vitro Antagonism Against Different Olive Pathogens, Including Verticillium dahliae
Results from in vitro antagonism tests showed that the three Bacillales strains originating from olive roots were able to inhibit the growth of different olive pathogens at variable levels (Table 1).The two Paenibacillus spp.strains (PIC73 and PIC167) effectively antagonized all pathogens tested (i.e., V. dahliae, Rosellinia necatrix Prill., P. cinnamomi, Colletotrichum nymphaeae Aa., Colletotrichum godetiae Neerg.and P. savastanoi) at least in one out the two culturing media assayed.However, Bacillus sp.PIC28 was effective neither against R. necatrix nor one P. savastanoi strain (Table 1).Inhibition of pathogen growth depended on the medium used (PDA or NA).When in vitro antagonism tests against V. dahliae were performed to calculate a relative inhibition index, the three olive rhizobacteria were able to significantly (p < 0.05) inhibit the growth of V. dahliae D pathotype, with strain PIC73 showing the highest antagonistic capacity.Relative inhibition indices differed among strains (PIC73 = 0.49 a, PIC167 = 0.47 ab, and PIC28 = 0.45 b).According to these results, the three olive rhizobacteria would be promising BCA candidates against VWO.

Bacillales from Olive Roots Showed In Vitro Antagonism Against Different Olive Pathogens, Including Verticillium dahliae
Results from in vitro antagonism tests showed that the three Bacillales strains originating from olive roots were able to inhibit the growth of different olive pathogens at variable levels (Table 1).The two Paenibacillus spp.strains (PIC73 and PIC167) effectively antagonized all pathogens tested (i.e., V. dahliae, Rosellinia necatrix Prill., P. cinnamomi, Colletotrichum nymphaeae Aa., Colletotrichum godetiae Neerg.and P. savastanoi) at least in one out the two culturing media assayed.However, Bacillus sp.PIC28 was effective neither against R. necatrix nor one P. savastanoi strain (Table 1).Inhibition of pathogen growth depended on the medium used (PDA or NA).When in vitro antagonism tests against V. dahliae were performed to calculate a relative inhibition index, the three olive rhizobacteria were able to significantly (p < 0.05) inhibit the growth of V. dahliae D pathotype, with strain PIC73 showing the highest antagonistic capacity.Relative inhibition indices differed among strains (PIC73 = 0.49 a, PIC167 = 0.47 ab, and PIC28 = 0.45 b).According to these results, the three olive rhizobacteria would be promising BCA candidates against VWO.

Strains PIC28, PIC73 and PIC167 Control VWO Under Greenhouse Conditions
Overall, results from independent biocontrol experiments ExpI and ExpII showed that 'Picual' plants pretreated with each of the three selected olive rhizobacteria reduced VWO development, as evidenced by a significant (p < 0.05) reduction of the AUDPC compared to non-bacterized plants, as well as other phytopathological parameters scored (final DI and DII, mortality and severity) (Table 2).In fact, the three Bacillales strains showed a good biocontrol ability against the D pathotype of V. dahliae as that displayed by P. fluorescens PICF7 and Pseudomonas sp.PICF141, two well-studied BCAs of VWO here used for comparative purposes.Non-inoculated plants (control treatment) and plants inoculated only with the bacterial strains showed normal growth in the two bioassays carried out.Neither abnormal appearance nor unexpected symptoms were observed during the bioassays.While all bacterial/V.dahliae treatments reduced the AUDPC compared with control (non-bacterized) plants, no significant difference in biocontrol performance was observed among olive rhizobacteria in ExpI (Table 2).Likewise, parameters such as final DI and DII, mortality, and severity were also reduced in all BCA treatments.However, in ExpII disease pressure was higher and Pseudomonas strains overall performed better than the Bacillales members.At the end of this bioassay, all BCA/Vd treatments showed significantly (p < 0.05) lower AUDPC values compared with the control (non-bacterized) treatment, although strains PIC73 and PIC167 performed slightly better than PIC28 for all disease parameters scored (Table 2).The three newly identified indigenous olive Gram-positive rhizobacteria were thus able to significantly reduce VWO progress under our experimental conditions, particularly when disease pressure was moderate. 1AUDPC, area under the disease progress curve (90 days).Final DI, final disease incidence (%).Final DII, disease intensity index (ranging 0-1) calculated with data on incidence and severity of symptoms (90 days).M, dead plants at the end of the experiment (%) (90 days).S, mean of disease severity symptoms (from 0 to 4) at the end of the experiments (90 days).Data are the average of three randomly distributed blocks each with three and five pots per treatment in Experiments I and II, respectively.Non-inoculated plants did not show any disease symptom and were excluded for statistical analysis.In each column, values followed by different letters are significantly different according to contrast Fisher's protected LSD test (p = 0.05).

Identification of Phenotypes Associated with Biological Control or Plant Growth Promotion
All strains were positive for protease, cellulase, β-glucosidase and catalase activities, as well as for siderophore(s) and 2,3-butanediol production.Strains PIC73 and PIC167 showed a higher number of positive activities than strain PIC28, including production of chitinase, amylase and xylanase.None of them were positive for the production of HCN or for pectinase and IAA activities (Table 3).Very weak phosphatase and glucanase activities were detected only for P. polymyxa PIC73.Strain PIC167 was the only strain showing phytase activity (Table 3).Results from the Biolog GEN III assays showed that Paenibacillus spp.strains (PIC73 and PIC167) shared most of the properties tested (85/94) (i.e., biochemical properties, utilization of carbon sources, chemical sensibilities) and differed only in a few metabolic characteristics (Supplementary Table S1).For instance, strain PIC167 was able to utilize L-aspartic acid, methyl piruvate, L-malic acid, bromo succinic acid and acetic acid, whereas strain PIC73 was unable to utilize any of these substrates but could grow, in contrast to PIC167, at pH 5 and in the presence of 4% NaCl, rifamycin, lincomycin, guanidine HCl and niaproof 4. The Biolog profile exhibited by strain PIC28 was clearly different from those of the Paenibacillus strains.For instance, only nine carbon sources were used by all three rhizobacteria.Moreover, PIC28 was the only strain showing tolerance to D-serine (Supplementary Table S1).

Whole-Genome General Features of the Newly Identified Biocontrol Agents
A general summary of the three genome sequencing projects is shown in Supplementary Table S2.The draft genomes of strains PIC28, PIC73 and PIC167 included 43, 31 and 18 large mapped contigs (largest contigs size were 1,384,757, 863,368 and 1,232,833 nucleotides, respectively) and 75, 22 and 7 non-mapped contigs, respectively.Additional relevant characteristics (genome size, GC content, predicted number of coding sequences, number of rRNA operons, etc.) are shown in Table 4.All these characteristics are within the range of previously sequenced genomes of Bacillus and Paenibacillus [49][50][51].A total of 5905, 5716 and 5413 predicted genes were identified for PIC28, PIC73 and PIC167, respectively.Predicted genes annotated for each strain are listed in Supplementary Table S3.Approximately 4012, 3757 and 3598 of them were predicted as protein-coding genes.Functional categories for each coding DNA sequences according to the COG (Clusters of Orthologous Groups) are summarized in Supplementary Table S4.

In Silico Identification of Bacterial Secretion Systems and Genetic Factors Involved in Plant-Bacteria Interactions
According to the 'T346Hunter' web-based tool, only flagellar T3SS was identified although differences were found among the three olive rhizobacteria.Strains PIC73 and PIC167 showed identical TSS profiles, and the flagellar T3SS1 and T3SS2 (PIC73 containing 33 related genes and PIC167 34 genes, including fla, flb, fle, fli, flg and flh genes) were detected in their genomes.However, only the flagellar T3SS1 was found in the genome of strain PIC28 (30 related genes including mot, fle, fli, flg, flh and flg).In addition to the detected T3SSs, genes putative coding for the T7SS or genes related to this TSS (not detectable by the 'T346hunter' tool) were found after genome annotation (WXG100 family type VII secretion target and different type VII secretion protein such as EsaA, EssB and EssC) (Supplementary Table S5).
The 'PIFAR' open-access web-based tool allowed the identification of several genetic factors involved in plant-bacteria interactions.A summary of the number of genes and percentage of genetic factors identified in the three newly identified Bacillales strains is displayed in Figure 1.Moreover, detailed information about the different genetic factors involved in plant-bacteria interactions is shown in Table 5.The number of annotated genes involved in plant-bacteria interaction present in the Paenibacillus strains (PIC73 and PIC167) was very similar, 86 and 82, respectively (Figure 2).The same gene clusters were identified in the genomes of the two strains except for genes related to metabolism of certain compounds enabling a BCA to proliferate inside the plants, EPSs and adhesion factors, confirming the close relatedness of these two olive rhizobacteria.In contrast, the profile of genetic factors related to plant-bacteria interaction differed significantly in strain PIC28 from that in strains PIC73 and PIC167, the former harboring a lower number of this type of factors ( 52) than the latter.For instance, a considerably small number of PCWDEs and antibiotic-related genes were found in strain PIC28.In contrast, more siderophore-(six genes) and biofilm-related (five genes) factors were identified in the PIC28 compared to the Paenibacillus strains (one siderophore and three biofilm genes) (Figure 2).The in silico analysis also confirmed that only the genome of Bacillus sp.PIC28 harbored genes encoding for the phytohormone IAA, pigments (Rubrifacine) and type III effectors (RipAP).Besides, the genome of strain PIC73 was the only one in which a gene related to adhesion factors was identified.
Agriculture 2018, 8, x FOR PEER REVIEW 12 of 24 detailed information about the different genetic factors involved in plant-bacteria interactions is shown in Table 5.The number of annotated genes involved in plant-bacteria interaction present in the Paenibacillus strains (PIC73 and PIC167) was very similar, 86 and 82, respectively (Figure 2).The same gene clusters were identified in the genomes of the two strains except for genes related to metabolism of certain compounds enabling a BCA to proliferate inside the plants, EPSs and adhesion factors, confirming the close relatedness of these two olive rhizobacteria.In contrast, the profile of genetic factors related to plant-bacteria interaction differed significantly in strain PIC28 from that in strains PIC73 and PIC167, the former harboring a lower number of this type of factors ( 52) than the latter.For instance, a considerably small number of PCWDEs and antibiotic-related genes were found in strain PIC28.In contrast, more siderophore-(six genes) and biofilm-related (five genes) factors were identified in the PIC28 compared to the Paenibacillus strains (one siderophore and three biofilm genes) (Figure 2).The in silico analysis also confirmed that only the genome of Bacillus sp.PIC28 harbored genes encoding for the phytohormone IAA, pigments (Rubrifacine) and type III effectors (RipAP).Besides, the genome of strain PIC73 was the only one in which a gene related to adhesion factors was identified.Table 5. Genetic factors involved in plant-bacteria interaction identified in the genomes of the Bacillales strains PIC28, PIC73 and PIC167, according to the web-based tool 'PIFAR' [47].Table 5. Genetic factors involved in plant-bacteria interaction identified in the genomes of the Bacillales strains PIC28, PIC73 and PIC167, according to the web-based tool 'PIFAR' [47].Finally, after manually examining the annotation of the three genomes here under study, at least 36 (PIC28), 22 (PIC73) and 24 (PIC167) genes potentially related to pathogenesis/virulence factors were detected (Supplementary Table S5).Among them, genes related to TSS7 (see above) were identified in the three strains, and three genes putatively coding for an enterotoxin were detected in the genome of strain PIC28.

In Vitro Antagonism among Effective BCAs of VWO
In order to assess whether the Bacillales strains here characterized could be combined (i.e., formulated as consortia) with other BCAs also effective against VWO (e.g., Pseudomonas spp.strains previously characterized by our group), in vitro dual tests with all possible strain combinations were performed (Table 6).Results varied depending on the culturing media used (PDA or LBA).In vitro antagonism assays for strains PIC25 and PIC105 were carried out only in LBA medium since these strains did not grow well in PDA medium.Only few antagonism events were observed and are summarized in Table 6.According to these results, we conclude that all Pseudomonas strains could be combined in a consortium since no antagonism was observed among them.Moreover, the two most effective BCAs from the Pseudomonas spp. group (i.e., PICF7 and PICF141), as well as P. indica PIC105, could also be combined with P. polymyxa PIC73.
PDA, Potato Dextrose Agar; LBA, Luria Bertani Agar.ND, Not determined since strains PIC25 and PIC105 did not grow well in PDA medium; +, positive antagonism; −, no antagonism.For experimental details see Material and Methods.

Discussion
The development and use of BCAs against V. dahliae is an attractive measure within an integrated management strategy of VWO.The early introduction of the BCA, for instance at the nursery propagation stage, could 'condition' the plant to better confront subsequent infections by the pathogen's propagules existing in the planting site.Moreover, selection of BCAs sharing the same ecological niche as V. dahliae is a plausible approach since these microbial antagonists can compete with the pathogen for infection sites, space and nutrients [5].For instance, in the tripartite interaction V. dahliae-olive-P.fluorescens PICF7, niche overlap between the BCA and the pathogen in planta was necessary for effective biocontrol [52].Among the bacteria thriving in the soil/rhizosphere, some spore-forming plant growth-promoting rhizobacteria (PGPR) and/or BCAs have gained attention due to their advantages over non-spore-formers in bioproduct formulation and stable maintenance in soil.Moreover, they are usually persistent to adverse conditions and can be preserved and distributed easily [53].
The holistic strategy previously designed to identify Pseudomonas spp.strains from the olive root/rhizosphere as BCAs against VWO [16], has also been successfully implemented to in-depth characterize Gram-positive bacteria effective against the D, highly virulent pathotype of V. dahliae.Eventually, three Bacillales members were selected (PIC28, PIC73 and PIC167) according to their ability to in vitro antagonize V. dahliae isolates and other olive pathogens.The MLSA performed enabled to identify strain PIC73 as P. polymyxa and strain PIC167 as P. terrae.However, strain PIC28 could only be assigned to the genus Bacillus and further analysis will be needed to determine the species level.Overall, results obtained from the phenotypic and metabolic analyses (e.g., Biolog Gen III) performed in our study confirmed the similarities between the two Paenibacillus strains.
The use of spores of Bacillus spp. as natural 'biopesticides' has previously been described as a safe and environmentally friendly alternative to traditional control methods.Bacillus cereus and B. subtilis (Ehrenberg) Cohn.[54] strains showed fungicidal and PGP ability when they were inoculated directly into soils as seed pretreatments.Moreover, the high survival rate of their spores in soils play an important role in its effectiveness ( [55], and references therein).Likewise, one advantage of treating fields with endospore-forming bacteria such as Paenibacillus spp. is their capability to persist for long periods of time in the soil under harsh environmental conditions [56].Representatives of the genus Bacillus has been demonstrated as BCAs of Verticillium wilts, highlighting Bacillus amyloliquefaciens Priest et al. and B. subtilis strains.For instance, B. amyloliquefaciens 5-127 reduced the percentage of diseased leaves in eggplants challenged with Verticillium dahliae in the greenhouse, and also decreased the DI in potato under field conditions [9].Bacillus strains have been also demonstrated to be effective against Verticillium wilt in oilseed rape, cotton and strawberry [5].Only a few studies have been focused on biological control of Verticillium wilt by Bacillales members in woody plants.Thus, some B. subtilis strains isolated from healthy maple trees were checked in a greenhouse against V. dahliae in this tree species, and disease reduction was reported [57].The genus Paenibacillus is currently gaining interest as BCA of a range of plant diseases [58,59].For instance, P. alvei K-165 obtained from tomato plants cultivated in solarized soil [9] showed effective biocontrol of V. dahliae in eggplant under greenhouse conditions [4,9,10,60].This strain was also able to control V. dahliae in potato plants under field conditions and in olive trees under both field and greenhouse conditions [9,17].In cotton, treatment with Paenibacillus xylanilyticus YUPP-1 and P. polymyxa YUPP-8 reduced DI and disease severity [61].
Results from in vitro antagonism assays showed differences depending on the culture medium used (PDA or NA).The influence of the culture media on the results of in vitro antagonism tests has been previously reported [16,62], a likely consequence of the effects on antagonistic bacteria of specific nutritional factors present in the used media.According to our in vitro antagonism tests, the most promising olive rhizobacteria to be used as BCAs against VWO would be the Paenibacillus spp.strains, particularly strain PIC73.Yet, these results are not enough to qualify them as effective BCAs since colonization of the target niche and biocontrol performance by a selected antagonist will be influenced by environmental factors (moisture, nutrients, pH, temperature, etc.) and its ability to compete with preexisting microbial communities present in the site where it will be deployed.It is well known that the capacity of a beneficial bacterium to act as a BCA could vary between laboratory and field conditions, and that the best in vitro antagonist is not always the best BCA [63][64][65][66].Therefore, in planta bioassays are important to qualify an in vitro antagonist as a true and effective BCA.Our independent VWO biocontrol experiments (performed under non-gnotobiotic conditions) showed that the three Bacillales members indigenous from olive roots were able to significantly reduce VWO development, and that their disease control ability was comparable to that observed for Pseudomonas spp.strains well characterized as VWO BCAs [15,16].Furthermore, all of the tested BCAs are well adapted to the olive root/rhizosphere since they originated from this niche.However, that biocontrol performance slightly decreased when the disease pressure was high, particularly in the case of Bacillus sp.PIC28.Therefore, it is highly recommended that different experiments are conducted to obtain an accurate picture of the biocontrol performance of any BCA candidate.
Biocontrol activity of PGPR relies, among other mechanisms, on the production of different cell wall-degrading enzymes (e.g., chitinase, protease/elastase, and β-1,3 glucanase, etc.) that play an important role in degrading the cell walls of phytopathogens [67,68].Some of these enzymes also degrade plant cell walls and some of their components can be utilized by BCAs, thereby enhancing their competitiveness during the root colonization process.The three Bacillales strains here studied were positive for protease, cellulase and β-glucosidase.Proteases play a significant role in cell wall lysis of phytopathogenic fungi, since proteolytic activity is a prerequisite to lyse whole fungal cells [69].Microorganisms displaying cellulase and protease activities do not only play a role in organic matter decomposition and PGP, but also contribute to disease suppression by inhibiting fungal growth [70].Only Paenibacillus spp.strains PIC73 and PIC167 strains showed chitinase, amylase and xylanase activities.Chitinases lyse fungal cell wall through degradation of the chitin polymer present in fungal cell walls [71].Strain PIC73 was the only one to show phosphatase and glucanase activities, although at very low levels.Phosphatases mediate the release of inorganic P from organically bound P returned to soil [72], increasing the root system and thus contributing to the growth of the plant.Glucanases participate in the degradation of cell walls in fungi, yeasts and plants [73].Complementing this phenotypic characterization, in silico analysis revealed that the genome of Bacillus sp.PIC28 harbored a considerably lower number of genes potentially coding for cell wall-degrading enzymes (as well as antibiotics) than the two Paenibacillus strains.The ability to produce fungal cell wall-degrading enzymes and a broad spectrum of antibiotics seem to be common to aerobic spore-forming bacteria [74].
Paenibacillus terrae PIC167 was the only strain displaying phytase activity.Phytases catalyze the hydrolysis of phytic acid, releasing a utilizable form of inorganic phosphorus (P) to the plants [75].Rhizosphere bacteria with phytase activity have been isolated and proposed to promote plant growth in soils with high content of organic P [76].Another trait traditionally associated with PGPR is the production of the phytohormone IAA.In silico analysis revealed the presence of two genes related to IAA biosynthesis in the genome of Bacillus sp.PIC28.However, production of IAA by this strain could not be demonstrated.Some PGPR elicit induced systemic resistance by emission of volatile organic compounds including short chain alcohols, acetoin, and 2,3-butanediol [77].In silico analysis predicted the ability to produce several volatile compounds, including 2,3-butanediol, for the three Bacillales strains.Indeed, production of this volatile was confirmed for all of them; 2,3-butanediol produced by Bacillus spp. is suggested to be involved in plant protection [78].Different types of siderophores promote the growth of several plant species and increase their yield by enhancing iron uptake in plants [79].Additionally, siderophore-producing microorganisms can suppress some soil-borne fungal pathogens through a siderophore-mediated iron competition mechanism.Genes putatively involved in siderophore(s) biosynthesis were detected by in silico analysis of the genomes of the three Bacillales strains.Moreover, siderophore activity was confirmed by the CAS (Chrome Azurol S) for all of them (Table 3), Bacillus sp.PIC28 always showing much bigger orange halos than those observed for the Paenibacillus strains.
Biofilm formation by beneficial strains of Bacillus spp.and Paenibacillus spp. is known to improve their biocontrol abilities [80][81][82].Biofilms facilitate the development of a microhabitat that protects bacteria against a/biotic stresses [83,84], enhancing root competitiveness of bacteria as well.The presence of genetic factors related to biofilm development has been predicted for strains PIC28, PIC73 and PIC167, although the number of genes related to this trait was higher in strain PIC28.The presence of TSSs in the newly identified Bacillales members of the olive roots was also investigated by in silico analysis.Secretion systems in biocontrol bacteria have been related to plant root colonization, rhizosphere competence, environmental competition or defense against pathogens [85][86][87].According to 'T346hunter' predictions, only the flagellar-T3SS was detected in the three genomes analyzed.It is worth mentioning that 'T346hunter' and 'PIFAR' web-based tools only detect positive matches for (nearly) complete gene clusters, and that factors defined by several genes are reported only if at least 90% of such genes are identified within a given genome [47,48].Many Gram-positive bacteria also produce a specialized protein secretion machinery termed the type VII secretion system (T7SS) that is not detected by the 'T346Hunter' tool.Genes related to T7SS have been detected during the examination of the annotated genomes.T7SS are widespread in representatives of the phlya Actinobacteria and Firmicutes, and affect a range of bacterial processes including virulence, sporulation, conjugation, and cell wall stability [88].Bioinformatics analyses have also predicted that this TSS is present in non-pathogenic bacteria [89].WXG-100-encoding genes have also been identified in different species of the genus Bacillus, both in the non-pathogenic species B. subtilis and in pathogenic representatives of B. cereus, B. thuringiensis and Bacillus anthracis Cohn.[90,91].Although some T7SSs can contribute to virulence, not all of these systems function in virulence or even in secretion of substrates.However, the biological functions of several type VII-secreted substrates and effector molecules are still being investigated.Finally, when the annotated genomes here reported were manually and thoroughly checked, genes putatively coding for pathogenicity/virulence factors were found, Bacillus sp.PIC28 being the strain harboring the largest number of genes potentially related with deleterious, non-desirable effects including three genes involved in enterotoxin synthesis.This fact must be carefully examined regarding the use of this strain in biocontrol frameworks.
Recent trends in biocontrol of plant diseases have focused on the development and use of consortia of beneficial microorganisms [25,92].This approach may improve the efficacy, reliability and consistency of individual microbe-based biocontrol tools [93].For instance, the mixture of two Bacillus strains increased protection in tomato and pepper against different pathogens [94].Also, the application of a consortium of three rhizobacteria, B. cereus AR156, B. subtilis SM21 and Serratia sp.XY21, resulted in higher biocontrol efficacy against Verticillium wilt in cotton compared to the outcomes observed with individual treatments [95].To be formulated in a consortium, microorganisms in the mixture must not show antagonism among them [96][97][98].Olive roots are a powerful reservoir of BCAs against VWO, mainly Pseudomonas spp.[13,16] and Bacillales (this study) strains.In order to check whether these BCAs could be formulated as consortia, the absence of antagonism among them was assessed.Results showed that, for instance, a consortium consisting of Pseudomonas spp.strains was feasible since no incompatibilities were detected among them.Moreover, a consortium based on the two most effective Pseudomonas strains (PICF7 and PICF141) [15,16] and P. polymyxa strain PIC73 would also be interesting to be tested in future biocontrol assays.These combinations would provide an array of BCAs colonizing different microhabitats or sites of the olive root/rhizosphere, thereby enhancing competition for the colonization sites also used by V. dahliae.Moreover, different biocontrol mechanisms (antagonism, competition, induction of host defense responses) [52,53,99,100] could be deployed by the olive indigenous BCAs in the consortia, improving biocontrol effectiveness.

Conclusions
Three novel indigenous Bacillales strains from the olive rhizosphere and effective against the (D) pathotype of V. dahliae under non-gnotobitic conditions were isolated and characterized.Among them, strain PIC73 (P.polymyxa) was the most promising BCA according to all parameters here analyzed.Biocontrol measures based on the introduction of a single BCA, combinations of them (consortia, microbiomes), or even formulated with organic amendments, are attractive and complementary actions to be implemented within an integrated disease management strategy of VWO [1] (Figure 3).Information here generated will be useful for the future development of effective formulations against this pathogen, including consortia between the newly identified Bacillales members and Pseudomonas spp.strains previously characterized.Most of the studies on VWO biocontrol are focused on preventive (before-planting) measures, mainly working with nursery-propagated olive young plants.In our opinion, it is interesting to introduce of microbial antagonists at the production stage (nurseries) due to easier application routines and to the fact that the plant could be transferred to fields potentially infested with the pathogen with its defense mechanisms 'activated' and/or harboring a cohort of effective antagonists.Our increasing knowledge on natural inhabitants of olive roots/rhizosphere displaying biocontrol activity, that in addition could effectively be combined, will be instrumental to also develop palliative (post-planting) measures to be implemented in a more effective and consistent way under field conditions.
Agriculture 2018, 8, x FOR PEER REVIEW 18 of 24 VWO [1] (Figure 3).Information here generated will be useful for the future development of effective formulations against this pathogen, including consortia between the newly identified Bacillales members and Pseudomonas spp.strains previously characterized.Most of the studies on VWO biocontrol are focused on preventive (before-planting) measures, mainly working with nurserypropagated olive young plants.In our opinion, it is interesting to introduce of microbial antagonists at the production stage (nurseries) due to easier application routines and to the fact that the plant could be transferred to fields potentially infested with the pathogen with its defense mechanisms 'activated' and/or harboring a cohort of effective antagonists.Our increasing knowledge on natural inhabitants of olive roots/rhizosphere displaying biocontrol activity, that in addition could effectively be combined, will be instrumental to also develop palliative (post-planting) measures to be implemented in a more effective and consistent way under field conditions.

Supplementary Materials:
The following tables are available online at www.mdpi.com/link.Table S1: Biochemical analysis of strains PIC28, PIC73 and PIC167 using Biolog GEN III Microplate including utilization of carbon sources and chemical sensibilities; Table S2: General information of the three genome sequence projects; Table S3: Predicted genes annotated in the genomes of strains PIC28, PIC73 and PIC167; Table S4: Number of genes associated with GO functional categories; and Table S5: Potential pathogenesis-related genes identified in annotated genomes of strains PIC28, PIC73 and PIC167.

Figure 1 .
Figure 1.Circular phylogenetic tree showing the taxonomic position of the three novel Bacillales strains (underlined red color) isolated from olive 'Picual' roots.The tree was constructed by the Neighbor-Joining method, based on the alignment of concatenated partial sequences of gyrB, atpA, nusA, pyrH and dnaJ (see text for details).Bar indicates sequence divergence.Bootstrap values (≥95%) based on 1000 re-sampled datasets are shown at branch nodes.Paenibacillus larvae DMS 25430 was used as out-group.

Figure 1 .
Figure 1.Circular phylogenetic tree showing the taxonomic position of the three novel Bacillales strains (underlined red color) isolated from olive 'Picual' roots.The tree was constructed by the Neighbor-Joining method, based on the alignment of concatenated partial sequences of gyrB, atpA, nusA, pyrH and dnaJ (see text for details).Bar indicates sequence divergence.Bootstrap values (≥95%) based on 1000 re-sampled datasets are shown at branch nodes.Paenibacillus larvae DMS 25430 was used as out-group.

Figure 2 .
Figure 2. Number of genes and percentage of genetic factors involved in plant-bacteria interaction detected in the genomes of the olive rhizobacteria PIC28, PIC73 and PIC167.

Figure 2 .
Figure 2. Number of genes and percentage of genetic factors involved in plant-bacteria interaction detected in the genomes of the olive rhizobacteria PIC28, PIC73 and PIC167.

Figure 3 .
Figure 3. Biological control measures within and integrated management strategy of Verticillium wilt of olive (based on López-Escudero and Mercado-Blanco [1]).

Table 2 .
Biocontrol performance of Bacillales members against the defoliating pathotype of Verticillium dahliae.

Table 3 .
Screening of activities associated with plant growth promotion or biocontrol in Bacillales strains.−, absence of activity; ±, ambiguous result.Results shown were obtained from Luria Bertani Agar (LBA) cultures amended with L-glycine.At least two biological replicates for each activity/strain were performed.

Table 4 .
General information of Bacillales genomes.
+, presence of the genetic factor; −, absence of the genetic factor.

Table 6 .
In vitro dual tests among olive rhizobacteria showing biocontrol activity against Verticillium wilt of olive.