Genome sequence analyses of Pseudomonas savastanoi pv. glycinea and subtractive hybridization-based comparative genomics with nine pseudomonads

Comparative Genomics of Prunus-Associated Members of the Pseudomonas syringae Species Complex Reveals Traits Supporting Co-evolution and Host Adaptation

Members of the Pseudomonas syringae species complex cause symptoms that are ranging from leaf spots to cankers on a multitude of plant species, including some of the genus Prunus. To date, a total of two species of the P. syringae species complex and six different pathovars have been associated with diseases on Prunus spp., which were shown to belong to different phylogenetic units (phylogroups, PG) based on sequence similarity of housekeeping genes or whole genomes, suggesting that virulence to Prunus spp. may be the result of convergent pathoadaptation. In this study, a comparative genomics approach was used to determine genes significantly associated with strains isolated from Prunus spp. across a phylogeny of 97 strains belonging to the P. syringae species complex. Our study revealed the presence of a set of orthologous proteins which were significantly associated with strains isolated from Prunus spp. than in strains isolated from other hosts or from non-agricultural environments. Among them, the type III effector HopAY predicted to encode for a C58 cysteine protease was found to be highly associated with strains isolated from Prunus spp. and revealed patterns supporting co-evolution and host adaptation.

Bacteriophage Control of Pseudomonas savastanoi pv. glycinea in Soybean

Bacterial viruses (bacteriophages) have been considered as potential agents for the biological control of bacterial phytopathogens due to their safety and host specificity. Pseudomonas savastanoi pv. glycinea (Psg) is a causative agent of the bacterial spotting of soybean (Glycine max Willd). The harm caused by this bacterium to crop production and the development of antibiotic resistance in Psg and other pathogenic microorganisms has led to the pursuit of alternative management strategies. In this study, three Psg-specific lytic bacteriophages were isolated from soybean field soil in geographically distant regions of Russia, and their potential for protective action on plants was assessed. Sequencing of phage genomes has revealed their close relatedness and attribution to the genus Ghunavirus, subfamily Studiervirinae, family Autographiviridae. Extensive testing of the biological properties of P421, the representative of the isolated phage group, has demonstrated a relatively broad host range covering closely related Pseudomonas species and stability over wide temperature (4-40 °C) and pH (pH 4-7) ranges, as well as stability under ultraviolet irradiation for 30 min. Application of the phages to prevent, and treat, Psg infection of soybean plants confirms that they are promising as biocontrol agents.

Symbiotic bacteria of the gall-inducing mite Fragariocoptes setiger (Eriophyoidea) and phylogenomic resolution of the eriophyoid position among Acari

Eriophyoid mites represent a hyperdiverse, phytophagous lineage with an unclear phylogenetic position. These mites have succeeded in colonizing nearly every seed plant species, and this evolutionary success was in part due to the mites' ability to induce galls in plants. A gall is a unique niche that provides the inducer of this modification with vital resources. The exact mechanism of gall formation is still not understood, even as to whether it is endogenic (mites directly cause galls) or exogenic (symbiotic microorganisms are involved). Here we (i) investigate the phylogenetic affinities of eriophyoids and (ii) use comparative metagenomics to test the hypothesis that the endosymbionts of eriophyoid mites are involved in gall formation. Our phylogenomic analysis robustly inferred eriophyoids as closely related to Nematalycidae, a group of deep-soil mites belonging to Endeostigmata. Our comparative metagenomics, fluorescence in situ hybridization, and electron microscopy experiments identified two candidate endosymbiotic bacteria shared across samples, however, it is unlikely that they are gall inducers (morphotype1: novel Wolbachia, morphotype2: possibly Agrobacterium tumefaciens). We also detected an array of plant pathogens associated with galls that may be vectored by the mites, and we determined a mite pathogenic virus (Betabaculovirus) that could be tested for using in biocontrol of agricultural pest mites.

Multiple virulence factors regulated by AlgU contribute to the pathogenicity of Pseudomonas savastanoi pv. glycinea in soybean

Pseudomonas savastanoi pv. glycinea (Psg) causes bacterial blight of soybean. To identify candidate virulence factors, transposon-mediated mutational analysis of Psg was carried out. We syringe-inoculated soybean leaves with Psg transposon mutants and identified 28 mutants which showed reduced virulence from 1,000 mutants screened. Next, we spray-inoculated soybean leaves with these mutants and demonstrated that the algU mutant showed significantly reduced virulence together with reduced bacterial populations in planta. Expression profiles comparison between the Psg wild-type (WT) and algU mutant in HSC broth revealed that expression of coronatine (COR)-related genes (including cmaA and corR) were down-regulated in the algU mutant compared with Psg WT. Moreover, we also showed that COR production were reduced in the algU mutant compared with WT. We also demonstrated that algD, which is related to alginate biosynthesis, showed reduced expression and biofilm formation was significantly suppressed in the algU mutant. Furthermore, hrpL also showed less expression in the algU mutant. These results indicate that AlgU plays a critical role in promoting Psg pathogenesis by regulating multiple virulence factors.

Analysis of biosorption and biotransformation mechanism of Pseudomonas chengduensis strain MBR under Cd(II) stress from genomic perspective

Microbial treatment of heavy metal-polluted sites is considered an environmentally friendly bioremediation technology with high potential. This study shows that Pseudomonas chengduensis strain MBR, a bacterium that can potentially be applied in the treatment of heavy metal pollution, is most affected by Cd(II) stress at the beginning of its growth. Up to 100% of total Cd(II) adsorption occurs in the first 48 h after treatment of stationary phase cells with Cd(II). A biofilm forms on the cell surface, Cd(II) adsorbs, and is reduced to Cd (0) in the form of nanoscale particles. The genome of strain MBR was sequenced, annotated and analyzed. We identified various genes potentially related to cadmium resistance, transport and metabolism. Analysis of the strain MBR genome is helpful to explore the mechanism of Cd(II) resistance, and can provide new ideas for cadmium pollution control.

Applications of Next-Generation Sequencing for Large-Scale Pathogen Diagnoses in Soybean

Soybean (Glycine max) has become an important crop in Manitoba, Canada, with a 10-fold increase in dedicated acreage over the past decade. Given the rapid increase in production, scarce information about foliar diseases present in the province has been recorded. In order to describe the foliar pathogens affecting this legume, we harnessed next-generation sequencing (NGS) to carry out a comprehensive survey across Manitoba in 2016. Fields were sampled during the V2/3 (33 fields) and R6 (70 fields) growth stages, with at least three symptomatic leaves per field collected and subjected to RNA sequencing. We successfully detected several bacteria, fungi, and viruses known to infect soybean, including Pseudomonas savastanoi pv. glycinea, Septoria glycines, and Peronospora manshurica, as well as pathogens not previously identified in the province (e.g., Pseudomonas syringae pv. tabaci, Cercospora sojina, and Bean yellow mosaic virus). For some microorganisms, we were able to disentangle the different pathovars present and/or assemble their genome sequence. Since NGS generates data on the entire flora and fauna occupying a leaf sample, we also identified residual pathogens (i.e., pathogens of crops other than soybean) and multiple species of arthropod pests. Finally, the sequence information produced by NGS allowed for the development of polymerase chain reaction-based diagnostics for some of the most widespread and important pathogens. Although there are many benefits of using NGS for large-scale plant pathogen diagnoses, we also discuss some of the limitations of this technology.

Many plant pathogenic Pseudomonas savastanoi pv glycinea isolates possess an inactive quorum sensing ahlR gene via a point mutation

Many plant bacterial pathogens monitor their group behaviour and their population density via production of N-acyl homoserine lactone signals which regulate the expression of several genes via the LuxI/R homologs. This regulatory network, termed quorum sensing (QS), is present in the soybean bacterial pathogen Pseudomonas savastanoi pv glycinea (Psg). The sequenced genomes of two strains of Psg, race 4 and B076, contain an N-acyl homoserine lactone (AHL) based LuxI/R QS system named AhlI/R. While studying the QS system of Psg strains race 4 and B076 isolated in USA, LMG5066 in New Zealand and IBSBF355 in Brazil, we found that B076, LMG5066 and IBSBF355 possess a point mutation in the ahlR gene that causes a frameshift resulting in a truncated AhlR protein. Psg race 4 does not possess the mutation in ahlR and the QS system is functional. The same mutation in the ahlR gene was found to be also present in 9 of 19 Psg strains isolated from diseased soybean in Illinois. Phenotypic analysis of strains showed that swarming motility is repressed whereas phosphate solubilisation was activated by QS in Psg. Analysing the secretome, we also found that four proteins were under QS regulation.

Comparative genomics of Pseudomonas syringae reveals convergent gene gain and loss associated with specialization onto cherry (Prunus avium)

Genome-wide analyses of the effector- and toxin-encoding genes were used to examine the phylogenetics and evolution of pathogenicity amongst diverse strains of Pseudomonas syringae causing bacterial canker of cherry (Prunus avium), including pathovars P. syringae pv morsprunorum (Psm) races 1 and 2, P. syringae pv syringae (Pss) and P. syringae pv avii. Phylogenetic analyses revealed Psm races and P. syringae pv avii clades were distinct and were each monophyletic, whereas cherry-pathogenic strains of Pss were interspersed amongst strains from other host species. A maximum likelihood approach was used to predict effectors associated with pathogenicity on cherry. Pss possesses a smaller repertoire of type III effectors but has more toxin biosynthesis clusters than Psm and P. syringae pv avii. Evolution of cherry pathogenicity was correlated with gain of genes such as hopAR1 and hopBB1 through putative phage transfer and horizontal transfer respectively. By contrast, loss of the avrPto/hopAB redundant effector group was observed in cherry-pathogenic clades. Ectopic expression of hopAB and hopC1 triggered the hypersensitive reaction in cherry leaves, confirming computational predictions. Cherry canker provides a fascinating example of convergent evolution of pathogenicity that is explained by the mix of effector and toxin repertoires acting on a common host.

The coronafacoyl phytotoxins: structure, biosynthesis, regulation and biological activities

Phytotoxins are secondary metabolites that contribute to the development and/or severity of diseases caused by various plant pathogenic microorganisms. The coronafacoyl phytotoxins are an important family of plant toxins that are known or suspected to be produced by several phylogenetically distinct plant pathogenic bacteria, including the gammaproteobacterium Pseudomonas syringae and the actinobacterium Streptomyces scabies. At least seven different family members have been identified, of which coronatine was the first to be described and is the best-characterized. Though nonessential for disease development, coronafacoyl phytotoxins appear to enhance the severity of disease symptoms induced by pathogenic microbes during host infection. In addition, the identification of coronafacoyl phytotoxin biosynthetic genes in organisms not known to be plant pathogens suggests that these metabolites may have additional roles other than as virulence factors. This review focuses on our current understanding of the structures, biosynthesis, regulation, biological activities and evolution of coronafacoyl phytotoxins as well as the different methods that are used to detect these metabolites and the organisms that produce them.

Coronafacoyl Phytotoxin Biosynthesis and Evolution in the Common Scab Pathogen Streptomyces scabiei

Coronafacoyl phytotoxins are an important family of plant toxins that are produced by several different phytopathogenic bacteria, including the gammaproteobacterium Pseudomonas syringae and the actinobacterium Streptomyces scabiei (formerly Streptomyces scabies). The phytotoxins consist of coronafacic acid (CFA) linked via an amide bond to different amino acids or amino acid derivatives. Previous work suggested that S. scabiei and P. syringae use distinct biosynthetic pathways for producing CFA, which is subsequently linked to its amino acid partner to form the complete phytotoxin. Here, we provide further evidence that the S. scabiei CFA biosynthetic pathway is novel by characterizing the role of CYP107AK1, a predicted cytochrome P450 that has no homologue in P. syringae Deletion of the CYP107AK1 gene abolished production of coronafacoyl-isoleucine (CFA-Ile), the primary coronafacoyl phytotoxin produced by S. scabiei Structural elucidation of accumulated biosynthetic intermediates in the ΔCYP107AK1 mutant indicated that CYP107AK1 is required for introducing the oxygen atom that ultimately forms the carbonyl group in the CFA backbone. The CYP107AK1 gene along with two additional genes involved in CFA-Ile biosynthesis in S. scabiei were found to be associated with putative CFA biosynthetic genes in other actinobacteria but not in other organisms. Analysis of the overall genetic content and organization of known and putative CFA biosynthetic gene clusters, together with phylogenetic analysis of the core biosynthetic genes, indicates that horizontal gene transfer has played an important role in the dissemination of the gene cluster and that rearrangement, insertion, and/or deletion events have likely contributed to the divergent biosynthetic evolution of coronafacoyl phytotoxins in bacteria.IMPORTANCE The ability of plants to defend themselves against invading pathogens relies on complex signaling pathways that are controlled by key phytohormones such as jasmonic acid (JA). Some phytopathogenic bacteria have evolved the ability to manipulate JA signaling in order to overcome host defenses by producing coronatine (COR), which functions as a potent JA mimic. COR and COR-like molecules, collectively referred to as coronafacoyl phytotoxins, are produced by several different plant-pathogenic bacteria, and this study provides supporting evidence that different biosynthetic pathways are utilized by different bacteria for production of these phytotoxins. In addition, our study provides a greater understanding of how coronafacoyl phytotoxin biosynthesis may have evolved in phylogenetically distinct bacteria, and we demonstrate that production of these compounds may be more widespread than previously recognized and that their role for the producing organism may not be limited to host-pathogen interactions.

Complete sequence and comparative genomic analysis of eight native Pseudomonas syringae plasmids belonging to the pPT23A family

BACKGROUND

The pPT23A family of plasmids appears to be indigenous to the plant pathogen Pseudomonas syringae and these plasmids are widely distributed and widely transferred among pathovars of P. syringae and related species. pPT23A-family plasmids (PFPs) are sources of accessory genes for their hosts that can include genes important for virulence and epiphytic colonization of plant leaf surfaces. The occurrence of repeated sequences including duplicated insertion sequences on PFPs has made obtaining closed plasmid genome sequences difficult. Therefore, our objective was to obtain complete genome sequences from PFPs from divergent P. syringae pathovars and also from strains of P. syringae pv. syringae isolated from different hosts.

RESULTS

The eight plasmids sequenced ranged in length from 61.6 to 73.8 kb and encoded from 65 to 83 annotated orfs. Virulence genes including type III secretion system effectors were encoded on two plasmids, and one of these, pPt0893-29 from P. syringae pv. tabaci, encoded a wide variety of putative virulence determinants. The PFPs from P. syringae pv. syringae mostly encoded genes of importance to ecological fitness including the rulAB determinant conferring tolerance to ultraviolet radiation. Heavy metal resistance genes encoding resistance to copper and arsenic were also present in a few plasmids. The discovery of part of the chromosomal genomic island GI6 from P. syringae pv. syringae B728a in two PFPs from two P. syringae pv. syringae hosts is further evidence of past intergenetic transfers between plasmid and chromosomal DNA. Phylogenetic analyses also revealed new subgroups of the pPT23A plasmid family and confirmed that plasmid phylogeny is incongruent with P. syringae pathovar or host of isolation. In addition, conserved genes among seven sequenced plasmids within the same phylogenetic group were limited to plasmid-specific functions including maintenance and transfer functions.

CONCLUSIONS

Our sequence analysis further revealed that PFPs from P. syringae encode suites of accessory genes that are selected at species (universal distribution), pathovar (interpathovar distribution), and population levels (intrapathovar distribution). The conservation of type IV secretion systems encoding conjugation functions also presumably contributes to the distribution of these plasmids within P. syringae populations.

Jasmonate signaling and manipulation by pathogens and insects

Plants synthesize jasmonates (JAs) in response to developmental cues or environmental stresses, in order to coordinate plant growth, development or defense against pathogens and herbivores. Perception of pathogen or herbivore attack promotes synthesis of jasmonoyl-L-isoleucine (JA-Ile), which binds to the COI1-JAZ receptor, triggering the degradation of JAZ repressors and induction of transcriptional reprogramming associated with plant defense. Interestingly, some virulent pathogens have evolved various strategies to manipulate JA signaling to facilitate their exploitation of plant hosts. In this review, we focus on recent advances in understanding the mechanism underlying the enigmatic switch between transcriptional repression and hormone-dependent transcriptional activation of JA signaling. We also discuss various strategies used by pathogens and insects to manipulate JA signaling and how interfering with this could be used as a novel means of disease control.

WHOP, a Genomic Region Associated With Woody Hosts in the Pseudomonas syringae Complex Contributes to the Virulence and Fitness of Pseudomonas savastanoi pv. savastanoi in Olive Plants

Bacteria from the Pseudomonas syringae complex belonging to phylogroups 1 and 3 (PG1 and PG3, respectively) isolated from woody hosts share a genomic region herein referred to as WHOP (from woody host and Pseudomonas spp.), which is absent in strains infecting herbaceous organs. In this work, we show that this region is also encoded in P. syringae pv. actinidifoliorum (PG1) and six additional members of PG3, namely, Pseudomonas savastanoi pv. retacarpa, three P. syringae pathovars, Pseudomonas meliae, and Pseudomonas amygdali. Partial conservation of the WHOP occurs in only a few PG2 strains. In P. savastanoi pv. savastanoi NCPPB 3335, the WHOP region is organized into four operons and three independently transcribed genes. While the antABC and catBCA operons mediate the catabolism of anthranilate and catechol, respectively, the ipoABC operon confers oxygenase activity to aromatic compounds. The deletion of antABC, catBCA, or ipoABC in NCPPB 3335 caused reduced virulence in woody olive plants without affecting knot formation in nonwoody plants; catBCA, dhoAB, and PSA3335_3206 (encoding a putative aerotaxis receptor) were also required for the full fitness of this strain exclusively in woody olive plants. Overall, this study sheds light on the evolution and adaptation of bacteria from the P. syringae complex to woody hosts and highlights the enzymatic activities encoded within the WHOP region that are essential for this process.

Comparative genomics reveals genes significantly associated with woody hosts in the plant pathogen Pseudomonas syringae

The diversification of lineages within Pseudomonas syringae has involved a number of adaptive shifts from herbaceous hosts onto various species of tree, resulting in the emergence of highly destructive diseases such as bacterial canker of kiwi and bleeding canker of horse chestnut. This diversification has involved a high level of gene gain and loss, and these processes are likely to play major roles in the adaptation of individual lineages onto their host plants. In order to better understand the evolution of P. syringae onto woody plants, we have generated de novo genome sequences for 26 strains from the P. syringae species complex that are pathogenic on a range of woody species, and have looked for statistically significant associations between gene presence and host type (i.e. woody or herbaceous) across a phylogeny of 64 strains. We have found evidence for a common set of genes associated with strains that are able to colonize woody plants, suggesting that divergent lineages have acquired similarities in genome composition that may form the genetic basis of their adaptation to woody hosts. We also describe in detail the gain, loss and rearrangement of specific loci that may be functionally important in facilitating this adaptive shift. Overall, our analyses allow for a greater understanding of how gene gain and loss may contribute to adaptation in P. syringae.

Complete Sequences of IncU Plasmids Harboring Quinolone Resistance Genes qnrS2 and aac(6')-Ib-cr in Aeromonas spp. from Ornamental Fish

The nucleotide sequences of three IncU plasmids from Aeromonas spp. isolated from ornamental fish are described. They had a typical IncU backbone for plasmid replication and maintenance functions, but conjugative transfer modules were disrupted. The gene qnrS2 was inserted into mpR as a mobile insertion cassette. Novel Tn3 family transposons carrying putative toxin-antitoxin and plasmid stability genes were identified. The study demonstrates high plasticity of IncU plasmids from aquatic environments.

Host target modification as a strategy to counter pathogen hijacking of the jasmonate hormone receptor

In the past decade, characterization of the host targets of pathogen virulence factors took a center stage in the study of pathogenesis and disease susceptibility in plants and humans. However, the impressive knowledge of host targets has not been broadly exploited to inhibit pathogen infection. Here, we show that host target modification could be a promising new approach to "protect" the disease-vulnerable components of plants. In particular, recent studies have identified the plant hormone jasmonate (JA) receptor as one of the common targets of virulence factors from highly evolved biotrophic/hemibiotrophic pathogens. Strains of the bacterial pathogen Pseudomonas syringae, for example, produce proteinaceous effectors, as well as a JA-mimicking toxin, coronatine (COR), to activate JA signaling as a mechanism to promote disease susceptibility. Guided by the crystal structure of the JA receptor and evolutionary clues, we succeeded in modifying the JA receptor to allow for sufficient endogenous JA signaling but greatly reduced sensitivity to COR. Transgenic Arabidopsis expressing this modified receptor not only are fertile and maintain a high level of insect defense, but also gain the ability to resist COR-producing pathogens Pseudomonas syringae pv. tomato and P. syringae pv. maculicola. Our results provide a proof-of-concept demonstration that host target modification can be a promising new approach to prevent the virulence action of highly evolved pathogens.

Antimicrobial activity of essential oils of Thymus vulgaris and Origanum vulgare on phytopathogenic strains isolated from soybean

The aim of this work was to study the antimicrobial activity of essential oils obtained from Thymus vulgaris (thyme) and Origanum vulgare (oregano) on phytopathogenic Pseudomonas species isolated from soybean. Strains with characteristics of P. syringae were isolated from leaves of soybean plants with blight symptoms. Ten of these could be identified in Group Ia of LOPAT as P. syringae. Six of these were confirmed as P. syringae using 16S rRNA, indicating the presence of these phytopathogenic bacteria in east and central Argentina. All the phytopathogenic bacteria were re-isolated and identified from the infected plants. MIC values for thyme were 11.5 and 5.7 mg·ml(-1) on P. syringae strains, while oregano showed variability in the inhibitory activity. Both essential oils inhibited all P. syringae strains, with better inhibitory activity than the antibiotic streptomycin. The oils were not bactericidal for all pseudomonads. Both oils contained high carvacrol (29.5% and 19.7%, respectively) and low thymol (1.5%). Natural products obtained from aromatic plants represent potential sources of molecules with biological activity that could be used as new alternatives for the treatment of phytopathogenic bacteria infections.

The plant growth-promoting bacteria Azospirillum amazonense: genomic versatility and phytohormone pathway

The rhizosphere bacterium Azospirillum amazonense associates with plant roots to promote plant growth. Variation in replicon numbers and rearrangements is common among Azospirillum strains, and characterization of these naturally occurring differences can improve our understanding of genome evolution. We performed an in silico comparative genomic analysis to understand the genomic plasticity of A. amazonense. The number of A. amazonense-specific coding sequences was similar when compared with the six closely related bacteria regarding belonging or not to the Azospirillum genus. Our results suggest that the versatile gene repertoire found in A. amazonense genome could have been acquired from distantly related bacteria from horizontal transfer. Furthermore, the identification of coding sequence related to phytohormone production, such as flavin-monooxygenase and aldehyde oxidase, is likely to represent the tryptophan-dependent TAM pathway for auxin production in this bacterium. Moreover, the presence of the coding sequence for nitrilase indicates the presence of the alternative route that uses IAN as an intermediate for auxin synthesis, but it remains to be established whether the IAN pathway is the Trp-independent route. Future investigations are necessary to support the hypothesis that its genomic structure has evolved to meet the requirement for adaptation to the rhizosphere and interaction with host plants.

Virulence determinants of Pseudomonas syringae strains isolated from grasses in the context of a small type III effector repertoire

Background

Pseudomonas syringae is pathogenic to a large number of plant species. For host colonization and disease progression, strains of this bacterium utilize an array of type III-secreted effectors and other virulence factors, including small secreted molecules such as syringolin A, a peptide derivative that inhibits the eukaryotic proteasome. In strains colonizing dicotyledonous plants, the compound was demonstrated to suppress the salicylic-acid-dependent defense pathway. Here, we analyze virulence factors of three strains colonizing wheat (Triticum aestivum): P. syringae pathovar syringae (Psy) strains B64 and SM, as well as P. syringae BRIP34876. These strains have a relatively small repertoire of only seven to eleven type III secreted effectors (T3Es) and differ in their capacity to produce syringolin A. The aim of this study was to analyze the contribution of various known virulence factors in the context of a small T3E repertoire.

Results

We demonstrate that syringolin A production enhances disease symptom development upon direct infiltration of strains into wheat leaves. However, it is not universally required for colonization, as Psy SM, which lacks syringolin biosynthesis genes, reaches cell densities comparable to syringolin A producer P. syringae BRIP34876. Next, we show that despite the small set of T3E-encoding genes, the type III secretion system remains the key pathogenicity determinant in these strains, and that phenotypic effects of deleting T3E-coding genes become apparent only when multiple effectors are removed.

Conclusions

Whereas production of syringolin A is not required for successful colonization of wheat leaves by P. syringae strains, its production results in increased lesion formation. Despite the small number of known T3Es encoded by the analyzed strains, the type III secretion system is essential for endophytic growth of these strains.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-014-0304-5) contains supplementary material, which is available to authorized users.

Ribosomally encoded antibacterial proteins and peptides from Pseudomonas

Members of the Pseudomonas genus produce diverse secondary metabolites affecting other bacteria, fungi or predating nematodes and protozoa but are also equipped with the capacity to secrete different types of ribosomally encoded toxic peptides and proteins, ranging from small microcins to large tailocins. Studies with the human pathogen Pseudomonas aeruginosa have revealed that effector proteins of type VI secretion systems are part of the antibacterial armamentarium deployed by pseudomonads. A novel class of antibacterial proteins with structural similarity to plant lectins was discovered by studying antagonism among plant-associated Pseudomonas strains. A genomic perspective on pseudomonad bacteriocinogeny shows that the modular architecture of S pyocins of P. aeruginosa is retained in a large diversified group of bacteriocins, most of which target DNA or RNA. Similar modularity is present in as yet poorly characterized Rhs (recombination hot spot) proteins and CDI (contact-dependent inhibition) proteins. Well-delimited domains for receptor recognition or cytotoxicity enable the design of chimeric toxins with novel functionalities, which has been applied successfully for S and R pyocins. Little is known regarding how these antibacterials are released and ultimately reach their targets. Other remaining issues concern the identification of environmental triggers activating these systems and assessment of their ecological impact in niches populated by pseudomonads.

The extent of genome flux and its role in the differentiation of bacterial lineages

Horizontal gene transfer (HGT) and gene loss are key processes in bacterial evolution. However, the role of gene gain and loss in the emergence and maintenance of ecologically differentiated bacterial populations remains an open question. Here, we use whole-genome sequence data to quantify gene gain and loss for 27 lineages of the plant-associated bacterium Pseudomonas syringae. We apply an extensive error-control procedure that accounts for errors in draft genome data and greatly improves the accuracy of patterns of gene occurrence among these genomes. We demonstrate a history of extensive genome fluctuation for this species and show that individual lineages could have acquired thousands of genes in the same period in which a 1% amino acid divergence accrues in the core genome. Elucidating the dynamics of genome fluctuation reveals the rapid turnover of gained genes, such that the majority of recently gained genes are quickly lost. Despite high observed rates of fluctuation, a phylogeny inferred from patterns of gene occurrence is similar to a phylogeny based on amino acid replacements within the core genome. Furthermore, the core genome phylogeny suggests that P. syringae should be considered a number of distinct species, with levels of divergence at least equivalent to those between recognized bacterial species. Gained genes are transferred from a variety of sources, reflecting the depth and diversity of the potential gene pool available via HGT. Overall, our results provide further insights into the evolutionary dynamics of genome fluctuation and implicate HGT as a major factor contributing to the diversification of P. syringae lineages.

The conserved upstream region of lscB/C determines expression of different levansucrase genes in plant pathogen Pseudomonas syringae

Background

Pseudomonas syringae pv. glycinea PG4180 is an opportunistic plant pathogen which causes bacterial blight of soybean plants. It produces the exopolysaccharide levan by the enzyme levansucrase. Levansucrase has three gene copies in PG4180, two of which, lscB and lscC, are expressed while the third, lscA, is cryptic. Previously, nucleotide sequence alignments of lscB/C variants in various P. syringae showed that a ~450-bp phage-associated promoter element (PAPE) including the first 48 nucleotides of the ORF is absent in lscA.

Results

Herein, we tested whether this upstream region is responsible for the expression of lscB/C and lscA. Initially, the transcriptional start site for lscB/C was determined. A fusion of the PAPE with the ORF of lscA (lscB_UpNA) was generated and introduced to a levan-negative mutant of PG4180. Additionally, fusions comprising of the non-coding part of the upstream region of lscB with lscA (lscB_UpA) or the upstream region of lscA with lscB (lscA_UpB) were generated. Transformants harboring the lscB_UpNA or the lscB_UpA fusion, respectively, showed levan formation while the transformant carrying lscA_UpB did not. qRT-PCR and Western blot analyses showed that lscB_UpNA had an expression similar to lscB while lscB_UpA had a lower expression. Accuracy of protein fusions was confirmed by MALDI-TOF peptide fingerprinting.

Conclusions

Our data suggested that the upstream sequence of lscB is essential for expression of levansucrase while the N-terminus of LscB mediates an enhanced expression. In contrast, the upstream region of lscA does not lead to expression of lscB. We propose that lscA might be an ancestral levansucrase variant upstream of which the PAPE got inserted by potentially phage-mediated transposition events leading to expression of levansucrase in P. syringae.

Genomics-Based Exploration of Virulence Determinants and Host-Specific Adaptations of Pseudomonas syringae Strains Isolated from Grasses

The Pseudomonas syringae species complex has recently been named the number one plant pathogen, due to its economic and environmental impacts, as well as for its role in scientific research. The bacterium has been repeatedly reported to cause outbreaks on bean, cucumber, stone fruit, kiwi and olive tree, as well as on other crop and non-crop plants. It also serves as a model organism for research on the Type III secretion system (T3SS) and plant-pathogen interactions. While most of the current work on this pathogen is either carried out on one of three model strains found on dicot plants with completely sequenced genomes or on isolates obtained from recent outbreaks, not much is known about strains isolated from grasses (Poaceae). Here, we use comparative genomics in order to identify putative virulence-associated genes and other Poaceae-specific adaptations in several newly available genome sequences of strains isolated from grass species. All strains possess only a small number of known Type III effectors, therefore pointing to the importance of non-Type III secreted virulence factors. The implications of this finding are discussed.

Commonalities and differences of T3SSs in rhizobia and plant pathogenic bacteria

Plant pathogenic bacteria and rhizobia infect higher plants albeit the interactions with their hosts are principally distinct and lead to completely different phenotypic outcomes, either pathogenic or mutualistic, respectively. Bacterial protein delivery to plant host plays an essential role in determining the phenotypic outcome of plant-bacteria interactions. The involvement of type III secretion systems (T3SSs) in mediating animal- and plant-pathogen interactions was discovered in the mid-80's and is now recognized as a multiprotein nanomachine dedicated to trans-kingdom movement of effector proteins. The discovery of T3SS in bacteria with symbiotic lifestyles broadened its role beyond virulence. In most T3SS-positive bacterial pathogens, virulence is largely dependent on functional T3SSs, while in rhizobia the system is dispensable for nodulation and can affect positively or negatively the mutualistic associations with their hosts. This review focuses on recent comparative genome analyses in plant pathogens and rhizobia that uncovered similarities and variations among T3SSs in their genetic organization, regulatory networks and type III secreted proteins and discusses the evolutionary adaptations of T3SSs and type III secreted proteins that might account for the distinguishable phenotypes and host range characteristics of plant pathogens and symbionts.

Complete nucleotide sequence of the self-transmissible TOL plasmid pD2RT provides new insight into arrangement of toluene catabolic plasmids

In the present study we report the complete nucleotide sequence of the toluene catabolic plasmid pD2RT of Pseudomonas migulae strain D2RT isolated from Baltic Sea water. The pD2RT is 129,894 base pairs in size with an average G+C content of 53.75%. A total of 135 open reading frames (ORFs) were predicted to encode proteins, among them genes for catabolism of toluene, plasmid replication, maintenance and conjugative transfer. ORFs encoding proteins with putative functions in stress response, transposition and site-specific recombination were also predicted. Analysis of the organization and nucleotide sequence of pD2RT backbone region revealed high degree of similarity to the draft genome sequence data of the plant-pathogenic pseudomonad Pseudomonas syringae pv. glycinea strain B076, exhibiting relatedness to pPT23A plasmid family. The pD2RT backbone is also closely related to that of pGRT1 of Pseudomonas putida strain DOT-T1E and pBVIE04 of Burkholderia vietnamiensis strain G4, both plasmids are associated with resistance to toluene. The ability of pD2RT to self-transfer by conjugation to P. putida recipient strain PaW340 was experimentally determined. Genetic organization of toluene-degrading (xyl) genes and flanking DNA segments resembles the structure of Tn1721-related class II transposon Tn4656 of TOL plasmid pWW53 of P. putida strain MT53. The complete sequence of the plasmid pD2RT extends the known range of xyl genes carriers, being the first completely sequenced TOL plasmid, which is not related to well-studied IncP plasmid groups. We also verified the functionality of the catabolic route encoded by pD2RT by monitoring the expression of the xylE gene in pD2RT bearing hosts along with bacterial strains containing TOL plasmid of IncP-9 group. The growth kinetics of plasmid-bearing strains was found to be affected by particular TOL plasmid.

Non contiguous-finished genome sequence of Pseudomonas syringae pathovar syringae strain B64 isolated from wheat

The Gram-negative gammaproteobacterium Pseudomonas syringae is one of the most wide-spread plant pathogens and has been repeatedly reported to cause significant damage to crop plantations. Research on this pathogen is very intensive, but most of it is done on isolates that are pathogenic to Arabidopsis, tomato, and bean. Here, we announce a high-quality draft genome sequence of Pseudomonas syringae pv. syringae B64 which is the first published genome of a P. syringae strain isolated from wheat up to date. The genome sequence will assist in gaining insights into basic virulence mechanisms of this pathogen which has a relatively small complement of type III effectors.

Structure of microcin B-like compounds produced by Pseudomonas syringae and species specificity of their antibacterial action

Escherichia coli microcin B (Ec-McB) is a posttranslationally modified antibacterial peptide containing multiple oxazole and thiazole heterocycles and targeting the DNA gyrase. We have found operons homologous to the Ec-McB biosynthesis-immunity operon mcb in recently sequenced genomes of several pathovars of the plant pathogen Pseudomonas syringae, and we produced two variants of P. syringae microcin B (Ps-McB) in E. coli by heterologous expression. Like Ec-McB, both versions of Ps-McB target the DNA gyrase, but unlike Ec-McB, they are active against various species of the Pseudomonas genus, including human pathogen P. aeruginosa. Through analysis of Ec-McB/Ps-McB chimeras, we demonstrate that three centrally located unmodified amino acids of Ps-McB are sufficient to determine activity against Pseudomonas, likely by allowing specific recognition by a transport system that remains to be identified. The results open the way for construction of McB-based antibacterial molecules with extended spectra of biological activity.

A novel l-isoleucine-4'-dioxygenase and l-isoleucine dihydroxylation cascade in Pantoea ananatis

A unique operon structure has been identified in the genomes of several plant- and insect-associated bacteria. The distinguishing feature of this operon is the presence of tandem hilA and hilB genes encoding dioxygenases belonging to the PF13640 and PF10014 (BsmA) Pfam families, respectively. The genes encoding HilA and HilB from Pantoea ananatis AJ13355 were cloned and expressed in Escherichia coli. The culturing of E. coli cells expressing hilA (E. coli-HilA) or both hilA and hilB (E. coli-HilAB) in the presence of l-isoleucine resulted in the conversion of l-isoleucine into two novel biogenic compounds: l-4′-isoleucine and l-4,4′-dihydroxyisoleucine, respectively. In parallel, two novel enzymatic activities were detected in the crude cell lysates of the E. coli-HilA and E. coli-HilAB strains: l-isoleucine, 2-oxoglutarate: oxygen oxidoreductase (4′-hydroxylating) (HilA) and l-4′-hydroxyisoleucine, 2-oxoglutarate: oxygen oxidoreductase (4-hydroxylating) (HilB), respectively. Two hypotheses regarding the physiological significance of C-4(4′)-hydroxylation of l-isoleucine in bacteria are also discussed. According to first hypothesis, the l-isoleucine dihydroxylation cascade is involved in synthesis of dipeptide antibiotic in P. ananatis. Another unifying hypothesis is that the C-4(4′)-hydroxylation of l-isoleucine in bacteria could result in the synthesis of signal molecules belonging to two classes: 2(5H)-furanones and analogs of N-acyl homoserine lactone.

Comparative genomics of multiple strains of Pseudomonas cannabina pv. alisalensis, a potential model pathogen of both monocots and dicots

Comparative genomics of closely related pathogens that differ in host range can provide insights into mechanisms of host-pathogen interactions and host adaptation. Furthermore, sequencing of multiple strains with the same host range reveals information concerning pathogen diversity and the molecular basis of virulence. Here we present a comparative analysis of draft genome sequences for four strains of Pseudomonas cannabina pathovar alisalensis (Pcal), which is pathogenic on a range of monocotyledonous and dicotyledonous plants. These draft genome sequences provide a foundation for understanding host range evolution across the monocot-dicot divide. Like other phytopathogenic pseudomonads, Pcal strains harboured a hrp/hrc gene cluster that codes for a type III secretion system. Phylogenetic analysis based on the hrp/hrc cluster genes/proteins, suggests localized recombination and functional divergence within the hrp/hrc cluster. Despite significant conservation of overall genetic content across Pcal genomes, comparison of type III effector repertoires reinforced previous molecular data suggesting the existence of two distinct lineages within this pathovar. Furthermore, all Pcal strains analyzed harbored two distinct genomic islands predicted to code for type VI secretion systems (T6SSs). While one of these systems was orthologous to known P. syringae T6SSs, the other more closely resembled a T6SS found within P. aeruginosa. In summary, our study provides a foundation to unravel Pcal adaptation to both monocot and dicot hosts and provides genetic insights into the mechanisms underlying pathogenicity.

Pseudomonas syringae pv. tomato DC3000 CmaL (PSPTO4723), a DUF1330 family member, is needed to produce L-allo-isoleucine, a precursor for the phytotoxin coronatine

Pseudomonas syringae pv. tomato DC3000 produces the phytotoxin coronatine, a major determinant of the leaf chlorosis associated with DC3000 pathogenesis. The DC3000 PSPTO4723 (cmaL) gene is located in a genomic region encoding type III effectors; however, it promotes chlorosis in the model plant Nicotiana benthamiana in a manner independent of type III secretion. Coronatine is produced by the ligation of two moieties, coronafacic acid (CFA) and coronamic acid (CMA), which are produced by biosynthetic pathways encoded in separate operons. Cross-feeding experiments, performed in N. benthamiana with cfa, cma, and cmaL mutants, implicate CmaL in CMA production. Furthermore, analysis of bacterial supernatants under coronatine-inducing conditions revealed that mutants lacking either the cma operon or cmaL accumulate CFA rather than coronatine, supporting a role for CmaL in the regulation or biosynthesis of CMA. CmaL does not appear to regulate CMA production, since the expression of proteins with known roles in CMA production is unaltered in cmaL mutants. Rather, CmaL is needed for the first step in CMA synthesis, as evidenced by the fact that wild-type levels of coronatine production are restored to a ΔcmaL mutant when it is supplemented with 50 μg/ml l-allo-isoleucine, the starting unit for CMA production. cmaL is found in all other sequenced P. syringae strains with coronatine biosynthesis genes. This characterization of CmaL identifies a critical missing factor in coronatine production and provides a foundation for further investigation of a member of the widespread DUF1330 protein family.

Identification of novel virulence genes and metabolic pathways required for full fitness of Pseudomonas savastanoi pv. savastanoi in olive (Olea europaea) knots

Comparative genomics and functional analysis of Pseudomonas syringae and related pathogens have mainly focused on diseases of herbaceous plants; however, there is a general lack of knowledge about the virulence and pathogenicity determinants required for infection of woody plants. Here, we applied signature-tagged mutagenesis (STM) to Pseudomonas savastanoi pv. savastanoi during colonization of olive (Olea europaea) knots, with the goal of identifying the range of genes linked to growth and symptom production in its plant host. A total of 58 different genes were identified, and most mutations resulted in hypovirulence in woody olive plants. Sequence analysis of STM mutations allowed us to identify metabolic pathways required for full fitness of P. savastanoi in olive and revealed novel mechanisms involved in the virulence of this pathogen, some of which are essential for full colonization of olive knots by the pathogen and for the lysis of host cells. This first application of STM to a P. syringae-like pathogen provides confirmation of functional capabilities long believed to play a role in the survival and virulence of this group of pathogens but not adequately tested before, and unravels novel factors not correlated previously with the virulence of other plant or animal bacterial pathogens.

Phylogenetic analysis of a gene cluster encoding an additional, rhizobial-like type III secretion system that is narrowly distributed among Pseudomonas syringae strains

Background

The central role of Type III secretion systems (T3SS) in bacteria-plant interactions is well established, yet unexpected findings are being uncovered through bacterial genome sequencing. Some Pseudomonas syringae strains possess an uncharacterized cluster of genes encoding putative components of a second T3SS (T3SS-2) in addition to the well characterized Hrc1 T3SS which is associated with disease lesions in host plants and with the triggering of hypersensitive response in non-host plants. The aim of this study is to perform an in silico analysis of T3SS-2, and to compare it with other known T3SSs.

Results

Based on phylogenetic analysis and gene organization comparisons, the T3SS-2 cluster of the P. syringae pv. phaseolicola strain is grouped with a second T3SS found in the pNGR234b plasmid of Rhizobium sp. These additional T3SS gene clusters define a subgroup within the Rhizobium T3SS family. Although, T3SS-2 is not distributed as widely as the Hrc1 T3SS in P. syringae strains, it was found to be constitutively expressed in P. syringae pv phaseolicola through RT-PCR experiments.

Conclusions

The relatedness of the P. syringae T3SS-2 to a second T3SS from the pNGR234b plasmid of Rhizobium sp., member of subgroup II of the rhizobial T3SS family, indicates common ancestry and/or possible horizontal transfer events between these species. Functional analysis and genome sequencing of more rhizobia and P. syringae pathovars may shed light into why these bacteria maintain a second T3SS gene cluster in their genome.

Extensive remodeling of the Pseudomonas syringae pv. avellanae type III secretome associated with two independent host shifts onto hazelnut

BACKGROUND

Hazelnut (Corylus avellana) decline disease in Greece and Italy is caused by the convergent evolution of two distantly related lineages of Pseudomonas syringae pv. avellanae (Pav). We sequenced the genomes of three Pav isolates to determine if their convergent virulence phenotype had a common genetic basis due to either genetic exchange between lineages or parallel evolution.

RESULTS

We found little evidence for horizontal transfer (recombination) of genes between Pav lineages, but two large genomic islands (GIs) have been recently acquired by one of the lineages. Evolutionary analyses of the genes encoding type III secreted effectors (T3SEs) that are translocated into host cells and are important for both suppressing and eliciting defense responses show that the two Pav lineages have dramatically different T3SE profiles, with only two shared putatively functional T3SEs. One Pav lineage has undergone unprecedented secretome remodeling, including the acquisition of eleven new T3SEs and the loss or pseudogenization of 15, including five of the six core T3SE families that are present in the other Pav lineage. Molecular dating indicates that divergence within both of the Pav lineages predates their observation in the field. This suggest that both Pav lineages have been cryptically infecting hazelnut trees or wild relatives for many years, and that the emergence of hazelnut decline in the 1970s may have been due to changes in agricultural practice.

CONCLUSIONS

These data show that divergent lineages of P. syringae can converge on identical disease etiology on the same host plant using different virulence mechanisms and that dramatic shifts in the arsenal of T3SEs can accompany disease emergence.

The molecular basis of host specialization in bean pathovars of Pseudomonas syringae

Biotrophic phytopathogens are typically limited to their adapted host range. In recent decades, investigations have teased apart the general molecular basis of intraspecific variation for innate immunity of plants, typically involving receptor proteins that enable perception of pathogen-associated molecular patterns or avirulence elicitors from the pathogen as triggers for defense induction. However, general consensus concerning evolutionary and molecular factors that alter host range across closely related phytopathogen isolates has been more elusive. Here, through genome comparisons and genetic manipulations, we investigate the underlying mechanisms that structure host range across closely related strains of Pseudomonas syringae isolated from different legume hosts. Although type III secretion-independent virulence factors are conserved across these three strains, we find that the presence of two genes encoding type III effectors (hopC1 and hopM1) and the absence of another (avrB2) potentially contribute to host range differences between pathovars glycinea and phaseolicola. These findings reinforce the idea that a complex genetic basis underlies host range evolution in plant pathogens. This complexity is present even in host-microbe interactions featuring relatively little divergence among both hosts and their adapted pathogens.

Genomic Distribution and Divergence of Levansucrase-Coding Genes in Pseudomonas syringae

In the plant pathogenic bacterium, Pseudomonas syringae, the exopolysaccharide levan is synthesized by extracellular levansucrase (Lsc), which is encoded by two conserved 1,296-bp genes termed lscB and lscC in P. syringae strain PG4180. A third gene, lscA, is homologous to the 1,248-bp lsc gene of the bacterium Erwinia amylovora, causing fire blight. However, lscA is not expressed in P. syringae strain PG4180. Herein, PG4180 lscA was shown to be expressed from its native promoter in the Lsc-deficient E. amylovora mutant, Ea7/74-LS6, suggesting that lscA might be closely related to the E. amylovora lsc gene. Nucleotide sequence analysis revealed that lscB and lscC homologs in several P. syringae strains are part of a highly conserved 1.8-kb region containing the ORF, flanked by 450-452-bp and 49-51-bp up- and downstream sequences, respectively. Interestingly, the 450-452-bp upstream sequence, along with the initial 48-bp ORF sequence encoding for the N-terminal 16 amino acid residues of Lsc, were found to be highly similar to the respective sequence of a putatively prophage-borne glycosyl hydrolase-encoding gene in several P. syringae genomes. Minimal promoter regions of lscB and lscC were mapped in PG4180 by deletion analysis and were found to be located in similar positions upstream of lsc genes in three P. syringae genomes. Thus, a putative 498-500-bp promoter element was identified, which possesses the prophage-associated com gene and DNA encoding common N-terminal sequences of all 1,296-bp Lsc and two glycosyl hydrolases. Since the gene product of the non-expressed 1,248-bp lscA is lacking this conserved N-terminal region but is otherwise highly homologous to those of lscB and lscC, it was concluded that lscA might have been the ancestral lsc gene in E. amylovora and P. syringae. Our data indicated that its highly expressed paralogs in P. syringae are probably derived from subsequent recombination events initiated by insertion of the 498-500-bp promoter element, described herein, containing a translational start site.

Pseudomonas syringae pv. actinidiae draft genomes comparison reveal strain-specific features involved in adaptation and virulence to Actinidia species

A recent re-emerging bacterial canker disease incited by Pseudomonas syringae pv. actinidiae (Psa) is causing severe economic losses to Actinidia chinensis and A. deliciosa cultivations in southern Europe, New Zealand, Chile and South Korea. Little is known about the genetic features of this pathovar. We generated genome-wide Illumina sequence data from two Psa strains causing outbreaks of bacterial canker on the A. deliciosa cv. Hayward in Japan (J-Psa, type-strain of the pathovar) and in Italy (I-Psa) in 1984 and 1992, respectively as well as from a Psa strain (I2-Psa) isolated at the beginning of the recent epidemic on A. chinensis cv. Hort16A in Italy. All strains were isolated from typical leaf spot symptoms. The phylogenetic relationships revealed that Psa is more closely related to P. s. pv. theae than to P. avellanae within genomospecies 8. Comparative genomic analyses revealed both relevant intrapathovar variations and putative pathovar-specific genomic regions in Psa. The genomic sequences of J-Psa and I-Psa were very similar. Conversely, the I2-Psa genome encodes four additional effector protein genes, lacks a 50 kb plasmid and the phaseolotoxin gene cluster, argK-tox but has acquired a 160 kb plasmid and putative prophage sequences. Several lines of evidence from the analysis of the genome sequences support the hypothesis that this strain did not evolve from the Psa population that caused the epidemics in 1984-1992 in Japan and Italy but rather is the product of a recent independent evolution of the pathovar actinidiae for infecting Actinidia spp. All Psa strains share the genetic potential for copper resistance, antibiotic detoxification, high affinity iron acquisition and detoxification of nitric oxide of plant origin. Similar to other sequenced phytopathogenic pseudomonads associated with woody plant species, the Psa strains isolated from leaves also display a set of genes involved in the catabolism of plant-derived aromatic compounds.

Use of low-coverage, large-insert, short-read data for rapid and accurate generation of enhanced-quality draft Pseudomonas genome sequences

Next-generation genomic technology has both greatly accelerated the pace of genome research as well as increased our reliance on draft genome sequences. While groups such as the Genomics Standards Consortium have made strong efforts to promote genome standards there is a still a general lack of uniformity among published draft genomes, leading to challenges for downstream comparative analyses. This lack of uniformity is a particular problem when using standard draft genomes that frequently have large numbers of low-quality sequencing tracts. Here we present a proposal for an “enhanced-quality draft” genome that identifies at least 95% of the coding sequences, thereby effectively providing a full accounting of the genic component of the genome. Enhanced-quality draft genomes are easily attainable through a combination of small- and large-insert next-generation, paired-end sequencing. We illustrate the generation of an enhanced-quality draft genome by re-sequencing the plant pathogenic bacterium Pseudomonas syringae pv. phaseolicola 1448A (Pph 1448A), which has a published, closed genome sequence of 5.93 Mbp. We use a combination of Illumina paired-end and mate-pair sequencing, and surprisingly find that de novo assemblies with 100x paired-end coverage and mate-pair sequencing with as low as low as 2–5x coverage are substantially better than assemblies based on higher coverage. The rapid and low-cost generation of large numbers of enhanced-quality draft genome sequences will be of particular value for microbial diagnostics and biosecurity, which rely on precise discrimination of potentially dangerous clones from closely related benign strains.

Information Management of Genome Enabled Data Streams for Pseudomonas syringae on the Pseudomonas-Plant Interaction (PPI) Website

Genome enabled research has led to a large and ever-growing body of data on Pseudomonas syringae genome variation and characteristics, though systematic capture of this information to maximize access by the research community remains a significant challenge. Major P. syringae data streams include genome sequence data, newly identified type III effectors, biological characterization data for type III effectors, and regulatory feature characterization. To maximize data access, the Pseudomonas-Plant Interaction (PPI) website [1] is primarily focused on categorization of type III effectors and curation of effector functional data represented in the Hop database and Pseudomonas-Plant Interaction Resource, respectively. The PPI website further serves as a conduit for incorporation of new genome characterization data into the annotation records at NCBI and other data repositories, and clearinghouse for additional data sets and updates in response to the evolving needs of the research community.

Identification of Genes Involved in the Glycosylation of Modified Viosamine of Flagellins in Pseudomonas syringae by Mass Spectrometry

Previously we revealed that flagellin proteins in Pseudomonas syringae pv. tabaci 6605 (Pta 6605) were glycosylated with a trisaccharide, modified viosamine (mVio)-rhamnose-rhamnose and that glycosylation was required for virulence. We further identified some glycosylation-related genes, including vioA, vioB, vioT, fgt1, and fgt2. In this study, we newly identified vioR and vioM in a so-called viosamine island as biosynthetic genes for glycosylation of mVio in Pta 6605 by the mass spectrometry (MS) of flagellin glycan in the respective mutants. Furthermore, characterization of the mVio-related genes and MS analyses of flagellin glycans in other pathovars of P. syringae revealed that mVio-related genes were essential for mVio biosynthesis in flagellin glycans, and that P. syringae pv. syringae B728a, which does not possess a viosamine island, has a different structure of glycan in its flagellin protein.

Sequence and role in virulence of the three plasmid complement of the model tumor-inducing bacterium Pseudomonas savastanoi pv. savastanoi NCPPB 3335

Pseudomonas savastanoi pv. savastanoi NCPPB 3335 is a model for the study of the molecular basis of disease production and tumor formation in woody hosts, and its draft genome sequence has been recently obtained. Here we closed the sequence of the plasmid complement of this strain, composed of three circular molecules of 78,357 nt (pPsv48A), 45,220 nt (pPsv48B), and 42,103 nt (pPsv48C), all belonging to the pPT23A-like family of plasmids widely distributed in the P. syringae complex. A total of 152 coding sequences were predicted in the plasmid complement, of which 38 are hypothetical proteins and seven correspond to putative virulence genes. Plasmid pPsv48A contains an incomplete Type IVB secretion system, the type III secretion system (T3SS) effector gene hopAF1, gene ptz, involved in cytokinin biosynthesis, and three copies of a gene highly conserved in plant-associated proteobacteria, which is preceded by a hrp box motif. A complete Type IVA secretion system, a well conserved origin of transfer (oriT), and a homolog of the T3SS effector gene hopAO1 are present in pPsv48B, while pPsv48C contains a gene with significant homology to isopentenyl-diphosphate delta-isomerase, type 1. Several potential mobile elements were found on the three plasmids, including three types of MITE, a derivative of IS801, and a new transposon effector, ISPsy30. Although the replication regions of these three plasmids are phylogenetically closely related, their structure is diverse, suggesting that the plasmid architecture results from an active exchange of sequences. Artificial inoculations of olive plants with mutants cured of plasmids pPsv48A and pPsv48B showed that pPsv48A is necessary for full virulence and for the development of mature xylem vessels within the knots; we were unable to obtain mutants cured of pPsv48C, which contains five putative toxin-antitoxin genes.

Application of high-throughput genome sequencing to intrapathovar variation in Pseudomonas syringae

One reason for the success of Pseudomonas syringae as a model pathogen has been the availability of three complete genome sequences since 2005. Now, at the beginning of 2011, more than 25 strains of P. syringae have been sequenced and many more will soon be released. To date, published analyses of P. syringae have been largely descriptive, focusing on catalogues of genetic differences among strains and between species. Numerous powerful statistical tools are now available that have yet to be applied to P. syringae genomic data for robust and quantitative reconstruction of evolutionary events. The aim of this review is to provide a snapshot of the current status of P. syringae genome sequence data resources, including very recent and unpublished studies, and thereby demonstrate the richness of resources available for this species. Furthermore, certain specific opportunities and challenges in making the best use of these data resources are highlighted.

Comparative Genomics of Erwinia amylovora and Related Erwinia Species-What do We Learn?

Erwinia amylovora, the causal agent of fire blight disease of apples and pears, is one of the most important plant bacterial pathogens with worldwide economic significance. Recent reports on the complete or draft genome sequences of four species in the genus Erwinia, including E. amylovora, E. pyrifoliae, E. tasmaniensis, and E. billingiae, have provided us near complete genetic information about this pathogen and its closely-related species. This review describes in silico subtractive hybridization-based comparative genomic analyses of eight genomes currently available, and highlights what we have learned from these comparative analyses, as well as genetic and functional genomic studies. Sequence analyses reinforce the assumption that E. amylovora is a relatively homogeneous species and support the current classification scheme of E. amylovora and its related species. The potential evolutionary origin of these Erwinia species is also proposed. The current understanding of the pathogen, its virulence mechanism and host specificity from genome sequencing data is summarized. Future research directions are also suggested.

Dynamic evolution of pathogenicity revealed by sequencing and comparative genomics of 19 Pseudomonas syringae isolates

Closely related pathogens may differ dramatically in host range, but the molecular, genetic, and evolutionary basis for these differences remains unclear. In many Gram- negative bacteria, including the phytopathogen Pseudomonas syringae, type III effectors (TTEs) are essential for pathogenicity, instrumental in structuring host range, and exhibit wide diversity between strains. To capture the dynamic nature of virulence gene repertoires across P. syringae, we screened 11 diverse strains for novel TTE families and coupled this nearly saturating screen with the sequencing and assembly of 14 phylogenetically diverse isolates from a broad collection of diseased host plants. TTE repertoires vary dramatically in size and content across all P. syringae clades; surprisingly few TTEs are conserved and present in all strains. Those that are likely provide basal requirements for pathogenicity. We demonstrate that functional divergence within one conserved locus, hopM1, leads to dramatic differences in pathogenicity, and we demonstrate that phylogenetics-informed mutagenesis can be used to identify functionally critical residues of TTEs. The dynamism of the TTE repertoire is mirrored by diversity in pathways affecting the synthesis of secreted phytotoxins, highlighting the likely role of both types of virulence factors in determination of host range. We used these 14 draft genome sequences, plus five additional genome sequences previously reported, to identify the core genome for P. syringae and we compared this core to that of two closely related non-pathogenic pseudomonad species. These data revealed the recent acquisition of a 1 Mb megaplasmid by a sub-clade of cucumber pathogens. This megaplasmid encodes a type IV secretion system and a diverse set of unknown proteins, which dramatically increases both the genomic content of these strains and the pan-genome of the species.

Breakthroughs in genomics have unleashed a new suite of tools for studying the genetic bases of phenotypic differences across diverse bacterial isolates. Here, we analyze 19 genomes of P. syringae, a pathogen of many crop species, to reveal the genetic changes underlying differences in virulence across host plants ranging from rice to maple trees. Surprisingly, a pair of strains diverged dramatically via the acquisition of a 1 Mb megaplasmid, which constitutes roughly 14% of the genome. Novel plasmids and horizontal genetic exchange have contributed extensively to species-wide diversification. Type III effector proteins are essential for pathogenicity, exhibit wide diversity between strains and are present in distinct higher-level patterns across the species. Furthermore, we use sequence comparisons within an evolutionary context to identify functional changes in multiple virulence genes. Overall, our data provide a unique overview of evolutionary pressures within P. syringae and an important resource for the phytopathogen research community.