Pseudomonas syringae is a widespread bacterial pathogen that causes disease on a broad range of economically important plant species. Pathogenicity of P. syringae strains is dependent on the type III secretion system, which secretes a suite of up to about thirty virulence 'effector' proteins into the host cytoplasm where they subvert the eukaryotic cell physiology and disrupt host defences. P. syringae pathovar tabaci naturally causes disease on wild tobacco, the model member of the Solanaceae, a family that includes many crop species as well as on soybean.
We used the 'next-generation' Illumina sequencing platform and the Velvet short-read assembly program to generate a 145X deep 6,077,921 nucleotide draft genome sequence for P. syringae pathovar tabaci strain 11528. From our draft assembly, we predicted 5,300 potential genes encoding proteins of at least 100 amino acids long, of which 303 (5.72%) had no significant sequence similarity to those encoded by the three previously fully sequenced P. syringae genomes. Of the core set of Hrp Outer Proteins that are conserved in three previously fully sequenced P. syringae strains, most were also conserved in strain 11528, including AvrE1, HopAH2, HopAJ2, HopAK1, HopAN1, HopI, HopJ1, HopX1, HrpK1 and HrpW1. However, the hrpZ1 gene is partially deleted and hopAF1 is completely absent in 11528. The draft genome of strain 11528 also encodes close homologues of HopO1, HopT1, HopAH1, HopR1, HopV1, HopAG1, HopAS1, HopAE1, HopAR1, HopF1, and HopW1 and a degenerate HopM1'. Using a functional screen, we confirmed that hopO1, hopT1, hopAH1, hopM1', hopAE1, hopAR1, and hopAI1' are part of the virulence-associated HrpL regulon, though the hopAI1' and hopM1' sequences were degenerate with premature stop codons. We also discovered two additional HrpL-regulated effector candidates and an HrpL-regulated distant homologue of avrPto1.
The draft genome sequence facilitates the continued development of P. syringae pathovar tabaci on wild tobacco as an attractive model system for studying bacterial disease on plants. The catalogue of effectors sheds further light on the evolution of pathogenicity and host-specificity as well as providing a set of molecular tools for the study of plant defence mechanisms. We also discovered several large genomic regions in Pta 11528 that do not share detectable nucleotide sequence similarity with previously sequenced Pseudomonas genomes. These regions may include horizontally acquired islands that possibly contribute to pathogenicity or epiphytic fitness of Pta 11528.
Pseudomonas syringae is a widespread bacterial pathogen that causes disease on a broad range of economically important plant species. The species P. syringae is sub-divided into about 50 pathovars, each exhibiting characteristic disease symptoms and distinct host-specificities. P. syringae pathovar tabaci (Pta) causes wild-fire disease in soybean and tobacco plants [1, 2], characterised by chlorotic halos surrounding necrotic spots on the leaves of infected plants. Formation of halos is dependent on the beta-lactam tabtoxin, which causes ammonia accumulation in the host cell by inhibition of glutamine synthetase . However, whether tabtoxin is an essential component of the disease process is unclear [4, 5].
Pathogenicity of P. syringae strains is dependent on the type III secretion system (T3SS). The T3SS secretes a suite of virulence 'effector' proteins into the host cytoplasm where they subvert the eukaryotic cell physiology and disrupt host defences [6–14]. Mutants lacking the T3SS do not secrete effectors, and as a consequence do not infect plants or induce disease symptoms. Thus, understanding effector action is central to understanding bacterial pathogenesis. A single P. syringae strain typically encodes about 30 different effectors . However, different P. syringae strains have different complements of effector genes. The emerging view is that of a core of common effectors encoded by most strains, augmented by a variable set. Individual effectors appear to act redundantly with each other and are individually dispensable with a small or no loss to pathogen virulence . Effectors are also thought to play an important role in deTermining host range. This is most clearly true when infections are restricted by host defences. Some plants have evolved specific mechanisms to recognise certain effectors; such recognition induces strong host defences which curtail infection. For example, expression of the T3SS effector HopQ1-1 from P. syringae pathovar tomato (Pto) DC3000 was sufficient to render Pta 11528 avirulent on Nicotiana benthamiana . The opposite situation, in which acquisition of a novel effector gene confers the ability to infect new host plants, has not been demonstrated and remains speculative. However, heterologous expression of the effector gene avrPtoB conferred a plasmid-cured strain of P. syringae pathovar phaseolicola (Pph) with increased virulence . We hope that further identification and characterisation of effector repertoires of particular strains will shine new light on their roles in deTermining host range. Finally, bacterial virulence is also likely to be influenced by other non-T3SS-dependent virulence factors such as toxins which are often co-regulated with the T3SS .
Complete genome sequences are available for strains representing three P. syringae pathovars: Pto, pathovar phaseolicola (Pph) and pathovar syringae (Psy) [18–20]. Comparisons of these have led to the identification of core effector gene sets and to explain some of the differences in host-specificity between pathovars. However, these three sequenced strains are representatives of three distinct phylogroups within the species P. syringae, and as such are phylogenetically quite distant [21, 22]. According to DNA-DNA hybridisation studies and ribotyping , P. syringae can be divided into 9 discrete genomospecies. Representative strains of Psy, Pph and Pto fell into genomospecies one, two and three respectively . Recently, a strain of pathovar oryzae (genomospecies four) was sequenced . A draft genome sequence was also published for Pto T1 , a strain closely related to Pto DC3000 but restricted to tomato hosts, whereas Pto DC3000 is able to cause disease on Arabidopsis. In the current study, we explore genetic differences at an inTermediate phylogenetic resolution; that is, we compared the genome sequences of Pta 11528 to that of P. phaseolicola (Pph) 1448A, which resides within the same phylogroup but possesses a distinct host range and causes different disease symptoms.
Pto DC3000 was the first plant-pathogenic pseudomonad to have its genome sequenced, helping to establish the Arabidopsis-Pto system as the primary model for plant-microbe interactions. However, Arabidopsis is not a natural host of Pto, and it is important to develop alternative systems given the genetic variability of P. syringae strains, particularly in regard to effectors. We work on the interaction between Pta and the wild tobacco plant N. benthamiana, which offers certain advantages over Arabidopsis. Firstly, N. benthamiana is an important model for the Solanaceae, which includes many important crop species. The Pta-N. benthamiana interaction is a natural pathosystem. Lastly, N. benthamiana is an important model plant that is more amenable to biochemistry-based approaches and facile manipulation of gene expression such as virus-induced gene silencing (VIGS). Thus N. benthamiana provides experimental options for understanding plant-bacterial interactions. Strains of Pta can cause disease on N. benthamiana, but relatively few genetic sequence data are available for this pathovar.
In this study we generated a draft complete genome sequence of Pta 11528 and used a functional screen for HrpL-dependent genes to infer its repertoire of T3SS effectors and associated Hrp Outer Proteins (Hops), which differs significantly from that of its closest relative whose complete genome has previously been published (Pph 1448A). Pta 11528 does not encode functional homologues of HopAF1 or HrpZ1. This was surprising since HopAF1 was conserved in the three previously sequenced pathovars [18–20]. HrpZ1 is conserved in most strains of P. syringae that have been investigated, albeit with differences in amino acid sequence . However, Pta strain 6605 and several other isolates from Japan, were previously shown to carry a major deletion leading to truncated HrpZ protein product . Pta 11528 encodes several novel potential T3SS effectors for which no close orthologues have been reported. We also discovered several large genomic regions in Pta 11528 that do not share detectable nucleotide sequence similarity with previously sequenced Pseudomonas genomes. These regions may be horizontally acquired islands that possibly contribute to pathogenicity or epiphytic fitness of Pta 11528.
Results and discussion
Sequencing and assembly of the Pta 11528 genome
The Illumina sequencing platform provides a cost-effective and rapid means to generate nucleotide sequence data [27–29]. Although this method generates very short sequence reads, several recent studies have demonstrated that it is possible to assemble these short reads into good quality draft genome sequences [30–41].
We generated 12,096,631 pairs of 36-nucleotide reads for a total of 870,957,432 nucleotides. This represents approximately 145X depth of coverage assuming a genome size of six megabases. We used Velvet 0.7.18  to assemble the reads de novo. Our resulting assembly had 71 supercontigs of mean length 85,604 nucleotides, an N50 number of eight, and N50 length of 317,167 nucleotides; that is, the eight longest supercontigs were all at least 317,167 nucleotides long and together covered more than 50% of the predicted genome size of six megabases. The largest supercontig was 606,547 nucleotides long. The total length of the 71 assembled supercontigs was 6,077,921 nucleotides. The G+C content of the assembly was 57.96%, similar to that of the previously sequenced P. syringae genomes (Table 1). The sequence data from this project have been deposited at DDBJ/EMBL/GenBank under the accession ACHU00000000. The version described in this paper is the first version, ACHU01000000. The data can also be accessed from the authors' website http://tinyurl.com/Pta11528-data and as Additional files submitted with this manuscript. In addition, an interactive genome browser is available from the authors' website http://tinyurl.com/Pta11528-browser.
Comparison of Pta 11528 genome properties with those of previously sequenced P. syringae genomes [18–20, 83–85], [86-93].
RefSeq accession number
G+C content (%)
Nucleic acid sequence identity to P. syringae pv tabaci 11528 draft assembly (%)
P. syringae pv. tabaci 11528 draft genome assembly
P. syringae pv. phaseolicola 1448A, chromosome
P. syringae pv. phaseolicola 1448A large plasmid
P. syringae pv. tomato str. DC3000 plasmid pDC3000B
P. syringae pv. syringae B728a, chromosome
P. syringae pv. tomato str. DC3000 plasmid pDC3000A
P. syringae pv. phaseolicola 1448A small plasmid
P. syringae pv. tomato str. DC3000
Percentage identities were calculated over the alignable portion of the genomes using MUMMER .
We aligned the 71 Pta supercontigs against published complete Pseudomonas genome sequences using MUMMER . The Pta 11528 genome was most similar to that of Pph 1448A, with 97.02% nucleotide sequence identity over the alignable portions. The next most similar genome was that of Pto DC3000, with less than 90% identity (Table 1). This pattern of sequence similarity is consistent with phylogenetic studies that placed strains of Pta in the same phylogroup as Pph and revealed a relatively distant relationship to Pto [21, 22].
Comparison of the protein complement of Pta 11528 versus Pph 1448A and other pseudomonads
Using the FgenesB annotation pipeline http://www.softberry.com, we identified 6,057 potential protein-coding genes, of which 5,300 were predicted to encode proteins of at least 100 amino acids long. Of 5,300 predicted Pta 11528 proteins, 575 (10.8%) had no detectable homology with Pph 1448A proteins (based on our criterion of an E-value less than 1e-10 using BLASTP). Of these 575 sequences, 303 had no detectable homologues in Psy B728a nor Pto DC3000. These 303 Pta-specific sequences had a median length of 198 amino acids whereas the median length of the 5,300 sequences was 216 amino acids. Automated gene prediction is not infallible and inevitably a subset of the predictions will be incorrect. The reliability of gene predictions is poorer for short sequences than for longer ones. This slight enrichment for very short sequences among the Pta-specific gene predictions might be explained by the inclusion of some open reading frames that are not functional genes among those 303. However, many of the predicted proteins showed significant similarity to other proteins in the NCBI NR databases (See Additional file 1: Table S1), confirming that these are likely to be genuine conserved genes.
Conservation of the T3SS apParatus and T3SS-dependent effectors
The Hop Database (HopDB, http://www.pseudomonas-syringae.org) provides a catalogue of confirmed and predicted hop genes . Figure 1 lists the hop genes in HopDB for the three previously fully sequenced P. syringae genomes. A 'core' set of hop genes are conserved in all three previously sequenced pathovars: avrE1, hopAF1, hopAH2, hopAJ2, hopAK1, hopAN1, hopI1, hopJ1, hopX1, hrpK1, hrpW1 and hrpZ1. In addition to this core set, each genome contains additional hop genes that are found in only a subset of the sequenced strains. The Pta 11528 homologues of hop genes are listed in Table 2. Figure 1 also indicates those hop genes for which a close homologue was found to be encoded in Pta 11528.
Homologues of known hop genes in Pta 11528. Homologues were detected by searching the Pta 11528 FgenesB-predicted protein sequences against HopDB http://www.pseudomonas-syringae.org using BLASTP
Gene in Pta 11528 genome (location)
Hrp-box HMM score (bioinformatic evidence)
HrpL-dependent (functional screen)
Homologue in Pph 1448A
(PSPPH is a truncated HopW1 homologue)
The locations of Pph 1448A homologous genes is indicated, including an indication of whether they are located on the chromosome or on the large plasmid. Also indicated is whether each gene appeared in the functional screen for Hrp-dependent transcription.
In sequenced strains of P. syringae, the gene cluster encoding the T3SS apParatus is flanked by collections of effector genes Termed the exchangeable effector locus (EEL) and the conserved effector locus (CEL). Together, these three genetic components comprise the Hrp pathogenicity island . A core set of hop genes is located in the Hrp pathogenicity island , which is highly conserved between Pta 11528 and Pph 1448A (Figure 2), except that in Pta 11528 there is a deletion in hrpZ1 and an insertion in the hrpV-hrcU intergenic region. The core hop genes avrE1, hopAH2, hopAJ2, hopAK1, hopAN1, hopI1, hopJ1, hopX1 and hrpK1 are conserved in Pta 11528 and encode intact full-length proteins. Pta 11528 encodes a full-length HrpW1 protein, albeit with insertions of 69 and 12 nucleotides relative to the Pph 1448A sequence. However, there is a large deletion in hrpZ1 that likely renders it non-functional and hopAF1 is completely absent.
Besides the core conserved hop genes, the Pta 11528 genome assembly contains full-length orthologues of hopR1, hopAS1, hopAE1 and hopV1, which are also found in Pph 1448A but are absent from Psy B728a and/or Pto DC3000.
The hrpZ1 gene encodes a harpin, which is not classified as a type III effector because it is not injected directly into host cells. Harpins are characteristically acidic, heat-stable and enriched for glycine, lack cysteine residues  and can induce defences in both host and non-host plants [45, 46]. HrpZ1 forms pores in the host membrane  suggesting a role in translocation of effectors across the host membrane. It also shows sequence-specific protein binding activity . HrpZ1 can induce defences in both host and non-host plants and tobacco has been extensively used as the non-host plant species [45, 46]. The inactivation of hrpZ1 in Pta 11528 and other strains of Pta  may be an adaptive strategy and have been an important process in the stepwise progression towards compatibility, allowing Pta 11528 to avoid detection by the tobacco host plant. This is reminiscent of the "black holes" and other processes that inactivate genes whose expressed products are detrimental to a pathogenic lifestyle [49, 50]. One excellent example is the inactivation of cadA in genomes of Shigella species as compared to the genome of their closely related but non-pathogenic Escherichia coli strain [51, 52].
Pta 11528 contains highly conserved homologues of hopAB2, hopW, hopO1-1, hopT1-1, hopAG1, hopAH1, hopF1 and hopAR1, which are absent in Pph 1448A. Although absent from the Pph 1448A genome, hopAR1 and hopF1 have been identified in other strains of Pph [53–57]. In Pph 1302A, hopAR1 is located on the pathogenicity island PPHGI-1, though its genomic location varies between strains [56, 57]. PPHGI-1 is absent from the Pph 1448A genome . The Pta 11528 genome (supercontig 1087) possesses a region of similarity to PTPHGI-1, but which contains a substantial number of insertions and deletions (Additional file 2: Figure S1). The Pta 11528 hopAR1 homologue (C1E_2036) is not located in the PPHGI-1 region; it falls on supercontig 672 about two kilobases upstream of a gene encoding a protein (C1E_2039) sharing 43% amino acid identity with Pto DC3000 avrPto1. In contrast to AvrPto1 from Pto DC3000, the AvrPto1 homologue (C1E_2039) from Pta 11528 is not recognised by the plant Pto/Prf system (S. Gimenez Ibanez and J. Rathjen, manuscript in preParation).
The homologues of hopAG1, hopAH1 and the degenerate hopAI1' are found within a region of the Pta 11528 genome that shares synteny with the chromosome of Psy B728a. This region is also conserved in Pto DC3000A, albeit with several deletions and insertions, suggesting that these effector genes are ancestral to the divergence of the pathovars and have been lost in Pph 1448A rather than having been laterally transferred laterally between Pta 11528 and Psy B728a. In Pto DC3000, hopAG1 (PSPTO_0901) has been disrupted by an insertion sequence (IS) element. This is consistent with a model of lineage-specific loss of certain ancestral effectors.
In Pto DC3000, hopO1-1 and hopT1-1 are located on the large plasmid pDC3000A; homologues of these effector-encoding genes are not found in Pph 1448A. The Pta 11528 genome contains a three kilobase region of homology to pDC3000 comprising homologues of these two effector genes and a homologue of the ShcO1 chaperone-encoding gene. These three genes are situated in a large (at least 50 kilobase) region of the Pta 11528 genome that has only limited sequence similarity with Pph 1448A. Two tRNA genes (tRNA-Pro and tRNA-Lys) are located at the boundary of this region (Figure 3), which would be consistent with this comprising a mobile island.
In plasmid pMA4326B from P. syringae pathovar maculicola (Pma), the hopW1 effector gene is immediately adjacent to a three-gene cassette comprising a resolvase, an integrase and exeA. This cassette is also found in plasmids and chromosomes of several human-pathogenic Gram-negative bacteria . We found a homologue of this cassette along with a hopW1 homologue on supercontig 955 of the Pta 11528 genome assembly. Stavrinides and Guttman  proposed that the boundaries of the cassette lay upstream of the resolvase and upstream of hopW1. The presence of this four-gene unit in a completely different location in Pta 11528 is indeed consistent with the hypothesis that it represents a discrete mobile unit.
Several hop genes are located on the large plasmid of Pph 1448A. We found no homologues of these genes in Pta 11528, suggesting that the plasmid is not present in Pta 11528. Consistent with this, only a small proportion of the plasmid was alignable to our 36-nucleotide Illumina sequence reads (Figure 4). This reveals that a large component of the pathogen's effector arsenal is deTermined by its complement of plasmids. However, simple loss or gain of a plasmid does not explain all of the differences in effector complement since Pta 11528 lacks homologues of several Pph 1448A chromosomally-located effector-encoding hop genes hopG1, hopAF1, avrB4, hopF3 and hopAT1 as well as the non-effector hopAJ1. It also lacks homologues of the Pph 1448A degenerate effector gene hopAB3'.
The regions of the Pph 1448A large plasmid that are apparently conserved in Pta 11528 include genes encoding the conjugal transfer system, suggesting the presence of one or more plasmids in this strain. We found an open reading frame (C1E_3950, located on supercontig 955 coordinates 59126-60394) encoding a protein with about 97% sequence identity to the RepA proteins characteristically encoded on pT23A-family plasmids (e.g. AAW01447; reviewed in ), suggesting that this 236 kilobase supercontig might represent a plasmid.
A functional screen for HrpL-regulated genes
We used a previously described functional screen  to complement our bioinformatics-based searches for type III effectors of Pta 11528. Our functional screen was based on two steps. The first step was employed to identify genes whose expression was regulated by the T3SS alternative sigma factor, HrpL. The second step was used to identify the subset of HrpL-regulated genes that encoded effectors. For Pta 11528, we employed only the first step to identify candidate effector genes based on induced expression by HrpL. A library was constructed from Pta 11528 into a broad-host range vector carrying a promoter-less GFP and mobilized into Pto lacking its endogenous hrpL but conditionally complemented with an arabinose-inducible hrpL. We used a fluorescence activated cell sorter (FACS) to select clones that carried HrpL-inducible promoters based on expression of GFP after growth in arabinose. Clones were sequenced and sequences were assembled. Clones representative of assembled supercontigs were verified again for HrpL regulation using FACS. Among the genes whose expression was confirmed to be HrpL-dependent were those encoding effectors hopAE1, hopI1, hopAR1, the avrPto1-like gene, hopF1, hopT1-1, hopO1-1, avrE1, hopX1, and the degenerate hopM1' and hopAI1' as well as known T3SS-associated genes hrpH (ORF1 of the CEL; ) and hrpW1. Interestingly, the screen also confirmed HrpL-dependent regulation of genes encoding a major facilitator superfamily (MFS) permease and a putative peptidase (Table 3).
Pta 11528 genes confirmed by the functional screen to be under the transcriptional control of HrpL
Gene in Pta 11528 gnome
Hrp-box HMM score (bioinformatic evidence)
Orthologue in Pph 1448A
Major facilitator superfamily permease
Putative M20 peptidase
Other differences in predicted proteomes of P. syringae strains
Host range and pathogenicity are likely to be further influenced by genes other than those associated with type III secretion. Virulence deTerminants in P. syringae include toxins as well as epiphytic fitness; that is, the ability to acquire nutrients and survive on the leaf surface . Epiphytic fitness depends on quorum-sensing , chemotaxis , osmo-protection, extracellular polysaccharides, glycosylation of extracellular structures  iron uptake  and the ability to form biofilms. Cell-wall-degrading hydrolytic enzymes play a role in virulence in at least some plant-pathogenic pseudomonads ). Secretion systems (including type I, type II, type IV, type V, type VI and twin arginine transporter) may also contribute to both virulence and epiphytic fitness , whilst multidrug efflux pumps may confer resistance to plant-derived antimicrobials .
To identify differences between Pta 11528 and the previously sequenced Pph 1448A, Psy B728a and Pto DC3000 with respect to their repertoires of virulence factors, we performed BLASTP searches between the predicted proteomes. We found no significant differences in the repertoires of secretion systems between the proteomes. However, we found that Pta 11528 lacks homologues of several Pph 1448A polysaccharide modifying enzymes (glycosyl transferase PSPPH_0951, polysaccharide lyase PSPPH_1510, glycosyl transferase PSPPH_3642). Conversely, Pta 11528 encodes two glycosyl transferases (C1E_0355 and C1E_0361) and a thermostable glycosylase (C1E_4802) that do not have homologues in any of the three fully sequenced P. syringae genomes. This may imply differences in the extracellular polysaccharide profiles. In contrast to Pph 1448A, Pta 11528 lacks homologues of RhsA insecticidal toxins (PSPPH_4042 and PSPPH_4043). However, a tabtoxin biosynthesis gene cluster is found in the Pta 11528 genome and shows a high degree of conservation with the previously sequenced Pta BR2 tabtoxin biosynthesis cluster .
Pta 11528 encodes several enzymes that do not have homologues in any of the three fully sequenced P. syringae genomes (Table 4), including a predicted gluconolactonase (C1E_2553), a predicted dienelactone hydrolase (C1E_2589), a predicted nitroreductase (C1E_6026), and a sulphotransferase (C1E_6026). C1E_0903 shares 71.4% amino acid sequence identity with a predicted epoxide hydrolase (YP_745600.1) from Granulibacter bethesdensis CGDNIH1  and has a significant match to the epoxide hydrolase N-Terminal domain in the Pfam database (PF06441) [71, 72]. Epoxide hydrolases are found in P. aeruginosa and P. fluorescens PfO-1, but not in any other pseudomonads. It is possible that this gene product has a function in detoxification of host-derived secondary metabolites.
Proteins encoded by the draft Pta 11528 genome that have no detectable homologues on three previously fully sequenced P. syringae genomes.
Length (amino acids)
Predicted function (FgenesB automated annotation)
Glycosyltransferase involved in cell wall biogenesis
Phosphoadenosine phosphosulfate reductase
Serine/threonine protein kinase
TetR family ranscriptional regulator. 49% amino acid sequence identity to G. bethesdensis GbCGDNIH1_1777 .
Hydrolases or acyltransferases (alpha/beta hydrolase superfamily)
Tabtoxin biosynthesis enzyme, TblA
Histone acetyltransferase HPA2
Amine oxidase, flavin-containing
Rhs family protein
Xenobiotic response element family of transcriptional regulator. 37% amino acid sequence identity to Xylella fastidiosa PD0954 .
Histone acetyltransferase HPA2
Similar to Mucin-1 precursor (MUC-1)
Site-specific recombinases, DNA invertase Pin homologs
LacI family transcriptional regulator. 57% amino acid sequence identity to Rhizobium leguminosarum plasmid-encoded pRL1201 .
Tfp pilus assembly protein, major pilin PilA. 42% amino acid sequence identity to P. aeruginosa UniProt:P17838 .
Histone acetyltransferase HPA2
NTPase (NACHT family)
Permeases of the major facilitator superfamily
Short-chain dehydrogenase/reductase SDR
ASPIC/UnbV domain-containing protein
Xenobiotic response element family of transcriptional regulator. 38% amino acid sequence identity to P. aeruginosa PACL_0260 .
Cro/CI family transcriptional regulator. 36% amino acid sequence identity to Pto DC3000 PSPTO_2855 .
TetR family transcriptional regulator. 56% amino acid sequence identity to Ralstonia solanacearum RSc0820 .
Protein-coding genes were predicted and automatically annotated using the FgenesB pipeline http://www.softberry.com. Only those proteins are shown for which a predicted function could be proposed.
Pta protein C1E_6026 has a significant match to the sulphotransferase domain (Pfam:PF00685). Examples of this protein domain have not been found in other pseudomonads except for P. fluorescens PfO-1. Sulphotransferase proteins include flavonyl 3-sulphotransferase, aryl sulphotransferase, alcohol sulphotransferase, estrogen sulphotransferase and phenol-sulphating phenol sulphotransferase. These enzymes are responsible for the transfer of sulphate groups to specific compounds. The sulphotransferase gene (C1E_6026, 82% amino acid identity to P. fluorescens Pfl01_0157) overlaps a two kilobase Pta 11528-specific genomic island that also encodes a phage tail collar-protein encoding gene (C1E_5461, 61% amino acid identity to P. fluorescens Pfl01_0155) and an acetyltransferase (C1E_5459, 76% amino acid identity to P. fluorescens Pfl01_0148). We speculate that this region has been horizontally acquired in the Pta 11528 lineage via a bacteriophage.
An 80 kilobase region of Pta 11528 supercontig 684 contains two open reading frames (ORFs) (C1E_2584 and C1E_2585) whose respective predicted protein products show 48 and 55% amino acid identity to the C- and N-Termini of a P. putida methyl-accepting chemotaxis protein (MCP) (PP_2643) and little similarity to any P. syringae protein. Since the N- and C-Termini are divided into seParate reading frames, this probably represents a degenerate pseudogene. Immediately downstream of these ORFs is a gene (C1E_2583) that specifies a MCP showing greatest sequence identity (70%) to PP_2643 from P. putida, whilst sharing only 65% identity to its closest homologue in P. syringae (PSPPH_4743). This region also encodes another MCP (C1E_2587) that shares only 50% amino acid identity with any previously sequenced P. syringae homologue. It remains to be tested whether these MCPs play a role in pathogenesis and/or epiphytic fitness.
Transcriptional regulators are not normally considered to be virulence factors. However, expression of virulence factors may be coordinated by and dependent on regulators. Moreover, heterologous expression of the RscS regulator was recently shown to be sufficient to transform a fish symbiont into a squid symbiont . Pta 11528 encodes several predicted transcriptional regulators that are not found in Pto DC3000, Psy B728a and Pph 1448A. These include two predicted TetR-like proteins (C1E_0901 and C1E_6027), two predicted xenobiotic response element proteins (C1E_2056 and C1E_2563), a LacI-like protein (C1E_2286), a Cro/CI family protein (C1E_2570) and an IclR family protein (C1E_5715).
Pta 11528 encodes a novel pilin (C1E_2329) not found in previously sequenced P. syringae strains but sharing significant sequence similarity with a type IV pilin from P. aeruginosa . Pilin is the major protein component of the type IV pili, which have functions in forming micro-colonies and biofilms, host-cell adhesion, signalling, phage-attachment, DNA uptake and surface motility, and have been implicated as virulence factors in animal-pathogenic bacteria . The precise function of the C1E_2329 pilin is unknown but it may be involved in epiphytic fitness or plant-pathogenesis or could even be involved in an interaction with an insect vector.
Pta-specific genomic islands
We identified 102 genomic regions of at least one kilobase in length which gave no BLASTN matches against previously sequenced Pseudomonas genomes (Additional file 3: Table S2). Ten of the Pta 11528-specific regions are longer than 10 kilobases, the longest being 37.7, 21.8, 18.7, 17.9 and 16.6 kilobases. The 16.6 kilobase region corresponds to the tabtoxin biosynthesis gene cluster . These regions will be good candidates for further study of the genetic basis for association of Pta with the tobacco host. For example, several of the islands encode MFS transporters and other efflux proteins that might be involved in protection from plant-derived antimicrobials (Additional file 3: Table S2).
We have generated a draft complete genome sequence for the Pta 11528 a pathogen that naturally causes disease in wild tobacco, an important model system for studying plant disease and immunity. From this sequence, combined with a functional screen, we were able to deduce the pathogen's repertoire of T3SS-associated Hop proteins. This has revealed some important differences between Pta and other pathovars with respect to the arsenal of T3SS effectors at their disposal for use against the host plant. We also revealed more than a hundred Pta-specific genomic regions that are not conserved in any other sequenced P. syringae, providing many potential leads for the further study of the Pta-tobacco disease system.
Solexa sequence data were assembled using Velvet 0.7.18 . We used Softberry's FgenesB pipeline http://www.softberry.com to predict genes encoding rRNAs, tDNAs and proteins. Annotation of protein-coding genes by FgenesB was based on the NCBI NR Proteins database.
Prediction of HrpL-binding sites (Hrp boxes)
We built a profile hidden Markov model (HMM) based on a multiple sequence alignment of 26 known Hrp boxes from Pto DC3000 using hmmb from the HMMER 1.8.5 package http://hmmer.janelia.org. DNA sequence was scanned against this profile-HMM using hmmls from HMMER 1.8.5 with a bit-score cut-off of 12.0.
Functional screen for candidate type III effectors
Library preParation and the Flow cytometric-based screen for HrpL-induced genes of Pta 11528 were done according to .
Visualisation of data
We generated graphical views of genome alignments using CGView . To visualise the annotation draft genome assembly of Pta11528, we used the 'gbrowse' Generic Genome Browser .
Library preParation for Illumina sequencing
DNA was prepared from bacteria grown in L-medium using the Puregene Genomic DNA Purification Kit (Gentra Systems, Inc., Minneapolis, USA) according to manufacturer's instructions. A library for Illumina Paired-End sequencing was prepared from 5 mg DNA using a Paired-End DNA Sample Prep Kit (Pe-102-1001, Illumina, Inc., Cambridge, UK). DNA was fragmented by nebulisation for 6 min at a pressure of 32 psi. For end-repair and phosphorylation, sheared DNA was purified using QIAquick Nucleotide Removal Kit (Quiagen, Crawley, UK). The end repaired DNA was A-tailed and ada Ptors were ligated according to manufacturer's instructions.
Size fractionation and purification of ligation products was performed using a 5% polyacrylamide gel run in TBE at 180V for 120 min. Gel slices were cut containing DNA in the 500 to 10 bp range. DNA was than extracted using 0.3 M sodium acetate and 2 mM EDTA [pH 8.0] followed by ethanol precipitation. Using 18 PCR cycles with primer PE1.0 and PE2.0 supplied by Illumina, 5' ada Ptor extension and enrichment of the library was performed. The library was finally purified using a QIAquick PCR Purification Kit and adjusted to a concentration of 10 nM in 0.1% Tween. The stock was kept at -20°C until used.
The flow cell was prepared according to manufacturer's instructions using a Paired-End Cluster Generation Kit (Pe-103-1001) and a Cluster Station. Sequencing reactions were performed on a 1G Genome Analyzer equipped with a Paired-End Module (Illumina, Inc., Cambridge, UK). 5 pM of the library were used to achieve ~20,000 to 25,000 clusters per tile. Capillary sequencing of avrE, HrpW1 and other individual genes was done on an ABI 3730. PCR products were directly sequenced after treatment with ExoI and SAP. Primer sequences are available upon request from JHC.
Verification of Illumina sequence data
Three of the core hop genes in Pta 11528 appeared to be degenerate, based on the de novo assembly of short Illumina sequence reads. The avrE1 gene appeared to have a 20-nucleotide deletion, hrpZ1 a 325-nucleotide deletion, whilst hrpW1 appeared to have three insertions of 22, 6 and 12 nucleotides. Currently, the reliability of de novo sequence assembly from short Illumina reads has not been fully characterised. In particular, repetitive and low-complexity sequence might generate artefacts in assembled supercontigs. Therefore, we checked these putative insertions and deletions by aligning the Illumina sequence reads against the relevant regions of both the Pph 1448A reference genome sequence and our Pta 11528 assembly. As an additional control, we also performed Velvet assembies on previously published Illumina short-read data from Psy B728a . We found that the B728a avreE1, hrpZ1, hrpW1 and hopAF1 were assembled intact [Additional file 4: Figure S2], indicating that there is nothing inherently 'un-assemble-able' about these gene sequences. Sequence alignment is much more robust than de novo assembly and is not subject to assembly artefacts. The alignments supported the presence of a large deletion in hrpZ1. However, the alignments were not consistent with the assembly for avrE1 and hrpW1. Therefore, we amplified the Pta 11528 avrE1 and hrpW1 genes by PCR and verified their sequences by capillary sequencing [Additional file 5: Table S3]. This confirmed that the apparent deletion in avrE1 was an artefact of the de novo assembly and that the avrE1 sequence encodes a full-length protein product. Furthermore, transient expression of avrE1 in N. benthamiana induces cell death (S. Gimenez Ibanez and J. Rathjen, unpublished). Capillary sequencing also confirmed that the de novo assembly of hrpW1 was incorrect and that Pta 11528 encodes a full-length HrpW1 protein, albeit with repetitive sequence insertions of 69 and 12 nucleotides relative to the Pph 1448A sequence.
The absence of hopAF1 from Pta 11528 is supported not only by the de novo assembly, but also by the absence of aligned (unassembled) reads. As an additional control for the degeneracy of hopAF1 and hrpZ1, we performed the same bioinformatics and sequencing protocols to Psy B728a  and recovered hopAF1 and hrpZ1 intact in the de novo assembly assembly (Additional file 4: Figure S1).
In addition to the data available from Genbank accession ACHU00000000, the Velvet assembly and predicted protein sequences are provided in FastA format in Additional files 6 and Additional file 7.
We used GenomeMatcher  for generating and visualising whole-genome alignments. For aligning short Illumina sequence reads against a reference genome, we used MAQ  and for other sequence alignments and searches we used BLAST . We used previously published complete genomes as reference sequences for comParative analyses [81–85].
List of abbreviations
conserved effector locus
exchangeable effector locus
hidden Markov model
methyl-accepting chemotaxis protein
polymerase chain reaction
Pseudomonas syringae pathovar maculicola
Pseudomonas syringae pathovar phaseolicola
Pph genomic island 1
Pseudomonas syringae pathovar tabaci
Pseudomonas syringae pathovar syringae
Pseudomonas syringae pathovar tomato
virus-induced gene silencing
We would like to thank Dr. Larry Arnold and Theresa Law for their assistance in the functional screen for candidate type III effector genes of Pta 11528, Jodie Pike for performing the Illumina sequencing, Eric Kemen for advice and help with Illumina sequencing and Kee Hoon Sohn for handling and maintaining the bacterial strains. We thank Ashley Chu and Caitlin Thireault for assistance with capillary sequencing and Dr. Jonathan Urbach for guidance on searching for hrp-boxes. We thank Robert Jackson, Gail Preston and two anonymous reviewers for very helpful comments on the manuscript. Financial support by the Gatsby Charitable Foundation is gratefully acknowledged.
The Sainsbury Laboratory
Department of Biology, CB# 3280, Coker Hall, The University of North Carolina at Chapel Hill
Department of Botany and Plant Pathology, Oregon State University
Center for Genome Research and Biocomputing, Oregon State University
Gasson MJ: Indicator Technique for Antimetabolic Toxin Production by Phytopathogenic Species ofPseudomonas. Appl Environ Microbiol 1980, 39:25–29.PubMed
Ribeiro R de LD, Hagedorn DJ, Durbin RD, Uchytil TF: Characterization of the Bacterium Inciting Bean Wildfire in Brazil.Phytopathology 1979, 69:208–212.View Article
Thomas MD, Langston-Unkefer PJ, Uchytil TF, Durbin RD: Inhibition of Glutamine Synthetase from Pea by Tabtoxinine-beta-lactam.Plant Physiol 1983, 71:912–915.View ArticlePubMed
Durbin RD: Bacterial phytotoxins: Mechanisms of action.Cell Mol Life Sci 1991, 47:776–783.View Article
Turner JG, Taha RR: Contribution of tabtoxin to the pathogenicity ofPseudomonas syringaepv.tabaci. Physiol plant pathol 1984, 25:55–69.View Article
Chang JH, Goel AK, Grant SR, Dangl JL: Wake of the flood: ascribing functions to the wave of type III effector proteins of phytopathogenic bacteria.Curr Opin Microbiol 2004, 7:11–18.View ArticlePubMed
Collmer A, Lindeberg M, Petnicki-Ocwieja T, Schneider DJ, Alfano JR: Genomic mining type III secretion system effectors inPseudomonas syringaeyields new picks for all TTSS prospectors.Trends Microbiol 2002, 10:462–469.View ArticlePubMed
Cornelis GR, Van Gijsegem F: Assembly and function of type III secretory systems.Annu Rev Microbiol 2000, 54:735–774.View ArticlePubMed
Cunnac S, Lindeberg M, Collmer A:Pseudomonas syringaetype III secretion system effectors: repertoires in search of functions.Curr Opin Microbiol 2009, 12:53–60.View ArticlePubMed
Grant SR, Fisher EJ, Chang JH, Mole BM, Dangl JL: Subterfuge and manipulation: type III effector proteins of phytopathogenic bacteria.Annu Rev Microbiol 2006, 60:425–449.View ArticlePubMed
Greenberg JT, Vinatzer BA: Identifying type III effectors of plant pathogens and analyzing their interaction with plant cells.Curr Opin Microbiol 2003, 6:20–28.View ArticlePubMed
Jin Q, Thilmony R, Zwiesler-Vollick J, He SY: Type III protein secretion inPseudomonas syringae. Microbes Infect 2003, 5:301–310.View ArticlePubMed
Lindeberg M, Cartinhour S, Myers CR, Schechter LM, Schneider DJ, Collmer A: Closing the circle on the discovery of genes encoding Hrp regulon members and type III secretion system effectors in the genomes of three modelPseudomonas syringaestrains.Mol Plant Microbe Interact 2006, 19:1151–1158.View ArticlePubMed
Lindeberg M, Myers CR, Collmer A, Schneider DJ: Roadmap to new virulence deTerminants inPseudomonas syringae: insights from comParative genomics and genome organization.Mol Plant Microbe Interact 2008, 21:685–700.View ArticlePubMed
Wei CF, Kvitko BH, Shimizu R, Crabill E, Alfano JR, Lin NC, Martin GB, Huang HC, Collmer A: APseudomonas syringaepv. tomato DC3000 mutant lacking the type III effector HopQ1–1 is able to cause disease in the model plantNicotiana benthamiana. Plant J 2007, 51:32–46.View ArticlePubMed
de Torres M, Mansfield JW, Grabov N, Brown IR, Ammouneh H, Tsiamis G, Forsyth A, Robatzek S, Grant M, Boch J:Pseudomonas syringaeeffector AvrPtoB suppresses basal defence in Arabidopsis.Plant J 2006, 47:368–382.View ArticlePubMed
Peñaloza-Vázquez A, Preston GM, Collmer A, Bender CL: Regulatory interactions between the Hrp type III protein secretion system and coronatine biosynthesis inPseudomonas syringaepv.tomatoDC3000.Microbiology 2000, 146:2447–2456.PubMed
Buell CR, Joardar V, Lindeberg M, Selengut J, Paulsen IT, Gwinn ML, Dodson RJ, Deboy RT, Durkin AS, Kolonay JF, Madupu R, Daugherty S, Brinkac L, Beanan MJ, Haft DH, Nelson WC, Davidsen T, Zafar N, Zhou L, Liu J, Yuan Q, Khouri H, Fedorova N, Tran B, Russell D, Berry K, Utterback T, Van Aken SE, Feldblyum TV, D'Ascenzo M, Deng WL, Ramos AR, Alfano JR, Cartinhour S, Chatterjee AK, Delaney TP, Lazarowitz SG, Martin GB, Schneider DJ, Tang X, Bender CL, White O, Fraser CM, Collmer A: The complete genome sequence of theArabidopsisand tomato pathogenPseudomonas syringaepv.tomatoDC3000.Proc Natl Acad Sci USA 2003, 100:10181–10186.View ArticlePubMed
Feil H, Feil WS, Chain P, Larimer F, DiBartolo G, Copeland A, Lykidis A, Trong S, Nolan M, Goltsman E, Thiel J, Malfatti S, Loper JE, Lapidus A, Detter JC, Land M, Richardson PM, Kyrpides NC, Ivanova N, Lindow SE: Comparison of the complete genome sequences ofPseudomonas syringaepv.syringaeB728a and pv.tomatoDC3000.Proc Natl Acad Sci USA 2005, 102:11064–11069.View ArticlePubMed
Joardar V, Lindeberg M, Jackson RW, Selengut J, Dodson R, Brinkac LM, Daugherty SC, Deboy R, Durkin AS, Giglio MG, Madupu R, Nelson WC, Rosovitz MJ, Sullivan S, Crabtree J, Creasy T, Davidsen T, Haft DH, Zafar N, Zhou L, Halpin R, Holley T, Khouri H, Feldblyum T, White O, Fraser CM, Chatterjee AK, Cartinhour S, Schneider DJ, Mansfield J, Collmer A, Buell CR: Whole-genome sequence analysis ofPseudomonas syringaepv.phaseolicola1448A reveals divergence among pathovars in genes involved in virulence and transposition.J Bacteriol 2005, 187:6488–6498.View ArticlePubMed
Gardan L, Shafik H, Belouin S, Broch R, Grimont F, Grimont PA: DNA relatedness among the pathovars ofPseudomonas syringaeand description ofPseudomonas tremae sp. nov. andPseudomonas cannabinasp.nov. (exSutic and Dowson 1959).Int J Syst Bacteriol 1999, 49:469–478.View ArticlePubMed
Sarkar SF, Guttman DS: Evolution of the core genome ofPseudomonas syringae, a highly clonal, endemic plant pathogen.Appl Environ Microbiol 2004, 70:1999–2012.View ArticlePubMed
Reinhardt JA, Baltrus DA, Nishimura MT, Jeck WR, Jones CD, Dangl JL:De novoassembly using low-coverage short read sequence data from the rice pathogenPseudomonas syringaepv.oryzae. Genome Res 2009, 19:294–305.View ArticlePubMed
Almeida NF, Yan S, Lindeberg M, Studholme DJ, Schneider DJ, Condon B, Liu H, Viana CJ, Warren A, Evans C, Kemen E, Maclean D, Angot A, Martin GB, Jones JD, Collmer A, Setubal JC, Vinatzer BA: A draft genome sequence ofPseudomonas syringaepv.tomatoT1 reveals a type III effector repertoire significantly divergent from that ofPseudomonas syringaepv.tomatoDC3000.Mol Plant Microbe Interact 2009, 22:52–62.View ArticlePubMed
Inoue Y, Takikawa Y: ThehrpZandhrpAgenes are variable, and useful for groupingPseudomonas syringaebacteria.J Gen Plant Pathol 2006, 72:26–33.View Article
Taguchi F, Tanaka R, Kinoshita S, Ichinose Y, Imura Y, Andi S, Toyoda K, Tomonori Shiraishi T, Yamada T: HarpinpstafromPseudomonas syringaepv.tabaciIs Defective and Deficient in Its Expression and HR-inducing Activity.Pseudomonas syringae 2001, 67:116–123.
Holt RA, Jones SJ: The new Paradigm of flow cell sequencing.Genome Res 2008, 18:839–846.View ArticlePubMed
MacLean D, Jones JDG, Studholme DJ: Application of 'next-generation' sequencing technologies to microbial genetics.Nature Rev Microbiol 2009, 7:287–296.
Aury JM, Cruaud C, Barbe V, Rogier O, Mangenot S, Samson G, Poulain J, Anthouard V, Scarpelli C, Artiguenave F, Wincker P: High quality draft sequences for prokaryotic genomes using a mix of new sequencing technologies.BMC Genomics 2008, 9:603.View ArticlePubMed
Butler J, MacCallum I, Kleber M, Shlyakhter IA, Belmonte MK, Lander ES, Nusbaum C, Jaffe DB: ALLPATHS: de novo assembly of whole-genome shotgun microreads.Genome Res 2008, 18:810–820.View ArticlePubMed
Chaisson M, Pevzner P, Tang H: Fragment assembly with short reads.Bioinformatics 2004, 20:2067–2074.View ArticlePubMed
Chaisson MJ, Pevzner PA: Short read fragment assembly of bacterial genomes.Genome Res 2008, 18:324–330.View ArticlePubMed
Dohm JC, Lottaz C, Borodina T, Himmelbauer H: SHARCGS, a fast and highly accurate short-read assembly algorithm forde novogenomic sequencing.Genome Res 2007, 17:1697–706.View ArticlePubMed
Farrer RA, Kemen E, Jones JDG, Studholme DJ:De novoassembly of thePseudomonas syringaepv.syringaeB728a genome using Illumina/Solexa short sequence reads.FEMS Microbiol Lett 2009, 291:103–111.View ArticlePubMed
Hernandez D, Franvßois P, Farinelli L, Ostervs M, Schrenzel J: De novo bacterial genome sequencing: millions of very short reads assembled on a desktop computer.Genome Res 2008, 18:802–809.View ArticlePubMed
Jeck WR, Reinhardt JA, Baltrus DA, Hickenbotham MT, Magrini V, Mardis ER, Dangl JL, Jones CD: Extending assembly of short DNA sequences to handle error.Bioinformatics 2007, 23:2942–2944.View ArticlePubMed
Salzberg SL, Sommer DD, Puiu D, Lee VT: Gene-boosted assembly of a novel bacterial genome from very short reads.PLoS Comput Biol 2008, 4:e1000186.View ArticlePubMed
Sundquist A, Ronaghi M, Tang H, Pevzner P, Batzoglou S: Whole-genome sequencing and assembly with high-throughput, short-read technologies.PLoS ONE 2007, 2:e484.View ArticlePubMed
Warren RL, Sutton GG, Jones SJ, Holt R: Assembling millions of short DNA sequences using SSAKE.Bioinformatics 2007, 23:500–501.View ArticlePubMed
Zerbino DR, Birney E: Velvet: algorithms forde novoshort read assembly using de Bruijn graphs.Genome Res 2008, 18:821–829.View ArticlePubMed
Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL: Versatile and open software for comparing large genomes.Genome Biol 2004, 5:R12.View ArticlePubMed
Lindeberg M, Stavrinides J, Chang JH, Alfano JR, Collmer A, Dangl JL, Greenberg JT, Mansfield JW, Guttman DS: Proposed guidelines for a unified nomenclature and phylogenetic analysis of type III Hop effector proteins in the plant pathogenPseudomonas syringae. Mol Plant Microbe Interact 2005, 18:275–282.View ArticlePubMed
Alfano JR, Charkowski AO, Deng WL, Badel JL, Petnicki-Ocwieja T, van Dijk K, Collmer A: ThePseudomonas syringaeHrp 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 USA 2000, 97:4856–4861.View ArticlePubMed
Tampakaki AP, Panopoulos NJ: Elicitation of hypersensitive cell death by extracellularly targeted HrpZPsph produced in planta.Mol Plant Microbe Interact 2000, 13:1366–74.View ArticlePubMed
Preston G, Huang HC, He SY, Collmer A: The HrpZ proteins ofPseudomonas syringaepvs.syringae, glycinea, andtomatoare encoded by an operon containingYersinia yschomologs and elicit the hypersensitive response in tomato but not soybean.Mol Plant Microbe Interact 1995, 8:717–732.View ArticlePubMed
Lee J, Klessig DF, Nürnberger T: A harpin binding site in tobacco plasma membranes mediates activation of the pathogenesis-related gene HIN1 independent of extracellular calcium but dependent on mitogen-activated protein kinase activity.Plant Cell 2001, 13:1079–1093.View ArticlePubMed
Li CM, Haapalainen M, Lee J, Nürnberger T, Romantschuk M, Taira S: Harpin ofPseudomonas syringaepv.phaseolicolaharbors a protein binding site.Mol Plant Microbe Interact 2005, 18:60–66.View ArticlePubMed
Maurelli AT: Black holes, antivirulence genes, and gene inactivation in the evolution of bacterial pathogens.FEMS Microbiol Lett 2007, 267:1–8.View ArticlePubMed
Maurelli AT, Fernández RE, Bloch CA, Rode CK, Fasano A: "Black holes" and bacterial pathogenicity: a large genomic deletion that enhances the virulence ofShigellaspp. and enteroinvasiveEscherichia coli. Proc Natl Acad Sci USA 1998, 95:3943–3948.View ArticlePubMed
Casalino M, Latella MC, Prosseda G, Ceccarini P, Grimont F, Colonna B: Molecular evolution of the lysine decarboxylase-defective phenotype inShigella sonnei. Int J Med Microbiol 2005, 294:503–512.View ArticlePubMed
Day WA Jr, Fernández RE, Maurelli AT: Pathoadaptive mutations that enhance virulence: genetic organization of thecadAregions ofShigellaspp.Infect Immun 2001, 69:7471–7480.View ArticlePubMed
Tsiamis G, Mansfield JW, Hockenhull R, Jackson RW, Sesma A, Athanassopoulos E, Bennett MA, Stevens C, Vivian A, Taylor JD, Murillo J: Cultivar-specific avirulence and virulence functions assigned toavrPphFinPseudomonas syringaepv.phaseolicola, the cause of bean halo-blight disease.EMBO J 2000, 19:3204–3214.View ArticlePubMed
Jackson RW, Athanassopoulos E, Tsiamis G, Mansfield JW, Sesma A, Arnold DL, Gibbon MJ, Murillo J, Taylor JD, Vivian A: Identification of a pathogenicity island, which contains genes for virulence and avirulence, on a large native plasmid in the bean pathogenPseudomonas syringaepathovarphaseolicola. Proc Natl Acad Sci USA 1999, 96:10875–10880.View ArticlePubMed
Jackson RW, Mansfield JW, Arnold DL, Sesma A, Paynter CD, Murillo J, Taylor JD, Vivian A: Excision from tRNA genes of a large chromosomal region, carrying avrPphB, associated with race change in the bean pathogen,Pseudomonas syringaepv.phaseolicola. Mol Microbiol 2000, 38:186–197.View ArticlePubMed
Pitman AR, Jackson RW, Mansfield JW, Kaitell V, Thwaites R, Arnold DL: Exposure to host resistance mechanisms drives evolution of bacterial virulence in plants.Curr Biol 2005, 15:2230–2235.View ArticlePubMed
Stavrinides J, Guttman DS: Nucleotide sequence and evolution of the five-plasmid complement of the phytopathogenPseudomonas syringaepv.maculicolaES4326.J Bacteriol 2004, 186:5101–5115.View ArticlePubMed
Sundin GW: Genomic insights into the contribution of phytopathogenic bacterial plasmids to the evolutionary history of their hosts.Annu Rev Phytopathol 2007, 45:129–151.View ArticlePubMed
Chang JH, Urbach JM, Law TF, Arnold LW, Hu A, Gombar S, Grant SR, Ausubel FM, Dangl JL: A high-throughput, near-saturating screen for type III effector genes fromPseudomonas syringae. Proc Natl Acad Sci USA 2005, 102:2549–2554.View ArticlePubMed
Oh HS, Kvitko BH, Morello JE, Collmer A:Pseudomonas syringaelytic transglycosylases coregulated with the type III secretion system contribute to the translocation of effector proteins into plant cells.J Bacteriol 2007, 189:8277–8289.View ArticlePubMed
Von Bodman SB, Bauer WD, Coplin DL: Quorum sensing in plant-pathogenic bacteria.Annu Rev Phytopathol 2003, 41:455–482.View Article
Brencic A, Winans SC: Detection of and response to signals involved in host-microbe interactions by plant-associated bacteria.Microbiol Mol Biol Rev 2005, 69:155–294.View ArticlePubMed
Taguchi F, Shibata S, Suzuki T, Ogawa Y, Aizawa S, Takeuchi K, Ichinose Y: Effects of glycosylation on swimming ability and flagellar polymorphic transformation inPseudomonas syringaepv.tabaci6605.J Bacteriol 2008, 190:764–768.View ArticlePubMed
Cha JY, Lee JS, Oh JI, Choi JW, Baik HS: Functional analysis of the role of Fur in the virulence ofPseudomonas syringaepv.tabaci11528: Fur controls expression of genes involved in quorum-sensing.Biochem Biophys Res Commun 2008, 366:281–287.View ArticlePubMed
Jakob K, Kniskern JM, Bergelson J: The role of pectate lyase and the jasmonic acid defense response inPseudomonas viridiflavavirulence.Mol Plant Microbe Interact 2007, 20:146–158.View ArticlePubMed
Preston GM, Studholme DJ, Caldelari I: Profiling the secretomes of plant pathogenic Proteobacteria.FEMS Microbiol Rev 2005, 29:331–360.PubMed
Barabote RD, Johnson OL, Zetina E, San Francisco SK, Fralick JA, San Francisco MJ:Erwinia chrysanthemi tolCis involved in resistance to antimicrobial plant chemicals and is essential for phytopathogenesis.J Bacteriol 2003, 185:5772–5778.View ArticlePubMed
Kinscherf TG, Willis DK: The biosynthetic gene cluster for the beta-lactam antibiotic tabtoxin inPseudomonas syringae. J Antibiot 2005, 58:817–821.View ArticlePubMed
Greenberg DE, Porcella SF, Zelazny AM, Virtaneva K, Sturdevant DE, Kupko JJ 3rd, Barbian KD, Babar A, Dorward DW, Holland SM: Genome sequence analysis of the emerging human pathogenic acetic acid bacteriumGranulibacter bethesdensis. J Bacteriol 2007, 189:8727–8736.View ArticlePubMed
Arand M, Hemmer H, Dürk H, Baratti J, Archelas A, Furstoss R, Oesch F: Cloning and molecular characterization of a soluble epoxide hydrolase fromAspergillus nigerthat is related to mammalian microsomal epoxide hydrolase.Biochem J 1999, 344:273–280.View ArticlePubMed
Finn RD, Tate J, Mistry J, Coggill PC, Sammut SJ, Hotz HR, Ceric G, Forslund K, Eddy SR, Sonnhammer EL, Bateman A: The Pfam protein families database.Nucleic Acids Res 2008, 36:D281–8.View ArticlePubMed
Mandel MJ, Wollenberg MS, Stabb EV, Visick KL, Ruby EG: A single regulatory gene is sufficient to alter bacterial host range.Nature 2009, 12:215–8.View Article
Pasloske BL, Sastry PA, Finlay BB, Paranchych W: Two unusual pilin sequences from different isolates ofPseudomonas aeruginosa. J Bacteriol 1988, 170:3738–3741.PubMed
Pizarro-Cerdá J, Cossart P: Bacterial adhesion and entry into host cells.Cell 2006, 124:715–727.View ArticlePubMed
Stothard P, Wishart DS: Circular genome visualization and exploration using CGView.Bioinformatics 2005, 21:537–539.View ArticlePubMed
Stein LD, Mungall C, Shu S, Caudy M, Mangone M, Day A, Nickerson E, Stajich JE, Harris TW, Arva A, Lewis S: The generic genome browser: a building block for a model organism system database.Genome Res 2002, 12:1599–1610.View ArticlePubMed
Ohtsubo Y, Ikeda-Ohtsubo W, Nagata Y, Tsuda M: GenomeMatcher: a graphical user interface for DNA sequence comparison.BMC Bioinformatics 2008, 9:376.View ArticlePubMed
Li H, Ruan J, Durbin R: Mapping short DNA sequencing reads and calling variants using mapping quality scores.Genome Res 2008, 18:1851–1858.View ArticlePubMed
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool.J Mol Biol 1990, 215:403–410.PubMed
Camargo LE, da Silva AC, Moon DH, Takita MA, Lemos EG, Machado MA, Ferro MI, da Silva FR, Goldman MH, Goldman GH, Lemos MV, El-Dorry H, Tsai SM, Carrer H, Carraro DM, de Oliveira RC, Nunes LR, Siqueira WJ, Coutinho LL, Kimura ET, Ferro ES, Harakava R, Kuramae EE, Marino CL, Giglioti E, Abreu IL, Alves LM, do Amaral AM, Baia GS, Blanco SR, Brito MS, Cannavan FS, Celestino AV, da Cunha AF, Fenille RC, Ferro JA, Formighieri EF, Kishi LT, Leoni SG, Oliveira AR, Rosa VE Jr, Sassaki FT, Sena JA, de Souza AA, Truffi D, Tsukumo F, Yanai GM, Zaros LG, Civerolo EL, Simpson AJ, Almeida NF Jr, Setubal JC, Kitajima JP: ComParative analyses of the complete genome sequences of Pierce's disease and citrus variegated chlorosis strains ofXylella fastidiosa. J Bacteriol(Edited by: Young JP, Crossman LC, Johnston AW, Thomson NR, Ghazoui ZF, Hull KH, Wexler M). 2003, 185:1018–26.View ArticlePubMed
Curson AR, Todd JD, Poole PS, Mauchline TH, East AK, Quail MA, Churcher C, Arrowsmith C, Cherevach I, Chillingworth T, Clarke K, Cronin A, Davis P, Fraser A, Hance Z, Hauser H, Jagels K, Moule S, Mungall K, Norbertczak H, Rabbinowitsch E, Sanders M, Simmonds M, Whitehead S, Parkhill J: The genome ofRhizobium leguminosarumhas recognizable core and accessory components.Genome Biol 2006, 7:R34.View ArticlePubMed
Hayden HS, Gillett W, Saenphimmachak C, Lim R, Zhou Y, Jacobs MA, Chang J, Rohmer L, D'Argenio DA, Palmieri A, Levy R, Haugen E, Wong GK, Brittnacher MJ, Burns JL, Miller SI, Olson MV, Kaul R: Large-insert genome analysis technology detects structural variation inPseudomonas aeruginosaclinical strains from cystic fibrosis patients.Genomics 2008, 91:530–537.View ArticlePubMed
Iacono M, Villa L, Fortini D, Bordoni R, Imperi F, Bonnal RJ, Sicheritz-Ponten T, De Bellis G, Visca P, Cassone A, Carattoli A: Whole-genome pyrosequencing of an epidemic multidrug-resistantAcinetobacter baumanniistrain belonging to the European clone II group.Antimicrob Agents Chemother 2008, 52:2616–2625.View ArticlePubMed
Salanoubat M, Genin S, Artiguenave F, Gouzy J, Mangenot S, Arlat M, Billault A, Brottier P, Camus JC, Cattolico L, Chandler M, Choisne N, Claudel-Renard C, Cunnac S, Demange N, Gaspin C, Lavie M, Moisan A, Robert C, Saurin W, Schiex T, Siguier P, Thébault P, Whalen M, Wincker P, Levy M, Weissenbach J, Boucher CA: Genome sequence of the plant pathogenRalstonia solanacearum. Nature 2002, 415:497–502.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the Terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.