Skip to main content

Unearthing the genomes of plant-beneficial Pseudomonas model strains WCS358, WCS374 and WCS417



Plant growth-promoting rhizobacteria (PGPR) can protect plants against pathogenic microbes through a diversity of mechanisms including competition for nutrients, production of antibiotics, and stimulation of the host immune system, a phenomenon called induced systemic resistance (ISR). In the past 30 years, the Pseudomonas spp. PGPR strains WCS358, WCS374 and WCS417 of the Willie Commelin Scholten (WCS) collection have been studied in detail in pioneering papers on the molecular basis of PGPR-mediated ISR and mechanisms of biological control of soil-borne pathogens via siderophore-mediated competition for iron.


The genomes of the model WCS PGPR strains were sequenced and analyzed to unearth genetic cues related to biological questions that surfaced during the past 30 years of functional studies on these plant-beneficial microbes. Whole genome comparisons revealed important novel insights into iron acquisition strategies with consequences for both bacterial ecology and plant protection, specifics of bacterial determinants involved in plant-PGPR recognition, and diversity of protein secretion systems involved in microbe-microbe and microbe-plant communication. Furthermore, multi-locus sequence alignment and whole genome comparison revealed the taxonomic position of the WCS model strains within the Pseudomonas genus. Despite the enormous diversity of Pseudomonas spp. in soils, several plant-associated Pseudomonas spp. strains that have been isolated from different hosts at different geographic regions appear to be nearly isogenic to WCS358, WCS374, or WCS417. Interestingly, all these WCS look-a-likes have been selected because of their plant protective or plant growth-promoting properties.


The genome sequences of the model WCS strains revealed that they can be considered representatives of universally-present plant-beneficial Pseudomonas spp. With their well-characterized functions in the promotion of plant growth and health, the fully sequenced genomes of the WCS strains provide a genetic framework that allows for detailed analysis of the biological mechanisms of the plant-beneficial traits of these PGPR. Considering the increasing focus on the role of the root microbiome in plant health, functional genomics of the WCS strains will enhance our understanding of the diversity of functions of the root microbiome.


Plants have evolved elaborate mechanisms to defend themselves against the plethora of pathogens by which they are attacked [1]. Although the interactions between plants and their attackers are often regarded as bilateral, a mesmerizing diversity of microorganisms in and around the plant influences the outcome of the battles between plants and their enemies [2, 3]. Whereas some microbes team up with ‘the bad guys’ [4], many others help to protect the plant and improve plant survival in ecological and agricultural settings [5]. Studies on structure and function of microbiota in and on plant roots have identified a wide array of bacteria and fungi from different taxonomic groups [69], many of which possess plant growth-promoting activities [10, 11]. Among the group of plant growth-promoting rhizobacteria (PGPR), the genus Pseudomonas is strongly represented. This genus comprises over one hundred species of aerobic bacteria that belong to the γ subclass of the Proteobacteria [12]. Although some Pseudomonas spp. are plant pathogenic, many have been found to protect plants by antagonizing soil-borne pathogens through competition for nutrients, production of antimicrobial compounds, or by eliciting a systemic immune response that is effective against a broad spectrum of pathogens, called induced systemic resistance (ISR) [11, 13]. Mutualistic root-colonizing Pseudomonas spp. emerged as important players in disease-suppressive soils [7, 1416], and served as model PGPR in research toward understanding how non-symbiotic root-associated bacteria protect plants against pests and diseases [11, 17].

In the past 30 years, over 300 publications described biological mechanisms involved in the ability of the PGPR strains Pseudomonas putida WCS358, Pseudomonas fluorescens WCS374, and P. fluorescens WCS417 to promote plant growth and health. WCS358, WCS374, WCS417 were isolated in the 1980’s at the Dutch Phytopathological Laboratory “Willie Commelin Scholten” (WCS; [18, 19]) from the rhizosphere of potato (WCS358 and WCS374; [20]) and wheat (WCS417; [21]). They efficiently colonize the rhizosphere of their host plants and have been shown to reduce plant diseases in many plant-pathogen combinations. Originally, the biological control activity of these WCS strains was linked to siderophore mediated competition for iron e.g. [2228]. Siderophores are small iron-chelating compounds that are secreted by microorganisms under low-iron conditions to enable sequestering and uptake of essential ferric iron from their environment. In most soils, ferric iron is scarcely available as it is present as poorly soluble ferric hydroxides [29]. Microorganisms possess siderophore receptors that specifically recognize and take up their cognate siderophore-iron complexes [30, 31]. In addition to receptors for their own siderophores, many rhizobacteria possess receptors for siderophores produced by other microbes, thereby enhancing their competitiveness under conditions of low iron availability [32, 33].

In addition to their siderophore-mediated antagonistic effect on soil-borne pathogens, the WCS strains also emerged as potent inducers of ISR. Strain WCS417 was among the Pseudomonas spp. strains with which the phenomenon of ISR was for the first time experimentally demonstrated [3436]. With WCS417, ISR was reported for the first time in carnation against Fusarium wilt disease, but since then all three WCS strains have been demonstrated to prime host immunity in different plant species providing enhanced protection against a broad spectrum of plant pathogens and even insect herbivores [13]. Although the three WCS strains are all capable of eliciting ISR, they show host specificity in terms of their ability to induce ISR in different plant species (Table 1). Also within a single plant species, the three WCS strains display differential effectiveness. For instance, in radish, WCS374 and WCS358 are potent elicitors of ISR, whereas WCS417 is not [25]. Arabidopsis thaliana (Arabidopsis) possesses natural genetic variation for the ability to express WCS417-ISR, a trait that could be mapped to a single genetic locus in the Arabidopsis genome, indicating that rhizobacteria-mediated ISR is genetically determined [37, 38].

Table 1 Induction of systemic resistance by Pseudomonas strains WCS358, WCS374, and WCS417

In the past 20 years, the Pseudomonas spp. WCS strains were highly instrumental in research on the molecular basis of rhizobacteria-mediated ISR signaling [13]. Using WCS417, it was demonstrated that the onset of ISR is regulated by the root-specific transcription factor MYB72 [39, 40]. Furthermore, it was shown that establishment of ISR in the leaves functions independently of the plant defense hormone salicylic acid (SA) but instead requires the plant hormones jasmonic acid (JA) and ethylene (ET) [41, 42]. These observations established that WCS417-ISR is mechanistically different from pathogen-induced systemic acquired resistance (SAR; [43]). Different bacterial determinants of the WCS strains have been implicated in the elicitation of ISR, including siderophores, lipopolysaccharides and flagella [4447]. However, while application of the purified determinants triggers ISR in the host plant, rhizobacterial knockout mutants that no longer produce the respective determinants were not impaired in their ISR-inducing ability [44, 48]. This suggests that bacterial determinants act redundantly in the elicitation of ISR.

The WCS strains can protect different plant species against infection by soil-borne and foliar pathogens, either directly by inhibiting pathogen growth via competition for iron, or indirectly by priming the plant immune system. For both mechanisms, sufficiently high population densities of these rhizosphere bacteria are essential [26]. The density of the WCS strains in the rhizosphere depends on their capacity to efficiently colonize the roots and effectively compete with other microbes in the rhizosphere for nutrients and space. Many of these nutrients are secreted by the plant roots and can manipulate the structure of the microbial community on the root system [2]. In turn, beneficial rhizobacteria can modulate the composition of root excretions and root immune responses, thereby influencing the structure of the root microbiome and the nature of the mutual benefits [40, 49, 50].

Despite the wealth of knowledge on the biology of the Pseudomonas spp. WCS strains in relation to their biological control and plant immunity-stimulating activity in host plants, their genomes were not yet elucidated. Here, we report on the whole genome sequences of the Pseudomonas strains WCS358, WCS374 and WCS417. We examined these three genomes in the light of the vast amount of functional data that have been published on these strains during the past 30 years. Moreover, we determined the taxonomic position of the WCS strains within the genus Pseudomonas and found that each of the three WCS strains is very closely related to one or more other plant-beneficial Pseudomonas strains with very different geographic origins. P. fluorescens WCS417 appeared to be taxonomically very similar to Pseudomonas simiae type strain R81 and was thus renamed into P. simiae WCS417. P. fluorescens WCS374 and P. putida WCS358 belong to the P. fluorescens and P. putida subgroups, respectively, but are taxonomically not within the 98 % identity range of any previously described type strain. Hence, WCS374 and WCS358 belong to thus far undescribed species, which we propose to name Pseudomonas defensor (type strain WCS374) and Pseudomonas capeferrum (type strain WCS358), respectively. Analysis of the genomes of WCS358, WCS374 and WCS417 in the light of their previously described biological functions revealed: 1) distinct strategies for siderophore-mediated iron acquisition among the three strains, 2) the structure of the WCS417 siderophore pyoverdine PVD417 and evidence that PVD417 is functionally different from the closely related pyoverdine of WCS374 PVD374; 3) insight into differences and similarities between the bacterial secretion systems and their relevance for rhizosphere competition and interactions with their host; and 4) bacterial determinants that may explain the differential abilities of the three strains to elicit ISR in different plants species. The whole genome sequences of WCS358, WCS374, and WCS417 in combination with the wealth of previously published knowledge on the biology of these PGPR strains will be highly instrumental for the research community working on the mode of action of plant-associated fluorescent Pseudomonas spp. in relation to their beneficial effects on plant growth and defense.


General genome characteristics

Sequencing of the genomes of WCS358 from the P. putida group, and WCS374 and WCS417 from the P. fluorescens group was carried out at the Beijing Genome Institute (Beijing, China). A summary of the general sequence characteristics is given in Table 2. The sizes of the genomes varied from approximately 5.94 Mb for WCS358 to 6.09 Mb for WCS374 and 6.17 Mb for WCS417 and they were predicted to contain 5188, 5351 and 5506 coding sequences, respectively. Their G + C contents ranged from 60 to 63.5 %. These characteristics are comparable to those described for taxonomically related and previously sequenced Pseudomonas spp. strains [51, 52].

Table 2 General sequencing and genome characteristics of the WCS Pseudomonas spp. strains

Phylogenetic analysis of WCS358, WCS374, and WCS417

The genus Pseudomonas is very diverse and is still undergoing taxonomic refinement [53]. Based on their nutritional and physiological characteristics, strains WCS374 and WCS417 were tentatively ascribed to the species Pseudomonas fluorescens, whereas WCS358 was described to belong to Pseudomonas putida [54, 55]. To obtain a comprehensive overview of the taxonomic position of the WCS strains, the sequences of the core housekeeping genes 16S rRNA, gyrB, rpoB and rpoD were compared with the corresponding genes of 107 Pseudomonas species type strains [53] and a selection of additional Pseudomonas spp. strains of which the full genome was available [51, 52, 5661]. A phylogenetic tree was generated based on multi-locus sequence analysis (MLSA) of the concatenated sequences of four core housekeeping genes (Fig. 1). This tree shows that both WCS374 and WCS417 are associated with the P. fluorescens subgroup as defined by Mulet and co-workers [53]. However, they are not most closely related to the P. fluorescens type strain. Similarly, WCS358 is related to the members of the P. putida group, but not most closely to the P. putida type strain. WCS417 has a 100 % nucleotide identity match with the concatenated sequences (NI) of the type strain of Pseudomonas simiae. WCS374 is most closely related to the type strain of Pseudomonas synxantha (96 % NI), whereas WCS358 shows highest homology with Pseudomonas monteilii (94 % NI). As the species boundary for MLSA was proposed to be 97 % [53], WCS417, which was previously ascribed to the P. fluorescens species, should be regarded as a P. simiae strain, whereas both WCS374 and WCS358 should be considered representatives of yet undescribed species. WCS417 appears to be closely related to the plant-beneficial P. simiae strain R81 (100 % NI), which was also isolated from wheat roots [60]. Furthermore, WCS374 appears to be closely related to two other plant-beneficial Pseudomonas spp. strains within its species boundary, i.e. strains A506 (99 % NI; [52]) and SS101 (95 % NI; [52]).

Fig. 1
figure 1

Phylogenetic tree showing the relationship of WCS358, WCS374, and WCS417 with other Pseudomonas spp. strains. Phylogenetic tree of the WCS strains (blue; also indicated by arrows) relative to 107 Pseudomonas sp. type strains (red), 24 selected Pseudomonas sp. strains of which the genomes were already sequenced (green), and Pseudomonas sp. strain CFBP2461 (purple). The tree is based on the alignment of concatenated sequences of four core housekeeping genes (16S rRNA, gyrB, rpoB, rpoD) of the strains. Bootstrap values from 1000 replicates are indicated at the nodes. Organization of Pseudomonas groups and subgroups is according to Mulet et al. [53]

To further establish the phylogenetic relationship of the WCS strains with other Pseudomonas spp. strains, we compared the whole genomes of our strains to available completely sequenced strains of the P. fluorescens subgroup and the P. putida group by calculating the Average Nucleotide Identity based on BLAST (ANIb) using JSpecies [62]. ANIb values of 94 % [63] and 95 % [62] have been proposed for defining the boundaries between species. The phylogeny based on ANIb values matched with that of MLSA, confirming the phylogenetic separation between strains of the P. putida group and the P. fluorescens subgroup (Fig. 2). None of the Pseudomonas strains for which the genome was available shared more than 94 % of its total genome sequence with WCS358, confirming that WCS358 belongs to a thus far undescribed Pseudomonas species within the P. putida group. WCS417 shared 99.8 % nucleotide identity with the draft sequence of Pseudomonas strain R81, indicating that these strains belong to the same species. Furthermore, Pseudomonas spp. strains A506 and SS101 shared 99.1 and 94.8 % of their respective genomes with WCS374 indicating that they belong to the same so far undescribed Pseudomonas species within the P. fluorescens subgroup (See also Additional file 1: Figure S1 and Additional file 2: Table S1). Together, our genomic analyses indicate that strain WCS417 belongs to the species P. simiae of the P. fluorescens subgroup, whereas strains WCS358 of the P. putida group and WCS374 of the P. fluorescens subgroup belong to sofar undescribed species. Here we designate these species as Pseudomonas capeferrum (type strain WCS358) and Pseudomonas defensor (type strain WCS374; see also Discussion section).

Fig. 2
figure 2

Average Nucleotide Identity based on BLAST for a selection of bacterial genomes. Black lines and bold numbers indicate putative species boundaries. Cell colors indicate similarity scaled from low (red) to high (green). ANIb values were calculated using Jspecies [62]

The whole-genome comparison unearthed that both WCS374 and WCS417 have a close relationship with other Pseudomonas spp. strains that were isolated based on their plant-beneficial properties. We aligned the genomes of WCS374 and WCS417 to their closest relative, i.e. Pseudomonas spp. strain A506 (isolated from pear phyllosphere, California, USA; [64]) and R81 (isolated from wheat rhizosphere in Bhawanipur, India; [60]), respectively, and found that large parts of the genomes were nearly identical, whereas only small, randomly dispersed genomic regions were missing in the corresponding relatives (Fig. 3). To further analyze the relationship between the pairs of strains, we used Islandviewer [65] to identify genomic islands in the compared genomes. Genomic islands are genomic regions that are derived from horizontal transfers [66]. Figure 3 shows a large overlap between genomic islands and the regions missing in the cognate closest relative, indicating that most of the differences between these closely related strains result from horizontal gene transfer. Although we did not find a close relative of WCS358 among the available sequenced Pseudomonas spp. strains, the plant-beneficial Pseudomonas sp. strain CFBP2461 (isolated from bean rhizosphere in Angers, France; [67]) was found to produce exactly the same pyoverdine siderophore as WCS358 [68]. Pyoverdines are iron-chelating molecules and their characteristics often closely reflect the phylogeny of the bacteria that produce them. Hence, in the past this feature was often used to characterize bacterial identity, a method known as siderotyping [69]. We amplified the housekeeping genes 16S rRNA, gyrB, rpoB and rpoD of CFBP2461 (Genbank accession numbers KM221193- KM221196) and included the concatenated sequence of this strain in an MLSA as described above (Fig. 1). Pseudomonas sp. CFBP2461 shared 99.5 % NI with WCS358, indicating that these strains are very similar and should be considered representatives of the same Pseudomonas species. Together, these results indicate that for each of the WCS strains other representatives of the same species have been found. The fact that these WCS “look-a-likes” have all been selected in independent searches for plant-beneficial microbes at very different geographic regions suggests that the WCS genome sequences can serve as representative genomes for plant-beneficial Pseudomonas spp.

Fig. 3
figure 3

Whole genome alignment of WCS374 and WCS417 with related Pseudomonas sp. strains. Whole genome alignments of WCS374 with its close relative Pseudomonas sp. strain A506 (a) and WCS417 with its close relative P. simiae R81 (b) were generated with Progressive MAUVE. Colored blocks indicate similar genome regions between the two strains. White gaps indicate genomic regions that are not shared between the compared strains. Genomic islands predicted by Islandviewer are indicated in purple above the black line for each genome. In R81, a large genomic region is represented indented from the rest of the genome as this region was reversely oriented in comparison to WCS417


The WCS strains produce siderophores that can inhibit plant pathogens directly via competition for iron, or indirectly via the onset of ISR. Application of purified pyoverdines of WCS358, WCS374 or WCS417 to the plant roots can elicit ISR. However, the pyoverdines of the WCS strains differentially trigger ISR in different plant species [25, 47] suggesting that their structure differs and that these pyoverdines are part of a specific host-microbe recognition system. Whilst the structures of the pyoverdines PVD358 and PVD374 of WCS358 and WCS374, respectively, have been elucidated, the structure of PVD417 of WCS417 is as yet unknown [68, 70]. Pyoverdines consist of a chromophore and a short peptide chain. Both are synthesized by the sequential action of multiple non-ribosomal peptide synthetases (NRPSs). NRPS are multi-modular enzymes in which each module is responsible for addition, attachment or modification of a specific amino acid onto a growing peptide chain [71]. AntiSMASH identified gene clusters that contained NRPSs in all three WCS strains. In each strain, two NRPS-containing gene clusters could be related to the production of pyoverdine. One of these gene clusters is likely to be involved in the synthesis of the chromophore as the only NRPS gene found in this cluster (data not shown) is an ortholog of NRPS pvdL of P. protegens Pf-5 [72]. The other pyoverdine biosynthesis gene cluster contains NRPS genes that are associated with synthesis of the respective pyoverdine peptide chains (Fig. 4a; indicated in green), pyoverdine transport (Fig. 4a; indicated in blue), or regulation (Fig. 4a; indicated in red). Bioinformatic analysis of the predicted NPRS-mediated peptide chains confirmed the previously elucidated composition of the pyoverdine peptide chains of PVD374 and PVD358 (Fig. 4b). In addition, the bioinformatics analysis predicts that the peptide chain of the WCS417 pyoverdine PVD417 contains the same amino acids as PVD374 in the same order (Ser-Lys-Gly-Orn-Lys-Orn-Ser; Fig. 4b).

Fig. 4
figure 4

Siderophore biosynthesis genes. (a) Siderophore biosynthetic gene clusters in the genomes of WCS358, WCS374, and WCS417 as identified by AntiSMASH. Colors represent different functional gene categories: biosynthetic genes (green); transport-related genes (blue); regulatory genes (red); and other genes (grey). (b) Bioinformatic analysis of the NRPSs that synthesize the peptide chain of pyoverdine in the WCS strains. For each NRPS, the domains recognized in the NRPS to function in condensation (C), adenylation (A), thiolation (T), epimerization (E), and epimerization and thioesterase (TE) are shown, as well as the amino acids predicted to be recognized by the adenylation domains

Recently, it was demonstrated that WCS374 produces exactly the same PVD as Pseudomonas strain A506 [32, 70]. The peptide chain of the WCS374/A506 PVD was demonstrated to be very similar to that of Pseudomonas sp. SBW25 as it contains the same amino acid residues in the same order. It was proposed that the SBW25 PVD differs slightly as its PvdJ contains an epimerization module that is lacking in the orthologous PvdD of the NRPS of WCS374 and A506 and that changed the stereo-isomeric configuration of the 6th amino acid in the chain from L-Orn to D-Orn [32]. Our bioinformatics analysis revealed that the PvdJ ortholog of WCS417 contains a similar epimerization domain that is lacking in WCS374 (Fig. 4b). This suggests that WCS417 produces the same PVD as SBW25, which is a close relative of WCS417 in our phylogenetic analysis (Figs. 1 and 2).

Pyoverdine knockout mutants WCS358-PVD and WCS417-PVD are completely abolished in siderophore activity in the universal siderophore assay on CAS agar medium (Fig. 5) [7375], indicating that both strains produce only one type of siderophore. However, pyoverdine knockout mutant WCS374-PVD still produced a halo on CAS medium (Fig. 5) [75], indicating that this strain produces one or more additional siderophores. Mercado-Blanco et al. [76], identified a second siderophore in WCS374, called pseudomonine (PSM). Pseudomonine is composed of salicylic acid (SA), cyclothreonine and histamine. Although microbially produced SA has iron-chelating properties and is implicated to function as a siderophore itself [77], the latter was questioned by other studies [78, 79]. Mercado-Blanco and co-workers [76] identified the pmsCEAB operon in WCS374 as responsible for production of SA and pseudomonine. Matthijs and co-workers [80] identified two NRPS genes upstream of pmsCEAB, i.e. basB and basAD, which are thought to function in the adenylation and cyclization of the amino acid threonine and subsequent assembly of pseudomonine [81]. All the genes required for pseudomonine biosynthesis were also detected in the genome sequence of WCS374 (Fig. 4a). Interestingly, mutants of WCS374 that do not produce PVD374, PSM or SA, still produce a halo in the CAS assay and are able to grow on KBA with 800 μM of the iron chelator bipyridyl, whereas the pyoverdine-deficient mutants WCS417-PVD and WCS358-PVD can only tolerate 400 μM of 2,2-bipyridyl (Fig. 5). This indicates that WCS374 produces yet another siderophore. Searching the WCS374 genome for additional siderophore biosynthesis genes using antiSMASH revealed a gene cluster (PD374_19610 - PD374_19665) involved in the production and transport of an aerobactin-like siderophore (Fig. 4a). This gene cluster is also present in the genomes of the WCS374 relatives SS101 and A506, although the ortholog of PD374_19640, which is related to IucC/IucA and is putatively required for the production of the siderophore, contains a frameshift mutation in A506.

Fig. 5
figure 5

Overview of siderophore production by the WCS strains and their mutant derivatives. Maximum amount of 2,2-bipyridyl in KBA that still allowed growth of strains after 48 h is indicated. Photographs display production of an orange halo by the bacterial strains on CAS agar, which is indicative for siderophore production. Abbreviations stand for Wild type (WT), pyoverdine (PVD), pseudomonine (PSM), pyochelin (PCH) or salicylic acid (SA). Production of pyochelin by WCS417 has not been demonstrated

Leeman et al. [25] investigated the production of SA by the three WCS strains and found that besides WCS374 also WCS417 is capable of producing SA. To identify the corresponding biosynthesis gene cluster in WCS417, the WCS374 genes pmsC and pmsB were used as bait in a BLASTp search of the WCS417 genome. This revealed two WCS417 orthologs (PS417_13200 and PS417_13205) in a genomic region that was identified by antiSMASH as being putatively responsible for the production of a pyochelin-like siderophore (Fig. 4a). Like pseudomonine, pyochelin is an iron-chelating siderophore found in Pseudomonas spp. that comprises a SA moiety [82]. However, pyoverdine knockout mutant WCS417-PVD does not show any siderophore activity on CAS agar at 28 °C, nor can it tolerate higher levels of bipyridyl than the pyoverdine-deficient mutant WCS358-PVD. Thus, it is unlikely that a pyochelin-like siderophore is indeed produced by WCS417 under the growth conditions tested.

Microbially produced siderophores are re-absorbed after they have complexed ferric iron from the environment. This is mediated by TonB-dependent proteins (TBDPs) that specifically recognize and transport siderophore-iron complexes into the periplasm [30]. We identified 40 TBDP-encoding genes in WCS358, 31 in WCS374 and 33 in WCS417. A minority of these 104 TBDPs also contain a short N-terminal domain, which is typical for TBDP-mediated transduction of environmental signals to the cytoplasm [83]. Most bacteria contain less than 14 TBDPs in their proteome, although some bacteria have larger numbers [32]. Recently, Hartney and co-workers [32] demonstrated that P. protegens Pf-5 contains 45 TBDPs. Six of these TBDPs (FpvU to FpvZ) function as ferric-pyoverdine receptors (FPVs) that facilitate the uptake of specific pyoverdine-iron complexes and enable this bacterium to use not only its native pyoverdine, but also heterologous pyoverdines from other strains. Previously it was reported that TBDPs cluster according to their substrate rather than to phylogeny [72]. Because, the ferric-pyoverdine receptors form a clear clade within the TBDPs [32, 84], we used the corresponding Fpv genes of strain Pf-5 to identify pyoverdine receptor genes in the three WCS strains. The amino acid sequences of the identified putative pyoverdine receptors were aligned together with the amino acid sequences of the six previously described pyoverdine receptors of Pf-5, after which a phylogenetic tree was build (Additional file 3: Figure S2). The six Pf-5 FPV pyoverdine receptors clustered together, confirming previous findings [72] and this cluster included 10 TBDPs of WCS358, five of WCS417 and four of WCS374.

Using combinations of deletion mutants of the six Fpv genes of Pf-5 it is possible to identify which FPV pyoverdine receptors are required for the uptake of specific ferric-pyoverdines [32]. Previously, it was demonstrated that Pf-5 can utilize ferric-pyoverdine of WCS374 and A506 through the receptors FpvU and FpvY, while ferric-pyoverdine of SBW25 can only be utilized through receptor FpvU. This discrepancy was attributed to the different stereo-isomeric configuration of the ornithine residue on position 6 of the pyoverdine peptide chain (Fig. 4b). Because from the genome sequence of WCS417 we predicted that PVD417 would be similar to pyoverdine of SBW25, we set out to functionally validate this using the set of Pf-5 fpv mutants [32]. To this end, the six fpv deletion mutants of Pf-5 were tested under low-iron conditions for their capacity to be cross-fed by WCS358, WCS374, or WCS417. In line with the descriptions of the specific nature of PVD358 [33, 68], none of the fpv mutants of Pf-5 could use the pyoverdine produced by WCS358 (Additional file 4: Figure S3). However, all mutants were able to grow in the presence of WCS374, confirming previous findings by Hartney and co-workers [32]. In the presence of WCS417, all fpv mutants of Pf-5 were able to grow, except mutant fpvU, indicating that pyoverdin receptor FpvU is required for the uptake of ferric-PVD417, which resembles the findings for SBW25. These results confirm the prediction from our bioinformatical analysis that WCS417 produces the same pyoverdin as SBW25 and highlights that slight differences in pyoverdine structure can have implications for specificity in heterologous uptake of iron-pyoverdin complexes.


Lipopolysaccharides (LPS) are molecules in the outermembrane of Gram negative bacteria, and are recognized as microbe-associated molecular patterns (MAMPs) that are able to elicit immune responses in plants and animals [85, 86]. Application of purified LPS of WCS417 to the roots of carnation plants was shown to trigger ISR, resulting in reduced Fusarium wilt disease when Fusarium was inoculated in the stem [87]. In radish, both WCS374 and WCS417 can trigger ISR, whereas WCS358 cannot. This differential effectiveness of ISR inducibility could be attributed to strain-specific differences in LPS, because in radish purified LPS of WCS417 and WCS374 triggered ISR, whereas that of WCS358 did not [88]. Bacterial LPS usually consists of three domains [86]. The first domain, lipid A, consists of a bisphosphorylated glucosamine-disaccharide backbone substituted with several fatty acids, which anchors the LPS to the outer membrane. Attached to lipid A is the LPS core, an oligosaccharide of about 9–10 sugars, which may be extended with the O-antigen or O-chain. The O-chain is a polysaccharide comprised of repeating units. The number and nature of sugars in the O-chain units is highly specific and can differ dramatically even between strains of the same species [86]. Under conditions of high iron availability, mutants of WCS374 and WCS417 that lack the O-antigen of the LPS (WCS374-ΔOA and WCS417-ΔOA) were no longer able to elicit ISR in radish, demonstrating that under these conditions the LPS of these strains is the only bacterial determinant implicated in ISR and that the O-antigen is the active component [88]. The importance of the O-chain of the LPS was also demonstrated in Arabidopsis. In this plant species, strains WCS417 and WCS358 are able to elicit ISR against Pseudomonas syringae, whereas strains WCS374 cannot [48]. Using LPS-containing cell envelopes of these strains, and mutants lacking the O-antigen it was shown that the O-chains of WCS417 and WCS358 are ISR-eliciting bacterial determinants, and that the O-antigen of the LPS of the three WCS strains are differentially recognized in Arabidopsis and radish.

In order to identify genes of the WCS strains involved in the biosynthesis of the highly variable O-antigen of LPS, we searched for putative O-antigen biosynthetic loci (OBL) using different OBL described for P. aeruginosa as bait. In P. aeruginosa, 20 distinct O-antigen serotypes are known for which the biosynthetic loci have been sequenced [89]. The WCS genomes were mined for orthologs of the predicted proteins in these 20 P. aeruginosa OBL. In each of the WCS strains, only one gene cluster was found with genes coding for more than four orthologs of the P. aeruginosa OBL proteins (Fig. 6). In P. aeruginosa, an OBL encoding the major enzymes for O-antigen biogenesis is usually found in between the highly conserved genes himD and wbpM. Orthologs for these two genes were also found in all the OBL of the WCS strains, although in the WCS strains genes encoding orthologs of the OBL proteins were also found upstream of himD (Fig. 6). The region upstream of himD is highly conserved between WCS417 and WCS374, but contains an ortholog of rspA that encodes an essential component of the protein synthesis machinery of Escherichia coli [90] and is most likely not involved in LPS biosynthesis.

Fig. 6
figure 6

LPS O-antigen biosynthetic loci. Putative LPS O-antigen biosynthetic loci (OBL) in the genomes of WCS358, WCS374, and WCS417 as identified in a BLASTp search using 20 OBL identified in P. aeruginosa as bait. Genes of which the gene tags are designated in red font are orthologs of genes found in at least one of the 20 P. aeruginosa OBL. Colors of the arrows indicate shared orthologs in the three WCS strains as determined with reciprocal BLASTp. Graphics under the arrows represent a measure of similarity as determined with Progressive MAUVE

Previously, the composition of the O-antigens of WCS374 and WCS358 were investigated biochemically and reported to consist of glucose and fucose (WCS374), and glucose and quinovosamine (WCS358), respectively [91]. However, the predicted function of the genes in the OBL of WCS374 and WCS358 are not clearly related to the biosynthesis of such subunits. Furthermore, the complexity of especially the OBL of WCS358 suggests that the O-antigen composition is not as clear-cut as described by De Weger and co-workers [91]. In contrast to the O-antigens of WCS358 and WCS374, the composition of the WCS417 O-antigen has not been previously explored. The genome sequence of the WCS417 OBL gives a number of clues towards the structure of the WCS417 LPS O-antigen. The WCS417 genes with locus tags PS417_08220 – PS417_08245 are homologs of ddhA-ddhD, yerE, and yerF, which are required for the synthesis of the O-antigen subunit yersiniose in Yersinia pseudotuberculosis [92, 93] a sugar that has been detected before in a Pseudomonas spp. O-antigen [94]. The presence of the genes with locus tags PS417_08165 – PS417_08175, which encode orthologs of WbjB, WbjC and WbjD that are required for the synthesis of the O-chain component N-acetyl-L-fucosamine (L-FucNAc) [95], suggests that L-FucNAc is also part of the WCS417 O-antigen. In addition, an ortholog of wbpM of P. aeruginosa is also found in the WCS417 OBL. WbpM catalyzes a reaction that ultimately results in either 2-acetamido-2-deoxy-D-fucose (D-FucNAc) or 2-acetamido-2-deoxy-D-quinovose (D- QuiNAc) [86]. It is therefore likely that the WCS417 LPS O-antigen contains yersiniose, L-FucNAc and D-QuiNAc or D-FucNAc. Within the WCS417 OBL, only four genes are shared between WCS417 and WCS358 (Locus tags PS417_08165, PS417_08175, PS417_08185 and PS417_08190) and only one of them (PS417_08190) is also present in WCS374 (Fig. 6). This lack of homology between the OBLs of the three WCS strains is likely to result in variation in O-chain composition, which may be causally related to the observed host-microbe specificity of LPS in eliciting ISR.

Antimicrobial compounds

Meziane and co-workers [44] investigated the redundancy of ISR-eliciting bacterial determinants in different plant species. They found that WCS358 could elicit ISR against Botrytis cinerea in tomato. This ISR could also be triggered by applying purified LPS or pyoverdine to the roots but not by applying flagella of WCS358. Nonetheless, a double knockout derivative of WCS358 that did not produce the LPS O-antigen or PVD358 still elicited ISR. This indicated that at least a third WCS358 determinant is able to elicit ISR in tomato [44]. In addition to siderophores, LPS and flagella, antimicrobial compounds are bacterial determinants that can induce systemic resistance [45, 9698]. However, to date no antimicrobial metabolites have been described for the three WCS strains under investigation. Tran and co-workers [96] demonstrated that the cyclic lipopeptide massetolide A, produced by Pseudomonas strain SS101, can induce ISR. Like siderophores, cyclic lipopeptides are non-ribosomally produced by NRPSs. All NRPS-containing gene clusters identified by antiSMASH in the genomes of the three WCS strains could be related to the production of the siderophores described above, except for one gene cluster in WCS358 (PC358_04000 - PC358_4180). This gene cluster putatively encodes three NRPS that share 100 % identity in a BLASTp search against PsoA, PsoB and PsoC. PsoA-C are responsible for the biosynthesis of the cyclo lipopeptides putisolvin I and II in P. putida PCL1445 [99]. In PCL1445, three additional genes have been shown to be important for putisolvin production. The LuxR-type regulator gene psoR is located upstream of psoA and is required for expression of the pso cluster, whereas macA and macB, located downstream of psoC, are involved in putisolvin production or export [99, 100]. Interestingly, also the entire psoR-macB cluster of PCL1445 shares 100 % nucleotide identity with WCS358. Tot test if WCS358 is indeed capable of producing these cyclic lipopeptides, a drop collapse assay was performed. In this assay, WCS358 was able to reduce surface tension, whereas WCS374 and WCS417 were not, demonstrating that WCS358 indeed produces a surfactant (Additional file 5: Figure S4).

Besides the putisolvin biosynthesis cluster, other genes that are potentially involved in the synthesis of compounds with broad-spectrum antibiotic activity could not be identified in the genomes of WCS358, WCS374 and WCS417. For instance, biosynthesis genes for the antimicrobial compounds 2,4-diacetylphloroglucinol, phenazines, hydrogen cyanide, and pyrrolnitrin, which are abundantly present in root-associated Pseudomonas strains [52, 82], are not present in the genomes of the sequenced WCS strains. This corroborates with early observations that the in vitro antagonistic activity of these strains was only apparent at low-iron conditions, which suggested that the observed antagonistic activity of the WCS strains is predominantly based on siderophore-mediated competition for iron [20, 21]. However, some genes were detected that encode putative bacteriocins. Bacteriocins are bacteriocidal proteins that are generally effective against a narrow taxonomic range of bacteria closely related to the producer [101]. We mined the WCS genomes with BAGEL3 [102] for orthologs of bacteriocin-encoding genes that were previously identified in taxonomically related strains [52]. Four putative bacteriocins were identified in WCS417, one related to S-type pyocins of P. aeruginosa (PS417_07930, pyocin AP41-like) and three to R-type pyocins (PS417_05796, PS417_10225 and PS417_05705). In WCS374, three bacteriocins were identified: again one related to S-type pyocins (PD374_22800, pyocin S6-like) and two to R-type pyocins (PD374_05705 and PD374_06430). In WCS358, no bacteriocins could be identified. Together, the genomes of WCS358, WCS374, and WCS417 indicate that these Pseudomonas strains produce a relatively small pallet of known antibiotics, confirming that their plant protective capacity is rather based on other mechanisms (e.g. siderophore-mediated competition for iron and ISR).

Protein secretion systems

In order to deal with different environments, competing microbes and accommodating hosts, bacteria need to secrete enzymes and other proteins into the extracellular environment [103]. Five different protein secretion systems (type I, II, III, V, and VI) are typically found within the Pseudomonas genus [104]. Differences in these protein secretion systems are likely to be important for the excretion of proteins involved in traits that influence rhizosphere competence of the WCS strains and for delivery of putative ISR-eliciting determinants.

Type I secretion systems consist of three components that together span the bacterial cell envelop and transport their products from the cytosol directly to the extracellular environment [104]. To identify gene clusters in the WCS genomes that encode the Type I secretion system (T1SS), a BLASTP search was performed with protein sequences encoded by the P. aeruginosa aprDEF gene cluster as bait. Two complete T1SS loci were found in the WCS417 genome (Additional file 6: Table S2). One T1SS gene cluster contained genes encoding a lipase and an ortholog of the alkaline protease AprA [105], whereas the second T1SS gene cluster contained a gene encoding an ortholog of the hemophore HasA [106]. It is likely that the T1SSs are dedicated to the secretion of these proteins that are encoded in their gene cluster. Also in WCS374, the genes for two complete T1SSs were found of which one also shared an operon with aprA and a lipase. For the second T1SS in WCS374 and for both T1SSs found in WCS358, the substrate was not so obvious, as the corresponding gene clusters did not provide strong clues (Additional file 6: Table S2). For instance, no ortholog of AprA could be found in the WCS358 genome. To experimentally confirm this, we used a milk powder assay to detect AprA protease activity. WCS374 and WCS417, which possess an aprA gene in one of their T1SS gene clusters, both showed a clearing zone around their colonies when grown on KBA supplemented with milk powder. For P. syringae pv. tomato DC3000, such a clearing zone could be attributed to AprA activity as its aprA mutant did not show a halo (Additional file 7: Figure S5). Also WCS358 lacked such extracellular protease activity, confirming previous observations [107].

Unlike the T1SS, proteins secreted by type II secretion systems (T2SS) do not directly transfer the secreted proteins over the entire envelop, but make a stopover in the periplasmic space [104]. To identify gene clusters that encode the Type II secretion system (T2SS) that transports proteins across the outermembrane, a BLASTp search was performed using protein sequences encoded by the P. aeruginosa xcpP-xcpZ and hxcP-hxcZ gene clusters. Clusters containing four or more genes with significant similarity to the bait genes (e-value <10−5) were considered T2SS loci. In this way, two T2SS loci were found in the WCS417 genome, whereas the WCS374 genome seemed to contain only one T2SS locus (Additional file 6: Table S2). Also in WCS358, only a single T2SS locus was found using the P. aeruginosa T2SSs as baits. This WCS358 T2SS locus resembled the Xcp T2SS of P. aeruginosa PA01, which was recently found to secrete the PhoX-type phosphatase UxpB under phosphate-limited conditions, thereby stimulating the phosphate uptake machinery of the bacterium [108]. A second T2SS of WCS358 was detected using the Xcm cluster described in P. putida GB-1 and KT2440 as bait (88 % nucleotide identity shared with the cluster of GB-1; [109]). Activity of this T2SS has not been confirmed in WCS358.

There is a strong resemblance between T2SS and the machinery required for type IV pili biogenesis [110]. In WCS358, type IV pilus genes were previously investigated, because type IV pili might have a role in attachment of bacteria to the root surface [111]. Even though type IV pili could not be detected on the surface of WCS358, a cluster containing orthologs of the pilA, pilC and pilD genes of P. aeruginosa PAO1 could be identified. However, a gene for an ortholog of the traffic ATPase PilB was lacking in this gene cluster. In the complete genome of WCS358, we could detect five genes encoding proteins with significant homology to P. aeruginosa PAO1 PilB. However, the PilB ortholog with the highest similarity (PC358-2543; 43 % AA identity) was found to be XcpR encoded in the WCS358 Type II secretion xcp gene cluster. A PilB ortholog is apparently not present in WCS358, which explains the lack of type IV pili in WCS358.

The type III secretion system (T3SS), which consists more than 20 proteins, is the most complex of all known protein secretion systems found in bacteria [104]. In the genomes of WCS417 and WCS374, ORFs encoding structural and regulatory components of the T3SS are organized in 26-kb and 18-kb clusters, respectively, and several ORFs within each cluster display a significant degree of similarity to the hrp/hrc cluster of pathogenic bacteria. Following the nomenclature initially proposed by Preston and co-workers [112] for the type three secretion (tts) gene cluster of P. fluorescens SBW25, we named these genes rsp (rhizosphere-expressed secretion protein) or rsc (rsp conserved). No tts gene cluster or individual ORFs encoding T3SS components could be identified in the WCS358 genome, suggesting that this strain lacks a typical T3SS. Immediately adjacent to the rspL regulatory gene of the WCS417 tts cluster, a large ORF is present that encodes RopE of the AvrE family of effectors. Effectors of the AvrE family can suppress the plant’s basal immune responses and promote cell death of the host [113]. To identify additional type III effectors in the WCS417 and WCS374 genomes, we employed bioinformatical analyses searching for conserved Hrp (Rsp) “box motifs” in the promoter regions of putative effectors and exploring the N-terminal protein sequence of candidate effectors for features typical of type III secreted effectors (i.e. abundance of Ser and polar residues, acidic residues in the first 12 positions, and an aliphatic amino acid in position 3 or 4) [114, 115]. Based on these criteria, we identified 11 putative effectors for WCS417 and 15 putative effectors for WCS374 (Additional file 8: Table S3), whereas not a single putative effector could be identified in WCS358. Interestingly, the majority of WCS417 and WCS374 effectors show no homology to known effector families (data not shown), suggesting that these proteins represent novel effectors.

Although recently the genes for a Type IV secretion system (T4SS) were found in a Pseudomonas genome [61], T4SSs are uncommon in Pseudomonas spp. and were not detected in the WCS strains.

The type V secretion system (T5SS) is the most simple of all Gram-negative bacterial secretion systems as the secreted proteins are transported across the outer membrane with the aid of a covalently connected translocator domain (autotransporters) or via a single dedicated outer membrane protein (called two-partner secretion or TPS) [104]. The C-terminal translocator domain of T5SS (PFAM03797) is conserved among different classical monomeric autotransporters (T5aSS). In total, we found 2, 10, and 5 T5aSS autotransporters with this domain in WCS358, WCS374, and WCS417, respectively (Additional file 6: Table S2). Although the functions of these autotransporters are unclear, some are homologous to autotransporters of P. aeruginosa PAO1 or other mammalian pathogens with functions in host immune activation [116, 117], peptidase activity [118], or biofilm formation and motility [119, 120].

In addition to the T5aSS, the WCS genomes were mined for orthologs of ORFs PA0692, PA4540 and PA4624, which code for the outer membrane components of the T5bSS two-partner secretion systems in P. aeruginosa PAO1. Nine orthologs were found in the WCS genomes: 2 in WCS358, 4 in WCS374, and 3 in WCS417. The orthologs found in this way all contained a periplasmic polypeptide transport-associated (POTRA) domain required for recognition of the substrate protein TpsA that is transported by the T5bSS. Although significant orthologs for P. aeruginosa PAO1 tpsA genes could not be found in the WCS genomes, six putative tpsA genes could be identified that contained the “haemagglutinin activity domain” (pfam05860), which specifically interacts with the POTRA domains in the TpsB component [104]. The tpsA genes were found in the same operon and adjacent to a tpsB gene that is likely dedicated to its secretion. Many TpsAs are toxins that play a role in contact-dependent growth inhibition of neighboring bacteria [121]. A clear example is PD374_00815, which carries a DUF637 and a pre-toxin Hint (PT-Hint) domain, which are characteristic for such toxic TpsAs. It is therefore likely that TpsAs enable the WCS strains to inhibit their competitors and successfully colonize the plant roots. For three of the putative TpsBs, no TpsA component could be detected. To our knowledge, such stand-alone TpsB proteins have not been previously described.

Using YadA and invasin of Yersinia enterocolitica as bait in BLAST searches, we did not detect members of the trimeric autotransporters (T5cSS) or intimin/invasin family (T5eSS) in the WCS strains. However, in each of the WCS genomes, a single ortholog of Patatin-like protein PlpD of P. aeruginosa PAO1 was found, which was recently suggested to represent a novel T5SS (T5dSS) [121]. Patatins form a group of glycoproteins with lipolytic activity found in potato tubers and have proposed functions against plant pathogens [104, 122, 123].

Type VI secretion systems (T6SS) are functionally similar to T3SS in that they deliver effector-like proteins into other organisms. Therefore, like T3SS, T6SS have been implied to play a role in manipulation of host immunity. They also play a prominent role in bacterial warfare by delivering toxic proteins into competing micro-organisms [124]. The machinery structurally resembles the contractile tails of bacteriophages and is used, in this case, for injecting toxins into target cells. Protein sequences corresponding to the genes encoded in the T6SS loci of P. aeruginosa PAO1 were used as bait in BLASTp searches of the three WCS genomes. Clusters containing five or more genes with significant similarity to the bait genes (e-value < 10−5) were considered a T6SS locus [125]. The genomes of WCS417 and WCS374 each contained two T6SS loci, while WCS358 contained one.

Overall, these data display the diversity of the protein secretion systems in WCS358, WCS374, and WCS417. Their presence is not surprising, but highlight the sophisticated mechanisms that these plant-beneficial rhizosphere bacteria have evolved in order to successfully compete with other soil microbiota and sustain a long-term mutualistic relationship with their host plants.


In the past 30 years, the plant-beneficial Pseudomonas model strains WCS358, WCS374 and WCS417 have been extensively studied for their plant-protective traits. In this study, we explored their genomes and established an overview of their phylogeny and traits supporting rhizosphere competence and plant functioning (Fig. 7).

Fig. 7
figure 7

Selected rhizosphere competence traits in WCS strains. Overview of selected traits that contribute to rhizosphere competence of WCS358, WCS374, and WCS417. Traits are divided in groups related to siderophore-mediated competition for iron, antimicrobial compound production, exoenzymes, and protein secretion systems. Colored boxes indicate the presence of a gene or gene cluster. Absence of a box means absence of the gene or gene cluster. Numbers within boxes indicate the number of copies of genes or gene clusters in the corresponding genome

WCS strains are representatives of universally-present plant-beneficial Pseudomonas spp.

The genus Pseudomonas is very diverse and its taxonomy is still subjected to refinement [53]. Ideally, the taxonomy of Pseudomonas spp. reflects their evolutionary history and is predictive for their functionality [126]. Previously, WCS358 was classified as P. putida and both WCS374 and WCS417 as P. fluorescens [54, 55]. Although these WCS strains are indeed closely related to the P. putida and P. fluorescens type strains, our genome analysis revealed that the level of similarity is not sufficiently high to group the WCS strains within these species. Instead, WCS417 was found to be very closely related to the P. simiae type strain and can be regarded as a representative of this species. The level of sequence homology of WCS358 and WCS374 did not fall within the boundary of any other Pseudomonas type strain. Hence, WCS358 and WCS374 should be regarded as a yet undescribed species. Genome information from Pseudomonas strains A506 and SS101 revealed that these strains are very similar to WCS374 and fall within the boundary of the same so far undescribed Pseudomonas species. WCS374, A506 and SS101 were all isolated from plants, albeit from different hosts (potato, pear and wheat, respectively), and different geographic locations in the Netherlands and Oregon [20, 52, 64, 127]. Since all three strains have been selected and studied for their abilities to protect plants against a wide range of pathogens, we propose the name Pseudomonas defensor for this species and strain WCS374 as its type strain. Strain WCS358 was also studied extensively because of its ability to protect plants against pathogen infection. Because the plant protective properties of WCS358 are predominantly linked to siderophore production and sequestration of iron, we propose the name Pseudomonas capeferrum (‘seizes iron’) for the species and WCS358 as its type strain. This species also includes Pseudomonas strain CFBP2461, which was isolated from roots of bean in France. Like WCS358, CFBP2461 has also been studied for its plant-protective capabilities and produces exactly the same siderophore as WCS358 [68].

For each of the WCS strains we found one or more other representatives of the same species, all of which surfaced in the Pseudomonas research community because of their plant-beneficial properties. The strains were isolated from different hosts and at different geographic locations, indicating that the sequenced WCS strains can be regarded as universal representatives of a group of Pseudomonas spp. with plant protective functions. The genomes of both WCS374 and WCS417 were so similar to the genomes of their respective closest relatives that they are nearly isogenic, even though these strains were isolated decades and continents apart. This is remarkable considering the enormous diversity of Pseudomonas species in the soil and the relatively small number of bacterial genomes that have been sequenced [128, 129]. The fact that different labs independently isolate very closely-related bacterial strains with similar plant protective properties at different locations across the World, suggests that plants have mechanism to specifically enrich their microbiome for such plant-beneficial microbes.

WCS genomes support different strategies for iron acquisition

Originally, the WCS strains were studied for their remarkable capacity to control plant pathogens through siderophore-mediated competition for iron. The WCS strains have evolved different strategies for the acquisition of iron. WCS358 seems to have limited options as it can only produce one siderophore, PVD358. However, very few Pseudomonas strains seem to be able to use this pyoverdine. Only 1 % of over 1000 randomly isolated Pseudomonas strains could utilize PVD358, whereas 30 % of these strains were able to use PVD374 [33]. In addition, WCS358 is able to utilize the pyoverdines produced by a wide range of other Pseudomonas strains [33], which is corroborated by the large number of pyoverdine receptors found in the genome of WCS358 (Fig. 7).

By contrast, WCS374 produces a less specific pyoverdine, and is not able to use as many heterologous siderophores as WCS358. The latter is reflected by the lower number of putative pyoverdine receptor genes in the genome of WCS374 (Fig. 7). This suggests that many other bacteria are able to utilize PVD374, while in turn WCS374 is not well capable of utilizing pyoverdine of competing bacteria. However, WCS374 can produce a different siderophore, pseudomonine. The competitive success of WCS358 and WCS374 has been investigated in the rhizosphere of radish [130]. After co-inoculation, WCS358 reduced population densities of wild-type WCS374 by up to 30-fold, whereas WCS374pMR, a WCS374 derivative harboring the PVD358 receptor pupA of WCS358 maintained its population density. Population densities of WCS358 were unaffected by co-inoculation with WCS374 or any of 20 random radish-rhizosphere isolates. This demonstrates that pyoverdine-mediated competition for iron is a major factor in the colonization of radish roots by these strains and that the iron acquisition strategy of WCS358 is more advantageous than that of WCS374.

The iron sequestration strategy of WCS417 is intermediate to that of WCS358 and WCS374. It has fewer ferric pyoverdine receptors than WCS358, but one more than WCS374 (Fig. 7). We discovered a gene cluster in the WCS417 genome that is involved in the putative production of a pyochelin-like siderophore, indicating it produces a siderophore additional to pyoverdine. Moreover, although the WCS417 siderophore PVD417 is very similar to PVD374, the epimerization module of the NRPS involved in biosynthesis of PVD417 is predicted to modify the configuration of one of the six amino acid residues in the peptide chain of PVD417. It is tempting to speculate that this is an evolutionary adjustment to decrease the number of other bacteria that are capable of using PVD417.

ISR elicitors encoded by the WCS genomes

The production of pyoverdines has been demonstrated to play an important role in the biological control activity of the WCS strains in the rhizosphere. Flax, potato, radish, eucalypt, carnation, bean and tomato were less diseased when their roots were colonized by one of the WCS strains, while this protection was no longer observed when pyoverdine knockout derivatives of the strains colonized the roots of these plants [22, 25, 27, 28, 44, 75, 131]. However, this protective effect is not always solely based on siderophore-mediated competition for iron, as pyoverdines were also implicated in the onset of ISR. PVD358 has been demonstrated to elicit ISR in Arabidopsis, bean, eucalypt and tomato [44, 131]. In tobacco, PVD358, PVD374 and PVD417 have all been demonstrated to trigger ISR [132]. In radish, only PVD374 is an effective ISR elicitor, whereas PVD358 and PVD417 are not [25]. This suggests that siderophores may play a role in host specificity with respect to the ability of the rhizobacterial strains to elicit ISR. Interestingly, our genome analysis predicts that PVD374 and PVD417 only differ in the isomeric configuration of one of the amino acids incorporated in the otherwise-identical peptide chain suggesting that plants can distinguish such minute differences between pyoverdines.

Besides pyoverdine, WCS374 also produces the siderophore pseudomonine, whereas the other WCS strains do not (Fig. 7). Interestingly, pseudomonine has a SA moiety in its molecule, the biosynthesis of which is encoded by the isochorismate synthase gene pmsC and the isochorismate-pyruvate lyase gene pmsB. In WCS374, both genes are part of the pseudomonine biosynthesis gene cluster pmsCEAB [76]. SA has been demonstrated to be a potent inducer of systemic acquired resistance (SAR). However, while exogenously applied SA induces systemic disease resistance in Arabidopsis, the SA producer WCS374 does not [48]. It was proposed that during colonization of plant roots by WC374, SA becomes swiftly incorporated in the siderophore pseudomonine, thereby preventing free SA from acting as an elicitor of plant immunity [70]. Interestingly, the close relative of WCS374, Pseudomonas strain SS101, is capable of inducing a systemic resistance in Arabidopsis, and it does so in an SA-dependent manner [133]. The SA biosynthesis genes of SS101 are highly similar to those of WCS374, but it may be possible that SS101 produces more SA in the rhizosphere or incorporates it less efficiently in its pseudomonine, thereby allowing the activation of the SAR pathway.

The nature of other potential ISR elicitors produced by the WCS strains has most extensively been investigated for WCS358 [44]. In Arabidopsis, purified pyoverdine, LPS and flagella of WCS358 were all demonstrated to be able to trigger ISR. Interestingly, knockout mutants of WCS358 that no longer produced one of these determinants also triggered ISR, which was concluded to be a logical consequence of the redundancy of these determinants. In tomato, the cyclic lipopeptide massetolide A of SS101 emerged as an additional elicitor of ISR [96]. The cyclic lipopeptide predicted to be produced by WCS358 might similarly be involved in ISR. Cyclic lipopeptides are also important in rhizosphere competence as they can inhibit other microbes and influence biofilm formation [71]. In this study, we were able to gain insight into the nature of the biosynthetic genes of all these ISR elicitors, which will be highly instrumental in future studies on the molecular mechanisms of the onset of ISR by plant-beneficial rhizosphere bacteria.

Protein secretion systems encoded by the WCS genomes in relation to rhizosphere competence

To protect plants from harmful organisms, PGPR need to reach sufficiently high population levels in the rhizosphere. For siderophore-mediated competition for iron or the onset of ISR, cell densities of at least 105 per gram of root have to be reached by the WCS strains [26]. Similar threshold densities have been reported for other PGPR [134, 135]. This requires that the PGPR strains are capable of out competing other microorganisms. To date, six types of protein secretion systems are recognized in bacteria, most of which play a role in host-microbe and microbe-microbe interactions. The Type I, II, III, V, and VI secretion systems were represented in the genomes of all three WCS strains (Fig. 7). The T1SSs and T2SSs of the WCS are likely to be involved in the secretion of extracellular enzymes, like extracellular protease AprA, lipase and phosphatase to facilitate nutrient acquisition [108, 134]. In addition, T1SS and T2SS may facilitate the secretion of antimicrobial compounds, such as cyclic lipopeptides [71], or bacteriocins [101] for which we identified biosynthetic genes in the genomes of the WCS strains. The genomes of WCS374 and WCS417 possess a gene cluster that encodes a putative T3SSs, as well as several genes encoding so far unknown putative effector proteins. T3SSs have been demonstrated to be functional in strains of the P. fluorescens group [112, 136, 137, 138]. However, their ecological significance in the rhizosphere is still unclear, as t3ss mutants of several Pseudomonas strains are not hampered in their ability to retain high population levels in the rhizosphere [112, 137, 139]. For the T5SS, a wide diversity of T5SS-related genes were identified in the WCS genomes. Many proteins that are secreted via the T5bSS are toxins that can play a role in contact-dependent growth inhibition of neighboring bacteria [121], thereby contributing to rhizosphere competence of the WCS strains. The same holds true for the T6SSs that were identified in the WCS genomes. T6SS have been demonstrated to function in competition with other microbes by delivering bactericidal effectors in a cell-contact dependent manner [103, 140, 141].

Suppression of root immune responses

To establish a long-lasting mutual relationship with their hosts, mutualistic microbes need to suppress or evade local host immune responses that are triggered in plant roots upon recognition of alien organisms [138]. For WCS417 it has been demonstrated that colonization of Arabidopsis roots by this rhizobacterium suppresses local root immune responses that are typically triggered by the flagellin epitope flg22 [50]. A possible mechanism by which WCS417 prevents the activation of this innate immune response is through the production of the alkaline protease AprA. AprA of P. aeruginosa and P. syringae was recently demonstrated to degrade flagellin monomers, thereby preventing flg22-triggered immunity in the leaves of Arabidopsis [142, 143]. The aprA gene found in the T1SS operon of the WCS417 genome will be instrumental in addressing this question. Another way by which plant-beneficial Pseudomonas strains may suppress local host immune responses is via T3SS-mediated injection of effector proteins. Pathogenic bacteria utilize this specialized secretion machinery to inject effector molecules into the cytoplasm of eukaryotic cells in order to suppress host immune responses and establish successful infections. In several non-pathogenic root-associated Pseudomonas bacteria, functional T3SSs or T3SS gene clusters have been found [52, 112, 137]. The T3SSs and putative effector proteins identified in WCS374 and WCS417 and the fact that we did not find T3SS-related genes in the genome of WCS358 may prove to be instrumental in investigating the role of T3SS-related effectors in suppression of root immunity.


In the past 30 years, a wealth of knowledge accumulated on the plant-beneficial functions of the PGPR WCS358, WCS374 and WCS417. The fully sequenced genomes of the WCS strains provide a genetic framework that allows for detailed analysis of the biological mechanisms of the plant-beneficial traits of the PGPR. The WCS genomes also clarified their taxonomy and revealed that very similar strains from the same species have been isolated elsewhere in the World because of their plant-beneficial properties. Hence, the sequenced WCS strains can be considered representatives of universally-present plant-beneficial Pseudomonas spp. with well-characterized functions in the stimulation of plant growth and health. Considering the increasing focus on the role of the root microbiome in plant health, functional genomics of the WCS strains will foster research toward understanding the diversity of functions of the root microbiome.


Sequencing, genome assembly, and annotation

The genomes of Pseudomonas spp. strains WCS358, WCS374, and WCS417 were sequenced following a whole-genome shotgun strategy using an Illumina HiSeq 2000 and two paired-end libraries with 500- and 2000-bp insert size at the Beijing Genomics Institute (BGI; Shenzhen, China). A hybrid genome assembly was prepared using SOAPdenovo (version 1.05) and ABySS (version 1.3.4). Multiple gaps were closed by merging overlapping contigs. The order of the remaining contigs was determined by alignment to close taxonomic relatives with MAUVE. Contigs could be combined in a single scaffold for WCS417 and WCS374. For WCS358, the order of the remaining contigs could not be determined as it differed depending on the relative it was compared to. The genome of WCS358 therefore remained distributed over eight scaffolds. Genome sequences were deposited at Genbank under the accessions CP007637 (WCS417), CP007638 (WCS374) and JMIT00000000 (WCS358). The genomes were annotated using NCBI’s Prokaryotic Genome Automatic Annotation Pipeline [144].

Bioinformatic analyses

Concatenated sequences of 16S rRNA, gyrB, rpoB and rpoD genes of WCS358, WCS374, and WCS417 were compared to the corresponding sequences of all Pseudomonas spp. type strains as described by Mulet et al. [53], and to the corresponding genes of plant-beneficial members of the P. fluorescens group of which the full genome was available [52, 145, 146]. 16S rRNA, gyrB, rpoB and rpoD genes of Pseudomonas sp. CFBP 2461 were amplified as described by Mulet et al. [53]. Throughout this study, alignments and phylogenetic trees were created with CLC main Workbench 6.7.2 (CLCbio, Aarhus, Denmark) using the neighbor joining algorithm and 1000 bootstrap replicates. Average Nucleotide Identity based on BLAST (ANIb) values were calculated using Jspecies [62]. Multiple genome alignments were performed with progressive MAUVE [147]. Genomic islands were identified using Islandviewer [65].

Mining for orthologs of genes in the genomes of WCS358, WCS374, and WCS417 was performed using reciprocal BLASTp analysis. Genes were considered orthologs when e-value <10−5. Characteristic protein domains were identified using the Pfam protein families database [148].

The genomes were screened for secondary metabolite biosynthetic clusters using AntiSMASH 2.0 [149]. Bacteriocin biosynthesis clusters were examined with BAGEL3 [102].

The amino acid composition of non-ribosomal peptide synthetase (NRPS) products was predicted using AntiSMASH, the PKS/NRPS Analysis Web-site [150] and the NRPS predictor 2 [151].

Cultivation of bacteria and media

The Pseudomonas spp. strains used in this study are listed in Table 3. All bacteria were routinely cultivated on King’s medium B agar (KBA; [152]) at 28 °C. Production of surfactants was tested in a drop-collapse assay as described by Kuiper et al. [100]. Production of siderophores by the WCS strains and their mutant derivatives was tested on CAS agar [73]. The abilities of the WCS strains and their mutant derivatives to sequester iron was tested by inoculating 10 μl of bacterial suspension (OD660 of 0.1) in 10 mM MgSO4 on KBA amended with a range of 2,2-bipyridyl concentrations (0, 400, 600, 800, 1000, 1200, 1400 or 2000 μM). Bacterial growth was assessed after 24 h of incubation at 28 °C. Cross feeding of P. fluorescens strain Pf-5 fpv mutants by the WCS strains was examined as described by Hartney et al. [32] on KBA. Briefly, bacterial cells were suspended in 10 mM MgSO4 at OD660 of 0.1 for the WCS strains and 10−2 dilutions thereof for the Pf-5 derivatives. Droplets of 10 μl were placed on the surface of KBA amended with 600 μM 2,2-bipyridyl. Receptor mutants were placed 1 cm apart from the WCS strains and growth of receptor mutants was examined after 48 h of incubation at 28 °C. Protease activity was assessed as described by Pel and co-workers [142] on KBA containing skimmed milk powder.

Table 3 Biocontrol strains and their mutant derivatives that were used in this study

Availability of supporting data

The genomes sequenced in this study are available at the website of the National Center for Biotechnology Information:

All other supporting data are included as additional files



Average nucleotide identity based on BLAST




Ferric-pyoverdine receptors


Induced systemic resistance


Jasmonic acid




Microbe-associated molecular patterns


Multi-locus sequence alignement


Nucleotide identity


Non-ribosomal peptide synthethase


O-antigen biosynthetic loci


Plant-growth promoting rhizobacteria






Salicylic acid


Systemic acquired resistance

T1SS, T2SS, etc.:

Type I secretion system, Type II secretion system, etc.


TonB dependent prtoeins


Willie Commelin Scholten


  1. 1.

    Pieterse CMJ, Van der Does D, Zamioudis C, Leon-Reyes A, Van Wees SCM. Hormonal modulation of plant immunity. Annu Rev Cell Dev Biol. 2012;28:489–521.

    CAS  PubMed  Google Scholar 

  2. 2.

    Berendsen RL, Pieterse CMJ, Bakker PAHM. The rhizosphere microbiome and plant health. Trends Plant Sci. 2012;17:478–86.

    CAS  PubMed  Google Scholar 

  3. 3.

    Mendes R, Garbeva P, Raaijmakers JM. The rhizosphere microbiome: significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol Rev. 2013;37:634–63.

    CAS  PubMed  Google Scholar 

  4. 4.

    Venturi V, Da Silva D. Incoming pathogens team up with harmless ‘resident’ bacteria. Trends Microbiol. 2012;20:160–4.

    CAS  PubMed  Google Scholar 

  5. 5.

    Bakker PAHM, Berendsen RL, Doornbos RF, Wintermans PCA, Pieterse CMJ. The rhizosphere revisited: root microbiomics. Front Plant Sci. 2013;4:165.

    PubMed Central  PubMed  Google Scholar 

  6. 6.

    Bulgarelli D, Rott M, Schlaeppi K, Ver Loren van Themaat E, Ahmadinejad N, Assenza F, et al. Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature. 2012;488:91–5.

    CAS  PubMed  Google Scholar 

  7. 7.

    Mendes R, Kruijt M, de Bruijn I, Dekkers E, van der Voort M, Schneider JHM, et al. Deciphering the rhizosphere microbiome for disease-suppressive bacteria. Science. 2011;332:1097–100.

    CAS  PubMed  Google Scholar 

  8. 8.

    Berg G, Grube M, Schloter M, Smalla K. Unraveling the plant microbiome: looking back and future perspectives. Front Microbiol. 2014;5:148.

    PubMed Central  PubMed  Google Scholar 

  9. 9.

    Lundberg DS, Lebeis SL, Paredes SH, Yourstone S, Gehring J, Malfatti S, et al. Defining the core Arabidopsis thaliana root microbiome. Nature. 2012;488:86–90.

    CAS  PubMed Central  PubMed  Google Scholar 

  10. 10.

    Kloepper JW, Ryu C-M, Zhang S. Induced systemic resistance and promotion of plant growth by Bacillus spp. Phytopathology. 2004;94:1259–66.

    CAS  PubMed  Google Scholar 

  11. 11.

    Lugtenberg B, Kamilova F. Plant-Growth-Promoting Rhizobacteria. Annu Rev Microbiol. 2009;63:541–56.

    CAS  PubMed  Google Scholar 

  12. 12.

    Mulet M, García-Valdés E, Lalucat J. Phylogenetic affiliation of Pseudomonas putida biovar A and B strains. Res Microbiol. 2013;164:351–9.

    CAS  PubMed  Google Scholar 

  13. 13.

    Pieterse CMJ, Zamioudis C, Berendsen RL, Weller DM, Van Wees SCM, Bakker PAHM. Induced systemic resistance by beneficial microbes. Ann Rev Phytopathol. 2014;52:347–75.

    CAS  Google Scholar 

  14. 14.

    Weller DM. Pseudomonas biocontrol agents of soilborne pathogens: looking back over 30 years. Phytopathology. 2007;97:250–6.

    PubMed  Google Scholar 

  15. 15.

    Mazurier S, Corberand T, Lemanceau P, Raaijmakers JM. Phenazine antibiotics produced by fluorescent pseudomonads contribute to natural soil suppressiveness to Fusarium wilt. ISME J. 2009;3:977–91.

    CAS  PubMed  Google Scholar 

  16. 16.

    Weller DM, Raaijmakers JM, Gardener B, Thomashow LS. Microbial populations responsible for specific soil suppressiveness to plant pathogens. Annu Rev Phytopathol. 2002;40:309–48.

    CAS  PubMed  Google Scholar 

  17. 17.

    Haas D, Défago G. Biological control of soil-borne pathogens by fluorescent pseudomonads. Nat Rev Microbiol. 2005;3:307–19.

    CAS  PubMed  Google Scholar 

  18. 18.

    Schippers B, Roosje GS. Hundred years of history and the future of the foundation ‘Willie Commelin Scholten Phytopathological Laboratory’. Eur J Plant Pathol. 1997;103:667–71.

    Google Scholar 

  19. 19.

    Faasse PE. In splendid isolation: A history of the Willie Commelin Scholten Phytopathology Laboratory 1894–1992. History of Science and Scholarship in the Netherlands. Amsterdam: KNAW Press; 2008.

    Google Scholar 

  20. 20.

    Geels FP, Schippers B. Selection of antagonistic fluorescent Pseudomonas spp. and their root colonization and persistence following treatment of seed potatoes. Phytopatholog Z. 1983;108:193–206.

    Google Scholar 

  21. 21.

    Lamers JG, Schippers B, Geels FP. Soil-borne diseases of wheat in the Netherlands and results of seed bacterization with Pseudomonas against Gaeumannomyces graminis var. tritici. In: Jorna ML, Slootmaker LAJ, editors. Cereal breeding related to integrated cereal production. Wageningen, the Netherlands: PUDOC; 1988. p. 134–9.

    Google Scholar 

  22. 22.

    Bakker PAHM, Lamers JG, Bakker AW, Marugg JD, Weisbeek PJ, Schippers B. The role of siderophores in potato tuber yield increase by Pseudomonas putida in a short rotation of potato. Neth J Plant Pathol. 1986;92:249–56.

    Google Scholar 

  23. 23.

    Duijff BJ, Bakker PAHM, Schippers B. Suppression of fusarium wilt of carnation by Pseudomonas putida WCS358 at different levels of disease incidence and iron availability. Biocontrol Sci Technol. 1994;4:279–88.

    Google Scholar 

  24. 24.

    De Vleesschauwer D, Djavaheri M, Bakker PAHM, Höfte M. Pseudomonas fluorescens WCS374r-induced systemic resistance in rice against Magnaporthe oryzae is based on pseudobactin-mediated priming for a salicylic acid-repressible multifaceted defense response. Plant Physiol. 2008;148:1996–2012.

    PubMed Central  PubMed  Google Scholar 

  25. 25.

    Leeman M, Den Ouden FM, Van Pelt JA, Dirkx FPM, Steijl H, Bakker PAHM, et al. Iron availability affects induction of systemic resistance to Fusarium wilt of radish by Pseudomonas fluorescens. Phytopathology. 1996;86:149–55.

    CAS  Google Scholar 

  26. 26.

    Raaijmakers JM, Leeman M, Van Oorschot MMP, Van der Sluis I, Schippers B, Bakker PAHM. Dose–response relationships in biological-control of Fusarium-wilt of radish by Pseudomonas spp. Phytopathology. 1995;85:1075–81.

    Google Scholar 

  27. 27.

    Lemanceau P, Bakker PAHM, De Kogel WJ, Alabouvette C, Schippers B. Effect of pseudobactin 358 production by Pseudomonas putida WCS358 on suppression of fusarium wilt of carnations by nonpathogenic Fusarium oxysporum Fo47. Appl Environ Microbiol. 1992;58:2978–82.

    CAS  PubMed Central  PubMed  Google Scholar 

  28. 28.

    Duijff BJ, Recorbet G, Bakker PAHM, Loper JE, Lemanceau P. Microbial antagonism at the root level is involved in suppression of Fusarium wilt by the combination of nonpathogenic Fusarium oxysporum Fo47 and Pseudomonas putida WCS358. Phytopathology. 1999;89:1073–9.

    CAS  PubMed  Google Scholar 

  29. 29.

    Lemanceau P, Expert D, Gaymard F, Bakker PAHM, Briat JF. Role of iron in plant-microbe interactions. In: Van Loon LC, editor. Plant Innate Immunity. Advances in Botanical Research. vol 51: Amsterdam: Elsevier; 2009. p. 491–549.

  30. 30.

    Cornelis P. Iron transport systems and iron homeostasis in Pseudomonas. In: Braun V, Hantke K, Cornelis P, Chakraborty R, editors. Iron uptake in bacteria with emphasis on E. coli and Pseudomonas. SpringerBriefs in Molecular Science. Netherlands: Springer; 2013. p. 67–89.

    Google Scholar 

  31. 31.

    Loper JE, Henkels MD. Utilization of Heterologous Siderophores Enhances Levels of Iron Available to Pseudomonas putida in the Rhizosphere. Appl Environ Microbiol. 1999;65:5357–63.

    CAS  PubMed Central  PubMed  Google Scholar 

  32. 32.

    Hartney S, Mazurier S, Girard M, Mehnaz S, Davis E, Gross H, et al. Ferric-pyoverdine recognition by Fpv outer membrane proteins of Pseudomonas protegens Pf-5. J Bacteriol. 2013;195:765–76.

    CAS  PubMed Central  PubMed  Google Scholar 

  33. 33.

    Bakker PAHM, Van Peer R, Schippers B. Specificity of siderophore receptors and biocontrol by Pseudomonas spp. In: Hornby D, editor. Biological control of soil-borne plant pathogens. Wallingford, UK: CAB international; 1990. p. 131–42.

    Google Scholar 

  34. 34.

    Van Peer R, Niemann GJ, Schippers B. Induced resistance and phytoalexin accumulation in biological control of Fusarium wilt of carnation by Pseudomonas strain WCS417r. Phytopathology. 1991;81:728–34.

    Google Scholar 

  35. 35.

    Alström S. Induction of disease resistance in common bean susceptible to halo blight bacterial pathogen after seed bacterization with rhizosphere Pseudomonads. J Gen Appl Microbiol. 1991;37:495–501.

    Google Scholar 

  36. 36.

    Wei G, Kloepper JW, Tuzun S. Induction of systemic resistance of cucumber to Colletrotichum orbiculare by select strains of plant-growth promoting rhizobacteria. Phytopathology. 1991;81:1508–12.

    Google Scholar 

  37. 37.

    Ton J, Davison S, Van Wees SCM, Van Loon LC, Pieterse CMJ. The Arabidopsis ISR1 locus controlling rhizobacteria-mediated induced systemic resistance is involved in ethylene signaling. Plant Physiol. 2001;125:652–61.

    CAS  PubMed Central  PubMed  Google Scholar 

  38. 38.

    Ton J, Pieterse CMJ, Van Loon LC. Identification of a locus in Arabidopsis controlling both the expression of rhizobacteria-mediated induced systemic resistance (ISR) and basal resistance against Pseudomonas syringae pv. tomato. Mol Plant-Microbe Interact. 1999;12:911–8.

    CAS  PubMed  Google Scholar 

  39. 39.

    Van der Ent S, Verhagen BWM, Van Doorn R, Bakker D, Verlaan MG, Pel MJC, et al. MYB72 is required in early signaling steps of rhizobacteria-induced systemic resistance in Arabidopsis. Plant Physiol. 2008;146:1293–304.

    PubMed Central  PubMed  Google Scholar 

  40. 40.

    Zamioudis C, Hanson J, Pieterse CMJ. β-Glucosidase BGLU42 is a MYB72-dependent key regulator of rhizobacteria-induced systemic resistance and modulates iron deficiency responses in Arabidopsis roots. New Phytol. 2014;204:368–79.

    CAS  PubMed  Google Scholar 

  41. 41.

    Pieterse CMJ, Van Wees SCM, Van Pelt JA, Knoester M, Laan R, Gerrits H, et al. A novel signaling pathway controlling induced systemic resistance in Arabidopsis. Plant Cell. 1998;10:1571–80.

    CAS  PubMed Central  PubMed  Google Scholar 

  42. 42.

    Pieterse CMJ, Van Wees SCM, Hoffland E, Van Pelt JA, Van Loon LC. Systemic resistance in Arabidopsis induced by biocontrol bacteria is independent of salicylic acid accumulation and pathogenesis-related gene expression. Plant Cell. 1996;8:1225–37.

    CAS  PubMed Central  PubMed  Google Scholar 

  43. 43.

    Fu ZQ, Dong X. Systemic acquired resistance: Turning local infection into global defense. Annu Rev Plant Biol. 2013;64:839–63.

    CAS  PubMed  Google Scholar 

  44. 44.

    Meziane H, Van der Sluis I, Van Loon LC, Höfte M, Bakker PAHM. Determinants of Pseudomonas putida WCS358 involved in inducing systemic resistance in plants. Mol Plant Pathol. 2005;6:177–85.

    PubMed  Google Scholar 

  45. 45.

    Weller DM, Mavrodi DV, van Pelt JA, Pieterse CMJ, van Loon LC, Bakker PAHM. Induced systemic resistance in Arabidopsis thaliana against Pseudomonas syringae pv. tomato by 2,4-diacetylphloroglucinol-producing Pseudomonas fluorescens. Phytopathology. 2012;102:403–12.

    CAS  PubMed  Google Scholar 

  46. 46.

    De Vleesschauwer D, Höfte M. Rhizobacteria-Induced Systemic Resistance. In: Van Loon LC, editor. Plant Innate Immunity. Advances in Botanical Research, vol 51: Amsterdam: Elsevier; 2009. p. 491–549.

  47. 47.

    Bakker PAHM, Pieterse CMJ, Van Loon LC. Induced systemic resistance by fluorescent Pseudomonas spp. Phytopathology. 2007;97:239–43.

    PubMed  Google Scholar 

  48. 48.

    Van Wees SCM, Pieterse CMJ, Trijssenaar A, Westende YAM V ‘t, Hartog F, Van Loon LC. Differential induction of systemic resistance in Arabidopsis by biocontrol bacteria. Mol Plant-Microbe Interact. 1997;10:716–24.

    PubMed  Google Scholar 

  49. 49.

    Rudrappa T, Czymmek KJ, Paré PW, Bais HP. Root-secreted malic acid recruits beneficial soil bacteria. Plant Physiol. 2008;148:1547–56.

    CAS  PubMed Central  PubMed  Google Scholar 

  50. 50.

    Millet YA, Danna CH, Clay NK, Songnuan W, Simon MD, Werck-Reichhart D, et al. Innate immune responses activated in Arabidopsis roots by microbe-associated molecular patterns. Plant Cell. 2010;22:973–90.

    CAS  PubMed Central  PubMed  Google Scholar 

  51. 51.

    Dos Santos V, Heim S, Moore E, Strätz M, Timmis K. Insights into the genomic basis of niche specificity of Pseudomonas putida KT2440. Environ Microbiol. 2004;6:1264–86.

    PubMed  Google Scholar 

  52. 52.

    Loper JE, Hassan KA, Mavrodi DV, Davis EW, Lim CK, Shaffer BT, et al. Comparative genomics of plant-associated Pseudomonas spp.: insights into diversity and inheritance of traits involved in multitrophic interactions. PLoS Genet. 2012;8:e1002784.

    CAS  PubMed Central  PubMed  Google Scholar 

  53. 53.

    Mulet M, Lalucat J, García-Valdés E. DNA sequence-based analysis of the Pseudomonas species. Environ Microbiol. 2010;12:1513–30.

    CAS  PubMed  Google Scholar 

  54. 54.

    De Weger LA, van Boxtel R, van der Burg B, Gruters R, Geels FP, Schippers B, et al. Siderophores and outer membrane proteins of antagonistic, plant-growth-stimulating, root-colonizing Pseudomonas spp. J Bacteriol. 1986;165:585–94.

    PubMed Central  PubMed  Google Scholar 

  55. 55.

    Duijff BJ, Erkelens A, Bakker PAHM, Schippers B. Influence of pH on suppression of fusarium wilt of carnation by Pseudomonas fluorescens WCS417r. J Phytopathol. 1995;143:217–22.

    Google Scholar 

  56. 56.

    Matilla MA, Pizarro-Tobias P, Roca A, Fernández M, Duque E, Molina L, et al. Complete genome of the plant growth-promoting rhizobacterium Pseudomonas putida BIRD-1. J Bacteriol. 2011;193:1290.

    CAS  PubMed Central  PubMed  Google Scholar 

  57. 57.

    Wu X, Monchy S, Taghavi S, Zhu W, Ramos J, van der Lelie D. Comparative genomics and functional analysis of niche-specific adaptation in Pseudomonas putida. FEMS Microbiol Rev. 2011;35:299–323.

    CAS  PubMed Central  PubMed  Google Scholar 

  58. 58.

    Nelson KE, Weinel C, Paulsen IT, Dodson RJ, Hilbert H, Martins dos Santos VAP, et al. Complete genome sequence and comparative analysis of the metabolically versatile Pseudomonas putida KT2440. Env Microbiol. 2002;4:799–808.

    CAS  Google Scholar 

  59. 59.

    Udaondo Z, Duque E, Fernández M, Molina L, De la Torre J, Bernal P, et al. Analysis of solvent tolerance in Pseudomonas putida DOT-T1E based on its genome sequence and a collection of mutants. FEBS Lett. 2012;586:2932–8.

    CAS  PubMed  Google Scholar 

  60. 60.

    Mathimaran N, Srivastava R, Wiemken A, Sharma A, Boller T. Genome sequences of two plant growth-promoting fluorescent Pseudomonas strains, R62 and R81. J Bacteriol. 2012;194:3272–3.

    CAS  PubMed Central  PubMed  Google Scholar 

  61. 61.

    Ye L, Hildebrand F, Dingemans J, Ballet S, Laus G, Matthijs S, et al. Draft genome sequence analysis of a Pseudomonas putida W15Oct28 Strain with antagonistic activity to Gram-positive and Pseudomonas sp. pathogens. PLoS One. 2014;9:e110038.

    PubMed Central  PubMed  Google Scholar 

  62. 62.

    Richter M, Rosselló-Móra R. Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci U S A. 2009;106:19126–31.

    CAS  PubMed Central  PubMed  Google Scholar 

  63. 63.

    Konstantinidis KT, Tiedje JM. Genomic insights that advance the species definition for prokaryotes. Proc Natl Acad Sci U S A. 2005;102:2567–72.

    CAS  PubMed Central  PubMed  Google Scholar 

  64. 64.

    Wilson M, Lindow S. Interactions between the biological control agent Pseudomonas fluorescens A506 and Erwinia amylovora in pear blossoms. Phytopathology. 1993;83:117–123.

  65. 65.

    Langille MGI, Brinkman FSL. IslandViewer: an integrated interface for computational identification and visualization of genomic islands. Bioinformatics. 2009;25:664–5.

    CAS  PubMed Central  PubMed  Google Scholar 

  66. 66.

    Wiedenbeck J, Cohan FM. Origins of bacterial diversity through horizontal genetic transfer and adaptation to new ecological niches. FEMS Microbiol Rev. 2011;35:957–76.

    CAS  PubMed  Google Scholar 

  67. 67.

    Lemanceau P, Samson R, editors. Relations entre quelques caractéristiques in vitro de 10 Pseudomonas fluorescents et leur effet sur la croissance du haricot (Phaseolus vulgaris). Les antagonismes microbiens, 24eme colloque de la SFP. Bordeaux, France: INRA; 1983.

  68. 68.

    Budzikiewicz H. Siderophores of the Pseudomonadaceae sensu stricto (fluorescent and non-fluorescent Pseudomonas spp.). Prog Ch Org Nat Prod. 2004;87:81–237.

    CAS  Google Scholar 

  69. 69.

    Fuchs R, Schafer M, Geoffroy V, Meyer J-M. Siderotyping a powerful tool for the characterization of pyoverdines. Curr Top Med Chem. 2001;1:31–57.

    CAS  PubMed  Google Scholar 

  70. 70.

    Djavaheri M, Mercado-Blanco J, Versluis C, Meyer JM, Van Loon LC, Bakker PAHM. Iron-regulated metabolites produced by Pseudomonas fluorescens WCS374r are not required for eliciting induced systemic resistance against Pseudomonas syringae pv. tomato in Arabidopsis. Microbiol Open. 2012;1:311–25.

    CAS  Google Scholar 

  71. 71.

    Raaijmakers JM, De Bruijn I, Nybroe O, Ongena M. Natural functions of lipopeptides from Bacillus and Pseudomonas: more than surfactants and antibiotics. FEMS Microbiol Rev. 2010;34:1037–62.

    CAS  PubMed  Google Scholar 

  72. 72.

    Hartney S, Mazurier S, Kidarsa T, Quecine M, Lemanceau P, Loper J. TonB-dependent outer-membrane proteins and siderophore utilization in Pseudomonas fluorescens Pf-5. Biometals. 2011;24:193–213.

    CAS  PubMed  Google Scholar 

  73. 73.

    Schwyn B, Neilands J. Universal chemical assay for the detection and determination of siderophores. Anal Biochem. 1987;160:47–56.

    CAS  PubMed  Google Scholar 

  74. 74.

    Marugg JD, van Spanje M, Hoekstra W, Schippers B, Weisbeek PJ. Isolation and analysis of genes involved in siderophore biosynthesis in plant-growth-stimulating Pseudomonas putida WCS358. J Bacteriol. 1985;164:563–70.

    CAS  PubMed Central  PubMed  Google Scholar 

  75. 75.

    Duijff BJ, Meijer JW, Bakker PAHM, Schippers B. Siderophore-mediated competition for iron and induced resistance in the suppression of fusarium wilt of carnation by fluorescent Pseudomonas spp. Neth J Plant Pathol. 1993;99:277–89.

    CAS  Google Scholar 

  76. 76.

    Mercado-Blanco J, van der Drift KMGM, Olsson PE, Thomas-Oates JE, van Loon LC, Bakker PAHM. Analysis of the pmsCEAB gene cluster involved in biosynthesis of salicylic acid and the siderophore pseudomonine in the biocontrol strain Pseudomonas fluorescens WCS374. J Bacteriol. 2001;183:1909–20.

    CAS  PubMed Central  PubMed  Google Scholar 

  77. 77.

    Ratledge C, Macham LP, Brown KA, Marshall BJ. Iron transport in Mycobacterium smegmatis: A restricted role for salicylic acid in the extracellular environment. Biochim Biophys. 1974;372:39–51.

    CAS  Google Scholar 

  78. 78.

    Chipperfield J, Ratledge C. Salicylic acid is not a bacterial siderophore: a theoretical study. Biometals. 2000;13:165–8.

    CAS  PubMed  Google Scholar 

  79. 79.

    Bakker PAHM, Ran LX, Mercado-Blanco J. Rhizobacterial salicylate production provokes headaches! Plant Soil. 2014;364:1–16.

    Google Scholar 

  80. 80.

    Matthijs S, Laus G, Meyer J-M, Abbaspour-Tehrani K, Schäfer M, Budzikiewicz H, et al. Siderophore-mediated iron acquisition in the entomopathogenic bacterium Pseudomonas entomophila L48 and its close relative Pseudomonas putida KT2440. Biometals. 2009;22:951–64.

    CAS  PubMed  Google Scholar 

  81. 81.

    Sattely ES, Walsh CT. A latent oxazoline electrophile for N − O − C bond formation in pseudomonine biosynthesis. J Am Chem Soc. 2008;130:12282–4.

    CAS  PubMed  Google Scholar 

  82. 82.

    Gross H, Loper JE. Genomics of secondary metabolite production by Pseudomonas spp. Nat Prod Rep. 2009;26:1408–46.

    CAS  PubMed  Google Scholar 

  83. 83.

    Ferguson A, Amezcua C, Halabi N, Chelliah Y, Rosen M, Ranganathan R, et al. Signal transduction pathway of TonB-dependent transporters. Proc Natl Acad Sci U S A. 2007;104:513–8.

    CAS  PubMed Central  PubMed  Google Scholar 

  84. 84.

    Mirus O, Strauss S, Nicolaisen K, von Haeseler A, Schleiff E. TonB-dependent transporters and their occurrence in cyanobacteria. BMC Biol. 2009;7:68.

    PubMed Central  PubMed  Google Scholar 

  85. 85.

    Erbs G, Newman M-A. The role of lipopolysaccharides in induction of plant defence responses. Mol Plant Pathol. 2003;4:421–5.

    CAS  PubMed  Google Scholar 

  86. 86.

    King J, Kocíncová D, Westman E, Lam J. Review: Lipopolysaccharide biosynthesis in Pseudomonas aeruginosa. Innate Immun. 2009;15:261–312.

    CAS  PubMed  Google Scholar 

  87. 87.

    Van Peer R, Schippers B. Lipopolysaccharides of plant-growth promoting Pseudomonas sp. strain WCS417r induce resistance in carnation to Fusarium wilt. Neth J Plant Pathol. 1992;98:129–39.

    Google Scholar 

  88. 88.

    Leeman M, Van Pelt JA, Den Ouden FM, Heinsbroek M, Bakker PAHM, Schippers B. Induction of systemic resistance against fusarium wilt of radish by lipopolysaccharides of Pseudomonas fluorescens. Phytopathology. 1995;85:1021–7.

    CAS  Google Scholar 

  89. 89.

    Raymond CK, Sims EH, Kas A, Spencer DH, Kutyavin TV, Ivey RG, et al. Genetic Variation at the O-Antigen Biosynthetic Locus in Pseudomonas aeruginosa. J Bacteriol. 2002;184:3614–22.

    CAS  PubMed Central  PubMed  Google Scholar 

  90. 90.

    Roberts MW, Rabinowitz JC. The effect of Escherichia coli ribosomal protein S1 on the translational specificity of bacterial ribosomes. J Biol Chem. 1989;264:2228–35.

    CAS  PubMed  Google Scholar 

  91. 91.

    De Weger LA, Jann B, Jann K, Lugtenberg B. Lipopolysaccharides of Pseudomonas spp. that stimulate plant growth: composition and use for strain identification. J Bacteriol. 1987;169:1441–6.

    PubMed Central  PubMed  Google Scholar 

  92. 92.

    Chen H, Guo Z, Liu H-w. Biosynthesis of yersiniose: attachment of the two-carbon branched-chain is catalyzed by a thiamine pyrophosphate-dependent flavoprotein. J Am Chem Soc. 1998;120:11796–7.

    CAS  Google Scholar 

  93. 93.

    Cunneen MM, Pacinelli E, Song WC, Reeves PR. Genetic analysis of the O-antigen gene clusters of Yersinia pseudotuberculosis O:6 and O:7. Glycobiology. 2011;21:1140–6.

    CAS  PubMed  Google Scholar 

  94. 94.

    Kondakova AN, Drutskaya MS, Shashkov AS, Nedospasov SA, Akimov VN, Arbatsky NP, et al. Structure of the O-polysaccharide of Pseudomonas mandelii CYar1 containing 3,6-dideoxy-4-C-[(S)-1-hydroxyethyl]-D-xylo-hexose (yersiniose A). Carbohydr Res. 2013;381:138–41.

    CAS  PubMed  Google Scholar 

  95. 95.

    Kneidinger B, O’Riordan K, Li J, Brisson J-R, Lee JC, Lam JS. Three highly conserved proteins catalyze the conversion of UDP-N-acetyl-d-glucosamine to precursors for the biosynthesis of O antigen in Pseudomonas aeruginosa O11 and capsule in Staphylococcus aureus type 5: Implications for the UDP-N-Acetyl-l-Fucosamine biosynthetic pathway. J Biol Chem. 2003;278:3615–27.

    CAS  PubMed  Google Scholar 

  96. 96.

    Tran H, Ficke A, Asiimwe T, Höfte M, Raaijmakers JM. Role of the cyclic lipopeptide massetolide A in biological control of Phytophthora infestans and in colonization of tomato plants by Pseudomonas fluorescens. New Phytol. 2007;175:731–42.

    CAS  PubMed  Google Scholar 

  97. 97.

    Audenaert K, Pattery T, Cornelis P, Höfte M. Induction of systemic resistance to Botrytis cinerea in tomato by Pseudomonas aeruginosa 7NSK2: role of salicylic scid, pyochelin, and pyocyanin. Mol Plant-Microbe Interact. 2002;15:1147–56.

    CAS  PubMed  Google Scholar 

  98. 98.

    Iavicoli A, Boutet E, Buchala A, Métraux J-P. Induced systemic resistance in Arabidopsis thaliana in response to root inoculation with Pseudomonas fluorescens CHA0. Mol Plant-Microbe Interact. 2003;16:851–8.

    CAS  PubMed  Google Scholar 

  99. 99.

    Dubern J-F, Coppoolse E, Stiekema W, Bloemberg G. Genetic and functional characterization of the gene cluster directing the biosynthesis of putisolvin I and II in Pseudomonas putida strain PCL1445. Microbiol-UK. 2008;154:2070–83.

    CAS  Google Scholar 

  100. 100.

    Kuiper I, Lagendijk EL, Pickford R, Derrick JP, Lamers GEM, Thomas-Oates JE, et al. Characterization of two Pseudomonas putida lipopeptide biosurfactants, putisolvin I and II, which inhibit biofilm formation and break down existing biofilms. Mol Microbiol. 2003;51:97–113.

    Google Scholar 

  101. 101.

    Parret A, De Mot R. Bacteria killing their own kind: novel bacteriocins of Pseudomonas and other gamma-proteobacteria. Trends Microbiol. 2002;10:107–12.

    CAS  PubMed  Google Scholar 

  102. 102.

    Van Heel A, De Jong A, Montalbán-López M, Kok J, Kuipers O. BAGEL3: Automated identification of genes encoding bacteriocins and (non-)bactericidal posttranslationally modified peptides. Nucleic Acids Res. 2013;41:53.

    Google Scholar 

  103. 103.

    Kapitein N, Mogk A. Deadly syringes: type VI secretion system activities in pathogenicity and interbacterial competition. Curr Opin Microbiol. 2013;16:52–8.

    PubMed  Google Scholar 

  104. 104.

    Filloux A. Protein Secretion Systems in Pseudomonas aeruginosa: An Essay on Diversity, Evolution, and Function. Front Microbiol. 2011;2:155.

    CAS  PubMed Central  PubMed  Google Scholar 

  105. 105.

    Guzzo J, Duong F, Wandersman C, Murgier M, Lazdunski A. The secretion genes of Pseudomonas aeruginosa alkaline protease are functionally related to those of Erwinia chrysanthemi proteases and Escherichia coli α-haemolysin. Mol Microbiol. 1991;5:447–53.

    CAS  PubMed  Google Scholar 

  106. 106.

    Létoffé S, Redeker V, Wandersman C. Isolation and characterization of an extracellular haem-binding protein from Pseudomonas aeruginosa that shares function and sequence similarities with the Serratia marcescens HasA haemophore. Mol Microbiol. 1998;28:1223–34.

    PubMed  Google Scholar 

  107. 107.

    De Groot A, Filloux A, Tommassen J. Conservation of xcp genes, involved in the two-step protein secretion process, in different Pseudomonas species and other gram-negative bacteria. Mol Gen Genet. 1991;229:278–84.

    PubMed  Google Scholar 

  108. 108.

    Putker F, Tommassen-van Boxtel R, Stork M, Rodríguez-Herva J, Koster M, Tommassen J. The type II secretion system (Xcp) of Pseudomonas putida is active and involved in the secretion of phosphatases. Environ Microbiol. 2013;15:2658–71.

    CAS  PubMed  Google Scholar 

  109. 109.

    De Vrind J, De Groot A, Brouwers GJ, Tommassen J, De Vrind-De JE. Identification of a novel Gsp-related pathway required for secretion of the manganese-oxidizing factor of Pseudomonas putida strain GB-1. Mol Microbiol. 2003;47:993–1006.

    PubMed  Google Scholar 

  110. 110.

    Filloux A, Bleves S, Van Ulsen P, Tommassen J. Protein secretion mechanisms in Pseudomonas. In: Ramos J-L, editor. Pseudomonas. New york: Kluwer; 2004. p. 749–91.

    Google Scholar 

  111. 111.

    De Groot A, Heijnen I, De Cock H, Filloux A, Tommassen J. Characterization of type IV pilus genes in plant growth-promoting Pseudomonas putida WCS358. J Bacteriol. 1994;176:642–50.

    PubMed Central  PubMed  Google Scholar 

  112. 112.

    Preston GM, Bertrand N, Rainey PB. Type III secretion in plant growth-promoting Pseudomonas fluorescens SBW25. Mol Microbiol. 2001;41:999–1014.

    CAS  PubMed  Google Scholar 

  113. 113.

    DebRoy S, Thilmony R, Kwack Y-B, Nomura K, He SY. A family of conserved bacterial effectors inhibits salicylic acid-mediated basal immunity and promotes disease necrosis in plants. Proc Natl Acad Sci U S A. 2004;101:9927–32.

    CAS  PubMed Central  PubMed  Google Scholar 

  114. 114.

    Guttman DS, Vinatzer BA, Sarkar SF, Ranall MV, Kettler G, Greenberg JT. A functional screen for the type III (Hrp) secretome of the plant pathogen Pseudomonas syringae. Science. 2002;295:1722–6.

    CAS  PubMed  Google Scholar 

  115. 115.

    Petnicki-Ocwieja T, Schneider DJ, Tam VC, Chancey ST, Shan L, Jamir Y, et al. Genome wide identification of proteins secreted by the Hrp type III protein secretion system of Pseudomonas syringae pv. tomato DC3000. Proc Natl Acad Sci U S A. 2002;99:7652–7.

    CAS  PubMed Central  PubMed  Google Scholar 

  116. 116.

    Kida Y, Taira J, Yamamoto T, Higashimoto Y, Kuwano K. EprS, an autotransporter protein of Pseudomonas aeruginosa, possessing serine protease activity induces inflammatory responses through protease-activated receptors. Cell Microbiol. 2013;1.

  117. 117.

    Campos IB, Darrieux M, Ferreira DM, Miyaji EN, Silva DA, Arêas APM, et al. Nasal immunization of mice with Lactobacillus casei expressing the Pneumococcal Surface Protein A: induction of antibodies, complement deposition and partial protection against Streptococcus pneumoniae challenge. Microbes Infect. 2008;10:481–8.

  118. 118.

    Luckett JCA, Darch O, Watters C, AbuOun M, Wright V, Paredes-Osses E et al. A Novel Virulence strategy for Pseudomonas aeruginosa mediated by an autotransporter with arginine-specific aminopeptidase activity. PLoS Pathog. 2012; 8:e1002854.

  119. 119.

    Wilhelm S, Tommassen J, Jaeger K-E. A novel lipolytic enzyme located in the outer membrane of Pseudomonas aeruginosa. J Bacteriol. 1999;181:6977–86.

    CAS  PubMed Central  PubMed  Google Scholar 

  120. 120.

    Wilhelm S, Gdynia A, Tielen P, Rosenau F, Jaeger K-E. The autotransporter esterase EstA of Pseudomonas aeruginosa is required for rhamnolipid production, cell motility, and biofilm formation. J Bacteriol. 2007;189:6695–703.

    CAS  PubMed Central  PubMed  Google Scholar 

  121. 121.

    Ruhe ZC, Low DA, Hayes CS. Bacterial contact-dependent growth inhibition. Trends Microbiol. 2013;21:230–7.

    CAS  PubMed Central  PubMed  Google Scholar 

  122. 122.

    Salacha R, Kovacic F, Brochier-Armanet C, Wilhelm S, Tommassen J, Filloux A, et al. The Pseudomonas aeruginosa patatin-like protein PlpD is the archetype of a novel Type V secretion system. Environ Microbiol. 2010;12:1498–512.

  123. 123.

    Banerji S, Flieger A. Patatin-like proteins: a new family of lipolytic enzymes present in bacteria? Microbiol-UK. 2004;150:522–5.

    CAS  Google Scholar 

  124. 124.

    Ho BT, Dong TG, Mekalanos JJ. A view to a kill: the bacterial type VI secretion system. Cell Host Microbe. 2014; 15:9-21.

  125. 125.

    Barret M, Egan F, Fargier E, Morrissey J, O’Gara F. Genomic analysis of the type VI secretion systems in Pseudomonas spp.: novel clusters and putative effectors uncovered. Microbiol-UK. 2011;157:1726–39.

    CAS  Google Scholar 

  126. 126.

    Silby M, Winstanley C, Godfrey S, Levy S, Jackson R. Pseudomonas genomes: diverse and adaptable. FEMS Microbiol Rev. 2011;35:652–80.

    CAS  PubMed  Google Scholar 

  127. 127.

    De Souza J, De Boer M, De Waard P, Van Beek T, Raaijmakers JM. Biochemical, genetic, and zoosporicidal properties of cyclic lipopeptide surfactants produced by Pseudomonas fluorescens. Appl Environ Microbiol. 2003;69:7161–72.

    PubMed Central  PubMed  Google Scholar 

  128. 128.

    Curtis TP, Sloan WT, Scannell JW. Estimating prokaryotic diversity and its limits. Proc Natl Acad Sci U S A. 2002;99:10494–9.

    CAS  PubMed Central  PubMed  Google Scholar 

  129. 129.

    Torsvik V, Øvreås L, Thingstad TF. Prokaryotic diversity - Magnitude, dynamics, and controlling factors. Science. 2002;296:1064–6.

    CAS  PubMed  Google Scholar 

  130. 130.

    Raaijmakers JM, Van der Sluis I, Koster M, Bakker PAHM, Weisbeek PJ, Schippers B. Utilization of heterologous siderophores and rhizosphere competence of fluorescent Pseudomonas spp. Can J Microbiol. 1995;41:126–35.

    CAS  Google Scholar 

  131. 131.

    Ran LC, Xiang ML, Zhou B, Bakker PAHM. Siderophores are the main determinants of fluorescent Pseudomonas strains in suppression of grey mould in Eucalyptus urophylla. Acta Phytopathol Sinica. 2005;35:6–12.

    Google Scholar 

  132. 132.

    Van Loon LC, Bakker PAHM, Van der Heijdt W, Wendehenne D, Pugin A. Early responses of tobacco suspension cells to rhizobacterial elicitors of induced systemic resistance. Mol Plant-Microbe Interact. 2008;21:1609–21.

    PubMed  Google Scholar 

  133. 133.

    Van de Mortel J, De Vos R, Dekkers E, Pineda A, Guillod L, Bouwmeester K, et al. Metabolic and transcriptomic changes induced in Arabidopsis by the rhizobacterium Pseudomonas fluorescens SS101. Plant Physiol. 2012;160:2173–88.

  134. 134.

    Bull CT, Weller DM, Thomashow LS. Relationship between root colonization and suppression of Gaeumannomyces graminis var. tritici by Pseudomonas fluorescens Strain 2–79. Phytopathology. 1991;81:954–9.

    Google Scholar 

  135. 135.

    Raaijmakers JM, Bonsall RE, Weller DM. Effect of population density of Pseudomonas fluorescens on production of 2,4-diacetylphloroglucinol in the rhizosphere of wheat. Phytopathology. 1999;89:470–5.

    CAS  PubMed  Google Scholar 

  136. 136.

    Anderson L, Stockwell VO, Loper JE. An extracellular protease of Pseudomonas fluorescens inactivates antibiotics of Pantoea agglomerans. Phytopathology. 2004;94:1228–34.

    CAS  PubMed  Google Scholar 

  137. 137.

    Rezzonico F, Binder C, Défago G, Moënne-Loccoz Y. The type III secretion system of biocontrol Pseudomonas fluorescens KD targets the phytopathogenic chromista Pythium ultimum and promotes cucumber protection. Mol Plant-Microbe Interact. 2005;18:991–1001.

    CAS  PubMed  Google Scholar 

  138. 138.

    Zamioudis C, Pieterse CMJ. Modulation of host immunity by beneficial microbes. Mol Plant-Microbe Interact. 2012;25:139–50.

    CAS  PubMed  Google Scholar 

  139. 139.

    Mavrodi DV, Joe A, Mavrodi OV, Hassan KA, Weller DM, Paulsen IT, et al. Structural and functional analysis of the type III secretion system from Pseudomonas fluorescens Q8r1-96. J Bacteriol. 2011;193:177–89.

  140. 140.

    Russell A, LeRoux M, Hathazi K, Agnello D, Ishikawa T, Wiggins P, et al. Diverse type VI secretion phospholipases are functionally plastic antibacterial effectors. Nature. 2013;496:508–12.

  141. 141.

    Hood R, Singh P, Hsu F, Güvener T, Carl M, Trinidad R, et al. A type VI secretion system of Pseudomonas aeruginosa targets a toxin to bacteria. Cell Host Microbe. 2010;7:25–37.

  142. 142.

    Pel MJC, Van Dijken AJH, Bardoel BW, Seidl MF, Van der Ent S, Van Strijp JAG, et al. Pseudomonas syringae evades host immunity by degrading flagellin monomers with alkaline protease AprA. Mol Plant-Microbe Interact. 2014;27:603–10.

  143. 143.

    Bardoel BW, Van der Ent S, Pel MJC, Tommassen J, Pieterse CMJ, Van Kessel KPM, et al. Pseudomonas evades immune recognition of flagellin in both mammals and plants. PLoS Pathog. 2011;7:e1002206.

  144. 144.

    Angiuoli S, Gussman A, Klimke W, Cochrane G, Field D, Garrity G, et al. Toward an online repository of Standard Operating Procedures (SOPs) for (meta)genomic annotation. Omics. 2008;12:137–41.

  145. 145.

    Paulsen IT, Press CM, Ravel J, Kobayashi DY, Myers GSA, Mavrodi DV, et al. Complete genome sequence of the plant commensal Pseudomonas fluorescens Pf-5. Nat Biotech. 2005;23:873–8.

  146. 146.

    Silby M, Cerdeno-Tarraga A, Vernikos G, Giddens S, Jackson R, Preston G, et al. Genomic and genetic analyses of diversity and plant interactions of Pseudomonas fluorescens. Genome Biol. 2009;10:R51.

  147. 147.

    Darling AE, Mau B, Perna NT. ProgressiveMauve: Multiple genome alignment with gene gain, loss and rearrangement. PLoS One. 2010;5:e11147.

    PubMed Central  PubMed  Google Scholar 

  148. 148.

    Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, et al. The Pfam protein families database. Nucleic Acids Res. 2012;40:290–301.

  149. 149.

    Blin K, Medema MH, Kazempour D, Fischbach MA, Breitling R, Takano E, et al. antiSMASH 2.0—a versatile platform for genome mining of secondary metabolite producers. Nucleic Acids Res. 2013;41:204–12.

  150. 150.

    Bachmann BO, Ravel J. In silico prediction of microbial polyketide and nonribosomal peptide biosynthetic pathways from DNA sequence data. Method Enzymol. 2009;458:181–217.

    CAS  Google Scholar 

  151. 151.

    Röttig M, Medema MH, Blin K, Weber T, Rausch C, Kohlbacher O. NRPSpredictor2—a web server for predicting NRPS adenylation domain specificity. Nucleic Acids Res. 2011;39:362–7.

    Google Scholar 

  152. 152.

    King EO, Ward MK, Raney DE. Two simple media for the demonstration of pyocyanin and fluorescin. J Lab Clin Med. 1954;44:301–7.

    CAS  PubMed  Google Scholar 

  153. 153.

    Verhagen BWM, Trotel-Aziz P, Couderchet M, Höfte M, Aziz A. Pseudomonas spp.-induced systemic resistance to Botrytis cinerea is associated with induction and priming of defence responses in grapevine. J Exp Bot. 2010; 61:249-60.

  154. 154.

    Nel B, Steinberg C, Labuschagne N, Viljoen A. The potential of nonpathogenic Fusarium oxysporum and other biological control organisms for suppressing fusarium wilt of banana. Plant Pathol. 2006; 55:217–23.

  155. 155.

    Bigirimana J, Höfte M. Induction of systemic resistance to Colletotrichum lindemuthianum in bean by a benzothiadiazole derivative and rhizobacteria. Phytoparasitica. 2002; 30:159-68.

  156. 156.

    Duijff BJ, Pouhair D, Olivain C, Alabouvette C, Lemanceau P. Implication of systemic induced resistance in the suppression of fusarium wilt of tomato by Pseudomonas fluorescens WCS417r and by nonpathogenic Fusarium. Eur J Plant Pathol. 1998; 104:903-10.

  157. 157.

    Weisbeek PJ, Van der Hofstad GAJM, Schippers B, Marugg JD. Genetic analysis of the iron-uptake system of two plant growth-promoting Pseudomonas strains. In: Swinburne TR, editor. Iron, siderophores and plant diseases. 2nd ed. New York: Plenum Press; 1986. p. 299–313.

    Google Scholar 

  158. 158.

    Howell C, Stipanovic R. Control of Rhizoctonia solani on cotton seedlings with Pseudomonas fluorescens and with an antibiotic produced by the bacterium. Phytopathology. 1979;69:480–2.

    CAS  Google Scholar 

  159. 159.

    Bitter W, Marugg JD, de Weger LA, Tommassen J, Weisbeek PJ. The ferric-pseudobactin receptor PupA of Pseudomonas putida WCS358: homology to TonB-dependent Escherichia coli receptors and specificity of the protein. Mol Microbiol. 1991;5:647–55.

    CAS  PubMed  Google Scholar 

  160. 160.

    Koster M, Ovaa W, Bitter W, Weisbeek PJ. Multiple outer membrane receptors for uptake of ferric pseudobactins in Pseudomonas putida WCS358. Mol Gen Genet. 1995;248:735–43.

    CAS  PubMed  Google Scholar 

  161. 161.

    Koster M, Van Klompenburg W, Bitter W, Leong J, Weisbeek PJ. Role for the outer membrane ferric siderophore receptor PupB in signal transduction across the bacterial cell envelope. EMBO J. 1994;13:2805–13.

    CAS  PubMed Central  PubMed  Google Scholar 

Download references


This work was financially supported by ERC Advanced Grant no. 269072 of the European Research Council. We would like to thank Dr. Joyce Loper for kindly providing Pseudomonas strains A506, Pf-5 and the mutant derivatives thereof, and Dr. Phillipe Lemanceau for providing Pseudomonas strain CFBP2461.

Author information



Corresponding author

Correspondence to Roeland L. Berendsen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RLB performed the experimental work. RLB, MCvV and IAS performed the bioinformatic analysis. RLB, CZ, JT, CMJP and PAHMB drafted the manuscript. All authors read and approved the final manuscript.

Additional files

Below is the link to the electronic supplementary material.

Additional file 1: Figure S1.

Venn diagrams showing the number of unique genes and genes shared between A) Pseudomonas strain WCS417 and R81 and between B) Pseudomonas strains WCS374, A506 and SS101.

Additional file 2: Table S1.

Lists of unique genes of Pseudomonas strain WCS374, not present in A506 and SS101 and vice versa and of WCS417 not present in R81 and vice versa.

Additional file 3: Figure S2.

Phylogenetic analysis of pyoverdine receptor genes in WCS358, WCS374, and WCS417. Neighbor-joining tree of all TonB-dependent proteins (TBDPs) identified in the WCS genomes and six confirmed ferric-pyoverdine receptors (FPVs) identified in the genome of P. protegens Pf-5 (in blue font). TBDPs with an N-terminal signaling domain characteristic of TonB-dependent transducers are indicated in green font; TBDPs without this domain are in red font. Bootstrap values from 1000 replicates are indicated at the nodes. Some of the WCS TBDPs clustered closely together with the FPVs of Pf-5, which indicates they have the same substrate. FpvU and FpvY, responsible for the uptake of PVD374 in Pf-5, clustered together with TBDP PC358-17710, PD374_15530 and PS417_11700 indicating that these are the receptors used for the uptake of PVD417 and PVD374 by the three strains. Previously described FPVs of WCS358 were among the 10 putative FPVs of WCS358 in the tree. PupA, which is required for the uptake of PVD358 [158, 159], does not seem to have a closely related FPV in the three other strains, which concurs with the fact that none of the other strains can be cross-fed by WCS358 on iron-limited medium. Again this receptor for PVD358 was found in the operon responsible for the biosynthesis of its peptide chain and immediately adjacent to the last NRPS gene. PupB clustered with the FpvV of Pf-5 and both TBDPs are responsible for the uptake of PVDBN7 of BN7 [160, 161]. Likewise, the TBDP RF3 clustered together with FpvW and both were demonstrated to function in the uptake of PVDB10 of B10. Although, RF2 and FpvZ clustered together, they are likely to have different substrates [32, 160]. This indicates that for definite conclusions on the substrates of the TBDPs found in this in silico analysis, further in vivo confirmation is required.

Additional file 4: Figure S3.

Identification of Pf-5 siderophore receptors involved in heterologous uptake of siderophores of WCS358, WCS374, and WCS417. Cross feeding of P. protegens Pf-5 and its siderophore receptor mutant derivatives LK032, LK148, LK150, LK151, LK153, LK154 and LK155 by siderophore donor strains WCS358, WCS374 and WCS417. WCS strains were placed in the center of the KBA plate that was supplemented with 600 μM 2,2-bipyridyl to create conditions of low iron availability. Letters indicate the TonB-dependent protein (TBDP) siderophore receptor mutants of Pf-5: A) LK032 (ΔpchA ΔpvdI), B) LK148 (ΔpchA ΔpvdI ΔfpvY), C), LK150 (ΔpchA ΔpvdI ΔfpvX), D) LK151 (ΔpchA ΔpvdI ΔfpvV), E) LK153 (ΔpchA ΔpvdI ΔfpvW), F) LK154 (ΔpchA ΔpvdI ΔfpvU), and G) LK155 (ΔpchA ΔpvdI ΔfpvZ). This set of mutants all lack biosynthesis genes for the siderophores pyochelin (ΔpchA) and pyoverdin (ΔpvdI), allowing to differentially test for the requirement of the TBDP receptors FpvU to FpvZ for the uptake of WCS siderophores.

Additional file 5: Figure S4.

Drop collapse assay for surfactant production by WCS strains. Bacterial cells of the WCS strains were suspended in a droplet of water placed on parafilm. Collapse of the droplet is an indication of surfactant production.

Additional file 6: Table S2.

Genes of WCS358, WCS374, and WCS417 putatively involved in protein secretion systems type I, II, III, IV, and VI. Gene clusters involved in protein secretion were identified in a BLASTp search with protein sequences of gene clusters of secretion systems in P. aeruginosa.

Additional file 7: Figure S5.

Protease activity of Pseudomonas strains on milk plates. Pseudomonas syringae pathovar tomato DC3000, its AprA-defective derivative [142] and the WCS strains WCS358, WCS374 and WCS417 were grown on KBA amended with skimmed milk powder. A clearing zone around the colony indicates extracellular protease activity.

Additional file 8: Table S3.

List of putative type III secreted effectors. Putative type III secreted effector genes were identified by searching for conserved Hrp (Rsp) “box motifs” in the promoter regions of putative effectors and exploring the N-terminal protein sequence of candidate effectors for features typical of type III secreted effectors (i.e. abundance of Ser and polar residues, acidic residues in the first 12 positions, and an aliphatic amino acid in position 3 or 4) [114, 115].

Rights and permissions

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Berendsen, R.L., van Verk, M.C., Stringlis, I.A. et al. Unearthing the genomes of plant-beneficial Pseudomonas model strains WCS358, WCS374 and WCS417. BMC Genomics 16, 539 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Plant growth-promoting rhizobacteria
  • Induced systemic resistance
  • Siderophores
  • Lipopolysaccharides
  • Antibiotics
  • Protein secretion system