- Open Access
Comprehensive analysis of draft genomes of two closely related pseudomonas syringae phylogroup 2b strains infecting mono- and dicotyledon host plants
© The Author(s). 2016
- Published: 28 December 2016
In recent years, the damage caused by bacterial pathogens to major crops has been increasing worldwide. Pseudomonas syringae is a widespread bacterial species that infects almost all major crops. Different P. syringae strains use a wide range of biochemical mechanisms, including phytotoxins and effectors of the type III and type IV secretion systems, which determine the specific nature of the pathogen virulence.
Strains 1845 (isolated from dicots) and 2507 (isolated from monocots) were selected for sequencing because they specialize on different groups of plants. We compared virulence factors in these and other available genomes of phylogroup 2 to find genes responsible for the specialization of bacteria. We showed that strain 1845 belongs to the clonal group that has been infecting monocots in Russia and USA for a long time (at least 50 years). Strain 1845 has relatively recently changed its host plant to dicots.
The results obtained by comparing the strain 1845 genome with the genomes of bacteria infecting monocots can help to identify the genes that define specific nature of the virulence of P. syringae strains.
- Pseudomonas syringae
- Core genome
- Virulence factors
In recent years, the damage caused by bacterial pathogens to major crops has been increasing worldwide. Pseudomonas syringae  is a widespread bacterial species that infects plants and causes many different diseases: leaf and fruit spots, cankers, blights, etc. P. syringae is a Gram-negative gamma-proteobacterium that can be isolated from more than 180 host plant species of different taxonomic groups, including nearly all major agricultural crops [2, 3]. Besides being a phytopathogen, this bacterium also occurs as an epiphyte on healthy plants or as a symbiont in phytophagous insects . It is also present in all phases of the natural water cycle: in clouds, rain water, snow, Antarctic ice, and-in association with algae-in streams and rivers. P. syringae is one of the most important objects for studying the molecular mechanisms of pathogenesis and plant response to infection .
The P. syringae species is divided into more than 50 pathovars indistinguishable by their physiological (microbiological) characteristics but infecting specific host plants or causing specific symptoms of a disease. Although P. syringae can infect a broad spectrum of plants, its individual strains are only virulent on a limited number of plant species; in the other plants, however, they either induce an immune response or cannot cause a disease at all . Based on the results of DNA–DNA hybridization, Gardan et. al. divided the P. syringae species into nine genotypes (genomospecies) . Taxonomic studies performed by multilocus sequence typing (MLST) revealed that the P. syringae species can be classified into 7 to 12 phylogenetic groups; some authors merge the P. cichorii and P. viridiflava with the P. syringae species [8–11]. The phylogenetic groups generally correspond to the genomospecies previously determined by DNA–DNA hybridization; this classification is confirmed by the genome-wide analysis of representatives of the species .
Different P. syringae strains use a wide range of biochemical mechanisms that determine their virulence to specific plants. These mechanisms include phytotoxins, ice-nucleation proteins , and effectors of the type III and type IV secretion systems, which determine the specific nature of the pathogen virulence . Bacteria of different P. syringae pathovars produce four main phytotoxins: coronatine, phaseolotoxin, syringomycin, and tabtoxin . Though these toxins are potentially important for the bacterial virulence, none of them is sufficient to cause a disease [13, 15].
To date, the genetic diversity of the P. syringae population in the Russian Federation remains poorly studied. There are several studies devoted to the comparative evaluation of 16S-23S rRNA intergenic spacer regions in P. syringae и P. fluorescens and to the rep-PCR analysis of a limited number of Pseudomonas strains isolated from cereal crops [16, 17]. Using MLST, Ignatov et. al. (in press, available at 2016) studied a collection of the P. syringae strains isolated from different crops and demonstrated that the strains of phylogroup 2 are predominant in Russia (more than 56% of all studied strains isolated from legumes, sunflower, cereal crops, brassicas, cucurbits, and grapes belong to phylogroup 2).
Phylogroup 2 of P. syringae comprises the majority of the strains isolated from the most different habitats. It comprises the pathovars P. s. pv aceris, P. s. pv aptata, P. s. pv atrofaciens, P. s. pv. avellanae, P. s. pv. coryli, P. s. pv dysoxyli, P. s. pv japonica, P. s. pv lapsa, P. s. pv papulans, P. s. pv pisi, P. s. pv solidagae, and P. s. syringae; it also comprises three previously described genetic subgroups (clades) 2a, 2b, and 2c . Subgroup 2b includes typical strains of such widespread phytopathogens as P. syringae pv. syringae, P. s. pv. aptata, and P. s. pv. atrofaciens, whereas nonpathogenic P. syringae strains, which are phenotypically similar to P. viridiflava, belong to subgroup 2c. In general, the strains of phylogroup 2 are characterized by the highest frequency of hypersensitivity reactions (HR) in tobacco plants, diseases in sprouts of wheat and sunflower, and the synthesis of syringotoxin .
To assess genetic features of the Russian P. syringae strains, we sequenced two P. syringae strains, 1845 и 2507, that infect dicots and monocots, respectively. Using phylogenetic analysis, we demonstrated that these strains are related and that they belong to phylogroup 2, clade 2b (according to Berge ) together with strains SM and B64, which infect wheat. Our aim was to compare the virulence factors in these and other available genomes of phylogroup 2 to identify genes responsible for the specialization of bacteria.
Strains and pathogenicity tests
Strains 1845 and 2507 were isolated in the Laboratory of Plant Bacterial Diseases at the Russian Research Institute for Phytopathology (Bolshie Vyazemy, Moscow region). The strains were stored at –80 °C in King’s B liquid medium  with 15% glycerol. Before use, a culture was replated on King’s B agar medium.
Strain 1845 was obtained from diseased leaves of a sunflower (cv. “Eklor”) collected in the Republic of North Ossetia-Alania in 2010. The strain was highly virulent to different dicots and inhibited the germination of their seeds. Strain 2507 was isolated from a winter wheat plant (cv. “Moscovskaya 39”) with symptoms of leaf blight collected in Krasnodar Krai in 2012. The strain was highly aggressive to monocots.
The pathogenicity and virulence of bacteria were assessed by the ability of the bacteria to induce the hypersensitivity reaction within 12 h after infiltrating leaves of tobacco (Nicotiana tabacum, Samsung cultivar) and pelargonium (Pelargonium × hortorum) with bacterial suspension. The virulence of bacteria was assessed using a number of dicot and monocot species (Additional file 1). Three methods of plants inoculation were used: (1) excised cotyledon assay-cotyledons were cut off using a scalpel dipped in a suspension of two pathogen concentrations 108 CFU/mL and 106 CFU/mL; (2) bacterial suspension (108 CFU/mL and 106 CFU/mL) was atomized into a true leaf of a plant; and (3) a hundred of seeds were soaked in a bacterial suspension (108 CFU/mL and 106 CFU/mL) for 1 h before their germination on filter paper. In all the inoculation methods, we used bacterial suspensions obtained from 48-h bacterial cultures on King’s B agar medium and diluted with Potassium Phosphate Buffer pH 7.4, 10 mM to the desirable concentrations.
Phenotypic and biochemical analysis
Morphological, physiological, and biochemical characteristics of bacterial cultures were determined by the methods for the phenotypic differentiation of the Pseudomonas genus: LOPAT (described in the manual for the identification of pathogenic bacteria ); GATTa tests (gelatin liquefaction, aesculin hydrolysis, tyrosinase activity, and L-tartrate utilization), and the analysis for the ice nucleation activity performed as described in .
DNA extraction and genome sequencing, assembly, and annotation
Bacteria were cultured on King’s B agar medium [19, 20]. Total DNA preparations were isolated from fresh cultures after 2–3 days of growth using the method of the DNA adsorption on magnetic particles (Miniprep kit, LLC Silex, Russia) according to the manufacturer’s instructions.
We sequenced 4–5 μg DNA of strains 1845 and 2507 on a 454 GS-FLX Titanium platform. The obtained readings were assembled into contigs using the GS De Novo Assembler software developed by Roche (http://www.454.com/products/analysis-software/). The sequence was annotated using the RAST Server .
Selection of strains for the joint analysis; phylogenetic analysis and calculation of the average nucleotide identity (ANI) values
To determine the P. syringae phylogroup membership, we used the phylogeny by the citrate synthase (cts) gene. It has previously been shown that the sequence of this gene is sufficiently informative to classify P. syringae strains by clades and phylogroups . For the analysis, we selected 87 P. syringae strains (Additional file 2) covering all the phylogroups . The tree was constructed by the maximum likelihood method using the RAxML package . Pseudomonas rhizosphaerae 6B4 was used as an outgroup.
To infer phylogenetic relationships between the strains of phylogroup 2, we selected 20 P. syringae strains with fully sequenced genomes, including 1845 and 2507 (Additional file 2). We then used multilocus sequence typing (MLST) to construct a phylogenetic tree based on the data on 7 household genes (RNA polymerase sigma factor – rpoD, Citrate synthase – gltA, Glyceraldehyde-3-phosphate dehydrogenase – gap1, DNA gyrase subunit – gyrB, potassium uptake protein – kup, Aconitate hydratase B – acnB, and Glucose-6-phosphate isomerase – pgi) . Strain DC3000 was used as an outgroup . As it has previously been shown, this method provides a sufficiently accurate reconstruction of the phylogenetic relationships in the Pseudomonas species .
Identification of the pan- and core genomes and unique genes; COG-analysis
The entire pan- and core genomes were identified by ortholog search using the blastp software . Genes were considered orthologous if the E-value was less than 1 × 10−10, the length of the aligned region was more than 60% of the entire gene length, and the identity of the aligned regions was more than 60%. To identify the entire pan- and core genomes, we selected 20 P. syringae strains with fully sequenced genomes, including 1845 and 2507 (Additional file 2). Among the selected strains, 8 strains infect monocots (hereinafter, group M) and 11 strains infect dicots (hereinafter, group D). Strain 1845 was used to identify pan- and core genomes but was not included into the group D because this strain, as we show in this study, has relatively recently changed the class of its host and may introduce errors into the study.
To construct the pan genome of the strains infecting monocots (Additional file 2), the genes found in the group M strains were selected from the entire pan-genome. To construct the pan-genome of the strains infecting dicots, the genes found in the group D strains were selected from the entire pan-genome (Additional file 2). Core genomes were obtained from the respective pan-genomes.
Clusters of Orthologous Groups (COG) analysis was conducted using the WebMGA service . COG class enrichment was calculated using Fisher’s exact test with the fdr correction for multiple testing (p-value < 0.05). The search for associations between the COG groups and the host plant class was performed using Fisher’s exact test with the fdr correction for multiple testing.
Finding the effectors of the type III secretion system (T3SS); T3SS, T4SS, and T6SS gene clusters; phytotoxin genes; quorum sensing genes; and mobile and CRISPR elements
We used the database on the effectors of the type III secretion system in P. syringae available on the site http://www.pseudomonas-syringae.org/. Its current version (04.10.2014) contains 128 unique effectors. The search for effectors was performed by the tblastn algorithm  with the E-value threshold of 1 × 10−5. If the alignment was incomplete, we used Baltrus’s algorithm: we searched for stop codons that are the nearest from the 5′ and 3′ termini of the aligned region. If a start codon is located between the 5′ stop codon and the start of the alignment and if this start codon belongs to the same reading frame as the aligned region and the 3′ stop codon, it is an active effector. Otherwise, the effector is truncated . The search for the clusters of the type III, IV, and VI secretion systems was conducted using the T346Hunter server .
To find the genes responsible for the phytotoxin synthesis, we used the blastn algorithm  with the E-value cutoff of 1 × 10−10 and the identity of 80%. To find the genes responsible for the quorum sensing, we used the blastn algorithm with the E-value threshold of 1 × 10−10 and the identity of 80%. Mobile elements were identified by the IS Finder service  with the E-value threshold of 1 × 10−70. Associations were found using Fisher’s test with the fdr correction for multiple testing (p-value < 0.05). CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) elements were identified using the CRISPR-finder server .
Characteristics of the strains
Main morphological, cultural, and biochemical properties of strains 1845 and 2507 of the P. syringae species are given in Additional file 3. Both strains belong to the LOPAT 1a group  and demonstrate a positive reaction for levan production (L); negative reaction for oxidase production (O); negative test result for pectinolytic activity on potato (P); negative reaction for arginine dihydrolase production (A); and negative reactions in the tobacco hypersensitivity test (T) on Nicotiana tabacum plants, ‘Samsung’ cultivar, and in the hypersensitivity test on Pelargonium × hortorum (pelargonium).
Within 12 h after the infiltration, strains 1845 and 2507 applied in the concentrations of 108 CFU/mL and 106 CFU/mL induced the hypersensitivity reaction in plants.
Moreover, they tested positive for gelatin liquefaction (G) and aesculin hydrolysis (A) and tested negative for tyrosinase activity (T) and tartrate utilization (Ta).
The results of the artificial inoculation of plants with Pseudomonas syringae strains 1845 and 2507 are given in Additional file 1.
Strains 1845 used in the concentrations of 108 CFU/mL and 106 CFU/mL induced water-soaked lesions within 48 h after inoculating the leaves of dicots by atomizer. The leaves of monocots did not show such symptoms.
When used in the concentration of 106 CFU/mL, strain 2507 did not induce severe disease symptoms in dicots but did affect monocots. At the concentration of 108 CFU/mL, the symptoms of strain 2507 infection manifested in dicots faster than in monocots and resembled the hypersensitivity reaction.
The inoculation of germinating dicots with strain 1845 used in the concentrations of 108 CFU/mL and 106 CFU/mL decreased the number of germinated seeds and inhibited the growth of seminal roots. No such symptoms were observed when inoculating the seeds of monocots. The inoculation of the plant seeds with strain 2507 slowed down the germination of monocots but did not affect the growth and development of dicots when the concentrations of 108 CFU/mL or 106 CFU/mL were used.
Sequencing, de novo assembly, and annotation of strains 1845 and 2507
Key parameters of the assembly and annotation of the genomes of P. syringae strains 1845 и 2507
Genome length (Mb)
Number of contigs
GC composition (%)
Total number of genes
rRNA (5S, 16S, 23S)
1, 1, 1
1, 1, 1
To clarify the relationship between the strains, out of the total of 87 strains, we selected 20 strains satisfying the following conditions: all selected strains belong to phylogroup 2, their genomes are fully sequenced, and the classes of their host plants are known (Additional file 2). The P. syringae pv. tomato strain DC3000 was used as an outgroup . A phylogenetic tree (Fig. 1b) was constructed for these strains based on the MLSA of seven household genes (rpoD, gltA, gap1, gyrB, kup, acnB, and pgi). The Russian strains, 1845 and 2507, belong to the same clade, namely clade 2b, as the strains B64 and SM isolated from wheat.
To confirm the correct clustering of the strains, we calculated pairwise ANI values between the genomes and performed the clustering based on these values (Fig. 1c, Additional file 4). This clustering confirmed that the Russian strains belong to clade 2b. Moreover, it showed that strains 1845, 2507, SM, and B64 cluster together.
Since strains 1845, 2507, and SM belong to the same clade (2b), we additionally compared their genomes. The contigs of genomes 1845 and 2507 were ordered according to the genome of strain SM and aligned by this genome using the progressive Mauve software . Additional file 5 shows that the genomes of the Russian strains are highly collinear with the SM genome. Long insertions and deletions in the alignment generally correspond to the regions of contig breaks. In strain SM, the chromosome region between PssSM_2181 and PssSM_2256 (57.3 Kb) is a prophage cluster. In strain 2507, the respective region changed its location and orientation (Additional file 6). Strain 1845 lacks this cluster altogether.
Another key indicator of the relationship between strains is the conservation of the hrp-cluster structure. We aligned the hrp clusters of 1845 and 2507 by strain SM (Additional file 7) . The figure shows that the structure of the cluster is conserved between the studied strains.
Analysis of the pan- and core genomes of phylogroup 2
To analyze the pan- and core genomes of phylogroup 2, we selected 20 strains with fully sequenced genomes, including two Russian strains (Additional file 2). These 20 selected strains were further divided into groups: the group M comprised 8 strains infecting monocots and the group D comprised 11 strains infecting dicots. We did not include strain 1845 into the group D because this strain has relatively recently changed the class of its host and, therefore, may contain both the genes inherent to monocot-infecting strains and the genes inherent to dicot-infecting strains.
Number of genes in the core and pan-genomes of P. syringae phylogroup 2, groups M (monocots) and D (dicots)
Second phylogroup (M + D, except strain 1845; 20 strains)
Group M (8 strains)
Group D (11 strains)
Unique genes of the group M
Unique genes of the group D
We annotated the pan-genome of all the 20 strains by COG-classes: 3633 out of the total of 6525 proteins in the pan genome were annotated. The COG-class enrichment was analyzed for the pan- and core genomes of groups M and D using the Fisher’s exact test with the correction for multiple testing. We found no significant differences between the pan- and core genomes of groups M and D (Additional files 9A and 4B). We compared the COG-class enrichment for unique genes of the pan-and core genomes of groups M and D (Additional files 9C and 4D). For the unique genes of the pan-genomes, there is a significant difference in the enrichment of three classes (D, K, and U). For the unique genes of the core genomes, there is a significant difference in the enrichment of seven classes (I, M, N, O, P, Q, and S) (Additional file 10).
COG-groups associated with the ability to infect dicots
Description of cog-group
Frequency in the group D
Frequency in the group M
Putative threonine efflux protein
Uncharacterized membrane-associated protein
Secretory systems of strains SM, 1845, and 2507 and toxins of the strains of P. syringae phylogroup 2
Number of genes of the main secretory systems in P. syringae strains 1845, 2507 and SM. Field notation: total number of the system genes in the strain/number of the system core genes in the strain/total number of the core genes in the system
T6SS was initially described in the Vibro cholera bacterium . This system was later found in almost a quarter of the species of gram-negative bacteria, mostly in known pathogens. This system is responsible for the transport of the effector proteins participating in the establishment of parasitic or symbiotic relationships between prokaryotes and eukaryotes and for the competition between prokaryotes . T6SS is generally represented by two gene clusters: T6SS-1 and T6SS-2. These two clusters were found in the Russian strains, but the core genes of these clusters are only partially present, which makes the full functioning of T6SS improbable (Table 4).
T4SS is involved in transporting proteins and DNA; particularly, it enables horizontal gene transfer between bacteria or between bacteria and plants [35, 36]. The genes of this system were found in all the three strains. SM strain comprises only 12 out of the total of 24 core components of the system. The Russian strains lost almost all T4SS. Apparently, T4SS is inessential for the virulence of these strains.
The main function of T3SS is the delivery of effector proteins into the cell body of a host [37–39]. In all the studied strains, there is a fully functional cluster T3SS-1 comprised of 15 genes and all of the 9 key genes are present. In the Russian strains, the cluster T3SS-2 only comprises 41 genes, but strain 2507 also lacks one of the 11 key genes of T3SS, whereas all of its components are present in strain 1845. The missing gene is fliP (flagellar biosynthesis protein), one of the nine mandatory membrane components . In strain SM, the cluster consists of 40 genes but it also lacks 1 of the 11 key genes of T3SS (Table 4).
T3SS effector repertoire of phylogroup 2 and the Russian strains
Mobile elements, quorum sensing, and CRISPR/Cas-system
The genomes of the second phylogroup were analyzed for the presence of mobile elements (IS). The results are shown in Fig. 3b. The total of 59 elements were found (Additional file 12): 15 of them were common, 14 elements were unique to the group M, and 40 elements were unique to the group D. Moreover, according to Fisher’s exact test, 4 of these 40 elements are significantly more common in the group D, namely ISPsy27, ISPsy28, ISPsy8, and ISPsy6. The presence of ISPs1 significantly determines the group M (p-value = 0.02). When correcting for the multiple testing using the fdr method, none of these mobile elements passes by the criterion of significance (p-value < 0.05). We found ten IS in strain 1845 and four IS in strain 2507. ISPs1 is present in both groups but none of the strains contain any of the four IS characteristic of the group D strains.
All genomes of the second clade were analyzed for the presence of genes responsible for the quorum sensing system. In gram-negative bacteria of the Pseudomonas family, this system is represented either by LuxI-LuxR-genes  or by ahlR-ahlI-locus . None of the LuxI system genes were found in the strains of phylogroup 2. The ahlR-ahlI system was found only in several of the D group strains (B728a, CC94, ISPaVe013, ISPaVe037, B301D, and CRAFRU12) that constitute one clade (Fig. 1b). These systems were not found in the Russian strains.
Using the CRISPR-finder service , we found that several regions in the Russian strains were identified as the spacer regions of the CRISPR/Cas system. However, the detailed analysis of these regions showed that these are parts of the ice nucleation protein; they also contain repetitions, which the system could mistake for a spacer region. The Cas-9 protein also wasn’t found in strains 1845 and 2507. Apparently, the absence of the CRISPR/Cas system is the feature of P. syringae strains .
Comparative genomics is often used to find the genes responsible for the virulence and specificity of phytopathogenic bacteria . We compared two previously sequenced groups of phylogroup 2 (clade 2b), infecting monocots and dicots, with two genetically similar P. syringae strains of phylogroup 2b; these two strains isolated in Russia in 2000s differ by their specialization against monocots and dicots.
Comparison of the Russian P. syringae strains
The expansion of crop infection areas and an increase in the harmfulness of P. syringae has been reported in Russia since 2004–2007 . The previous epiphytotics of basal bacteriosis and leaf spots in crops (P. syringae pv. atrofaciens and P. syringae pv. syringae) in Russia was described in the early 1970s , which is a bit later than the epiphytotics of similar diseases in 1967–1974 in USA and Canada [50–52], and seminal infection played the major role in the spread of pathogens . Taking into account the earlier emergence of strains SM (isolated in 1990s ) and B64 (isolated before 1976 [54–56]), which are the most similar to strains 1845 and 2507, it is possible that the Russian strains have common ancestors with strains B64 and SM. We can assume that the emergence of strains 1845 and 2507 in the Russian Federation is associated with the import of grain from Canada and the United States that took place from the beginning of the 1960s until the 1990s.
To the best of our knowledge, this paper is the first to provide data on sequenced P. syringae strains isolated on the territory of the Russian Federation.
The strains isolated in the Russian Federation possess the canonical type III secretion system and a very small set of T3SS effectors. Moreover, the strains contain no functional secretion systems of type IV and VI. We also found no genes of the quorum sensing system, which are crucial to inactivate bactericidal substances synthesized by the host plant. The strains contain a limited number of mobile elements (IS) and have no genes of the CRISPRs system to protect them from bacteriophages and foreign plasmids.
Comparative phylogenetic analysis
The phylogenetic analysis of strains 1845 and 2507 conducted using the sequence of the cts gene fragment  and the total of seven household genes (rpoD, gltA, gap1, gyrB, kup, acnB, and pgi) showed that Russian strains 1845 and 2507 belong to clade 2b of phylogroup 2 of the P. syringae species. Moreover, both strains are in the same clade with strains B64 and SM, which also infect monocots . The relationship of these four strains is confirmed by the conservative structure of the hrp cluster, similar ANI values, and high homology of the genomes of strains SM, B64, 1845 and 2507.
Based on the structure of the phylogenetic tree, we can conclude that strain 1845 isolated from a dicotyledonous crop is evolutionarily the youngest and that the change of its host class has occurred recently. We conducted a comparative analysis of strains infecting monocots (including SM, B64, and 2507) and dicots (except 1845) to identify the changes on the genomic level that could lead to the change of the specialization in strain 1845.
Search for the features of the M and D groups of strains
As expected, the pan- and core genomes of the strain groups M and D (pathogens to monocots and dicots) do not differ at the level of COG-classes (Additional files 9A and 4B). The difference is only observed in unique genes of the genomes of these groups (Additional files 9C and 4D). We also tried to identify associations between the host class and the COG-groups. We identified seven COG-groups (YP_235010.1, YP_234265.1, YP_234264.1, YP_234628.1, YP_235987.1, YP_237386.1, and YP_235613.1) associated with the strains that infect dicots (Table 3, Additional file 13).
Three unique genes (YP_235010.1 YP_237386.1, and YP_235613.1) belong to the LysR group of regulatory proteins. It has previously been shown that the genes of the LysR group play the regulatory role in the expression of the rovA gene, which is responsible for the virulence of Yersinia pseudotuberculosis enterobacteria . The LysR group genes, rovM, are homologous to the virulence regulators PecT/HexA of the Erwinia phytopathogens, which is another enterobacterial genus. Unique proteins of the pseudomonas of this group are most similar to the LysR proteins of gram-positive bacteria Spirosoma linguale and proteobacteria Burkholderia ubonensis and Brevundimonas diminuta.
The unique gene YP_234265.1 of the lysin exporter group LysE/YggA also participates in transporting other proteins of bacterial metabolism . Regulatory protein YP_234264.1 of the AsC/Lrp group (leucine-responsive regulatory protein/asparagine synthase C products) is one of the essential bacterial transcription regulators, which determine the metabolism intensity . Interestingly, the groups YP_234264.1 and YP_234265.1 are present in strain 1845, in all strains of the group D, and only in one strain of the group M, namely strain SM, which is the closest to strain 1845 on the phylogenetic tree.
The unique gene YP_234628.1 encodes membrane proteins of the DedA group that participate in the protection of the bacterial membrane in the human and animal pathogens Salmonella enterica and Neisseria meningitides from cationic cytolytic peptides. These proteins are also necessary for the functioning of the type III secretory system . The unique gene YP_235987.1 encodes the protein of the AraC family of transcription regulators, which control the expression of virulence genes in pathogenic bacteria . The correct functioning of the genes of the third transport system is achieved by the interaction of several regulatory systems affecting the central gene expression regulator of the AraC family . It should be noted that most of the unique genes have the closest homologues outside the Pseudomonas genus, in the genomes of pathogenic enterobacteria.
When studying the repertoire of the T3SS effectors, we showed that only 3 out of the total of 81 identified effectors are unique for the group M, whereas the group D contains 37 unique effectors (Additional file 11). However, the average number of the effectors does not differ significantly in the strains of groups M and D (18 ± 8 for the group M and 26 ± 11 for the group D). This fact indicates greater individual variety of the effectors in the strains of the group D, which might partially be explained by the host range. While for the group M hosts came from the same plant family (Poaceae), the range of host organisms of the group D consists of several families.
The analysis of statistical significance of the representation of mobile elements in the strains of groups M and D showed that the mobile element ISPs1 is inherent in the group M, whereas four elements (ISPsy27, ISPsy28, ISPsy8, and ISPsy6) are inherent in the group D. Paper  describes ISPs1 as an element unique to the strains that infect monocots . Interestingly, strain 1845 contains ISPs1 but does not contain the four mobile elements characteristic of the group D.
The genomes of Pseudomonas syringae strains 2507 (wheat) and 1845 (sunflower) isolated on the territory of the Russian Federation were determined by pyrosequencing and compared with previously published genome sequences of 18 genomes of the strains belonging to the same phylogroup and affecting dicots and monocots. We analyzed seven informative genes used in MLST genotyping of P. syringae, calculated the average nucleotide identity (ANI), studied the synteny of the hrp-gene clusters, and examined the compositions of the type III secretion system (T3SS) effectors and of the elements of insertion sequences (IS). Based on the obtained data, we found that strains 2507 and 1845 and strains SM and B64 (strains SM and B64 were isolated from wheat in the USA in 1990 and before 1976, respectively) form a subgroup that is stable among the other strains of phylogroup 2b. Within this subgroup, the strains 1845 and 2507 demonstrated the greatest similarity in the number of common unique genes. Moreover, the analysis of the genome of strain 1845 indicated the recent loss of several genetic elements (the cluster of genes responsible for the synthesis of syringolin and the prophage cluster) that are present in strains 2507, B64, and SM. We found three genes (YP_234264.1, YP_234265.1, and YP_237386.1), the acquisition of which by strain 1845 could lead to the change in its host class. The obtained results make it possible to perform a detailed study on the role of the identified genes in the specialization of P. syringae.
We thank V.A. Polityko and E.S. Pekhtereva from Russian Research Institute of Phytopathology for their contribution with isolation, identification, and virulence analysis of the bacterial strains.
This article has been published as part of BMC Genetics Vol 17 Suppl 14, 2016: Selected articles from BGRS\SB-2016: genomics. The full contents of the supplement are available online at http://bmcgenomics.biomedcentral.com/articles/supplements/volume-17-supplement-14.
Sequencing of strains 1845 and 2507 was supported by the International Science and Technology Center (project No. 3431). Bioinformatics analysis and the publication costs for this article were funded by the Russian Science Foundation (project No. 14-50-00131). Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture. USDA is an equal opportunity provider and employer.
Availability of data and materials
The genomes of P. syringae strains 1845 and 2507 are deposited at NCBI under the accession numbers LYUP00000000 and LYUO00000000, respectively.
AI, DL and SV conceived the study, supervised the experiments, performed most of the experiments; RS, GA and VG analyzed the data; RS, AI, GA and SV wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hall CJJ. Bijdragen tot de kennis der bakteriëele plantenziekten. Nederland: Proefschrift Universiteit van Amsterdam, De Gemeente Universiteit; 1902. Available from: http://library.wur.nl/WebQuery/clc/1765393.
- Hirano SS. Upper CD. Population Biology and Epidemiology of Pseudomonas. Syringae Annu Rev Phytopathol. 1990;28:155–77.View ArticleGoogle Scholar
- Young JM. Taxonomy of pseudomonas syringae. J Plant Pathol. 2010;92:S5–14. Società Italiana di Patologia Vegetale (SIPaV).Google Scholar
- Stavrinides J, McCloskey JK, Ochman H. Pea aphid as both host and vector for the phytopathogenic bacterium Pseudomonas syringae. Appl Environ Microbiol Am Soc Microbiol. 2009;75:2230–5.View ArticleGoogle Scholar
- Mansfield J, Genin S, Magori S, Citovsky V, Sriariyanum M, Ronald P, et al. Top 10 plant pathogenic bacteria in molecular plant pathology. Mol Plant Pathol. 2012;13:614–29. Wiley Online Library.View ArticlePubMedGoogle Scholar
- Lindeberg M, Cunnac S, Collmer A. The evolution of Pseudomonas syringae host specificity and type III effector repertoires. Mol Plant Pathol. 2009;10:767–75. Wiley Online Library.View ArticlePubMedGoogle Scholar
- Gardan L, Shafik H, Belouin S, Broch R, Grimont F, Grimont P. DNA relatedness among the pathovars of Pseudomonas syringae and description of Pseudomonas tremae sp. nov. and Pseudomonas cannabina sp. nov. (ex Sutic and Dowson 1959. Int J Syst Evol Microbiol Microbiol Soc. 1999;49:469–78.Google Scholar
- Sarkar SF, Guttman DS. Evolution of the core genome of Pseudomonas syringae, a highly clonal, endemic plant pathogen. Appl Environ Microbiol Am Soc Microbiol. 2004;70:1999–2012.View ArticleGoogle Scholar
- Hwang MSH, Morgan RL, Sarkar SF, Wang PW, Guttman DS. Phylogenetic characterization of virulence and resistance phenotypes of Pseudomonas syringae. Appl Environ Microbiol Am Soc Microbiol. 2005;71:5182–91.View ArticleGoogle Scholar
- Kaluzna M, Ferrante P, Sobiczewski P, Scortichini M. Characterization and genetic diversity of pseudomonas syringae from stone fruits and hazelnut using repetitive-pcr and mlst. J Plant Pathol. 2010;92:781–7. Società Italiana di Patologia Vegetale (SIPaV).Google Scholar
- Berge O, Monteil CL, Bartoli C, Chandeysson C, Guilbaud C, Sands DC, et al. A User’s Guide to a Data Base of the Diversity of Pseudomonas syringae and Its Application to Classifying Strains in This Phylogenetic Complex. PLoS One. 2014 [cited 2015 Sep 21];9. Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4153583/
- Scortichini M, Marcelletti S, Ferrante P, Firrao G. A Genomic Redefinition of Pseudomonas avellanae species. PLoS One. 2013;8:e75794. Ceнтябpь 25.View ArticlePubMedPubMed CentralGoogle Scholar
- Bender CL, Alarcón-Chaidez F, Gross DC. Pseudomonas syringae phytotoxins: mode of action, regulation, and biosynthesis by peptide and polyketide synthetases. Microbiol Mol Biol Rev Am Soc Microbiol. 1999;63:266–92.Google Scholar
- Abramovitch RB, Kim YJ, Chen S, Dickman MB, Martin GB. Pseudomonas type III effector AvrPtoB induces plant disease susceptibility by inhibition of host programmed cell death. EMBO J. 2003;22:60–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Bultreys A, Gheysen II. Biological and molecular detection of toxic lipodepsipeptide-producing pseudomonas syringae strains and PCR identification in plants. Appl Environ Microbiol Am Soc Microbiol. 1999;65:1904–9.Google Scholar
- Smirnova AV, Wang L, Rohde B, Budde I, Weingart H, Ullrich MS. Control of temperature-responsive synthesis of the phytotoxin coronatine in Pseudomonas syringae by the unconventional two-component system CorRPS. J Mol Microbiol Biotechnol. 2002;4:191–6. horizonpress.com.PubMedGoogle Scholar
- Agner G, Kaulin YA, Schagina LV, Takemoto JY, Blasko K. Effect of temperature on the formation and inactivation of syringomycin E pores in human red blood cells and bimolecular lipid membranes. Biochim Biophys Acta. 2000;1466:79–86. Elsevier.View ArticlePubMedGoogle Scholar
- Fatmi MB, Collmer A, Sante Iacobellis N, Mansfield JW. Pseudomonas syringae Pathovars and Related Pathogens. In: Fatmi ‘barek M, Collmer A, Iacobellis NS, Mansfield JW, Murillo J, Schaad NW, et al., editors. Identification, Epidemiology and Genomics. Netherlands: Springer; 2008.Google Scholar
- Lelliott RA, Billing E, Hayward AC. A determinative scheme for the fluorescent plant pathogenic pseudomonads. J Appl Bacteriol. 1966;29:470–89.View ArticlePubMedGoogle Scholar
- Braun-Kiewnick A, Sands DC. II. Gram-negative bacteria. Pseudomonas. In: N.W. Schaad, J.B. Jones, and W. Chun, editor. Laboratory Guide for Identification of Plant Pathogenic Bacteria. American Phytopathological Society Press; 2001.Google Scholar
- Brettin T, Davis JJ, Disz T, Edwards RA, Gerdes S, Olsen GJ, et al. RASTtk: a modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. Sci Rep. 2015;5:8365.View ArticlePubMedPubMed CentralGoogle Scholar
- Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, et al. Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci USA National Acad Sci. 1998;95:3140–5.View ArticleGoogle Scholar
- Buell CR, Joardar V, Lindeberg M, Selengut J, Paulsen IT, Gwinn ML, et al. The complete genome sequence of the Arabidopsis and tomato pathogen Pseudomonas syringae pv. tomato DC3000. Proc Natl Acad Sci U S A. 2003;100:10181–6.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu S, Zhu Z, Fu L, Niu B, Li W. WebMGA: a customizable web server for fast metagenomic sequence analysis. BMC Genomics. 2011;12:444.View ArticlePubMedPubMed CentralGoogle Scholar
- Baltrus DA, Nishimura MT, Romanchuk A, Chang JH, Mukhtar MS, Cherkis K, et al. Dynamic Evolution of Pathogenicity Revealed by Sequencing and Comparative Genomics of 19 Pseudomonas syringae Isolates. PLoS Pathog. 2011;7:e1002132. Июль 14.View ArticlePubMedPubMed CentralGoogle Scholar
- Martínez-García PM, Ramos C, Rodríguez-Palenzuela P. T346Hunter: a novel web-based tool for the prediction of type III, type IV and type VI secretion systems in bacterial genomes. PLoS One. 2015;10:e0119317.View ArticlePubMedPubMed CentralGoogle Scholar
- Siguier P, Perochon J, Lestrade L, Mahillon J, Chandler M. ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006;34:D32–6.View ArticlePubMedGoogle Scholar
- Grissa I, Vergnaud G, Pourcel C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007;35:W52–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Darling AE, Mau B, Perna NT. progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One. 2010;5:e11147.View ArticlePubMedPubMed CentralGoogle Scholar
- Pukatzki S, Ma AT, Sturtevant D, Krastins B, Sarracino D, Nelson WC, et al. Identification of a conserved bacterial protein secretion system in Vibrio cholerae using the Dictyostelium host model system. Proc Natl Acad Sci U S A. 2006;103:1528–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Jani AJ, Cotter PA. Type VI secretion: not just for pathogenesis anymore. Cell Host Microbe. 2010;8:2–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Juhas M, Crook DW, Hood DW. Type IV secretion systems: tools of bacterial horizontal gene transfer and virulence. Cell Microbiol. 2008;10:2377–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Voth DE, Broederdorf LJ, Graham JG. Bacterial Type IV secretion systems: versatile virulence machines. Future Microbiol. 2012;7:241–57.View ArticlePubMedPubMed CentralGoogle Scholar
- Alfano JR, Charkowski AO, Deng WL, Badel JL, Petnicki-Ocwieja T, van Dijk K, et al. The Pseudomonas syringae Hrp pathogenicity island has a tripartite mosaic structure composed of a cluster of type III secretion genes bounded by exchangeable effector and conserved effector loci that contribute to parasitic fitness and pathogenicity in plants. Proc Natl Acad Sci U S A. 2000;97:4856–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Coburn B, Sekirov I, Finlay BB. Type III secretion systems and disease. Clin Microbiol Rev. 2007;20:535–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Gazi AD, Sarris PF, Fadouloglou VE, Charova SN, Mathioudakis N, Panopoulos NJ, et al. Phylogenetic analysis of a gene cluster encoding an additional, rhizobial-like type III secretion system that is narrowly distributed among Pseudomonas syringae strains. BMC Microbiol. 2012;12:188.View ArticlePubMedPubMed CentralGoogle Scholar
- Minamino T, Macnab RM. Components of the Salmonella flagellar export apparatus and classification of export substrates. J Bacteriol. 1999;181:1388–94.PubMedPubMed CentralGoogle Scholar
- Carrión VJ, Gutiérrez-Barranquero JA, Arrebola E, Bardaji L, Codina JC, de Vicente A, et al. The Mangotoxin Biosynthetic Operon (mbo) Is Specifically Distributed within Pseudomonas syringae Genomospecies 1 and Was Acquired Only Once during Evolution. Appl Environ Microbiol. 2013;79:756–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Scholz-Schroeder BK, Hutchison ML, Grgurina I, Gross DC. The Contribution of Syringopeptin and Syringomycin to Virulence of Pseudomonas syringae pv. syringae strain B301D on the Basis of sypA and syrB1 Biosynthesis Mutant Analysis. Mol Plant Microbe Interact. 2001;14:336–48.View ArticlePubMedGoogle Scholar
- Carrión VJ, Arrebola E, Cazorla FM, Murillo J, de Vicente A. The mbo Operon Is Specific and Essential for Biosynthesis of Mangotoxin in Pseudomonas syringae. PLoS One. 2012;7:e36709.View ArticlePubMedPubMed CentralGoogle Scholar
- Miller MB, Bassler BL. Quorum sensing in bacteria. Annu Rev Microbiol. 2001;55:165–99. annualreviews.org.View ArticlePubMedGoogle Scholar
- Dumenyo CK, Mukherjee A, Chun W, Chatterjee AK. Genetic and physiological evidence for the production of N-acyl homoserine lactones by Pseudomonas syringae pv. syringae and other fluorescent plant pathogenic Pseudomonas species. Eur J Plant Pathol. 1998;104:569–82. Springer.View ArticleGoogle Scholar
- Dudnik A, Dudler R. Genomics-Based Exploration of Virulence Determinants and Host-Specific Adaptations of Pseudomonas syringae Strains Isolated from Grasses. Pathogens. 2014;3:121–48.View ArticlePubMedPubMed CentralGoogle Scholar
- Ravindran A, Jalan N, Yuan JS, Wang N, Gross DC. Comparative genomics of Pseudomonas syringae pv. syringae strains B301D and HS191 and insights into intrapathovar traits associated with plant pathogenesis. Microbiologyopen. 2015;4:553–73. Aвгуcт 1.View ArticlePubMedPubMed CentralGoogle Scholar
- Matveeva EV, Ignatov AN, Bobrova VK, Milyutina IA, Troitsky AV, Polityko VA, et al. Genetic Diversity Among Pseudomonad Strains Associated with Cereal Diseases in Russian Federation. Pseudomonas syringae Pathovars and Related Pathogens-Identification, Epidemiology and Genomics. Netherlands: Springer; 2008. p. 337–45.View ArticleGoogle Scholar
- Shneider YI, Ilyukhina MK. Brown bacteriosis of winter wheat. Zashchita Rastenii. 1975;10:22–3.Google Scholar
- Otta JD. Others. Occurrence and characteristics of isolates of Pseudomonas syringae on winter wheat. Phytopathology. 1977;67:22–6. apsnet.org.View ArticleGoogle Scholar
- Hagborg W. Others. Notes on bacterial diseases of cereals and some other crop plants. Can Plant Dis Surv. 1974;54:129–51. phytopath.ca.Google Scholar
- Wilkie JP. Basal glume rot of wheat in New Zealand. N Z J Agric Res. 1973;16:155–60. Taylor & Francis.View ArticleGoogle Scholar
- Dudnik A, Dudler R. High-Quality Draft Genome Sequence of Pseudomonas syringae pv. Syringae Strain SM, Isolated from Wheat. Genome Announc. 2013;1:e00610–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Denny TP. Phenotypic Diversity in Pseudomonas syringae pv. tomato. Microbiology. Microbiol Soc. 1988;134:1939–48.View ArticleGoogle Scholar
- Sundin GW, Demezas DH, Bender CL. Genetic and plasmid diversity within natural populations of Pseudomonas syringae with various exposures to copper and streptomycin bactericides. Appl Environ Microbiol Am Soc Microbiol. 1994;60:4421–31.Google Scholar
- Dudnik A, Dudler R. Non contiguous-finished genome sequence of Pseudomonas syringae pathovar syringae strain B64 isolated from wheat. Stand Genomic Sci. 2013;8:420–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Heroven AK, Dersch P. RovM, a novel LysR-type regulator of the virulence activator gene rovA, controls cell invasion, virulence and motility of Yersinia pseudotuberculosis. Mol Microbiol. 2006;62:1469–83. Wiley Online Library.View ArticlePubMedGoogle Scholar
- Eggeling L, Sahm H. New ubiquitous translocators: amino acid export by Corynebacterium glutamicum and Escherichia coli. Arch Microbiol . 2003;180:155–60. Springer.View ArticlePubMedGoogle Scholar
- Deng W, Wang H, Xie J. Regulatory and pathogenesis roles of Mycobacterium Lrp/AsnC family transcriptional factors. J Cell Biochem. 2011;112:2655–62. Wiley Online Library.View ArticlePubMedGoogle Scholar
- Doerrler WT, Sikdar R, Kumar S, Boughner LA. New functions for the ancient DedA membrane protein family. J Bacteriol Am Soc Microbiol. 2013;195:3–11.View ArticleGoogle Scholar
- Yang J, Tauschek M, Robins-Browne RM. Control of bacterial virulence by AraC-like regulators that respond to chemical signals. Trends Microbiol . 2011;19:128–35. Elsevier.View ArticlePubMedGoogle Scholar
- Francis MS, Wolf-Watz H, Forsberg A. Regulation of type III secretion systems. Curr Opin Microbiol . 2002;5:166–72. Elsevier.View ArticlePubMedGoogle Scholar